Diapiric Ascent Research Articles

Extreme Mantle Heterogeneity Revealed by Geochemical Investigation of In Situ Lavas at the Central Mohns Ridge, Arctic Mid‐Ocean Ridges

AbstractMid‐ocean ridge basalts reflect the mantle’s composition and reveal processes from melting to eruption. The Mohns and Knipovich Ridges have ultraslow spreading rates, low magma budgets and erupted lavas indicating various mantle domains. Here, we use geochemistry and isotope systematics of in situ samples from two axial volcanic ridges (AVRs) to study mantle heterogeneity and melt production. By linking chemical variations to high‐resolution bathymetry and age data, we document systematic changes over time in the mantle source of the volcanic sequence. At Mohns Ridge AVR‐M10 (72.3°N), we observed significant variations in chemistry (e.g., (La/Sm)N from 0.7 to 2.9) and isotope systematics in basaltic samples from a small area (∼1 km2), suggesting the emplacement of multiple small‐volume lava flows. Pb isotope variations, for example, 206Pb/204Pb (17.91–18.76), are comparable with the observed range along the entire Mohns and Knipovich Ridges. Temporal constraints document that erupted basalts have changed from highly radiogenic Pb compositions to a more depleted signature within 30 ka. To explain the extreme variations in the erupted lavas at the Mohns Ridge, the mantle would need to be highly heterogeneous in composition with effective melt extraction and limited mixing prior to eruption. We use the highly heterogenous mantle underneath the Mohns Ridge to understand the melt extraction processes and mixing of melts and propose a two‐stage melting model: continuous generation of enriched melts from a deep and fertile source in the first stage, while depleted melts from a shallower and more refractory mantle occur sporadically and simultaneously with the intermittent ascent of diapirs.

The fate of ultramafic-rich mélanges in cold to hot subduction zones: Implications for diapirism (or not) and chemical geodynamics

Buoyant ultramafic-rich (serpentine- or chlorite-rich) mélange diapirs in sediment-starved subduction zones can transport slab material to arc sources. While the buoyancy of chlorite-rich mélanges was previously investigated, serpentine-rich mélanges were never explored. Thus, the overall contribution of ultramafic-rich mélanges to buoyancy, the conditions for diapir formation, and their fate in subduction zones are not well constrained. Here, we investigate the partial melting behavior and the associated density transformations of a serpentine-rich matrix (5–10 wt.% H2O) with minor sediments (9:1 ratio) at fore-arc (∼65 km) to sub-arc (∼95 km) depths (2–3 GPa and 800–1250 °C) and compare to that of chlorite-rich mélanges from the literature. Our results show that the solidus of serpentine-rich matrices is between 1050 and 1100 °C and requires either diapiric rise of the mélange into the hotter mantle wedge or interactions with a hotter asthenosphere through slab tears to partially melt and produce basaltic melts, whether in hot or cold slab channels. Chlorite-rich mélanges may account for the sources of some arc lavas, but partial melting of serpentine-rich mélanges produce melts depleted in CaO, TiO2, alkalis, and are highly enriched in MgO compared to basaltic arc lavas. Both serpentine-rich and chlorite-rich matrices dehydrate to form denser peridotite and lose buoyancy at ∼800 °C and ≥1000 °C, respectively. Even if diapirism initiates in such mélanges near the slab-mantle interface, they would likely lose buoyancy upon ascent into the hotter mantle wedge resulting in stalled or failed diapirs. Diapir growth (τa) is controlled by the interplay of density, thickness and viscosity of the mélange, as well as the timescale of slab subduction (τs) and thermal structure of the subduction zone. We observe that the onset of diapirs in cold subduction zones requires mélanges that may sometimes be thicker than that observed by field and geophysical studies, while hot subduction zones overall require thinner mélanges. Thus, ultramafic-rich mélange diapirs may occur but only under specific conditions and when the diapiric ascent timescale is faster than the thermal equilibration barrier of ∼800–1000 °C (especially at the core of the mélange). Dehydration or partial melting of ultramafic-rich mélanges can affect the large ion lithophile element (LILE), volatiles, and high-field strength element (HFSE) budgets in the mantle wedge. Partial melting (caused by a diapiric rise or slab tear) does not fractionate LILEs from HFSEs at T ≥ 1100 °C and if the mélange has a lower LILE/HFSE to begin with, that signature is transferred to arc sources. Dehydration releases aqueous fluids rich in fluid-mobile elements (LILE and volatiles) relative to HFSE. Thus, the characteristic high LILE/HFSE signature of aqueous fluids is transferred to arc magma sources. Given high LILE/HFSE ratio is a ubiquitous arc magma signature, but slab tears are not, and diapirism in ultramafic-rich mélanges is highly conditional, this study corroborates that aqueous fluids released from sediment-starved mélanges are the predominant mass transfer agents rather than diapirs.

