Search

The Eastern Mediterranean is the most seismically active region in Europe due to the complex interactions of the Arabian, African, and Eurasian tectonic plates. Deformation is achieved by faulting in the brittle crust, distributed flow in the viscoelastic lower-crust and mantle, and Hellenic subduction, but the long-term partitioning of these mechanisms is still unknown. We exploit an extensive suite of geodetic observations to build a kinematic model connecting strike-slip deformation, extension, subduction, and shear localization across Anatolia and the Aegean Sea by mapping the distribution of slip and strain accumulation on major active geological structures. We find that tectonic escape is facilitated by a plate-boundary-like, trans-lithospheric shear zone extending from the Gulf of Evia to the Turkish-Iranian Plateau that underlies the surface trace of the North Anatolian Fault. Additional deformation in Anatolia is taken up by a series of smaller-scale conjugate shear zones that reach the upper mantle, the largest of which is located beneath the East Anatolian Fault. Rapid north–south extension in the western part of the system, driven primarily by Hellenic Trench retreat, is accommodated by rotation and broadening of the North Anatolian mantle shear zone from the Sea of Marmara across the north Aegean Sea, and by a system of distributed transform faults and rifts including the rapidly extending Gulf of Corinth in central Greece and the active grabens of western Turkey. Africa–Eurasia convergence along the Hellenic Arc occurs at a median rate of 49.8 mm yr–1 in a largely trench-normal direction except near eastern Crete where variably oriented slip on the megathrust coincides with mixed-mode and strike-slip deformation in the overlying accretionary wedge near the Ptolemy–Pliny–Strabo trenches. Our kinematic model illustrates the competing roles the North Anatolian mantle shear zone, Hellenic Trench, overlying mantle wedge, and active crustal faults play in accommodating tectonic indentation, slab rollback and associated Aegean extension. Viscoelastic flow in the lower crust and upper mantle dominate the surface velocity field across much of Anatolia and a clear transition to megathrust-related slab pull occurs in western Turkey, the Aegean Sea and Greece. Crustal scale faults and the Hellenic wedge contribute only a minor amount to the large-scale, regional pattern of Eastern Mediterranean interseismic surface deformation. [ABSTRACT FROM AUTHOR]

Constitutive models of fault friction form the basis of physics-based simulations of seismic activity. A generally accepted framework for the slip-rate and state dependence of friction involves a thermally activated process, whereby the probability of slip along microasperities adheres to an Arrhenius law. This model, which has become widely adopted among experimentalists and theoreticians, predicts a continuous increase of the direct effect with absolute temperature, but is it observed experimentally? Leveraging comprehensive laboratory data across diverse hydrothermal, barometric, and lithological conditions, we demonstrate that, contrary to the classical view, the direct effect for a given deformation mechanism remains largely temperature-independent. Instead, the incremental shifts in the direct effect often coincide with the brittle to semi-brittle transition, across which distinct deformation mechanisms operate. These considerations challenge the validity of the classical model. Realistic constitutive laws for rock failure within the lithosphere must incorporate the contributions of multiple deformation mechanisms within active fault zones., Competing Interests: Competing interests statement:The author declares no competing interest.

Although the surface deformation of tectonic plate boundaries is well determined by geological and geodetic measurements, the pattern of flow below the lithosphere remains poorly constrained. We use the crustal velocity field of the Plate Boundary Observatory to illuminate the distribution of horizontal flow beneath the California margin. At lower-crustal and upper-mantle depths, the boundary between the Pacific and North American plates is off-centered from the San Andreas fault, concentrated in a region that encompasses the trace of nearby active faults. A major step is associated with return flow below the Eastern California Shear Zone, leading to the extrusion of the Mojave block and a re-distribution of fault activity since the Pleistocene. Major earthquakes in California have occurred above the regions of current plastic strain accumulation. Deformation is mechanically coupled from the crust to the asthenosphere, with mantle flow overlaid by a kinematically consistent network of faults in the brittle crust.

