SEARCH · Search NASA

Results for “temperature dependence of ocean viscosity”

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.

SourceResearch connectionFrom yearThrough year

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

17 records

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Reduced viscosity steadily weakens oceanic currents

The viscosity of both air and water is temperature dependent. A rising temperature leads to an increased viscosity for air but a decreased viscosity for water. As climate becomes warmer, this increased air viscosity can partly inhibit the reduction of wind stress over the ocean, and the reduced water viscosity causes less downward momentum and heat transport. As these opposing effects of warming on air and water viscosity are not included in the state-of-the-art climate models, the understanding of their potential impacts on the response of the climate system to the anthropogenic warming is lacking. Here, via analyzing the Simple Ocean Data Assimilation oceanic reanalysis dataset, we show that the ocean heat content increases at a rate of ~1.3 × 10 22 J/yr over 35 years, which leads to a continuous reduction of oceanic viscosity. As a result, the ocean vertical shear enhances with a shoaling of the mixed layer depth and a reduced vertical linkage in the ocean. Our calculations show a reduction of the oceanic kinetic energy at a rate of ~2.4 × 10 16 J/yr. Potentially, this could generate far-reaching impacts on the energy storage of the climate system and, hence, could pace the global warming. Thus, it is important to include the temperature-dependent viscosity in our climate models. Freshwater discharged from polar ice sheets and mountain glaciers also contributes to the reduction in oceanic viscosity but, at present, to a lesser extent than that in oceanic warming. Reduced oceanic viscosity, therefore, is an important, but hitherto overlooked, response to a warming climate and contributes to many recent weather extremes including heavier rainfall rates in hurricanes, slackening of the polar vortex, and oceanic heat waves.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1981

Long wavelength gravity and topography anomalies

It is shown that gravity and topography anomalies on the earth's surface may provide new information about deep processes occurring in the earth, such as those associated with mantle convection. Two main reasons are cited for this. The first is the steady improvement that has occurred in the resolution of the long wavelength gravity field, particularly in the wavelength range of a few hundred to a few thousand km, mainly due to increased coverage of terrestrial gravity measurements and the development of radar altimeters in orbiting satellites. The second reason is the large number of numerical and laboratory experiments of convection in the earth, including some with deformable upper and lower boundaries and temperature-dependent viscosity. The oceans are thought to hold the most promise for determining long wavelength gravity and topography anomalies, since their evolution has been relatively simple in comparison with that of the continents. It is also shown that good correlation between long wavelength gravity and topography anomalies exists over some portions of the ocean floor

Watts, A. B.↗

NASA NTRS2014

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↗

NASA NTRS1976

Oceanic lithosphere and asthenosphere: The thermal and mechanical structure

A coupled thermal and mechanical solid state model of the oceanic lithosphere and asthenosphere is presented. The model includes vertical conduction of heat with a temperature dependent thermal conductivity, horizontal and vertical advection of heat, viscous dissipation or shear heating, and linear or nonlinear deformation mechanisms with temperature and pressure dependent constitutive relations between shear stress and strain rate. A constant horizontal velocity u sub 0 and temperature t sub 0 at the surface and zero horizontal velocity and constant temperature t sub infinity at great depth are required. In addition to numerical values of the thermal and mechanical properties of the medium, only the values of u sub 0, t sub 0 and t sub infinity are specified. The model determines the depth and age dependent temperature horizontal and vertical velocity, and viscosity structures of the lithosphere and asthenosphere. In particular, ocean floor topography, oceanic heat flow, and lithosphere thickness are deduced as functions of the age of the ocean floor.

Schubert, G.↗

NASA NTRS1976

Oceanic lithosphere and asthenosphere - Thermal and mechanical structure

