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 163 records · Page 9

Does Mercury have a molten core

A feasible nonseismic observational experiment is proposed for determining the existence and extent of a conducting molten core within Mercury. This experiment would utilize the effects of a liquid core on the dynamics of Mercury's rotation; two necessary conditions for performing it are that the core must not follow the mantle's forced librations in longitude but must follow the mantle on the timescale of the 250,000-yr precession. A method is developed by assuming these conditions to be satisfied, and bounds are established on the core viscosity for which they are satisfied. It is shown that the value of the ratio of the moment of inertia of the mantle to the largest principal moment of inertia of the entire planet would indicate whether the core is most probably solid, partially fluid, or entirely fluid. Techniques are suggested for determining the unknowns required to compute the necessary ratio.

Peale, S. J.↗

The Tethered Moon

A reasonable initial condition on Earth after the Moonforming impact is that it begins as a hot global magma ocean1,2. We therefore begin our study with the mantle as a liquid ocean with a surface temperature on the order of 3000- 4000 K at a time some 100-1000 years after the impact, by which point we can hope that early transients have settled down. A 2nd initial condition is a substantial atmosphere, 100-1000 bars of H2O and CO2, supplemented by smaller amounts of CO, H2, N2, various sulfur-containing gases, and a suite of geochemical volatiles evaporated from the magma. Third, we start the Moon with its current mass at the relevant Roche limit. The 4th initial condition is the angular momentum of the Earth-Moon system. Canonical models hold this constant, whilst some recent models begin with considerably more angular momentum than is present today. Here we present a ruthlessly simplified model of Earth's cooling magmasphere based on a full-featured atmosphere and including tidal heating by the newborn Moon. Thermal blanketing by H2O-CO2 atmospheres slows cooling of a magma ocean. Geochemical volatiles - chiefly S, Na, and Cl - raise the opacity of the magma ocean's atmosphere and slow cooling still more. We assume a uniform mantle with a single internal (potential) temperature and a global viscosity. The important "freezing point" is the sharp rheological transition between a fluid carrying suspended crystals and a solid matrix through which fluids percolate. Most tidal heating takes place at this "freezing point" in a gel that is both pliable and viscous. Parameterized convection links the cooling rate to the temperature and heat generation inside the Earth. Tidal heating is a major effect. Tidal dissipation in the magma ocean is described by viscosity. The Moon is entwined with Earth by the negative feedback between thermal blanketing and tidal heating that comes from the temperature-dependent viscosity of the magma ocean. Because of this feedback, the rate that the Moon's orbit evolves is limited by the modest radiative cooling rate of Earth's atmosphere, which in effect tethers the Moon to the Earth. Consequently the Moon's orbit evolves orders of magnitude more slowly than in conventional models. Slow orbital evolution promotes capture by orbital resonances that may have been important in the Earth-Moon system

Earth-Moon system↗

Volatile solubities in magmas - Transport of volatiles from mantles to planet surfaces

The solubility of CO2, CO, CH4, H2, and H2O in melts of NaAlSi3O8 composition was measured at the mantle conditions of pressure (10-20 Kb) and temperature (1200 C). The melt volatiles were found to have dramatically lower H2/H2O and higher CO/CO2 ratios than the fluid, due to large variations in relative solubilities of volatile species in the melt. Partial melting of source regions of either high or low oxygen fugacity (fO2) will result in magmas with intermediate fO2. This means that volcanic gases can be either more, or less reducing than their source regions and that a volcanic gas composition cannot be used to directly estimate either the fO2 or the volatile composition of the source region. The results suggest that volcanic gases will usually lie in the 'neutral' range, with the fO2 values near those of the quartz-fayalite-magnetite buffer. These gases are predominately H2O with minor CO2, CO, CH4, and H2. These conclusions should apply to earth, Mars, and Venus mantles, in which the magmas produced by partial melting have moderate silica contents, but will probably not apply to very low silica magmas, such as kimberlites, because of the very high solubility of CO2 in those magmas.

Holloway, J. R.↗

Surface topography due to convection in a variable viscosity fluid - Application to short wavelength gravity anomalies in the central Pacific Ocean

Finite difference calculations of thermal convection in a fluid layer with a viscosity exponentially decreasing with temperature are performed in the context of examining the topography and gravity anomalies due to mantle convection. The surface topography and gravity anomalies are shown to be positive over regions of ascending flow and negative over regions of descending flow; at large Rayleigh numbers the amplitude of surface topography is inferred to depend on Rayleigh number to the power of 7/9. Compositional stratifications of the mantle is proposed as a mechanism for confining small-scale convection to a thin layer. A comparative analysis of the results with other available models is included.

Lin, J.↗

The pole tide in deep oceans

The fluid-dynamical theory of the pole tide is examined by describing the oceanic response to the Chandler wobble and assessing its implications for mantle anelasticity and low-frequency ocean dynamics. The Laplace tide equations accounting for bottom friction are given, and a spherical harmonic approach is delineated in which the time-independent portion of the tide height is expanded. Pole-tide height and related inertia products are linearly proportional to wobble amplitude, and the final equations are modified to account for mantle elasticity and oceanic loading. Results for pole tide effects are given for various earth models with attention to the role of boundary constraints. A dynamic effect is identified which lengthens the Chandler period by about 1 day more than static lengthening, a contribution that suggests a vigorous low-frequency response. The values derived are shown to agree with previous models that do not incorporate the effects of the pole tide.

Dickman, S. R.↗

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↗

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.↗

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.↗

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.↗

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↗

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.↗

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.↗

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↗

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↗

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.↗

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.↗

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.↗