Return to search

Two-dimensional, viscous flow modeling of roll-back subduction : numerical investigation into the role of slab density in subduction dynamics

Thesis: S.B., Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences, 2004. / Cataloged from PDF version of thesis. / Includes bibliographical references (pages 37-38). / Observations of retreating subduction systems in the Mediterranean region suggest the density of subducting lithosphere is dynamically related to trench retreat rate and upper-plate deformation. Most numerical and analog studies of retreating subduction systems have not explored the effects of lithospheric density variations on subduction processes. This study is a preliminary effort to construct a two-dimensional, viscous flow model of "roll-back" subduction to explicitly examine how slab density influences retreat rate, mantle flow, and slab geometry. For a given lithosphere-mantle density contrast, the model computes the evolution of a viscous, thermal slab using a finite element code for incompressible convection (ConMan). Imposed velocity boundary conditions guide lithospheric material into a uniformly weak "subduction zone" and out into the mantle below, generating stable, asymmetric subduction. Slabs driven faster than the "intrinsic" (dynamically consistent), steady-state retreat rate of the system (vr) are characteristically arcuate, pushed upward from the base of the mantle layer by strong horizontal "return flow" beneath the descending lithosphere. Slabs driven slower than vr are sigmoidal: the slabs steepen at depth, where vertical buoyancy forces overcome lateral viscous forces set up by weak surface velocities. The diagnostic behaviors of slabs driven faster and slower than vr define a set of qualitative criteria (slab geometry, mantle flow patterns) for converging on the consistent, steady-state retreat rate of the system. For slab-mantle density contrast [delta]p = 198 kg/m³ (defined as the density difference between lithosphere at surface of the system and mantle material at the base of the system), vr ~~ 16 mm/yr. The slab is roughly planar, with 500 dip. For [delta]p = 168 kg/m³, vr is slightly slower (14 mm/yr), and steady-state slab geometry is nearly identical (moderately-dipping planar surface). It is found that the angle at which lithospheric material is forced into the mantle does not significantly affect either steady-state retreat rate or slab geometry. / by Jessica L. Haurin. / S.B.

Identiferoai:union.ndltd.org:MIT/oai:dspace.mit.edu:1721.1/114112
Date January 2004
CreatorsHaurin, Jessica L. (Jessica Lyn)
ContributorsLeigh Royden., Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences., Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences.
PublisherMassachusetts Institute of Technology
Source SetsM.I.T. Theses and Dissertation
LanguageEnglish
Detected LanguageEnglish
TypeThesis
Format38 pages, application/pdf
RightsMIT theses are protected by copyright. They may be viewed, downloaded, or printed from this source but further reproduction or distribution in any format is prohibited without written permission., http://dspace.mit.edu/handle/1721.1/7582

Page generated in 0.0021 seconds