Magma Extraction from the Mantle Wedge at Convergent Margins Through Dikes: A Parametric Sensitivity Analysis

We address under what conditions a magma generated by partial melting at 100 km depth in the mantle wedge above a subduction zone can reach the crust in dikes before stalling. We also address under what conditions primitive basaltic magma (Mg # >60) can be delivered from this depth to the crust. We employ linear elastic fracture mechanics with magma solidification theory and perform a parametric sensitivity analysis. All dikes are initiated at a depth of 100 km in the thermal core of the wedge, and the Moho is fixed at 35 km depth. We consider a range of melt solidus temperatures (800-1100 degrees C), viscosities (10-100 Pa s), and densities (2400-2700 kg m(-3)). We also consider a range of host rock fracture toughness values (50-300 MPa m(1/2)) and dike lengths (2-5 km) and two thermal structures for the mantle wedge (1260 and 1400 degrees C at 100 km depth and 760 and 900 degrees C at 35 km depth). For the given parameter space, many dikes can reach the Moho in less than a few hundred hours, well within the time constraints provided by U series isotope disequilibria studies. Increasing the temperature in the mantle wedge, or increasing the dike length, allows additional dikes to propagate to the Moho. We conclude that some dikes with vertical lengths near their critical lengths and relatively high solidus temperatures will stall in the mantle before reaching the Moho, and these may be returned by corner flow to depths where they can melt under hydrous conditions. Thus, a chemical signature in arc lavas suggesting partial melting of slab basalts may be partly influenced by these recycled dikes. Alternatively, dikes with lengths well above their critical lengths can easily deliver primitive magmas to the crust, particularly if the mantle wedge is relatively hot. Dike transport remains a viable primary mechanism of magma ascent in convergent tectonic settings, but the potential for less rapid mechanisms making an important contribution increases as the mantle temperature at the Moho approaches the solidus temperature of the magma.

Recent contribution of sediments and fluids to the mantle's volatile budget

Subduction modifies the cycling of Earth's volatile elements. Fluid-rich sediments and hydrated oceanic lithosphere enter the convecting mantle at subduction zones. Some of the sediments and volatile components are released from the subducting slab, promote mantle melting and are returned to the surface by volcanism. The remainder continue into the deeper mantle. Quantification of the fate of these volatiles requires an understanding of both the nature and timing of fluid release and mantle melting(1). Here we analyse the trace element and isotopic geochemistry of fragments of upper mantle rocks that were transported to the surface by volcanic eruptions above the Batan Island subduction zone, Philippines. We find that the mantle fragments exhibit extreme disequilibrium between their U-Th-Ra isotopic ratios, which we interpret to result from the interaction of wet sediment melts and slab-derived fluids with rocks in the overlying mantle wedge. We infer that wet sediments were delivered from the slab to the mantle wedge between 8,000 and 10,000 years ago, whereas aqueous fluids were delivered separately much later. We estimate that about 625 ppm of water is retained in the wedge. A significant volume of water could therefore be delivered to the mantle transition zone at the base of the upper mantle, or even to the deeper mantle.

Fluid-related inclusions in Alpine high-pressure peridotite reveal trace element recycling during subduction-zone dehydration of serpentinized mantle (Cima di Gagnone, Swiss Alps).

Serpentinites release at sub-arc depths volatiles and several fluid-mobile trace elements found in arc magmas. Constraining element uptake in these rocks and defining the trace element composition of fluids released upon serpentinite dehydration can improve our understanding of mass transfer across subduction zones and to volcanic arcs. The eclogite-facies garnet metaperidotite and chlorite harzburgite bodies embedded in paragneiss of the subduction melange from Cima di Gagnone derive from serpentinized peridotite protoliths and are unique examples of ultramafic rocks that experienced subduction metasomatism and devolatilization. In these rocks, metamorphic olivine and garnet trap polyphase inclusions representing the fluid released during high-pressure breakdown of antigorite and chlorite. Combining major element mapping and laser-ablation ICP-MS bulk inclusion analysis, we characterize the mineral content of polyphase inclusions and quantify the fluid composition. Silicates, Cl-bearing phases, sulphides, carbonates, and oxides document post-entrapment mineral growth in the inclusions starting immediately after fluid entrapment. Compositional data reveal the presence of two different fluid types. The first (type A) records a fluid prominently enriched in fluid-mobile elements, with Cl, Cs, Pb, As, Sb concentrations up to 10(3) PM (primitive mantle), similar to 10(2) PM Tit Ba, while Rb, B, Sr, Li, U concentrations are of the order of 10(1) PM, and alkalis are similar to 2 PM. The second fluid (type B) has considerably lower fluid-mobile element enrichments, but its enrichment patterns are comparable to type A fluid. Our data reveal multistage fluid uptake in these peridotite bodies, including selective element enrichment during seafloor alteration, followed by fluid-rock interaction along with subduction metamorphism in the plate interface melange. Here, infiltration of sediment-equilibrated fluid produced significant enrichment of the serpentinites in As, Sb, B, Pb, an enriched trace element pattern that was then transferred to the fluid released at greater depth upon serpentine dehydration (type A fluid). The type B fluid hosted by garnet may record the composition of the chlorite breakdown fluid released at even greater depth. The Gagnone study-case demonstrates that serpentinized peridotites acquire water and fluid-mobile elements during ocean floor hydration and through exchange with sediment-equilibrated fluids in the early subduction stages. Subsequent antigorite devolatilization at subarc depths delivers aqueous fluids to the mantle wedge that can be prominently enriched in sediment-derived components, potentially triggering arc magmatism without the need of concomitant dehydration/melting of metasediments or altered oceanic crust.

Focused fluid transfer through the mantle above subduction zones

Volcanic arcs above subduction zones are enriched in volatiles and fluid-mobile elements with respect to mid-oceanic ridge basalts. There is general consensus that this particular subduction zone signature is generated by fluid-induced extraction of these elements from subducted oceanic crust and its sedimentary cover. However, how these fluids are transferred through the mantle wedge to the locus of partial melting and what modification the fluids will experience is unresolved. Here we investigate the interaction of slab fluids with the mantle wedge through a series of high-pressure experiments. We explore two end-member processes of focused and porous reactive flow of hydrous slab melts through the mantle. Transfer by porous flow leads to the formation of hydrous minerals that sequester fluid-mobile elements and residual fluids characterized by trace element patterns inconsistent with typical arc lavas. In contrast, no hydrous minerals are formed in the reaction zone of experiments mimicking focused flow, and the typical trace element signature acquired during fluid extraction from the slab is preserved, indicating that this is an efficient process for element transfer through the mantle wedge.

Tracking the magmatic evolution of island arc volcanism : insights from a high-precision Pb isotope record of Montserrat, Lesser Antilles

The volcanic succession on Montserrat provides an opportunity to examine the magmatic evolution of island arc volcanism over a ∼2.5 Ma period, extending from the andesites of the Silver Hills center, to the currently active Soufrière Hills volcano (February 2010). Here we present high-precision double-spike Pb isotope data, combined with trace element and Sr-Nd isotope data throughout this period of Montserrat's volcanic evolution. We demonstrate that each volcanic center; South Soufrière Hills, Soufrière Hills, Centre Hills and Silver Hills, can be clearly discriminated using trace element and isotopic parameters. Variations in these parameters suggest there have been systematic and episodic changes in the subduction input. The SSH center, in particular, has a greater slab fluid signature, as indicated by low Ce/Pb, but less sediment addition than the other volcanic centers, which have higher Th/Ce. Pb isotope data from Montserrat fall along two trends, the Silver Hills, Centre Hills and Soufrière Hills lie on a general trend of the Lesser Antilles volcanics, whereas SSH volcanics define a separate trend. The Soufrière Hills and SSH volcanic centers were erupted at approximately the same time, but retain distinctive isotopic signatures, suggesting that the SSH magmas have a different source to the other volcanic centers. We hypothesize that this rapid magmatic source change is controlled by the regional transtensional regime, which allowed the SSH magma to be extracted from a shallower source. The Pb isotopes indicate an interplay between subduction derived components and a MORB-like mantle wedge influenced by a Galapagos plume-like source.

Pulling apart the mid to late Cenozoic magmatic record of the Gulf of California : no room for the Comondú Arc

Understanding the link between tectonic-driven extensional faulting and volcanism is crucial from a hazard perspective in active volcanic environments, while ancient volcanic successions provide records on how volcanic eruption styles, compositions, magnitudes and frequencies can change in response to extension timing, distribution and intensity. Significantly, incorrect tectonic interpretations can be made when the spatial-temporal-compositional trends of, and source contributions to magmatism are not properly considered. This study draws on intimate relationships of volcanism and extension preserved in the Sierra Madre Occidental (SMO) and Gulf of California (GoC) regions of western Mexico. Here, a major Oligocene rhyolitic ignimbrite “flare-up” (>300,000 km3) switched to a dominantly bimodal and mixed effusive-explosive volcanic phase in the Early Miocene (~100,000 km3), associated with distributed extension and opening of numerous grabens. Rhyolitic dome fields were emplaced along graben edges and at intersections of cross-graben and graben-parallel structures during early stages of graben development. Concomitant with this change in rhyolite eruption style was a change in crustal source as revealed by zircon chronochemistry with rapid rates of rhyolite magma generation due to remelting of mid- to upper crustal, highly differentiated igneous rocks emplaced during earlier SMO magmatism. Extension became more focused ~18 Ma resulting in volcanic activity being localised along the site of GoC opening. This localised volcanism (known as the Comondú “arc”) was dominantly effusive and andesite-dacite in composition. This compositional change resulted from increased mixing of basaltic and rhyolitic magmas rather than fluid flux melting of the mantle wedge above the subducting Guadalupe Plate. A poor understanding of space-time relationships of volcanism and extension has thus led to incorrect past tectonic interpretations of Comondú-age volcanism.

Pulling apart the mid to late Cenozoic magmatic record of the Gulf of California : is there a Comondú arc?

The composition of the lithosphere can be fundamentally altered by long-lived subduction processes such that subduction-modified lithosphere can survive for 100's Myrs. Incorrect petrotectonic interpretations result when spatial-temporal-compositional trends of, and source contributions to, magmatism are not properly considered. Western Mexico has had protracted Cenozoic magmatism developed mostly in-board of active oceanic plate subduction beneath western North America. A broad range of igneous compositions from basalt to high-silica rhyolite were erupted with intermediate to silicic compositions in particular, showing calc-alkaline and other typical subduction-related geochemical signatures. A major Oligocene rhyolitic ignimbrite “flare-up” (>300,000 km3) switched to a bimodal volcanic phase in the Early Miocene (~100,000 km3), associated with distributed extension and opening of numerous grabens. Extension became more focussed ~18 Ma resulting in localised volcanic activity along the future site of the Gulf of California. This localised volcanism (known as the Comondú “arc”) was dominantly effusive and andesite-dacite in composition. Past tectonic interpretations of Comondú-age volcanism may have been incorrect as these regional temporal-compositional changes are alternatively interpreted as a result of increased mixing of mantle-derived basaltic and crust-derived rhyolitic magmas in an active rift environment rather than fluid flux melting of the mantle wedge above the subducting Guadalupe Plate.

Consequences of extensional tectonics on volcanic eruption style, compositions and source regions : new insights from the southern Sierra Madre Occidental and Gulf of California regions, western Mexico

Understanding the link between tectonic-driven extensional faulting and volcanism is crucial from a hazard perspective in active volcanic environments, while ancient volcanic successions provide records on how volcanic eruption styles, compositions, magnitudes and frequencies can change in response to extension timing, distribution and intensity. This study draws on intimate relationships of volcanism and extension preserved in the Sierra Madre Occidental (SMO) and Gulf of California (GoC) regions of western Mexico. Here, a major Oligocene rhyolitic ignimbrite “flare-up” (>300,000 km3) switched to a dominantly bimodal and mixed effusive-explosive volcanic phase in the Early Miocene (~100,000 km3), associated with distributed extension and opening of numerous grabens. Rhyolitic dome fields were emplaced along graben edges and at intersections of cross-graben and graben-parallel structures during early stages of graben development. Concomitant with this change in rhyolite eruption style was a change in crustal source as revealed by zircon chronochemistry with rapid rates of rhyolite magma generation due to remelting of mid- to upper crustal, highly differentiated igneous rocks emplaced during earlier SMO magmatism. Extension became more focused ~18 Ma resulting in volcanic activity being localised along the site of GoC opening. This localised volcanism (known as the Comondú “arc”) was dominantly effusive and andesite-dacite in composition. This compositional change resulted from increased mixing of basaltic and rhyolitic magmas rather than fluid flux melting of the mantle wedge above the subducting Guadalupe Plate. A poor understanding of space-time relationships of volcanism and extension has thus led to incorrect past tectonic interpretations of Comondú-age volcanism.

Neoarchean sanukitoid series in the Karelian Province, Finland

Sanukitoid series intrusions can be found throughout the Archean Karelian Province of the Fennoscandian shield. All sanukitoids share the same controversial elemental characteristics: they have high content of incompatible elements such as K, Ba, and Sr as well as high content of the compatible elements Mg, Cr, and Ni, and high Mg#. This composition is explained by an enriched mantle wedge origin in a Neoarchean subduction setting. This study concentrates on sanukitoid intrusions and tonalite-trondhjemite-granodiorite series (TTGs) from Finnish part of the Karelian Province. The collected rock samples have been studied in the field and under microscope as well as for their whole-rock (including isotopes) and mineral compositions. The new data together with previously published analyses help us to better understand the petrogenesis, tectonic setting and reworking of the Archean rock units. TTGs from the Karelian Province form a voluminous series of granitoids and reworked migmatites. This study divides TTG series into two subgroups based on their elemental composition: low-HREE (heavy rare earth element) TTGs and high-HREE TTGs indicating pressure differences in their source. Sanukitoid series is a minor, divergent group of intrusions. These intrusions are variable sized, and the texture varies from even-grained to K-feldspar porphyritic. The elemental composition differentiates sanukitoids from more voluminous TTG groups, the SiO2 in sanukitoids varies to include series of gabbro, diorite, and granodiorite. U Pb age determinations from sanukitoid series show temporally limited emplacement between ~ 2745 2715 Ma after the main crust forming period in the area. Hafnium, neodymium, common lead, and oxygene isotopes indicate well homogenized characteristics. Recycled crust has made a variable, yet minor, contribution to sanukitoids, as evidenced by oxygene isotopes and inherited zircon cores. A proposed tectonic setting for the formation of the sanukitoid series is slab breakoff of oceanic lithosphere in subduction setting, with sanukitoids deriving from an enriched mantle wedge. The proposed setting explains some of the peculiar features of sanukitoids, such as their temporally limited occurrence and controversial elemental composition. Sanukitoids would occur after cessation of the regional growth of Archean crust, and they could be derived from mantle wedge previously enriched by melts and fluids from oceanic crust and sediments. A subsequent event during the Paleoproterozoic Svecofennian orogeny at ~1.9 Ga affected the appearance and microstructures of the rocks as well as caused redistribution of lead between minerals and whole rock. However, the deformation was not able to obliterate the original geochemical characteristics of these sanukitoids.

Chromite-silicate chemistry of the Neoarchean Sittampundi Complex, southern India: Implications for subduction-related arc magmatism

The Neoarchean layered anorthositic complex at Sittampundi in southern India is known for its chromitite layers that are mostly associated with anorthosite (An(90-100)). The chromitites contain FeAl-rich chromites concentrated in layers between amphibole-rich layers with a dominant mineralogy of amphibole-spinel-plagiocase+/-sapphirine. The chromite-rich layers contain only amphibole and plagioclase. Mineral compositions illustrated by X-ray composition maps and profiles show subtle chemical differences. The chrome spinels are of refractory grade with Cr2O3 and Al2O3 contents varying between 34-40 wt.% and 23-28 wt.%. The chromite compositions are noticeably different from those in layered igneous intrusions of the Bushveld-Stillwater type. The existence of original highly calcic plagioclase, FeAl-rich chromite, and magmatic amphibole is consistent with derivation from a parental magma of hydrous tholeiitic composition that was most likely generated in a supra-subduction zone arc setting. In terms of mineralogy and field relations, the Sittampundi chromitites are remarkably similar to anorthosite-hosted chromitites in the Neoarchean Fiskensset anorthositic complex, Greenland. We propose that the Sittampundi chromitites formed by partial melting of unusually aluminous harzburgite in a hydrated mantle wedge above a subduction zone. This melting process produced hydrous, aluminous basalt, which fractionated at depth to give rise to a variety of high-alumina basalt compositions from which the anorthositic complex with its cumulate chromite-rich and amphibole-rich layers formed within the magma chamber of a supra-subduction zone arc. (C) 2011 Elsevier B.V. All rights reserved.

Imaging the earth with ambient noise and earthquakes

In this thesis, I develop the velocity and structure models for the Los Angeles Basin and Southern Peru. The ultimate goal is to better understand the geological processes involved in the basin and subduction zone dynamics. The results are obtained from seismic interferometry using ambient noise and receiver functions using earthquake- generated waves. Some unusual signals specific to the local structures are also studied. The main findings are summarized as follows:

(1) Los Angeles Basin

The shear wave velocities range from 0.5 to 3.0 km/s in the sediments, with lateral gradients at the Newport-Inglewood, Compton-Los Alamitos, and Whittier Faults. The basin is a maximum of 8 km deep along the profile, and the Moho rises to a depth of 17 km under the basin. The basin has a stretch factor of 2.6 in the center decreasing to 1.3 at the edges, and is in approximate isostatic equilibrium. This "high-density" (~1 km spacing) "short-duration" (~1.5 month) experiment may serve as a prototype experiment that will allow basins to be covered by this type of low-cost survey.

(2) Peruvian subduction zone

Two prominent mid-crust structures are revealed in the 70 km thick crust under the Central Andes: a low-velocity zone interpreted as partially molten rocks beneath the Western Cordillera – Altiplano Plateau, and the underthrusting Brazilian Shield beneath the Eastern Cordillera. The low-velocity zone is oblique to the present trench, and possibly indicates the location of the volcanic arcs formed during the steepening of the Oligocene flat slab beneath the Altiplano Plateau.

The Nazca slab changes from normal dipping (~25 degrees) subduction in the southeast to flat subduction in the northwest of the study area. In the flat subduction regime, the slab subducts to ~100 km depth and then remains flat for ~300 km distance before it resumes a normal dipping geometry. The flat part 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 western half of the flat slab, which indicates the lack of melting and thus explains the cessation of the volcanism above. The velocity turns to normal values before the slab steepens again, indicating possible resumption of dehydration and ecologitization.

(3) Some unusual signals

Strong higher-mode Rayleigh waves due to the basin structure are observed in the periods less than 5 s. The particle motions provide a good test for distinguishing between the fundamental and higher mode. The precursor and coda waves relative to the interstation Rayleigh waves are observed, and modeled with a strong scatterer located in the active volcanic area in Southern Peru. In contrast with the usual receiver function analysis, multiples are extensively involved in this thesis. In the LA Basin, a good image is only from PpPs multiples, while in Peru, PpPp multiples contribute significantly to the final results.

Trace element composition of peridotites from the southern Mariana forearc: Insights into the geochemical effects of serpentinization and/or seafloor weathering

Peridotites from the southern Mariana forearc were sampled on the landward trench slope of the Izu-Bonin-Mariana (IBM) subduction zone by dredging. These mantle wedge peridotites underwent hydration by fluid derived from a dehydrated descending slab, and later interacted with seawater after emplacement at or near the seafloor. This study investigates how these two different rock-fluid interaction processes influenced trace element distribution in the southern Mariana forearc peridotites. We measured trace element concentrations of peridotites from the southern Mariana forearc. The southern Mariana forearc peridotites are characterized by a distinct seawater-like REE pattern with an obvious negative Ce anomaly, and La shows good correlation with other REEs (except Ce). In addition, there is a great enrichment of U, Pb, Sr and Li elements, which show a distinct positive anomaly relative to adjacent elements in the multi-element diagram. For the seawater-like REE pattern, we infer that REEs are mainly influenced by seawater during peridotite-seawater interactions after their emplacement at or near the seafloor, by serpentinization or by marine weathering. Furthermore, the anomalous behavior of Ce, compared with other rare earth elements in these samples, may indicate that they have undergone reactions involving Ce (IV) when the peridotites interacted with seawater. Positive U, Pb, Sr and Li anomalies are inferred to be related to seawater and/or fluids released during dehydration of the subducting slab.

Geochemical properties of the variolitic andesite from Susong in the Dabie orogenic belt, China

The variolitic andesite from the Susong County in the Dabie Mountains implies that it was erupted in water. The mineralogy of the varioles is primarily radiate plagioclase (albite sind oligoclase), with little pyroxene, hornblende and quartz (derived from alteration). The pyroxene, hornblende and quartz are in the interstices between plagiocalse. The matrix consists of glass, hornblende, chlorite, epidote and zoisite. It is clearly subjected an extensive alteration. The andesite has an uncommon chemical composition. The SiO2 content is about 56.8%, TiO2 = 0.9%, MgO = 6.4%, Fe2O3 (tot) = 6.7%similar to 7.6%, 100Mg/(Mg + Fe) = 64.1 similar to 66.2. Mg-# is significantly high. The andesite has high abundances of large-lithophile trace elements (e.g. K, Ba. Sr, LREE), e.g. La/Nb = 5.56 similar to 6.07, low abundances of high-strength-field elements (HFSE e.g. Ta, Nb, P, Ti), particularly Ta and Nb strongly depleted. These are consistent with the characteristics of subduction-related magmas. In the spider diagram of trace elements, from Ce to right hand, the abundances of elements decrease quickly, showing a character of the continental margins. There has a strong punishment of light-rare-earth elements, with a significant diffraction of REEs (the mean value of (La/Yb)(N) is 32.84). No Eu anomaly, but there are anomaly high (La/Yb)(N) = 28.63 similar to 36.74, (La/Y)(N) = 70.33 similar to 82.4. The elements Y and Yb are depleted greatly, Y<20 g/g, Y-N = 2.74 similar to 2.84, Yb-N = 2.18 similar to 2.35. From the La-(La/Sm) diagram, the andesite is derived from partial melting. But the epsilone value of Nd is - 18.7 similar to -19.2, so that the material source may be the mantle materials affected by the crustal materials. The Nd model age is 1.9 Ga indicating that the variolitic basaltic andesite was resulted from the mantle wedge of North China block, which had the Nd model age of 2.5Ga, when the Yangze block which had the Nd model age of 1.7Ga subducted beneath it. So the variolitic andesite has characteristics of the island-are volconic rocks oil a continental basement in the vicinity of the destructive continental margin.