scispace - formally typeset
Search or ask a question

Showing papers on "Slab published in 2009"


Journal ArticleDOI
TL;DR: In this article, the authors developed numerical thermal models with a nonlinear mantle rheology for seventeen subduction zones, spanning a large range of slab age, descent rate, and geometry.
Abstract: Processes in subduction zones such as slab and mantle-wedge metamorphism, intraslab earthquakes, and arc volcanism vary systematically with the age-dependent thermal state of the subducting slab. In contrast, the configuration of subduction zones is rather uniform in that the arc is typically situated where the slab is ∼100 km deep. Toward reconciling the diversity and uniformity, we developed numerical thermal models with a nonlinear mantle rheology for seventeen subduction zones, spanning a large range of slab age, descent rate, and geometry. Where there are adequate observations, such as in Cascadia, northeast Japan, and Kamchatka, we find that surface heat flows can be explained if the interface between the slab and the mantle wedge is decoupled to a depth of 70–80 km. Models with this common decoupling depth predict that the region of high mantle temperatures and optimal fluid supply from the dehydrating slab, both required for melt generation for arc volcanism, occurs where the slab is ∼100 km deep. These models also reproduce the variations of the metamorphic, seismic, and volcanic processes with the thermal state of the slab. The shallow decoupling results in a stagnant fore arc whose thermal regime is controlled mainly by the subducting slab. The deeper coupling leads to a sudden onset of mantle wedge flow that brings heat from greater depths and the back arc, and its thermal effect overshadows that of the slab in the arc region. Our results serve to recast the research of subduction zone geodynamics into a quest for understanding what controls the common depth of decoupling.

388 citations


Journal ArticleDOI
TL;DR: A stagnant slab is a subducted slab of oceanic lithosphere subhorizontally deflected above, across, or below the 660 km discontinuity in the seafloor as mentioned in this paper.
Abstract: A stagnant slab is a subducted slab of oceanic lithosphere subhorizontally deflected above, across, or below the 660 km discontinuity. This phenomenon has now been widely recognized beneath subduction zones around the circum-Pacific and in the Mediterranean. Collaboration of seismic and electromagnetic observations, mineral physics measurements, and geodynamic modeling has begun to provide a consistent picture of stagnant slab.

326 citations


Journal ArticleDOI
TL;DR: In this article, a series of tests were conducted to investigate the blast resistances of slabs constructed with both plain ultra-high performance fiber concrete (UHPFC) and reinforced UHPFC, and slabs reinforced with externally bonded (EB) fibre reinforced polymer (FRP) plates.

303 citations


Journal ArticleDOI
22 Oct 2009-Nature
TL;DR: It is shown that the crystal-preferred orientation of serpentine, the main hydrous mineral in the upper mantle, can produce the strong trench-parallel seismic anisotropy observed in subduction systems.
Abstract: Seismic anisotropy is a powerful tool for detecting the geometry and style of deformation in the Earth's interior, as it primarily reflects the deformation-induced preferred orientation of anisotropic crystals. Although seismic anisotropy in the upper mantle is generally attributed to the crystal-preferred orientation of olivine, the strong trench-parallel anisotropy (delay time of one to two seconds) observed in several subduction systems is difficult to explain in terms of olivine anisotropy, even if the entire mantle wedge were to act as an anisotropic source. Here we show that the crystal-preferred orientation of serpentine, the main hydrous mineral in the upper mantle, can produce the strong trench-parallel seismic anisotropy observed in subduction systems. High-pressure deformation experiments reveal that the serpentine c-axis tends to rotate to an orientation normal to the shear plane during deformation; consequently, seismic velocity propagating normal to the shear plane (plate interface) is much slower than that in other directions. The seismic anisotropy estimated for deformed serpentine aggregates is an order of magnitude greater than that for olivine, and therefore the alignment of serpentine in the hydrated mantle wedge results in a strong trench-parallel seismic anisotropy in the case of a steeply subducting slab. This hypothesis is also consistent with the presence of a hydrous phase in the mantle wedge, as inferred from anomalously low seismic-wave velocities.