A coupled thermomechanical subsolidus model of the oceanic lithosphere and asthenosphere is developed which includes vertical heat conduction, a temperature-dependent thermal conductivity, heat advection by a horizontal and vertical mass flow that depends on depth and age, contributions of viscous dissipation or shear heating, a linear or nonlinear deformation law relating shear stress and strain rate, as well as a temperature- and pressure-dependent viscosity. The model requires a constant horizontal velocity and temperature at the surface, but zero horizontal velocity and constant temperature at great depths. The depth- and age-dependent temperature, horizontal and vertical velocities, and viscosity structure of the lithosphere and asthenosphere are determined along with the age-dependent shear stress in those two zones. The ocean-floor topography, oceanic heat flow, and lithosphere thickness are deduced as functions of ocean-floor age; seismic velocity profiles which exhibit a marked low-velocity zone are constructed from the age-dependent geotherms and assumed values of the elastic parameters. It is found that simple boundary-layer cooling determines the thermal structure at young ages, while effects of viscous dissipation become more important at older ages.

Schubert, G.↗

NASA NTRS1982

The evolution of the moon - A finite element approach

The present lunar evolution model analyzes the thermal history of a self-gravitating spherical planetary body, including the effects of viscous dissipation, internal melting, adiabatic gradient, core formation, variable viscosity, radioactive nucleide decay, and a depth-dependent initial temperature profile, together with physical parameters corresponding to the moon. Although no initial basalt ocean is assumed, partial melting is observed early in the model moon's history. This is suggested to be related to the formation of the basalt maria. The model's present lithospheric thickness is 600 km, with core-mantle temperatures close to 1600 K and surface heat flux of 25.3 mW/sq m. The finite element method is judged to be applicable to the problem of planetary evolution, although faster solution algorithms will be required for the examination of a sufficient number of models.

Chacko, S.↗

NASA NTRS1978

Shear flow beneath oceanic plates - Local nonsimilarity boundary layers for olivine rheology

The principle of local similarity, which has been used to model the two-dimensional boundary layers in the oceanic upper mantle, permits calculation of the temperature, velocity, and stress fields with essentially analytic techniques. Finite difference numerical methods are hard pressed to resolve the detail required by the large variation of viscosity between the lithosphere and the asthenosphere. In this paper the local similarity approximation has been justified by quantitatively evaluating the effect of nonsimilarity due to viscous heating, nonlinear temperature- and pressure-dependent rheology, buoyancy, adiabatic cooling, etc. Nonsimilar effects produce only small modifications of the locally similar boundary layers; important geophysical observables such as surface heat flux and ocean floor topography are given to better than 10 percent by the locally similar solution. A posteriori evaluations of the terms neglected in the boundary layer simplification of the complete equations have been conducted on the locally similar temperature and velocity profiles close to the spreading ridge. The boundary layer models are valid to depths of 100 km at 3 m.y. and 10 km at 0.3 m.y.

Yuen, D. A.↗

NASA NTRS1990

Evolution of midplate hotspot swells: Numerical solutions

The evolution of midplate hotspot swells on an oceanic plate moving over a hot, upwelling mantle plume is numerically simulated. The plume supplies a Gaussian-shaped thermal perturbation and thermally-induced dynamic support. The lithosphere is treated as a thermal boundary layer with a strongly temperature-dependent viscosity. The two fundamental mechanisms of transferring heat, conduction and convection, during the interaction of the lithosphere with the mantle plume are considered. The transient heat transfer equations, with boundary conditions varying in both time and space, are solved in cylindrical coordinates using the finite difference ADI (alternating direction implicit) method on a 100 x 100 grid. The topography, geoid anomaly, and heat flow anomaly of the Hawaiian swell and the Bermuda rise are used to constrain the models. Results confirm the conclusion of previous works that the Hawaiian swell can not be explained by conductive heating alone, even if extremely high thermal perturbation is allowed. On the other hand, the model of convective thinning predicts successfully the topography, geoid anomaly, and the heat flow anomaly around the Hawaiian islands, as well as the changes in the topography and anomalous heat flow along the Hawaiian volcanic chain.

Liu, Mian↗

NASA NTRS1985

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

NASA NTRS1994

Mantle devolatilization and rheology in the framework of planetary evolution