Serpentinite Diapirs and the Evolution of Oceanic Core Complexes

Serpentinites and peridotites are the predominant lithological components of Oceanic Core Complexes (OCCs), located commonly at triple junctions of slow-spreading oceanic accreting ridges, fracture zones and transform faults. These lithologies differ from the surrounding oceanic brittle lithosphere, built of basalt and gabbro, and the structural evolution of these OCCs is enigmatic. The present investigation suggests that the tectonics of OCCs is derived from the ascent of serpentinic diapirs generated by the unique proximity of shallow asthenosphere, faulted lithosphere and the juxtaposition of oceanic crusts of contrasting densities. Such setting initiated two structural stages in the evolution of the OCC, the first is spontaneous subduction of old and dense oceanic crust under the fresh and lighter basalt at the edge of the MOR across its intersection with transform fault - fracture zone. The sub ducted slab would be affected by the combined effect of the steep thermal gradient and the availability of volatiles there to enhance the alteration of pyroxene into serpentine. Analog and numeric experiments show that spontaneous subduction can initiate spontaneously if the density contrast between the juxtaposed slabs is significant, as is the case between fresh basalt (ρ~2.7 g/cm3 ) and older basalt (ρ~2.9 g/cm3 ). Since the average thermal gradient under the MOR is ca. 130o /km, then at depth of 4-5 km the pyroxene and plagioclase would start their alteration, mostly to serpentine. The abundant faults of MOR normal rifting and strike-slip transform faulting could enable the light and malleable serpentinite to form diapirs, which would ascent to the seafloor at the rift-fracture zone intersections. The friction at the top of the diapirs during their ascent would probably generate breccia at their tops. Highlights a) OCCs are built mostly of serpentinite diapirs that carry peridotite inclusions. b) Serpentinization occurs under low pressure-high temperature conditions. c) OCC source rock is slab of cold oceanic crust sub ducted under fresh crust. d) Such subduction occurs mostly where large density contrast between slabs exists. e) The link between OCCs and detachment normal faults requires support.

Mechanism of the historical and the ongoing Vulcanian eruptions of Ebeko volcano, Northern Kuriles.

Ebeko is one of the most active volcanoes of the Kurile island arc, producing frequent mild Vulcanian explosions with eruption clouds up to 5 km high. The volcano poses a serious threat to the Severo-Kurilsk town with a population of around 2500 inhabitants, located at a distance of only 7 km on a fan of the volcano’s laharic deposits. Here, we report an overview of the activity of the volcano in the 20th–21st centuries and the results of our geological and petrological investigations of the ongoing Vulcanian eruption that started in 2016. We have found that eruptions of Ebeko span a range of mechanisms from purely magmatic to phreatic/hydrothermal. Three of its historical eruptions (the 1934–1935, 1987–1991, and the 2016–ongoing) involved fresh magma, while during the others (1967–1971, 2009–2011) fresh magma was not erupted. Juvenile material of the ongoing eruption represents highly crystalline and highly viscous (more than 108 pa s) low-silica (56–58 wt% SiO2) andesite. Historical data and our observations of the ongoing eruption allowed us to suggest a functional model of the volcano where Vulcanian explosions are caused by shallow intrusions of small diapir-like batches of strongly crystallized and highly viscous andesitic magma ascending into water-saturated, hydrothermally altered rocks composing the volcano summit. We suggest that the diapir’s ascent is governed by their positive buoyancy. Some of the diapirs reach and breach the ground surface producing magmatic eruptions of Ebeko, while the others are stuck at the shallow subsurface level and feed intensive hydrothermal activity as well as phreatic eruptions of the volcano. Positive buoyancy of the diapirs is too weak to allow them to extrude high above the ground surface to form lava domes.Supplementary InformationThe online version contains supplementary material available at 10.1007/s00445-020-01426-z.

Structural evolution of continental and marine Permian rock salt of the North German Basin: constraints from microfabrics, geochemistry and U\u2013Pb ages