216 citations


Journal ArticleDOI
TL;DR: In this paper, the temperature of slab fluids indicate relatively hot conditions, and hint at a shallow and fast return path for ocean water, which is the case in many arc volcanism.
Abstract: Slab fluids drive mantle melting and return ocean water to the Earth's surface through arc volcanism. New ways of estimating the temperature of slab fluids indicate relatively hot conditions, and hint at a shallow and fast return path for ocean water.

211 citations


Journal ArticleDOI
TL;DR: Shirey et al. as mentioned in this paper showed that the major and trace element composition of sanukitoids is intermediate between typical Archaean TTG and modern arc granitoids, and two groups can be distinguished on the basis of the Ti content of the less differentiated rocks of the suite.
Abstract: A specific type of granitoid, referred to as sanukitoid (Shirey & Hanson 1984), was emplaced mainly across the Archaean–Proterozoic transition. The major and trace element composition of sanukitoids is intermediate between typical Archaean TTG and modern arc granitoids. However, among sanukitoids, two groups can be distinguished on the basis of the Ti content of the less differentiated rocks of the suite: high- and low-Ti sanukitoids. Melting experiments and petrogenetic modelling show that they may have formed by either (1) melting of mantle peridotite previously metasomatised by felsic melts of TTG composition, or (2) by reaction between TTG melts and mantle peridotite (assimilation). Rocks of the sanukitoid suite were emplaced at the Archaean–Proterozoic boundary, possibly marking the time when TTG-dominated granitoid magmatism changed to a more modern-style, arc-dominated magmatism. Consequently, the intermediate character of sanukitoids is not only compositional but chronological. The succession of granitoid magmatism with time is integrated in a plate tectonic model where it is linked to the thermal evolution of subduction zones, reflecting the progressive cooling of Earth: (1) the Archaean Earth’s heat production was high enough to allow the production of large amounts of TTG granitoids formed by partial melting of recycled basaltic crust (‘slab melting’); (2) at the end of the Archaean, due to the progressive cooling of the Earth, the extent of slab melting was reduced, resulting in lower melt:rock ratios. In such conditions the slab melts can be strongly contaminated by assimilation of mantle peridotite, thus giving rise to low-Ti sanukitoids. It is also possible that the slab melts were totally consumed in reactions with mantle peridotite, subsequent melting of this ‘melt-metasomatised mantle’ producing the high-Ti sanukitoid magmas; (3) after 2·5 Ga, Earth heat production was too low to allow slab melting, except in relatively rare geodynamic circumstances, and most modern arc magmas are produced by melting of the mantle wedge peridotite metasomatised by fluids from dehydration of the subducted slab. Of course, such changes did not take place exactly at the same time all over the world. The Archaean mechanisms coexisted with new processes over a relatively long time period, even if they were subordinate to the more modern processes.

201 citations


Journal ArticleDOI
TL;DR: The geometry of the Patagonia slab window is a subsurface tectonic feature resulting from subduction of the Nazca-Antarctic spreading-ridge system (Chile Rise) beneath southern South America as mentioned in this paper.

183 citations