We investigate the thermal history of an Earth-like planet with the help of a parameterized mantle convection model including the volatile exchange between mantle and surface reservoirs. The weakening of mantle silicates by dissolved volatiles is described by a functional relationship between creep rate and water fugacity. We use flow law parameters of diffusion creep in olivine under dry and wet conditions. The mantle degassing rate is considered as directly proportional to the seafloor spreading rate, which is also dependent on the mantle heat flow. To calculate the spreading rate, we assume that the heat flow under the mid-ocean ridges is double the average mantle heat flow. The rate of regassing also depends on the seafloor spreading rate as well as on other factors like the efficiency of volatile recycling through island arc volcanism. Both mechanisms (de- and regassing) are coupled self-consistently with the help of the parameterized convection model under implementation of a temperature and volatile-content-dependent mantle viscosity.

Franck, S.↗

NASA NTRS1984

The effects of the variations in sea surface temperature and atmospheric stability in the estimation of average wind speed by SEASAT-SASS

The average wind speeds from the scatterometer (SASS) on the ocean observing satellite SEASAT are found to be generally higher than the average wind speeds from ship reports. In this study, two factors, sea surface temperature and atmospheric stability, are identified which affect microwave scatter and, therefore, wave development. The problem of relating satellite observations to a fictitious quantity, such as the neutral wind, that has to be derived from in situ observations with models is examined. The study also demonstrates the dependence of SASS winds on sea surface temperature at low wind speeds, possibly due to temperature-dependent factors, such as water viscosity, which affect wave development.

Liu, W. T.↗

NASA NTRS1992

Estimation of the viscosities of silicate liquids at high pressure from measurements of oxygen diffusivity

The dynamics and evolution of a magma ocean depend significantly on the structure and physical properties of silicate liquids at high pressure. Viscosity is particularly important because of its effect on the dynamics of convection and crystal setting. Traditionally, the viscosities of silicate liquids have been measured at high pressure by falling sphere viscometry but, due to some limitations of this technique, few measurements have been made at pressures exceeding 2.5 GPa. An alternative approach, which has previously been proven up to 2 GPa, is to use high pressure measurements of oxygen self-diffusivity (D) to estimate viscosity (eta) from the Eyring relationship eta = kT/D lambda (where k is the Boltzmann constant, T is absolute temperature, and lambda is the diffusive jump distance). We have performed preliminary experiments on Na2Si4O9 liquid in order to explore the feasibility of this method up to pressures in excess of 10 GPa. Diffusion couples were prepared from glass starting materials with one half of the sample enriched in O-18 and the other half containing the natural abundance (0.2 wt percent). Diffusion experiments were performed at 1600 to 1800 C at pressures in the range of 2.5 to 10 GPa for times up to 6 minutes using a 1200 ton multianvil apparatus. Oxygen diffusivities were derived from the resulting O-18 concentration profiles, which were analyzed using an ion microprobe. At 1800 C, oxygen diffusivities increase continuously from 1 x 10(exp -10) m(exp -2) s(exp -1) at 2.5 Gpa to 5 x 10(exp -10) m(exp -2) s(exp -1) at 10 GPa. From the Eyring relationship, the viscosity of Na2SiO9 liquid is predicted to decrease by approx. 0.7 log units over this pressure range at 1800 C. These trends, which can be related to changes in the structure of the liquid at high pressure, are in agreement with results of molecular dynamics calculations. These results suggest that it should be possible to estimate the viscosities of silicate liquids up to at least 15 GPa and 2200 C by this method.

Rubie, D. C.↗

NASA NTRS1988

Dynamic compensation in the central Pacific Ocean

