Siberian flood basalt magmatism and Mongolia-Okhotsk slab dehydration

Experimental data combined with numerical calculations suggest that fast subducting slabs are cold enough to carry into the deep mantle a significant portion of the water in antigorite, which transforms with increasing depth to phase A and then to phase E and/or wadsleyite by solid-solid phase transition. Clathrate hydrates and ice VII are also stable at PT conditions of cold slabs and represent other potential phases for water transport into the deep mantle. Some cold slabs are expected to deflect while crossing the 410 km and stagnate in transition zone being unable to penetrate through 660 km discontinuity. In this way slabs can move a long way beneath continents after long-lived subduction. With time, the stagnant slabs are heated to the temperature of the ambient transition zone and release free H~2~O-bearing fluid. Combining with transition zone water filter model this may cause voluminous melting of overlying upper mantle rocks. If such process operates in nature, magmas geochemically similar to island-arc magmas are expected to appear in places relatively remote from active arcs at the time of their emplacement. Dolerites of the south-eastern margin of the Siberian flood basalt province, located about 700 km from suggested trench, were probably associated with fast subduction of the Mongolia-Okhotsk slab and originated by dehydration of the stagnant slab in the transition zone. We show that influence of the subduction-related deep water cycle on Siberian flood basalt magmatism gradually reduced with increasing distance from the subduction zone

Sulfur loss from subducted altered oceanic crust and implications for mantle oxidation

© The Author(s), [year]. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Walters, J. B., Cruz-Uribe, A. M., & Marschall, H. R. Sulfur loss from subducted altered oceanic crust and implications for mantle oxidation. Geochemical Perspectives Letters, 13, (2020): 36-41, doi:10.7185/geochemlet.2011.Oxygen fugacity (fO2) is a controlling factor of the physics of Earth’s mantle; however, the mechanisms driving spatial and secular changes in fO2 associated with convergent margins are highly debated. We present new thermodynamic models and petrographic observations to predict that oxidised sulfur species are produced during the subduction of altered oceanic crust. Sulfur loss from the subducting slab is a function of the protolith Fe3+/ΣFe ratio and subduction zone thermal structure, with elevated sulfur fluxes predicted for oxidised slabs in cold subduction zones. We also predict bi-modal release of sulfur-bearing fluids, with a low volume shallow flux of reduced sulfur followed by an enhanced deep flux of sulfate and sulfite species, consistent with oxidised arc magmas and associated copper porphyry deposits. The variable SOx release predicted by our models both across and among active margins may introduce fO2 heterogeneity to the upper mantle.We thank James Connolly for modelling support and Peter van Keken for providing updated P–T paths for the Syracuse et al. (2010) models. The manuscript benefited from the editorial handling by Helen Williams and from constructive reviews of Maryjo Brounce, Katy Evans, and an anonymous reviewer. JBW acknowledges Fulbright and Chase Distinguished Research fellowships. This work was supported by NSF grant EAR1725301 awarded to AMC

Water, oceanic fracture zones and the lubrication of subducting plate boundaries - insights from seismicity

We investigate the relationship between subduction processes and related seismicity for the Lesser Antilles Arc using the Gutenberg-Richter law. This power lawdescribes the earthquakemagnitude distribution, with the gradient of the cumulative magnitude distribution being commonly known as the b-value. The Lesser Antilles Arc was chosen because of its alongstrike variability in sediment subduction and the transition from subduction to strike-slip movement towards its northern and southern ends. The data are derived from the seismicity catalogues from the Seismic Research Centre of The University of the West Indies and the Observatoires Volcanologiques et Sismologiques of the Institut de Physique du Globe de Paris and consist of subcrustal events primarily from the slab interface. The b-value is found using a Kolmogorov-Smirnov test for a maximum-likelihood straight line-fitting routine. We investigate spatial variations in b-values using a grid-search with circular cells as well as an along-arc projection. Tests with different algorithms and the two independent earthquake cataloges provide confidence in the robustness of our results. We observe a strong spatial variability of the b-value that cannot be explained by the uncertainties. Rather than obtaining a simple north-south b-value distribution suggestive of the dominant control on earthquake triggering being water released from the sedimentary cover on the incoming American Plates, or a b-value distribution that correlates with on the obliquity of subduction, we obtain a series of discrete, high b-value 'bull's-eyes' along strike. These bull's-eyes, which indicate stress release through a higher fraction of small earthquakes, coincide with the locations of known incoming oceanic fracture zones on the American Plates. We interpret the results in terms of water being delivered to the Lesser Antilles subduction zone in the vicinity of fracture zones providing lubrication and thus changing the character of the related seismicity. Our results suggest serpentinization around mid-ocean ridge transform faults, which go on to become fracture zones on the incoming plate, plays a significant role in the delivery of water into the mantle at subduction zones

