Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mantle fluids”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

The Edge of Wetness: The Case for Dry Magmatism on Mars

The issue of whether martian magmas are wet or dry is an important one. The answer to this basic question has profound consequences for how we think about Mars as a planet. Recently, several lines of evidence have been presented that collectively suggest that shergottite parent magmas were once wet. These include: (i) phase equilibria studies that indicate that the Shergotty parent magma required ~2 wt.% water in order to be co-saturated with both pigeonite and augite, (ii) reverse zoning of light lithophile elements (Li and B) in shergottite pyroxenes, suggesting the exsolution and removal of an aqueous fluid, and (iii) measurement of D/H ratios in SNC minerals that are much lower than atmospheric, suggesting that there may be juvenile (primordial) mantle water. Below I will review the evidence for the diametrically opposite case, that shergottite magmas were effectively dry (<< 1 wt.% H 2 O).

J.H. Jones↗
NASA NTRS2016

Apatite-Melt Partitioning at 1 Bar: An Assessment of Apatite-Melt Exchange Equilibria Resulting from Non-Ideal Mixing of F and Cl in Apatite

The mineral apatite [Ca5(PO4)3(F,Cl,OH)] is present in a wide range of planetary materials. Due to the presence of volatiles within its crystal structure (X-site), many recent studies have attempted to use apatite to constrain the volatile contents of planetary magmas and mantle sources. In order to use the volatile contents of apatite to precisely determine the abundances of volatiles in coexisting silicate melt or fluids, thermodynamic models for the apatite solid solution and for the apatite components in multi-component silicate melts and fluids are required. Although some thermodynamic models for apatite have been developed, they are incomplete. Furthermore, no mixing model is available for all of the apatite components in silicate melts or fluids, especially for F and Cl components. Several experimental studies have investigated the apatite-melt and apatite-fluid partitioning behavior of F, Cl, and OH in terrestrial and planetary systems, which have determined that apatite-melt partitioning of volatiles are best described as exchange equilibria similar to Fe-Mg partitioning between olivine and silicate melt. However, McCubbin et al. recently reported that the exchange coefficients may vary in portions of apatite compositional space where F, Cl, and OH do not mix ideally in apatite. In particular, solution calorimetry data of apatite compositions along the F-Cl join exhibit substantial excess enthalpies of mixing. In the present study, we conducted apatite-melt partitioning experiments in evacuated, sealed silica-glass tubes at approximately 1 bar and 950-1050 degrees Centigrade on a synthetic Martian basalt composition equivalent to the basaltic shergottite Queen Alexandria Range (QUE) 94201. These experiments were conducted dry, at low pressure, to assess the effects of temperature and apatite composition on the partitioning behavior of F and Cl between apatite and basaltic melt along the F-Cl apatite binary join, where there is non-ideal mixing of F and Cl in apatite.

McCubbin, F. M.↗
NASA NTRS2017

Partitioning of F and Cl Between Apatite and a Synthetic Shergottite Liquid (QUE 94201) at 4 Gpa from 1300 TO 1500 C

Apatite [Ca5(PO4)3(F,Cl,OH)] is present in a wide range of planetary materials. Due to the presence of volatiles within its crystal structure (Xsite), many recent studies have attempted to use apatite to constrain the volatile contents of planetary magmas and mantle sources. In order to use the volatile contents of apatite to accurately determine the abundances of volatiles in coexisting silicate melt or fluids, thermodynamic models for the apatite solid solution and for the apatite components in multicomponent silicate melts and fluids are required. Although some thermodynamic models for apatite have been developed, they are incomplete. Furthermore, no mixing model is available for all of the apatite components in silicate melts or fluids, especially for F and Cl components. Several experimental studies have investigated the apatite-melt and apatite-fluid partitioning behavior of F, Cl, and OH in terrestrial and planetary systems, which have determined that apatite-melt partitioning of volatiles are best described as exchange equilibria similar to Fe-Mg partitioning between olivine and silicate melt. However, McCubbin et al., recently reported that the exchange coefficients vary in portions of apatite compositional space where F, Cl, and OH do not mix ideally in apatite. In particular, solution calorimetry data of apatite compositions along the F-Cl join exhibit substantial excess enthalpies of mixing, and McCubbin et al. reported substantial deviations in the Cl-F exchange Kd along the F-Cl apatite join that could be explained by the preferential incorporation of F into apatite. In the present study, we assess the effect of apatite crystal chemistry on F-Cl exchange equilibria between apatite and melt at 4 GPa over the temperature range of 1300-1500 C. The goal of these experiments is to assess the variation in the Ap-melt Cl-F exchange Kd over a broad range of F:Cl ratios in apatite. The results of these experiments could be used to understand at what composition apatite shifts from a hexagonal unit cell with space group P63/m to a unit cell with monoclinic symmetry within space group P21/b. We anticipate that this transition occurs at >70% chlorapatite based on solution calorimetry data.

