Glacial isostatic uplift of the European Alps

Following the last glacial maximum (LGM), the demise of continental ice sheets induced crustal rebound in tectonically stable regions of North America and Scandinavia that is still ongoing. Unlike the ice sheets, the Alpine ice cap developed in an orogen where the measured uplift is potentially attributed to tectonic shortening, lithospheric delamination and unloading due to deglaciation and erosion. Here we show that ∼90% of the geodetically measured rock uplift in the Alps can be explained by the Earth's viscoelastic response to LGM deglaciation. We modelled rock uplift by reconstructing the Alpine ice cap, while accounting for postglacial erosion, sediment deposition and spatial variations in lithospheric rigidity. Clusters of excessive uplift in the Rhône Valley and in the Eastern Alps delineate regions potentially affected by mantle processes, crustal heterogeneity and active tectonics. Our study shows that even small LGM ice caps can dominate present-day rock uplift in tectonically active regions

Linked 3-D modelling of megathrust earthquake-tsunami events: from subduction to tsunami run up

How does megathrust earthquake rupture govern tsunami behaviour? Recent modelling advances permit evaluation of the influence of 3-D earthquake dynamics on tsunami genesis, propagation, and coastal inundation. Here, we present and explore a virtual laboratory in which the tsunami source arises from 3-D coseismic seafloor displacements generated by a dynamic earthquake rupture model. This is achieved by linking open-source earthquake and tsunami computational models that follow discontinuous Galerkin schemes and are facilitated by highly optimized parallel algorithms and software. We present three scenarios demonstrating the flexibility and capabilities of linked modelling. In the first two scenarios, we use a dynamic earthquake source including time-dependent spontaneous failure along a 3-D planar fault surrounded by homogeneous rock and depth-dependent, near-lithostatic stresses. We investigate how slip to the trench influences tsunami behaviour by simulating one blind and one surface-breaching rupture. The blind rupture scenario exhibits distinct earthquake characteristics (lower slip, shorter rupture duration, lower stress drop, lower rupture speed), but the tsunami is similar to that from the surface-breaching rupture in run-up and length of impacted coastline. The higher tsunami-generating efficiency of the blind rupture may explain how there are differences in earthquake characteristics between the scenarios, but similarities in tsunami inundation patterns. However, the lower seafloor displacements in the blind rupture result in a smaller displaced volume of water leading to a narrower inundation corridor inland from the coast and a 15 per cent smaller inundation area overall. In the third scenario, the 3-D earthquake model is initialized using a seismo-thermo-mechanical geodynamic model simulating both subduction dynamics and seismic cycles. This ensures that the curved fault geometry, heterogeneous stresses and strength and material structure are consistent with each other and with millions of years of modelled deformation in the subduction channel. These conditions lead to a realistic rupture in terms of velocity and stress drop that is blind, but efficiently generates a tsunami. In all scenarios, comparison with the tsunamis sourced by the time-dependent seafloor displacements, using only the time-independent displacements alters tsunami temporal behaviour, resulting in later tsunami arrival at the coast, but faster coastal inundation. In the scenarios with the surface-breaching and subduction-initialized earthquakes, using the time-independent displacements also overpredicts run-up. In the future, the here presented scenarios may be useful for comparison of alternative dynamic earthquake-tsunami modelling approaches or linking choices, and can be readily developed into more complex applications to study how earthquake source dynamics influence tsunami genesis, propagation and inundation

Plio-Quaternary tectonic evolution of the southern margin of the Alboran Basin (Western Mediterranean)

