Measuring fault offsets preserved at the ground surface is of primary importance to recover earthquake and long‐term slip distributions and understand fault mechanics. The recent explosion of high‐resolution topographic data, such as Lidar and photogrammetric digital elevation models, offers an unprecedented opportunity to measure dense collections of fault offsets. We have developed a new Matlab code, 3D_Fault_Offsets, to automate these measurements. In topographic data, 3D_Fault_Offsets mathematically identifies and represents nine of the most prominent geometric characteristics of common sublinear markers along faults (especially strike slip) in 3‐D, such as the streambed (minimum elevation), top, free face and base of channel banks or scarps (minimum Laplacian, maximum gradient, and maximum Laplacian), and ridges (maximum elevation). By calculating best fit lines through the nine point clouds on either side of the fault, the code computes the lateral and vertical offsets between the piercing points of these lines onto the fault plane, providing nine lateral and nine vertical offset measures per marker. Through a Monte Carlo approach, the code calculates the total uncertainty on each offset. It then provides tools to statistically analyze the dense collection of measures and to reconstruct the prefaulted marker geometry in the horizontal and vertical planes. We applied 3D_Fault_Offsets to remeasure previously published offsets across 88 markers on the San Andreas, Owens Valley, and Hope faults. We obtained 5,454 lateral and vertical offset measures. These automatic measures compare well to prior ones, field and remote, while their rich record provides new insights on the preservation of fault displacements in the morphology.

The Mw 7.8 2016 Kaikoura earthquake ruptured the Kekerengu-Needle fault resulting in the loading of its eastern continuation, the Wairarapa fault. Since the most recent earthquake on Wairarapa occurred in 1855 and is one of the strongest continental earthquakes ever observed, it is critical to assess the seismic potential of the Wairarapa fault, which might be prone to break. Using Lidar data, we examine its bare-earth morphology and reveal ~650 mostly undiscovered offset geomorphic markers. Using a code we developed in earlier work, we automatically measure the lateral and vertical offsets of these markers providing more than 7000 well constrained measurements. The data document the lateral and vertical slip profiles of the 1855 earthquake for the first time and show its total slip reached ~20 m at surface. Modeling the entire offset dataset reveals 7 prior earthquakes ruptured the entire fault, each similarly producing 16.9 ± 1.4 m dextral slip and ~0.6 m vertical slip at surface in the same central bend zone of the fault. Thus, the Wairarapa fault repeatedly produced giant earthquakes and is likely able to produce a similarly strong forthcoming event. The extreme large size of the Wairarapa earthquakes questions our understanding of earthquake physics.