Journal ArticleDOI
04 Jun 2009-Nature
TL;DR: It is demonstrated that the location of arc volcanoes is controlled by a combination of conditions: melting in the wedge is induced at the overlap of regions in the wedges that are hotter than the melting curve (solidus) of vapour-saturated peridotite and regions where hydrous minerals both inThe wedge and in the subducting slab break down.
Abstract: The processes that give rise to arc magmas at convergent plate margins have long been a subject of scientific research and debate. A consensus has developed that the mantle wedge overlying the subducting slab and fluids and/or melts from the subducting slab itself are involved in the melting process. However, the role of kinematic variables such as slab dip and convergence rate in the formation of arc magmas is still unclear. The depth to the top of the subducting slab beneath volcanic arcs, usually approximately 110 +/- 20 km, was previously thought to be constant among arcs. Recent studies revealed that the depth of intermediate-depth earthquakes underneath volcanic arcs, presumably marking the slab-wedge interface, varies systematically between approximately 60 and 173 km and correlates with slab dip and convergence rate. Water-rich magmas (over 4-6 wt% H(2)O) are found in subduction zones with very different subduction parameters, including those with a shallow-dipping slab (north Japan), or steeply dipping slab (Marianas). Here we propose a simple model to address how kinematic parameters of plate subduction relate to the location of mantle melting at subduction zones. We demonstrate that the location of arc volcanoes is controlled by a combination of conditions: melting in the wedge is induced at the overlap of regions in the wedge that are hotter than the melting curve (solidus) of vapour-saturated peridotite and regions where hydrous minerals both in the wedge and in the subducting slab break down. These two limits for melt generation, when combined with the kinematic parameters of slab dip and convergence rate, provide independent constraints on the thermal structure of the wedge and accurately predict the location of mantle wedge melting and the position of arc volcanoes.

170 citations


Journal ArticleDOI
TL;DR: In this article, travel-time tomography is carried out to estimate three-dimensional seismic velocity structures around Kanto using 735,520 P wave and 444,049 S wave arrival times from 6508 local earthquakes.
Abstract: [1] We first determine the configuration of the upper surface of the Pacific (PAC) slab beneath Kanto, Japan, from the distribution of interplate earthquakes relocated by an appropriate 1-D velocity model. Then, traveltime tomography is carried out to estimate three-dimensional seismic velocity structures around Kanto using 735,520 P wave and 444,049 S wave arrival times from 6508 local earthquakes. The obtained results suggest that the Philippine Sea (PHS) slab is subducting to depths of 130–140 km without a gap, even to the northwest of the Izu collision zone. We subsequently define the lateral extent of the contact zone between the bottom of the PHS slab and the upper surface of the PAC slab (PHS-PAC interface) and reveal that the slab contact zone underlies a wider area beneath Kanto in harmony with the Kanto plain. The downdip limit of interplate (thrust-type) earthquakes on the PAC slab is deepened by ∼30 km locally under the slab contact zone. This deepening is probably caused by a lower-temperature environment in the PAC slab, resulting from the overlap with the PHS slab subducting above and consequent thermal shielding by the PHS slab from the hot mantle wedge. We detect an extremely low-velocity anomaly in the easternmost portion of the PHS slab, which is probably attributable to serpentinization of mantle peridotite. Interplate earthquakes are almost absent along the PHS-PAC interface overlain by the serpentinized mantle in the PHS slab, suggesting that ductile deformation takes place along the interface because of low viscosity of the serpentine.

168 citations


Journal ArticleDOI
TL;DR: In this paper, an alternating-direction, implicit, finite-difference scheme is used to compute the thermal models taking into account all heating sources and phase boundaries, with the appropriate spreading rates and dip angles, compared with different island arc systems.
Abstract: Summary The thermal evolution of the lithospheric slab at subduction zones and its geophysical effects are numerically calculated. An alternating-direction, implicit, finite-difference scheme is used to compute the thermal models taking into account all heating sources and phase boundaries. These models, with the appropriate spreading rates and dip angles, are compared with different island arc systems. Temperatures inside the slab are strongly controlled by the conductivity and by the time elapsed since the initiation of descent. The depth to which temperature anomalies persist is generally about 700 km or less. The thermal results are used to construct the seismic velocities and ray paths, density anomalies, and the resulting stress distribution. Comparing the theoretical stress distribution and the focal mechanism studies of intermediate- and deep-focus earthquakes indicates the importance of both the mantle's rheology and the temperature dependence of the slab's elastic properties. The intermediate and deep focus earthquakes are located along the coolest region of the slab. The theoretical results explain the source mechanisms and the orientation of principal stresses under major island arcs.

