{"id":{"repo_id":"mit","oai_identifier":"oai:dspace.mit.edu:1721.1/114112"},"canonical_url":"https://search.dev.ndltd.org/etd/mit/oai:dspace.mit.edu:1721.1/114112","repository":{"repo_id":"mit","name":"MIT","base_url":"https://dspace.mit.edu/oai/request"},"display":{"title":"Two-dimensional, viscous flow modeling of roll-back subduction : numerical investigation into the role of slab density in subduction dynamics","abstract":"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.","abstract_html":"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.","abstract_has_math":false,"creators":["Haurin, Jessica L. (Jessica Lyn)"],"institution":"Massachusetts Institute of Technology","degree_name":null,"degree_level":null,"degree_discipline":null,"degree_department":"Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences.","school":null,"contributors":[],"advisors":["Leigh Royden."],"committee_chairs":[],"committee_members":[],"year":2004,"date_issued":"2004","date_published":"2004","updated_at":"2026-07-22T22:21:55Z","subjects":["Earth, Atmospheric, and Planetary Sciences."],"languages":["eng"],"rights":["MIT 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."],"rights_urls":["http://dspace.mit.edu/handle/1721.1/7582"],"identifier_entries":[]},"links":{"outbound_url":"http://hdl.handle.net/1721.1/114112","outbound_label":"Handle","outbound_source":"dc:identifier.uri"},"metadata_groups":[{"id":"people","label":"People","entries":[{"key":"dc:contributor.advisor","label":"Advisor","values":["Leigh Royden."]},{"key":"dc:contributor.department","label":"Department","values":["Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences."]},{"key":"dc:contributor.other","label":"Dc Contributor Other","values":["Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences."]},{"key":"dc:creator","label":"Author","values":["Haurin, Jessica L. (Jessica Lyn)"]}]},{"id":"academic_context","label":"Academic Context","entries":[{"key":"dc:date.accessioned","label":"Dc Date Accessioned","values":["2018-03-12T19:30:23Z"]},{"key":"dc:date.available","label":"Dc Date Available","values":["2018-03-12T19:30:23Z"]},{"key":"dc:date.issued","label":"Date","values":["2004"]},{"key":"dc:publisher","label":"Institution","values":["Massachusetts Institute of Technology"]},{"key":"dc:type","label":"Dc Type","values":["Thesis"]}]},{"id":"subjects_keywords","label":"Subjects and Keywords","entries":[{"key":"dc:subject","label":"Dc Subject","values":["Earth, Atmospheric, and Planetary Sciences."]}]},{"id":"language_rights","label":"Language and Rights","entries":[{"key":"dc:language.iso","label":"Language (ISO)","values":["eng"]},{"key":"dc:rights","label":"Dc Rights","values":["MIT 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."]},{"key":"dc:rights.uri","label":"Rights URI","values":["http://dspace.mit.edu/handle/1721.1/7582"]}]},{"id":"identifiers","label":"Identifiers","entries":[{"key":"dc:identifier.uri","label":"Identifier URI","values":["http://hdl.handle.net/1721.1/114112"]}]},{"id":"additional","label":"Additional Metadata","entries":[{"key":"dc:description","label":"Description","values":["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)."]},{"key":"dc:description.abstract","label":"Abstract","values":["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."]},{"key":"dc:description.degree","label":"Dc Description Degree","values":["S.B."]},{"key":"dc:title","label":"Title","values":["Two-dimensional, viscous flow modeling of roll-back subduction : numerical investigation into the role of slab density in subduction dynamics"]}]}],"canonical_facts":{"dc:contributor.advisor":["Leigh Royden."],"dc:contributor.department":["Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences."],"dc:contributor.other":["Massachusetts Institute of Technology. Department of Earth, Atmospheric, and Planetary Sciences."],"dc:creator":["Haurin, Jessica L. (Jessica Lyn)"],"dc:date.accessioned":["2018-03-12T19:30:23Z"],"dc:date.available":["2018-03-12T19:30:23Z"],"dc:date.issued":["2004"],"dc:description":["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)."],"dc:description.abstract":["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."],"dc:description.degree":["S.B."],"dc:identifier.uri":["http://hdl.handle.net/1721.1/114112"],"dc:language.iso":["eng"],"dc:publisher":["Massachusetts Institute of Technology"],"dc:rights":["MIT 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."],"dc:rights.uri":["http://dspace.mit.edu/handle/1721.1/7582"],"dc:subject":["Earth, Atmospheric, and Planetary Sciences."],"dc:title":["Two-dimensional, viscous flow modeling of roll-back subduction : numerical investigation into the role of slab density in subduction dynamics"],"dc:type":["Thesis"]},"updated_at":"2026-07-22T22:21:55Z"}