Geochemical evidence for the melting of subducting oceanic lithosphere at plate edges

Vertical tearing of subducting plates controlled by geometry and rheology of oceanic plates

Lateral non-uniform subduction is impacted by continuous plate segmentation owing to vertical tearing of the subducting plate. However, the dynamics and physical controls of vertical tearing remain controversial. Here, we employed 3D numerical models to investigate the effects of trench geometry (offset by a transform boundary) and plate rheology (plate age and the magnitude of brittle/plastic strain weakening) on the evolution of shear stress-controlled vertical tearing within a homogenous subducting oceanic plate. Numerical results suggest that the trench offset geometry could result in self-sustained vertical tearing as a narrow shear zone within the intact subducting oceanic plate, and that this process of tearing could operate throughout the entire subduction process. Further, the critical trench offset length for the maturation of vertical tearing is impacted by plate rheology. Comparison between numerical modelling results and natural observations suggests that vertical tearing attributed to trench offset geometry is broadly developed in modern subduction and collision systems worldwide.

Climatic control on the location of continental volcanic arcs

Orogens and volcanic arcs at continental plate margins are primary surface expressions of convergent plate tectonics. Although it is established that climate affects the shape, size, and architecture of orogens via orographic erosion gradients, the ascent of magma through the crust and location of volcanoes along magmatic arcs have been considered insensitive to erosion. However, available data reveal westward migration of late-Cenozoic volcanic activity in the Southern Andes and Cascade Range where orography drives an eastward migration of the topographic water divide by increased precipitation and erosion along west-facing slopes. Thermomechanical numerical modeling shows that orographic erosion and the associated leeward topographic migration may entail asymmetric crustal structures that drive the magma ascent toward the region of enhanced erosion. Despite the different tectonic histories of the Southern Andes and the Cascade Range, orographic erosion is a shared causal mechanism that can explain the late-Cenozoic westward migration of the volcanic front along both magmatic arcs.

Sources of dehydration fluids underneath the Kamchatka arc

Fluids mediate the transport of subducted slab material and play a crucial role in the generation of arc magmas. However, the source of subduction-derived fluids remains debated. The Kamchatka arc is an ideal subduction zone to identify the source of fluids because the arc magmas are comparably mafic, their source appears to be essentially free of subducted sediment-derived components, and subducted Hawaii-Emperor Seamount Chain (HESC) is thought to contribute a substantial fluid flux to the Kamchatka magmas. Here we show that Tl isotope ratios are unique tracers of HESC contribution to Kamchatka arc magma sources. In conjunction with trace element ratios and literature data, we trace the progressive dehydration and melting of subducted HESC across the Kamchatka arc. In succession, serpentine (<100 km depth), lawsonite (100–250 km depth) and phengite (>250 km depth) break down and produce fluids that contribute to arc magmatism at the Eastern Volcanic Front (EVF), Central Kamchatka Depression (CKD), and Sredinny Ridge (SR), respectively. However, given the Tl-poor nature of serpentine and lawsonite fluids, simultaneous melting of subducted HESC is required to explain the HESC-like Tl isotope signatures observed in EVF and CKD lavas. In the absence of eclogitic crust melting processes in this region of the Kamchatka arc, we propose that progressive dehydration and melting of a HESC-dominated mélange offers the most compelling interpretation of the combined isotope and trace element data.

Fluids released from progressive breakdown of minerals at increasing pressure within a mélange may explain the trace element systematics and stable thallium isotope data of the Kamchatka arc lavas from volcanic front to back arc.

Petrogenesis of isotopically enriched Quaternary magma with adakitic affinity associated with subduction of old lithosphere beneath central Myanmar

