Search

The variation of stress on faults is important for our understanding of fault friction and the dynamics of earthquake ruptures. However, we still have little observational constraints on their absolute magnitude, or their variations in space and in time over the seismic cycle. Here we use a new geodetic imaging technique to measure the 3D coseismic slip vectors along the 2019 Ridgecrest surface ruptures and invert them for the coseismic stress state. We find that the coseismic stresses show an eastward rotation that becomes increasingly transtensional from south‐to‐north along the rupture, that matches the known background stress state. We find that the main fault near the Mw 7.1 mainshock hypocenter was critically stressed. Coseismic slip was maximum there and decreased gradually along strike as the fault became less optimally oriented due its curved geometry. The variations of slip and stress along the curved faults are used to infer the static and dynamic fault friction assuming Mohr‐Coulomb failure. We find shear stresses of 4–9 MPa in the shallow crust (∼1.3 km depth) and that fault friction drops from a static, Byerlee‐type, value of 0.61 ± $\pm $ 0.14 to a dynamic value of 0.29 ± $\pm $ 0.04 during seismic slip. These values explain quantitatively the slip variations along a transpressional fault bend. Plain Language Summary: Understanding the orientation and magnitude of stresses within the crust are important because they can affect the location, size and spatial extent of earthquake rupture. However, measuring the absolute magnitude and orientation of stresses as well as the frictional properties of the fault surface (i.e., how strongly the fault resists the applied driving forces) is very difficult. Here we use optical images acquired by satellites to measure how the surface deformed in 3D during the 2019 Ridgecrest event. These 3D measurements allow us to extract the direction of fault slip movement along the entire rupture length which we use to estimate the direction of stresses by assuming the shear stress is parallel to the direction of the observed fault slip motion. We find that the main fault near the mainshock epicenter was the most optimally aligned for failure, which could be one contributing reason for the location of rupture initiation. By deriving a relation between how much a fault slips with how well aligned it is to the stress field we can estimate the absolute magnitude of stresses, the frictional resistance at initial fault sliding (finding a static friction = 0.61) and during sliding (a dynamic friction = 0.29). Key Points: Inverting surface coseismic slip vectors show a variable stress state that matches the background stresses constrained by seismicityFaults at the mainshock epicenter were the most critically stressed; we find slip increases linearly as faults become more optimally alignedWe find absolute stress magnitudes of 8–26 MPa in the upper crust, a static frictional coefficient of 0.61 and a dynamic value of 0.29 [ABSTRACT FROM AUTHOR]

We report Global Positioning System (GPS) measurements of postseismic deformation following the 2015 Mw7.8 Gorkha (Nepal) earthquake, including previously unpublished data from 13 continuous GPS stations installed in southern Tibet shortly after the earthquake. We use variational Bayesian Independent Component Analysis (vbICA) to extract the signal of postseismic deformation from the GPS time series, revealing a broad displacement field extending >150 km northward from the rupture. Kinematic inversions and dynamic forward models show that these displacements could have been produced solely by afterslip on the Main Himalayan Thrust (MHT) but would require a broad distribution of afterslip extending similarly far north. This would require the constitutive parameter (a − b)σ to decrease northward on the MHT to ≤0.05 MPa (an extreme sensitivity of creep rate to stress change) and seems unlikely in light of the low interseismic coupling and high midcrustal temperatures beneath southern Tibet. We conclude that the northward reach of postseismic deformation more likely results from distributed viscoelastic relaxation, possibly in a midcrustal shear zone extending northward from the seismogenic MHT. Assuming a shear zone 5–20 km thick, we estimate an effective shear‐zone viscosity of ~3·1016–3·1017 Pa·s over the first 1.12 postseismic years. Near‐field deformation can be more plausibly explained by afterslip itself and implies (a − b)σ ~ 0.5–1 MPa, consistent with other afterslip studies. This near‐field afterslip by itself would have re‐increased the Coulomb stress by ≥0.05 MPa over >30% of the Gorkha rupture zone in the first postseismic year, and deformation further north would have compounded this reloading. Key Points: New GPS data show that postseismic deformation of the 2015 Mw7.8 Gorkha (Nepal) earthquake extended far northward into southern TibetKinematic and dynamic models suggest that the northward reach of deformation likely reflects the influence of viscoelastic relaxationAfterslip likely dominated postseismic deformation in the near field, transferring stress back to the (relocked) Gorkha rupture zone [ABSTRACT FROM AUTHOR]

