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.
During 1983 to 1986 emphasis in the study of mantle convection shifted away from fluid mechanical analysis of simple systems with uniform material properties and simple geometries, toward analysis of the effects of more complicated, presumably more realistic models. The important processes related to mantle convection are considered. The developments in seismology are discussed.
Model calculations are described for a compressible fluid in a three-dimensional spherical shell with 80 percent of the surface heat flow generated within the model mantle. The numerical solutions are strongly chaotic, with surface planforms dominated by long curvilinear downflows reminiscent of the descending slabs in the earth's mantle. The results suggest that descending slabs play an important part in driving mantle convection, and that their chaotic evolution may influence the spatial and temporal behavior of plates and thus the dispersal and aggregation of continents.
Dehydration reactions within subducted oceanic crust are important for fluid-mediated element transfer within the subducting plate and potentially to the mantle wedge. The effects of metamorphic reactions and fluid flow on element recycling that occur during retrogression and exhumation of subducted continental crust from mantle depths are poorly understood. We study two metabasite pods with fresh eclogite cores and retrogressed amphibolite-facies rims and surrounding host gneiss within the Western Gneiss Region (WGR), Norway, to better understand element mobility and mass transfer during exhumation of subducted continental crust. Bulk-rock data were collected from samples taken across the pod and into the host gneiss. Phengite breakdown in eclogite and epidote recrystallization in veins and/or gneiss within pod cores contributed large ion lithophile elements and REE to retrogressed eclogite closest to the pod cores. In gneiss hosting the pods, phengite and epidote breakdown provided fluid that mediated elemental transfer and redistribution to the pod rim or tail. Compared to the studied pod in the southern WGR, the pod in the northern WGR underwent higher P-T conditions, partial melting and higher strain rates. This resulted in the infiltration of external fluid farther into the pod interior from the rim and facilitated larger mass gain in trace elements in the amphibolite tail of the pod relative to fresh eclogite in the core. The results show clear evidence for retrogression dehydration reactions driving significant fluid-mediated element redistribution as observed on the outcrop scale during exhumation following ultrahigh-pressure metamorphism, which directly impacts element signatures within the exhuming crust.
Oxygen fugacity and water content are crucial parameters for many chemical and physical properties of the Earth's mantle, for example bearing on fluid type, melting initiation, and deformation. However, the exact behaviour of Fe3+ and H during melting and metasomatism is still under debate. Here, the Fe3+/ΣFe ratio (Mössbauer and EMP) and water content (FTIR) of peridotite minerals are examined in mantle xenoliths from Kilbourne Hole (KH), NM, and Dish Hill (DH), CA (USA). These spinel peridotites have compositions consistent with partial melting with variable degrees of metasomatism (undetectable to cryptic to modal). Pyroxenites also allow to examine melt-rock reactions. Bulk-rock Fe2O3 content of the KH peridotites correlates with indices of melting (positive with bulk-rock Al2O3 and Cpx Yb contents, and negative with spinel Cr#) confirming that Fe3+ behaves as an incompatible element during melting. Correlations of the Fe3+/ΣFe ratio of minerals with these indices, however, indicates that Fe3+ is incompatible in Cpx but compatible in Opx and spinel during melting. Water contents in olivine, Cpx and Opx from most KH peridotites can be explained by partial melting and correlate negatively with the Fe3+/ΣFe ratio of spinel and Opx but positively with that of Cpx. This indicates partial control of Fe3+ on the incorporation of H in pyroxene, but not related to a redox equilibrium in Cpx. The higher Fe3+/ΣFe ratio of spinel in the metasomatized KH and DH peridotites, and in the pyroxenites confirms that oxidation characterizes modal metasomatism. Metasomatism, however, is not necessarily accompanied by water addition.
Large-scale mass transports in the geophysical fluids of the Earth system excite Earth's rotational variations in both length-of-day and polar motion. The excitation process is via the conservation of angular momentum. Therefore Earth rotation observations contain information about the integrated angular momentum (consisting of both the mass term and the motion term) of the geophysical fluids, which include atmosphere, hydrosphere, mantle, and the outer and inner cores. Such global information is often important and otherwise unattainable depending on the nature of the mass transport, its magnitude and time scale. The last few years have seen great advances in VLBI measurement of Earth rotation in precision and temporal resolution. These advances have opened new. areas in geophysical fluid studies, such as oceanic tidal angular momentum, atmospheric tides, Earth librations, and rapid atmospheric angular momentum fluctuations. Precision of 10 microseconds in UTI and 200 microarcseconds in polar motion can now be achieved on hourly basis. Building upon this heritage, the multi-network geodetic VLBI project, Continuous Observation of the Rotation of the Earth (CORE), promises to further these studies and to make possible studies on elusive but tell-tale geophysical processes such as oscillatory modes in the core and in the atmosphere. Currently the early phase of CORE is underway. Within a few years into the new mellinnium, the upcoming space gravity missions (such as GRACE) will measure the temporal variations in Earth's gravitational field, thus providing complementary information to that from Earth rotation study for a better understanding of global geophysical fluid processes.
Numerical calculations in two dimensions show that solitary wave disturbances can propagate along thermal plumes in a homogeneous, viscous fluid with a thermally activated rheology similar to the mantle. Comparison of the numerical results with analytic results from an idealized two-fluid model indicates that the two-fluid model can be applied to solitary waves in the thermally activated plume. A two-fluid model of solitary waves on cylindrical conduits is accordingly applied to mantle plumes to estimate the propagation speeds, time durations, and pulse lengths of solitary waves in mantle plumes as a function of background mantle viscosity, plume flux, plume density deficit, plume viscosity, and the volume of material transported by the solitary wave. It is inferred that mantle plumes may be eposodic via solitary waves and that these disturbances might account for observed 0(10 m.y.) variations in the durations and spacings of episodes of enhanced hotspot volcanic activity. Solitary waves in mantle plumes could be generated by interactions among plumes and between plumes and the large-scale, time-dependent mantle circulation.
Abstract Subduction‐related lavas have higher Fe 3+ /∑Fe than midocean ridge basalts (MORB). Hypotheses for this offset include imprint from subducting slabs and differentiation in thickened crust. These ideas are readily tested through examination of the time‐dependent evolution of slab‐derived signatures, thickening crust of the overriding plate, and evolving redox during subduction initiation. Here, we present Fe 3+ /ΣFe and volatile element abundances of volcanic glasses recovered from International Ocean Discovery Program (IODP) Expedition 352 to the Izu‐Bonin‐Mariana (IBM) forearc. The samples include forearc basalts (FAB) that are stratigraphically overlain by low‐ and high‐silica boninite lavas. The FAB glasses have 0.18–0.85 wt% H 2 O, 75–233 ppm CO 2 , S contents controlled by saturation with a sulfide phase (602–1,386 ppm), Ba/La from 3.9‐10, and Fe 3+ /ΣFe ratios from 0.136 to 0.177. These compositions are similar to MORB and suggest that decompression melting of dry and reduced mantle dominates the earliest stages of subduction initiation. Low‐ and high‐silica boninite glasses have 1.51–3.19 wt% H 2 O, CO 2 below detection, S contents below those required for sulfide saturation (5–235 ppm), Ba/La from 11 to 29, and Fe 3+ /∑Fe from 0.181 to 0.225. The compositions are broadly similar to modern arc lavas in the IBM arc. These data demonstrate that the establishment of fluid‐fluxed melting of the mantle, which occurs in just 0.6–1.2 my after subduction initiation, is synchronous with the production of oxidized, mantle‐derived magmas. The coherence of high Fe 3+ /∑Fe and Ba/La ratios with high H 2 O contents in Expedition 352 glasses and the modern IBM arc rocks strongly links the production of oxidized arc magmas to signatures of slab dehydration.
Variations in rotation and orientation of the Moon are sensitive to solid-body tidal dissipation, dissipation due to relative motion at the fluid-core/solid-mantle boundary, and tidal Love number k2 [1,2]. There is weaker sensitivity to flattening of the core-mantle boundary (CMB) [2-5] and fluid core moment of inertia [1]. Accurate Lunar Laser Ranging (LLR) measurements of the distance from observatories on the Earth to four retroreflector arrays on the Moon are sensitive to lunar rotation and orientation variations and tidal displacements. Past solutions using the LLR data have given results for dissipation due to solid-body tides and fluid core [1] plus Love number [1-5]. Detection of CMB flattening has been improving [3,5] and now seems significant. This strengthens the case for a fluid lunar core.
Satellite systems do the best job of resolving the long wavelength components of the Earth's gravity field. Over the oceans, satellite-borne radar altimeters such as SEASAT provide the best resolution observations of the intermediate wavelength components. Satellite observations of gravity contributed to the understanding of the dynamics of subduction. Large, long wavelength geoidal highs generally occur over subduction zones. These highs are attributed to the superposition of two effects of subduction: (1) the positive mass anomalies of subducting slabs themselves; and (2) the surface deformations such as the trenches convectively inducted by these slabs as they sink into the mantle. Models of this subduction process suggest that the mantle behaves as a nonNewtonian fluid, its effective viscosity increases significantly with depth, and that large positive mass anomalies may occur beneath the seismically defined Benioff zones.
The predicted orientation of Venus' rotation axis relative to its orbit ran be uniquely determined given knowledge of its J2 gravity coefficient and polar moment of inertia C if its free obliquity is fully damped. This assumption seems warranted given the dominant damping mechanism: turbulent fluid friction at a core mantle boundary (CMB). This skin friction results from differential obliquity of mantle and core spin axes, and the associated damping rate could be as short as 1/10(exp 6) year. However, the observed pole orientation indicates a free obliquity amplitude epsilon approx. = 2.1 deg compared with a nominal forced amplitude of 0.5 deg. There are two plausible explanations. The most likely is that the observed obliquity is a tidally evolved end state in which core friction, modulated by CMB ellipticity and core obliquity amplitude, counterbalances solid and atmospheric tidal torques. This concept is similar to the explanation for the retrograde spin omega as an end state in which solid and atmospheric thermal tidal torques balance at the present spin rate because of the omega(sup -1) dependence of the axial thermal torque. Weaker effects such as solid and thermal tides can then compete with core friction and for plausible models, their sum tends to increase free obliquity. The obliquity balance is controlled by the nonlinear (and nearly quadratic) dependence of the CMB turbulent "skin friction" torque on obliquity.Core ellipticity also has a profound effect on tidal evolution of Venus obliquity.
The predicted orientation of Venus' rotation axis relative to its orbit can be uniquely determined given knowledge of its J2 gravity coefficient and polar moment of inertia C if its free obliquity is fully damped. This assumption seems warranted given the dominant damping mechanism: turbulent fluid friction at a core mantle boundary. The determination of the free obliquity is examined.
It is noted that dissipation in the moon produces a small offset (approximately 0.23 arcsec) of the moon's rotation axis from the plane defined by the ecliptic and lunar orbit normals. Both solid body tidal friction and viscous fluid friction at a core-mantle interface are thought to be plausible mechanisms. The merits of both are discussed, and it is found that solid friction requires a low lunar tidal Q (approximately 28), whereas turbulent fluid friction requires a core with a radius of approximately 330 km to cause the signature observed by lunar laser ranging. Large (approximately 0.4-8.0 arcsec) free librations of the lunar figure have also been detected through laser ranging. Both a very recent impact on the moon and fluid turbulence in the lunar core are considered plausible mechanisms for generating these librations.
Combined inferences from seismology, high-pressure experiment and theory, geomagnetism, fluid dynamics, and current views of terrestrial planetary evolution lead to models of the earth's core with five basic properties. These are that core formation was contemporaneous with earth accretion; the core is not in chemical equilibrium with the mantle; the outer core is a fluid iron alloy containing significant quantities of lighter elements and is probably almost adiabatic and compositionally uniform; the more iron-rich inner solid core is a consequence of partial freezing of the outer core, and the energy release from this process sustains the earth's magnetic field; and the thermodynamic properties of the core are well constrained by the application of liquid-state theory to seismic and labroatory data.
Variations in rotation and orientation of the Moon are sensitive to solid-body tidal dissipation, dissipation due to relative motion at the fluid-core/solid-mantle boundary, and tidal Love number k2 [1,2]. There is weaker sensitivity to flattening of the core-mantle boundary (CMB) [2,3,4] and fluid core moment of inertia [1]. Accurate Lunar Laser Ranging (LLR) measurements of the distance from observatories on the Earth to four retroreflector arrays on the Moon are sensitive to lunar rotation and orientation variations and tidal displacements. Past solutions using the LLR data have given results for dissipation due to solid-body tides and fluid core [1] plus Love number [1-5]. Detection of CMB flattening, which in the past has been marginal but improving [3,4,5], now seems significant. Direct detection of the core moment has not yet been achieved.
Seismic studies suggest that the bulk of the mantle is heterogeneous, with density variations in depth as well as in horizontal directions (latitude and longitude). This density variation produces a three- dimensional gravity field throughout the Earth. On the other hand, the core density also varies in both time and space, due to convective core flow. Consequently, the fluid outer core and the solid mantle interact gravitationally due to the mass anomalies in both regions. This gravitational core-mantle interaction could play a significant role in exchange of angular momentum between the core and the mantle, and thus the change in Earth's rotation on time scales of decades and longer. Aiming at estimating the significance of the gravitational core-mantle interaction on Earth's rotation variation, we introduce in our MoSST core dynamics model a heterogeneous mantle, with a density distribution derived from seismic results. In this model, the core convection is driven by the buoyancy forces. And the density variation is determined dynamically with the convection. Numerical simulation is carried out with different parameter values, intending to extrapolate numerical results for geophysical implications.
Most known porphyry Cu deposits formed in the Phanerozoic and are exclusively associated with moderately oxidized, sulfur-rich, hydrous arc-related magmas derived from partial melting of the asthenospheric mantle metasomatized by slab-derived fluids. Yet, whether similar metallogenic processes also operated in the Precambrian remains obscure. Here we address the issue by investigating the origin, fO 2 , and S contents of calc-alkaline plutonic rocks associated with the Haib porphyry Cu deposit in the Paleoproterozoic Richtersveld Magmatic Arc (southern Namibia), an interpreted mature island-arc setting. We show that the ca. 1886–1881 Ma ore-forming magmas, originated from a mantle-dominated source with minor crustal contributions, were relatively oxidized (1–2 log units above the fayalite-magnetite-quartz redox buffer) and sulfur-rich. These results indicate that moderately oxidized, sulfur-rich arc magma associated with porphyry Cu mineralization already existed in the late Paleoproterozoic, probably as a result of recycling of sulfate-rich seawater or sediments from the subducted oceanic lithosphere at that time.
Observational results for the earth's C21 and S21 gravity coefficients can be used to constrain the mean equatorial rotation of the core with respect to the mantle. Current satellite gravity solutions suggest the equatorial rotation rate is no larger than 1 x 10 to the -7th times the earth's diurnal spin rate, a limit more than one order of magnitude smaller than the polar rotation rate inferred from the westward drift of the earth's magnetic field. The next generation gravity solutions should improve this constraint by more than one order of magnitude. Implications for the fluid pressure at the core-mantle boundary and for the shape of that boundary are discussed.
It is noted that the trace element concentration patterns of continental and ocean island basalts and of mid-ocean ridge basalts are complementary. Estimates of the relative sizes of the source regions for these fundamentally different basalt types can be arrived at from the trace element enrichment-depletion patterns. Their combined volume occupies the greater part of the mantle above the 670 km discontinuity. It is pointed out that the source regions separated as a result of early mantle differentiation and crystal fractionation from the resulting melt. The mid-ocean ridge basalts source evolved from an eclogite cumulate that gave up its late-stage enriched fluids at various times to the shallower mantle and continental crust. The mid-ocean ridge basalts source is rich in garnet and clinopyroxene, while the continental and ocean island basalt source is a garnet peridotite that has experienced secondary enrichment. These relationships are found to be consistent with the evolution of a terrestrial magma ocean.