Showing papers on "Slab published in 2016"

Philippine Sea and East Asian plate tectonics since 52 Ma constrained by new subducted slab reconstruction methods

Abstract: We reconstructed Philippine Sea and East Asian plate tectonics since 52 Ma from 28 slabs mapped in 3-D from global tomography, with a subducted area of ~25% of present-day global oceanic lithosphere Slab constraints include subducted parts of existing Pacific, Indian, and Philippine Sea oceans, plus wholly subducted proto-South China Sea and newly discovered “East Asian Sea” Mapped slabs were unfolded and restored to the Earth surface using three methodologies and input to globally consistent plate reconstructions Important constraints include the following: (1) the Ryukyu slab is ~1000 km N-S, too short to account for ~20° Philippine Sea northward motion from paleolatitudes; (2) the Marianas-Pacific subduction zone was at its present location (±200 km) since 48 ± 10 Ma based on a >1000 km deep slab wall; (3) the 8000 × 2500 km East Asian Sea existed between the Pacific and Indian Oceans at 52 Ma based on lower mantle flat slabs; (4) the Caroline back-arc basin moved with the Pacific, based on the overlapping, coeval Caroline hot spot track These new constraints allow two classes of Philippine Sea plate models, which we compared to paleomagnetic and geologic data Our preferred model involves Philippine Sea nucleation above the Manus plume (0°/150°E) near the Pacific-East Asian Sea plate boundary Large Philippine Sea westward motion and post-40 Ma maximum 80° clockwise rotation accompanied late Eocene-Oligocene collision with the Caroline/Pacific plate The Philippine Sea moved northward post-25 Ma over the northern East Asian Sea, forming a northern Philippine Sea arc that collided with the SW Japan-Ryukyu margin in the Miocene (~20–14 Ma)

Philippine Sea and East Asian plate tectonics since 52 Ma constrained by new subducted slab reconstruction methods

Abstract: We reconstructed Philippine Sea and East Asian plate tectonics since 52 Ma from 28 slabs mapped in 3-D from global tomography, with a subducted area of ~25% of present-day global oceanic lithosphere. Slab constraints include subducted parts of existing Pacific, Indian, and Philippine Sea oceans, plus wholly subducted proto-South China Sea and newly discovered “East Asian Sea.” Mapped slabs were unfolded and restored to the Earth surface using three methodologies and input to globally consistent plate reconstructions. Important constraints include the following: (1) the Ryukyu slab is ~1000 km N-S, too short to account for ~20° Philippine Sea northward motion from paleolatitudes; (2) the Marianas-Pacific subduction zone was at its present location (±200 km) since 48 ± 10 Ma based on a >1000 km deep slab wall; (3) the 8000 × 2500 km East Asian Sea existed between the Pacific and Indian Oceans at 52 Ma based on lower mantle flat slabs; (4) the Caroline back-arc basin moved with the Pacific, based on the overlapping, coeval Caroline hot spot track. These new constraints allow two classes of Philippine Sea plate models, which we compared to paleomagnetic and geologic data. Our preferred model involves Philippine Sea nucleation above the Manus plume (0°/150°E) near the Pacific-East Asian Sea plate boundary. Large Philippine Sea westward motion and post-40 Ma maximum 80° clockwise rotation accompanied late Eocene-Oligocene collision with the Caroline/Pacific plate. The Philippine Sea moved northward post-25 Ma over the northern East Asian Sea, forming a northern Philippine Sea arc that collided with the SW Japan-Ryukyu margin in the Miocene (~20–14 Ma).

Continuity of the Alpine slab unraveled by high-resolution P wave tomography