Analyzing the dynamics of microstructural response on natural deformation in rock salt, we present microfabric, EBSD, geochemical and U–Pb data, obtained from Permian salt formations of the Kiel-Honigsee salt wall in Northern Germany. The samples were recovered from deep drillings, which penetrated through an overturned rock salt sequence of both Rotliegend and Zechstein deposits. The bromide concentration in halite indicates a continental and marine origin for the Rotliegend and Zechstein deposits, respectively. Despite intense deformation, relics of early diagenetic fabrics are still preserved. Deformation of the impure Rotliegend rock salt was accommodated by pressure solution and hydrofracturing as is indicated by the microfabrics and bromide concentration in halite. Fractures in siliciclastic domains were filled with fibrous halite and deformed by subgrain rotation recrystallization (SGR). Fluid-rich Zechstein rock salt, on the other hand, was deformed by formation of subgrains and grain boundary migration (GBM). The distribution of mineral phases and fluids had a significant impact on the fabric evolution and on strain localization. U–Pb dating of carbonate phases of the Rotliegend sequence yielded Permian depositional ages and Jurassic to Cretaceous deformation ages, the latter related to diapiric ascent. The combination of results traces a dynamic evolution of the rock fabric inside the diapir structure driven by locally active deformation processes that can be correlated with early stages of halite deposition and diagenesis and syntectonic fabric reorganization related to diapirism in an extensional setting.

Granitoids of the Khara Complex in the Eastern Bureya Continental Massif of the Central Asian Orogenic Belt: Age and Geodynamic Settings of Formation

Petrogeochemical and geochronological studies were conducted on the main varieties of granitoids from the Sedelga and Berezovsky batholiths in the eastern Bureya continental Massif of the Central Asian Fold Belt. SHRIMP-II U–Pb zircon dating of granitoids constrained the time of formation to the Late Triassic and Early Jurassic between 208 to 193 Ma. Our studies revealed that the batholiths contain peraluminous granitoids of high K calc-alkaline series clearly enriched in Rb, K, Th, U, and Pb and relatively depleted in Nb, Ta, Y, HREE, Sr, Ba, Ti, and P and belong to I-type granites. Synthesis of the U–Pb geochronological data indicates that magmatism of tectonic blocks constituting the Bureya Massif at 213–183 Ma was accompanied by the emplacement of granitoid rocks of the Khara Complex. Strike-slip tectonics as a reflection of oblique collision of continental blocks and transform movement might have been an important mechanism of granitoid derivation in the Bureya Massif. Strike-slip faults were likely to penetrate deep into the crust initiating the ascent of an asthenospheric diapir and its fluids into the higher crustal levels thus transferring the heat and triggering melting with appreciable concentrations of Th, U, Rb, Pb and losses of Sr, P, Ba, and Ti. The depletion in Nb, Ta, Ti, and P and the enrichment in Pb make these granitoids similar to those of subduction zones. Granitoids of the Khara Complex display subduction-related (associated with the preceding subduction) and within-plate (associated with an asthenospheric diapir and its fluids) geochemical features, which is characteristic of transform margin environments. The most likely settings for granitoid magmatism to occur throughout the eastern margin of the Bureya Massif in the Late Triassic and Early Jurassic were transform settings associated with asthenospheric diapirism following the cessation of subduction.

Melt Segregation and Depletion During Ascent of Buoyant Diapirs in Subduction Zones

AbstractCold, low‐density diapirs arising from hydrated mantle and/or subducted sediments on the top of subducting slabs have been invoked to transport key chemical signatures to the source region of arc magmas. However, to date there have been few quantitative models to constrain melting in such diapirs. Here we use a two‐phase Darcy‐Stokes‐energy model to investigate thermal evolution, melting, and depletion in a buoyant sediment diapir ascending through the mantle wedge. Using a simplified 2‐D circular geometry, we investigate diapir evolution in three scenarios with increasing complexity. In the first two scenarios we consider instantaneous heating of a diapir by thermal diffusion with and without the effect of the latent heat of melting. Then, these simplified calculations are compared to numerical simulations that include melting, melt segregation, and the influence of depletion on the sediment solidus along pressure‐temperature‐time (P‐T‐t) paths appropriate for ascent through the mantle wedge. The high boundary temperature induces a rim of high porosity, into which new melts are focused and then migrate upward. The rim thus acts like an annulus melt channel, while the effect of depletion buffers additional melt production. Solid matrix flow combined with recrystallization of melt pooled near the top of the diapir can result in large gradients in depletion across the diapir. These large depletion gradients can either be preserved if the diapir leaks melt during ascent, or rehomogenized in a sealed diapir. Overall our numerical simulations predict less melt production than the simplified thermal diffusion calculations. Specifically, we show that diapirs whose ascent paths favor melting beneath the volcanic arc will undergo no more than ~40–50% total melting.

