Search

During continental collision, considerable amounts of buoyant continental crust subduct to depth and subsequently exhume. Whether various exhumation paths contribute to contrasting styles of magmatism across modern collision zones is unclear. Here we present 2D thermomechanical models of continental collision combined with petrological databases to investigate the effect of the main contrasting buoyancy forces, in the form of continental crustal buoyancy versus oceanic slab age (i.e., its thickness). We specifically focus on the consequences for crustal exhumation mechanisms and magmatism. Results indicate that it is mainly crustal density that determines the degree of steepening of the subducting continent and separates the models' parameter space into two regimes. In the first regime, high buoyancy values (∆ρ > 500 kg/m3) steepen the slab most rapidly (to 45–58°), leading to opening of a gap in the subduction channel through which the subducted crust exhumes (“subduction channel crustal exhumation”). A shift to a second regime (“underplating”) occurs when the density contrast is reduced by 50 kg/m3. In this scenario, the slab steepens less (to 37–50°), forcing subducted crust to be placed below the overriding plate. Importantly, the magmatism changes in the two cases: Crustal exhumation through the subduction channel is mainly accompanied by a narrow band of mantle melts, while underplating leads to widespread melting of mixed sources. Finally, we suggest that the amount (or density) of subducted continental crust, and the resulting buoyancy forces, could contribute to contrasting collision styles and magmatism in the Alps and Himalayas/Tibet., Key Points We analyse buoyancy forces in 2D thermomechanical‐petrological models of continental collision and resulting postcollisional magmatismCrustal buoyancy controls whether subducting crust exhumes between plates with only mantle melts or underplates with mixed meltingTiming and distribution of crust and melting for two endmembers fits Alps and Himalaya/Tibet

\ud Crystallization temperatures of primitive olivine crystals have been widely used as both a proxy for, or an intermediate step in calculating, mantle temperatures. The olivine-spinel aluminum-exchange thermometer has been applied to samples from mid-ocean ridges and large igneous provinces, yielding considerable variability in olivine crystallization temperatures. We supplement the existing data with new crystallization temperature estimates for Hawaii, between 1282 ± 21 and 1375 ± 19°C. Magmatic temperatures may be linked to mantle temperatures if the thermal changes during melting can be quantified. The magnitude of this temperature change depends on melt fraction, itself controlled by mantle temperature, mantle composition and lithosphere thickness. Both mantle composition and lithosphere thickness vary spatially and temporally, with systematic differences between mid-ocean ridges, ocean islands and large igneous provinces. For crystallization temperatures to provide robust evidence of mantle temperature variability, the controls of lithosphere thickness and mantle lithology on crystallization temperature must be isolated. We develop a multi-lithology melting model for predicting crystallization temperatures of magmas in both intra-plate volcanic provinces and mid-ocean ridges. We find that the high crystallization temperatures seen at mantle plume localities do require high mantle temperatures. In the absence of further constraints on mantle lithology or melt productivity, we cannot robustly infer variable plume temperatures between ocean-islands and large igneous provinces from crystallization temperatures alone; for example, the extremely high crystallization temperatures obtained for the Tortugal Phanerozoic komatiite could derive from mantle of comparable temperature to modern-day Hawaii. This work demonstrates the limit of petrological thermometers when other geodynamic parameters are poorly known.\ud \ud Plain Language Summary\ud The temperature inside the Earth varies a lot. There are many ways of measuring the mantle's temperature in the present-day, but to understand how our planet has changed through time, we need to know how hot its interior was in the past. One of the ways we can estimate mantle temperatures from ancient and modern rocks is from their crystal chemistry. By measuring the aluminum content in crystals of olivine, we can estimate their crystallization temperature. We do this for crystals from Hawaii. To turn the crystallization temperatures into the mantle temperature we need to know the proportions of minerals the mantle is made of. However, we often don't know all of this information. Using a new model of mantle melting we can calculate how uncertain the mantle temperature is when we only have a crystallization temperature. We find that the mantle under Hawaii is 1582 ± 65°C, much hotter than normal mantle, which has a temperature 1364 ± 23°C. We also apply our method to crystallization temperatures from other locations, including ancient volcanic rocks. We find that crystallization temperatures from large igneous provinces, formed by unusually hot mantle, are consistent with their mantle having a similar temperature to Hawaii

