Core Mantle Research Articles

Lost City’s Plumbing Exposed by the Longest Mantle Core Ever Drilled

Lost City’s Plumbing Exposed by the Longest Mantle Core Ever Drilled

Superionic iron hydride shapes ultralow-velocity zones at Earth’s core–mantle boundary

Seismological studies have exposed numerous ultralow velocity zones (ULVZs) exhibiting extraordinary physical attributes at Earth's core-mantle boundary, yet their compositions and origins remain controversial. Water-iron reaction can generate unique phases under lowermost-mantle conditions and likely plays a crucial role in forming ULVZs. Through first-principles molecular dynamic simulations with machine learning techniques, we determine that iron hydride, the product of water-iron reaction, is stable as a superionic phase at the core-mantle boundary. This superionic iron hydride has much slower velocities and a higher density than the ambient mantle under lowermost-mantle conditions. Accumulation of iron hydride, created through either a chemical reaction between subducted water and iron or solidification of core material entrained in the lower mantle by convection, can explain the seismic observations of ULVZs particularly those associated with subduction. This work suggests that water may have a substantial role in creating seismic heterogeneities at the core-mantle boundary.

Consequences of Continental Core Complexes on Rifting Patterns: Insights From Multichannel Seismic Data From the Northern Margin of the South China Sea

AbstractCore complex type structures are common in continental margins. Recent studies have proposed exhumation of deeper crustal and mantle levels during extension of the continental margin along the Northern South China Sea (NSCS) in a core complex mode. Here, we present high‐resolution seismic profiles across different sections of the NSCS that allow us to observe two types of core complex systems. The Kaiping Rift and eastern Baiyun Rift are characterized by a sub‐horizontal Moho and exhumation of ductile deeper crust beneath an extremely thinned brittle upper crust. The central Baiyun Rift in contrast is marked by homogeneous thinning of the entire crust that is floored by mantle unroofing with a concave‐down Moho geometry. While the former type, referred to as crustal core complexes, coincides with relatively intense syn‐rift magmatism, the latter, referred to as mantle core complexes, is marked by limited magmatism. Consistently, basin subsidence rates are lower over crustal core complexes than mantle core complexes. Compared to mantle core complexes at the magma‐poor Iberian margin, the widely distributed crustal core complexes in the NSCS indicate higher geotherms and more intense magmatism during extension of the continental lithosphere. We propose that the continental core complexes observed at the NSCS represent a typical outcome of extension of a gravitationally unstable lithosphere.

A New Database of Giant Impacts over a Wide Range of Masses and with Material Strength: A First Analysis of Outcomes

In the late stage of terrestrial planet formation, planets are predicted to undergo pairwise collisions known as giant impacts. Here, we present a high-resolution database of giant impacts for differentiated colliding bodies of iron–silicate composition, with target masses ranging from 1 × 10−4 M ⊕ up to super-Earths (5 M ⊕). We vary the impactor-to-target mass ratio, core–mantle (iron–silicate) fraction, impact velocity, and impact angle. Strength in the form of friction is included in all simulations. We find that, due to strength, the collisions with bodies smaller than about 2 ×10−3 M ⊕ can result in irregular shapes, compound-core structures, and captured binaries. We observe that the characteristic escaping velocity of smaller remnants (debris) is approximately half of the impact velocity, significantly faster than currently assumed in N-body simulations of planet formation. Incorporating these results in N-body planet formation studies would provide more realistic debris–debris and debris–planet interactions.

Crustal and Uppermost Mantle Azimuthal Seismic Anisotropy of Antarctica From Ambient Noise Tomography

AbstractSeismic anisotropy provides essential information for characterizing the orientation of deformation and flow in the crust and mantle. The isotropic structure of the Antarctic crust and upper mantle has been determined by previous studies, but the azimuthal anisotropy structure has only been constrained by mantle core phase (SKS) splitting observations. This study determines the azimuthal anisotropic structure of the crust and mantle beneath the central and West Antarctica based on 8—55 s Rayleigh wave phase velocities from ambient noise cross‐correlation. An anisotropic Rayleigh wave phase velocity map was created using a ray—based tomography method. These data are inverted using a Bayesian Monte Carlo method to obtain an azimuthal anisotropy model with uncertainties. The azimuthal anisotropy structure in most of the study region can be fit by a two‐layer structure, with one layer at depths of 0–15 km in the shallow crust and the other layer in the uppermost mantle. The azimuthal anisotropic layer in the shallow crust of West Antarctica, where it coincides with strong positive radial anisotropy quantified by the previous study, shows a fast direction that is subparallel to the inferred extension direction of the West Antarctic Rift System. Fast directions of upper mantle azimuthal anisotropy generally align with teleseismic shear wave splitting fast directions, suggesting a thin lithosphere or similar lithosphere‐asthenosphere deformation. However, inconsistencies in this exist in Marie Byrd Land, indicating differing ancient deformation patterns in the shallow mantle lithosphere sampled by the surface waves and deformation in the deeper mantle and asthenosphere sampled more strongly by splitting measurements.

