Sub-arc mantle enrichment in the Sunda rear-arc inferred from HFSE systematics in high-K lavas from Java

Many terrestrial silicate reservoirs display a characteristic depletion in Nb, which has been explained in some studies by the presence of reservoirs on Earth with superchondritic Nb/Ta. As one classical example, K-rich lavas from the Sunda rear-arc, Indonesia, have been invoked to tap such a high-Nb/Ta reservoir. To elucidate the petrogenetic processes active beneath the Java rear-arc and the causes for the superchondritic Nb/Ta in some of these lavas, we studied samples from the somewhat enigmatic Javanese rear-arc volcano Muria, which allow conclusions regarding the across-arc variations in volcanic output, source mineralogy and subduction components. We additionally report some data for an along-arc sequence of lavas from the Indonesian part of the Sunda arc, extending from Krakatoa in the west to the islands of Bali and Lombok in the east. We present major and trace element concentrations, Sr–Nd–Hf–Pb isotope compositions, and high-field-strength element (HFSE: Nb, Ta, Zr, Hf, W) concentrations obtained via isotope dilution and MC-ICP-MS analyses. The geochemical data are complemented by melting models covering different source compositions with slab melts formed at variable P–T conditions. The radiogenic isotope compositions of the frontal arc lavas in combination with their trace element systematics confirm previously established regional variations of subduction components along the arc. Melting models show a clear contribution of a sediment-derived component to the HFSE budget of the frontal arc lavas, particularly affecting Zr–Hf and W. In contrast, the K-rich rear-arc lavas tap more hybrid and enriched mantle sources. The HFSE budget of the rear-arc lavas is in particular characterized by superchondritic Nb/Ta (up to 25) that are attributed to deep melting involving overprint by slab melts formed from an enriched garnet–rutile-bearing eclogitic residue. Sub-arc slab melting was potentially triggered along a slab tear beneath the Sunda arc, which is the result of the forced subduction of an oceanic basement relief ~ 8 Myr ago as confirmed by geophysical studies. The purported age of the slab tear coincides with a paucity in arc volcanism, widespread thrusting of the Javanese basement crust as well as the short-lived nature of the K-rich rear-arc volcanism at that time.

Early Results of Eastern Indonesia P-wave Tomography Study Using Regional Events

The tectonic system of Eastern Indonesia is controlled by several major and minor plates, such as Indo-Australian, Australian plate, and Pacific plates. This area is known for its complexity, and high seismic activity. This study tries to image the complex structures beneath this region by employing regional events data and seismic tomography methods. We used five years of regional events catalog provided by the Indonesian Agency of Meteorology, Climatology, and Geophysics. We have sorted 7336 events recorded between 120° – 136° longitude and 0° – 13°(-) latitude consisting of 46446 P and 15467 S wave arrival data. Relocated hypocenter map shows a better constrain location on seismicity along outer Bandar Arc. A dipping pattern of seismicity is seen that is going deeper to the Banda Sea. The seismicity map also images a steep angle pattern of seismicity that could be related to the subduction slab roll-back model at North of Wetar island. Interestingly, we spotted a seismicity gap in West Seram that could be linked with slab tear zone. The checker-board test suggests a proper resolution is still reliable to a depth of 200 km with a less interpretable model at a depth of 300 km. P-wave tomographic models image the high velocity dipping down going slab. The Banda slab is seen to subduct from south Timor Island to the north, from east Tanimbar and Aru Island to west part, and from north Seram Island to south. We observed the down-going slab meet from all directions at about 300 km beneath the Banda sea. P wave tomogram also shows the Timor Island slab has a steeper dip that agrees with the seismicity pattern. Near the Seram island, we identify a low-velocity anomaly zone infiltrate the Banda slab beneath the shallow part of West Seram, which was previously interpreted as slab tear zone. This study also noticed a higher velocity tomogram model at North of Wetar island that might indicate a back-arc thrust. Lastly, a low-velocity band is also exposed at a shallow depth close to the volcano chain along that Banda volcanic arc.

Timescales of continental subduction: Constraints from ultrahigh-pressure metapelites in the Western Gneiss Region, Norway