Tourmaline records hydrothermal events related to Zn-Pb mineralization around the Murguía diapi.

The chemical composition of tourmaline has been used as a host environment register as well as a potential exploration tool for mineral deposits. In this study, the textural and chemical composition of tourmalines associated with Zn-Pb mineralizations around the Murguia diapir (Basque Cantabrian Basin, N Spain) are examined to verify if they record the mineralizing events in the area. Petrographically, tourmalines have been differentiated between inherited and authigenic. Colorless authigenic tourmalines are present as halos partially around green and pleochroic detrital grains or as individual crystals. Inherited and authigenic tourmalines are also chemically distinct. Authigenic tourmalines show different X-site occupancies, a Mg/(Mg+Fe) ratio above 0.77, and are aluminum rich and plot to the right of the povondraite-oxidravite join, above the schorl-dravite join. Inherited tourmalines plot within the alkaline (Na+K) group field, and have a Mg/(Mg+Fe) ratio below 0.77. These data suggest that authigenic tourmalines grew under reducing conditions, compatible with the hydrothermal event responsible for the ore deposition and caprock formation during the diapir ascent.

Archean diapirism recorded by vertical sheath folds in the core of the Yalgoo Dome, Yilgarn Craton

In Archean granite-greenstone terrains, both vertical and horizontal tectonic processes can generate the typical dome-and-keel structural geometries. One of the classic Archean domes interpreted to result from a horizontal tectonic regime is the Yalgoo Dome in the Yilgarn Craton of Western Australia. In this contribution, we reassess the structural evolution of the Yalgoo Dome by detailing the structures exposed in the migmatitic core. The core contains multi-scale domes-and-basins. These structures were previously interpreted as the result of superposed folding. Our structural mapping shows that the dome-and-basin patterns result from vertical sheath folds developed in a single, progressive deformation event (D1). This first event most likely resulted from diapiric ascent of a buoyant partially molten tonalitic crust and was overprinted by two contractional episodes (D2–D3) both reflecting bulk EW shortening. The structural evolution proposed here points towards a shift from an early history of vertical tectonics dominated by diapirism to a regional regime dominated by horizontal tectonics.

The nature of the heat source of mafic magmatism during the formation of the Vilyui rift based on the ages of dike swarms and results of numerical modeling

Possible mechanisms of rifting and the thermal regime of the lithosphere beneath the rift zone of the Vilyui sedimentary basin are considered based on the available isotopic ages of dike swarms, rates of sedimentation, and results of numerical modeling. Temporal correlations between the intrusion of mafic magma and a sharp increase in the rate of subsidence and sedimentation in the rift basin prove the contribution of both plate-tectonic and magmatic factors to the formation of the Vilyui rift. The results show a relationship between the rapid extension of the lithosphere and the formation of mafic dike swarms in the Yakutsk-Vilyui Large Igneous Province of the Siberian Platform at the Frasnian-Famennian boundary, with a peak at ~374.1 Ma, and at the end of the Late Devonian, with a peak at ~363.4 Ma. There were two pulses of dike formation during rapid subsidence of the basin basement in the period 380-360 Ma, with a sedimentation rate of 100-130 m/Myr, at a background rate of 10-20 m/Myr. Analysis of numerical thermomechanical models revealed that the best-fit model is that combining the mechanisms of intraplate extension (passive rifting) and the ascent of a mantle magmatic diapir (active rifting). A conclusion about the nature of the heat source of trap magmatism has been drawn: The plume-driven regime of the lithosphere can better explain the dynamics of extension during rifting than the decompression melting mechanism.

Temporal correlation between dyke swarms and crustal extension in the middle Palaeozoic Vilyui rift basin, Siberian platform