We present a model for the petrogenesis of magma with adakitic affinity in an old subduction zone, which does not involve slab melting and is constrained by new geochronological and geochemical data for Mt. Popa, the largest of three Quaternary volcanoes in central Myanmar (Popa, Monywa and Singu). The edifice is composed of Popa Plateau (0.8–0.6 Ma) with high-K rocks and a stratovolcano (< 0.33 Ma) predominantly composed of medium-K rocks with adakitic affinity (Mg# 45–63, Sr/Y > 40). The distinct K contents indicate that the adakitic magmas cannot be derived from Popa high-K rocks, but they share trace-element signatures and Sr–Nd isotope ratios with medium-K basalts from Monywa volcano. Our estimation of water contents in Popa magma reveals that primary magma for medium-K basalts was generated by partial melting of wedge mantle with normal potential temperature (T_P 1330–1340 °C) under wet conditions (H₂O 0.25–0.54 wt%). Its melting was probably induced by asthenospheric upwelling that is recognized by tomographic images. Mafic adakitic magma (Mg# ~ 63, Sr/Y ~ 64) was derived from the medium-K basaltic magma in fractional crystallization of a garnet-bearing assemblage at high pressure, and felsic adakitic rocks (Mg# ~ 45, Sr/Y ~ 50) were produced by assimilation-fractional crystallization processes at mid-crustal depths.

Adakites, High-Nb Basalts and Copper–Gold Deposits in Magmatic Arcs and Collisional Orogens: An Overview

Adakites are Y- and Yb-depleted, SiO2- and Sr-enriched rocks with elevated Sr/Y and La/Yb ratios originally thought to represent partial melts of subducted metabasalt, based on their association with the subduction of young (<25 Ma) and hot oceanic crust. Later, adakites were found in arc segments associated with oblique, slow and flat subduction, arc–transform intersections, collision zones and post-collisional extensional environments. New models of adakite petrogenesis include the melting of thickened and delaminated mafic lower crust, basalt underplating of the continental crust and high-pressure fractionation (amphibole ± garnet) of mantle-derived, hydrous mafic melts. In some cases, adakites are associated with Nb-enriched (10 ppm < Nb < 20 ppm) and high-Nb (Nb > 20 ppm) arc basalts in ancient and modern subduction zones (HNBs). Two types of HNBs are recognized on the basis of their geochemistry. Type I HNBs (Kamchatka, Honduras) share N-MORB-like isotopic and OIB-like trace element characteristics and most probably originate from adakite-contaminated mantle sources. Type II HNBs (Sulu arc, Jamaica) display high-field strength element enrichments in respect to island-arc basalts coupled with enriched, OIB-like isotopic signatures, suggesting derivation from asthenospheric mantle sources in arcs. Adakites and, to a lesser extent, HNBs are associated with Cu–Au porphyry and epithermal deposits in Cenozoic magmatic arcs (Kamchatka, Phlippines, Indonesia, Andean margin) and Paleozoic-Mesozoic (Central Asian and Tethyan) collisional orogens. This association is believed to be not just temporal and structural but also genetic due to the hydrous (common presence of amphibole and biotite), highly oxidized (>ΔFMQ > +2) and S-rich (anhydrite in modern Pinatubo and El Chichon adakite eruptions) nature of adakite magmas. Cretaceous adakites from the Stanovoy Suture Zone in Far East Russia contain Cu–Ag–Au and Cu–Zn–Mo–Ag alloys, native Au and Pt, cupriferous Ag in association witn barite and Ag-chloride. Stanovoy adakites also have systematically higher Au contents in comparison with volcanic arc magmas, suggesting that ore-forming hydrothermal fluids responsible for Cu–Au(Mo–Ag) porphyry and epithermal mineralization in upper crustal environments could have been exsolved from metal-saturated, H2O–S–Cl-rich adakite magmas. The interaction between depleted mantle peridotites and metal-rich adakites appears to be capable of producing (under a certain set of conditions) fertile sources for HNB melts connected with some epithermal Au (Porgera) and porphyry Cu–Au–Mo (Tibet, Iran) mineralized systems in modern and ancient subduction zones.

Iron-Copper-Zinc Isotopic Compositions of Andesites from the Kueishantao Hydrothermal Field off Northeastern Taiwan

The studies of iron (Fe), copper (Cu), and zinc (Zn) isotopic compositions in seafloor andesites are helpful in understanding the metal stable isotope fractionation during magma evolution. Here, the Fe, Cu, and Zn isotopic compositions of andesites from the Kueishantao hydrothermal field (KHF) off northeastern Taiwan, west Pacific, have been studied. The majority of δ56Fe values (+0.02‰ to +0.11‰) in the KHF andesites are consistent with those of MORBs (mid-ocean ridge basalts). This suggests that the Fe in the KHF andesites is mainly from a MORB-type mantle. The Fe-Cu-Zn isotopic compositions (δ56Fe +0.22‰, δ65Cu +0.16‰ to +0.64‰, and δ66Zn +0.29‰ to +0.71‰) of the KHF andesites, which are significantly different from those of the MORBs and the continental crust (CC), have a relatively wide range of Cu and Zn isotopic compositions. This is most likely to be a result of the entrainment of the sedimentary carbonate-derived components into an andesitic magma. The recycled altered rocks (higher δ56Fe, lower δ66Zn) could preferentially incorporate isotopically light Fe and heavy Zn into the magma, resulting in relative enrichment of the lighter Fe and heavier Zn isotopes in the andesites. The majority of the δ56Fe values in the KHF andesites are higher than those of the sediments and the local CC and lower than those of the subducted altered rocks, while the reverse is true for δ66Zn, suggesting that the subseafloor sediments and CC materials (lower δ56Fe, higher δ66Zn) contaminating the rising andesitic magma could preferentially incorporate isotopically heavy Fe and light Zn into the magma, resulting in relative enrichment of the heavier Fe and lighter Zn isotopes in the andesites. Thus, the characteristics of the Fe and Zn isotopes in back-arc and island-arc volcanic rocks may also be influenced by the CC and plate subduction components.

Multistage Magmatism in Ophiolites and Associated Metavolcanites of the Ulan-Sar’dag Mélange (East Sayan, Russia)

We present new whole-rock major and trace element, mineral chemistry, and U-Pb isotope data for the Ulan-Sar’dag mélange, including different lithostratigraphic units: Ophiolitic, mafic rocks and metavolcanites. The Ulan-Sar’dag mélange comprises of a seafloor and island-arc system of remnants of the Paleo-Asian Ocean. Detailed studies on the magmatic rocks led to the discovery of a rock association that possesses differing geochemical signatures within the studied area. The Ulan-Sar’dag mélange includes blocks of mantle peridotite, podiform chromitite, cumulate rocks, deep-water siliceous chert, and metavolcanic rocks of the Ilchir suite. The ophiolitic unit shows overturned pseudostratigraphy. The nappe of mantle tectonites is thrusted over the volcanic-sedimentary sequence of the Ilchir suite. The metavolcanic series consist of basic, intermediate, and alkaline rocks. The mantle peridotite and cumulate rocks formed in a supra-subduction zone environment. The mafic and metavolcanic rocks belong to the following geochemical types: (1) Ensimatic island-arc boninites; (2) island-arc calc-alkaline andesitic basalts, andesites, and dacites; (3) tholeiitic basalts of mid-ocean ridges; and (4) oceanic island basalts. U–Pb dating of zircons from the trachyandesite, belonging to the second geochemical type, yielded a date of 833 ± 4 Ma which is interpreted as the crystallization age during mature island-arc and intra-arc rifting stages. The possible influence of later plume magmatic-hydrothermal activities led to the appearance of moderately alkaline igneous rocks (monzogabbro, trachybasalt, trachyandesite, subalkaline gabbro, and metasomatized peridotites) with a significant subduction geochemical fingerprint.

A multi-elemental approach to assess potential contamination in tourist beaches: The case of Loreto Bay (Marine Protected Area), NW Mexico

Metal/metalloid concentrations in water sediment and commercial fishes of Loreto Maritime National Park (MNP), Baja California Sur, Mexico were determined for a comprehensive geochemical study. In-situ physical characteristics (pH, conductivity, redox potential, dissolved oxygen, turbidity) of water clearly indicated the unique oceanographic properties of the Gulf of California. Likewise, the distribution pattern of metals/metalloid in water, sediments and fishes denoted the influences of local geology, longshore currents, upwelling process, natural hydrothermal vents and the 100-year old mining activities of Santa Rosalia region, situated to the north of Loreto. Calculated carcinogenic indices in commercial fish species showed safe human consumption. Thus, the present research validates a comprehensive geochemical study of protected areas upholding the need for continuous monitoring for a better conservation of coastal ecosystems.

Growth of, and diffusion in, olivine in ultra-fast ascending basalt magmas from Shiveluch volcano

Complex core-rim zoning of Mg-Fe-Ni-Ca-Cr-Al-P in high-Mg olivine crystals from a tuff ring of Shiveluch volcano, Kamchatka, enables reconstruction of the entire olivine crystallization history from mantle conditions to eruption. Bell-shaped Fo_86–92 and Ni profiles in crystal cores were formed by diffusion after mixing with evolved magma. Diffusion proceeded to the centres of crystals and completely equilibrated Fo and Ni in some crystals. Diffusion times extracted from Fo and Ni core profiles range from 100 to 2000 days. During subsequent mixing with mafic mantle-equilibrated melt, the cores were partially dissolved and overgrown by Fo₉₀ olivine. Times extracted from Fo and Ni diffusion profiles across the resorption interface between the core and its overgrowth range within 1–10 days, which corresponds to the time of magma ascent to the surface. The overgrowth shows identical smooth Fo-Ni decreasing zoning patterns for all crystals towards the margin, indicating that all crystals shared the same growth history after last mixing event prior to eruption. At the same time, Ca, and to an even greater extent Cr, Al, and P have oscillatory growth patterns in the crystals overgrowth. Our data show that magma ascent can be extremely short during maar/tuff ring eruption.

Hepatocellular Carcinoma in the Cirrhotic Liver: Evaluation Using Computed Tomography and Magnetic Resonance Imaging

Hepatocellular carcinoma is the fifth most common tumor in patients worldwide and the third most common cause of cancer-related death, after lung and stomach cancer. Cirrhosis of the liver is the strongest predisposing factor for hepatocellular carcinoma, with approximately 80% of cases of hepatocellular carcinoma developing in a cirrhotic liver. The annual incidence of hepatocellular carcinoma is 2.0% to 6.6% in patients with cirrhosis compared with 0.4% in patients without cirrhosis. The 5-year survival rates of patients undergoing curative therapies for hepatocellular carcinoma, including liver transplant, hepatic resection, and percutaneous ablative techniques, range between 40% and 75%. Orthotropic liver transplant offers the prima facie cure for both hepatocellular carcinoma and liver cirrhosis. In hepatocellular carcinoma confined to the liver without macrovascular invasion, patients with a single tumor ≤ 5 cm or up to 3 tumors ≤ 3 cm each had a 5-year survival rate of 75% and a disease-free survival rate of 83%. In the adult population, liver transplant for hepatocellular carcinoma yields good results for patients whose tumor masses do not exceed the Milan criteria. The diagnosis of hepatocellular carcinoma using imaging tests has had a substantial impact on transplant decisions. Radiologists should be aware of this responsibility and exercise the utmost scrutiny before making a diagnosis of hepatocellular carcinoma. Erroneous diagnosis of hepatocellular carcinoma based on imaging tests could deny deserving patients the opportunity of a life-saving liver transplant and result in unnecessary liver transplants for others. Contrast-enhanced magnetic resonance imaging and helical computed tomography are the best imaging techniques currently available for the noninvasive diagnosis of hepatocellular carcinoma. With technological advances in hardware and software, diffusion-weighted imaging can be readily applied to the liver with resulting improved image quality.

Genesis of ultra-high-Ni olivine in high-Mg andesite lava triggered by seamount subduction

The Kamchatka Peninsula is a prominent and wide volcanic arc located near the northern edge of the Pacific Plate. It has highly active volcanic chains and groups, and characteristic lavas that include adakitic rocks. In the north of the peninsula adjacent to the triple junction, some additional processes such as hot asthenospheric injection around the slab edge and seamount subduction operate, which might enhance local magmatism. In the forearc area of the northeastern part of the peninsula, monogenetic volcanic cones dated at <1 Ma were found. Despite their limited spatiotemporal occurrence, remarkable variations were observed, including primitive basalt and high-Mg andesite containing high-Ni (up to 6300 ppm) olivine. The melting and crystallization conditions of these lavas indicate a locally warm slab, facilitating dehydration beneath the forearc region, and a relatively cold overlying mantle wedge fluxed heterogeneously by slab-derived fluids. It is suggested that the collapse of a subducted seamount triggered the ascent of Si-rich fluids to vein the wedge peridotite and formed a peridotite–pyroxenite source, causing the temporal evolution of local magmatism with wide compositional range.

Geochemical evidence for mélange melting in global arcs

In subduction zones, sediments and hydrothermally altered oceanic crust, which together form part of the subducting slab, contribute to the chemical composition of lavas erupted at the surface to form volcanic arcs. Transport of this material from the slab to the overlying mantle wedge is thought to involve discreet melts and fluids that are released from various portions of the slab. We use a meta-analysis of geochemical data from eight globally representative arcs to show that melts and fluids from individual slab components cannot be responsible for the formation of arc lavas. Instead, the data are compatible with models that first invoke physical mixing of slab components and the mantle wedge, widely referred to as high-pressure mélange, before arc magmas are generated.

Partial melting of deeply subducted eclogite from the Sulu orogen in China

We report partial melting of an ultrahigh pressure eclogite in the Mesozoic Sulu orogen, China. Eclogitic migmatite shows successive stages of initial intragranular and grain boundary melt droplets, which grow into a three-dimensional interconnected intergranular network, then segregate and accumulate in pressure shadow areas and then merge to form melt channels and dikes that transport magma to higher in the lithosphere. Here we show, using zircon U–Pb dating and petrological analyses, that partial melting occurred at 228–219 Myr ago, shortly after peak metamorphism at 230 Myr ago. The melts and residues are complimentarily enriched and depleted in light rare earth element (LREE) compared with the original rock. Partial melting of deeply subducted eclogite is an important process in determining the rheological structure and mechanical behaviour of subducted lithosphere and its rapid exhumation, controlling the flow of deep lithospheric material, and for generation of melts from the upper mantle, potentially contributing to arc magmatism and growth of continental crust.

Natural examples of eclogite-derived migmatite were previously unknown. Here, the authors show evidence of high-pressure melting of eclogite in the Sulu orogen of China, and suggest the melts represent an exhumed example of flow channels in the lower crust, possibly a source of syn-exhumation magmas.

Laboratory models of the thermal evolution of the mantle during rollback subduction

The subduction of oceanic lithosphere plays a key role in plate tectonics, the thermal evolution of the mantle and recycling processes between Earth's interior and surface. Information on mantle flow, thermal conditions and chemical transport in subduction zones come from the geochemistry of arc volcanoes, seismic images and geodynamic models. The majority of this work considers subduction as a two-dimensional process, assuming limited variability in the direction parallel to the trench. In contrast, observationally based models increasingly appeal to three-dimensional flow associated with trench migration and the sinking of oceanic plates with a translational component of motion (rollback). Here we report results from laboratory experiments that reveal fundamental differences in three-dimensional mantle circulation and temperature structure in response to subduction with and without a rollback component. Without rollback motion, flow in the mantle wedge is sluggish, there is no mass flux around the plate and plate edges heat up faster than plate centres. In contrast, during rollback subduction flow is driven around and beneath the sinking plate, velocities increase within the mantle wedge and are focused towards the centre of the plate, and the surface of the plate heats more along the centreline.

Back-arc spreading and mantle flow in the East Scotia Sea

The East Scotia Ridge exhibits systematic variations in axial morphology and basalt geochemistry. Central segments have morphology typical of intermediate-rate spreading centres and erupt mainly normal mid-ocean ridge basalt (N-MORB). Segments near the ridge ends exhibit anomalous, inflated, axial morphology and erupt more evolved basalts, influenced by the Bouvet plume in the north. As the end segments lie closer to the volcanic arc, these variations could be caused by coupled flow within the mantle wedge, as inferred from similar studies in the Lau Basin. Three of the four zones of crustal accretion defined from the Lau Basin may be identified in the East Scotia Sea, although there is no counterpart to a zone of diminished magma supply observed at the East Lau Spreading Centre. Superimposed on the pattern of plate-driven flow is a ridge-parallel flow related to inflow of Atlantic mantle into the East Scotia Sea back-arc region at both ends of the South Sandwich slab. The inflow causes enhanced magmatism and propagation of the end segments towards the middle of the back-arc region, and may be related to trench-parallel flow beneath the rapidly retreating slab. Alternatively, it may be driven by buoyancy flux from Atlantic hot spots. There is no evidence that retreat was ever driven by escape flow of Pacific mantle.

Seismic evidence for catastrophic slab loss beneath Kamchatka

In the northwest Pacific Ocean, a sharp corner in the boundary between the Pacific plate and the North American plate joins a subduction zone running along the southern half of the Kamchatka peninsula with a region of transcurrent motion along the western Aleutian arc. Here we present images of the seismic structure beneath the Aleutian-Kamchatka junction and the surrounding region, indicating that: the subducting Pacific lithosphere terminates at the Aleutian-Kamchatka junction; no relict slab underlies the extinct northern Kamchatka volcanic arc; and the upper mantle beneath northern Kamchatka has unusually slow shear wavespeeds. From the tectonic and volcanic evolution of Kamchatka over the past 10 Myr (refs 3-5) we infer that at least two episodes of catastrophic slab loss have occurred. About 5 to 10 Myr ago, catastrophic slab loss shut down island-arc volcanic activity north of the Aleutian-Kamchatka junction. A later episode of slab loss, since about 2 Myr ago, seems to be related to the activity of the world's most productive island-arc volcano, Klyuchevskoy. Removal of lithospheric mantle is commonly discussed in the context of a continental collision, but our findings imply that episodes of slab detachment and loss are also important agents in the evolution of oceanic convergent margins.