The intermediate-wavelength geoid (lambda similar to 2000 km) and sea-floor topography fields in the central Pacific Ocean were studied in terms of static and dynamic compensation models. Topographic features on the sea-floor with lambda less than 1000 km were found to be compensated both regionally, by the elastic strength of the lithosphere, and locally, by displacing mantle material to reach isostatic adjustment. The larger-scale sea-floor topography and the corresponding geoid anomalies with lambda similar to 2000 km cannot be explained by either local or regional compensation. The topography and the resulting geoid anomaly at this wavelength were modeled by considering the dynamic effects arising from viscous stresses in a layer of fluid with a highly temperature-dependent viscosity for the cases of: (1) surface cooling, and (2) basal heating. In this model, the mechanical properties of the elastic part of the lithosphere were taken into account by considering an activation energy of about 520 kJ/mol in the Arrhenius law for the viscosity. Numerical predictions of the topography, total geoid anomaly, and admittance were obtained, and the results show that the thermal perturbation in the layer, which accounts for the mass deficit, must be located close to the surface to compensate the gravitational effect of the surface deformation. For the case of basal heating, the temperature dependence of viscosity results in a separation of the upper, quasi-rigid lid from the lower mobile fluid, hence inhibiting the development of a compensating thermal perturbation at shallow depths. The results clearly rule out small-scale, upper-mantle convection as the source of these anomalies. Instead, the geophysical observables can be well explained by a shallow, transient thermal perturbation.

Hinojosa, Juan Homero↗

NASA NTRS2005

Lunar and Planetary Science XXXVI, Part 17

The following topics were discussed: A Model for the Formation of Paterae on Io; LIBS-based Detection of As, Br, C, Cl, P, and S in the VUV Spectral Region in a Mars Atmosphere; Mass Independent Sulfur in Achondrites: Possible Evidence of Photochemistry in the Solar Nebula; Grain Size-dependent Viscosity and Oceans in Icy Satellites; Claritas Paleolake Studied from the MEX HRSC Data; Mars Express HRSC Colors of White Rock, Arabia, Mars; Lava and Flows of the Arcadia Region of Mars; Isotopic Composition of Lunar Soils and the Early Differentiation of the Moon; Trace Element Analysis of Lunar Soils by ICP-MS; Highly Siderophile Elements and Osmium Isotope Systematics in Ureilites: Are the Carbonaceous Veins Primary Components?; Evaporative Evolution of Martian Brines Based on Halogens in Nakhlites and MER Samples; Io from High-Resolution Galileo PPR Data Taken Simultaneously with SSI or NIMS Observations; Loki, Io: Groundbased Observations and a Model for Periodic Overturn; Deconstructing a Few Myths in the Interpretation of Satellite-Altitude Crustal Magnetic Field: Examples from Mars Global Surveyor; Semi-Autonomous Rover Operations: A Mars Technology Program Demonstration; Rotational Studies of Asteroids with Small Telescopes; Mineralogy and Temperature-induced Spectral Investigations of A-type Asteroids 246 Asporina and 446 Aeternitas; and Thermal History Calculations Versus Full Convection Models: Application to the Thermal Evolution of Mercury. Recent Solar-Proton Fluxes

Source record↗

NASA NTRS1994

Subduction and volatile recycling in Earth's mantle

The subduction of water and other volatiles into the mantle from oceanic sediments and altered oceanic crust is the major source of volatile recycling in the mantle. Until now, the geotherms that have been used to estimate the amount of volatiles that are recycled at subduction zones have been produced using the hypothesis that the slab is rigid and undergoes no internal deformation. On the other hand, most fluid dynamical mantle flow calculations assume that the slab has no greater strength than the surrounding mantle. Both of these views are inconsistent with laboratory work on the deformation of mantle minerals at high pressures. We consider the effects of the strength of the slab using two-dimensional calculations of a slab-like thermal downwelling with an endothermic phase change. Because the rheology and composition of subducting slabs are uncertain, we consider a range of Clapeyron slopes which bound current laboratory estimates of the spinel to perovskite plus magnesiowustite phase transition and simple temperature-dependent rheologies based on an Arrhenius law diffusion mechanism. In uniform viscosity convection models, subducted material piles up above the phase change until the pile becomes gravitationally unstable and sinks into the lower mantle (the avalanche). Strong slabs moderate the 'catastrophic' effects of the instabilities seen in many constant-viscosity convection calculations; however, even in the strongest slabs we consider, there is some retardation of the slab descent due to the presence of the phase change.