This paper presents results from new 40Ar/39Ar isotope dating of nine mafic dykes from three large dyke swarms (Vilyui–Markha, Kontai–Dzherba, and Chara–Sinsk) of the Yakutsk–Vilyui large igneous province (LIP), in addition to a reconstruction of the subsidence history of the middle Palaeozoic Vilyui paleorift basin (eastern Siberian platform). All previously published 40Ar/39Ar and U–Pb dates are summarized. Statistical analysis of the dyke ages reveals repeated magmatic events in the study area. Two major pulses of mafic magmatism are identified: one at the Frasnian–Famennian boundary, with a main peak at ca. 374.1Ma, and another in the latest Devonian with a peak at ca. 363.4Ma. The time of maximum intensity of dyke intrusion coincides (within uncertainty) with rapid subsidence in the Vilyui basin. The minimum total volume of middle Palaeozoic magmatism produced in the Yakutsk–Vilyui LIP is 100–215Kkm3, which is much less than earlier estimates. Most of the mafic material within the Yakutsk–Vilyui LIP is related to the Vilyui basin and associated dyke swarms. Backstripping analysis of sedimentation in depressions of the Vilyui basin was carried out. Estimates were obtained for the spatial distribution of the stretching factor of the crust and mantle lithosphere, averaging 1.17 and 1.44, respectively. The amount of extension due to dyke intrusion is estimated to be ~6%. Backstripping analysis of sedimentation in the Vilyui basin was used to assess the effect of both intraplate far-field forces and upwelling magma flows initiated by a mantle plume. A numerical thermomechanical model was developed to investigate the relations between two possible mechanisms by which the Vilyui rift was initiated: intraplate extension (passive rifting) and the ascent of a mantle magmatic diapir (active rifting). A model considering both of these mechanisms shows the contribution of the far-field extension forces and the effect of convective flows around the mantle plume, assuming a temperature of 1500°C for a heat source operating for ≥10Myr. The best-fit model predicts the formation of the Vilyui basin under extensional conditions with a depth-integrated horizontal driving force of 2.0×1013N/m, which corresponds to a lithostatic pressure of 0.73.

Restitic or not? Insights from trace element content and crystal — Structure of spinels in African mantle xenoliths

The lithospheric architecture of Africa consists of several Archean cratons and smaller cratonic fragments, stitched together and flanked by polycyclic fold belts. Here we investigate the structure and chemistry of spinels from lithospheric mantle xenoliths from distinct tectonic settings, i.e. from the Saharan metacraton in Libya (Waw-En-Namus) which could show archaic chemical features, Cameroon (Barombi Koto and Nyos Lakes) where the Sub Continental Lithospheric Mantle was modified during the Pan-African event and fluxed by asthenospheric melts of the Tertiary Cameroon Volcanic Line and Morocco (Tafraoute, Bou-Ibalrhatene maars) in the Middle Atlas where different metasomatic events have been recorded. From a structural point of view it is to notice that the Libyan spinels can be divided into two groups having different oxygen positional parameter (u>0.2632 and u<0.2627, respectively), while those from Cameroon are in between those values as the Moroccan ones already studied by other authors. The intracrystalline closure temperature (Tc) of the here studied spinels is different among the different samples with one Libyan group (LB I) showing Tc in the range 490–640°C and the other 680–950°C (LB II). Cameroon and Morocco spinels show a Tc in the range 630–760°C. About 150 different spinels have been studied for their trace element content and it can be seen that many of them are related to Cr content, while Zn and Co are not and clearly distinguish the occurrences. Differences in the trace element chemistry, in the structural parameters and in the intracrystalline closure temperatures suggest that a different history should be considered for Cameroon, Morocco and LB I and LB II spinels. Even if it was not considered for this purpose, we tentatively used the Fe2+/Fe3+ vs. TiO2 diagram that discriminate between peridotitic and the so-called “magmatic” spinels, i.e. spinel crystallized from melts. LB I and LB II spinels plot in the peridotitic field while Cameroon and Morocco spinels fall in the magmatic one. Consequently, the xenoliths sampled from a probably juvenile SCLM at the edge of the most important lithospheric roots (i.e. Cameroon and Morocco) apparently have spinels possibly fractionated in situ from percolating melts and do not represent a real spinel-peridotite facies. On the contrary mantle xenoliths from Libya exhibit spinels with peridotitic features compatible with a slow ascent of a mantle diapir (plume).

The mechanism of magma ascent through the solid lithosphere and relation between mantle and crustal diapirism: numerical modeling and natural examples