Carbonate veins in the igneous basement and in the lithified sedimentary cover of the Cocos Ridge at International Ocean Discovery Program (IODP) Hole 344‐U1414A reveal the hydrologic system and fluid‐rock interactions. IODP Hole 344‐U1414A was drilled on the northern flank of the Cocos Ridge and is situated 1 km seaward from the Middle America Trench offshore Costa Rica. Isotopic and elemental compositions were analyzed to constrain the fluid source of the carbonate veins and to reveal the thermal history of Hole 344‐U1414A. The formation temperatures (oxygen isotope thermometer) of the carbonate veins in the lithified sedimentary rocks range from 70 to 92 °C and in the basalt from 32 to 82 °C. 87Sr/86Sr ratios of the veins in the altered Cocos Ridge basalt range between 0.707307 and 0.708729. The higher ratios are similar to seawater strontium ratios in the Neogene. 87Sr/86Sr ratios lower 0.7084 indicate exchange of Sr with the igneous host rock. The calcite veins hosted by the sedimentary rocks are showing more primitive 87Sr/86Sr ratios, Key Points Oxygen and strontium isotope ratios of carbonate veins suggest seawater and hydrothermally modified seawater as fluid sourceIsotopic and elemental compositions reveal fluid‐rock interactions in the Cocos Ridge basalt and in the lithified sedimentary rocksIntraplate seamount volcanism in the area between the Galapagos hot spot and the Cocos Island acted as an additional heating source

The architecture and mechanical properties of the subduction interface impact large‐scale subduction processes, including mass and volatile recycling, upper‐plate orogenesis, and seismic behavior. The nature of the deep subduction interface, where a dominantly frictional megathrust likely transitions to a distributed ductile shear zone, is poorly understood, due to a lack of constraints on rock types, strain distribution, and interface thickness in this depth range. We characterized these factors in the Condrey Mountain Schist, a Late Jurassic to Early Cretaceous subduction complex in northern California that consists of an upper and lower unit. The Lower Condrey unit is predominantly pelagic and hemipelagic metasediment with m‐to km‐scale metamafic and metaserpentinitic ultramafic lenses all deformed at epidote blueschist facies (0.7–1.1 GPa, 450°C). Major and trace element geochemistry suggest tectonic erosion of the overriding plate sourced all ultramafic and some mafic lenses. We identified two major ductile thrust zones responsible for Lower Condrey unit assembly, with earlier strain distributed across the structural thickness between the ductile thrusts. The Lower Condrey unit records distributed deformation across a sediment‐dominated, 2+ km thick shear zone, possibly consistent with low velocity zones observed in modern subduction zones, despite subducting along a sediment poor, tectonically erosive margin. Periodic strain localization occurred when rheological heterogeneities (i.e., km‐scale ultramafic lenses) entered the interface, facilitating underplating that preserved 10%–60% of the incoming sediment. Modern mass and volatile budgets do not account for erosive margin underplating, so improved quantification is crucial for predicting mass and volatile net flux to Earth′s interior., Geochemistry, Geophysics, Geosystems, 22 (3), ISSN:1525-2027