Abstract: The question of lateral and/or vertical continuity of subducted slabs in active orogens is a hot topic partly due to poorly resolved tomographic data. The complex slab structure beneath the Alpine region is only partly resolved by available geophysical data, leaving many geological and geodynamical issues widely open. Based upon a finite-frequency kernel method, we present a new high-resolution tomography model using P wave data from 527 broadband seismic stations, both from permanent networks and temporary experiments. This model provides an improved image of the slab structure in the Alpine region and fundamental pinpoints for the analysis of Cenozoic magmatism, (U)HP metamorphism, and Alpine topography. Our results document the lateral continuity of the European slab from the Western Alps to the central Alps, and the downdip slab continuity beneath the central Alps, ruling out the hypothesis of slab break off to explain Cenozoic Alpine magmatism. A low-velocity anomaly is observed in the upper mantle beneath the core of the Western Alps, pointing to dynamic topography effects. A NE dipping Adriatic slab, consistent with Dinaric subduction, is possibly observed beneath the Eastern Alps, whereas the laterally continuous Adriatic slab of the Northern Apennines shows major gaps at the boundary with the Southern Apennines and becomes near vertical in the Alps-Apennines transition zone. Tear faults accommodating opposite-dipping subductions during Alpine convergence may represent reactivated lithospheric faults inherited from Tethyan extension. Our results suggest that the interpretations of previous tomography results that include successive slab break offs along the Alpine-Zagros-Himalaya orogenic belt might be proficiently reconsidered.

Fresnel coefficients of a two-dimensional atomic crystal

Abstract: In general the experiments on the linear optical properties of a single-layer two-dimensional atomic crystal are interpreted by modeling it as a homogeneous slab with an effective thickness. Here I fit the most remarkable experiments in graphene optics by using the Fresnel coefficients, fixing both the surface susceptibility and the surface conductivity of graphene. It is shown that the Fresnel coefficients and the slab model are not equivalent. Experiments indicate that the Fresnel coefficients are able to simulate the overall experiments here analyzed, while the slab model fails to predict absorption and the phase of the reflected light.

Snow fracture in relation to slab avalanche release: critical state for the onset of crack propagation

Abstract: . The failure of a weak snow layer buried below cohesive slab layers is a necessary, but insufficient, condition for the release of a dry-snow slab avalanche. The size of the crack in the weak layer must also exceed a critical length to propagate across a slope. In contrast to pioneering shear-based approaches, recent developments account for weak layer collapse and allow for better explaining typical observations of remote triggering from low-angle terrain. However, these new models predict a critical length for crack propagation that is almost independent of slope angle, a rather surprising and counterintuitive result. Based on discrete element simulations we propose a new analytical expression for the critical crack length. This new model reconciles past approaches by considering for the first time the complex interplay between slab elasticity and the mechanical behavior of the weak layer including its structural collapse. The crack begins to propagate when the stress induced by slab loading and deformation at the crack tip exceeds the limit given by the failure envelope of the weak layer. The model can reproduce crack propagation on low-angle terrain and the decrease in critical length with increasing slope angle as modeled in numerical experiments. The good agreement of our new model with extensive field data and the ease of implementation in the snow cover model SNOWPACK opens a promising prospect for improving avalanche forecasting.

Zinc isotope evidence for sulfate-rich fluid transfer across subduction zones

Abstract: Subduction zones modulate the chemical evolution of the Earth's mantle. Water and volatile elements in the slab are released as fluids into the mantle wedge and this process is widely considered to result in the oxidation of the sub-arc mantle. However, the chemical composition and speciation of these fluids, which is critical for the mobility of economically important elements, remain poorly constrained. Sulfur has the potential to act both as oxidizing agent and transport medium. Here we use zinc stable isotopes (δ66Zn) in subducted Alpine serpentinites to decipher the chemical properties of slab-derived fluids. We show that the progressive decrease in δ66Zn with metamorphic grade is correlated with a decrease in sulfur content. As existing theoretical work predicts that Zn-SO42- complexes preferentially incorporate heavy δ66Zn, our results provide strong evidence for the release of oxidized, sulfate-rich, slab serpentinite-derived fluids to the mantle wedge.

Ground‐Motion Prediction Equations for Subduction Slab Earthquakes in Japan Using Site Class and Simple Geometric Attenuation Functions