Diapirism can be regarded as the main mechanism of transport through the lithosphere for both felsic and mafic/ultramafic magmas. However, the lack of field observations makes it difficult to identify the key mechanism responsible for the formation of dome-shaped structures. In this study, emplacement of natural diapirs is reconstructed by numerical experiments handling realistic rheological and petrological models for the crust and mantle lithosphere. Three different regimes of diapiric ascent were established depending on the chosen model rheology: (1) single-stage diapir ascent; (2) pulsating ascent of successive batches of mantle-derived magma to the base of the crust with a periodicity of 2-3 Myr; (3) emplacement of extensive magma bodies in the form of sills either beneath the base of the crust (underplating) or to deeper mantle levels. The timescale of 30 Myr for a heat source at the base of the lithosphere is sufficient to initiate the ascent of a diapir through the mantle and crust. The study provides the estimates of rheological properties of the lithosphere and partially molten material at which diapiric ascent through the mantle and crust can occur.

The lithospheric shear-wave velocity structure of Saudi Arabia: Young volcanism in an old shield

We investigate the lithospheric shear-wave velocity structure of Saudi Arabia by conducting H-κ stacking analysis and jointly inverting teleseismic P-receiver functions and fundamental-mode Rayleigh wave group velocities at 56 broadband stations deployed by the Saudi Geological Survey (SGS). The study region, the Arabian plate, is traditionally divided into the western Arabian shield and the eastern Arabian platform: The Arabian shield itself is a complicated mélange of crustal material, composed of several Proterozoic terrains separated by ophiolite-bearing suture zones and dotted by outcropping Cenozoic volcanic rocks (locally known as harrats). The Arabian platform is primarily covered by 8 to 10km of Paleozoic, Mesozoic and Cenozoic sedimentary rocks. Our results reveal high Vp/Vs ratios in the region of Harrat Lunayyir, which are interpreted as solidified magma intrusions from old magmatic episodes in the shield. Our results also indicate slow velocities and large upper mantle lid temperatures below the southern and northern tips of the Arabian shield, when compared with the values obtained for the central shield. We argue that our inferred patterns of lid velocity and temperature are due to heating by thermal conduction from the Afar plume (and, possibly, the Jordan plume), and that volcanism in western Arabia may result from small-scale adiabatic ascent of magma diapirs.

Timing of UHP exhumation and rock fabric development in gneiss domes containing the world's youngest eclogite Facies rocks, southeastern Papua New Guinea

AbstractThe youngest known ultrahigh‐pressure (UHP) rocks in the world occur in the Woodlark Rift of southeastern Papua New Guinea. Since their crystallization in the Late Miocene to Early Pliocene, these eclogite facies rocks have been rapidly exhumed from mantle depths to the surface and today they remain in the still‐active geodynamic setting that caused this exhumation. For this reason, the rocks provide an excellent opportunity to study rates and processes of (U)HP exhumation. We present New Rb–Sr results from 12 rock samples from eclogite‐bearing gneiss domes in the D'Entrecasteaux Islands, and use those results to examine the time lag between (U)HP metamorphism and later ductile thinning, penetrative fabric development and accompanying metamorphic retrogression at amphibolite facies conditions during their exhumation. A Rb–Sr age for a sample of mafic eclogite (with no preserved coesite) from the core zone of the Mailolo gneiss dome (Fergusson Island) provides a new estimate of the timing of HP metamorphism (5.6 ± 1.6 Ma). The strongly deformed quartzofeldspathic and granitic gneisses (90–95% by volume) that enclose variably retrogressed relict blocks of mafic eclogite (5–10% by volume) yield Rb–Sr isochron ages from 4.4 to 2.4 Ma. For the UHP‐bearing gneisses of Mailolo dome, previously published U–Pb ages on zircon and our Rb–Sr isochron ages are consistent with a mean time lag of 2.2 ± 1.5 Ma (~95% c.i.) for passage of the rock between eclogite and amphibolite facies conditions. New thermobarometric data indicate that the main syn‐exhumational foliation developed at amphibolite facies conditions of 630–665 °C and 12.1–14.4 kbar. These pressure estimates indicate that the lower crust of the Woodlark Rift was unusually thick (&gt;40 km) at the time of the amphibolite facies overprint, possibly as a result of accumulation and underplating of UHP‐derived material from below. Our data imply a minimum unroofing rate of 10 ± 7 mm year−1 (~95% c.i.) for the (U)HP body from minimum HP depths (73 ± 7 km) to lower crustal depths. This minimum unroofing rate reinforces previous inferences that the exhumation from the mantle to the surface of the gneiss domes in the D'Entrecasteaux Islands took place at plate tectonic rates. On the basis of previous structural studies and the new thermobarometry, we attribute the high (cm year−1) exhumation to diapiric ascent of the partially molten terrane from mantle depths, with a secondary contribution from pure shear thinning of the terrane after its arrival in the crust.

