Analyzing the horizontal orientation of the crustal stress adjacent to plate boundaries

The spatial analysis of horizontal stress orientation is important to study stress sources and understand tectonics and the deformation of the lithosphere. Additional to the stress sources, the geometry of stress fields depends on the underlying coordinate reference system, which causes spatial distortions that bias the analysis and interpretation of stresses. The bias can be avoided when the stress field is decomposed and transformed into the reference frame of its first-order stress source. We present a modified and extended theory based on the empirical link between the orientation of first-order stresses and the trajectories of lateral plate boundary forces. This link is applied to analyze the orientation of horizontal stresses, their patterns, and tectonic structures from the perspective of their first-order source or cause. By using only parameters for the relative motion between two neighboring plates, we model the first-order orientation of the maximum horizontal stress that statistically fits the orientation of [Formula: see text]80% of the global stress data adjacent to plate boundaries. Considerable deviations of the observed stress from the predicted first-order stress direction can reveal the geometry of second-order stresses and confine areas where other stress sources dominate. The model's simple assumptions, independence from the sample size, potential application to regional to global scale analysis, and compatibility with other spatial interpolation algorithms make it a powerful method for analyzing stress fields. For immediate use, the presented method is implemented in the free and open-source software package tectonicr, which is written in the computer language R.

Relations between earthquake distributions, geological history, tectonics and rheology on the continents

This paper is concerned with the distribution of earthquakes, particularly their depths, with the temperature of the material in which they occur, and with the significance of both for the rheology and deformation of the continental lithosphere. Earthquakes on faults are generated by the sudden release of elastic energy that accumulates during slow plate motions. The nonlinear high-temperature creep that localizes such energy accumulation is, in principle, well understood and can be described by rheological models. But the same is not true of seismogenic brittle failure, the main focus of this paper, and severely limits the insights that can be obtained by simulations derived from geodynamical modelling of lithosphere deformation. Through advances in seismic tomography, we can now make increasingly detailed maps of lithosphere thickness on the continents. The lateral variations are dramatic, with some places up to 300 km thick, and clearly relate to the geological history of the continents as well as their present-day deformation. Where the lithosphere thickness is about 120 km or less, continental earthquakes are generally confined to upper crustal material that is colder than about 350°C. Within thick lithosphere, and especially on its edges, the entire crust may be seismogenic, with earthquakes sometimes extending into the uppermost mantle if the Moho is colder than 600°C, but the continental mantle is generally aseismic. Earthquakes in the continental lower crust at 400-600°C require the crust to be anhydrous and so are a useful guide or proxy to both composition and strength. These patterns and correlations have important implications for the geological evolution of the continents. They can be seen to have influenced features as diverse as the location of post-collisional rifting; cratonic basin formation; the location, origin and timing of granulite-facies metamorphism; and the formation, longevity and strength of cratons. In addition, they have important consequences for earthquake hazard assessment in the slowly deforming edges and interiors of continental shields or platforms, where the large seismogenic thickness can host very large earthquakes. This article is part of a discussion meeting issue 'Understanding earthquakes using the geological record'.

Strike-slip 23 January 2018 MW 7.9 Gulf of Alaska rare intraplate earthquake: Complex rupture of a fracture zone system

Large intraplate earthquakes in oceanic lithosphere are rare and usually related to regions of diffuse deformation within the oceanic plate. The 23 January 2018 M_W 7.9 strike-slip Gulf of Alaska earthquake ruptured an oceanic fracture zone system offshore Kodiak Island. Bathymetric compilations show a muted topographic expression of the fracture zone due to the thick sediment that covers oceanic basement but the fracture zone system can be identified by offset N-S magnetic anomalies and E-W linear zones in the vertical gravity gradient. Back-projection from global seismic stations reveals that the initial rupture at first propagated from the epicenter to the north, likely rupturing along a weak zone parallel to the ocean crustal fabric. The rupture then changed direction to eastward directed with most energy emitted on Aka fracture zone resulting in an unusual multi-fault earthquake. Similarly, the aftershocks show complex behavior and are related to two different tectonic structures: (1) events along N-S trending oceanic fabric, which ruptured mainly strike-slip and additionally, in normal and oblique slip mechanisms and (2) strike-slip events along E-W oriented fracture zones. To explain the complex faulting behavior we adopt the classical stress and strain partitioning concept and propose a generalized model for large intra-oceanic strike-slip earthquakes of trench-oblique oriented fracture zones/ocean plate fabric near subduction zones. Taking the Kodiak asperity position of 1964 maximum afterslip and outer-rise Coulomb stress distribution into account, we propose that the unusual 2018 Gulf of Alaska moment release was stress transferred to the incoming oceanic plate from co- and post-processes of the nearby great 1964 M_W 9.2 megathrust earthquake.

Global Tectonics and Space Geodesy

Much of the success of plate tectonics can be attributed to the near rigidity of tectonic plates and the availability of data that describe the rates and directions of motion across narrow plate boundaries \m=~\1 to 60 kilometers wide. Nonetheless, many plate boundaries in both continental and oceanic lithosphere are not narrow but are hundreds to thousands of kilometers wide. Wide plate boundary zones cover \m=~\15 percent of Earth's surface area. Space geodesy, which includes very long baseline radio interferometry, satellite laser ranging, and the global positioning system, is providing the accurate long-distance measurements needed to estimate the present motion across and within wide plate boundary zones. Space geodetic data show that plate velocities averaged over years are remarkably similar to velocities averaged over millions of years.