Evolution and progressive geomorphic manifestation of surface faulting: A comparison of the Wairau and Awatere faults, South Island, New Zealand

Field mapping and lidar analysis of surface faulting patterns expressed in flights of geologically similar fluvial terraces at the well-known Branch River and Saxton River sites along the Wairau (Alpine) and Awatere strike-slip faults, South Island, New Zealand, reveal that fault-related deformation patterns expressed in the topography at these sites are markedly less structurally complex along the higher-displacement (hundreds of kilometers), structurally mature Wairau fault than along the Awatere fault (∼13–20 km total slip). These differences, which are generally representative of the surface traces of these faults, provide direct evidence that surface faulting becomes structurally simpler with increasing cumulative fault offset. We also examine the degree to which off-fault deformation (OFD) is expressed in the landscape at the Saxton River site along the less structurally mature Awatere fault. Significantly greater amounts of OFD are discernible as a wide damage zone (∼460 m fault-perpendicular width) in older (ca. 15 ka), more-displaced (64–74 m) fluvial terraces than in younger (ca. 1–7 ka), less-displaced (<55 m) terraces; no OFD is discernible in the lidar data on the least-displaced (<35 m) terraces. From this, we infer that OFD becomes progressively more geomorphically apparent with accumulating displacement. These observations imply that (1) the processes that accommodate OFD are active during each earthquake, but may not be evident in deposits that have experienced relatively small displacements; (2) structures accommodating OFD will become progressively geomorphically clearer with increasing displacement; (3) geomorphic measurements of overall fault zone width taken in deposits that have experienced small displacements will be underestimates; and (4) fault slip rates based on geomorphic surface offsets will be underestimates for immature faults if based solely on measurements along the high-strain fault core

One tune, many tempos: faults trade off slip in time and space to accommodate relative plate motions

Analysis of incremental slip rates from the four major strike-slip faults of the Marlborough fault system (MFS) of northern South Island, New Zealand, provides a first-ever record at the scale of an entire plate-boundary fault system of how relative plate motions are accommodated in time and space. This record, which spans the past 350–450 m of relative plate motion and ca. 12–14 ky, demonstrates that the fault system as a whole accommodates a steady plate-boundary slip rate, with the MFS faults “keeping up” with the overall rate of relative Pacific-Australia plate motion at relatively short displacement (10 s of meters) and time (102–103 yr) scales. These results affirm the often-assumed but until now unproven assumption that the relative plate-motion rate provides a robust basic constraint on both geodynamical models and analyses of system-level seismic hazard at these scales. In marked contrast, the incremental slip rates of each of the four main Marlborough faults are highly variable through time, marked by coordinated accelerations and decelerations spanning 4–6 earthquakes and several millennia as the faults trade off slip to accommodate a steady relative plate motion rate. These results suggest that (a) the weakest fault in the system will slip faster than average while adjacent mechanically complementary faults slip more slowly, and (b) that these patterns switch back and forth through time, likely reflecting reversible changes in the strength (i.e., resistance to shear) of the individual faults as they collectively accommodate relative plate motion. Interestingly, the periods of fast slip on the MFS faults exhibit ∼20–25 m of displacement, suggesting that these may record periods of fast slip on a weakened fault/ductile shear zone that continues until it uses up all locally stored elastic strain energy, thus potentially approaching local complete stress drop, albeit during a few tens of meters of rapid fault slip during multiple earthquakes, rather than during a single event. This hypothesis is consistent with typical earthquake stress drops of ∼1–10 MPa and estimates of depth-averaged crustal shear stress of a few 10 s of MPa, such as might be released in clusters of 4–6 earthquakes. These results emphasize the need to analyze the collective behavior of the entire fault systems, rather than just individual faults, to understand the mechanics of the system. Moreover, these patterns suggest a potential path forward for more accurate estimation of time-dependent seismic hazard, with the possible incorporation of current position of a fault within a fast- or slow/no-slip period into the probability analysis, as well as a means of potentially estimating crustal shear stress