Volatile (H<Subscript>2</Subscript>O,CO<Subscript>2</Subscript>, Cl,S) budget of the Central American subduction zone

After more than a decade of multidisciplinary studies of the Central American subduction zone mainly in the framework of two large research programmes, the US MARGINS program and the German Collaborative Research Center SFB 574, we here review and interpret the data pertinent to quantify the cycling of mineral-bound volatiles (H₂O, CO₂, Cl, S) through this subduction system. For input-flux calculations, we divide the Middle America Trench into four segments differing in convergence rate and slab lithological profiles, use the latest evidence for mantle serpentinization of the Cocos slab approaching the trench, and for the first time explicitly include subduction erosion of forearc basement. Resulting input fluxes are 40–62 (53) Tg/Ma/m H₂O, 7.8–11.4 (9.3) Tg/Ma/m CO₂, 1.3–1.9 (1.6) Tg/Ma/m Cl, and 1.3–2.1 (1.6) Tg/Ma/m S (bracketed are mean values for entire trench length). Output by cold seeps on the forearc amounts to 0.625–1.25 Tg/Ma/m H₂O partly derived from the slab sediments as determined by geochemical analyses of fluids and carbonates. The major volatile output occurs at the Central American volcanic arc that is divided into ten arc segments by dextral strike-slip tectonics. Based on volcanic edifice and widespread tephra volumes as well as calculated parental magma masses needed to form observed evolved compositions, we determine long-term (10⁵ years) average magma and K₂O fluxes for each of the ten segments as 32–242 (106) Tg/Ma/m magma and 0.28–2.91 (1.38) Tg/Ma/m K₂O (bracketed are mean values for entire Central American volcanic arc length). Volatile/K₂O concentration ratios derived from melt inclusion analyses and petrologic modelling then allow to calculate volatile fluxes as 1.02–14.3 (6.2) Tg/Ma/m H₂O, 0.02–0.45 (0.17) Tg/Ma/m CO₂, and 0.07–0.34 (0.22) Tg/Ma/m Cl. The same approach yields long-term sulfur fluxes of 0.12–1.08 (0.54) Tg/Ma/m while present-day open-vent SO₂-flux monitoring yields 0.06–2.37 (0.83) Tg/Ma/m S. Input–output comparisons show that the arc water fluxes only account for up to 40 % of the input even if we include an “invisible” plutonic component constrained by crustal growth. With 20–30 % of the H₂O input transferred into the deeper mantle as suggested by petrologic modeling, there remains a deficiency of, say, 30–40 % in the water budget. At least some of this water is transferred into two upper-plate regions of low seismic velocity and electrical resistivity whose sizes vary along arc: one region widely envelopes the melt ascent paths from slab top to arc and the other extends obliquely from the slab below the forearc to below the arc. Whether these reservoirs are transient or steady remains unknown.     

Magma genesis beneath Northeast Japan arc: A new perspective on subduction zone magmatism

It is being accepted that earthquakes in subducting slab are caused by dehydration reactions of hydrous minerals. In the context of this “dehydration embrittlement” hypothesis, we propose a new model to explain key features of subduction zone magmatism on the basis of hydrous phase relations in peridotite and basaltic systems determined by thermodynamic calculations and seismic structures of Northeast Japan arc revealed by latest seismic studies. The model predicts that partial melting of basaltic crust in the subducting slab is an inevitable consequence of subduction of hydrated oceanic lithosphere. Aqueous fluids released from the subducting slab also cause partial melting widely in mantle wedge from just above the subducting slab to just below overlying crust at volcanic front. Hydrous minerals in the mantle wedge are stable only in shallow (< 120 km) areas, and are absent in the layer that is dragged into deep mantle by the subducting slab. The position of volcanic front is not restricted by dehydration reactions in the subducting slab but is controlled by dynamics of mantle wedge flow, which governs the thermal structure and partial melting regime in the mantle wedge.     

Nitrogen recycling in subducted mantle rocks and implications for the global nitrogen cycle

