Locating Fully Locked Asperities Along the South America Subduction Megathrust: A New Physical Interseismic Inversion Approach in a Bayesian Framework

The largest earthquakes in subduction zones occur where significant interseismic slip deficit has accumulated on the plate interface. Slip deficit accumulates most quickly in mechanically locked regions, and these also cause the regions around them to accumulate slip deficit; therefore, large earthquakes are typically expected to rupture in and around locked areas. The locations and dimensions of these locked zones have been difficult to resolve using standard techniques and available data sets. We develop a new statistical interseismic inversion approach that incorporates the physical interactions between nearby fault areas to directly determine the distribution of locking on the subduction plate interface (simultaneously with rigid forearc motions) from interseismic surface velocities. Because we include physical prior information in the inversion procedure, this approach reduces uncertainties in the rate of slip deficit accumulation, even in locations (such as near the trench) where kinematic inversions of onshore data have relatively low resolution. Applying the inversion to the South America subduction zone, we find that the pattern of locking and corresponding slip deficit rates correlate well with recent and historical large earthquake ruptures. Locked patch dimensions are <40 km and account for no more than 30% of the area of the plate interface. The small size of the imaged locked zones is a natural outcome of our physical assumptions and implies that mechanical locking is caused by correspondingly small geological features. Despite their small dimensions, locked zones generate substantial slip deficit on the surrounding plate interface, consistent with the slip patterns of large megathrust earthquakes.