Abstract: The frequency content of strong ground motions from subduction slab earthquakes differs significantly from that of ground motions produced by other categories (tectonic locations: shallow crustal, upper mantle, and subduction interface) of earthquakes in subduction zones. In the last two decades, a large number of records from subduction slab events have been obtained in Japan. We present a ground‐motion prediction equation (GMPE) for this category of earthquakes. We used a large dataset from reliably identified slab events up to the end of 2012. The GMPEs were based on a set of simple geometric attenuation functions, site classes were used as site terms, and nonlinear site amplification ratios were adopted. A bilinear magnitude‐scaling function was adopted for large earthquakes with moment magnitude M w≥7.1, with the scaling rates for large events being much smaller than for the smaller events. A magnitude‐squared term was used for events with M w<7.1 as well as the bilinear magnitude‐scaling function. We also modeled the effect of volcanic zones using an anelastic attenuation coefficient applied to a horizontal portion of the seismic‐wave travel distance within possible volcanic zones. We found that excluding the records from sites with inferred site classes improved the model goodness of fit. The within‐event residuals were approximately separated into within‐site and between‐site residuals, and the corresponding standard deviations were calculated using a random effects model. The separation of within‐event residuals into within‐site and between‐site components allows for the possibility of adopting different standard deviations for different site classes in a probabilistic seismic‐hazard analysis if desired. Online Material: Figures showing the distribution of between‐event residuals with respect to magnitude and fault‐top depth and the distribution of within‐event residuals with respect to magnitude and source distance.

Magnesium isotope geochemistry in arc volcanism

Abstract: Incorporation of subducted slab in arc volcanism plays an important role in producing the geochemical and isotopic variations in arc lavas. The mechanism and process by which the slab materials are incorporated, however, are still uncertain. Here, we report, to our knowledge, the first set of Mg isotopic data for a suite of arc lava samples from Martinique Island in the Lesser Antilles arc, which displays one of the most extreme geochemical and isotopic ranges, although the origin of this variability is still highly debated. We find the δ26Mg of the Martinique Island lavas varies from −0.25 to −0.10, in contrast to the narrow range that characterizes the mantle (−0.25 ± 0.04, 2 SD). These high δ26Mg values suggest the incorporation of isotopically heavy Mg from the subducted slab. The large contrast in MgO content between peridotite, basalt, and sediment makes direct mixing between sediment and peridotite, or assimilation by arc crust sediment, unlikely to be the main mechanism to modify Mg isotopes. Instead, the heavy Mg isotopic signature of the Martinique arc lavas requires that the overall composition of the mantle wedge is buffered and modified by the preferential addition of heavy Mg isotopes from fluids released from the altered subducted slab during fluid−mantle interaction. This, in turn, suggests transfer of a large amount of fluid-mobile elements from the subducting slab to the mantle wedge and makes Mg isotopes an excellent tracer of deep fluid migration.

Mantle transition zone structure beneath the Changbai volcano: Insight into deep slab dehydration and hot upwelling near the 410 km discontinuity

Abstract: We study the detailed mantle transition zone structure beneath the active Changbai intraplate volcano in Northeast China by using a receiver-function method. A total of 3005 teleseismic receiver functions recorded by 70 broadband stations are obtained by using a common-conversion-point stacking method. For conducting the time-to-depth conversion, we use a three-dimensional velocity model of the study region so as to take into account the influence of structural heterogeneities. Our results reveal significant depth variations of the 410, 520, and 660 km discontinuities. A broad depression of the 410 km discontinuity and a low-velocity anomaly are revealed beneath the Changbai volcano, which may reflect a large-scale hot mantle upwelling around the 410 km discontinuity with a positive Clapeyron slope. The 520 km discontinuity is identified clearly, and its uplift occurs above the stagnant Pacific slab. We also find a prominent depression of the 660 km discontinuity, which is elongated along the trend of deep earthquake clusters in a range of 39°N–44°N latitude, and the depression area has a lateral extent of about 400 km. Because the 520 and 660 km discontinuities correspond to positive and negative Clapeyron slopes, respectively, we think that the 520 uplift and the 660 depression are caused by the cold subducting Pacific slab. A part of the Pacific slab may have penetrated into the lower mantle and so caused the large-scale 660 depression in front of the deep earthquake clusters. Our results also reveal a part of the upper boundary of the subducting Pacific slab in the mantle transition zone.