Origins of chemical diversity of back-arc basin basalts: a segment-scale study of the Eastern Lau Spreading Center

We report major, trace, and volatile element data on basaltic glasses from the northernmost segment of the Eastern Lau Spreading Center (ELSC1) in the Lau back-arc basin to further test and constrain models of back-arc volcanism. The zero-age samples come from 47 precisely collected stations from an 85 km length spreading center. The chemical data covary similarly to other back-arc systems but with tighter correlations and well-developed spatial systematics. We confirm a correlation between volatile content and apparent extent of melting of the mantle source but also show that the data cannot be reproduced by the model of isobaric addition of water that has been broadly applied to back-arc basins. The new data also confirm that there is no relationship between mantle temperature and the wet melting productivity. Two distinct magmatic provinces can be identified along the ELSC1 axis, a southern province influenced by a “wet component” with strong affinities to arc volcanism and a northern province influenced by a “damp component” intermediate between enriched mid-ocean ridge basalts (E-MORB) and arc basalts. High–field strength elements and rare earth elements are all mobilized to some extent by the wet component, and the detailed composition of this component is determined. It differs in significant ways from the Mariana component reported by E. Stolper and S. Newman (1994), particularly by having lower abundances of most elements relative to H_(2)O. The differences can be explained if the slab temperature is higher for the Mariana and the source from which the fluid is derived is more enriched. The ELSC1 damp component is best explained by mixing between the wet component and an E-MORB-like component. We propose that mixing between water-rich fluids and low-degree silicate melts occurs at depth in the subduction zone to generate the chemical diversity of the ELSC1 subduction components. These modified sources then rise independently to the surface and melt, and these melts mix with melts of the background mantle from the ridge melting regime to generate the linear data arrays characteristic of back-arc basalts. The major and trace element framework for ELSC1, combined with different slab temperatures and compositions for difference convergent margins, may be able to be applied to other back-arc basins around the globe

Two-sided asymmetric subduction; implications for tectonomagmatic and metallogenic evolution of the Lut Block, Eastern Iran

West directed subduction zones show common characteristics, such as low structural elevation, deep trench, steep slab and a conjugate back-arc basin that are opposite to those of the east directed subduction zones. The tectonomagmatic and metallogenic setting of the Lut Block is still a matter of debate and several hypotheses have been put forward. Despite some authors denying the influence of the operation of Benioff planes, the majority propose that it occurred beneath the Afghan Block, while others consider that oceanic lithosphere was dragged under the Lut Block. Cu-Au porphyry deposits seem to occur in an island arc geotectonic setting during the middle Eocene while Mo-bearing deposits are coincident with the crustal thickening during Oligocene. We introduce new trace element and isotope geochemical data for granitoids and structural evidences testifying the two-sided asymmetric subduction beneath both Afghan and Lut Blocks, with different rates of consumption of oceanic lithosphere

Boninite and Harzburgite from Leg 125 (Bonin-Mariana Forearc): A Case Study of Magma Genesis during the Initial Stages of Subduction

