Hager, B. H. | Search NASA

NASA NTRS1989

Long-wavelength variations in earth's geoid - Physical models and dynamical implications

Seismic tomography and fluid dynamical models are extended to include a detailed and realistic parameterization of the viscosity structure of the mantle. The model also includes the effects of heterogeneity inferred from surface-wave studies of upper-mantle structure. It is found that about 90 percent of the observed geoid may be explained by density anomalies inferred from tomography, a model of subducted slabs, and an earth model with a low-viscosity asthenosphere over a moderate-viscosity lower mantle.

Hager, B. H.↗

NASA NTRS1988

Subduction zone earthquakes and stress in slabs

Simple viscous fluid models of subducting slabs are used to explain observations of the distribution of earthquakes as a function of depth and the orientation of stress axes of deep (greater than 300 km) and intermediate (70-300 km) earthquakes. Results suggest the following features in the distribution of earthquakes with depth: (1) an exponential decrease from shallow depths down to 250 to 300 km, (2) a minimum near 250 to 300 km, and (3) a deep peak below 300 km. Many shallow subducting slabs show only the first characteristic, while deeper extending regions tend to show all three features, with the deep peak varying in position and intensity. These data, combined with the results on the stress orientations of various-depth earthquakes, are consistent with the existence of a barrier of some sort at 670-km depth and a uniform viscosity mantle above this barrier.

Vassiliou, M. S.↗

NASA NTRS1987

Solving Finite-Element Problems on a Concurrent Processor

By use of "conjugate gradients" technique, concurrent efficiency greater than 90 percent. Algorithm applies method of conjugate gradients to iterative solution of finite-element problems on concurrent processor. With algorithm, iteration rates nearly proportional to number of processors. For sufficiently large problems, fraction of proportional speedup achieved, called concurrent efficiency, exceeds 90 percent. Results indicate future application of this and related algorithms to large finite-element problems depend primarily upon applicability of iteratiave techniques, not upon issues of concurrency or efficiency.

Lyzenga, G. A.↗

NASA NTRS1986

Dynamically supported geoid highs over hotspots: Observation and theory

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.

Richards, M. A.↗

NASA NTRS1986

A dynamic model of Venus's gravity field

Unlike Earth, long wavelength gravity anomalies and topography correlate well on Venus. Venus's admittance curve from spherical harmonic degree 2 to 18 is inconsistent with either Airy or Pratt isostasy, but is consistent with dynamic support from mantle convection. A model using whole mantle flow and a high viscosity near surface layer overlying a constant viscosity mantle reproduces this admittance curve. On Earth, the effective viscosity deduced from geoid modeling increases by a factor of 300 from the asthenosphere to the lower mantle. These viscosity estimates may be biased by the neglect of lateral variations in mantle viscosity associated with hot plumes and cold subducted slabs. The different effective viscosity profiles for Earth and Venus may reflect their convective styles, with tectonism and mantle heat transport dominated by hot plumes on Venus and by subducted slabs on Earth. Convection at degree 2 appears much stronger on Earth than on Venus. A degree 2 convective structure may be unstable on Venus, but may have been stabilized on Earth by the insulating effects of the Pangean supercontinental assemblage.

Kiefer, W. S.↗

NASA NTRS1985

A critical assessment of viscous models of trench topography and corner flow

Stresses for Newtonian viscous flow in a simple geometry (e.g., corner flow, bending flow) are obtained in order to study the effect of imposed velocity boundary conditions. Stress for a delta function velocity boundary condition decays as 1/R(2); for a step function velocity, stress goes as 1/R; for a discontinuity in curvature, the stress singularity is logarithmic. For corner flow, which has a discontinuity of velocity at a certain point, the corresponding stress has a 1/R singularity. However, for a more realistic circular-slab model, the stress singularity becomes logarithmic. Thus the stress distribution is very sensitive to the boundary conditions, and in evaluating the applicability of viscous models of trench topography it is essential to use realistic geometries. Topography and seismicity data from northern Hoshu, Japan, were used to construct a finite element model, with flow assumed tangent to the top of the grid, for both Newtonian and non-Newtonian flow (power law 3 rheology). Normal stresses at the top of the grid are compared to the observed trench topography and gravity anomalies. There is poor agreement. Purely viscous models of subducting slables with specified velocity boundary conditions do not predict normal stress patterns compatible with observed topography and gravity. Elasticity and plasticity appear to be important for the subduction process.

Zhang, J.↗

NASA NTRS1985

The source of the Earth's long wavelength geoid anomalies: Implications for mantle and core dynamics