Friction Experiments on Lunar Analog Gouges and Implications for the Mechanism of the Apollo 17 Long Runout Landslide

AbstractThe Light Mantle landslide is a hypermobile landslide on the Moon. Apollo 17 astronauts collected a core sample of the top 60 cm of the Light Mantle deposit, which is currently being analyzed as part of the NASA's Apollo Next Generation Sample Analysis program. The origin of its hypermobility remains undetermined, as the proposed mechanisms are difficult to prove because of the lack of theoretical and experimental support and the scarcity of field data related to the internal structures of its deposit. Regardless of the emplacement mechanisms, it has been proposed that localized dynamic frictional weakening is responsible for the early stage instability that leads to catastrophic failure. Here, we conduct friction experiments under vacuum to investigate the viability of dynamic friction weakening in lunar analog anorthosite‐bearing gouges (i.e., rock powders). Our results show that localized dynamic friction weakening does not occur in these gouges at loading conditions where, instead, weakening is observed in other materials on Earth. Therefore, possibly other fluidization‐related mechanisms contributed to the initiation of the hypermobile Light Mantle landslide. Finally, we describe the microstructures formed in the experiments, including the presence of clast cortex aggregates. Preliminary investigation of the Light Mantle core samples (73001/73002) shows the presence of similar microstructures. Therefore, our microstructural observations will help the analysis and interpretation of the Apollo 17 core samples, as keys to insights about internal processes occurring during the emplacement of the landslide.

Kilometer-scale structure on the core–mantle boundary near Hawaii

The lowermost mantle right above the core-mantle boundary is highly heterogeneous containing multiple poorly understood seismic features. The smallest but most extreme heterogeneities yet observed are ‘Ultra-Low Velocity Zones’ (ULVZ). We exploit seismic shear waves that diffract along the core-mantle boundary to provide new insight into these enigmatic structures. We measure a rare core-diffracted signal refracted by a ULVZ at the base of the Hawaiian mantle plume at unprecedentedly high frequencies. This signal shows remarkably longer time delays at higher compared to lower frequencies, indicating a pronounced internal variability inside the ULVZ. Utilizing the latest computational advances in 3D waveform modeling, here we show that we are able to model this high-frequency signal and constrain high-resolution ULVZ structure on the scale of kilometers, for the first time. This new observation suggests a chemically distinct ULVZ with increasing iron content towards the core-mantle boundary, which has implications for Earth’s early evolutionary history and core-mantle interaction.

Fractional crystallization of a basal lunar magma ocean: A dense melt-bearing garnetite layer above the core?