164 citations


Journal ArticleDOI
TL;DR: In this article, the authors present new pieces of evidence from seismology and mineral physics for the existence of low-velocity zones in the deep part of the upper mantle wedge and the mantle transition zone that are caused by fluids from the deep subduction and deep dehydration of the Pacific and Philippine Sea slabs under western Pacific and East Asia.

Journal ArticleDOI
TL;DR: In this paper, the authors demonstrate that the timing, occurrence, and surface effects resulting from the detachment of subducted lithosphere are strongly dependent upon the stiffness (i.e. stress supported viscously) of the lithosphere and rheology of the upper mantle.

Journal ArticleDOI
TL;DR: In this article, the authors used several different plate motion models to describe plate and trench motions and evaluated the differences among the models and found that only one parameter, namely, the trench migration velocity in a Pacific hot spot reference frame, appears to correlate well with the strength of the subslab splitting signal.
Abstract: [1] The character of the mantle flow field in subduction zones remains poorly understood, despite its importance for our understanding of subduction dynamics. In particular, little attention has been paid to mantle flow beneath subducting slabs. In order to identify processes that make first-order contributions to the global pattern of subslab mantle flow, we have compiled shear wave splitting measurements from subduction zones worldwide from previously published studies and estimated average splitting parameters for the subslab region. Globally, the subslab region is overwhelmingly dominated by trench-parallel fast splitting directions. We tested for relationships between splitting delay time, a measure of the strength of anisotropy, and parameters that are indicators of tectonic processes, such as trench migration velocity, convergence velocity and obliquity, age of subducting lithosphere, and slab dip, curvature, seismicity, and thickness. We used several different plate motion models to describe plate and trench motions and evaluated the differences among the models. We find that only one parameter, namely, the trench migration velocity in a Pacific hot spot reference frame, appears to correlate well with the strength of the subslab splitting signal. This finding supports a model in which the mantle beneath subducting slabs is dominated by three-dimensional flow induced by trench migration. We explore several implications of our model for various aspects of mantle dynamics, including the choice of a suitable reference frame(s) for mantle flow, the existence of a thin decoupling zone between slabs and the subslab mantle, and consequences for mass transfer between the upper and lower mantle.

Journal ArticleDOI
TL;DR: In this paper, the authors investigated the bending strength of the downgoing plate in the subduction of a viscoelastic lithosphere driven solely by the upgoing plate's buoyancy in a passive mantle, represented by drag forces.
Abstract: [1] The bending strength of subducting lithosphere plays a critical role in the Earth's plate tectonics and mantle convection, modulating the amount of slab pull transmitted to the surface and setting the boundary conditions under which plates move and deform. However, it is the subject of a lively debate how much of the potential energy of the downgoing plate is consumed in bending the plate and how the lithospheric strength is defined during this process. We model the subduction of a viscoelastic lithosphere, driven solely by the downgoing plate's buoyancy, freely sinking in a passive mantle, represented by drag forces. To investigate the dynamics of bending, (1) we vary the density and the viscosity profile within the plate from isoviscous, where strength is distributed, to strongly layered, where strength is concentrated in a thin core, and (2) we map the stress, strain, and dissipation along the downgoing plate. The effective plate strength during bending is not a simple function of average plate viscosity but is affected by rheological layering and plate thinning. Earth-like layered plates allow for the transmission of large fractions of slab pull (∼75–80%) through the bend and yield a net slab pull of FSPnet = 1 to 6 × 1012 N m−1, which varies with the buoyancy of plates. In all models, only a minor portion of the energy is dissipated in the bending. Surprisingly, bending dissipation hardly varies with lithospheric viscosity because in our dynamic system, the plates minimize overall dissipation rate by adjusting their bending curvature. As a result, bending dissipation, ΦB, is mainly controlled by the bending moment work rate exerted by slab pull. We propose a new analytical formulation that includes a viscosity-dependent bending radius, which allows for assessment of the relative bending dissipation in the Earth's subduction zones using parameters from a recent global compilation. This yields estimates of ΦB/ΦTOT < 20%. These results suggest that plates on Earth weakly resist bending, yet are able to propagate a large amount of slab pull.