The relative motion of tectonic plates is accommodated at boundary faults through slow and fast ruptures that encompass a wide range of source properties. Near the Parkfield segment of the San Andreas fault, low-frequency earthquakes and slow-slip events take place deeper than most seismicity, at temperature conditions typically associated with stable sliding. However, laboratory experiments indicate that the strength of granitic gouge decreases with increasing temperature above 350°C, providing a possible mechanism for weakening if temperature is to vary dynamically. Here, we argue that recurring low-frequency earthquakes and slow-slip transients at these depths may arise because of shear heating and the temperature dependence of frictional resistance. Recurring thermal instabilities can explain the recurrence pattern of the mid-crustal low-frequency earthquakes and their correlative slip distribution. Shear heating associated with slow slip is sufficient to generate pseudotachylyte veins in host rocks even when fault slip is dominantly aseismic., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).)

Deformation associated with plate convergence at subduction zones is accommodated by a complex system involving fault slip and viscoelastic flow. These processes have proven difficult to disentangle. The 2010 M w 8.8 Maule earthquake occurred close to the Chilean coast within a dense network of continuously recording Global Positioning System stations, which provide a comprehensive history of surface strain. We use these data to assemble a detailed picture of a structurally controlled megathrust fault frictional patchwork and the three-dimensional rheological and time-dependent viscosity structure of the lower crust and upper mantle, all of which control the relative importance of afterslip and viscoelastic relaxation during postseismic deformation. These results enhance our understanding of subduction dynamics including the interplay of localized and distributed deformation during the subduction zone earthquake cycle., (Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY).)

The strength of the lithosphere controls tectonic evolution and seismic cycles, but how rocks deform under stress in their natural settings is usually unclear. We constrain the rheological properties beneath the Taiwan orogenic belt using the stress perturbation following the 1999 Chi-Chi earthquake and fourteen-year postseismic geodetic observations. The evolution of stress and strain rate in the lower crust is best explained by a power-law Burgers rheology with rapid increases in effective viscosities from ~10 17 to ~10 19 Pa s within a year. The short-term modulation of the lower-crustal strength during the seismic cycle may alter the energy budget of mountain building. Incorporating the laboratory data and associated uncertainties, inferred thermal gradients suggest an eastward increase from 19.5±2.5°C/km in the Coastal Plain to 32±3°C/km in the Central Range. Geodetic observations may bridge the gap between laboratory and lithospheric scales to investigate crustal rheology and tectonic evolution.

Plate tectonics is one mode of mantle convection that occurs when the surface layer (the lithosphere) is relatively weak. When plate tectonics operates on a terrestrial planet, substantial exchange of materials occurs between planetary interior and its surface. This is likely a key in maintaining the habitable environment on a planet. Therefore it is essential to understand under which conditions plate tectonics operates on a terrestrial planet. One of the puzzling observations in this context is the fact that plate tectonics occurs on Earth but not on Venus despite their similar size and composition. Factors such as the difference in water content or in grain-size have been invoked, but these models cannot easily explain the contrasting tectonic styles between Earth and Venus. We propose that strong dynamic weakening in friction is a key factor. Fast unstable fault motion is found in cool Earth, while slow and stable fault motion characterizes hot Venus, leading to substantial dynamic weakening on Earth but not on Venus. Consequently, the tectonic plates are weak on Earth allowing for their subduction, while the strong plates on Venus promote the stagnant lid regime of mantle convection.

Surveillance of clandestine nuclear tests relies on a global seismic network, but the potential of spaceborne monitoring has been underexploited. We used satellite radar imagery to determine the complete surface displacement field of up to 3.5 meters of divergent horizontal motion with 0.5 meters of subsidence associated with North Korea's largest underground nuclear test. Combining insight from geodetic and seismological remote sensing, we found that the aftermath of the initial explosive deformation involved subsidence associated with subsurface collapse and aseismic compaction of the damaged rocks of the test site. The explosive yield from the nuclear detonation with best-fitting source parameters for 450-meter depth was 191 kilotonnes of TNT equivalent. Our results demonstrate the capability of spaceborne remote sensing to help characterize large underground nuclear tests., (Copyright © 2018 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Understanding the rheological properties of the upper mantle is essential to develop a consistent model of mantle dynamics and plate tectonics. However, the spatial distribution and temporal evolution of these properties remain unclear. Here, we infer the rheological properties of the asthenosphere across multiple great megathrust earthquakes between 2004 and 2014 along the Sumatran subduction zone, taking advantage of decade-long continuous GPS and tide-gauge measurements. We observe transient mantle wedge flow following these earthquakes, and infer the temporal evolution of the effective viscosity. We show that the evolution of stress and strain rate following these earthquakes is better matched by a bi-viscous than by a power-law rheology model, and we estimate laterally heterogeneous transient and background viscosities on the order of ~10 17 and ~10 19  Pa s, respectively. Our results constitute a preliminary rheological model to explain stress evolution within earthquake cycles and the development of seismic hazard in the region.

