Rheologic controls on slab dynamics

Author Posting. © American Geophysical Union, 2007. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 8 (2007): Q08012, doi:10.1029/2007GC001597.Several models have been proposed to relate slab geometry to parameters such as plate velocity or plate age. However, studies on the observed relationships between slab geometry and a wide range of subduction parameters show that there is not a simple global relationship between slab geometry and any one of these other subduction parameters for all subduction zones. Numerical and laboratory models of subduction provide a method to explore the relative importance of different physical processes in determining subduction dynamics. Employing 2-D numerical models with a viscosity structure constrained by laboratory experiments for the deformation of olivine, we show that the observed range in slab dip and the observed trends between slab dip and convergence velocity, subducting plate age, and subduction duration can be reproduced without trench motion (i.e., slab roll-back) for locations away from slab edges. Successful models include a stiff slab that is 100–1000 times more viscous than previous estimates from models of plate bending, the geoid, and global plate motions. We find that slab dip in the upper mantle depends primarily on slab strength and plate boundary coupling, with a small dependence on subducting plate age. Once the slab sinks into the lower mantle the primary processes controlling slab evolution are (1) the ability of the stiff slab to transmit stresses up dip, (2) resistance to slab descent into the higher-viscosity lower mantle, and (3) subduction-induced flow in the mantle-wedge corner.This research was partially supported by NSF award EAR0125919

Morphology and origin of the Osbourn Trough

The Osbourn Trough is a 900 km long, east-west trending gravity low, visible in satellite altimetry data, which intersects the Tonga Trench at 25°30′S. A recent survey collected gravity, magnetic, echo sounder, and swath bathymetry data on three adjacent, north-south trending ship tracks centered on the trough. The linear gravity low is 20–30 mGal less than the regional value and is accompanied by a flat-lying, 200–500 m deep sediment-filled valley. Swath bathymetry images reveal several parallel, east-west trending linear ridges and valleys on either side of the main trough as well as other morphologic features characteristic of relict spreading centers, including a prominent inside corner high and possible pseudo-fault trace. The presence of magnetic anomalies (possibly anomalies 33 and 32) suggests that the seafloor here was formed after the end of the Cretaceous Normal Superchron (anomaly 34). These data support the conclusion that this trough is a spreading center, which stopped spreading in Late Cretaceous time. The existence of this feature has important implications for tectonic reconstructions in this region. The Osbourn Trough could be part of the fossil spreading center between the Pacific Plate and a fragment of the Phoenix Plate, the Bellingshausen Plate. This would require the seafloor between the Osbourn Trough and the Chatham Rise to the south to be a remnant fragment of the Bellingshausen Plate. This remnant may have joined to the Pacific Plate when the Hikurangi Plateau entered the Gondwana subduction zone at the Chatham Rise possibly causing the cessation of spreading on the Osbourn Trough

Young solid Earth researchers of the world unite!

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/95272/1/eost14667.pd

Newtonian versus non-Newtonian upper mantle viscosity : implications for subduction initiation

Author Posting. © American Geophysical Union, 2005. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 32 (2005): L19304, doi:10.1029/2005GL023457.The effect of rheology on the evolution of the slab-tip during subduction initiation is analyzed using 2-D numerical flow models. Experimentally determined flow laws have both strong temperature- and stress-dependence, which leads to large local variations in viscosity with direct consequences for subduction initiation. We find that models with Newtonian viscosity lead to flat or coupled subduction due to hydrodynamic stresses that pull the slab-tip up towards the overriding plate. Non-Newtonian rheology reduces these hydrodynamic stresses by decreasing the wedge viscosity and the slab coupling to wedge-corner flow, rendering the small negative-slab buoyancy of the slab-tip sufficient to maintain its dip during the early stages of subduction

Influence of cratonic lithosphere on the formation and evolution of flat slabs : insights from 3-D time-dependent modeling.