Journal ArticleDOI
TL;DR: In this paper, the authors describe field evidence for brittle hydrofracture in previously subducted rocks from the Western Alps, with a free aqueous fluid phase produced by dehydration reactions in the host blueschists and serpentinites.

Book ChapterDOI
TL;DR: Humphreys et al. as mentioned in this paper showed that the Laramide orogeny was caused by the combined effects of oceanic plateau subduction and unusually great suction in the mantle wedge.
Abstract: Flat subduction of the Farallon plate beneath the western United States during the Laramide orogeny was caused by the combined effects of oceanic plateau subduction and unusually great suction in the mantle wedge, the latter of which was caused by the shallowing slab approaching the North American craton root. Once in contact with basal North America, the slab cooled and hydrated the lithosphere. Upon removal, asthenospheric contact with lithosphere resulted in magma production that was especially intense where the basal lithosphere was fertile (in what now is the Basin and Range Province), and this heating weakened the lithosphere and made it convectively unstable. Small-scale convection has since affected many areas. With slab sinking and unloading, the western United States elevated into a broad plateau, and the weak portion gravitationally collapsed. With development of a transform plate boundary, the western part of the weak zone became entrained with the Pacifi c plate, and deformation there is dominated by shear. 2 Humphreys mwr204-04 1st pages E. Surface heat flow P-wave velocity D. Seismic velocity at 100 km Background: S-wave velocity/2 100-175 km C. S-wave velocity

Journal ArticleDOI
TL;DR: In this article, a hot rolling process of AA5083 aluminum alloy is simulated using the finite element method, where the temperature distribution in the roll and the slab, the stress, strain and strain rate fields, are extracted throughout a steady-state analysis of the process.

Journal ArticleDOI
TL;DR: In this article, the authors focus on the dynamic effects on the basis of 2D thermomechanical modelling incorporating effects of slab dehydration, mantle-wedge melting and surface topography development, and conclude that subduction of a moderate-size ridge does not typically result in strong slab flattening and related decrease of magmatic activity.
Abstract: The subduction of oceanic aseismic ridges, oceanic plateaus and seamount chains is a common process that takes place in a variety of tectonic settings and seems to coincide spatially and temporally with a gap of volcanic activity, shallow or even horizontal slab angles, enhanced seismic activity and various topographic features. In the present study we focus on these dynamic effects on the basis of 2D thermomechanical modelling incorporating effects of slab dehydration, mantle-wedge melting and surface topography development. In order to ascertain the impact of a moderate-size (200 × 18 km) aseismic ridge, 12 pairs of experiments (one for the case with a ridge, the other without) were carried out varying slab density and subducting- and overriding-plate velocities. By analysing pairs of experiments we conclude that subduction of a moderate-sized ridge does not typically result in strong slab flattening and related decrease of magmatic activity. This, in turn, suggests that, when slab flattening is indeed associated with the ridge subduction in nature, the slab itself should be in a nearly critical ( i.e ., transient from inclined to flat) state so that any local addition of positive buoyancy may strongly affect overall slab dynamics. Therefore, subducting ridges may serve as indicators of transient slab states in nature. Another important result from our study is the numerical quantification of strongly decreased magma production associated with flat slabs that may explain gaps in recent active volcanism at low-angle subduction margins. Lowering of magmatic rock production is caused by the absence of a hot mantle wedge above the flat slabs and does not directly depend on the mechanism responsible for the triggering of slab flattening. Finally we document several very distinct surface effects associated with the moderate-size ridge subduction such as local increase in elevation of overriding margin, enhancement of subduction erosion and landward trench displacement. Surface uplift may exceed the original ridge height due to additional uplift resulting from the overriding plate shortening. Topographic perturbations within the accretionary wedge domain are transient and have a tendency to relax after the ridge passes the trench. In contrast, the topographic high created in the continental portion of the overriding plate relaxes more slowly and may even be sustained for several millions of year after the ridge subduction.