<p>The Western Gneiss Region (WGR), Norway is an archetypal continental ultrahigh-pressure (U)HP terrane with an extensive metamorphic history, recording the subduction and subsequent exhumation of continental crust to depths exceeding 120 km. The vast bulk of past work within the WGR has focused on mafic eclogites. In this study, data from rare garnet-kyanite metapelites in (UHP) domains of the WGR is presented. U–Pb geochronology and trace element compositions in zircon, monazite, apatite, rutile and garnet were acquired, and P–T conditions were calculated by mineral equilibria forward modelling and Zr-in-rutile thermometry. The Ulsteinvik metapelite defines a prograde path that traverses through ~600–710 °C and ~11–14 kbar. Minimum peak conditions are ~750 °C and ~2.9 GPa in an inferred garnet-kyanite-coesite-omphacite-muscovite-rutile-quartz-H<sub>2</sub>O assemblage. Plagioclase-biotite-quartz intergrowths developed after omphacite-phengite-rutile breakdown on the early retrograde path, followed by cordierite-spinel-plagioclase symplectites after garnet-kyanite-biotite, defining a retrograde P–T point at ~740 °C and ~7 kbar. Late Ordovician-Early Silurian (~470–440 Ma) zircon and rutile age data in Ulsteinvik pre-dates the major Scandian UHP subduction episode in the WGR, interpreted as recording early Caledonian subduction within the Blåhø nappe. Monazite and apatite U-Pb geochronology and trace element data suggest exhumation occurred at ~400 Ma. The Fjørtoft metapelite is a constituent of the Blåhø nappe. Minimum peak P–T conditions are ~1.8 GPa and ~750 °C, with poor peak mineral fidelity attributed to extensive retrograde deformation. Negative Eu anomalies in ~423 Ma monazite suggest retrograde conditions were reached [RJT1] by ~423 Ma. Ulsteinvik and Fjørtoft may have experienced pre-Scandian subduction together within the Blåhø nappe, but record dissimilar histories after this. Two potential scenarios are presented: (1) Ulsteinvik resided within the mantle for 20 million-years longer than Fjørtoft during Scandian subduction, or (2), the samples were exhumed at different times during pre-Scandian subduction of the Blåhø nappe. The preservation of prograde zoning within Ulsteinvik garnets precludes a long-term residence within the mantle and suggests the latter option. In this scenario, the subducting Blåhø nappe experienced a degree of slab tear and partial underplating of the upper plate during the early stages of continental underthrusting. Discrete pieces may have later reattached to the lower plate at different times, partially exhumed, and then subducted to mantle-depths during the Scandian.</p>

A Laurentian margin subduction perspective: Geodynamic constraints from phase equilibria modeling of barroisite greenstones, northern USA Appalachians

Phase equilibria modeling of sodic-calcic amphibole-epidote assemblages in greenstones in the northern Appalachians, USA, is compatible with relatively shallow subduction of the early Paleozoic Laurentian margin along the Laurentia-Gondwana suture zone during closure of a portion of the Iapetus Ocean basin. Pseudosection and isopleth calculations demonstrate that peak metamorphic conditions ranged between 0.65 GPa, 480 °C and 0.85 GPa, 495 °C down-dip along the subducted Laurentian continental margin between ∼20 km and ∼30 km depth. Quantitative petrological data are explained in the context of an Early Ordovician geodynamic model involving shallow subduction of relatively young, warm, and buoyant Laurentian margin continental-oceanic lithosphere and Iapetus Ocean crust beneath a relatively warm and wet peri-Gondwanan continental arc. A relatively warm subduction zone setting may have contributed to the formation of a thin, ductile metasedimentary rock-rich channel between the down-going Laurentian slab and the overriding continental arc. This accretionary channel accommodated metamorphism and tectonization of continental margin sediments and mafic volcanic rocks (greenstones) of the Laurentian margin and provided a pathway for exhumation of serpentinite slivers and rare eclogite blocks. Restricted asthenospheric flow in the forearc mantle wedge provides one explanation for the lack of ophiolites and absence of a well-preserved ultra-high-pressure terrane in central and northern Vermont. Exhumation of the subducted portion of the Laurentian margin may have been temperature triggered due to increased asthenospheric flow following a slab tear at relatively shallow depths.

Image of the Iberian Tethys paleomargin beneath the eastern Betic mountain range

<p>We obtain P-wave receiver functions from recordings at a dense seismic broadband transect, deployed along 170 km across the eastern Betic orogen in south Spain. Migrated images show the crustal structure of the orogen in detail. In particular, they reveal the situation of the subducted Iberian paleomargin, with full preservation of the proximal domain and the ~50 km wide necking domain. Crustal thinning affects the lower continental crust. The Variscan crust of the Tethys margin is bending downward beneath the Betics, reaching ~45 km depth, and terminates abruptly at a major slab tear fault. The distal domain of the paleomargin cannot be reconstructed, but the migrated section suggests that material has been exhumed through the subduction channel and integrated into the Betic Orogene. This supports an origin of the HP-LT Nevado-Filabride units from subducted, hyperextended Variscan crust.</p>

Chapter 23: Alteration, Mineralization, and Age Relationships at the Kışladağ Porphyry Gold Deposit, Turkey