The nitrogen concentrations [N] and isotopic compositions of ultramafic mantle rocks that represent various dehydration stages and metamorphic conditions during the subduction cycle were investigated to assess the role of such rocks in deep-Earth N cycling. The samples analyzed record low-grade serpentinization on the seafloor and/or in the forearc wedge (low-grade serpentinites from Monte Nero/Italy and Erro Tobbio/Italy) and two successive stages of metamorphic dehydration at increasing pressures and temperatures (high-pressure (HP) serpentinites from Erro Tobbio/Italy and chlorite harzburgites from Cerro del Almirez/Spain) to allow for the determination of dehydration effects in ultramafic rocks on the N budget. In low-grade serpentinites, δ¹⁵N_air values (?3.8 to +3.5 ‰) and [N] (1.3–4.5 μg/g) are elevated compared to the pristine depleted MORB mantle (δ¹⁵N_air ~ ?5 ‰, [N] = 0.27 ± 0.16 μg/g), indicating input from sedimentary organic sources, at the outer rise during slab bending and/or in the forearc mantle wedge during hydration by slab-derived fluids. Both HP serpentinites and chlorite harzburgites have δ¹⁵N_air values and [N] overlapping with low-grade serpentinites, indicating no significant loss of N during metamorphic dehydration and retention of N to depths of 60–70 km. The best estimate for the δ¹⁵N_air of ultramafic rocks recycled into the mantle is +3 ± 2 ‰. The global N subduction input flux in serpentinized oceanic mantle rocks was calculated as 2.3 × 10⁸ mol N₂/year, assuming a thickness of serpentinized slab mantle of 500 m. This is at least one order of magnitude smaller than the N fluxes calculated for sediments and altered oceanic crust. Calculated global input fluxes for a range of representative subducting sections of unmetamorphosed and HP-metamorphosed slabs, all incorporating serpentinized slab mantle, range from 1.1 × 10¹⁰ to 3.9 × 10¹⁰ mol N₂/year. The best estimate for the δ¹⁵N_air of the subducting slab is +4 ± 1 ‰, supporting models that invoke recycling of subducted N in mantle plumes and consistent with general models for the volatile evolution on Earth. Estimates of the efficiency of arc return of subducted N are complicated further by the possibility that mantle wedge hydrated in forearcs, then dragged to beneath volcanic fronts, is capable of conveying significant amounts of N to subarc depths.     

The Pushtashan juvenile suprasubduction zone assemblage of Kurdistan (northeastern Iraq): A Cretaceous (Cenomanian) Neo-Tethys missing link

The Pushtashan suprasubduction zone assemblage of volcanic rocks, gabbros, norites and peridotites occurs in the Zagros suture zone, Kurdistan region, northeastern Iraq. Volcanic rocks are dominant in the assemblage and consist mainly of basalt and basaltic andesite flows with interlayered red shale and limestone horizons. Earlier lavas tend to be MORB-like, whereas later lavas display island arc tholeiite to boninitic geochemical characteristics. Tholeiitic gabbros intrude the norites and display fractionation trends typical of crystallisation under low-pressure conditions, whereas the norites display calc-alkaline traits, suggesting their source included mantle metasomatised by fluids released from subducted oceanic crust. Enrichment of Rb, Ba, Sr, Th and the presence of negative Nb anomalies indicate generation in a suprasubduction zone setting. Trondhjemite and granodiorite intrusions are present in the volcanic rocks, gabbros and norites. SHRIMP U-Pb dating of magmatic zircons from a granodiorite yields a mean~(206)Pb/~(238)U age of 96.0 ±2.0 Ma(Cenomanian). The initial ε_(Hf) value for the zircons show a narrow range from +12.8 to+15.6, with a weighted mean of + 13.90±0.96. This initial value is within error of model depleted mantle at 96 Ma or slightly below that, in the field of arc rocks with minimal contamination by older continental crust. The compositional bimodality of the Pushtashan suprasubduction sequence suggests seafloor spreading during the initiation of subduction, with a lava stratigraphy from earlyerupted MORB transitioning into calc-alkaline lavas and finally by 96 Ma intrusion of granodioritic and trondhjemitic bodies with juvenile crustal isotopic signatures. The results confirm another Cretaceous arc remnant preserved as an allochthon within the Iraqi segment of the Cenozoic Zagros suture zone. Implications for the closure of Neo-Tethys are discussed.     

Tomographic evidence for the subducting oceanic crust and forearc mantle serpentinization under Kyushu, Japan

We determined high-resolution three-dimensional P- and S-wave velocity (Vp, Vs) structures beneath Kyushu in Southwest Japan using 177,500 P and 174,025 S wave arrival times from 8515 local earthquakes. A Poisson's ratio structure was derived from the obtained Vp and Vs values. Our results show that significant low-Vp, low-Vs and high Poisson's ratio zones are extensively distributed along the volcanic front in the uppermost mantle, which extend and dip toward the back-arc side in the mantle wedge. In the crust, low-Vp, low-Vs and high Poisson's ratio anomalies exist beneath the active volcanoes. The subducting Philippine Sea slab is clearly imaged as a high-Vp, high-Vs and low Poisson's ratio zone from the Nankai Trough to the back-arc. A thin low-velocity zone is detected above the subducting Philippine Sea slab in the mantle wedge, and earthquakes in the upper mantle are distributed along the transition zone between this thin low-velocity zone and the high-velocity Philippine Sea slab, which may imply that oceanic crust exists on the top of the slab and the forearc mantle wedge is serpentinized due to the slab dehydration. The seismic velocity of the subducting oceanic crust with basaltic or gabbroic composition is lower than that of the mantle according to the previous studies. The serpentinization process could also dramatically reduce the seismic velocity in the forearc mantle wedge.     

