SEARCH · Engineering Papers
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.
The Earth is a dynamic system-it has a fluid, mobile atmosphere and oceans, a continually changing distribution of ice, snow, and groundwater, a fluid core undergoing hydromagnetic motion, a mantle undergoing both thermal convection and rebound from glacial loading of the last ice age, and mobile tectonic plates.
Subduction zones are presently the dominant sites on Earth for recycling and mass transfer between the crust and mantle; they feed hydrated basaltic oceanic crust into the upper mantle, where dehydration reactions release aqueous fluids and/or hydrous melts. The loci for fluid and/or melt generation will be determined by the intersection of dehydration reaction boundaries of primary hydrous minerals within the subducted lithosphere with slab geotherms. For metabasalt of the oceanic crust, amphibole is the dominant hydrous mineral. The dehydration melting solidus, vapor-absent melting phase relationships; and amphibole-out phase boundary for a number of natural metabasalts have been determined experimentally, and the pressure-temperature conditions of each of these appear to be dependent on bulk composition. Whether or not the dehydration of amphibole is a fluid-generating or partial melting reaction depends on a number of factors specific to a given subduction zone, such as age and thickness of the subducting oceanic lithosphere, the rate of convergence, and the maturity of the subduction zone. In general, subduction of young, hot oceanic lithosphere will result in partial melting of metabasalt of the oceanic crust within the garnet stability field; these melts are characteristically high-Al2O3 trondhjemites, tonalites and dacites. The presence of residual garnet during partial melting imparts a distinctive trace element signature (e.g., high La/Yb, high Sr/Y and Cr/Y combined with low Cr and Y contents relative to demonstrably mantle-derived arc magmas). Water in eclogitized, subducted basalt of the oceanic crust is therefore strongly partitioned into melts generated below about 3.5 GPa in 'hot' subduction zones. Although phase equilibria experiments relevant to 'cold' subduction of hydrated natural basalts are underway in a number of high-pressure laboratories, little is known with respect to the stability of more exotic hydrous minerals (e.g., ellenbergite) and the potential for oceanic crust (including metasediments) to transport water deeper into the mantle.
Abstract The physical processes during planet formation span a large range of pressures and temperatures. Giant impacts, such as the one that formed the Moon, achieve peak pressures of 100s of GPa. The peak shock states generate sufficient entropy such that subsequent decompression to low pressures intersects the liquid‐vapor phase boundary. The entire shock‐and‐release thermodynamic path must be calculated accurately in order to predict the post‐impact structures of planetary bodies. Forsterite (Mg 2 SiO 4 ) is a commonly used mineral to represent the mantles of differentiated bodies in hydrocode models of planetary collisions. Here, we performed shock experiments on the Sandia Z Machine to obtain the density and temperature of the liquid branch of the liquid‐vapor phase boundary of forsterite. This work is combined with previous work constraining pressure, density, temperature, and entropy of the forsterite principal Hugoniot. We find that the vapor curves in previous forsterite equation of state models used in giant impacts vary substantially from our experimental results, and we compare our results to a recently updated equation of state. We have also found that due to under‐predicted entropy production on the principal Hugoniot and elevated temperatures of the liquid vapor phase boundary of these past models, past impact studies may have underestimated vapor production. Furthermore, our results provide experimental support to the idea that giant impacts can transform much of the mantles of rocky planets into supercritical fluids.
Summary Serpentinization is the hydration and oxidation of ultramafic rock, which occurs as oceanic lithosphere is emplaced onto continental margins (ophiolites), and along the seafloor as faulting exposes this mantle‐derived material to circulating hydrothermal fluids. This process leads to distinctive fluid chemistries as molecular hydrogen (H 2 ) and hydroxyl ions (OH − ) are produced and reduced carbon compounds are mobilized. Serpentinizing ophiolites also serve as a vector to transport sulfur compounds from the seafloor onto the continents. We investigated hyperalkaline, sulfur‐rich, brackish groundwater in a serpentinizing continental ophiolite to elucidate the role of sulfur compounds in fuelling in situ microbial activities. Here we illustrate that key sulfur‐cycling taxa, including Dethiobacter , Desulfitispora and ‘ Desulforudis ’, persist throughout this extreme environment. Biologically catalysed redox reactions involving sulfate, sulfide and intermediate sulfur compounds are thermodynamically favourable in the groundwater, which indicates they may be vital to sustaining life in these characteristically oxidant‐ and energy‐limited systems. Furthermore, metagenomic and metatranscriptomic analyses reveal a complex network involving sulfate reduction, sulfide oxidation and thiosulfate reactions. Our findings highlight the importance of the complete inorganic sulfur cycle in serpentinizing fluids and suggest sulfur biogeochemistry provides a key link between terrestrial serpentinizing ecosystems and their submarine heritage.
Time-dependent numerical models of two-dimensional convection which simulate mantle convection are presented. Starting from a fluid initially at rest with a purely conductive temperature profile, it is found that for a Rayleigh number of ten million, the model experiences a transient period of double-layer convection which lasts on the order of hundreds of millions of years when scaled to the earth's mantle. It is suggested that transient periods of double-layer convection in the earth's mantle may have existed in the past, whether or not such a period exists today.
Hotspots are associated with long wavelength geoid highs, an association that is even stronger when the geoid highs associated with subduction zones are removed. These associations are quantified by expanding the hotspot distribution in spherical harmonics and calculating correlation coefficients as a function of harmonic degree. The hotspot distribution spectrum is essentially white, with peaks at degrees 2 and 6. It is correlated positively with the slab residual geoid for degrees 2 to 6, with low seismic velocity in the lower mantle at degree 2, and with low seismic velocity in the upper mantle at degree 6. A variety of fluid mechanical models were tested for hotspots, including lithospheric delamination and hot plumes, by calculating their predicted dynamic geoid responses and comparing them to the observations. These models include the effects of temperature dependent rheology. The preferred hotspot model, based on observations of the geoid and seismic tomography, has plumes preferentially occurring in regions of large scale background temperature highs in a mantle with substantial viscosity increase with depth, although other models are possible.
Hotspots are associated with long wavelength geoid highs, an association that is even stronger when the geoid highs associated with subduction zones are removed. These associations are quantified by expanding the hotspot distribution in spherical harmonics and calculating correlation coefficients as a function of harmonic degree. The hotspot distribution spectrum is essentially white, with peaks at degrees 2 and 6. It is correlated positively with the slab residual geoid for degrees 2 to 6, with low seismic velocity in the lower mantle at degree 2, and with low seismic velocity in the upper mantle at degree 6. A variety of fluid mechanical models were tested for hotspots, including lithospheric delamination and hot plumes, by calculating their predicted dynamic geoid responses and comparing them to the observations. These models include the effects of temperature dependent rheology. The preferred hotspot model, based on observations of the geoid and seismic tomography, has plumes preferentially occurring in regions of large scale background temperature highs in a mantle with substantial viscosity increase with depth, although other models are possible.
An alternative method for determining the fluid motion immediately beneath the core-mantle boundary is presented which is based on solving the full nonlinear core motions problem. This method is used to examine three dynamical hypotheses about the flow: (1) the steady motions hypothesis; (2) the geostrophic hypothesis; and (3) the toroidal flow hypothesis. Better fits to the field are obtained with the toroidal flows than with geostrophic flows, casting considerable doubt on the validity of the geostrophic hypothesis. Additionally, some indication is found that failure of the frozen-flux approximation, a concomitant assumption, may be a serious obstacle to obtaining reliable maps of the core fluid flow.
Select models of the main geomagnetic field and its secular variation are extrapolated to the base of an insulating mantle and used to estimate the adjacent fluid motion of a perfectly conducting outer core. The assumption of steady motion provides formally unique solutions and is tested along with that of no upwelling. The hypothesis of no upwelling is found to be substantially worse than that of steady motion. Although the actual motion is not thought to be steady, the large-scale secular variation at the top of the core can be adequately described by a large-scale, combined toroidal-poloidal circulation which is steady for intervals of at least a decade or two. The derived flows include a bulk westward drift but are complicated by superimposed jets, gyres, and surface divergence indicative of vigorous vertical motion at depth. The circulation pattern and key global properties including rms speed, upwelling, and westward drift are found to be fairly insensitive to variations in modeling parameters.
Direct numerical simulations of two-dimensional high Rayleigh (Ra) number, base-heated thermal convection in large aspect-ratio boxes are presented for infinite Prandtl number fluids, as applied to the earth's mantle. A transition is characaterized in the flow structures in the neighborhood of Ra between 10 to the 7th and 10 to the 8th. These high Ra flows consist of large-scale cells with strong intermittent, boundary-layer instabilities. For Ra exceeding 10 to the 7th it is found that the heat-transfer mechanism changes from one characterized by mushroom-like plumes to one consisting of disconnected ascending instabilities, which do not carry with them all the thermal anomaly from the bottom boundary layer. Plume-plume collisions become much more prominent in high Ra situations and have a tendency of generating a pulse-like behavior in the fixed plume. This type of instability represents a distinct mode of heat transfer in the hard turbulent regime. Predictions of this model can be used to address certain issues concerning the mode of time-dependent convection in the earth's mantle.
The main findings of a study of earth rotation changes on the longer interannual and decade range of time scales are summarized. The earth rotation data provide important constraints on models of the internal structure and dynamics of the core and lower mantle. Torques at the core-mantle interface are due to time varying fluid motions in the core. A torque of about 0.2 Hadley units is found in the axial component during the 1980 period.
A numerical investigation is conducted of the effect of compressibility on 3D thermal convection in a basally heated highly viscous fluid spherical shell (representing the earth's mantle). It was found that compressibility has little effect on the spatial structure of steady convection when the superadiabatic temperature drop across the shell, Delta T(sa), is comparable to a characteristic adiabatic temperature. When the Delta T(sa) is about an order of magnitude smaller than the adiabatic temperature, compressibility effect is significant. It is shown that compressibility affects the pattern evolution of the irregular solutions; fewer upwelling plumes are produced with increasing compressibility.
Pressure load Love numbers are presented for the mantle deformation induced by the variation of the pressure field at the core mantle boundary (CNB). We find that the CMB geostrophic pressure fields, derived from 'frozen-flux' core surface flow estimates at epochs 1965 and 1975, produce a relative radial velocity (RRV) field in the range of 3mm/decade with uplift near the equator and subsidence near the poles. The contribution of this mechanism to the change in the length of day (l.o.d) is small --- about 2.3 x 10(exp -2) ms/decade. The contribution to the time variation of the ellipticity coefficient is more important --- -1.3 x 10(exp -11)/yr.
Calcium isotope ratios are readily mass-fractionated by low-temperature depositional processes that can be recorded in sedimentary rocks and, therefore, have the potential to track distinct geochemical signatures of recycled materials in mantle-derived igneous rocks. In this study, I report calcium isotopic compositions for well-characterized lavas from the Central American volcanic arc that exhibit a range of trace element and radiogenic isotope signatures that have been used to indicate variable amounts of subducted marine carbonate and hemipelagic sedimentary rocks, fresh carbonatite lavas from Oldoinyo Lengai Volcano, Tanzania, and an intrusive carbonatite clast erupted 12.8 ka from the Laacher See Volcano, Germany. I observed no calcium isotope evidence for recycled sedimentary rocks in the Central American arc magmas or in the Oldoinyo Lengai carbonatite volcanic rocks. They all exhibit isotopic compositions similar to rocks dominated by a primitive, mantle-like, bulk silicate Earth (BSE) composition. The exception in this work is the calcium isotope composition measured in the intrusive Laacher See carbonatite that is resolvable from BSE (δ^44 Ca = 0.0). Although this calcium isotopic signature might be related to ancient carbonate recycling, the magnitude and relatively light calcium isotopic composition (δ^44 Ca = –0.4) of this carbonatite could also reflect an origin that involved partial melting of altered lithospheric mantle, from which isotopically heavy Ca-bearing fluids have been lost. The decoupled signatures between trace elements and their radiogenic isotopes and those recorded by calcium isotope data are important because they likely reflect different sources and processes, and demonstrate that sediment subduction is not a bulk mixing process.
It was established that the total absolute magnetic flux crossing the core- mantle boundary has been a constant of the core motion for the last 50 years. This provides a scalar constraint that could be added to the geometric modelling procedure. The GSFC 8 8/80 model is being evaluated. The absolute magnetic flux linking the CMB to that model was plotted as a function of time during the span covered by the data, and increasing truncation level. The inclusion of the standard error of each Gauss coefficient derived from the statistics of fit in the GSFC 9/80 model is useful. The magnitude and sense (upwelling or downe welling) of vertical fluid motion adjacent to the core-mantle boundary was calculated using the model. Standard errors were found to be sufficiently small at all but one or two of the 40 or more critical points of B sub r. They do not nearly overlap the value gamma u/gamma r = 0. It is concluded that the core is upwelling and downwelling at an observationally detectable level.
We report on the mineralogy, petrography, and O and Al-Mg isotopic systematics of secondary mineralization in the metasomatically altered igneous Ca,Al-rich inclusions (CAIs) [compact type A (CTA), B1, B2, forsterite-bearing B (FoB), and C] from the CV3 carbonaceous chondrite Allende. This alteration affected mainly melilite, and to a lesser degree anorthite, and resulted in the formation of a variety of secondary minerals, including adrianite, Al-diopside, andradite, anorthite, calcite, celsian, clintonite, corundum, dmisteinbergite, ferroan olivine, ferroan monticellite, ferroan Al-diopside, forsterite, grossular, heazlewoodite, hedenbergite, hutcheonite, kushiroite, margarite, monticellite, Na-melilite, nepheline, pentlandite, pyrrhotite, sodalite, spinel, tilleyite, wadalite, and wollastonite. The secondary mineral assemblages are mainly defined by chemical compositions of the primary melilite replaced and elements introduced by an aqueous fluid. Gehlenitic melilite (Åk <35 ) in CTAs and mantles of B1s is mainly replaced by anorthite + grossular; clintonite, corundum, spinel, and Al-diopside are minor. Åkermanitic melilite (Åk 35-90 ) in type B2s, FoBs, and cores of B1s is replaced by the grossular + monticellite + wollastonite, grossular + monticellite, and grossular + Al-diopside assemblages; forsterite, spinel, clintonite, and Na-melilite are minor. In type Cs, lacy melilite (åkermanitic melilite with rounded inclusions of anorthite) is pseudomorphically replaced by the grossular + forsterite + monticellite and grossular + Al-diopside assemblages; Na-melilite is minor. Primary and secondary anorthites in the peripheral portions of CAIs are replaced by nepheline, sodalite, and ferromagnesian olivine. Some CAIs contain voids and cracks filled by andradite, hedenbergite, wollastonite, ±sodalite, ±grossular, ±monticellite, ±tilleyite, and ±calcite. All CAIs studied are surrounded by Wark-Lovering rims, fine-grained matrix-like rims composed of lath-shaped ferroan olivine and abundant nepheline grains, and a layer of salite-hedenbergite pyroxenes + andradite + wollastonite. Grossular associating with monticellite, Al-diopside, and forsterite and replacing åkermanitic melilite ( 27 Al/ 24 Mg ~ 2) has high 27 Al/ 24 Mg ratios (30–100) and shows no resolvable excess of radiogenic 26 Mg (26Mg*). The 27 Al/ 24 Mg ratios (7–10) and 26 Mg* (2–3‰) in the nearly monomineralic grossular veins crosscutting gehlenitic melilite are similar to those of the host melilite and plot along a regression line with 26 Al/ 27 Al ratio of ~5×10 –5 . Oxygen isotopic compositions of secondary minerals in the most Type Bs measured in situ with the UH Cameca ims-1280 and matrix-matched standards plot along mass-dependent fractionation line with Δ 17 O of ~ –3±2‰ with δ 18 O ranging from ~0 to ~10‰. Primary melilite and anorthite in the host CAIs are similarly 16 O-depleted, whereas spinel, forsterite, and most Al,Ti-diopside grains have 16 O-rich compositions (Δ 17 O ~ –25±2‰). Secondary grossular and forsterite in type Cs and type B1 CAI TS-34 show a range of Δ 17 O, from ~ –15 to ~ –1‰; the 16 O-enriched compositions of grossular and forsterite plot along the carbonaceous chondrite anhydrous mineral line. The similar ranges of Δ 17 O and positions on the three-isotope oxygen diagram are observed for primary anorthite; melilite is generally 16 O-depleted compared to anorthite (Δ 17 O ~ –5 to –1±2‰); spinel and fassaite are 16 O-rich (except very Ti-rich fassaite in TS-34 and CTA CAIs). We conclude that Allende CAIs experienced an open-system in situ metasomatic alteration at relatively high temperatures (200-250 °C) in the presence of CO 2 - and H 2 O-bearing fluid with Δ 17 O of ~ –3±2‰ followed by thermal metamorphism at ~ 500 °C on the CV chondrite parent asteroid. During the alteration, most elements were mobile: Si, Na, Cl, K, Fe, S, and Ni were introduced; Al, Ti, Mg, and Ba were locally mobilized; Ca and some Mg and Al were lost from the host inclusions. The alteration occurred after nearly complete decay of 26 Al, >3 Ma after crystallization of CAIs with the canonical ( 26 Al/ 27 Al) 0 of (5.25±0.02)×10 -5 ; 26 Mg* in grossular was inherited from the primary melilite and provide no chronological significance. Oxygen isotopic heterogeneity of primary minerals in the Allende CAIs at least partly is due to isotopic exchange with an aqueous fluid that largely affected melilite, anorthite, perovskite, Zr- and Sc-rich oxides and silicates, and possibly very Ti-rich fassaite.
The surface of Ganymede, as seen in Voyager images, was greatly modified by tectonic activity. Geologic evidence is consistent with the formation of younger bright terrain by the emplacement of icy material into rift zones formed in the older dark terrain. To better understand the formation of both terrestrial rift zones and bright terrain bands on Ganymede, the mechanics of rift initiation by the growth of necking instabilities in an extending lithosphere are studied. In particular, how an initial, small perturbation in lithosphere thickness at the lithosphere/asthenosphere boundary is amplified by uniform horizontal extension is studied. The lithosphere is treated as a layer of power-law viscous fluid of thickness H overlying a weaker, viscous mantle substrate. In order to assess the relative importance of buoyant upwelling and mechanical instability in rift formation, both density and strength stratification were incorporated into the model.
The work done under NASA grant NAG5-485 included modelling the deformation of the earth caused by variations in atmospheric pressure. The amount of deformation near coasts is sensitive to the nature of the oceanic response to the pressure. The PSMSL (Permanent Service for Mean Sea Level) data suggest the response is inverted barometer at periods greater than a couple months. Green's functions were constructed to describe the perturbation of the geoid caused by atmospheric and oceanic loading and by the accompanying load-induced deformation. It was found that perturbation of up to 2 cm are possible. Ice mass balance data was used for continental glaciers to look at the glacial contributions to time-dependent changes in polar motion, the lod, the earth's gravitational field, the position of the earth's center-of-mass, and global sea level. It was found that there can be lateral, non-hydrostatic structure inside the fluid core caused by gravitational forcing from the mantle, from the inner core, or from topography at the core/mantle or inner core/outer core boundaries. The nutational and tidal response of a non-hydrostatic earth with a solid inner core was modeled. Monthly, global tide gauge data from PSMSL was used to look at the 18.6-year ocean tide, the 14-month pole tide, the oceanic response to pressure, the linear trend and inter-annual variability in the earth's gravity field, the global sea level rise, and the effects of post glacial rebound. The effects of mantle anelasticity on nutations, earth tides, and tidal variation in the lod was modeled. Results of this model can be used with Crustal Dynamics observations to look at the anelastic dissipation and dispersion at tidal periods. The effects of surface topography on various components of crustal deformation was also modeled, and numerical models were developed of post glacial rebound.