Are Adakites Slab Melts or High-pressure Fractionated Mantle Melts?

Abstract: Adakites are unusual felsic igneous rocks commonly associated with asthenospheric slab window opening or fast subduction of a young (<25 Ma) oceanic plate that may allow slab melting at shallow depths. Their genesis has been extensively debated, as they are also observed in other geodynamic settings where thermal models do not predict slab melting in the fore-arc region. Here, we present a new approach that provides new constraints on adakite petrogenesis in hot subduction zones (e.g. the Philippines) and above an asthenospheric window (e.g. Baja California, Mexico). We use amphibole compositions to estimate magma storage depths and the composition of the host melts to test the hypothesis that adakites are pristine slab melts. We find that adakites from the Philippines and Baja California fore-arcs formed in two distinct petrogenetic scenarios: in the Philippines, water-rich mantle melts stalled and crystallized within lower and upper crustal magma storage regions to produce silica-rich melts with an adakitic signature; in Baja California, slab melts that percolated through the mantle wedge mixed or mingled with water-rich mantle melts within a lower crustal (∼30 km depth) magma storage region before stalling in the upper arc crust (∼7–15 km depth). Alternatively, the Baja California adakites may represent mixing products between high-pressure differentiated mantle melts and mantle melts in a lower crustal magma reservoir, periodically refluxed by mantle melts. Thus, slab melting is not necessarily required to produce an adakitic geochemical fingerprint in hot subduction zones. The hot downgoing plate may cross the ‘adakitic window’ and melt in specific geodynamic settings such as the opening of a slab tear, as beneath Baja California.

Constraints from GPS measurements on the dynamics of deformation in Anatolia and the Aegean

Abstract: We estimate the strength of the lithosphere in Anatolia and the Aegean, and the boundary forces acting upon it, using a dynamical model that treats the lithosphere as a thin fluid sheet deforming in response to variations in gravitational potential energy. This model has one free material parameter, the power law exponent, n, of the vertically averaged rheology of the lithosphere, and two parameters that specify the forces per unit length applied to its edges. Solutions to this model that best fit the velocities of 346 reliable GPS sites require an effective viscosity of the lithosphere of 1022 to 1021 Pa s at strain rates of 10 to 100 nanostrain per year. The best-fitting force at the Arabia-Anatolia boundary is consistent with the lithostatic pressure due to the high topography there, and the force at the Nubia-Aegean boundary is consistent with the contrast in lithostatic pressure across that boundary. No additional force, from “slab rollback” or basal tractions due to convection in the mantle, is required to explain the observations. These results are supported by scaling relations derived from approximate analytical solutions. The inverse relationship between the viscosity of the lithosphere and deviatoric stress produces strong slowly deforming regions in the Southern Aegean and Central Anatolia whose motions resemble those of microplates. The distribution of geodetic strain rates within the region, and the partitioning between normal and strike-slip faulting, are explained by the interplay between boundary conditions, internal variations in gravitational potential energy, and the power law rheology of the lithosphere.

Origins of ultralow velocity zones through slab-derived metallic melt

Abstract: Understanding the ultralow velocity zones (ULVZs) places constraints on the chemical composition and thermal structure of deep Earth and provides critical information on the dynamics of large-scale mantle convection, but their origin has remained enigmatic for decades. Recent studies suggest that metallic iron and carbon are produced in subducted slabs when they sink beyond a depth of 250 km. Here we show that the eutectic melting curve of the iron−carbon system crosses the current geotherm near Earth’s core−mantle boundary, suggesting that dense metallic melt may form in the lowermost mantle. If concentrated into isolated patches, such melt could produce the seismically observed density and velocity features of ULVZs. Depending on the wetting behavior of the metallic melt, the resultant ULVZs may be short-lived domains that are replenished or regenerated through subduction, or long-lasting regions containing both metallic and silicate melts. Slab-derived metallic melt may produce another type of ULVZ that escapes core sequestration by reacting with the mantle to form iron-rich postbridgmanite or ferropericlase. The hypotheses connect peculiar features near Earth's core−mantle boundary to subduction of the oceanic lithosphere through the deep carbon cycle.