King, S. D.↗

NASA NTRS2018

High-resolution Wave Propagation Method for Stratified Flows

The implementation of the multidimensional f-waves Riemann solver for the time-dependent, three-dimensional, nonhydrostatic, meso- and microscale atmospheric flows is described in detail. The Riemann solver employs flux-based wave decomposition (f-waves) for the calculation of Godunov fluxes in which the flux differences are written directly as the linear combination of the right eigenvectors of the hyperbolic system. The scheme incorporates the source term due to gravity without introducing discretization errors which is an important property in the context of atmospheric flows. The resulting flow solver is conservative, accurate, stable, and well-balanced. The implementation of the solver is evaluated using benchmark test cases for atmospheric dynamics.

Riemann problem↗

NASA NTRS2010

Stability of Basalt plus Anhydrite plus Calcite at HP-HT: Implications for Venus, the Earth and Mars

"Canali" observed at Venus surface by Magellan are evidence for very long melt flows, but their composition and origin remain uncertain. The hypothesis of water-rich flow is not reasonable regarding the temperature at Venus surface. The length of these channels could not be explained by a silicate melt composition but more likely, by a carbonate-sulfate melt which has a much lower viscosity (Kargel et al 1994). One hypothesis is that calcite CaCO3 and anhydrite CaSO4 which are alteration products of basalts melted during meteorite impacts. A famous example recorded on the Earth (Chicxulub) produced melt and gas rich in carbon and sulfur. Calcite and sulfate evaporites are also present on Mars surface, associated with basalts. An impact on these materials might release C- and S-rich melt or fluid. Another type of planetary phenomenon (affecting only the Earth) might provoke a high pressure destabilization of basalt+anhydrite+calcite. Very high contents of C and S are measured in some Earth s magmas, either dissolved or in the form of crystals (Luhr 2008). As shown by the high H content and high fO2 of primary igneous anhydrite-bearing lavas, the high S content in their source may be explained by subduction of an anhydrite-bearing oceanic crust, either directly (by melting followed by eruption) or indirectly (by release of S-rich melt or fluid that metasomatize the mantle) . Calcite is a major product of oceanic sedimentation and alteration of the crust. Therefore, sulfate- and calcite-rich material may be subducted to high pressures and high temperatures (HP-HT) and release S- and C-rich melts or fluids which could influence the composition of subduction zone lavas or gases. Both phenomena - meteorite impact and subduction - imply HP-HT conditions - although the P-T-time paths are different. Some HP experimental/theoretical studies have been performed on basalt/eclogite, calcite and anhydrite separately or on a combination of two. In this study we performed piston-cylinder experiments at 1 GPa between 950 and 1700 C using a mixture of 70wt% tholeiitic basalt + 15wt% anhydrite + 15wt% calcite. Up to ~1440 C, an ultracalcic (CaO>13.5 wt%; CaO/Al2O3>1 wt%) picrobasaltic (SiO2~40-45 wt%; Na2O+K2O<2 wt%) melt containing up to 5 wt% SO3 and up to 5.3 wt% CO2+H2O (calculated by difference) is present in equilibrium with clinopyroxene, anhydrite, spinel-chromite, a CAS-phase and a gas composed mainly of CO and an aliphatic thiol (CH2)4SH. Hydrogen was incorporated either by contact between the starting material and air or by diffusion through the capsule during the experiments. The S content in the gas increases with temperature and run duration, implying that gases with various C/S ratios might be released during an impact or at subduction zones, depending on the P-T-t path and on the H content. Above approx.1440 C, a Ca-rich carbonate-sulfate melt forms (in equilibrium with the picrobasaltic melt) which contains a few percents of Na and K. Such melt is not expected to form at Earth s subduction temperatures. If it forms by meteorite impact, it might crystallize too fast to explain long flows like Venus canali. A different basalt/anhydrite/calcite ratio might, however, decrease its formation temperature.

Martin, A. M.↗