The long wavelength components of the Earth's gravity field result mainly from density contrasts associated with convection in the mantle. Direct interpretation of the geoid for mantle convection is complicated by the fact that convective flow results in dynamically maintained deformation of the surface of the Earth, the core mantle boundary (CMB), and any interior chemical boundaries which might exist. These boundary deformations effect the geoid opposite in sign and are comparable in magnitude to those of the interior density contrasts driving the flow. The total difference of two relatively large quantities.

Hager, B. H.↗

NASA NTRS1985

Constraints on the rheological structure of the mantle

Rheological models of the mantle are at present limited to radial symmetry, usually with homogeneous linearly viscous or viscoelastic incompressible layers. While such models are probably overly simple, they readily allow calculation of geophysical effects, such as post glacial rebound and related changes in the Earth's shape. They are also directly applicable to problems of global mantle flow and plate motions.

Oconnell, R. J.↗

NASA NTRS1985

Small-scale convection beneath the transverse ranges, California: Implications for interpretation of gravity anomalies

Tomographic inversion of upper mantle P wave velocity heterogeneities beneath southern California shows two prominent features: an east-west trending curtain of high velocity material (up to 3% fast) in the upper 250 km beneath the Transverse Ranges and a region of low velocity material (up to 4% slow) in the 100 km beneath the Salton Trough. These seismic velocity anomalies were interpreted as due to small scale convection in the mantle. Using this hypothesis and assuming that temperature and density anomalies are linearly related to seismic velocity anomalies through standard coefficients of proportionality, leads to inferred variations of approx. + or - 300 C and approx. + or - 0.03 g/cc.

Humphreys, E. D.↗

NASA NTRS1985

Lower mantle heterogeneity, dynamic topography and the geoid

Density contrasts in the lower mantle, recently imaged using seismic tomography, drive convective flow which results in kilometers of dynamically maintained topography at the core-mantle boundary and at the earth's surface. The total gravity field due to interior density constrasts and boundary topography predicts the largest wavelength components of the geoid remarkably well. Neglecting dynamic surface deformation leads to geoid anomalies of opposite sign than are observed.

Hager, B. H.↗

NASA NTRS1985

Dynamic compensation of Venus's geoid: A comparison with Earth

Unlike Earth, on Venus long wavelength geoid anomalies correlate well with topography. Venus's admittance curve between harmonic degrees 3 and 18 is inconsistent with Airy isostasy but is consistent with dynamic support from convection being the dominant mechanism of compensation on Venus. We model dynamic compensation on Venus using simple flow models which assume a spherically symmetric Newtonian mantle viscosity profile. Preliminary models parameterize the viscosity variation with depth as a 2 layer model with a boundary at 720 km depth. A model in which viscosity in the lower mantle is a factor of 10 lower than in the upper mantle can explain Venus's observed admittance curve for degrees 3 through 18. Dynamic models which include a chemical boundary between the upper and lower mantle do not successfully explain the observed admittance curve, indicating that Venus does not have a chemically layered mantle.

Kiefer, W. S.↗

NASA NTRS1984

A dynamic model of Venus's gravity field

Unlike Earth, long wavelength gravity anomalies and topography correlate well on Venus. Venus's admittance curve from spherical harmonic degree 2 to 18 is inconsistent with either Airy or Pratt isostasy, but is consistent with dynamic support from mantle convection. A model using whole mantle flow and a high viscosity near surface layer overlying a constant viscosity mantle reproduces this admittance curve. On Earth, the effective viscosity deduced from geoid modeling increases by a factor of 300 from the asthenosphere to the lower mantle. These viscosity estimates may be biased by the neglect of lateral variations in mantle viscosity associated with hot plumes and cold subducted slabs. The different effective viscosity profiles for Earth and Venus may reflect their convective styles, with tectonism and mantle heat transport dominated by hot plumes on Venus and by subducted slabs on Earth. Convection at degree 2 appears much stronger on Earth than on Venus. A degree 2 convective structure may be unstable on Venus, but may have been stabilized on Earth by the insulating effects of the Pangean supercontinental assemblage.

Kiefer, W. S.↗

NASA NTRS1984

Geoid anomalies in a dynamic earth

Potential and surface deformation Love numbers for internal loading have been calculated in order to obtain a dynamically consistent relationship between the geoid and the earth's response to internal buoyancy forces. These quantities depend on the depth and harmonic degree of loading, and can be integrated as Green functions to obtain the dynamic response due to an arbitrary distribution of internal density contrasts. Constructing a series of spherically symmetric, self-gravitating flow models for a variety of radial Newtonian viscosity variations and flow configurations, and calculating relaxation times for spherically symmetric viscous earth models, it is demonstrated that boundary deformation due to internal loading reaches its equilibrium value on the same time scale as postglacial rebound; this is much less time than the time scale for significant change in the convective flow pattern.

