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

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