McCubbin, F. M.↗
NASA NTRS1999

A View of the Lunar Interior Through Lunar Laser Range Analysis

Laser ranges between observatories on the Earth and retroreflectors on the Moon started in 1969 and continue to the present. Recent range accuracies are 2 cm while earliest ranges are an C, order of magnitude less certain. Four retroreflectors are ranged: three located at the Apollo 11, 14, and 15 sites and one on the Lunakhod 2 rover. Accurate analysis of the range data determines a number of lunar science parameters. The lunar interior variables include a fluid core parameter. The Lunar Laser Ranging effort is reviewed elsewhere. Many parameters are detected through their influence on rotation. Also detected are solid-body tides and accurate selenocentric reflector locations. Determined through the rotation are moment-of-inertia differences, gravitational harmonics, potential Love number, and dissipation effects due to tides and molten core. The rotation of the Moon is not at its minimum energy state; some recently active process has caused free librations. The moment differences contributed to the recent improvement of the Moon's moment of inertia from the Lunar Prospector gravity field. The Love numbers provide bulk elastic properties. Future possibilities for measurement include oblateness of the core-mantle boundary and core moment. A study of dissipation signatures in the rotation determines tidal Q vs. frequency and concludes that the Moon has a molten core. At 1 month the tidal Q is 37 and at 1 yr it is 60. The core radius is < or = 352 km for Fe and < or = 374 km for the Fe-FeS eutectic. The core detection exceeds 3x its uncertainty. The spin of the core is not aligned with the spin of the mantle and torque arises from the velocity difference at the boundary. Yoder's turbulent boundary layer theory is used to compute the radii. The present heat generation from tides and core interaction is minor compared to radiogenic heating. The heating for ancient times is more interesting. Peale and Cassen investigated lunar tidal heating while the lunar orbit expanded due to tides on Earth. Their calculations predate the measurement of Q and should be multiplied by 3.45 to match the lunar-laser-determined Love number and monthly Q. Tidal-heating computations depend on how fast the lunar orbit evolved and whether the tidal dissipation is localized. Neither is known, but plausible assumptions lead to early central region temperature increases of several hundred degrees. Most of the energy is deposited early in the Moon's history. The turbulent boundary layer theory allows a prediction of energy dissipated at the core-mantle boundary during orbit evolution. Under the assumption that the properties of the early core are the same as at present, the energy dissipated by core-mantle interaction is about the same as for tidal dissipation, but it is deposited in a smaller volume. This source of energy is capable of promoting convection in an early fluid core and driving a dynamo. This is a transient phase with a duration depending on the rate of orbit evolution. Plausible assumptions lead to a duration of a few hundred million years. Thus, the remnant magnetization of many lunar rocks is compatible with a brief global magnetic field powered by dynamical energy dissipation. Analysis of the lunar laser ranges is providing information on lunar geophysics. Future data will improve accuracies of present solution parameters, and several more interior effects should be detectable.

Williams, J. G.↗
NASA NTRS1990

On the consequences of strong stable stratification at the top of earth's outer core

The consequences of strong stable stratification at the top of the earth's fluid outer core are considered, concentrating on the generation of the geomagnetic secular variation. It is assumed that the core near the core-mantle boundary is both strongly stably stratified and free of Lorentz forces: it is found that this set of assumptions severely limits the class of possible motions, none of which is compatible with the geomagnetic secular variation. Relaxing either assumption is adequate: tangentially geostrophic flows are consistent with the secular variation if the assumption that the core is strongly stably stratified is relaxed (while retaining the assumption that Lorentz forces are negligible); purely toroidal flows may explain the secular variation if Lorentz forces are included.