We thank the members of the SARAS and Marlboro cruises in 2011 and 2012. We thank Emanuele Lodolo, Jacques Déverchère, Guillermo Booth-Rea for their helpful comments and discussion. We also thank the editor, Federico Rossetti, for the attention provided to this article. This work was funded by the French program Actions Marges, the EUROFLEETS program (FP7/2007-2013; no. 228344) and project FICTS-2011-03-01. The French program ANR-17-CE03-0004 also supported this work. Seismic reflection data were processed using the Seismic UNIX SU and Geovecteur software. The processed seismic data were interpreted using Kingdom IHS Suite© software. This work also benefited from the Fauces Project (Ref CTM2015-65461-C2-R; MINCIU/FEDER) financed by Ministerio de Economía y Competitividad y al Fondo Europeo de Desarrollo Regiona (FEDER).Progress in the understanding and dating of the sedimentary record of the Alboran Basin allows us to propose a model of its tectonic evolution since the Pliocene. After a period of extension, the Alboran Basin underwent a progressive tectonic inversion starting around 9–7.5 Ma. The Alboran Ridge is a NE–SW transpressive structure accommodating the shortening in the basin. We mapped its southwestern termination, a Pliocene rhombic structure exhibiting series of folds and thrusts. The active Al-Idrissi Fault zone (AIF) is a Pleistocene strike-slip structure trending NNE– SSW. The AIF crosses the Alboran Ridge and connects to the transtensive Nekor Basin and the Nekor Fault to the south. In the Moroccan shelf and at the edge of a submerged volcano we dated the inception of the local subsidence at 1.81– 1.12 Ma. The subsidence marks the propagation of the AIF toward the Nekor Basin. Pliocene thrusts and folds and Quaternary transtension appear at first sight to act at different tectonic periods but reflect the long-term evolution of a transpressive system. Despite the constant direction of Africa– Eurasia convergence since 6 Ma, along the southern margin of the Alboran Basin, the Pliocene–Quaternary compression evolves from transpressive to transtensive along the AIF and the Nekor Basin. This system reflects the logical evolution of the deformation of the Alboran Basin under the indentation of the African lithosphere.This research has been supported by the CNRSINSU-TOTAL-BRGM-IFREMER Actions Marges program, EUROFLEETS program FP7/2007-2013 (grant no. 228344), EU Regional Structural Fund (grant no. FICTS-2011-03-01) and DAMAGE Project (project no. FEDER/CGL2016-80687-R AEI), Fauces Project (project no. FEDER/CTM2015-65461-C2-R; MINCI), ALBAMAR Project (project no. ANR/ANR-17-CE03-0004)

Plume‐Induced Sinking of Intracontinental Lithospheric Mantle: An Overlooked Mechanism of Subduction Initiation?

Although many different mechanisms for subduction initiation have been proposed, only few of them are viable in terms of consistency with observations and reproducibility in numerical experiments. In particular, it has recently been demonstrated that intra‐oceanic subduction triggered by an upwelling mantle plume could greatly contribute to the onset and operation of plate tectonics in the early and, to a lesser degree, modern Earth. On the contrary, the initiation of intra‐continental subduction still remains underappreciated. Here we provide an overview of 1) observational evidence for upwelling of hot mantle material flanked by downgoing proto‐slabs of sinking continental mantle lithosphere, and 2) previously published and new numerical models of plume‐induced subduction initiation. Numerical modeling shows that under the condition of a sufficiently thick (>100 km) continental plate, incipient downthrusting at the level of the lowermost lithospheric mantle can be triggered by plume anomalies of moderate temperatures and without significant strain‐ and/or melt‐related weakening of overlying rocks. This finding is in contrast with the requirements for plume‐induced subduction initiation within oceanic or thinner continental lithosphere. As a result, plume‐lithosphere interactions within continental interiors of Paleozoic‐Proterozoic‐(Archean) platforms are the least demanding (and thus potentially very common) mechanism for initiation of subduction‐like foundering in the Phanerozoic Earth. Our findings are supported by a growing body of new geophysical data collected in various intra‐continental areas. A better understanding of the role of intra‐continental mantle downthrusting and foundering in global plate tectonics and, particularly, in the initiation of “classic” ocean‐continent subduction will benefit from more detailed follow‐up investigations

Low-Temperature Plasticity in Olivine: Grain Size, Strain Hardening, and the Strength of the Lithosphere

