magali i. billen u.c. davis dept. of geology margins successor program workshop, feb. 15-17, 2010...
TRANSCRIPT
![Page 1: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/1.jpg)
Magali I. Billen
U.C. DavisDept. of Geology
MARGINS Successor Program Workshop, Feb. 15-17, 2010
Geometry & Rates of 3D Mantle Flow in Subduction Zones
![Page 2: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/2.jpg)
MARGINS & Geodynamic Modeling
• Models of wedge convection– Rheology (deformation mechanisms, fabrics, LPO directions, dynamics)
– Fluids, petrology...– Mostly kinematic slabs & mostly 2D
![Page 3: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/3.jpg)
How will Geodynamics fit into a MARGINS Successor
Program?
1. Develop better tools for...– 3D & time-dependent models– Dynamic slabs (evolving trench & slab
geometry)– Coupling & tracking fluid & melt migration
flow– Understanding of special processes
• ie., subduction initiation, slab detachment, flat slabs...
We’re making progress here but it takes time to develop and test the required numerical methods.
![Page 4: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/4.jpg)
How will Geodynamics fit into a MARGINS Successor
Program?
2. Integrate modeling with all stages of MARGINS research– Guide deployment of seismic
stations, sample collection, etc...• Region specific models
– Analyze/interpret results from various focus sites• Generic (process-related) & regional
models
– Integrate & interpret multi-disciplinary observations
![Page 5: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/5.jpg)
Two Illustrative Examples
1. Ridge-Trench Interaction– PhD candidate Erin Burkett
2. 3D Mantle Flow at a Slab Edge– Margarete Jadamec (PhD 2009)
... illustrate two ways in which geodynamic modeling can be even better integrated into a MARGINS successor program.
![Page 6: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/6.jpg)
Ex. 1: Ridge-Trench Interaction
QuickTime™ and a decompressor
are needed to see this picture.
Burkett & Billen, JGR 2009
![Page 7: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/7.jpg)
Detachments & Plate Strength
• Detachment: integrated strength of subducted lithosphere => less than stress from sinking slab– plate age & rock yield strength.
![Page 8: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/8.jpg)
Regions With Slab Detachment?
• Costa Rica (continued sub.) & Baja Calif. (halted sub.)
• What are effects of 3D geometry?
![Page 9: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/9.jpg)
3D Ridge-Trench Interaction
Temperature
isosurface
Slabviscosityisosurface
ridgetrench
![Page 10: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/10.jpg)
3D Ridge-Trench Interaction
• Side view
• Front view
![Page 11: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/11.jpg)
3D Ridge-Trench Interaction
• Slab sinking induces complex 3D flow & interaction with approaching ridge & small-scale instabilities.
![Page 12: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/12.jpg)
Ex. 2: 3D Flow Models of Alaska
• Detailed regional model (2 km resolution).• Slab shape constructed from seismic observations.
![Page 13: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/13.jpg)
Geometry of 3D Flow at a Slab Edge
• Corner-flow dominates away from slab edge.• Slab is steepening (sinking back & down).• Toroidal flow around slab edge (slab-parallel flow).
![Page 14: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/14.jpg)
Decoupling of Plate & Mantle Flow
• Pacific plate motion matches observations. – Speed and direction.
• Mantle flows at rates of up to 90 cm/yr.– Slab-parallel component near slab edge ~ 10 cm/yr.
• Significant decoupling of mantle flow from plates.
![Page 15: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/15.jpg)
Evidence For Fast Mantle Flow
• Costa Rica: tracking isotopic signature transport along arc.– 6.5 - 19.0 cm/yr– Sub. Rate: 8.5 cm/yrHoernle et al., Nature 2008.
• If slab-parallel component is fraction (10 %) of mantle flow, predicts mantle flow rates of > 65 cm/yr
![Page 16: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/16.jpg)
ISA orientation, LPO & SKS Fast-Axis
• ISA can be non-parallel to mantle flow – wedge, slab edge. -- need B-type fabric in wedge nose.
• ISA match observations of SKS fast-axis orientations (from Christensen & Abers, 2009).
![Page 17: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/17.jpg)
ISA Sensitive to Rheology & Geometry
• Need broad (strategic) distribution of observations• Can distinguish successful models from unsuccessful
![Page 18: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/18.jpg)
3D Geometry of ISA Orientation
• Highly variable orientations in the mantle wedge: shallow horizontal, dipping slab-parallel, middle dipping and...
![Page 19: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/19.jpg)
3D Geometry of ISA Orientations
• Need: Better calculation of LPO from flow (A,B...)– 3D analysis of seismic anisotropy data & model results.
Slab-parallelstretching
![Page 20: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/20.jpg)
Conclusions
• Many opportunities to use dynamic modeling – to integrate observations & test hypothesis,– to help plan other experiments & observations.
• Need to create a strategy for development of better numerical methods for future MARGINS sceince.– What tools do we need most now?– How do we create these tool in tandem with collection & interpretation of data (field or laboratory-based)?
– How do we leverage work being done by CIG (Computational Infrastructure for Geodynamics)?
![Page 21: Magali I. Billen U.C. Davis Dept. of Geology MARGINS Successor Program Workshop, Feb. 15-17, 2010 Geometry & Rates of 3D Mantle Flow in Subduction Zones](https://reader034.vdocument.in/reader034/viewer/2022051416/56649e695503460f94b65da8/html5/thumbnails/21.jpg)