Bloxham, Jeremy↗
NASA NTRS1999

Mass Transport in Global Geophysical Fluids

Mass transports occurring in the atmosphere-hydrosphere-solid Earth-core system (the "global geophysical fluids") are important geophysical phenomena. They occur on all temporal and spatial scales. Examples include air mass and ocean circulations, tides, hydrological water redistribution, mantle processes such as post-glacial rebound, earthquakes and tectonic motions, and core geodynamo activities. With only a few exceptions on the Earth surface, the temporal history and spatial pattern of such mass transport are often not amenable to direct observations. Space geodesy techniques, however, have the capability of monitoring certain direct consequences of the mass transport, including Earth's rotation variations, gravitational field variations, and the geocenter motion. These techniques include the very-long-baseline interferometry, satellite laser ranging and Doppler tracking, and the Global Positioning System, all entail global observational networks. While considerable advances have been made in observing and understanding of the dynamics of Earth's rotation, only the lowest-degree gravitational variations have been observed and limited knowledge of geocenter motion obtained. New space missions, projects and initiatives promise to further improve the measurements and hence our knowledge about the global mass transports. The latter contributes to our understanding and modeling capability of the geophysical processes that produce and regulate the mass transports, as well as the solid Earth's response to such changes in constraining the modeling of Earth's mechanical properties.

Chao, B. F.↗
NASA NTRS2012

Electrical Conductivity of Rocks and Dominant Charge Carriers: Thermally Activated Positive Holes - Part 1

The prevailing view in the geophysics community is that the electrical conductivity structure of the Earth's continental crust over the 5-35 km depth range can best be understood by assuming the presence of intergranular fluids and/or of intragranular carbon films. Based on single crystal studies of melt-grown MgO, magma-derived sanidine and anorthosite feldspars and upper mantle olivine, we present evidence for the presence of electronic charge carriers, which derive from peroxy defects that are introduced during cooling, under non-equilibrium conditions, through a redox conversion of pairs of solute hydroxyl arising from dissolution of H2O.The peroxy defects become thermally activated in a 2-step process, leading to the release of defect electrons in the oxygen anion sublattice. Known as positive holes and symbolized by h(dot), these electronic charge carriers are highly mobile. Chemically equivalent to O(-) in a matrix of O(2-) they are highly oxidizing. Being metastable they can exist in the matrix of minerals, which crystallized in highly reduced environments. The h(dot) are highly mobile. They appear to control the electrical conductivity of crustal rocks in much of the 5-35 km depth range.

Freund, Friedemann T.↗
NASA NTRS2005

Motion of the Mantle in the Translational Modes of the Earth and Mercury

Slichter modes refer to the translational motion of the inner core with respect to the outer core and the mantle. The polar Slichter mode is the motion of the inner core along the axis of rotation. Busse presented an analysis of the polar mode which yielded an expression for its period. Busse's analysis included the assumption that the mantle was stationary. This approximation is valid for planets with small inner cores, such as the Earth whose inner core is about 1/60 of the total planet mass. On the other hand, many believe that Mercury's core may be enormous. If so, the motion of the mantle should be expected to produce a significant effect. We present a formal framework for including the motion of the mantle in the analysis of the translational motion of the inner core. We analyze the effect of the motion of the mantle on the Slichter modes for a non-rotating planet with an inner core of arbitrary size. We omit the effects of viscosity in the outer core, magnetic effects, and solid tides. Our approach is perturbative and is based on a linearization of Euler's equations for the motion of the fluid and Newton's second law for the motion of the inner core. We find an analytical expression for the period of the Slichter mode. Our result agrees with Busse's in the limiting case of small inner core. We present the unexpected result that even for Mercury the motion of the mantle does not significantly change the period of oscillation.

Grinfeld, Pavel↗
NASA NTRS2017

Water and Metasomatism in the Slave Cratonic Lithosphere (Canada): An FTIR Study