In the wake of its accretion, the Moon was likely partially or fully molten, forming a lunar magma ocean (LMO). A fully molten Moon may crystallize neutrally buoyant olivine, whose accumulation may form a barrier leading to a separation into two independently evolving melt reservoirs. This study investigates the crystallization of such a putative basal lunar magma ocean, ranging from the core mantle boundary (with radius r = 380 km) to the level of neutral olivine buoyancy (r = 600 km). Consecutive crystallization experiments determine liquidus temperatures and crystallize 10–30 wt% minerals, before conceptually segregating these according to their buoyancy, and then stepping to a new bulk composition that corresponds to the residual melt after removal of the cumulate minerals.The first olivine to crystallize from a Taylor Whole Moon composition (twm, XMg = Mg/(Mg + Fe2+) = 0.83) has XMg = 0.94 and is neutrally buoyant at 3.8 GPa according to the melt density model of Lange and Carmichael (1990). Crystallization begins at the core mantle boundary, but early olivine floats and re-dissolves until the magma ocean cools to the liquidus temperature at the depth of neutral buoyancy (1850 °C). At this point a > 500 km thick olivine-only layer could form, mainly fed from the upper magma shell. The crystallization sequence in the basal magma ocean is olivine-only at 1850–1675 °C, 0–26 pcs (wt% percent solidified) → olivine + opx (to 1600 °C, 42 pcs) → opx + cpx + garnet (to 1580 °C, 60 pcs) → cpx + garnet (to 1520 °C, 73 pcs) → cpx + garnet + olivine leaving 11 wt% residual melt at 1450 °C.Cooling this last melt to 1300 °C leads to 80% garnet + cpx + olivine + Ti-spinel, the residual liquid corresponding to 2.2% of the basal magma ocean. Olivine, opx and cpx remain buoyant over the entire crystallization interval and various pyroxenite layers are added to the olivine layer. Garnet and the final cotectic cpx + garnet + olivine + FeTi-oxide assemblage instead form a 70 km thick basal layer on the core-mantle boundary. Such a layer provides a gravitationally stable high-density lowermost lunar mantle, which would concur with a recent re-analysis of the lunar seismic data. The lowest temperature experiment at 1300 °C, not much above the 1250 °C proposed for the core-mantle boundary from inversion of geophysical data, had about 20% of a highly evolved, dense Fe-rich melt. It is feasible that this melt has remained in its liquid state to present day, providing an explanation for the proposed low vp and vs layer and high dissipation just above the core-mantle boundary.

Earth and Mars – Distinct inner solar system products

