Search NASA⌕ Search

Engineering topics

Anderson, Charles

Publications and source records attributed to Anderson, Charles.

Mantle plume interaction with an endothermic phase change

High spatial resolution numerical simulations of mantle plumes impinging from below on the endothermic phase change at 660-km depth are used to investigate the effects of latent heat release on the plume-phase change interaction. Both axisymmetric and planar upflows are considered, and the strong temperature dependence of mantle viscosity is taken into account. For plume strengths considered, a Clapeyron slope of -4 MPa/K prevents plume penetration of the phase change. Plumes readily penetrate the phase change for a Clapeyron slope of -2 MPa/K and arrive in the upper mantle considerably hotter than if they had not traversed the phase change. For the same amount of thermal drive, i.e., the same excess basal temperature, axisymmetric plumes are hotter upon reaching the upper mantle than are planar upwellings. Heating of plumes by their passage through the spinel-perovskite endothermic phase change can have important consequences for the ability of the plume to thermally thin the lithosphere and cause melting and volcanism.

Schubert, Gerald↗

Modeling of periodic great earthquakes on the San Andreas fault: Effects of nonlinear crustal rheology

We analyze the cycle of great earthquakes along the San Andreas fault with a finite element numerical model of deformation in a crust with a nonlinear viscoelastic rheology. The viscous component of deformation has an effective viscosity that depends exponentially on the inverse absolute temperature and nonlinearity on the shear stress; the elastic deformation is linear. Crustal thickness and temperature are constrained by seismic and heat flow data for California. The models are for anti plane strain in a 25-km-thick crustal layer having a very long, vertical strike-slip fault; the crustal block extends 250 km to either side of the fault. During the earthquake cycle that lasts 160 years, a constant plate velocity v(sub p)/2 = 17.5 mm yr is applied to the base of the crust and to the vertical end of the crustal block 250 km away from the fault. The upper half of the fault is locked during the interseismic period, while its lower half slips at the constant plate velocity. The locked part of the fault is moved abruptly 2.8 m every 160 years to simulate great earthquakes. The results are sensitive to crustal rheology. Models with quartzite-like rheology display profound transient stages in the velocity, displacement, and stress fields. The predicted transient zone extends about 3-4 times the crustal thickness on each side of the fault, significantly wider than the zone of deformation in elastic models. Models with diabase-like rheology behave similarly to elastic models and exhibit no transient stages. The model predictions are compared with geodetic observations of fault-parallel velocities in northern and central California and local rates of shear strain along the San Andreas fault. The observations are best fit by models which are 10-100 times less viscous than a quartzite-like rheology. Since the lower crust in California is composed of intermediate to mafic rocks, the present result suggests that the in situ viscosity of the crustal rock is orders of magnitude less the rock viscosity determined in the laboratory.

Reches, Ze'ev↗

Structure of axisymmetric mantle plumes

The structure of axisymmetric subsolidus thermal plumes in the earth's lower mantle is inferred from calculations of axisymmetric thermal plumes in an infinite Prandtl number fluid with thermally activated viscosity. The velocity and temperature distribution is determined for axisymmetric convection above a heated disk in an incompressible fluid cylinder 2,400 km in height and 1,200 km in diameter. Several calculations of plumes with heat transport in the range 100-400 GW, similar to the advective heat transport at the Hawaiian hotspot, are presented. Hotspot formation by plumes originating at the base of the mantle requires both large viscosity variations and a minimum heat transport.

Olson, Peter↗

Plume formation and lithosphere erosion - A comparison of laboratory and numerical experiments

The mechanics of thermal plume formation and intrusion into the lithosphere are investigated using a combination of laboratory and numerical simulations. The sequence of events leading to lithospheric thinning and uplift by thermal plumes is established, and some numerical estimates of the time scales for each stage in this process are derived that are applicable to the mantle. It is demonstrated that the two-dimensional finite element computations successfully reproduce the qualitative features seen in the experiments, with a quantitative discrepancy of typically 30 percent or less. The results of some calculations on plume formation and intrusion into model lithospheres with a variety of rheologies are presented.

Olson, Peter↗

Plume formation in the D-double prime-layer and the roughness of the core-mantle boundary

A series of very-high-resolution finite element calculations of plume formation in the D-double prime-layer has been performed for several plausible rheologies and boundary conditions in order to study both the early and later stages of boundary layer development. The results show that plumes are initiated by coalescence of small-scale convective instabilities within the low-viscosity region immediately above the core-mantle boundary (CMB). These instabilities support topographic roughness on the CMB having horizontal scales of 20-50 km and provide a source for scattered P-waves seen as precursors to the phases PKIKP and PKKP. The calculated structure of fully developed plumes emerging from the D-double prime-layer consists of 5-50 cm/yr flow confined to 50-100 km thick vertical conduits. With strongly temperature-dependent viscosity, plumes exhibit time-dependent behavior, including upward propagating solitary conduit waves, which may contribute to episodicity in hotspot volcanism.

Olson, Peter↗