Coseismic slip variation assessed from terrestrial lidar scans of the El Mayor–Cucapah surface rupture

We analyze high-resolution (>10³points/m²) terrestrial lidar surveys of the 4 April 2010 El Mayor–Cucapah earthquake rupture (Baja California, Mexico), collected at three sites 12–18 days after the event. Using point cloud-based tools in an immersive visualization environment, we quantify coseismic fault slip for hundreds of meters along strike and construct densely constrained along-strike slip distributions from measurements of offset landforms. Uncertainty bounds for each offset, determined empirically by repeatedly measuring offsets at each site sequentially, illuminate measurement uncertainties that are difficult to quantify in the field. These uncertainties are used to define length scales over which variability in slip distributions may be assumed to reflect either recognizable earthquake mechanisms or measurement noise. At two sites characterized by 2–3m of concentrated right-oblique slip, repeat measurements yield 2σ uncertainties of ±11–12%. Each site encompasses ∼200m along strike, and a smoothed linear slip gradient satisfies all measurement distributions, implying along-fault strains of ∼10⁻³. Conversely, the common practice of defining the slip curve by the local slip maxima distorts the curve, overestimates along-fault strain, and may overestimate actual fault slip by favoring measurements with large, positive, uncertainties. At a third site characterized by 1–2.5m of diffuse normal slip, repeat measurements of fault throw summed along fault-perpendicular profiles yield 2σ uncertainties of ±17%. Here, a low order polynomial fit through the measurement averages best approximates surface slip. However independent measurements of off-fault strain accommodated by hanging wall flexure suggest that over the ∼200m length of this site, a linear interpolation through the average values for the slip maxima at either end of this site most accurately represents subsurface displacement. In aggregate, these datasets show that given uncertainties of greater than ±11% (2σ), slip distributions over shorter scales are likely to be less uneven than those derived from a single set of field- or lidar-based measurements. This suggests that the relatively smooth slip curves we obtain over ∼10²m distances reflect real physical phenomena, whereas short wavelength variability over ∼10⁰–10¹m distances can be attributed to measurement uncertainty.