Plastic deformation of olivine at relatively low temperatures (i.e., low-temperature plasticity) likely controls the strength of the lithospheric mantle in a variety of geodynamic contexts. Unfortunately, laboratory estimates of the strength of olivine deforming by low-temperature plasticity vary considerably from study to study, limiting confidence in extrapolation to geological conditions. Here we present the results of deformation experiments on olivine single crystals and aggregates conducted in a deformation-DIA at confining pressures of 5 to 9 GPa and temperatures of 298 to 1473 K. These results demonstrate that, under conditions in which low-temperature plasticity is the dominant deformation mechanism, fine-grained samples are stronger at yield than coarse-grained samples, and the yield stress decreases with increasing temperature. All samples exhibited significant strain hardening until an approximately constant flow stress was reached. The magnitude of the increase in stress from the yield stress to the flow stress was independent of grain size and temperature. Cyclical loading experiments revealed a Bauschinger effect, wherein the initial yield strength is higher than the yield strength during subsequent cycles. Both strain hardening and the Bauschinger effect are interpreted to result from the development of back stresses associated with long-range dislocation interactions. We calibrated a constitutive model based on these observations, and extrapolation of the model to geological conditions predicts that the strength of the lithosphere at yield is low compared to previous experimental predictions but increases significantly with increasing strain. Our results resolve apparent discrepancies in recent observational estimates of the strength of the oceanic lithosphere.Support for this research was provided by Natural Environment Research Council (NERC) grant NE/M000966/1 and NSF Division of Earth Sciences grants 1255620, 1464714, and 1550112. D.E.J.A. acknowledges funding from the Royal Academy of Engineering through a research fellowship

The impact of water on slip system activity in olivine and the formation of bimodal crystallographic preferred orientations

Crystallographic preferred orientations (CPOs) in olivine are widely used to infer the mechanisms, conditions, and kinematics of deformation of mantle rocks. Recent experiments on water-saturated olivine were the first to produce a complex CPO characterised by bimodal orientation distributions of both [100] and [001] axes and inferred to form by combined activity of (001)[100], (100)[001], and (010)[100] slip. This result potentially provides a new microstructural indicator of deformation in the presence of elevated concentrations of intracrystalline hydrous point defects and has implications for the interpretation of seismic anisotropy. Here, we document a previously unexplained natural example of this CPO type in a xenolith from Lesotho and demonstrate that it too may be explained by elevated concentrations of hydrous point defects. We test and confirm the hypothesis that combined (001)[100], (100)[001], and (010)[100] slip were responsible for formation of this CPO by (1) using high-angular resolution electron backscatter diffraction to precisely characterise the dislocation types present in both the experimental and natural samples and (2) employing visco-plastic self-consistent simulations of CPO evolution to assess the ability of these slip systems to generate the observed CPO. Finally, we utilise calculations based on effective-medium theory to predict the anisotropy of seismic wave velocities arising from the CPO of the xenolith. Maxima in S-wave velocities and anisotropy are parallel to both the shear direction and shear plane normal, whereas maxima in P-wave velocities are oblique to both, adding complexity to interpretation of deformation kinematics from seismic anisotropy.D. Wallis, L.N. Hansen, and A.J. Wilkinson acknowledge support from the Natural Environment Research Council Grant NE/M000966/1. M. Tasaka acknowledges support through a JSPS Research Fellowship for Young Scientists (26-4879) and the Japan Society for the Promotion of Science (16K17832). D.L. Kohlstedt acknowledges support through NASA Grant NNX15AL53G. K.M. Kumamoto acknowledges support through NSF Division of Earth Science grants 1255620 and 1625032

The role of water and compression in the genesis of alkaline basalts: Inferences from the Carpathian-Pannonian region

We present a new model for the formation of Plio-Pleistocene alkaline basalts in the central part of the Carpathian-Pannonian region (CPR). Based on the structural hydroxyl content of clinopyroxene megacrysts, the ‘water’ content of their host basalts is 2.0–2.5 wt.%, typical for island arc basalts. Likewise, the source region of the host basalts is ‘water’ rich (290–660 ppm), akin to the source of ocean island basalts. This high ‘water’ content could be the result of several subduction events from the Mesozoic onwards (e.g. Penninic, Vardar and Magura oceans), which have transported significant amounts of water back to the upper mantle, or hydrous plumes originating from the subduction graveyard beneath the Pannonian Basin. The asthenosphere with such a relatively high ‘water’ content beneath the CPR may have been above the ‘pargasite dehydration’ (90 km) solidi. This means that neither decompressional melting nor the presence of voluminous pyroxenite and eclogite lithologies are required to explain partial melting. While basaltic partial melts have been present in the asthenosphere for a long time, they were not extracted during the syn-rift phase, but were only emplaced at the onset of the subsequent tectonic inversion stage at ~8–5 Ma. We propose that the extraction has been facilitated by evolving vertical foliation in the asthenosphere as a response to the compression between the Adriatic indenter and the stable European platform. The vertical foliation and the prevailing compression effectively squeezed the partial basaltic melts from the asthenosphere. The overlying lithosphere may have been affected by buckling in response to compression, which was probably accompanied by formation of deep faults and deformation zones. These zones formed conduits towards the surface for melts squeezed out of the asthenosphere. This implies that basaltic partial melts could be present in the asthenosphere in cases where the bulk ‘water’ content is relatively high (>~200 ppm) at temperatures exceeding ~1000–1100 °C. These melts could be extracted even under a compressional tectonic regime, where the combination of vertical foliation in the asthenosphere and deep fractures and deformation zones in the folded lithosphere provides pathways towards the surface. This model is also valid for deep seated transpressional or transtensional fault zones in the lithosphere