A large mantle water source for the northern San Andreas fault system: a ghost of subduction past

Recent research indicates that the shallow mantle of the Cascadia subduction margin under near-coastal Pacific Northwest, USA is cold and partially serpentinized, storing large quantities of water in this wedge-shaped region. Such a wedge probably formed to the south in California during an earlier period of subduction. We show by numerical modeling that after subduction ceased with the creation of the San Andreas Fault System (SAFS), the mantle wedge warmed, slowly releasing its water over a period of more than 25 Ma by serpentine dehydration into the crust above. This deep, long-term water source could facilitate fault slip in San Andreas System at low shear stresses by raising pore pressures in a broad region above the wedge. Moreover, the location and breadth of the water release from this model gives insights into the position and breadth of the SAFS. Such a mantle source of water also likely plays a role in the occurrence of non-volcanic tremor (NVT) that has been reported along the SAFS in central California. This process of water release from mantle depths could also mobilize mantle serpentinite from the wedge above the dehydration front, permitting upward emplacement of serpentinite bodies by faulting or by diapiric ascent. Specimens of serpentinite collected from tectonically emplaced serpentinite blocks along the SAFS show mineralogical and structural evidence of high fluid pressures during ascent from depth. Serpentinite dehydration may also lead to tectonic mobility along other plate boundaries that succeed subduction, such as other continental transforms, collision zones, or along present-day subduction zones where spreading centers are subducting.

Intrusion of lamprophyre dyke and related deformation effects in the host rock salt: A case study from the Loulé diapir, Portugal

A rock salt–lamprophyre dyke contact zone (sub-vertical, NE–SW strike) was investigated for its petrographic, mechanic and physical properties by means of anisotropy of magnetic susceptibility (AMS) and rock magnetic properties, coupled with quantitative microstructural analysis and thermal mathematical modelling. The quantitative microstructural analysis of halite texture and solid inclusions revealed good spatial correlation with AMS and halite fabrics. The fabrics of both lamprophyre and rock salt record the magmatic intrusion, “plastic” flow and regional deformation (characterized by a NW–SE trending steep foliation). AMS and microstructural analysis revealed two deformation fabrics in the rock salt: (1) the deformation fabrics in rock salt on the NW side of the dyke are associated with high temperature and high fluid activity attributed to the dyke emplacement; (2) On the opposite side of the dyke, the emplacement-related fabric is reworked by localized tectonic deformation. The paleomagnetic results suggest significant rotation of the whole dyke, probably during the diapir ascent and/or the regional Tertiary to Quaternary deformation.

Zircon U–Pb ages of the Mianning–Dechang syenites, Sichuan Province, southwestern China: Constraints on the giant REE mineralization belt and its regional geological setting

The Himalayan Mianning–Dechang (MD) rare earth element (REE) belt in western Sichuan Province, southwestern China, is approximately 270km long and 15km wide, and contains total reserves of more than 3Mt of light REEs (LREEs), comprising one giant (Maoniuping), one large (Dalucao), two small–medium-sized (Muluozhai and Lizhuang), and numerous smaller REE deposits. The belt occurs within the eastern Indo-Asian collision zone (EIACZ), where its location is controlled by large-scale strike-slip faults and tensional fissure zones. Himalayan carbonatite–syenite complexes consist predominantly of alkaline syenite stocks and carbonatite sills or dikes that host REE mineralization. Previous studies have reported inconsistent ages for alkaline magmatism syenite formation and REE mineralization. Here, we present new results of sensitive high-resolution ion micro-probe U–Pb dating of zircons from syenites from the Dalucao, Maoniuping, Lizhuang and Diaoloushan areas, the first systematic and precise age determinations for these rocks in the MD belt. The new data give concordant ages of 12.13±0.19 and 11.32±0.23Ma for the Dalucao deposit, 22.81±0.31 and 21.3±0.4Ma for Maoniuping, 26.77±0.32Ma for Muluozhai, and 27.41±0.35Ma for Lizhuang. These ages, which should be regarded as maximum ages for the REE mineralization in the study area, can be split into two groups, i.e. 11–12Ma in the southern part of the MD belt and 12–27Ma in the northern part, suggesting a progression of magmatism from north to south. These data suggest that the majority of carbonatite–syenite magmatism within the EIACZ occurred during the main stage of Himalayan metallogenesis. The ages presented in this study suggest that strike-slip shear along the MD belt was initiated at ca. 27Ma and ended ca. 12Ma. This timing is consistent with movements along the adjacent Ailaoshan–Red River strike-slip fault in southeastern Tibet (to the south of the MD belt) and one of the three Cenozoic strike-slip faults in eastern Tibet. Ascent of an asthenospheric mantle diapir beneath the EIACZ in the Cenozoic may have provided a thermal mechanism for the generation of magmas that formed the carbonatite–syenite complexes in the study area. Alternatitvely, the magmas may have been generated by decompression melting associated with the transition from a transpressional to a transtensional regime at 38–40Ma. The precise age results for syenite magmatism in the study area indicate that this transition occurred prior to carbonatite–syenite magmatism and the formation of the MD REE belt, which is consistent with the regional tectonic model.