Journal ArticleDOI
TL;DR: In this paper, the environmental impacts of two types of slab systems are studied and the assessment shows that the environmental impact of a building structure with precast concrete floors is 12.2% lower than that with in situ cast floors for the defined functional unit.

Journal ArticleDOI
TL;DR: In this article, the authors use 3-D petrological-thermomechanical numerical simulations to investigate small-scale convection and hydrous, partially molten, cold plumes formed in the mantle wedge in response to slab dehydration.
Abstract: Hydration and partial melting along subducting slabs can trigger Rayleigh-Taylor-like instabilities. We use 3-D petrological-thermomechanical numerical simulations to investigate small-scale convection and hydrous, partially molten, cold plumes formed in the mantle wedge in response to slab dehydration. The simulations were carried out with the I3ELVIS code, which is based on a multigrid approach combined with marker-in-cell methods and conservative finite difference schemes. Our numerical simulations show that three types of plumes occur above the slab-mantle interface: (1) finger-like plumes that form sheet-like structure parallel to the trench, (2) ridge-like structures perpendicular to the trench, and (3) flattened wave-like instabilities propagating upward along the upper surface of the slab and forming zigzag patterns parallel to the trench. The viscosity of the plume material is the main factor controlling the geometry of the plumes. Our results show that lower viscosity of the partially molten rocks facilitates the Rayleigh-Taylor-like instabilities with small wavelengths. In particular, in low-viscosity models (1018–1019 Pa s) the typical spacing of finger-like plumes is about 30–45 km, while in high-viscosity models (1020–1021 Pa s) plumes become rather sheet-like, and the spacing between them increases to 70–100 km. Water released from the slab forms a low-viscosity wedge with complex 3-D geometries. The computed spatial and temporal pattern of melt generation intensity above the slab is compared to the distribution and ages of volcanoes in the northeast Japan. Based on the similarity of the patterns we suggest that specific clustering of volcanic activity in this region could be potentially related to the activity of thermal-chemical plumes.

Journal ArticleDOI
TL;DR: In this article, the authors review geophysical and geological data from 23 active and ancient plate boundaries to document a compelling spatial and temporal relationship between collision, plate boundary curvature, rapid tectonic rotations, and the occurrence of back-arc rifting.
Abstract: [1] We review geophysical and geological data from 23 active and ancient plate boundaries to document a compelling spatial and temporal relationship between collision, plate boundary curvature, rapid tectonic rotations, and the occurrence of back-arc rifting Our observations support a conceptual model where the change from subduction to collision causes rapid fore-arc rotation (when viewed in a reference frame relative to the bounding tectonic plates), leading to marked plate boundary curvature Our global compilation reveals that most active back-arc rifts are associated with rapidly rotating fore arcs and nearby collisions Thus, we propose that collision and rapid fore-arc rotations play a major role in the evolution and kinematics of back-arc basins We conduct numerical modeling to better understand the physical processes behind these observed rapid tectonic rotations at convergent margins Our results suggest that the presence of an indentor or choke point in the subduction system (eg, collision) can generate rapid rotation of the fore arc about a nearby pole and lead to back-arc rifting The rate of fore-arc rotation and back-arc rifting depends on the incoming indentor velocity and can be greatly enhanced by slab rollback and the presence of a low-viscosity back arc Where viscosity of the back arc is low, fore-arc rotation dominates; where back-arc viscosity is high, the formation of strike-slip faults and tectonic escape dominates Our observational and model-derived results illustrate that shallow crustal forces produced by the entry of buoyant features into subduction zones are a fundamental mechanism for the generation of fore-arc block rotations, plate boundary curvature, back-arc rift evolution, and tectonic escape in subduction systems Rollback of the subducting slab within the mantle will certainly enhance the processes we describe but it is not the only driving force