International Ocean Discovery Program Expeditions 352 and 351 drilled into the Western Pacific Izu‐Bonin forearc and rear arc. The drill cores revealed that the forearc is composed of forearc basalts (FAB) and boninites and the rear arc consists of FAB‐like rocks. These rocks are pervaded by calcite veins. Blocky vein microtextures enclosing host rock fragments dominate in all locations and suggest hydrofracturing and advective fluid flow. Significant diffusion‐fed and crystallization pressure‐driven antitaxial veining is restricted to the rear arc. The lack of faults and presence of an Eocene sedimentary cover in the rear arc facilitated antitaxial veining. Rare earth element and isotopic (δ18O, δ13C, 87Sr/86Sr, and Δ47) tracers indicate varying parental fluid compositions ranging from pristine to variably modified seawater. The most pristine seawater signatures are recorded by FAB‐hosted low‐T (95% seawater and, Geochemistry, Geophysics, Geosystems, 21 (2), ISSN:1525-2027

In seismic tomography we map the wave speed structure inside the Earth, but we ultimately seek to interpret those images in terms of physical parameters. This is challenging because many parameters can trade‐off with each other to produce a given wave speed. The problem is compounded by the convention of mapping seismic structures as perturbations relative to a 1‐D reference model, rather than absolute wave speeds. Using a full waveform tomography model of Europe as a case study, we quantify the extent to which thermochemical and dynamic properties can be constrained using only S wave speed, expressed in absolute values. The wave speed distributions of this tomography model are compared with 4 million thermochemical models, whose seismic properties are computed via thermodynamic modeling. These models sample the full range of realistic mantle compositions, including variable water and melt contents, and mineral intrinsic anelasticity is taken into account. Intrinsic anelasticity causes waves to travel more slowly at higher temperatures, leading to seismic attenuation, but the sensitivity of the wave speed reduction to temperature is, in turn, controlled by the wave frequency. Global studies of surface waves indicate an anticorrelation between S wave speed and attenuation. We therefore only retain thermochemical models satisfying this anticorrelation. Our study indicates that the frequency dependence of anelasticity, α, depends on temperature or rheology, with α ≈ 0.1 being most appropriate in cold or lithospheric mantle and α ≈ 0.3 in warmer regions (i.e., the asthenosphere). Additionally, the slowest regions require specific compositions and/or a velocity‐weakening mechanism, such as partial melting, elastically accommodated grain boundary sliding, or water., Key Points Knowing the S wave speed precisely with a loose constraint on mantle attenuation places strong constraints on the temperatureThe frequency dependence of anelasticity α is variable: ~0.1 in cold/lithospheric mantle and ~0.3 in warm/asthenospheric mantleEither partial melting or elastically accommodated grain boundary sliding is required to explain the slowest wave speeds at 130 km

Viscosity has been determined during isothermal crystallization of an andesite from Tungurahua volcano (Ecuador). Viscosity was continuously recorded using the concentric cylinder method and employing a Pt‐sheathed alumina spindle at 1 bar and from 1400°C to subliquidus temperatures to track rheological changes during crystallization. The disposable spindle was not extracted from the sample but rather left in the sample during quenching thus preserving an undisturbed textural configuration of the crystals. The inspection of products quenched during the crystallization process reveals evidence for heterogeneous crystal nucleation at the spindle and near the crucible wall, as well as crystal alignment in the flow field. At the end of the crystallization, defined when viscosity is constant, plagioclase is homogeneously distributed throughout the crucible (with the single exception of experiment performed at the lowest temperature). In this experiments, the crystallization kinetics appear to be strongly affected by the stirring conditions of the viscosity determinations. A TTT (Time‐Temperature‐Transformation) diagram illustrating the crystallization “nose” for this andesite under stirring conditions and at ambient pressure has been constructed. We further note that at a given crystal content and distribution, the high aspect ratio of the acicular plagioclase yields a shear‐thinning rheology at crystal contents as low as 13 vol %, and that the relative viscosity is higher than predicted from existing viscosity models. These viscosity experiments hold the potential for delivering insights into the relative influences of the cooling path, undercooling, and deformation on crystallization kinetics and resultant crystal morphologies, as well as their impact on magmatic viscosity., Key Points: Viscosity determination at subliquidus temperatures of an andesitic sampleUse of disposable Pt‐sheathed alumina spindlesDescription of viscosity and crystallization kinetics under constant stirring