Influence of the precollisional stage on subduction dynamics and the buried crust thermal state: Insights from numerical simulations

At continental subduction initiation, the continental crust buoyancy may induce, first, a convergence slowdown, and second, a compressive stress increase that could lead to the forearc lithosphere rupture. Both processes could influence the slab surface P–T conditions, favoring on one side crust partial melting or on the opposite the formation of ultra-high pressure/low temperature (UHP-LT) mineral. We quantify these two effects by performing numerical simulations of subduction. Water transfers are computed as a function of slab dehydration/overlying mantle hydration reactions, and a strength decrease is imposed for hydrated mantle rocks. The model starts with an old oceanic plate ( 100 Ma) subducting for 145.5 Myr with a 5 cm/yr convergence rate. The arc lithosphere is thermally thinned between 100 km and 310 km away from the trench, due to small-scale convection occuring in the water-saturated mantle wedge. We test the influence of convergence slowdown by carrying on subduction with a decreased convergence rate (≤ 2 cm/yr). Surprisingly, the subduction slowdown yields not only a strong slab warming at great depth (> 80 km), but also a significant cooling of the forearc lithosphere at shallower depth. The convergence slowdown increases the subducted crust temperature at 90 km depth to 705 ± 62 °C, depending on the convergence rate reduction, and might thus favor the oceanic crust partial melting in presence of water. For subduction velocities ≤ 1 cm/yr, slab breakoff is triggered 20–32 Myr after slowdown onset, due to a drastic slab thermal weakening in the vicinity of the interplate plane base. At last, the rupture of the weakened forearc is simulated by imposing in the thinnest part of the overlying lithosphere a dipping weakness plane. For convergence with rates ≥ 1 cm/yr, the thinned forearc first shortens, then starts subducting along the slab surface. The forearc lithosphere subduction stops the slab surface warming by hot asthenosphere corner flow, and decreases in a first stage the slab surface temperature to 630 ± 20 °C at 80 km depth, in agreement with P–T range inferred from natural records of UHP-LT metamorphism. The subducted crust temperature is further reduced to 405 ± 10 °C for the crust directly buried below the subducting forearc. Such a cold thermal state at great depth has never been sampled in collision zones, suggesting that forearc subduction might not be always required to explain UHP-LT metamorphsim.     

Late Eocene to Early Miocene Andean uplift inferred from detrital zircon fission track and U–Pb dating of Cenozoic forearc sediments (15–18°S)

Timing, amount, and mechanisms of uplift in the Central Andes have been a matter of debate in the last decade. Our study is based on the Cenozoic Moquegua Group deposited in the forearc basin between the Western Cordillera and the Coastal Cordillera in southern Peru from ∼50 to ∼4 Ma. The Moquegua Group consists mainly of mud-flat to fluvial siliciclastic sediments with upsection increasing grain size and volcanic intercalations. Detrital zircon U–Pb dating and fission track thermochronology allow us to refine previous sediment provenance models and to constrain the timing of Late Eocene to Early Miocene Andean uplift. Uplift-related provenance and facies changes started around 35 Ma and thus predate major voluminous ignimbrite eruptions that started at ∼25 by up to 10 Ma. Therefore magmatic addition to the crust cannot be an important driving factor for crustal thickening and uplift at Late Eocene to Early Oligocene time. Changes in subduction regime and the subducting plate geometry are suggested to control the formation of significant relief in the area of the future Western Cordillera which acts as an efficient large-scale drainage divide between Altiplano and forearc from at least 15.5 to 19°S already at ∼35 Ma. The model integrates the coincidence of (i) onset of provenance change no later than 35 Ma, (ii) drastic decrease in convergence rates at ∼40, (iii) a flat-subduction period at around ∼40 to ∼30 Ma leading to strong interplate coupling, and (iv) strong decrease in volcanic activity between 45 and 30 Ma.     

Early Cretaceous continental arc-related volcanic rocks in the Duobuzha area,northern Tibet: implications for evolution history of the Bangong–Nujiang Ocean