Richards, M. A.↗

NASA NTRS1984

Subducted slabs and the geoid - Constraints on mantle rheology and flow

The total geoid anomaly which is the result of a given density contrast in a convecting viscous earth is affected by the mass anomalies associated with the flow induced deformation of the upper surface and internal compositional boundaries, as well as by the density contrast itself is discussed. If the internal density contrasts can be estimated, the depth and variation of viscosity with depth of the convecting system can be constrained. The observed long wavelength geoid is highly correlated with that predicted by a density model for seismically active subducted slabs. The amplitude of the correlation is explained if the density contrasts associated with subduction extend into the lower mantle or if subducted slabs exceeding 350 km in thickness are piled up over horizontal distances of thousands of km at the base of the upper mantle. Mantle wide convection in a mantle that has a viscosity increasing with depth provides the explanation of the long-wavelength geoid anomalies over subduction zones. Previously announced in STAR as N83-22874

Hager, B. H.↗

NASA NTRS1984

A tomographic image of mantle structure beneath southern California

The variations in seismic structure beneath southern California were determined by using a tomographic method of inversion on teleseismic P delays recorded with the Southern California Array. The inversion reveals two prominent features beneath the region. The first is a thin, vertical wedge directly beneath the Transverse Ranges that is 2-3 percent faster than the surrounding region. This feature deepens to the east, attaining a maximum depth of about 250 km beneath the San Bernardino Mountains. The second feature is a major zone of low velocity material that is 2-4 percent slow under the Salton Trough rift valley, extending to a depth of about 125 km. Two possible explanations for the spatial association of the Transverse Ranges with the velocity anomaly below are lithospheric subduction or small-scale sublithospheric convection in the region of the Big Bend of the San Andreas Fault. The low velocity anomaly beneath the Salton Trough is consistent with convective upwelling there.

Humphreys, E.↗

NASA NTRS1984

Lower mantle heterogeneity, dynamic topography and the geoid

Density contrasts in the lower mantle, recently imaged using seismic tomography, drive convective flow which results in kilometers of dynamically maintained topography at the core-mantle boundary and at the Earth's surface. The total gravity field due to interior density contrasts and boundary topography predicts the largest wavelength components of the geoid remarkably well. Neglecting dynamic surface deformation leads to geoid anomalies of opposite sign than are observed.

Hager, B. H.↗

NASA NTRS1984

A critical assessment of viscous models of trench topography and corner flow

Stresses for Newtonian viscous flow in a simple geometry (e.g., corner flow, bending flow) are obtained in order to study the effect of imposed velocity boundary conditions. Stress for a delta function velocity boundary condition decays as 1/R(2); for a step function velocity, stress goes as 1/R; for a discontinuity in curvature, the stress singularity is logarithmic. For corner flow, which has a discontinuity of velocity at a certain point, the corresponding stress has a 1/R singularity. However, for a more realistic circular-slab model, the stress singularity becomes logarithmic. Thus the stress distribution is very sensitive to the boundary conditions, and in evaluating the applicability of viscous models of trench topography it is essential to use realistic geometries. Topography and seismicity data from northern Hoshu, Japan, were used to construct a finite element model, with flow assumed tangent to the top of the grid, for both Newtonian and non-Newtonian flow (power law 3 rheology). Normal stresses at the top of the grid are compared to the observed trench topography and gravity anomalies. There is poor agreement. Purely viscous models of subducting slables with specified velocity boundary conditions do not predict normal stress patterns compatible with observed topography and gravity. Elasticity and plasticity appear to be important for the subduction process.

Zhang, J.↗

NASA NTRS1984

Convection experiments in a centrifuge and the generation of plumes in a very viscous fluid

In this paper, experiments are described for which inertial effects are negligible. A small aspect-ratio tank filled with a very viscous fluid (Pr = 10 to the 6th) is used to observe the behavior of convection for Rayleigh numbers up to 6.3 x 10 to the 5th. These high values are reached by conducting the experiment in a centrifuge which provides a 130-fold increase in apparent gravity. Rotational effects are small, but cannot be totally dismissed. In this geometry, thermal boundary layer instabilities are indeed observed, and are found to be very similar to their lower Prandtl number counterparts. It is tentatively concluded that once given a certain degree of 'vulnerability' convection can develop 'plume' like instabilities, even when the Prandtl number is infinite. The concept is applied to the earth's mantle and it is speculated that 'plumes' could well be the dominant mode of small-scale convection under the lithospheric plates.

Nataf, H.-C.↗