Coseismic surface deformation in large earthquakes is typically measured using field mapping and with a range of geodetic methods (e.g., InSAR, lidar differencing, and GPS). Current methods, however, either fail to capture patterns of near-field coseismic surface deformation or lack preevent data. Consequently, the characteristics of off-fault deformation and the parameters that control it remain poorly understood. We develop a standardized method to fully measure the surface, near-field, coseismic deformation patterns at high resolution using the COSI-Corr program by correlating pairs of aerial photographs taken before and after the 1992 M w 7.3 Landers earthquake. COSI-Corr offers the advantage of measuring displacement across the entire zone of surface deformation and over a wider aperture than that available to field geologists. For the Landers earthquake, our measured displacements are systematically larger than the field measurements, indicating the presence of off-fault deformation. We show that 46% of the total surface displacement occurred as off-fault deformation, over a mean deformation width of 154 m. The magnitude and width of off-fault deformation along the rupture is primarily controlled by the macroscopic structural complexity of the fault system, with a weak correlation with the type of near-surface materials through which the rupture propagated. Both the magnitude and width of distributed deformation are largest in stepovers, bends, and at the southern termination of the surface rupture. We find that slip along the surface rupture exhibits a consistent degree of variability at all observable length scales and that the slip distribution is self-affine fractal with dimension of 1.56.

Supershear earthquakes with rupture velocity exceeding shear-wave speeds, previously observed in laboratory experiments and large strike-slip events, often have an initial sub-shear stage before they transition to supershear. In this study, integrated geophysical observations of the 2018 Mw 7.5 Palu, Indonesia earthquake, provide robust evidence of an early and persistent supershear rupture speed. Slowness-enhanced back-projection (SEBP) of teleseismic data provides a sharp image of the rupture process, consistently across multiple arrays. The inferred rupture path agrees with the surface rupture trace inferred from the net surface displacement field derived by sub-pixel InSAR image correlation. The SEBP results indicate a sustained rupture velocity of 4.1 km/s from the rupture initiation to the end, despite large fault bends. The persistent supershear speed is further validated by evidence of far-field Rayleigh Mach waves in regional seismograms. The short or absent supershear transition distance can be caused by high initial shear stress or short critical slip-weakening distance, and promoted by fault roughness near the hypocenter. Steady rupture propagation at a supershear speed considered to be unstable, lower than the Eshelby speed, could result from the presence of a damaged fault zone.

S U M M A R YGeodetic slip inversions for three major (M w > 7) strike-slip earthquakes (1992( Landers, 1999 Hector Mine and 2010 El Mayor-Cucapah) show a 15-60 per cent reduction in slip near the surface (depth < 2 km) relative to the slip at deeper depths (4-6 km). This significant difference between surface coseismic slip and slip at depth has been termed the shallow slip deficit (SSD). The large magnitude of this deficit has been an enigma since it cannot be explained by shallow creep during the interseismic period or by triggered slip from nearby earthquakes. One potential explanation for the SSD is that the previous geodetic inversions lack data coverage close to surface rupture such that the shallow portions of the slip models are poorly resolved and generally underestimated. In this study, we improve the static coseismic slip inversion for these three earthquakes, especially at shallow depths, by: (1) including data capturing the near-fault deformation from optical imagery and SAR azimuth offsets; (2) refining the interferometric synthetic aperture radar processing with non-boxcar phase filtering, modeldependent range corrections, more complete phase unwrapping by SNAPHU (Statistical Nonlinear Approach for Phase Unwrapping) assuming a maximum discontinuity and an onfault correlation mask; (3) using more detailed, geologically constrained fault geometries and (4) incorporating additional campaign global positioning system (GPS) data. The refined slip models result in much smaller SSDs of 3-19 per cent. We suspect that the remaining minor SSD for these earthquakes likely reflects a combination of our elastic model's inability to fully account for near-surface deformation, which will render our estimates of shallow slip minima, and potentially small amounts of interseismic fault creep or triggered slip, which could 'make up' a small percentages of the coseismic SSD during the interseismic period. Our results indicate that it is imperative that slip inversions include accurate measurements of near-fault surface deformation to reliably constrain spatial patterns of slip during major strike-slip earthquakes.

, a M w 6.5 earthquake, followed 34 h later by a M w 7.1 event, struck Searles Valley, California. These events are part of a long-lived cluster of historical earthquakes along the Eastern California Shear Zone (ECSZ) which started in 1872 and are associated with temporarily elevated strain rates. We find that the M w 6.5 event initiated on a right-lateral NW striking fault and then ruptured a left-lateral fault to the surface. This event triggered right-lateral slip during the M w 7.1 earthquake. It started as a bilateral, crack-like rupture on a segment brought closer to failure by the M w 6.5 event. The rupture evolved to pulse-like as it propagated at a relatively slow velocity (2 km/s) along a segment that was unloaded by the M w 6.5 event. It stopped abruptly at the Coso volcanic area and at the Garlock Fault and brought some neighbouring faults closer to failure.