Analyzing the displacement time series from continuous GPS (cGPS) with an Independent Component Analysis, we detect a transient deformation signal that correlates both in space and time with a seismic swarm activity (maximum M w=3.69 ± 0.09) occurred in the hanging wall of the Altotiberina normal fault (Northern Apennines, Italy) in 2013-2014. The geodetic transient lasted ∼6 months and produced a NW-SE trending extension of ∼5.3 mm, consistent with the regional tectonic regime. The seismicity and the geodetic signal are consistent with slip on two splay faults in the Altotiberina fault (ATF) hanging wall. Comparing the seismic moment associated with the geodetic transient and the seismic events, we observe that seismicity accounts for only a fraction of the measured geodetic deformation. The combined seismic and aseismic slip decreased the Coulomb stress on the locked shallow portion of the ATF, while the transition region to the creeping section has been loaded. [ABSTRACT FROM AUTHOR]

Measurement of interseismic strain along subduction zones reveals the location of both locked asperities, which might rupture during megathrust earthquakes, and creeping zones, which tend to arrest such seismic ruptures. The heterogeneous pattern of interseismic coupling presumably relates to spatial variations of frictional properties along the subduction interface and may also show up in the fore-arc morphology. To investigate this hypothesis, we compiled information on the extent of earthquake ruptures for the last 500 years and uplift rates derived from dated marine terraces along the South American coastline from central Peru to southern Chile. We additionally calculated a new interseismic coupling model for that same area based on a compilation of GPS data. We show that the coastline geometry, characterized by the distance between the coast and the trench; the latitudinal variations of long-term uplift rates; and the spatial pattern of interseismic coupling are correlated. Zones of faster and long-term permanent coastal uplift, evidenced by uplifted marine terraces, coincide with peninsulas and also with areas of creep on the megathrust where slip is mostly aseismic and tend to arrest seismic ruptures. We conclude that spatial variations of frictional properties along the megathrust dictate the tectono-geomorphological evolution of the coastal zone and the extent of seismic ruptures along strike. [ABSTRACT FROM AUTHOR]

We model evolution of a continent-continent collision and draw some parallels with the tectonic evolution of the Himalaya. We use a large-scale visco-plasto-elastic thermomechanical model that has a free upper surface, accounts for erosion and deposition and allows for all modes of lithospheric deformation. For quartz/olivine rheology and 60 mm/yr convergence rate, the continental subduction is stable, and the model predicts three distinct phases. During the phase 1 (120 km or 6% of shortening), deformation is characterized by back thrusting around the suture zone. Some amount of delaminated lower crust accumulates at depth. During phase 2 (120 km-420 km or 6%-22% of shortening), this crustal root is exhumed (medium- to high-grade rocks) along a newly formed major thrust fault. This stage bears similarities with the period of coeval activity of the Main Central thrust and of the South Tibetan Detachment between 20-16 Myr ago. During phase 3 (>420 km or 22% of shortening), the crust is scraped off from the mantle lithosphere and is incorporated into large crustal wedge. Deformation is localized around frontal thrust faults. This kinematics should produce only low- to medium-grade exhumation. This stage might be compared with the tectonics that has prevailed in the Himalaya over the last 15 Myr allowing for the formation of the Lesser Himalaya. The experiment is conducted at constant convergence rate, which implies increasing compressive force. Considering that this force is constant in nature, this result may be equivalent to a slowing down of the convergence rate as was observed during the India-Asia collision. [ABSTRACT FROM AUTHOR]