New whole-rock major and trace elements data, zircon laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) U–Pb ages, and zircon Hf isotope compositions were analysed for Early Cretaceous volcanic rocks, also called Meiriqieco Formation (MF) in the Duobuzha area of the Southern Qiangtang–Baoshan Block (SQBB), northern Tibet. Our aim is to clarify their petrogenesis and tectonic setting, and constrain the evolution process on the northern margin of Bangong–Nujiang suture zone (BNSZ) during Early Cretaceous time. The MF volcanic rocks are mainly composed of andesites with subordinate basalts and rhyolites with high-K calc-alkaline affinity. Zircon LA-ICP-MS U–Pb dating for two andesite and one rhyolite samples give uniform ages within error of ca.113, 114, and 118 Ma, respectively, indicating they were erupted on the Early Cretaceous. The MF andesites have variable zircon ε_Hf(t) values (+0.5 to +10.5), which is different from those of MF rhyolites (+7.9 to +10.7). All the MF rocks are enriched in large ion lithophile elements, and depleted in high field strength elements, yielding the affinity of arc rocks. The MF basalts were most likely derived from the mantle wedge that was metasomatized by fluids released from subducting slab with the involvement of subducted sediments. The MF rhyolites were generated by partial melting of the juvenile mafic lower crust. The MF andesites are interpreted to have formed by mixing of the magmas that parental of the MF basalts and the MF rhyolites. In addition, a couple of distinctly magmatic sources are identified in the SQBB, and this may be related to mantle components injected into the continental crust. Combined with published geological data in the BNSZ and SQBB, we consider that the MF volcanic rocks are formed in a continental arc setting, suggesting that BNO were subducting during the Early Cretaceous time in the Duobuzha area.     

Zircon U‐Pb Ages and Geochemistry of Permo‐Carboniferous Mafic Intrusions in the Xilinhot Area,Inner Mongolia: Constraints on the Northward Subduction of the Paleo‐Asian Ocean

The Central Asian Orogenic Belt(CAOB) resulted from accretion during the Paleozoic subduction of the PaleoAsian Ocean. The Xilinhot area in Inner Mongolia is located in the northern subduction zone of the central-eastern CAOB and outcropped a large number of late Paleozoic mafic intrusions. The characteristics of magma source and tectonic setting of the mafic intrusions and their response to the closure process of the Paleo-Asian Ocean are still controversial. This study presents LA-ICPMS zircon U-Pb ages and geochemical features of mafic intrusions in the Xilinhot area to constrain the northward subduction of the Paleo-Asian Ocean. The mafic intrusions consist of gabbro, hornblende gabbro, and diabase. Their intrusion times can be divided into three stages of 326–321 Ma, 276 Ma and 254 Ma by zircon U-Pb ages. The first two stages of the 326–276 Ma intrusions mostly originated from subduction-modified continental lithospheric mantle sources that underwent a variable degree partial melting(5–30%), recording the subduction of oceanic crust. The third stage of the 254 Ma mafic rocks also show arc-related features. The primary magma compositions calculated by PRIMELT2 modeling on three samples of ～326 Ma and two samples of ～254 Ma show that these mafic samples are characterized by a variable range in SiO_2(47.51–51.47 wt%), Al_2O_3(11.46–15.55 wt%), ΣFeO(8.27–9.61 wt%), MgO(13.01–15.18 wt%) and CaO(9.13–11.67 wt%), consisting with the features between enriched mantle and lower continental crust. The source mantle melting of mafic intrusions occurred under temperatures of 1302–1351°C and pressures of 0.92–1.30 GPa. The magmatic processes occurred near the crust-mantle boundary at about 33–45 km underground. Combined with previous studies, it is concluded that Carboniferous to early Permian(～326–275 Ma) northward subduction of the Paleo-Asian oceanic crust led to the formation of the mafic magmatism in the Baolidao arc zone. The whole region had entered the collision environment at ～254 Ma, but with subduction-related environments locally. The final collision between the North China craton and the South Mongolian microcontinent may have lasted until ca. 230 Ma.     

Effect of time-evolving age and convergence rate of the subducting plate on the Cenozoic adakites and boninites

Partial melting of subducting oceanic crust expressed as high-Mg volcanic rocks such as adakites and boninites has been actively studied for decades, and Lee and King (2010) reported that time-evolving subduction parameters such as the age and the subduction rate of the converging oceanic plate play important roles in transient partial melting of the subducting oceanic crust (e.g., Aleutians). However, few subduction model experiments have considered time-evolving subduction parameters, posing problems for studies of transient partial melting of subducting oceanic crust in many subduction zones. Therefore, we constructed two-dimensional kinematic–dynamic subduction models for the Izu–Bonin, Mariana, Northeast Japan, Kuril, Tonga, Java–Sunda, and Aleutian subduction zones that account for the last 50 Myr of their evolution. The models include the time-evolving age and convergence rate of the incoming oceanic plate, so the effect of time-evolving subduction parameters on transient partial melting of oceanic crust can be evaluated. Our model calculations revealed that adakites and boninites in the Izu–Bonin and Aleutian subduction zones resulted from transient partial melting of oceanic crust. However, the steady-state subduction model using current subduction parameters did not produce any partial melting of oceanic crust in the aforementioned subduction zones, indicating that time-evolving subduction parameters are crucial for modeling transient eruption of adakites and boninites. Our model calculations confirm that other geological processes such as forearc extension, back-arc opening, mantle plumes and ridge subduction are required for partial melting of the oceanic crust in the Mariana, Northeast Japan, Tonga, and southeastern Java–Sunda subduction zones.     

