Properties of iron at the Earth's core conditions

Summary. The phase diagram of iron up to 330 GPa is solved using the experimental data of static high pressure (up to 11 GPa) and the experimental data of shock wave data (up to 250 GPa). A solution for the highest triple point is found ( P = 280 GPa and T = 5760 K) by imposing the thermodynamic constraints of triple points. This pressure of the triple point is less than the pressure of the inner core–outer core boundary of the Earth. These results indicate that the density of iron at the inner core–outer core boundary pressure is close to 13 g cm⁻³, which lies close to the seismic solutions of the Earth at that pressure. It is thus concluded that the Earth's inner core is very likely to be virtually pure iron in its hexagonal close packed (hcp) phase.
It is shown that four properties of the Earth's inner core determined from seismology are close in value to the corresponding properties of hcp iron at inner core conditions: density, bulk modulus, longitudinal velocity, and Poisson's ratio. The density–pressure profile of hcp iron at inner core conditions matches the density–pressure profile of the inner core as determined by seismic methods, within the spread of values given by recent seismic models.
This indicates that the Earth is slowly cooling, the Earth's inner core is growing by crystallization, and the impurities of the core are concentrated in the outer core. The calculated temperature at the Earth's centre is 6450 K.     

Coercive force of single crystals of magnetite at low temperatures

The temperature dependence of coercive force H _c was studied on well-characterized and stoichiometric millimetre-sized single crystals of magnetite at a series of 16 temperatures from 300 to 10 K using a SQUID magnetometer. H _c decreases gradually with cooling to the isotropic temperature, T _i = 130 K, where the first magnetocrystalline anisotropy constant K ₁ becomes zero. H _c exhibits a sharp increase at the Verwey transition, T _v = 120 K, where the structure changes from cubic to monoclinic. In crossing the Verwey transition, H _c increases by more than two orders of magnitude, from 20 μT to 2.4 mT, and the shape of the hysteresis loops becomes wasp-waisted.
Observed coercivity between 300 K and 170 K varies with temperature as λ _s / M _s , where λ _s is the magnetostriction constant and M _s is the saturation magnetization, indicating that the coercivity in MD magnetite is controlled mainly by internal stress associated with dislocations or other crystal defects. It seems likely that the stable single-domain-like magnetic memory observed in large MD magnetite crystals is due to magnetoelastically pinned domain walls. The discontinuous change in H _c at the Verwey transition is controlled by abrupt changes in magnetocrystalline and magnetostriction constants due to crystal deformation from cubic to monoclinic structure.     

The deformation and compaction of partial molten zones

Summary. The segregation of melt from a partially molten source region requires a corresponding deformation of the unmelted residue ('matrix'). The role of matrix deformation during melt segregation is examined using simple one-dimensional models, for which the deformation consists only of bulk compression or 'compaction'. In model I, a volume fraction φ₀ of ascending mantle material undergoes pressure-release melting at a depth z = 0 (localized melting). Compaction of the matrix occurs in a boundary layer whose thickness (reduced compaction length δ_R) is proportional to the square root of the matrix viscosity. In the Earth's mantle, δ_R∼ 10–100 m, indicating that compaction cannot be important over large distances. Model II examines the case in which melting occurs over a depth range of order h (distributed melting). In the limit h ≪δ_R, the solution is the same as for the case of localized melting, except in a 'melting layer' of thickness ∼ h near z = 0. In the more realistic limit h ≫δ_R, compaction makes a negligible contribution to the balance of forces associated with melt segregation. This result is also valid for the more general case of two-dimensional flow. Compaction is therefore likely to be of negligible importance in the Earth's mantle, with the consequence that melt segregation can be accurately described by Darcy's law.     

A diurnal resonance in the ocean tide and in the Earth's load response due to the resonant free 'core nutation'