Holes drilled into the volcanic and ultrabasic basement of the Izu-Ogasawara and Mariana forearc terranes during Leg 125 provide data on some of the earliest lithosphere created after the start of Eocene subduction in the Western Pacific. The volcanic basement contains three boninite series and one tholeiite series. (1) Eocene low-Ca boninite and low-Ca bronzite andesite pillow lavas and dikes dominate the lowermost part of the deep crustal section through the outer-arc high at Site 786. (2) Eocene intermediate-Ca boninite and its fractionation products (bronzite andesite, andesite, dacite, and rhyolite) make up the main part of the boninitic edifice at Site 786. (3) Early Oligocene intermediate-Ca to high-Ca boninite sills or dikes intrude the edifice and perhaps feed an uppermost breccia unit at Site 786. (4) Eocene or Early Oligocene tholeiitic andesite, dacite, and rhyolite form the uppermost part of the outer-arc high at Site 782. All four groups can be explained by remelting above a subduction zone of oceanic mantle lithosphere that has been depleted by its previous episode of partial melting at an ocean ridge. We estimate that the average boninite source had lost 10-15 wt% of melt at the ridge before undergoing further melting (5-10%) shortly after subduction started. The composition of the harzburgite (<2% clinopyroxene, Fo content of about 92%) indicates that it underwent a total of about 25% melting with respect to a fertile MORB mantle. The low concentration of Nb in the boninite indicates that the oceanic lithosphere prior to subduction was not enriched by any asthenospheric (OIB) component. The subduction component is characterized by (1) high Zr and Hf contents relative to Sm, Ti, Y, and middle-heavy REE, (2) light REE-enrichment, (3) low contents of Nb and Ta relative to Th, Rb, or La, (4) high contents of Na and Al, and (5) Pb isotopes on the Northern Hemisphere Reference Line. This component is unlike any subduction component from active arc volcanoes in the Izu-Mariana region or elsewhere. Modeling suggests that these characteristics fit a trondhjemitic melt from slab fusion in amphibolite facies. The resulting metasomatized mantle may have contained about 0.15 wt% water. The overall melting regime is constrained by experimental data to shallow depths and high temperatures (1250°C and 1.5 kb for an average boninite) of boninite segregation. We thus envisage that boninites were generated by decompression melting of a diapir of metasomatized residual MORB mantle leaving the harzburgites as the uppermost, most depleted residue from this second stage of melting. Thermal constraints require that both subducted lithosphere and overlying oceanic lithosphere of the mantle wedge be very young at the time of boninite genesis. This conclusion is consistent with models in which an active transform fault offsetting two ridge axes is placed under compression or transpression following the Eocene plate reorganization in the Pacific. Comparison between Leg 125 boninites and boninites and related rocks elsewhere in the Western Pacific highlights large regional differences in petrogenesis in terms of mantle mineralogy, degree of partial melting, composition of subduction components, and the nature of pre-subduction lithosphere. It is likely that, on a regional scale, the initiation of subduction involved subducted crust and lithospheric mantle wedge of a range of ages and compositions, as might be expected in this type of tectonic setting

Segmented Hellenic slab rollback driving Aegean deformation and seismicity

The NE dipping slab of the Hellenic subduction is imaged in unprecedented detail using teleseismic receiver function analysis on a dense 2-D seismic array. Mapping of slab geometry for over 300 km along strike and down to 100 km depth reveals a segmentation into dipping panels by along-dip faults. Resolved intermediate-depth seismicity commonly attributed to dehydration embrittlement is shown to be clustered along these faults. Large earthquakes occurrence within the upper and lower plate and at the interplate megathrust boundary show a striking correlation with the slab faults suggesting high mechanical coupling between the two plates. Our results imply that the general slab rollback occurs here in a differential piecewise manner imposing its specific stress and deformation pattern onto the overriding Aegean plate

Flat slab deformation caused by interplate suction force

We image the structure at the southern end of the Peruvian flat subduction zone, using receiver function and surface wave methods. The Nazca slab subducts to ~100 km depth and then remains flat for ~300 km distance before it resumes the dipping subduction. The flat slab closely follows the topography of the continental Moho above, indicating a strong suction force between the slab and the overriding plate. A high-velocity mantle wedge exists above the initial half of the flat slab, and the velocity resumes to normal values before the slab steepens again, indicating the resumption of dehydration and ecologitization. Two prominent midcrust structures are revealed in the 70 km thick crust under the Central Andes: molten rocks beneath the Western Cordillera and the underthrusting Brazilian Shield beneath the Eastern Cordillera