Water in the mantle influences melting, viscosity, seismic velocity, and electrical conductivity. The role played by water in the long-term stabilization of cratonic roots is currently being debated. This study focuses on water contents of mantle minerals (olivine, pyroxene and garnet) from xenoliths found in kimberlites of the Archean Slave craton. 19 mantle xenoliths from central Lac de Gras, and 10 from northern Jericho were analyzed by FTIR for water, and their equilibration depths span the several compositional layers identified beneath the region. At both locations, the shallow peridotites have lower water contents in their olivines (11-30 ppm H2O) than those from the deeper layers (28-300 ppm H2O). The driest olivines, however, are not at the base of the cratonic lithosphere (>6 GPa) as in the Kaapvaal craton. Instead, the deepest olivines are hydrous (31-72 ppm H2O at Lac de Gras and 275 ppm H2O at Jericho). Correlations of water in clinopyroxene and garnet with their other trace element contents are consistent with water being added by metasomatism by melts resembling kimberlite precursors in the mantle approx.0.35 Ga ago beneath Lac de Gras. The northern Jericho xenoliths are derived from a region of the Slave craton that is even more chemically stratified, and was affected at depth by the 1.27 Ga Mackenzie igneous events. Metasomatism at Jericho may be responsible for the particularly high olivine water contents (up to 300 ppm H2O) compared to those at Lac de Gras, which will be investigated by acquiring trace-element data on these xenoliths. These data indicate that several episodes of metasomatic rehydration occurred in the deep part of the Slave craton mantle lithosphere, with the process being more intense in the northern part beneath Jericho, likely related to a translithospheric suture serving as a channel to introduce fluids and/or melts in the northern region. Consequently, rehydration of the lithosphere does not necessarily cause cratonic root delamination and these peridotites may represent localized metasomatic zones - the wall rocks to kimberlite magma passage.

Kilgore, McKensie↗
NASA NTRS1997

Thermal History and Mantle Dynamics of Venus

One objective of this research proposal is to develop a 3-D thermal history model for Venus. The basis of our study is a finite-element computer model to simulate thermal convection of fluids with highly temperature- and pressure-dependent viscosities in a three-dimensional spherical shell. A three-dimensional model for thermal history studies is necessary for the following reasons. To study planetary thermal evolution, one needs to consider global heat budgets of a planet throughout its evolution history. Hence, three-dimensional models are necessary. This is in contrasts to studies of some local phenomena or local structures where models of lower dimensions may be sufficient. There are different approaches to treat three-dimensional thermal convection problems. Each approach has its own advantages and disadvantages. Therefore, the choice of the various approaches is subjective and dependent on the problem addressed. In our case, we are interested in the effects of viscosities that are highly temperature dependent and that their magnitudes within the computing domain can vary over many orders of magnitude. In order to resolve the rapid change of viscosities, small grid spacings are often necessary. To optimize the amount of computing, variable grids become desirable. Thus, the finite-element numerical approach is chosen for its ability to place grid elements of different sizes over the complete computational domain. For this research proposal, we did not start from scratch and develop the finite element codes from the beginning. Instead, we adopted a finite-element model developed by Baumgardner, a collaborator of this research proposal, for three-dimensional thermal convection with constant viscosity. Over the duration supported by this research proposal, a significant amount of advancements have been accomplished.

Hsui, Albert T.↗
NASA NTRS1997

Planetary Geophysics and Tectonics

Research supported by grant NAGW-1928 has addressed a variety of problems related to planetary evolution. One important focus has been on questions related to the role of chemical buoyancy in planetary evolution with application to both Venus and the Moon. We have developed a model for the evolution of the Moon (Hess and Parmentier, 1995) in which dense, highly radioactive, late stage magma ocean cumulates sink forming a core. This core heats the overlying, chemically layered mantle giving rise to a heated, chemically well-mixed layer that thickens with time. This Mixed layer eventually becomes hot enough and thick enough that its top begins to melt at a pressure low enough that melt is buoyant, thus creating mare basalts from a high pressure source of the correct composition and at an appropriate time in lunar evolution. In work completed during the last year, numerical experiments on convection in a chemically stably stratified fluid layer heated from below have been completed. These results show us how to calculate the evolution of a mixed layer in the Moon, depending on the heat production in the ilmenite- cumulate core and the chemical stratification of the overlying mantle. Chemical stratification of the mantle after its initial differentiation is would trap heat in the deep interior and prevent the rapid rise of plumes with accompanying volcanism. This trapping of heat in the interior can explain the thickness of the lunar lithosphere as a function of time as well as the magmatic evolution. We show that heat transported to the base of the lithosphere at a rate determined by current estimates of radioactivity in the Moon would not satisfy constraints on elastic lithosphere thickness from tectonic feature associated with basin loading. Trapping heat at depth by a chemically stratified mantle may also explain the absence of global compressional features on the surface that previous models predict for an initially hot lunar interior. For Venus, we developed a model in which the chemical buoyancy of crust and a depleted mantle layer stabilizes the lithosphere for long periods of time and provides a mechanism of episodic planetary evolution (Parmentier and Hess, 1992). Continued thickening of a residual depleted mantle layer eventually suppresses pressure release melting and the creation of depleted mantle. Continued cooling then allows the lithosphere to become heavier than the underlying hotter, undepleted mantle. This repeated instability can occur on time scales appropriate for episodic global resurfacing on Venus. We have also examined the role of the gabbro-eclogite phase transformation on crust and lithosphere stability and as a mechanism of crustal recycling in the absence of plate tectonics. Our work thus far concentrates on the scale of instability that would occur due to cooling or crustal thickening associated with horizontal shortening. Whether repeated overturn can explain the evolution of Venus depends in part on whether sufficient heat transfer can occur between overturns and on constraints provided by understanding observed surface features and evolution.