Journal ArticleDOI
TL;DR: In this article, a petrological model (pHMELTS) coupled with a 2D thermal and variable viscosity flow model (ConMan) is used to describe and compare fundamental processes occurring within subduction zones.

Patent
13 Feb 2009
TL;DR: In this article, a planar slab of piezoelectric material having a top surface is provided and a plurality of holes are made through the top surface and into the slab.
Abstract: A method comprises providing a substantially planar slab of piezoelectric material having a top surface. The method further comprises drilling a plurality of holes through the top surface and into the slab. The method further comprises making a plurality of cuts through the top surface and into the slab. The cuts form a plurality of polygons that are generally centered about one of the holes. The method further comprises plating the slab with an electrically conductive material. The method further comprises removing the electrically conductive material from the top surface of the slab. The method further comprises cutting the slab substantially parallel to the top surface.

Journal ArticleDOI
TL;DR: The Arc Basalt Simulator (ABS2) as mentioned in this paper uses an Excel® spreadsheet-based calculator to predict the partitioning of incompatible element and Sr-Nd-Pb isotopic composition in a slab-derived fluid and in arc basalt magma generated by an open system fluid-fluxed melting of mantle wedge peridotite.
Abstract: [1] Convergent margin magmas typically have geochemical signatures that include elevated concentrations of large-ion lithophile elements; depleted heavy rare earth elements and high field strength elements; and variously radiogenic Sr, Pb, and Nd isotopic compositions. These have been attributed to the melting of depleted mantle peridotite by the fluxing of fluids or melts derived from subducting oceanic crust. High Mg # basalts and high Mg # andesites are inferred to make up the bulk of subduction-related primary magmas and may be generated by fluid or melt fluxing of mantle peridotite. The difference in contributions from the subducted slab found among various arcs appears to be mostly controlled by thermal structure. Cold slabs yield fluids, and hot slabs yield melts. Recent experimental studies and thermodynamic models better constrain the phase petrology of the slab components during prograde metamorphism and melting, mantle wedge melting, and mantle slab melt reaction. Experimental results also constrain the behavior of many elements in these processes. In addition, geodynamic models allow increasingly realistic, quantitative modeling of the temperature and pressure in the subducted slab and mantle wedge. These developments together enable generation of forward models to explain arc magma geochemistry. The Arc Basalt Simulator (ABS) version 2 (ABS2) uses an Excel® spreadsheet-based calculator to predict the partitioning of incompatible element and Sr-Nd-Pb isotopic composition in a slab-derived fluid and in arc basalt magma generated by an open system fluid-fluxed melting of mantle wedge peridotite. The ABS2 model is intended to simulate high Mg # basalt geochemistry in relatively cold subduction zones. The modeling scheme of ABS2 is presented and is applied to primitive arc magmas.

Journal ArticleDOI
TL;DR: In this paper, the authors present a small-scale laboratory model of oceanic subduction in which plates motion is imposed by lateral boundary conditions and the long-term process of subduction is not steady-state.

Journal ArticleDOI
TL;DR: In this paper, the authors calculated the mineralogical properties of the Earth's mantle and the lithological units constituting the subducting oceanic slabs within a wide range of P-T conditions within the CaO-FeO-MgO-Al2O3-SiO2 system, except for the basalt top-layer of a slab.