Composition of terrestrial planets records planetary accretion, core–mantle and crust–mantle differentiation, and surface processes. Here we compare the compositional models of Earth and Mars to reveal their characteristics and formation processes. Earth and Mars are equally enriched in refractory elements (1.9 × CI), although Earth is more volatile-depleted and less oxidized than Mars. Their chemical compositions were established by nebular fractionation, with negligible contributions from post-accretionary losses of moderately volatile elements. The degree of planetary volatile element depletion might correlate with the abundances of chondrules in the accreted materials, planetary size, and their accretion timescale, which provides insights into composition and origin of Mercury, Venus, the Moon-forming giant impactor, and the proto-Earth. During its formation before and after the nebular disk's lifetime, the Earth likely accreted more chondrules and less matrix-like materials than Mars and chondritic asteroids, establishing its marked volatile depletion. A giant impact of an oxidized, differentiated Mars-like (i.e., composition and mass) body into a volatile-depleted, reduced proto-Earth produced a Moon-forming debris ring with mostly a proto-Earth's mantle composition. Chalcophile and some siderophile elements in the silicate Earth added by the Mars-like impactor were extracted into the core by a sulfide melt (∼0.5% of the mass of the Earth's mantle). In contrast, the composition of Mars indicates its rapid accretion of lesser amounts of chondrules under nearly uniform oxidizing conditions. Mars’ rapid cooling and early loss of its dynamo likely led to the absence of plate tectonics and surface water, and the present-day low surface heat flux. These similarities and differences between the Earth and Mars made the former habitable and the other inhospitable to uninhabitable.

Collisions between ice-covered silica grains: An atomistic study

We use molecular dynamics simulation to study collisions of silica grains covered by an ice mantle. Such heterogeneous particles may belong to debris disks of planetary systems below the snow line. For hydroxylated silica spheres, our results of the bouncing velocity obtained for nm-sized grains extrapolate well to experimental results for μm-sized particles. However, with increasing thickness of the ice mantle, the bouncing velocity strongly increases while the coefficient of restitution is reduced, meaning that the inelastic losses increase. An analysis of the processes occurring during the collision of two core–mantle grains shows that the ice mantle in the collision zone is strongly heated during the collision-induced compression. Temperatures surpass the triple point of water such that the ice may melt and become soft enough to yield sideways during the collision, dissipating collision energy. The ice thus acts both as a ‘cushion’ to soften the collision dynamics and as a glue to bind the two grains together. As a consequence of the energy dissipation, the bouncing velocity increases with increasing mantle thickness. The consequences of our findings on collisions of larger ice-coated grains are discussed. We conclude that collision-induced heating constitutes an important process in the collision of ice-covered grains, since it changes the collision physics from a ‘dry’ to a ‘wet’ contact.

Grain boundary diffusion of W in lower mantle phase with implications for isotopic heterogeneity in oceanic island basalts by core-mantle interactions

Tungsten isotopes provide important constraints on the ocean-island basalt (OIB) source regions. Recent analyses of μ182W in modern basalts with high 3He/4He originating from the core-mantle boundary region reveal two distinct features: positive μ182W in Phanerozoic flood basalts indicating the presence of primordial reservoir, and negative μ182W in modern OIBs. One possibility to produce large variations in μ182W is interaction between the mantle and outer core. Here, we report grain boundary diffusion of W in lower mantle phases. High pressure experimental results show that grain boundary diffusion of W is fast and strongly temperature dependent. Over Earth's history, diffusive transport of W from the core to the lowermost mantle may have led to significant modification of the W isotopic composition of the lower mantle at length scales exceeding one kilometer. Such grain boundary diffusion can lead to large variations in μ182W in modern basalts as a function of the distance of their source regions from the core mantle boundary. Modern oceanic island basalts from Hawaii, Samoa and Iceland exhibit negative μ182W and likely originated from the modified isotope region just above the core-mantle boundary, whereas those with positive μ182W could be derived from the thick Large Low Shear Velocity Provinces (LLSVPs) far from the core-mantle boundary (CMB). When highly-oxidized slabs accumulate at the CMB oxidizing the outer core at the interface, a large W flux with negative μ182W can be added to the silicate mantle. As a result, the source region of the OIB would be effectively modified to a negative μ182W.

Hypothetical Physics and Chemistry of Volcanic Eruptions: The Doorway to Their Prediction

This article presents a further development of the hypotheses concerning the possibility of predicting (“tectonic”) earthquakes [1]. Those hypotheses are based on the conversion of all types of released energy into heat and active chemical substances. One of the important sources of this phenomenon is the release of the latent energy trapped and stored during the Earth’s accretion. The latent energy of primordial hydrogen and helium escaping from the Earth’s core and lower mantle causes degassing processes [2] [3]. This latent energy converts into totally different types of chemical, electromagnetic and thermal energies of active compounds that are responsible for the major endogenic terrestrial processes. The dominating theories in seismology and volcanology are that an earthquake results from a sudden slip of a tectonic fault and that only magma and the gases contained in magma supply the volcanic energy resulting in the conclusions that earthquakes and eruptions are unpredictable. Volcanic eruption is considered herein to be a special case of the earthquake-process in which earthquake hypocenters rise to the Earth’s surface. A possible solution is proposed ([1] and herein) based on the analyses of the physicochemical processes as participants in earthquake and eruption preparations (foreshocks - major shock - aftershocks - volcanic eruptions) and on the characteristic rates of reflection of these processes on the Earth’s surface. Influences of Sun-Moon-tides and volcanic (“harmonic”) tremors are analyzed from physical-chemical point of view. The case of the 1980 eruption of Mount St. Helens and the proposed monitoring of the recommended additional data provides a way of selecting a complex of reliable earthquake and volcanic eruption precursors.

Prospects for an ancient dynamo and modern crustal remanent magnetism on Venus

Venus lacks an internally generated magnetic field today. Whether one existed in the past is unknown, but critical to atmospheric evolution and potential habitability. Canonical models assume the core of Venus is cooling too slowly for convection and thus a magnetic dynamo to occur today. Core/mantle heat flow is suppressed in these models after a putative transition in mantle dynamics associated with widespread, volcanic resurfacing. However, recent studies of impact craters and other surface features support more steady heat loss over geologic time. Precipitation of MgO and/or SiO2 from the core can also drive compositional convection even with slow cooling. Here we reevaluate the likelihood that Venus has an “Earth-like” (at least partially liquid and chemically homogeneous) core using numerical simulations of the coupled atmosphere–surface–mantle–core evolution. An Earth-like core is only compatible with the modern lack of a dynamo if the thermal conductivity of core material is towards the higher end of modern estimates (i.e., >100 W m−1 K−1). If lower estimates like ∼40–50 W m−1 K−1 are actually correct, then we favor recent proposals that the core has completely solidified or preserved primordial stratification. Any simulation initialized with a homogeneous, liquid core predicts a global magnetic field with Earth-like surface strength for >2–3 billion years after accretion—consistent with all available observations—and also sporadic activity within the surface age while temperatures remain below the Curie point of magnetite. Therefore, future spacecraft missions should prioritize the first-ever magnetometer measurements below the ionosphere to search for crustal remanent magnetism.

The Thermal Conductivity of Earth's Core: A Key Geophysical Parameter's Constraints and Uncertainties

The thermal conductivity of iron alloys at high pressures and temperatures is a critical parameter in governing ( a) the present-day heat flow out of Earth's core, ( b) the inferred age of Earth's inner core, and ( c) the thermal evolution of Earth's core and lowermost mantle. It is, however, one of the least well-constrained important geophysical parameters, with current estimates for end-member iron under core-mantle boundary conditions varying by about a factor of 6. Here, the current state of calculations, measurements, and inferences that constrain thermal conductivity at core conditions are reviewed. The applicability of the Wiedemann-Franz law, commonly used to convert electrical resistivity data to thermal conductivity data, is probed: Here, whether the constant of proportionality, the Lorenz number, is constant at extreme conditions is of vital importance. Electron-electron inelastic scattering and increases in Fermi-liquid-like behavior may cause uncertainties in thermal conductivities derived from both first-principles-associated calculations and electrical conductivity measurements. Additional uncertainties include the role of alloying constituents and local magnetic moments of iron in modulating the thermal conductivity. Thus, uncertainties in thermal conductivity remain pervasive, and hence a broad range of core heat flows and inner core ages appear to remain plausible.

Contamination tracer testing with seabed drills: IODP Expedition 357

Abstract. IODP Expedition 357 utilized seabed drills for the first time in the history of the ocean drilling program, with the aim of collecting intact sequences of shallow mantle core from the Atlantis Massif to examine serpentinization processes and the deep biosphere. This novel drilling approach required the development of a new remote seafloor system for delivering synthetic tracers during drilling to assess for possible sample contamination. Here, we describe this new tracer delivery system, assess the performance of the system during the expedition, provide an overview of the quality of the core samples collected for deep biosphere investigations based on tracer concentrations, and make recommendations for future applications of the system.

The great sulfur depletion of Earth’s mantle is not a signature of mantle–core equilibration

The extreme depletion of the Earth’s mantle in sulfur is commonly seen as a signature of metal segregation from Earth’s mantle to Earth’s core. However, in addition to S, the mantle contains other elements as volatile as S that are hardly depleted relative to the lithophile volatility trend although they are potentially as siderophile as sulfur. We report experiments in metal-sulfide–silicate systems to show that the CI normalized abundances of S, Pb, and Sn in Earth’s mantle cannot be reproduced by element partitioning in Fe ± S–silicate systems, neither at low nor at high pressure. Much of the volatile inventory of the Earth’s mantle must have been added late in the accretion history, when metal melt segregation to the core had become largely inactive. The great depletion in S is attributed to the selective segregation of a late sulfide matte from an oxidized and largely crystalline mantle. Apparently, the volatile abundances of Earth’s mantle are not in redox equilibrium with Earth’s core.

Complete synthetic seismograms based on a spherical self-gravitating Earth model with an atmosphere–ocean–mantle–core structure

Complete synthetic seismograms based on a spherical self-gravitating Earth model with an atmosphere–ocean–mantle–core structure

Concentration of Platinum Group Elements during the Early Earth Evolution: A Review*

Numerous unique geological processes [1] took place during the early Earth evolution; several of them, especially those occurring in the Hadean—Early Archean and later, are reflected in the modern geological (geophysical, geochemical, etc.) pattern. One such significant enigmatic feature is the preservation of extremely dense and heavy platinum group elements (PGEs): Pt, Pd, Rh, Ru, Ir, Os. Concentration of PGEs during this period could have taken place in two ways: 1) presence of particular matter capable of preserving PGEs near the earth's surface, 2) transportation of PGEs by magma flows from deep lithospheric (asthenospheric) layers (slabs) to the subsurface. Clearly, much of the dense and heavy PGEs did not sink through to the Earth’s mantle (core) at the time of the magma-ocean, and occur near Earth’s surface in abundances for formation of ore deposits with PGE concentrations found to be 2 - 3 orders of magnitude greater than those in their host media. Their enrichments are associated in numerous cases with such enigmatic phenomena as formation of anorthosites and anorthosite-bearing layered magmatic intrusions. PGE deposits and mineralization zones are also found in associations with chromitites, dunites and serpentinites. In this review, problems related to the initial concentration and preservation of PGEs, their association with anorthosites, and formation of layered intrusions are discussed in detail. The main aim of this article is analysis of the requirements—initial concentration and preservation of PGE and PGM (Platinum Group Minerals) during the early Earth evolution, as well as examination of the distribution behavior of some PGEs in different ore deposits and meteorites. It is supposed that meteoritic bombardment of Earth has played a significant role in formation of PGEs deposits. Some conclusions made in this article may be useful for developing and enhancing strategies of prospecting for PGEs deposits.