Parmentier, E. M.↗
NASA NTRS1979

MANTLE: A finite element program for the thermal-mechanical analysis of mantle convection. A user's manual with examples

A finite element computer code for the analysis of mantle convection is described. The coupled equations for creeping viscous flow and heat transfer can be solved for either a transient analysis or steady-state analysis. For transient analyses, either a control volume or a control mass approach can be used. Non-Newtonian fluids with viscosities which have thermal and spacial dependencies can be easily incorporated. All material parameters may be written as function statements by the user or simply specified as constants. A wide range of boundary conditions, both for the thermal analysis and the viscous flow analysis can be specified. For steady-state analyses, elastic strain rates can be included. Although this manual was specifically written for users interested in mantle convection, the code is equally well suited for analysis in a number of other areas including metal forming, glacial flows, and creep of rock and soil.

Thompson, E.↗
NASA NTRS2004

Magnetohydrodynamic Convection in the Outer Core and its Geodynamic Consequences

The Earth's fluid outer core is in vigorous convection through much of the Earth's history. In addition to generating and maintaining Earth s time-varying magnetic field (geodynamo), the core convection also generates mass redistribution in the core and a dynamical pressure field on the core-mantle boundary (CMB). All these shall result in various core-mantle interactions, and contribute to surface geodynamic observables. For example, electromagnetic core-mantle coupling arises from finite electrically conducting lower mantle; gravitational interaction occurs between the cores and the heterogeneous mantle; mechanical coupling may also occur when the CMB topography is aspherical. Besides changing the mantle rotation via the coupling torques, the mass-redistribution in the core shall produce a spatial-temporal gravity anomaly. Numerical modeling of the core dynamical processes contributes in several geophysical disciplines. It helps explain the physical causes of surface geodynamic observables via space geodetic techniques and other means, e.g. Earth's rotation variation on decadal time scales, and secular time-variable gravity. Conversely, identification of the sources of the observables can provide additional insights on the dynamics of the fluid core, leading to better constraints on the physics in the numerical modeling. In the past few years, our core dynamics modeling efforts, with respect to our MoSST model, have made significant progress in understanding individual geophysical consequences. However, integrated studies are desirable, not only because of more mature numerical core dynamics models, but also because of inter-correlation among the geophysical phenomena, e.g. mass redistribution in the outer core produces not only time-variable gravity, but also gravitational core-mantle coupling and thus the Earth's rotation variation. They are expected to further facilitate multidisciplinary studies of core dynamics and interactions of the core with other components of the Earth.

Kuang, Weijia↗
NASA NTRS2004

Space Geodesy Monitoring Mass Transport in Global Geophysical Fluids