Origin of Basalts in a Hot Subduction Setting: Petrological and Geochemical Insights from Mt. Baker, Northern Cascade Arc

The northern Cascade arc is an end-member ‘hot’ subduction zone where slab dehydration may be essentially complete prior to reaching sub-arc depths, presenting a potential problem for a flux melting origin for arc basalts. Nevertheless, mafic lavas from the Mt. Baker volcanic field, the most magmatically productive volcano in the northern Cascade arc, record subduction input and compositions typical of calc-alkaline magmas. Relative to normal mid-ocean ridge basalt, the most primitive lavas at Mt. Baker show elevated abundances of large ion lithophile elements (LILE), Pb and light rare earth elements (LREE). High field strength elements (HFSE) and heavy REE (HREE) are depleted relative to LILE. Reconstructed primary magmas define three groups: Group I calcalkaline basalts (Sulphur Creek, Lake Shannon, Hogback) have lower SiO2, (La/Yb)N and LILE/HFSE, and are olivineand plagioclase-phyric; Group II high-Mg basaltic andesites (Tarn Plateau, Cathedral Crag) also contain clinopyroxene phenocrysts and have higher H2O, SiO2, (La/Yb)N and LILE/HFSE; Group III low-K olivine tholeiite (Park Butte) has the lowest (La/Yb)N but highest LILE/HFSE.The redox states of all groups lie between QFM (quartz^fayalite^magnetite) and QFM þ1·4. Pb, Sr and Nd isotope compositions are similar among all the analysed samples and are consistent with a depleted mantle source modified by input from the subducting slab. Trace element^isotope modeling indicates a subduction component composed predominantly of metabasaltderived fluid with lesser amounts of sediment melt and metabasalt melt. Group I basalts record the smallest melt fractions ( 5^7%), lowest water contents (1·5^2·1wt %), and highest temperatures and pressures of mantle segregation (up to 13548C, 1·5 GPa) from lherzolitic ( spinel) residues. Group II basaltic andesites show the greatest extents of mantle metasomatism, the highest water contents (2·7^3·7 wt %) and partial melt fractions (10^12%), and segregated from harzburgite at 12708C, 1GPa, consistent with pooling of melts at the Moho. Group III records P^Tconditions similar to Group I ( 1·4 GPa, 13268C) but melt fractions (12%) and mantle residues (harzburgite) are more similar to Group II, and H2O contents ( 2·1wt %) are intermediate. Melting beneath Mt. Baker was initiated by dehydration melting of amphibole peridotite at 95 km and 10208C, within the stability field of garnet lherzolite. Initial melt fractions were small ( 1^2%) and near water-saturation. Phase equilibria and trace element modeling show no evidence for garnet-bearing mantle residues, indicating that progressive melting during ascent of diapirs through the hot core of the convecting mantle wedge reduced H2O contents and erased any residual garnet signature. Because fluid release from the slab is restricted to the forearc, mantle hydrated at shallow depths in the serpentine and/or chlorite stability fields must be down-dragged to the region of amphibole stability to initiate dehydration melting.

A quasi-asymptotic model for the ascent of magmatic diapirs

Recent asymptotic analysis for the slow, steady ascent of a zero-traction sphere in a power-law fluid with temperature-dependent viscosity is applied to the derivation of a model for magmatic diapirism in the Earth’s mantle and continental crust. As an example, the model is applied to the case of a mantle diapir rising through country rock that is treated as a strongly shear-thinning fluid. It is found that, within the geologically relevant timescale of years, the diapir would only encounter two of the four possible asymptotic regimes that the earlier analysis uncovered, before stalling at a considerably greater depth below the Earth’s surface than was predicted by an earlier model. The implications of this result for the modelling of magmatic diapirism are discussed.