HYDROGEN IN THE EARTH’S OUTER CORE, AND ITS ROLE IN THE DEEP EARTH GEODYNAMICS

The content of hydrogen in the outer core of the Earth is roughly quantified from the dependence of the density of iron (viewed as the main component of the core) on the amount of hydrogen dissolved in the core, with account of the most likely presence of iron hydrogen in the outer core, and the matter’s density jumps at the boundaries between the outer liquid core and the internal solid core (that is devoid of hydrogen) and the mantle. Estimations for the outer liquid core show that the hydrogen content varies from 0.67 wt. % at the boundary with the solid inner core to 3.04 wt. % at the boundary with the mantle. Iron occlusion is viewed as the most likely mechanism for the iron–nickel core to capture such a significant amount of hydrogen. Iron occlusion took place at the stage of the young sun when the metallic core emerged in the cooling protoplanetary cloud containing hydrogen in high amounts, and non-volatile hydrogen was accumulated. Absorption (occlusion) of molecular hydrogen was preceded by dissociation of molecules into atoms and ionization of the atoms, as proved by results of studies focused on Fe–H2 system, and hydrogen dissipation was thus prevented. The core matter was subject to gravitational compression at high pressures that contributed to the forced rapprochement of protons and electrons which interaction resulted by the formation of hydrogen atoms. Highly active hydrogen atoms reacted with metals and produced hydrides of iron and nickel, FeH and NiH. While the metallic core and then the silicate mantle were growing and consolidating, the stability of FeH and NiH was maintained due to pressures that were steadily increasing. Later on, due to the impacts of external forces on the Earth, marginal layers at the mantle–core boundary were detached and displaced, pressures decreased in the system, and iron and nickel hydrides were decomposed to produce molecular hydrogen. Consequences of the hydrides transformation into molecular hydrogen are important in terms of petrology, mineralogy and geodynamics. At high temperatures, molecular hydrogen can be involved in redox reactions with iron silicates and carbonaceous gases (CO and CO2), and the synthesis of water becomes possible throughout the entire mantle. It is known that water can significantly reduce the temperature of rock melting, which leads to partial melting of the rocks and pluming in the asthenosphere (in the D” layer) at the bottom of the mantle, and causes the hydrolysis of magnesium silicates, which results in the chemically bound state (hydroxyl ions). Highly ductile hydroxyl-containing magnesium silicates can alter rheological properties of the rocks. A combination of rheologically weak areas in the mantle rocks and the external cosmic effects can cause significant impacts on the tectonic activity and facilitate its manifestation throughout the entire mantle.

