Three‐dimensional near‐fault coseismic deformation fields from high‐resolution differential topography provide new information on the behavior of the shallow fault zone in large surface‐rupturing earthquakes. Our work focuses on the 16 April 2016M_w7.0 Kumamoto, Japan, earthquake, which ruptured ~40 km of the Futagawa‐Hinagu Fault Zone on Kyushu Island with an oblique strike‐slip mechanism and surface offset exceeding 2 m. Our differential lidar analysis constrains the structural style of strain accommodation along the primary fault trace and the surrounding damage zone. We show that 36 ± 29% and 62 ± 32% of the horizontal and vertical deformation, respectively, was accommodated off the principal fault trace. The horizontal strains of up to 0.03 suggest that the approximate elastic strain limit was exceeded over a ~250 m width in many locations along the rupture. The inelastic deformation of the fault volume produced the observed distributed deformation at the Earth's surface. We demonstrate a novel approach for calculating 3‐D displacement uncertainties, indicating errors of centimeters to a few decimeters for displacements computed over 50 m horizontal windows. Errors correlate with land cover and relief, with flatter agricultural land associated with the highest displacement uncertainty. These advances provide a framework for future analyses of shallow earthquake behavior using differential topography.

Observations of surface deformation within 1–2 km of a surface rupture contain invaluable information about the coseismic behavior of the shallow crust. We investigate the oblique strike‐slip 2016 M7 Kumamoto, Japan, earthquake, which ruptured the Futagawa‐Hinagu Fault. We solve for variable fault slip in an inversion of differential lidar topography, satellite optical image correlation, and Interferometric Synthetic Aperture Radar (InSAR)‐derived surface displacements. The near‐fault differential lidar pose several challenges. The model fault geometry must follow the surface trace at the sub‐kilometer scale. Integration of displacement datasets with different sensitivities to the 3D deformation field and varying spatial distribution permits additional complexity in the inferred slip but introduces ambiguity that requires careful selection of the regularization. We infer a M_w earthquake. The maximum slip of 6.9 m occurred at 4.5‐km depth, suggesting an on‐fault slip deficit in the upper several kilometers of the crust that likely reflects distributed and inelastic deformation within the shallow fault zone.

We present high-resolution mapping and surface faulting measurements along the Lost River fault (Idaho-USA), a normal fault activated in the 1983 (M_w6.9) earthquake. The earthquake ruptured ~35 km of the fault with a maximum throw of ~3 m. From new 5 to 30 cm-pixel resolution topography collected by an Unmanned Aerial Vehicle, we produce the most comprehensive dataset of systematically measured vertical separations from ~37 km of fault length activated by the 1983 and prehistoric earthquakes. We provide Digital Elevation Models, orthophotographs, and three tables of: (i) 757 surface rupture traces, (ii) 1295 serial topographic profiles spaced 25 m apart that indicate rupture zone width and (iii) 2053 vertical separation measurements, each with additional textual and numerical fields. Our novel dataset supports advancing scientific knowledge about this fault system, refining scaling laws of intra-continental faults, comparing to other earthquakes to better understand faulting processes, and contributing to global probabilistic hazard approaches. Our methodology can be applied to other fault zones with high-resolution topographic data.

Following observations made in a survey campaign along the Lost River Fault (Idaho, USA) in 2019, we integrate both original and previously published data to obtain a detailed segmentation of the fault sections that failed in the 1983 Borah Peak earthquake (⁠ $M$_w 6.9). The earthquake ruptured the topographic surface with an oblique-normal faulting mechanism, activating two SW-dipping fault segments (Thousand Springs and Warm Springs) and a branching SSW-dipping fault (Arentson Gulch Fault) and producing coseismic surface ruptures with up to 3 m of vertical separation. We augment the 1983 earthquake description by interpreting high-resolution topography and fault mapping. We use quality vertical separation data, rupture zone width measurements, and fault slip data to analyze major and minor structural-geometric complexities, highlighting a partition of the deformation and a fault segmentation into four detail levels (i.e., segments, sections, subsections, and sectors). Our work provides new details of the 1983 Borah Peak earthquake, constraints for paleoseismic and seismotectonic studies, and a methodological approach applicable in other areas of the world. Our fault-slip data show variations along fault-strike that we interpret as kinematic partitioning. In 1983, the main southern segment had a large rupture zone width, while the northern segment localized the deformation. The distributed ruptures accommodate a large portion of the rupture length (~19.5 km versus 31 km for the main rupture) and displacement (~66%). 83% of the surface faulting and 80% of the displacement are located at the hanging wall of the main rupture. There is a strong correlation between vertical separation, rupture zone width, rupture position (footwall or hanging wall), and fault geometry. We highlight the control of the obliquity and kinematic partitioning in the surface expression of the earthquake propagation. We interpret the coseismic (i.e., 1983) and long-term (i.e., Quaternary) behavior, showing that the two activated segments had similar cumulated behaviors in distributing the deformation between synthetic and antithetic ruptures, despite the different geometries. Our results have implications for fault rupture behavior with application to rupture hazard.

Imaging tectonic creep along active faults is critical for measuring strain accumulation and ultimately understanding the physical processes that guide creep and the potential for seismicity. We image tectonic deformation along the central creeping section of the San Andreas Fault at the Dry Lake Valley paleoseismic site (36.468°N, 121.055°W) using three data sets with varying spatial and temporal scales: (1) an Interferometric Synthetic Aperture Radar (InSAR) velocity field with an ~100‐km footprint produced from Sentinel‐1 satellite imagery, (2) light detection and ranging (lidar) and structure‐from‐motion 3‐D topographic differencing that resolves a decade of deformation over a 1‐km aperture, and (3) surface fractures that formed over the 3‐ to 4‐m wide fault zone during a drought from late 2012 to 2014. The InSAR velocity map shows that shallow deformation is localized to the San Andreas Fault. We demonstrate a novel approach for differencing airborne lidar and structure‐from‐motion topography that facilitates resolving deformation along and adjacent to the San Andreas Fault. The 40‐m resolution topographic differencing resolves a 2.5 ± 0.2 cm/yr slip rate localized to the fault. The opening‐mode fractures accommodate cm/yr of fault slip. A 90% ± 30% of the 1‐km aperture deformation is accommodated over the several meter‐wide surface trace of the San Andreas Fault. The extension direction inferred from the opening‐mode fractures and topographic differencing is 40°–48° from the local trend of the San Andreas Fault. The localization of deformation likely reflects the well‐oriented and mature fault.

Fault creep reduces seismic hazard and serves as a window into plate boundary processes; however, creep rates are typically constrained with sparse measurements. We use differential lidar topography (11–13 year time span) to measure a spatially dense surface deformation field along a 150 km section of the Central San Andreas and Calaveras faults. We use an optimized windowed‐iterative‐closest‐point approach to resolve independent creep rates every 400 m at 1–2 km apertures. Rates vary from <10 mm/year along the creeping fault ends to over 30 mm/year along much of the central 100 km of the fault. Creep rates are 3–8 mm/year higher than most rates from alignment arrays and creepmeters, likely due to the larger aperture of the topographic differencing. Creep is often focused along discrete fault traces, but strain is sometimes distributed in areas of complex fault geometry, such as Mustang Ridge. Our observations constrain shallow seismic moment accumulation and the location of the creeping fault trace.