Dislocation interactions during low-temperature plasticity of olivine and their impact on the evolution of lithospheric strength

The strength of the lithosphere is typically modelled based on constitutive equations for steady-state flow. However, strain hardening may cause significant evolution of strength in the colder load-bearing portion of the lithosphere. Recent rheological data from low-temperature deformation experiments on olivine suggest that strain hardening occurs due to the presence of temperature-independent back stresses generated by long-range elastic interactions among dislocations. These interpretations provided the basis for a flow law that incorporates hardening by the development of back stress. Here, we test this dislocation-interaction hypothesis by examining the microstructures of olivine samples deformed plastically at room temperature either in a deformation-DIA apparatus at differential stresses of ≤4.3GPa or in a nanoindenter at applied contact stresses of ≥10.2GPa. High-angular resolution electron backscatter diffraction maps reveal the presence of geometrically necessary dislocations with densities commonly above 1014m−2 and intragranular heterogeneities in residual stress on the order of 1 GPa in both sets of samples. Scanning transmission electron micrographs reveal straight dislocations aligned in slip bands and interacting with dislocations of other types that act as obstacles. The resulting accumulations of dislocations in their slip planes, and associated stress heterogeneities, are consistent with strain hardening resulting from long-range back-stresses acting among dislocations and thereby support the form of the flow law for low-temperature plasticity. Based on these observations, we predict that back stresses among dislocations will impart significant mechanical anisotropy to deformed lithosphere by enhancing or reducing the effective stress. Therefore, strain history, with associated microstructural and micromechanical evolution, is an important consideration for models of lithospheric strength. The microstructural observations also provide new criteria for identifying the operation of back-stress induced strain hardening in natural samples and therefore provide a means to test the applicability of the flow law for low-temperature plasticity.This research was supported by Natural Environment Research Council grants NE/M000966/1 to LNH, AJW, and DW and 1710DG008/JC4 to LNH and AJW; European Plate Observing System Transnational Access grant EPOS-TNA-MSL 2018-022 to LNH; Advanced Photon Source General User Proposal 55176 to LNH, DLG, and WBD; and National Science Foundation Awards EAR-1361319 to WBD, EAR-1625032 to JMW, and EAR-1806791 to KMK

A Tsunami Generated by a Strike-Slip Event: Constraints From GPS and SAR Data on the 2018 Palu Earthquake

A devastating tsunami struck Palu Bay in the wake of the 28 September 2018 Mw = 7.5 Palu earthquake (Sulawesi, Indonesia). With a predominantly strike-slip mechanism, the question remains whether this unexpected tsunami was generated by the earthquake itself, or rather by earthquake-induced landslides. In this study we examine the tsunami potential of the co-seismic deformation. To this end, we present a novel geodetic data set of Global Positioning System and multiple Synthetic Aperture Radar-derived displacement fields to estimate a 3D co-seismic surface deformation field. The data reveal a number of fault bends, conforming to our interpretation of the tectonic setting as a transtensional basin. Using a Bayesian framework, we provide robust finite fault solutions of the co-seismic slip distribution, incorporating several scenarios of tectonically feasible fault orientations below the bay. These finite fault scenarios involve large co-seismic uplift (>2 m) below the bay due to thrusting on a restraining fault bend that connects the offshore continuation of two parallel onshore fault segments. With the co-seismic displacement estimates as input we simulate a number of tsunami cases. For most locations for which video-derived tsunami waveforms are available our models provide a qualitative fit to leading wave arrival times and polarity. The modeled tsunamis explain most of the observed runup. We conclude that co-seismic deformation was the main driver behind the tsunami that followed the Palu earthquake. Our unique geodetic data set constrains vertical motions of the sea floor, and sheds new light on the tsunamigenesis of strike-slip faults in transtensional basins