The Miocene Kışladağ deposit (~17 Moz), located in western Anatolia, Turkey, is one of the few global examples of Au-only porphyry deposits. It occurs within the West Tethyan magmatic belt that can be divided into Cretaceous, Cu-dominant, subduction-related magmatic arc systems and the more widespread Au-rich Cenozoic magmatic belts. In western Anatolia, Miocene magmatism was postcollisional and was focused in extension-related volcanosedimentary basins that formed in response to slab roll back and a major north-south slab tear. Kışladağ formed within multiple monzonite porphyry stocks and dikes at the contact between Menderes massif metamorphic basement and volcanic rocks of the Beydağı stratovolcano in the Uşak-Güre basin. The mineralized magmatic-hydrothermal system formed rapidly (<400 kyr) between ~14.75 and 14.36 Ma in a shallow (<1 km) volcanic environment. Volcanism continued to at least 14.26 ± 0.09 Ma based on new age data from a latite lava flow at nearby Emiril Tepe. Intrusions 1 and 2 were the earliest (14.73 ± 0.05 and 14.76 ± 0.01 Ma, respectively) and best mineralized phases (average median grades of 0.64 and 0.51 g/t Au, respectively), whereas younger intrusions host progressively less Au (Intrusion 2A: 14.60 ± 0.06 Ma and 0.41 g/t Au; Intrusion 2 NW: 14.45 ± 0.08 Ma and 0.41 g/t Au; Intrusion 3: 14.39 ± 0.06 and 14.36 ± 0.13 Ma and 0.19 g/t Au). A new molybdenite age of 14.60 ± 0.07 Ma is within uncertainty of the previously published molybdenite age (14.49 ± 0.06 Ma), and supports field observations that the bulk of the mineralization formed prior to the emplacement of Intrusion 3. Intrusions 1 and 2 are altered to potassic (biotite-K-feldspar-quartz ± magnetite) and younger but deeper sodic-calcic (feldspar-amphibole-magnetite ± quartz ± carbonate) assemblages, both typically pervasive with disseminated to veinlet-hosted pyrite ± chalcopyrite ± molybdenite and localized quartz-feldspar stockwork veinlets and sodic-calcic breccias. Tourmaline-white mica-quartz-pyrite alteration surrounds the potassic core both within the intrusions and outboard in the volcanic rocks. Tourmaline was most strongly developed on the inner margins of the tourmaline-white mica zone, particularly along the Intrusion 1 volcanic contact where it formed breccias and veins, including Maricunga-style veinlets. Field relationships show that the early magmatic-hydrothermal events were cut by Intrusion 2A, which was then overprinted by Au-bearing argillic (kaolinite-pyrite ± quartz) alteration, followed by Intrusion 3 and late-stage, low-grade to barren argillic and advanced argillic alteration (quartz-pyrite ± alunite ± dickite ± pyrophyllite). Gold deportment changes with each successive hydrothermal event. The early potassic and sodic-calcic alteration controls much of the original Au distribution, with the Au dominantly deposited with feldspar and lesser quartz and pyrite. Tourmaline-white mica and argillic alteration events overprinted and altered the early Au-bearing feldspathic alteration and introduced additional Au that was dominantly associated with pyrite. Analogous Au-only deposits such as Maricunga, Chile, La Colosa, Colombia, and Biely Vrch, Slovakia, are characterized by similar alteration styles and Au deportment. The deportment of Au in these Au-only porphyry deposits differs markedly from that in Au-rich porphyry Cu deposits where Au is typically associated with Cu sulfides.

Surface and mantle records reveal an ancient slab tear beneath Gondwana

Vertical slab-tearing has been widely reported in modern convergent settings profoundly influencing subduction and mantle dynamics. However, evaluating a similar impact in ancient convergent settings, where oceanic plates have been subducted and the geological record is limited, remains challenging. In this study, we correlate the lower mantle structure, which retained the past subduction configuration, with the upper-plate geological record to show a deep slab rupture interpreted as a large-scale tearing event in the early Mesozoic beneath southwestern Gondwana. For this purpose, we integrated geochronological and geological datasets with P-wave global seismic tomography and plate-kinematic reconstructions. The development of a Late Triassic-Early Jurassic slab-tearing episode supports (i) a slab gap at lower mantle depths, (ii) a contrasting spatiotemporal magmatic evolution, (iii) a lull in arc activity, and (iv) intraplate extension and magmatism in the Neuquén and Colorado basins. This finding not only has implications for identifying past examples of a fundamental process that shapes subduction zones, but also illustrates an additional mechanism to trigger slab-tearing in which plate rupture is caused by opposite rotation of slab segments.

Comparison of 3-D Raytracing and Finite Frequency Tomography

We performed 3-D seismic tomography using teleseismic arrival time at Southwest Mexico. The Mexican subduction zone results from successive fragmentation events that affected the ancient Farallon plate as various segments of the East Pacific rise approached the paleo-trench off western North America. The complexity in this region is related to two subducting oceanic plates, the Rivera and Cocos plates, that have different ages, compositions, convergence velocities and subduction dip angles. In this study, we compared the 3-D raytracing tomography model with finite frequency tomography model. Final models show the differences in amplitude and pattern between the raytracing and finite frequency. 3D raytracing models produced sharper images of fast velocity structures in the mantle. The deeper slabs are more coherent and show less broadening with depth than using 1D finite frequency kernels. However, although the finite frequency and 3-D ray tracing models show some differences in amplitude and pattern, the overall agreement of the models supports the interpretation of Yang et al. (2009) that slab rollback is occurring in South Western Mexico. One possible different interpretation between the raytracing and finite frequency theory results concerns the deep structure of the Rivera slab. The finite frequency models show that the Rivera slab is clearly observable at a depth of about 300km but fades away at greater depths. However, the 3-D ray tracing model shows a clear fast velocity band down to a depth of 400 km and thus our model does not support a slab tear of the Rivera plate above 400 km depth