Origin of Tertiary to Recent EM- and subduction-like chemical and isotopic signatures in Auca Mahuida region (37°–38°S) and other Patagonian plateau lavas

The alkaline volcanic rocks of the 1.8–0.9 Ma Auca Mahuida and post-mid-Pliocene Rio Colorado backarc volcanic fields east of the Andean Southern Volcanic Zone at ~37°–38°S have pronounced intraplate-like chemical signatures with some striking similarities to oceanic DM-EM1-like lavas of the south Atlantic Tristan da Cunha type. These backarc lavas are considered to have formed as a series of mantle batches typified by 4–7 % melting, with decompression melting initiating in a garnet-bearing mantle above a steepening subduction zone, and final equilibration occurring near the base of a ~65- to 70-km-thick lithosphere at temperatures of ~1,350–1,380 °C. Evolved Auca Mahuida mugearite to trachytic samples are best explained by crystal fractionation with limited mixing of partial melts of recently underplated basalts, in line with isotopic signatures that preclude significant radiogenic contamination in a preexisting refractory crust. Higher Ba/La and subtly higher La/Ta ratios than in nearby ~24–20 Ma primitive basalts or oceanic (OIB) lavas are attributed to the residual effects of slab fluids introduced during a shallow subduction episode recorded in the arc-like chemistry of the adjacent 7–4 Ma Chachahuén volcanic complex. Positive Sr, K and Ba spikes on mantle-normalized patterns of both primitive Auca Mahuida and ~24–20 Ma basalts, like those in EM-like OIB basalts, are attributed to mixing of continental lithosphere into the asthenosphere. In Patagonia, this mixing is suggested to have peaked as the South America continent accommodated to major late Oligocene plate convergence changes, as similar Sr, K and Ba spikes and DM-EM1 signatures are absent in ~50–30 Ma backarc lavas north of 51°S, and all of those south of 51°S. Introduction of an EM1-like component associated with lateral mantle flow of a Tristan da Cunha source is largely precluded by its Cretaceous age and distance to Patagonia.     

The dynamics of big mantle wedge,magma factory,and metamorphic–metasomatic factory in subduction zones

The East Asian continental margin is underlain by stagnant slabs resulting from subduction of the Pacific plate from the east and the Philippine Sea plate from the south. We classify the upper mantle in this region into three major domains: (a) metasomatic–metamorphic factory (MMF), subduction zone magma factory (SZMF), and the ‘big mantle wedge’ (BMW). Whereas the convection pattern is anticlockwise in the MMF domain, it is predominantly clockwise in the SZMF and BMW, along a cross section from the south. Here we define the MMF as a small wedge corner which is driven by the subducting Pacific plate and dominated by H₂O-rich fluids derived by dehydration reactions, and enriched in large ion lithophile elements (LILE) which cause the metasomatism. The SZMF is a zone intermediate between MMF and BMW domains and constitutes the main region of continental crust production by partial melting through wedge counter-corner flow. Large hydrous plume generated at about 200 km depth causes extensive reduction in viscosity and the smaller scale hydrous plumes between 60 km and 200 km also bring about an overall reduction in the viscosity of SZMF. More fertile and high temperature peridotites are supplied from the entrance to this domain. The domain extends obliquely to the volcanic front and then swings back to the deep mantle together with the subducting slab. The BMW occupies the major portion of upper mantle in the western Pacific and convects largely with a clockwise sense removing the eastern trench oceanward. Sporadic formation of hydrous plume at the depth of around 410 km and the curtain flow adjacent to the trench cause back arc spreading. We envisage that the heat source in BMW could be the accumulated TTG (tonalite–trondhjemite–granodiorite) crust on the bottom of the mantle transition zone. The ongoing process of transportation of granitic crust into the mantle transition zone is evident from the deep subduction of five intra-oceanic arcs on the subducting Philippine Sea plate from the south, in addition to the sediment trapped subduction by the Pacific plate and Philippine Sea plate. The dynamics of MMF, SZMF and BMW domains are controlled by the angle of subduction; a wide zone of MMF in SW Japan is caused by shallow angle subduction of the Philippine Sea plate and the markedly small MMF domain in the Mariana trench is due to the high angle subduction of Pacific plate. The domains in NE Japan and Kyushu region are intermediate between these two. During the Tertiary, a series of marginal basins were formed because of the nearly 2000 km northward shift of the subduction zone along the southern margin of Tethyan Asia, which may be related to the collision of India with Asia and the indentation. The volume of upper mantle under Asia was reduced extensively on the southern margin with a resultant oceanward trench retreat along the eastern margin of Asia, leading to the formation of a series of marginal basins. The western Pacific domain in general is characterized by double-sided subduction; from the east by the oldest Pacific plate and from the south by the oldest Indo-Australian plate. The old plates are hence hydrated extensively even in their central domains and therefore of low temperature. The cracks have allowed the transport of water into the deeper portions of the slab and these domains supply hydrous fluids even to the bottom of the upper mantle. Thus, a fluid dominated upper mantle in the western Pacific drives a number of microplates and promote the plate boundary processes.     