Thermodynamic Predictions of Thermal Expansivity and Elastic Compliances at High Temperatures and Pressures Applied to Perovskite Crystals

The possibility of near zero thermal expansion coefficients at very high pressures is explored for application to the Earth’s core materials and mantle dynamics. The pressures in the Earth are large enough to effectively reduce thermal expansion coefficients to values which will decouple heat from mechanical work. It is shown that at pressures below the bulk modulus the thermal expansion coefficient will approach zero in all simple linear-elastic crystalline models. Advanced models of crystalline elastic solids based on interatomic potentials and density functional theory are shown to violate Gibb’s potential for a solid, crystalline material described by three elastic matrix compliance entries; it is established that the temperature dependence of S 11 and S 12 are thermodynamically identical; it is also established that the pressure dependence of S 11 and S 12 are thermodynamically identical. The basis for thermal energy in materials is the phonon energy in solids. However, it is noted that heat capacity measurements which are obtained from constant pressure heat capacity conditions converted to constant volume values on isobars are not in the correct state when compared to theoretical models; at atmospheric pressure there may be very little difference between these states but at very high pressures the effect may be major. Very large pressures always reduce thermal expansion coefficients; the importance of very small thermal expansion coefficients is discussed in relation to physical processes deep in the core and mantle of the Earth.