The Lesser Himalaya (LH) consists of metasedimentary rocks that have been scrapped off from the underthrusting Indian crust and accreted to the mountain range over the last ∼20 Myr. It now forms a significant fraction of the Himalayan collisional orogen. We document the kinematics and thermal metamorphism associated with the deformation and exhumation of the LH, combining thermometric and thermochronological methods with structural geology. Peak metamorphic temperatures estimated from Raman spectroscopy of carbonaceous material decrease gradually from 520°-550°C below the Main Central Thrust zone down to less than 330°C. These temperatures describe structurally a 20°-50°C/km inverted apparent gradient. The Ar muscovite ages from LH samples and from the overlying crystalline thrust sheets all indicate the same regular trend; i.e., an increase from about 3-4 Ma near the front of the high range to about 20 Ma near the leading edge of the thrust sheets, about 80 km to the south. This suggests that the LH has been exhumed jointly with the overlying nappes as a result of overthrusting by about 5 mm/yr. For a convergence rate of about 20 mm/yr, this implies underthrusting of the Indian basement below the Himalaya by about 15 mm/yr. The structure, metamorphic grade and exhumation history of the LH supports the view that, since the mid-Miocene, the Himalayan orogen has essentially grown by underplating, rather than by frontal accretion. This process has resulted from duplexing at a depth close to the brittle-ductile transition zone, by southward migration of a midcrustal ramp along the Main Himalayan Thrust fault, and is estimated to have resulted in a net flux of up to 150 m2/yr of LH rocks into the Himalayan orogenic wedge. The steep inverse thermal gradient across the LH is interpreted to have resulted from a combination of underplating and post metamorphic shearing of the underplated units. [ABSTRACT FROM AUTHOR]

The combination of detailed topographic leveling on the southwest segment of the El Asnam thrust fault with existing seismic and geologic data implies that the geometry of this fault involves shallow dipping flats and steep ramps. The fault appears to be growing along strike toward the southwest end, where the main shock initiated in 1980. From a depth of about 10 km, the main thrust appears to ramp to the basement-Cenozoic cover interface on a plane striking N40°E and dipping 50°-55° to the northwest. Along the southwest segment where folding has not yet developed, the thrust continues steeply through the Cenozoic cover to the near surface where it flattens, causing normal faulting. Along the central and northeast segments, which display a more evolved fold structure, the deep thrust probably flattens at a depth of 5-6 km, into a decollement along the Cenozoic-Jurassic interface before ramping to the surface. The Sara El Marouf and Kef El Mes anticlines have thus formed as fault propagation folds. Normal faults at Beni Rached probably branch with the thrust to maintain kinematic compatibility between the deep ramp and decollement. The greater separation (∼7 km) between the normal faults at Beni Rached and the thrust where it crosses Oued Cheliff than along the southwest segment (∼1 km) reflects the greater depth of the ramp to flat bend. We infer that the September 9, 1954, earthquake activated only the central deep segment of the main thrust together with the Beni Rached normal faults, while that of October 10, 1980, activated the whole system of flat decollements, ramp thrusts and compatibility normal faults. Further complexities of the faulting in map view are related to changes of strike of the thrust (in particular north of Oued Cheliff). [ABSTRACT FROM AUTHOR]

Deformation of the Earth's surface associated with redistributions of continental water mass explains, to first order, the seasonal signals observed in geodetic position time series. Discriminating these seasonal signals from other sources of deformation in geodetic measurements is essential to isolate tectonic signals and to monitor spatio‐temporal variations in continental water storage. We propose a new methodology to identify and extract these seasonal signals. The approach uses a variational Bayesian Independent Component Analysis (vbICA) to extract the seasonal signals and a gravity‐based deformation model to identify which of these signals are caused by surface loading. We test the procedure on two study areas, the Arabian Peninsula and the Nepal Himalaya, and find that the technique successfully extracts the seasonal signals with one or two independent components, depending on whether the load is stationary or moving. The approach is robust to spatial heterogeneities inherent to geodetic measurements and can help extract systematic errors in geodetic products (e.g., draconitic errors). We also discuss how to handle the degree‐1 deformation field present in the geodetic data set but not captured by the gravity‐based model. Key Points: We propose a procedure to isolate seasonal geodetic signals due to surface loads using a GRACE‐based deformation model and an independent component analysisGNSS and GRACE‐derived displacements yield consistent seasonal time evolutions in both the Arabian Peninsula and the Nepal HimalayaThe method is robust to topographic and elastic heterogeneities, simple moving loads and geodetic products errors [ABSTRACT FROM AUTHOR]