Incipient boninitic arc crust built on denudated mantle: the Khantaishir ophiolite (western Mongolia)

The ~ 570 Ma old Khantaishir ophiolite is built by up to 4 km harzburgitic mantle with abundant pyroxenites and dunites followed by ~ 2 km of hornblende-gabbros and gabbronorites and by a ~ 2.5 km thick volcanic unit composed of a dyke + sill complex capped by pillow lavas and some volcanoclastics. The volcanics are mainly basaltic andesites and andesites (or boninites) with an average of 58.2 ± 1.0 wt% SiO₂, X _Mg = 0.61 ± 0.03 (X _Mg = molar MgO/(MgO + FeO^tot), TiO₂ = 0.4 ± 0.1 wt% and CaO = 7.5 ± 0.6 wt% (errors as 2σ). Normalized trace element patterns show positive anomalies for Pb and Sr, a negative Nb-anomaly, large ion lithophile elements (LILE) concentrations between N- and E-MORB and distinctly depleted HREE. These characteristics indicate that the Khantaishir volcanics were derived from a refractory mantle source modified by a moderate slab-component, similar to boninites erupted along the Izu-Bonin-Mariana subduction system and to the Troodos and Betts Cove ophiolites. Most strikingly and despite almost complete outcrops over 260 km², there is no remnant of any pre-existing MORB crust, suggesting that the magmatic suite of this ophiolite formed on completely denudated mantle, most likely upon subduction initiation. The architecture of this 4–5 km thick early arc crust resembles oceanic crust formed at mid ocean ridges, but lacks a sheeted dyke complex; volcanic edifices are not observed. Nevertheless, low melting pressures combined with moderate H₂O-contents resulted in high-Si primitive melts, in abundant hornblende-gabbros and in a fast enrichment in bulk SiO₂. Fractional crystallization modeling starting from the observed primitive melts (56.6 wt% SiO₂) suggests that 25 wt% pyroxene + plagioclase fractionation is sufficient to form the average Khantaishir volcanic crust. Most of the fractionation happened in the mantle, the observed pyroxenite lenses and layers in and at the top of the harzburgites account for the required cumulate volumes. Finally, the multiply documented occurrence of highly depleted boninites during subduction initiation suggests a causal relationship of subduction initiation and highly depleted mantle. Possibly, a discontinuity between dense fertile and buoyant depleted mantle contributes to the sinking of the future dense subducting plate, while the buoyant depleted mantle of the future overriding plate forms the infant mantle wedge.     

Enrichments of the mantle sources beneath the Southern Volcanic Zone (Andes) by fluids and melts derived from abraded upper continental crust

Mafic basaltic-andesitic volcanic rocks from the Andean Southern Volcanic Zone (SVZ) exhibit a northward increase in crustal components in primitive arc magmas from the Central through the Transitional and Northern SVZ segments. New elemental and Sr–Nd-high-precision Pb isotope data from the Quaternary arc volcanic centres of Maipo (NSVZ) and Infernillo and Laguna del Maule (TSVZ) are argued to reflect mainly their mantle source and its melting. For the C-T-NSVZ, we identify two types of source enrichment: one, represented by Antuco in CSVZ, but also present northward along the arc, was dominated by fluids which enriched a pre-metasomatic South Atlantic depleted MORB mantle type asthenosphere. The second enrichment was by melts having the characteristics of upper continental crust (UCC), distinctly different from Chile trench sediments. We suggest that granitic rocks entered the source mantle by means of subduction erosion in response to the northward increasingly strong coupling of the converging plates. Both types of enrichment had the same Pb isotope composition in the TSVZ with no significant component derived from the subducting oceanic crust. Pb–Sr–Nd isotopes indicate a major crustal compositional change at the southern end of the NSVZ. Modelling suggests addition of around 2 % UCC for Infernillo and 5 % for Maipo.     

地壳与弱化岩石圈地幔的相互作用:以燕山造山带为例

燕山造山带中生代发育4期钙碱性火山活动,它们的源区组成都是受壳幔相互作用的制约,其中髫髻山组和义县组分布广泛,具有代表性．髫髻山组岩性比较单一,地球化学参数变化范围小,岩浆的AFC作用不强烈,源区成分不复杂．依据Kay et al．(1991)的方法,估算了早-中侏罗世燕山地区的地壳厚度为40-45 km.髫髻山组粗安岩是在加厚的地壳 (40-45 km)条件下,源区是含角闪石的石榴石麻粒岩 底侵的基性岩的壳幔过渡带熔融形成．义县组火山岩的源区为下地壳 岩石圈地幔,地幔组分较髫髻山组增加．研究区中生代早期地壳开始加厚,发生下地壳拆沉,进入流变学性质改变了的“弱化的岩石圈地幔”,二者发生作用．岩石圈地幔在中生代晚期受到流体、熔体、地幔矿物中活化的分子水、剪切构造作用,以及温、压条件改变的影响,导致岩石圈地幔发生不均一的局部弱化,为容纳拆沉的下地壳提供了优化场所．推测弱化岩石圈地幔出现于135 Ma以后燕山地区发育的小型拉伸盆地之下,以及对应的小型软流圈底辟体之上．上述模型可以与俯冲带的楔形地幔与俯冲洋壳的相互作用相对比．     

岛弧火山岩成因和地球化学特征

岛弧火山岩主要为俯冲带的俯冲板片脱水形成的富大离子亲石元素流体交代地幔楔,并使其发生部分熔融,产生岛弧岩浆作用而形成的,岩石组合通常为玄武岩—安山岩—英安岩—流纹岩及相应侵入岩组合。它以Al2O3、K2O高,低Ti O2,且K2ONa2O为特征,相对富集LILE,亏损HFSE,特别是Ti、Nb、Ta等。本文主要从岛弧岩浆作用的起因着手,分析流体和熔体对地幔楔的交代作用,以及岛弧岩浆作用过程,进而分析岛弧火山岩的地球化学特征。     

西太平洋雅浦海沟地幔演化与岩浆作用研究进展

俯冲带地幔演化与岩浆作用是地球各固体圈层之间发生物质和能量交换的重要地质过程.西太平洋雅浦海沟因其极短的沟-弧距离和洋脊碰撞等独特的地质构造特征成为研究复杂条件下俯冲带演化的理想场所.为了探究雅浦海沟地幔演化与岩浆作用,本文将前人对雅浦海沟火成岩的研究数据进行整合,分析了雅浦海沟火成岩的成因,并根据火成岩形成的制约条件,对卡罗琳板块俯冲到菲律宾海板块的地幔演化与岩浆作用过程进行了讨论.结果显示雅浦海沟火成岩均具有与俯冲相关火成岩的典型特征.橄榄岩地球化学特征指示雅浦海沟地幔熔融程度为20%~25%,地幔在部分熔融过程中受到了流体与熔体的双重交代作用.Re-Os同位素特征指示雅浦海沟地幔中存在约1.16 Ga非常古老的残余地幔,表明地幔可能经历过多期熔融事件,从而导致雅浦海沟地幔非常亏损.雅浦岛弧成因至今仍存争议,主要包括:(1)现今雅浦岛弧为帕里希维拉海盆洋壳的一部分,在中新世因卡罗琳洋脊的碰撞导致帕里希维拉海盆洋壳逆冲到原雅浦岛弧之上.(2)雅浦岛弧在不同构造时期经历过多期岛弧岩浆作用,包括俯冲初始阶段(~52 Ma)的弧前玄武岩、俯冲开始后的岛弧玄武岩(~25 Ma)、与卡罗琳洋脊碰撞(21 Ma)后的岛弧拉斑玄武岩(7~11 Ma).其中7~11 Ma的岛弧拉斑玄武岩指示雅浦岛弧岩浆活动并未因卡罗琳洋脊的碰撞完全停止,很有可能在晚中新世短暂恢复活动.      

High- and low-Cr chromitite and dunite in a Tibetan ophiolite: evolution from mature subduction system to incipient forearc in the Neo-Tethyan Ocean

The microstructures, major- and trace-element compositions of minerals and electron backscattered diffraction (EBSD) maps of high- and low-Cr# [spinel Cr# = Cr³⁺/(Cr³⁺ + Al³⁺)] chromitites and dunites from the Zedang ophiolite in the Yarlung Zangbo Suture (South Tibet) have been used to reveal their genesis and the related geodynamic processes in the Neo-Tethyan Ocean. The high-Cr# (0.77–0.80) chromitites (with or without diopside exsolution) have chromite compositions consistent with initial crystallization by interaction between boninitic magmas, harzburgite and reaction-produced magmas in a shallow, mature mantle wedge. Some high-Cr# chromitites show crystal-plastic deformation and grain growth on previous chromite relics that have exsolved needles of diopside. These features are similar to those of the Luobusa high-Cr# chromitites, possibly recycled from the deep upper mantle in a mature subduction system. In contrast, mineralogical, chemical and EBSD features of the Zedang low-Cr# (0.49–0.67) chromitites and dunites and the silicate inclusions in chromite indicate that they formed by rapid interaction between forearc basaltic magmas (MORB-like but with rare subduction input) and the Zedang harzburgites in a dynamically extended, incipient forearc lithosphere. The evidence implies that the high-Cr# chromitites were produced or emplaced in an earlier mature arc (possibly Jurassic), while the low-Cr# associations formed in an incipient forearc during the initiation of a new episode of Neo-Tethyan subduction at ~130–120 Ma. This two-episode subduction model can provide a new explanation for the coexistence of high- and low-Cr# chromitites in the same volume of ophiolitic mantle.     

大陆俯冲隧道板片-地幔楔界面反向流体交代作用:西阿尔卑斯造山带超高压变质白片岩的地球化学证据

前人对深俯冲板片释放熔/流体交代地幔楔形成弧岩浆源区的过程和机制已得到充分认识,然而对地幔楔岩石能否脱水交代深俯冲地壳并不清楚.在对欧洲西阿尔卑斯造山带Dora-Maira地体白片岩的地球化学研究中,首次发现地幔楔交代岩能够脱水反向交代深俯冲地壳岩石,从而极大影响俯冲地壳的地球化学组成.结合白片岩和围岩的全岩地球化学特征以及锆石学结果,查明了白片岩的原岩为S型花岗岩,澄清了关于Dora-Maira白片岩原岩属性的长期争议.在此基础上,发现白片岩中变质锆石相对残留岩浆锆石δ18O值显著降低,指示原岩形成后受到低δ18O变质流体的交代作用.白片岩具有高温岩石中最高的δ26Mg值达0.75‰,显著高于围岩变花岗岩,指示交代流体具有重Mg同位素组成.基于地球主要岩石储库的Mg同位素组成,推测交代流体来自俯冲隧道中富滑石地幔楔蛇纹岩在弧下深度的脱水分解,而这些地幔楔蛇纹岩是新特提斯洋壳在弧前深度变质脱水产生的流体与地幔楔浅部橄榄岩反应形成.这些结果不仅提供了利用Mg-O同位素示踪俯冲隧道中流体来源的新思路,也提供了地幔楔蛇纹岩来源流体反向交代深俯冲地壳岩石的首个典型实例.这种反向交代不仅极大改变了深俯冲地壳的地球化学组成,而且对弧岩浆岩重Mg同位素成因具有重要意义.      

硼的地球化学性质及其在俯冲带的循环与成矿初探

硼是广泛应用于化工、农业、材料科学及核工业领域的重要元素。硼与氢的核聚变反应是未来具备运用潜力的清洁能源。硼作为典型的亲石元素,是高度不相容元素。硼元素容易富集于蚀变洋壳及蛇纹石化地幔橄榄岩中。而在板块俯冲过程中,由于硼具有强的流体活动性,会优先赋存于流体中。因此,当蛇纹石化的大洋岩石圈及覆于其上的沉积物在俯冲过程中发生脱水,这使得弧前地幔楔发生大规模的蛇纹石化。此时大量硼元素很可能随俯冲流体释放并封存于弧前地幔楔中。目前已发现的超大型硼矿床主要位于聚合型板块边缘,尤其土耳其拥有世界上最大的硼酸盐储量。我们推测这些矿床的形成基础条件之一可能与弧前高度蛇纹石化的地幔楔有关。尤其是在洋 陆俯冲环境,弧前蛇纹岩或蛇绿混杂岩首先通过俯冲侵蚀再循环到火山弧岩浆中,使得岩浆更富集硼。随后弧火山喷发大量富硼的火山岩、岩浆热液及水气。在岩浆冷却过程中,硼元素析出、沉淀于火山表面,并伴随风化、侵蚀过程汇聚至碰撞造山带的封闭湖盆之中。此外,干冷的气候条件下也进一步促进了硼的成矿。我国具有形成大型、超大型硼矿的地质条件,应加大研究及探勘力度,并适当购买硼作为战略储备。