Summary. The luni-solar forced nutations and body tide are believed to be resonant at frequencies near (1 + 1/460) cycle sidereal day⁻¹ as seen from the rotating Earth. This resonance is due to the Earth's rotating, elliptical fluid core. We show here that tides in the open ocean and the Earth's response to those tides must also be resonant at (1 + 1/460) cycle day⁻¹. We examine these resonant oceanic effects on the Earth's nutational motion and on the body tide. Effects on the forced nutations might be as large as 0.002 arcsec at 18.6 yr. The effects on the observed resonance in the body tide are more important. For tidal gravity, for example, the difference between K ₁ and 0 ¹ which is usually used to determine the resonance, can be perturbed by 30 per cent or more due to the oceanic resonance effects.     

A model for energy dissipation at the mantle—core boundary

Summary. An existing experimentally verified model for energy dissipation in a processing spherical cavity filled with liquid assumed to be in a semirigidized state except for a viscous Ekman boundary layer is applied to the Earth's liquid core to assess energy dissipation magnitudes. Application of the model to the best available Earth data occurs at the derived energy dissipation maximum for the model. Other existing research showing that the Earth's atmosphere appears to adjust to a state of maximum dissipation led to generic models for systems of maximum dissipation. The maximum dissipation mantle—core model with core motion driven by Earth precession alone, coupled to the mantle only by viscous shear stresses, and with a spherical mantle—core boundary leads to energy dissipation rates on the order of 10⁴ times those necessary for an Earth dynamo. The maximum dissipation model also leads to excessive magnetic field drift rates and to excessive retardation of the Earth's rotation rate. Effects of the mantle—core ellipticity and of magnetic field coupling are briefly discussed and are used to help develop a less than maximum dissipation model also driven by precession alone but using the additional coupling to yield a model more consistent with observed phenomena.     

A linear approximation to the solution of a one-dimensional Stefan problem and its geophysical implications

Summary. The motion of a phase boundary in the Earth caused by temperature and pressure excitations at the Earth's surface is determined under a linear approximation. The solution is found as a sum of convolutions of pressure and temperature Green's functions with the corresponding excitations. The Green's functions are given under the form of Laplace transforms that can be inverted either by numerical evaluation of a branch cut integral or by inversion of a series expansion. This solution is a generalization of a solution previously derived by Gjevik. This latter solution is the first term in the series expansion. The relaxation times associated with the phase boundary motion are of the order of 10⁵–10⁷yr for the olivine—spinel phase transition and of 10⁶–10⁷yr for the basalt—eclogite transition. The linear approximation remains valid for long times only if the phase boundary moves slowly.     

Dynamic behaviour of a phase boundary under non-uniform surface loads

