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.
As faults accumulate greater amounts of displacement, strain progressively localizes into a relatively narrow, structurally simple zone (e.g., Wesnousky, 1988; Chester et al., 2004; Dolan and Haravitch, 2014; Hatem, 2014). This process, known as structural maturation, decreases the proportion of total shear accommodated as off-fault deformation (OFD), which includes secondary faulting, warping, rotation, and distributed granular flow, that occur outside the fault core (e.g., Dolan and Haravitch, 2014). Additionally, the width and expression of OFD in geologic and geomorphic features will increase as deformation accumulates over multiple earthquakes, especially along structurally immature faults. In this way, the style and structural complexity of surface deformation recorded in recent deposits along a fault will reflect both the structural maturity of the underlying fault as well as the cumulative deformation resulting from earthquakes experienced by the deposit. Understanding how these processes become manifest in the landscape is essential to the proper geomorphic interpretation of surface faulting patterns.
We examine the relationships between structural maturity, strain localization, and the geomorphic expression of OFD by characterizing surface faulting patterns along two strike-slip faults with vastly different cumulative displacements. Specifically, we combine field observations with aerial lidar topographic data to compare surface faulting patterns in flights of fluvial terraces at the Branch River and Saxton River sites along the Wairau and Awatere faults in the Marlborough fault system of South Island, New Zealand (Fig. 1). These observations provide key insights into the interplay between overall fault structural maturity and the progressive geomorphic manifestation of OFD with increasing fault offset in young deposits.
COMPARISON OF SURFACE DEFORMATION PATTERNS ALONG FAULTS WITH DIFFERENT STRUCTURAL MATURITIES
The Wairau and Awatere faults have similar slip rates, lithologies, and tectonic settings (see the GSA Data Repository1; e.g., Wallace et al., 2012, and references therein), but have accommodated considerably different amounts of slip. Whereas the Wairau fault has accommodated much of the total ∼460 km of right-lateral offset along the Alpine fault system, the Awatere fault has accommodated only ∼13–20 km of strike slip (Fig. 1; Fig. DR1 in the Data Repository; Sutherland, 1999; Rattenbury et al., 2006; this study). As a first-order observation, the Wairau fault is geomorphically and structurally simpler over its length than the Awatere fault (e.g., Fig. DR2; Rattenbury et al., 2006). Fault surface expressions are, however, subject to numerous variables, including variations in the age, composition, fabric, and thickness of the near-surface geologic material through which the fault propagates (e.g., Zinke et al., 2014b; Teran et al., 2015). The Branch River and Saxton River sites are ideal for comparing surface deformation patterns on these faults, as both consist of flights of fluvial terraces that are similar in depositional age, geologic material, grain size, fault slip experienced by the deposits, and climate. We adopt the names and interpretations of the terraces used by previous workers at each site (Fig. 2; Lensen, 1968; Mason et al., 2006). Seven terraces (A–F and W in Fig. 2A) compose the flight at Branch River, whereas the Saxton River site contains six terraces (T1–T6 in Fig. 2B) and a bedrock promontory (“bedrock spur” of Mason et al., 2006). The terrace deposits at both of these sites consist of greywacke-dominated cobble gravels. Luminescence dating and cobble weathering-rind analyses indicate that the terrace ages at Branch River (7–18 ka) overlap with those at Saxton River (0–15 ka) (Knuepfer, 1992; Mason et al., 2006). Moreover, the oldest terraces at each site have experienced similar amounts of displacement (53 ± 2 m at Branch River; 69 ± 5 m at Saxton River [see the Data Repository; Lensen, 1968; Mason et al., 2006; Grenader et al., 2014; Zinke et al., 2014a]).
We used high-resolution (≥12 shots m–2) aerial lidar topographic data (see details in the Data Repository) to map fault-related deformation expressed in the micro-topography (Fig. 2; Fig. DR3). Field observations of surface deformation at these sites support our lidar mapping. At both sites, subtle (≥20 cm) intra-terrace channels and geomorphically expressed minor secondary faults are well preserved in the landscape. Throughout both study sites, secondary faulting expressed by vertical relief on this order is the only form of OFD discernible in the lidar data, and in the following we use the distribution of secondary faults as a proxy for other, geomorphically undetectable forms of OFD (e.g., warping, distributed granular flow).
The trace of the Wairau fault at Branch River (Fig. 2A) is straight and continuous over most of the 1.2-km-long study area, though we note the presence of a sag basin and a pop-up ridge near the eastern edge of the study area, each expressed by 2–3 m of topographic relief in terrace A. No secondary fault strands are geomorphically evident in the youngest four terraces (F–C). Two secondary fault strands, which are subparallel to the main trace, are discernible in the older W and A terraces. The fault-perpendicular width of geomorphically evident deformation at this site is <100 m at its widest and is effectively discrete (<3 m wide) along ∼80% of the length of the study site in all but the oldest, most-displaced terraces.
In contrast, the surface expression of the Awatere fault at Saxton River (Fig. 2B) is much more structurally complex. The main fault trace is multi-stranded throughout most of its length, and the main fault strands exhibit curvature on the half-kilometer scale. A 160-m-long, 2-m-deep sag basin is evident in the oldest T1 surface. The most striking difference in fault trace complexity between the two study sites is the abundance of secondary fault strands along the Awatere fault. These secondary faults are 14–342 m in length and are expressed by 0.2- to 1.1-m-tall scarps. Secondary fault strands are discernible to a maximum fault-perpendicular width of ∼460 m in the oldest terrace, T1.
The Branch River and Saxton River sites are generally representative of surface deformation patterns along the Wairau and Awatere faults, as similar trends in structural complexity are observed along the entire 40 km and 80 km fault lengths covered by our lidar data (Fig. DR2). Although sections of relatively structurally simple surface faulting occur at some locations along the Awatere fault, such examples are not representative of the overall structural complexity of this structurally less-mature fault. We also note that the wider zone of OFD at Saxton River is not a function of thicker sediments at that site, as the terrace gravels are only a few meters thick there (Mason, 2004), whereas they are ≥10 m thick at the structurally simpler Branch River site (observed in nearby cut banks).
We suggest that the dearth of OFD expressed as secondary faulting at the Branch River site, relative to that at the Saxton River site, results from the greater cumulative offset and consequent higher structural maturity of the Wairau fault. This inference is consistent with previous theoretical and observational studies that demonstrate that whereas faulting may be complex and relatively diffuse along lower-displacement, structurally immature faults, strain progressively localizes into a relatively discrete, structurally simple zone with increasing amounts of displacement (e.g., Sibson, 1985; Wesnousky, 1988; Stirling et al., 1996; Frost et al., 2009; Dolan and Haravitch, 2014).
Changes in fault structure will lead to the co-evolution of many important fault behaviors. Structural complexities such as steps and bends exert strong controls on fault rupture arrest, so structurally simpler faults will tend to rupture in larger events (e.g., Wesnousky, 1988, 2006; Stirling et al., 1996; Elliott et al., 2009). Additionally, the frequency content of seismic energy released in earthquakes likely depends on fault complexity (Dolan, 2006). Furthermore, fault structural maturity exerts a primary control on percentages of OFD in large (MW > 7) earthquakes (Dolan and Haravitch, 2014). Percentages of OFD can also affect measurements of slip in past earthquakes, with implications for the proper interpretation of fault slip rates based on offset geomorphic features and for the use of such data in seismic hazard assessments (Dolan and Haravitch, 2014).
PROGRESSIVE GEOMORPHIC MANIFESTATION OF OFD WITH INCREASING FAULT SLIP
The distribution of OFD expressed as secondary fault strands is strongly heterogeneous across the Saxton River site (Fig. 2B). In the oldest T1 surface, secondary faults discernible in the lidar define a zone of deformation across a total fault-perpendicular width of ∼460 m. Secondary faults evident in the younger T2 surface define a narrower deformation zone with a fault-perpendicular width of ∼100 m. OFD is also evident in the bedrock spur at fault-perpendicular distances of ∼220 m north of the main fault trace (the same distance as OFD in the adjacent T1 terrace) and ∼30 m south of the main fault trace, though some evidence for OFD south of the main fault trace may be obscured by colluvium or destroyed by slumping. Regardless, the total width of OFD within the bedrock spur evident in the lidar is ≥230 m. In contrast to these wide zones of secondary faulting in the older deposits, no geomorphic evidence of OFD is visible in either the lidar data or in our field mapping in the T3–T6 surfaces. In these younger terraces, the fault is expressed as a single, discrete, through-going feature.
Given that flights of fluvial terraces such as those at Saxton River represent time-transgressive sequences of deposits, the observations above indicate that OFD becomes more geomorphically discernible in progressively older terraces. The lack of topographically expressed OFD in the T6–T3 terraces is not likely to reflect changes in the structural maturity of the bedrock Awatere fault beneath the terrace deposits that occurred during the additional ∼19 m of displacement accommodated between formation of the T2 and T3 terraces, as tens of kilometers of fault slip is necessary to cause such an increase in structural maturity (e.g., Wesnousky, 1988; Stirling et al., 1996; Dolan and Haravitch, 2014). Rather, the processes that accommodate OFD are equally active during each earthquake recorded at the Saxton River site. Variations in the thicknesses of the gravel deposits composing each terrace could influence the distribution of fault-related deformation, such that thicker gravel deposits would cause more OFD (e.g., Van Dissen et al., 2011; Zinke et al., 2014b; Teran et al., 2015). Paleoseismic trenching revealed bedrock as shallow as ∼1.2 m below the T1 surface (Mason, 2004), whereas a road cut along the T2–T6 terrace riser exposes bedrock cropping out at 2–3 m below the T2 surface, suggesting that the T1 terrace gravels are similar in thickness to, or somewhat thinner than, the deposits composing the younger T2–T6 terraces. The fact that we observe more OFD in the T1 terrace than in the slightly thicker T2–T6 terraces reinforces our argument that cumulative slip, not sediment thickness, is the primary control on the manifestation of OFD in the landscape at Saxton River. Thus, whereas OFD may not be immediately evident in the landscape, it will become geomorphically better expressed as a function of increasing displacement (Fig. 3). In gravel-dominated deposits such as those at Saxton River, along faults of similar structural maturity to the Awatere fault, 35–50 m of displacement may be necessary for OFD to become geomorphically evident.
We observe a similar though much less-pronounced relationship between OFD and terrace displacement on the Wairau fault at Branch River. Whereas no OFD is observable in the lidar along the five youngest, least-offset terraces (F–B), we observe minor secondary faulting, out to ∼90 m fault-perpendicular width, in the two oldest terraces (W and A) which have been displaced by ≥∼55 m (Fig. 2). We reemphasize, however, that the dearth of geomorphically expressed OFD at the Branch River site relative to that at the Saxton River site in similarly offset terraces supports our inference that strain is much more localized along the structurally mature Wairau fault than along the much less-mature Awatere fault.
Analysis of surface deformation patterns associated with recent large-magnitude earthquakes using remote sensing techniques and field studies of offset linear cultural features reveals that surface ruptures are commonly surrounded by zones of OFD that extend as much as several hundred meters perpendicular to the fault (e.g., Rockwell et al., 2002; Van Dissen et al., 2013; Zinke et al., 2014b; Milliner et al., 2015). As much as ∼50% of total coseismic surface displacement can be accommodated as OFD along structurally immature faults (Dolan and Haravitch, 2014), and numerous geologic studies report evidence for off-fault coseismic strain (e.g., cracking, en echelon shearing) associated with recent surface ruptures (e.g., Zinke et al., 2014b; Teran et al., 2015). Yet quantifiable evidence of OFD is seldom manifest in the surrounding geomorphology, as many of the structures that accommodate OFD during surface ruptures are not likely to be preserved in the landscape. Thus, despite the importance of OFD in accommodating coseismic surface strain, well-defined geomorphic evidence for these zones that is likely to remain discernible for centuries to millennia is commonly lacking.
Geomorphic measurements of overall fault zone width may thus underestimate the true width of the fault zone, especially in deposits that have experienced small displacements along structurally immature faults. The ability to predict the width of fault-related deformation in advance of future events is important for understanding fault rupture mechanics and mitigation of seismic hazards to the built environment. For instance, improved constraints on expected surface deformation patterns could aid in the potential development of earthquake hazard microzonation maps detailing the expected magnitude and width of OFD, which in turn could be incorporated into standards for more effective earthquake-resistant design. Furthermore, failure to account for OFD in geomorphic measurements of fault slip may lead to systematic underestimation of the total fault slip, and consequently to underestimation of paleoearthquake magnitudes and fault slip rates. Because fault slip rates are a basic input into most probabilistic seismic hazard analyses, this will lead to underestimation of the seismic hazard associated with a fault. Proper interpretation of fault-related deformation patterns in the landscape therefore rests on the recognition that although OFD may not be fully geomorphically expressed, the processes that accommodate OFD nevertheless play a crucial role in accommodating coseismic strain in the near surface, especially along structurally immature faults.
This research was funded by the U.S. National Science Foundation (grant EAR-1321914 to Dolan). Lidar data were collected by the U.S. National Center for Airborne Laser Mapping (NCALM) and New Zealand Aerial Mapping. Thanks to Michael Sartori, Juan Fernandez, and Heezin Lee for help with lidar collection and processing. Special thanks to Simon and Jane Fowler, Muller Station (and the Satterthwaites), and Molesworth Station (and the Wards) for graciously allowing us access. This study benefited from reviews by Kate Scharer, Olaf Zielke, and Nic Barth.