Mass transports occurring in the atmosphere-hydrosphere-cryosphere-solid Earth-core system (the 'global geophysical fluids') are important geophysical phenomena. They occur on all temporal and spatial scales. Examples include air mass and ocean circulations, oceanic and solid tides, hydrological water and idsnow redistribution, mantle processes such as post-glacial rebound, earthquakes and tectonic motions, and core geodynamo activities. The temporal history and spatial pattern of such mass transport are often not amenable to direct observations. Space geodesy techniques, however, have proven to be an effective tool in monitorihg certain direct consequences of the mass transport, including Earth's rotation variations, gravitational field variations, and the geocenter motion. Considerable advances have been made in recent years in observing and understanding of these geodynamic effects. This paper will use several prominent examples to illustrate the triumphs in research over the past years under a 'Moore's law' in space geodesy. New space missions and projects promise to further advance our knowledge about the global mass transports. The latter contributes to our understanding of the geophysical processes that produce and regulate the mass transports, as well as of the solid Earth's response to such changes in terms of Earth's mechanical properties.

Chao, Benjamin F.↗
NASA NTRS1995

Rheological transition in mantle convection with a composite temperature-dependent, non-Newtonian and Newtonian rheology

Numerical simulations of mantle convection with a composite temperature-dependent, Newtonian and non-Newtonian creep law have revealed a transition in the dominant creep mechanism with the increasing vigour of convection. Newtonian creep is found to dominate in the low Rayleigh number regime. With sufficiently high effective Rayleigh number, the overall creep mechanism in the convective flow becomes non-Newtonian. The transitional Rayleigh number increases strongly with the activation energy. These results would suggest a scenario that in the early epochs of Earth the flow in the mantle would have been governed by non-Newtonian rheology and would have exhibited both strong spatial and temporal fluctuations. With time the flow mechanism would behave like a Newtonian fluid and would have a different time-dependent character. In time-dependent Newtonian-dominated flows there are still localized features with distinctly non-Newtonian character. Our analysis of the relative contributions to the lateral viscosity field supports the idea that the inference of the nature of lateral viscosity heterogeneities by seismic tomography may be strongly contaminated by the dominant non-Newtonian contributions to the total lateral viscosity field.

Van Den Berg, Arie P.↗
NASA NTRS2014

Dynamic Responses of the Earth's Outer Core to Assimilation of Observed Geomagnetic Secular Variation

Assimilation of surface geomagnetic observations and geodynamo models has advanced very quickly in recent years. However, compared to advanced data assimilation systems in meteorology, geomagnetic data assimilation (GDAS) is still in an early stage. Among many challenges ranging from data to models is the disparity between the short observation records and the long time scales of the core dynamics. To better utilize available observational information, we have made an effort in this study to directly assimilate the Gauss coefficients of both the core field and its secular variation (SV) obtained via global geomagnetic field modeling, aiming at understanding the dynamical responses of the core fluid to these additional observational constraints. Our studies show that the SV assimilation helps significantly to shorten the dynamo model spin-up process. The flow beneath the core-mantle boundary (CMB) responds significantly to the observed field and its SV. The strongest responses occur in the relatively small scale flow (of the degrees L is approx. 30 in spherical harmonic expansions). This part of the flow includes the axisymmetric toroidal flow (of order m = 0) and non-axisymmetric poloidal flow with m (is) greater than 5. These responses can be used to better understand the core flow and, in particular, to improve accuracies of predicting geomagnetic variability in future.

secular variation↗
NASA NTRS1988

Driving forces: Slab subduction and mantle convection

Mantle convection is the mechanism ultimately responsible for most geological activity at Earth's surface. To zeroth order, the lithosphere is the cold outer thermal boundary layer of the convecting mantle. Subduction of cold dense lithosphere provides tha major source of negative buoyancy driving mantle convection and, hence, surface tectonics. There are, however, importnat differences between plate tectonics and the more familiar convecting systems observed in the laboratory. Most important, the temperature dependence of the effective viscosity of mantle rocks makes the thermal boundary layer mechanically strong, leading to nearly rigid plates. This strength stabilizes the cold boundary layer against small amplitude perturbations and allows it to store substantial gravitational potential energy. Paradoxically, through going faults at subduction zones make the lithosphere there locally weak, allowing rapid convergence, unlike what is observed in laboratory experiments using fluids with temperature dependent viscosities. This bimodal strength distribution of the lithosphere distinguishes plate tectonics from simple convection experiments. In addition, Earth has a buoyant, relatively weak layer (the crust) occupying the upper part of the thermal boundary layer. Phase changes lead to extra sources of heat and bouyancy. These phenomena lead to observed richness of behavior of the plate tectonic style of mantle convection.

Hager, Bradford H.↗