Several mechanisms have been suggested for the formation of flat slabs including buoyant features on the subducting plate, trenchward motion and thermal or cratonic structure of the overriding plate. Analysis of episodes of flat subduction indicate that not all flat slabs can be attributed to only one of these mechanisms and it is likely that multiple mechanisms work together to create the necessary conditions for flat slab subduction. In this study we examine the role of localized regions of cratonic lithosphere in the overriding plate in the formation and evolution of flat slabs. We explicitly build on previous models, by using time-dependent simulations with three-dimensional variation in overriding plate structure. We find that there are two modes of flat subduction: permanent underplating occurs when the slab is more buoyant (shorter or younger), while transient flattening occurs when there is more negative buoyancy (longer or older slabs). Our models show how regions of the slab adjacent to the subcratonic flat portion continue to pull the slab into the mantle leading to highly contorted slab shapes with apparent slab gaps beneath the craton. These results show how the interpretation of seismic images of subduction zones can be complicated by the occurrence of either permanent or transient flattening of the slab, and how the signature of a recent flat slab episode may persist as the slab resumes normal subduction. Our models suggest that permanent underplating of slabs may preferentially occur below thick and cold lithosphere providing a built-in mechanism for regeneration of cratons

I. Seafloor Morphology of the Osbourn Trough and Kermadec Trench and II. Multiscale DynamiCS of Subduction Zones

This thesis aims to demonstrate that integration of detailed observations of deformation at short to long length scales with carefully formulated numerical modeling is an effective method for simulating the complex multiscale nature of mantle-lithosphere dynamics. In Part I, marine geophysical observations are used to determine the origin of the Osbourn Trough, a long linear depression within the Pacific Plate seaward of the Tonga-Kermadec Trench, and to determine the elastic strength of the subducting plate within the Kermadec Trench. Based on the morphology of the seafloor from swath bathymetry mapping and modeling of magnetic data, we conclude that the Osbourn Trough is an extinct spreading center which stopped spreading about 72 million years ago. Swath bathymetry mapping within the Kermadec Trench reveals extensive faulting within the trench on the subducting plate, with oblique grabens aligned perpendicular to the absolute plate motion direction. Using isostatic flexural response methods, we find that the flexural rigidity (1e19-1e20 Nm) is smaller than normally found for old oceanic lithosphere reflecting a local reduction in the strength of the plate. In Part II, regional 3-D dynamic models of the Tonga-Kermadec and Aleutian subduction zones are used to constrain lateral variations in viscosity in the upper mantle. Modeling of the dynamic topography of the overriding plate for the Tonga-Kermadec subduction zone requires a low viscosity and low density (-20 kg/m³) region within the wedge above the slab to decouple the slab-induced flow from. These efforts lead to a good fit to the observed shallow bathymetry on the overriding plate for a model with a slab density anomaly due to temperature of ~80 kg/m³. However, the geoid anomaly above the subduction zone is too large by 20-40 m at length scales of 100-1000 km. A reduction of the slab density by a factor of 1.5 is needed to match both the geoid and topography, suggesting the density anomaly of the slab due to temperature is compensated within the upper mantle (~100-300 km). Similar modeling for the Aleutians including a narrower low viscosity region and smaller density anomaly (-10 kg/m³) in the wedge is able to fit the geoid and topography without reducing the slab density.</p

Deep slab seismicity limited by rate of deformation in the transition zone.

Deep earthquakes within subducting tectonic plates (slabs) are enigmatic because they appear similar to shallow earthquakes but must occur by a different mechanism. Previous attempts to explain the depth distribution of deep earthquakes in terms of the temperature at which possible triggering mechanisms are viable, fail to explain the spatial variability in seismicity. In addition to thermal constraints, proposed failure mechanisms for deep earthquakes all require that sufficient strain accumulates in the slab at a relatively high stress. Here, I show that simulations of subduction with nonlinear rheology and compositionally dependent phase transitions exhibit strongly variable strain rates in space and time, which is similar to observed seismicity. Therefore, in addition to temperature, variations in strain rate may explain why there are large gaps in deep seismicity (low strain rate), and variable peaks in seismicity (bending regions), and, possibly, why there is an abrupt cessation of seismicity below 660 km