Summary. A method is outlined to determine the dynamic behaviour of a phase boundary in the Earth when non-uniform time-varying pressure and temperature conditions are assumed at the Earth's surface. An integral equation describing the phase boundary motion is derived and it is solved under a linearizing assumption. The solution is obtained in the form of a double integral transform. Short and long time-expansions of the solution can be obtained from series expansion and integration of the Laplace transform along a branch cut. The method is illustrated by considering a stepwise change in surface pressure conditions.
For short times, the solution exhibits the same type of time dependence (i.e. the first-order term is in t ^1/2) as the one obtained in the one-dimensional case (i.e. uniform pressure perturbation at the Earth's surface).
For long times, it is shown that the time dependence of the phase boundary motion is almost identical to the one derived for the one- dimensional case if the wavenumber k _L of the surface excitation is such that κ k ²_Lτ≤ 1 (where τ is the relaxation time associated with the one-dimensional phase boundary motion and κ is the thermal diffusivity). If κ k ²_Lτ > 1, then the relaxation time for the phase boundary motion in two dimensions is of the order of κ⁻¹ k ⁻²_L.
When considering parameters that would be appropriate for a basalt to eclogite phase transition at Moho depth, the latter situation is met only when the load wavelength is smaller than 35 km.     

Magnetic and viscous coupling at the core–mantle boundary: inferences from observations of the Earth's nutations

Dissipative core–mantle coupling is evident in observations of the Earth's nutations, although the source of this coupling is uncertain. Magnetic coupling occurs when conducting materials on either side of the boundary move through a magnetic field. In order to explain the nutation observations with magnetic coupling, we must assume a high (metallic) conductivity on the mantle side of the boundary and a rms radial field of 0.69 mT. Much of this field occurs at short wavelengths, which cannot be observed directly at the surface. High levels of short-wavelength field impose demands on the power needed to regenerate the field through dynamo action in the core. We use a numerical dynamo model from the study of Christensen & Aubert (2006) to assess whether the required short-wavelength field is physically plausible. By scaling the numerical solution to a model with sufficient short-wavelength field, we obtain a total ohmic dissipation of 0.7–1 TW, which is within current uncertainties. Viscous coupling is another possible explanation for the nutation observations, although the effective viscosity required for this is 0.03 m² s⁻¹ or higher. Such high viscosities are commonly interpreted as an eddy viscosity. However, physical considerations and laboratory experiments limit the eddy viscosity to 10⁻⁴ m² s⁻¹, which suggests that viscous coupling can only explain a few percent of the dissipative torque between the core and the mantle.     

Temperature Dependence of Single-Crystal Spinel (MgAl2O4) Elastic Constants from 293 to 423°K Measured by Light-Sound Scattering in the Raman-Nath Region

Summary. The temperature dependence of single-crystal elastic constants of synthetic stoichiometric MgAl₂O₄ spinel has been measured by the light-sound scattering technique in the Raman-Nath region. The crystal is set into forced vibration by a single crystal LiNbO₃ transducer coupled to one crystal face. A He-Ne Laser beam is diffracted by the stress-induced birefringence inside the crystal. The diffraction angle is determined from the distance between two spots exposed on a photographic plate by the first order diffracted beams as measured by a microdensitometer. The sound wavelength inside the crystal is then inferred from the laser diffraction angle. Combining the sound wavelength with the measured transducer frequency, the velocity inside the crystal is determined typically to a precision of 0·05 per cent. In this method, the measurement of velocity is not dependent on either the determination of sample length or on phase shifts at sample-transducer interface. Velocities of four pure modes, L //[001], T //[001], L //[110], and T //[110]( P //[1 1 0] are measured in the temperature range between 293 and 423 °K. A linear temperature dependence is fit to the data by a least square method. Values obtained at 25 °C from this linear fit are
The temperature dependence of the adiabatic elastic constants and bulk and shear (VRH average) moduli is computed using the density and literature value of thermal expansion coefficient. Values obtained are:
A comparison with previous measurements by pulse superposition and ultrasonic interferometry methods is made. Disagreement, when present, is discussed in terms of the separate measuring techniques. Finally, the present method, with its possibility for further improvement, is evaluated as a new method to measure temperature and pressure dependence of elastic constants.     

Physical oceanography studies in the Weddell Sea during the Norwegian Antarctic Research Expedition 1978/79

Hydrographic and current measurements obtained during the Norwegian Antarctic Research Expedition 1978/79 to the southern Weddell Sea are presented. Cold, dense Ice Shelf Water circulating under the floating ice shelves is observed to leave the shelf as a concentrated bottom flow. From moored current metres this discharge is estimated at 0.7 10⁶ m³/s at -2.0°C (one year average) and with no appreciable seasonal variation. This contribution to the Weddell Sea Bottom Water is clearly identified through extreme temperature gradients at our deepest stations (below 2500 m). The core of Weddell Deep Water shows a considerable (T ∼ 0.5°C) warming up since 1977, presumably due to the lack of polynya activity in the intervening period. Measurements in the coastal current at the ice shelf (70°S, 2°W) show step structures which are probably due to cooling and melting at the vertical ice barrier. Slight supercooling due to circulation under the ice shelf is also seen. The net effect of the ice shelf boundary seems to be a deep reaching cooling and freshening of the coastal current providing the low salinity, freezing point Eastern Shelf Water. This process is considered a preconditioning which enhances production of the saline Western Shelf Water which in turn is transformed to Ice Shelf Water.     

The subseismic approximation in core dynamics—II. Love numbers and surface gravity

This paper extends our earlier examinations of the utility of various approximations for treating the dynamics of the Earth's liquid core on time-scales of the order of 10⁴ to 10⁸ s. We discuss the effects of representing the response of the mantle and inner core by static (versus dynamic) Love numbers, and of invoking the subseismic approximation for treating core flow, used either only in the interior of the liquid core (SSA-1) or also at the boundaries (SSA-2). The success of each approximation (or combinations thereof) is measured by comparing the resulting surface gravity effects (computed for a given earthquake excitation), and (for the Slichter mode) the distribution of translational momentum, with reference calculations in which none of these approximations is made. We conclude that for calculations of the Slichter triplet, none of the approximations is satisfactory, i.e. a full solution (using dynamic Love numbers at elastic boundaries and no core flow approximation) is required in order to avoid spurious eigenfrequencies and to yield correct eigenfunctions (e.g. conserving translational momentum) and surface gravity. For core undertones, the use of static Love numbers at rigid boundaries is acceptable, along with SSA-1 (i.e. provided the subseismic approximation is not invoked at the core boundaries). Although the calculations presented here are for a non-rotating earth model, we argue that the principal conclusions should be applicable to the rotating Earth. Shortcomings of the subseismic approximation appear to arise because both SSA-1 and SSA-2 lower the order of the governing system of differential equations (giving rise to a singular perturbation problem), and because SSA-2 overdetermines the boundary conditions (making it impossible for solutions to satisfy all continuity requirements at core boundaries).     

Transient and impulse responses of a one-dimensional linearly attenuating medium—II. A parametric study

Summary. We investigate one-dimensional waves in a standard linear solid for geophysically relevant ranges of the parameters. The critical parameters are shown to be T*= t_u/Q_m where t _u is the travel time and Q_m the quality factor in the absorption band, and τ⁻¹ _m , the high-frequency cut-off of the relaxation spectrum. The visual onset time, rise time, peak time, and peak amplitude are studied as functions of T* and τ _m. For very small τ _m , this model is shown to be very similar to previously proposed attenuation models. As τ _m grows past a critical value which depends on T* , the character of the attenuated pulse changes. Seismological implications of this model may be inferred by comparing body wave travel times with a'one second'earth model derived from long-period observations and corrected for attenuation effects assuming a frequency independent Q over the seismic band. From such a comparison we speculate that there may be a gap in the relaxation spectrum of the Earth's mantle for relaxation times shorter than about one second. However, observational constraints from the attenuation of body waves suggest that such a gap might in fact occur at higher frequencies. Such a hypothesis would imply a frequency dependence of Q in the Earth's mantle for short-period body waves.     

Numerical models of mantle convection with secular cooling

A 2-D time-dependent finite-difference numerical model is used to investigate the thermal character and evolution of a convecting layer which is cooling as it convects. Two basic cooling modes are considered: in the first, both upper and lower boundaries are cooled at the same rate, while maintaining the same temperature difference across the layer; in the second, the lower boundary temperature decreases with time while the upper boundary temperature is fixed at 0°C. The first cooling mode simulates the effects of internal heating while the second simulates planetary cooling as mantle convection extracts heat from, and thereby cools, the Earth's core. The mathematical analogue between the effects of cooling and internal heating is verified for finite-amplitude convection. It is found that after an initial transient period the central core of a steady but vigorous convection cell cools at a constant rate which is governed by the rate of cooling of the boundaries and the viscosity structure of the layer. For upper-mantle models the transient stage lasts for about 30 per cent of the age of the Earth, while for the whole mantle it lasts for longer than the age of the Earth. Consequently, in our models the bulk cooling of the mantle lags behind the cooling of the core-mantle boundary. Models with temperature-dependent viscosity are found to cool in the same manner as models with depth-dependent viscosity; the rate of cooling is controlled primarily by the horizontally averaged variation of viscosity with depth. If the Earth's mantle cools in a similar fashion, secular cooling of the planet may be insensitive to lateral variations of viscosity.     

Travel-time anomalies in the mantle under the North Atlantic

Summary. Lateral heterogeneity exists in the Earth's mantle, and may result in seismic velocity anomalies up to several per cent. If convection cells and plumes extend down to the core, then these features may be associated with local inhomogeneities observed in the lower mantle.
Published data for direct and core-reflected P -wave residuals are used to delineate velocity anomalies in the middle—lower mantle under the North Atlantic. Differential ( PcP — P ) residuals indicate travel-time anomalies near the core—mantle transition, and may be due to core topography or lateral variations in velocity. It is assumed that the anomalies occur near the midpoints of the ray paths. The main source of error in the data set may arise from phases which have been identified incorrectly. Hence trend-surfaces are fitted to the residual data to show only the large-scale trends in anomaly values, with wavelengths of the order of 1000 km.
The Azores and Colorado hot spots occur in a region covered by the data. A possible interpretation of the trend maps is that an anomalous zone extends from a relatively fast region at the core boundary at 35° N, 50° W up to these hot spots, at about 30 degrees from the vertical. This may agree with the suggestion of Anderson that plumes are chemical rather than thermal in origin. If inclⁱned plumes do exist, the deviation from the ideal vertical plume or convection cell boundary may imply that lateral shear or other distortion effects exist in the mantle.     

The causes of low-temperature demagnetization of remanence in multidomain magnetite

In this paper, the mechanisms of low-temperature demagnetization of remanence in multidomain magnetite are considered. New experimental observations of the behaviour of the saturation isothermal remanence, thermoremanence and partial thermoremanence at low temperatures are presented. The results show that there are two main contributions to this low-temperature demagnetization. The first and predominant contribution (type-1 demagnetization) is due to 'kinematic' domain state reorganization and occurs throughout cooling from room temperature to the Verwey transition, T _v , at 120–124  K. The second contribution arises from the change in anisotropy from cubic to monoclinic at T _v , which changes the overall domain structure of the grain. On warming in zero field, some domain walls will not return to their original positions but will take up a position that leads to a lower net remanence (type-2 demagnetization). In stoichiometric magnetite, demagnetization does not occur at 130  K at the isotropic point, T _k , contrary to some previous predictions. In non-stoichiometric magnetite, the influence of the Verwey transition is greatly reduced, and anomalous behaviour is observed at T _k .     

The secular acceleration of the Earth's spin

Summary. The power spectrum of the Earth's spin has important components with periods ranging from a few days to at least a few thousand years, and probably to the age of the Earth. The secular acceleration, as the term is used here, refers to the components with periods longer than three centuries. In the year 600, the secular acceleration was —19.9 ± 0.8 parts in 10⁹ per century, while the value at the present time is less than half this size. The spin acceleration has important contributions from tidal friction and from an effect that is proportional to the square of the magnetic dipole moment. When these contributions are subtracted from the observed acceleration, we are left with a contribution that amounts to +41 parts in 10⁹ per century. This amount probably results from an unknown combination of changes in the size of the core, in the amount of glaciation, and in the size of the gravitational constant.     

Experiments on precessing flows in the Earth's liquid core

Experiments simulating flow in the Earth's liquid core induced by luni-solar precession of the solid mantle indicate, to a first approximation, that the core behaves like a rigidized fluid sphere spinning slower than the mantle and with its spin axis lagging the mantle spin axis in precession. Secondary flow patterns are always present. At low precession rates the fluid sphere is subdivided into a set of cylinders coaxial with the fluid spin axis, the cylinders rotating alternately at slightly faster and slower rates relative to the net retrograde motion of the fluid as a whole. Slow non-axisymmetric columnar wave patterns develop between the differentially rotating cylinders. Axial flows between the spheroidal cavity boundary and the interior are observed. Fluid motion becomes turbulent only at precession rates large enough to cause the fluid spin axis to align nearly with the precession axis. There is no evidence that the Earth's liquid spin axis direction departs more than a fraction of a degree from geographic north. Our observations suggest precession induces a complex variety of laminar flows, including slowly varying and/or periodic patterns, in the Earth's liquid core.     

On the upper limit of undertone periods

Summary. The severity of the effect of the Earth's rotation on the transverse motions within the fluid core of the Earth is examined. It is shown that, for any model of the fluid core, the formalism adapted from free oscillation theory is applicable only for periods less than 12 hr. In the limiting case, as the frequency is decreased to 2 cycle day^-1, all modes of response in the fluid are excited to the same order.