Response of precast foamed concrete sandwich panels to flexural loading

Abstract: This paper presents the results of an experimental and analytical investigation of a total of six precast foamed concrete sandwich panels (PFCSPs) as one-way acting slabs tested under flexural loads. Foamed concrete of 25.73 MPa was used to produce the PFCSP concrete wythes. The results obtained from the tests have been discussed in terms of ultimate flexural strength capacity, moment-vertical deflection profile, load–strain relationship, strain variation across the slab depth, influence of aspect ratio, cracking patterns, and ultimate flexural load at failure. An analytical study of finite element analysis (FEA) as a one-way slab model was then conducted. The increase in aspect ratio (L/d) from 18.33 to 26.67 shows a reduction of 50% and 69.6% on the ultimate flexural strength capacity as obtained experimentally and in FEA models, respectively. Theoretical analyses on the extremes of fully composite and non-composite actions were also determined. The experimental results showed that cracking patterns were observed in one direction only, similar to those reported on a reinforced concrete solid slab, as well as precast concrete sandwich panels, when both concrete wythes act in a single composite manner. The experimental results were compared with FEA model data, and a significant degree of accuracy was obtained. Therefore, the PFCSP slab can serve as an alternative to the normal concrete slab system in buildings.

Global correlation of lower mantle structure and past subduction

Abstract: Advances in global seismic tomography have increasingly motivated identification of subducted lithosphere in Earth's deep mantle, creating novel opportunities to link plate tectonics and mantle evolution. Chief among those is the quest for a robust subduction reference frame, wherein the mantle assemblage of subducted lithosphere is used to reconstruct past surface tectonics in an absolute framework anchored in the deep Earth. However, the associations heretofore drawn between lower mantle structure and past subduction have been qualitative and conflicting, so the very assumption of a correlation has yet to be quantitatively corroborated. Here we show that a significant, time-depth progressive correlation can be drawn between reconstructed subduction zones of the last 130 Myr and positive S wave velocity anomalies at 600-2300 km depth, but that further correlation between greater times and depths is not presently demonstrable. This correlation suggests that lower mantle slab sinking rates average between 1.1 and 1.9 cmyr(-1).

Full-waveform inversion of the Japanese Islands region

Abstract: We present a full-waveform tomographic model of the crust and upper mantle beneath the Japanese Islands region. This is based on the combination of GPU-accelerated spectral-element wavefield simulations, adjoint techniques, and nonlinear optimization. Our model explains complete seismic waveforms of events not used in the inversion in the period range from 20 to 80 s. Quantitative resolution analysis indicates that resolution lengths within the well-covered areas are around 150 km in the horizontal and around 30 km in the vertical directions. In addition to the high-velocity signatures of known lithospheric slabs in the region, our model reveals a pronounced low-velocity anomaly beneath the volcanic island of Ulleung in the Sea of Japan, reaching −19% around 100 km depth. The Ulleung anomaly originates at or above the Pacific slab, rises vertically upward to the base of the Philippine Sea slab at ∼200 km depth, circumvents it in NW direction, and then significantly strengthens in the uppermost mantle above the Philippine Sea slab. Among the numerous hypotheses for the generation of low-velocity anomalies in subduction systems, those invoking instabilities before or when a slab enters the transition zone seem most likely. The age and fast subduction of the Pacific slab may facilitate the transport of fluids into the transition zone. This may promote the reduction in viscosity and the onset of convective upwelling, aided by ambient mantle flow, such as return flow within the mantle wedge.

Seismic Simulation of an Integrated Ceiling-Partition Wall-Piping System at E-Defense. I: Three-Dimensional Structural Response and Base Isolation