The deformation of mantle and crustal rocks in response to stress plays a crucial role in the distribution of seismic and volcanic hazards, controlling tectonic processes ranging from continental drift to earthquake triggering. However, the spatial variation of these dynamic properties is poorly understood as they are difficult to measure. We exploited the large stress perturbation incurred by the 2016 earthquake sequence in Kumamoto, Japan, to directly image localized and distributed deformation. The earthquakes illuminated distinct regions of low effective viscosity in the lower crust, notably beneath the Mount Aso and Mount Kuju volcanoes, surrounded by larger-scale variations of viscosity across the back-arc. This study demonstrates a new potential for geodesy to directly probe rock rheology in situ across many spatial and temporal scales., (Copyright © 2017, American Association for the Advancement of Science.)

The deformation at well-defined, narrow plate boundaries depends on the relative plate motion, but how the deformation takes place within a distributed plate boundary zone remains a conundrum. This was confirmed by the seismological analyses of the 2012 great Wharton Basin earthquakes [moment magnitude ( M w ) 8.6], which suggested the rupture of several faults at high angles to one another. Using high-resolution bathymetry and seismic reflection data, we report the discovery of new N294°E-striking shear zones, oblique to the plate fabric. These shear zones are expressed by sets of normal faults striking at N335°E, defining the direction of the principal compressional stress in the region. Also, we have imaged left-lateral strike-slip faults along reactivated N7°E-oriented oceanic fracture zones. The shear zones and the reactivated fracture zones form a conjugate system of faults, which accommodate present-day intraplate deformation in the Wharton Basin.

The deep extension of the San Andreas Fault is believed to be creeping, but the recent observations of tectonic tremors from these depths indicate a complex deformation style. In particular, an isolated tremor source near Parkfield has been producing a sequence of low-frequency earthquakes that indicates an uncommon mechanism of stress accumulation and release. The tremor pattern regularly oscillated between three and six days from mid-2003 until it was disrupted by the 2004 magnitude 6.0 Parkfield earthquake. After that event, the tremor source ruptured only about every three days, but over the next two years it gradually returned to its initial alternating recurrence pattern. The mechanism that drives this recurrence pattern is unknown. Here we use physics-based models to show that the same tremor asperity--the region from which the low-frequency earthquakes radiate--can regularly slip in slow and fast ruptures, naturally resulting in recurrence intervals alternating between three and six days. This unusual slip behaviour occurs when the tremor asperity size is close to the critical nucleation size of earthquakes. We also show that changes in pore pressure following the Parkfield earthquake can explain the sudden change and gradual recovery of the recurrence intervals. Our findings suggest a framework for fault deformation in which the same asperity can release tectonic stress through both slow and fast ruptures.

Advances in observational, laboratory, and modeling techniques open the way to the development of physical models of the seismic cycle with potentially predictive power. To explore that possibility, we developed an integrative and fully dynamic model of the Parkfield segment of the San Andreas Fault. The model succeeds in reproducing a realistic earthquake sequence of irregular moment magnitude (M(w)) 6.0 main shocks--including events similar to the ones in 1966 and 2004--and provides an excellent match for the detailed interseismic, coseismic, and postseismic observations collected along this fault during the most recent earthquake cycle. Such calibrated physical models provide new ways to assess seismic hazards and forecast seismicity response to perturbations of natural or anthropogenic origins.