Journal ArticleDOI
29 May 2009-Science
TL;DR: It is shown that the Izu-Bonin and Japan slabs are torn apart at their junction when they bend to flatten over the 660-kilometer discontinuity, as is expected from a simple geometric argument.
Abstract: Subducted slabs of oceanic lithosphere below the western Pacific tend to be stagnant in the transition zone with poorly known mechanical properties Typical examples are the Izu-Bonin and Japan slabs that meet each other to form a cusplike junction beneath southwest Japan Here, we show that these two slabs are torn apart at their junction when they bend to flatten over the 660-kilometer discontinuity, as is expected from a simple geometric argument We present three lines of evidence for this ongoing slab tear

Journal ArticleDOI
TL;DR: In this article, a new method for modelling floating-slab tracks with discontinuous slabs in underground railway tunnels is presented, where the track is subjected to a harmonic load moving with a constant velocity.

Book ChapterDOI
01 Jan 2009
TL;DR: In this paper, the authors review models of subduction dynamics that have been studied over the last decade by means of numerical and analog experiments and provide a better understanding of slab strength which is still uncertain since large-scale subduction zone observables and laboratory results do not put strong constraints on slab rheology.
Abstract: Subduction of oceanic lithosphere and deep slabs control several aspects of plate tectonics. We review models of subduction dynamics that have been studied over the last decade by means of numerical and analog experiments. Regional models indicate that trench rollback, trench curvature, and back-arc deformation may be explained by fl uid slabs that are 250–500 times stiffer than the upper mantle. Slab width and, more importantly, rheology determine the role of viscous bending, poloidal-sinking fl ow and toroidal-rollback stirring, and interactions of the slab with the higher viscosity lower mantle. Several of these contributions can be represented by a local sinking veloCity. Back-arc deformation may then result from an imbalance if larger-scale plate forcing leads to deviations of the convergence rate from the local equilibrium. Lateral viscosity variations (LVVs) are also key for understanding plate driving forces. The realism of global circulation computations has advanced and such models with weak zones and other LVVs have lead to an improved match to observed plate tectonic scores. Those include the correlation with plate motions, the magnitude of intraplate deformation, and oceanic to continental plate veloCity ratios. Net rotation of the lithosphere with respect to the lower mantle may be caused jointly by regional slab forcing and the stirring effect of cratonic keels. However, slab models have so far only produced net rotations that are small compared to recent hotspot reference-frame models. Progress in the next years will likely come from a better understanding of slab strength, which is still uncertain since large-scale subduction zone observables and laboratory results do not put strong constraints on slab rheology. Importantly, circulation models with an improved representation of convergent margins will help to close the gap between regional and global approaches to subduction, and to better understand the potential role of the overriding plate.

Journal ArticleDOI
TL;DR: A regularization procedure to define electrostatic energies and forces in a slab system of thickness h that is periodic in two dimensions and carries a net charge is introduced and it is shown that the method does indeed provide consistent long-range electrostatics.
Abstract: We introduce a regularization procedure to define electrostatic energies and forces in a slab system of thickness h that is periodic in two dimensions and carries a net charge. The regularization corresponds to a neutralization of the system by two charged walls and can be viewed as the extension to the two-dimensional (2D)+h geometry of the neutralization by a homogeneous background in the standard three-dimensional Ewald method. The energies and forces can be computed efficiently by using advanced methods for systems with 2D periodicity, such as MMM2D or P3M/ELC, or by introducing a simple background-charge correction to the Yeh-Berkowitz approach of slab systems. The results are checked against direct lattice sum calculations on simple systems. We show, in particular, that the Madelung energy of a 2D square charge lattice in a uniform compensating background is correctly reproduced to high accuracy. A molecular dynamics simulation of a sodium ion close to an air/water interface is performed to demonstrate that the method does indeed provide consistent long-range electrostatics. The mean force on the ion reduces at large distances to the image-charge interaction predicted by macroscopic electrostatics. This result is used to determine precisely the position of the macroscopic dielectric interface with respect to the true molecular surface.