Abstract: The seismic response of a full-scale, 5-story steel moment frame building in base-isolated and fixed-base configurations with an integrated suspended ceiling-partition wall-sprinkler piping system that was shaken at E-Defense is critically assessed. Horizontal floor accelerations were constrained by the isolation systems to relatively low levels, which allowed observation of damage to the integrated system that was directly related to the vertical component of input acceleration. The floor slabs exhibited single mode vibration at their natural periods with widely varying effective damping. Peak vertical accelerations were amplified by an average factor ranging from 3 to 6 from the table to the middle of the floor slabs, at which amplification factors increased as slab vibration periods lengthened. Damage to the ceiling-partion-piping components initiated at slab accelerations of approximately 2 g and became extensive for slab accelerations exceeding 5 g. These metrics establish target vertical acc...

Estimating the effective elastic modulus and specific fracture energy of snowpack layers from field experiments

Abstract: Measurements of the mechanical properties of snow are essential for improving our understanding and the prediction of snow failure and hence avalanche release. We performed fracture mechanical experiments in which a crack was initiated by a saw in a weak snow layer underlying cohesive snow slab layers. Using particle tracking velocimetry (PTV), the displacement field of the slab was determined and used to derive the mechanical energy of the system as a function of crack length. By fitting the estimates of mechanical energy to an analytical expression, we determined the slab effective elastic modulus and weak layer specific fracture energy for 80 different snowpack combinations, including persistent and nonpersistent weak snow layers. The effective elastic modulus of the slab ranged from 0.08 to 34 MPa and increased with mean slab density following a power-law relationship. The weak layer specific fracture energy ranged from 0.08 to 2.7 J m−2 and increased with overburden. While the values obtained for the effective elastic modulus of the slab agree with previously published low-frequency laboratory measurements over the entire density range, the values of the weak layer specific fracture energy are in some cases unrealistically high as they exceeded those of ice. We attribute this discrepancy to the fact that our linear elastic approach does not account for energy dissipation due to non-linear parts of the deformation in the slab and/or weak layer, which would undoubtedly decrease the amount of strain energy available for crack propagation.

Imaging the transition from flat to normal subduction: variations in the structure of the Nazca slab and upper mantle under southern Peru and northwestern Bolivia

Abstract: S U M M A R Y Two arrays of broad-band seismic stations were deployed in the north central Andes between 8◦ and 21◦S, the CAUGHT array over the normally subducting slab in northwestern Bolivia and southern Peru, and the PULSE array over the southern part of the Peruvian flat slab where the Nazca Ridge is subducting under South America. We apply finite frequency teleseismic Pand S-wave tomography to data from these arrays to investigate the subducting Nazca plate and the surrounding mantle in this region where the subduction angle changes from flat north of 14◦S to normally dipping in the south. We present new constraints on the location and geometry of the Nazca slab under southern Peru and northwestern Bolivia from 95 to 660 km depth. Our tomographic images show that the Peruvian flat slab extends further inland than previously proposed along the projection of the Nazca Ridge. Once the slab re-steepens inboard of the flat slab region, the Nazca slab dips very steeply (∼70◦) from about 150 km depth to 410 km depth. Below this the slab thickens and deforms in the mantle transition zone. We tentatively propose a ridge-parallel slab tear along the north edge of the Nazca Ridge between 130 and 350 km depth based on the offset between the slab anomaly north of the ridge and the location of the re-steepened Nazca slab inboard of the flat slab region, although additional work is needed to confirm the existence of this feature. The subslab mantle directly below the inboard projection of the Nazca Ridge is characterized by a prominent low-velocity anomaly. South of the Peruvian flat slab, fast anomalies are imaged in an area confined to the Eastern Cordillera and bounded to the east by well-resolved low-velocity anomalies. These low-velocity anomalies at depths greater than 100 km suggest that thick mantle lithosphere associated with underthrusting of cratonic crust from the east is not present. In northwestern Bolivia a vertically elongated fast anomaly under the Subandean Zone is interpreted as a block of delaminating lithosphere.