Fault zone architecture of a large plate-bounding strike-slip fault: a case study from the Alpine Fault, New Zealand

New Zealand’s Alpine Fault is a large, platebounding strike-slip fault, which ruptures in large (Mw > 8) earthquakes. We conducted field and laboratory analyses of fault rocks to assess its fault zone architecture. Results reveal that the Alpine Fault Zone has a complex geometry, comprising an anastomosing network of multiple slip planes that have accommodated different amounts of displacement. This contrasts with the previous perception of the Alpine Fault Zone, which assumes a single principal slip zone accommodated all displacement. This interpretation is supported by results of drilling projects and geophysical investigations. Furthermore, observations presented here show that the young, largely unconsolidated sediments that constitute the footwall at shallow depths have a significant influence on fault gouge rheological properties and structure.


Introduction
A fault, a planar discontinuity in rock where one side has moved relative to the other parallel to the discontinuity plane, constitutes a rheological and mechanical manifestation of localized deformation (Twiss and Moores, 2007;Ben-Zion, 2008;Fossen, 2016;Fossen and Cavalcante, 2017). Large fault zones are typically composed of networks of smaller, individual, but related and interacting, faults of self-similar geometry (Ben-Zion and Sammis, 2003;Twiss and Moores, 2007;Peacock et al., 2016). The structure, composition, hydrological properties and seismo-mechanical behavior of faults are typically intimately related. These interactions also govern strain distribution and depend on various factors, such as lithology (Faulkner et al., 2003;Schleicher et al., 2010;Holdsworth et al., 2011;Rybacki et al., 2011), fluid pressure (Hickman et al., 1995;Janssen et al., 1998;Fagereng et al., 2010), stress field and stress magnitudes (Sibson, 1985;Faulkner et al., 2006;Lindsey et al., 2014).
Faults control the strength of the Earth's lithosphere (Townend and Zoback, 2000;Bürgman and Dresen, 2008) and govern substantially fluid flow (Caine et al., 1996;Wibberley et al., 2008). Hydrocarbon production from faultcompartmentalized reservoirs (Van Eijs et al., 2006), exploitation of fault-hosted mineral deposits (Cox et al., 1986) and long-term integrity of potential nuclear waste repositories (Laurich et al., 2018) are practical examples demonstrating how important it is to understand fault zone properties and their spatial and temporal evolution. Furthermore, large, plate-bounding strike-slip faults such as the Alpine Fault (New Zealand), the North Anatolian Fault Zone (Turkey) and the San Andreas Fault (USA) rupture in large (M w > 7) earthquakes (Toppozada et al., 2002;Sutherland et al., 2007;Bohnhoff et al., 2016). Many of these faults are located in densely populated areas so they pose a significant geohazard (Eguchi et al., 1998;Sahin and Tari, 2000;Martínez-Garzón et al., 2015). Thus, from a georesource, seismic hazard and  Caine et al. (1996) a fault is a relatively simple structure, where most of the strain is accommodated at a single, quite narrow fault core hosting a principal slip zone. (b) Faults described by the model of Faulkner et al. (2003) are more complex and consist of a damage zone hosting multiple, anastomosing principal slip zones forming a complex network. risk perspective, it is important to characterize the seismomechanical properties of faults (Aksu et al., 2000;Zoback et al., 2007;Hollingsworth et al., 2017). Caine et al. (1996) presented a conceptual model of the structure of fault zones that has three primary components (Fig. 1a). The protolith (I) hosts a damage zone (II) that is characterized by a fracture density significantly higher than the background values of the surrounding host rock (Chester et al., 1993;Biegel and Sammis, 2004;Faulkner et al., 2010). The damage zone contains the fault core (III) where most of the displacement has been accommodated (Caine et al., 1996). This conceptual framework has been found to be applicable to faults across the full range of natural and experimental spatial scales (Anders and Wiltschko, 1994). Biegel and Sammis (2004) suggested that a fault core should be defined as the zone within a fault where strain leads to granulation of rocks, distinct from a damage zone within which fracture density is high but fracturing has "not [been] sufficient to produce distinct particles". A fault core can be structurally and lithologically heterogeneous, and most of the displacement it accommodates may be localized along one or more principal slip zones (PSZs) defined by ultracataclasites or fault gouges (Sibson, 2003;Janssen et al., 2014;Toy et al., 2015). Although most commonly less than 10 cm thick, PSZs may be up to 1 m thick (Sibson, 2003) and thicknesses tend to increase with increasing displacement (Evans, 1990;Faulkner et al., 2010;Ben-Zion and Sammis, 2003, and references therein). In general, fault zone thicknesses may range from millimeters (e.g., Antonellini and Aydin, 1994) to a few hundred meters (e.g., Bruhn et al., 1994).
The Punchbowl Fault, an inactive, exhumed part of the San Andreas Fault system, typifies this conceptual model of fault zone architecture. It has a single PSZ embedded in a fault core and the surrounding damage zone. A single, continuous gouge layer with 30 cm thickness on average hosts the ∼ 1 mm thick PSZ. The core of the Punchbowl Fault is surrounded by an approximately 15 m thick damage zone (Chester and Logan, 1986;Chester et al., 2005).
However, other faults show a more complex structure with changing properties along strike or down dip (e.g., Wibberley et al., 2008;Faulkner et al., 2010). For example, detailed studies of the Carboneras Fault, Spain, yielded a conceptual model that is suitable for broader, typically phyllosilicaterich fault zones, which tend to contain multiple high-strain zones (Fig. 1b;Faulkner et al., 2003). The Carboneras Fault is a predominantly strike-slip structure that comprises an ∼ 1 km thick tabular zone of continuous and anastomosing strands of phyllosilicate-rich fault gouge containing lenses of fractured protolith surrounded by an approximately 100 m wide damage zone. Scaly clays, which typically contain anastomosing shear planes, are examples of this distributed fault zone model on the hand specimen scale (Vannucchi et al., 2003;Laurich et al., 2017).
The fault zone architecture of New Zealand's Alpine Fault, a large, transpressional plate-bounding fault and a significant geohazard, has attracted increasing attention in the last 10 years. Considered to be appropriately described by the conceptual model of Caine et al. (1996), a single PSZ is commonly supposed to accommodate displacement along the Alpine Fault (e.g., Barth et al., 2012;Sutherland et al., 2012;Barth et al., 2013;Norris and Toy, 2014;Toy et al., 2015;Williams et al., 2016;Townend et al., 2017;Schuck et al., 2018;Williams et al., 2018). Here, by combining results of these previous studies on the Alpine Fault's structure with new field observations, microstructural, mineralogical and geochemical analyses, we show that the Alpine Fault has a complex fault geometry manifested by PSZ thicknesses differing substantially between investigated outcrops. This complexity is not controlled by lithology, which is unexpected given that the Alpine Fault is hosted in a quartzofeldspathic protolith, which commonly fosters localization compared to phyllosilicate host rocks (Faulkner et al., 2003). This implies that the Alpine Fault does not fit our paradigmatic models of fault zone architecture.
from Milford Sound in the SW to Hokitika in the NE, where it transfers displacement onto the four main fault strands of the Marlborough Fault System (the Wairau, Awatere, Clarence and Hope Faults; Fig. 2a). The Alpine Fault maintains a constant NE-SW strike for its entire length, but the fault dip changes from 80-90 • SE SW of the Dun Mountain Ophiolite Belt (DMOB) to 30-50 • SE in the central segment (Barth et al., 2013). The shallowest 1-2 km of the fault NE of Haast display "serial partitioning", i.e., northerly striking oblique thrust sections alternate with easterly striking dextral strike slip sections of 1-10 km length (Norris and Cooper, 1995;Barth et al., 2012).
A 470 km right-lateral offset of the DMOB defines the minimum cumulative displacement along the fault (Sutherland et al., 2007). Displacement-normal shortening is on the order of 90 ± 20 km (Little et al., 2005). Strike-slip displacement rates are between 21 and 29 ± 6 mm yr −1 (Norris and Cooper, 2000) corresponding to 60 %-80 % of the total relative velocity between the bounding plates (DeMets et al., 2010). Seismic investigations indicate a maximum exhumation of ∼ 35 km from a deep crustal, subhorizontal, NE-dipping detachment Little et al., 2005). Long-term exhumation rates of 6-9 mm yr −1 are inferred from 40 Ar/ 39 Ar ages encountered between Franz Josef Glacier and Fox Glacier (Little et al., 2005). A combination of these high exhumation rates and meteoric fluid circulation driven by steep topographic gradients results in a very high geothermal gradient of up to 125 • C km −1 encountered in valleys (Menzies et al., 2016;Sutherland et al., 2017).
Currently, the Alpine Fault does not exhibit creep and is thought to be seismically locked to a depth of 12-18 km . It is known to rupture in large earthquakes (M w > 8), while generating up to 8-9 m of lateral and up to 1 m of normal displacement (Sutherland et al., 2007;Nicol et al., 2016). Offset and deformed Quaternary features demonstrate that Alpine Fault earthquakes rupture the ground surface (e.g., Cooper and Norris, 1990;Berryman et al., 2012;Schuck et al., 2018). While events like the most recent one in 1717 AD might rupture the entire fault, differing recurrence intervals of 263±68 years for the central segment and 291±23 years for the southern segment demonstrate that individual sections of the fault might fail independently ( Fig. 2a; Sutherland et al., 2007;Howarth et al., 2018 and references therein). Considering the potential to produce large-magnitude earthquakes and the time passed since the last event, the Alpine Fault is late in its seismic cycle and thus constitutes one of the South Island's major geohazards.

Lithology
The Australian Plate footwall assemblage encompasses Paleozoic to Cretaceous plutonic rocks intruded into metasediments (Fig. 2b). As these units are mostly overlain by Quaternary fluvioglacial sediments, footwall rocks are poorly ex-posed (e.g., Toy et al., 2015). In the Pacific Plate hanging wall, a narrow, 12-25 km wide elongate belt of metamorphosed sediments, the Alpine Schist, is exposed from the SW of Jackson Bay to the fault's NE termination (e.g., Sibson et al., 1979;Grapes and Watanabe, 1992;Scott et al., 2015). The Alpine Schist originates mostly from the Torlesse Terrane. This terrane is a polygenetic metamorphic suite with the Alpine Schist being its high grade and most recently formed (Late Cretaceous) part (Roser and Scott et al., 2015).
The amphibolite-greenschist facies rocks of the Alpine Schist have dominantly metapelitic to metapsammitic compositions, with rare metabasite and metachert. Metamorphic grade decreases from K-feldspar and oligoclase amphibolite facies through garnet, biotite and chlorite greenschist facies to pumpellyite-actinolite and prehnite-pumpellyite facies with increasing SE distance from the fault plane (e.g., Grapes and Watanabe, 1992;Scott et al., 2015). However, a combination of significant right-lateral displacement and high exhumation rates resulted in 200-300 m wide inverted metamorphic sequences cropping out structurally above brittle fault rocks in the central segment (Cooper and Norris, 2011). Irrespective of local variations, typical mineral phases encountered in Alpine Schist-derived fault rocks are quartz, feldspar, garnet, muscovite, biotite, other minor phases and, in the case of metabasites, hornblende and epidote . At the Waikukupa Thrust there are remnants of highly disrupted, intensely strained and sheared granite pegmatites indicating simple shear strains of more than 150 and coeval pure shear stretches of ∼ 3.5 (Norris and Cooper, 2003;Toy et al., 2013).
Whereas the Alpine Schist constitutes the protolith of fault rocks in the central portion of the Alpine Fault, a comparatively small (∼ 550 m thick) sliver of the Brook Street Terrane (BST) is the protolith of samples from the southern segment (Fig. 2b). The sliver is surrounded by the DMOB and bound to the NE by the Alpine Fault. The BST is the remnant of a Permian island-arc system, mainly composed of basaltic to andesitic volcaniclastic and sedimentary rocks of prehnite to pumpellyite and locally greenschist metamorphic grade (Spandler et al., 2005).

Fault rocks of the Alpine Fault
Exhumation from deep-crustal levels resulted in the formation of an approximately 1 km wide characteristic sequence of hanging wall fault rocks (e.g., Reed, 1964;Sibson et al., 1979;Cooper and Norris, 2011;Toy et al., 2015, Fig. 2c). Progressive northwestward-increasing strain of the Alpine Schist yielded protomylonites to ultramylonites that are overprinted by a cataclastic fault zone displaying increasing damage towards the PSZ (Williams et al., 2016). Cataclastic shears in the ultramylonites are filled with phyllosilicates and exhibit strike-slip, normal and reverse kinematics. Cataclasites are composed of comminuted mylonite fragments   (Edbrooke et al., 2015) and the digital elevation model of Columbus et al. (2011). (c) Typical shallow-depth sequence of Alpine Fault rocks. The figure was modified from Norris and Cooper (2007). Extents of alteration zone were derived from Sutherland et al. (2012) and extents of inner and outer damage zone from Townend et al. (2017). within a matrix of pulverized host rock, authigenic chlorite, muscovite and illite (Sibson et al., 1979;Norris and Cooper, 2007;. Intense chloritization imparted a typical pale green color to these 10-30 m wide cataclasites (Warr and Cox, 2001), which are cemented by authigenic phases, dominantly phyllosilicates and carbonates (Sutherland et al., 2012;Toy et al., 2015;Williams et al., 2016). The fault's PSZ cross-cuts all other structures, which indicates that it is the last active part of the fault system and accommodates the coseismic displacements that manifest as offset Quaternary features (e.g., Cooper and Norris, 1990;Berryman et al., 2012;Schuck et al., 2018). It is an incohesive and fine-grained gouge up to ∼ 50 cm in thickness . Authigenic phyllosilicates cement the PSZ converting it into an impermeable hydraulic barrier preventing cross-fault fluid flow (Sutherland et al., 2012;Menzies et al., 2016). There is no location providing a continuous section from hard rock hanging wall to hard rock footwall (Townend et al., 2009).
Another two boreholes were drilled in DFDP phase 2 (DFDP-2A: 212.6 m MD -measured depth; DFDP-2B: 893.1 m MD) about 7.5 km ENE of the location of DFDP-1 (Fig. 2b, Toy et al., 2017). Phase 2 yielded cuttings for petrographic investigations  and provided various insights into fault zone architecture by wireline logs (e.g., Massiot et al., 2018) but was not able to penetrate and sample the fault core due to scientific and technical difficulties . DFDP boreholes from both phases provide the opportunities for longterm monitoring of the Alpine Fault (e.g., Sutherland et al., 2015).

Study locations
For this study, all known outcrops with accessible PSZ between Martyr River in the SW and Kokatahi River, SSE of Hokitika in the NE were investigated in the austral summer 2015/2016 (Fig. 2b). Furthermore, cataclasite and the fault gouge samples were available from the DFDP-1A core. Four of the five investigated locations, Havelock Creek, Waikukupa Thrust, Gaunt Creek and the core of DFDP-1A, are in the central segment of the Alpine Fault with the Alpine Schist as hanging wall host rock. Rocks of the BST constitute the host of fault rocks encountered at Martyr River, the fifth and southwestern-most location. All outcrops are located along river banks providing good outcrop conditions. Waikukupa Thrust is the only studied site that is considered to be inactive. Incision of the Waikukupa River resulted in a geomechanically unfavorable geometry leading to the abandonment of this thrust segment and re-localization of deformation approximately 700 m to the NE at Hare Mare Creek (Norris and Cooper, 1997).

Sampling
Fault rocks at Martyr River and Havelock Creek were sampled using hammer and chisel. At Waikukupa Thrust and Gaunt Creek samples were recovered using a hand-held chainsaw equipped with a silicon carbide chain. After recovery, all outcrop samples were immediately wrapped in aluminum and plastic foil to slow down drying. This procedure allowed for continuous transects across the PSZ, which juxtapose hanging wall cataclasites and footwall gravels, to be sampled from Martyr River and Waikukupa Thrust. However, the samples broke apart along the PSZ during shipping. An intact contact of PSZ on footwalls material was only preserved from Gaunt Creek.
In summary, hanging wall fault rocks are available from Martyr River, Waikukupa Thrust and DFDP-1A. PSZ samples originate from all five locations and footwall rocks were taken at Martyr River, Waikukupa Thrust and Gaunt Creek.

Samples from Martyr River, Waikukupa Thrust and Gaunt
Creek were cut dry with a low-speed saw perpendicular and parallel to the fault trace. From all locations, subsamples were selected for microstructural analyses and subsequently embedded in resin prior to the preparation of 29 dry-polished thin sections using otherwise standard techniques.
Microstructural investigations were conducted using optical microscopy, scanning electron microscopy (SEM), transmission electron microscopy (TEM) and cathodoluminescence analyses (CL). Thin sections were studied with a Leica DM RX optical microscope and a FEI Quanta 3D SEM with focused ion beam (FIB; dual-beam machine). The SEM, equipped with a field emission gun, operated at 20 kV in backscatter electron mode (BSE) and allowed semiquantitative geochemical analyses with its EDAX energy dispersive X-ray analyzer (EDX). For TEM analyses, a platinum strip was deposited on the sites selected with the FIB of the SEM to enable subsequent preparation of 14 thin foils (10 × 8 × 0.15 µm) with a FEI FIB 200 following the procedure outlined by Wirth (2004Wirth ( , 2009. A FEI Tecnai G2 F20 X-Twin TEM with Gatan Tridiem energy filter, Fischione high-angle annular dark field detector (HAADF) and EDX operated at 200 kV and allowed nanoscale investigation. Cathodoluminescence analyses on calcite veins and cement were performed on thin sections using an Olympus polarizing microscope. The Lumic HC1-LM hot cathode CL microscope operated at 14 kV, 0.0001 mbar and ∼ 0.6 mA electron beam and 2.5 A filament current, respectively.
Microstructural analyses are mostly based on qualitative observations, but some quantitative information was obtained by image analysis using the open source software FIJI/ImageJ (Schindelin et al., 2012;Schneider et al., 2012).

Mineralogical and geochemical investigations
18 bulk rock powder samples for mineralogical and geochemical analyses were prepared with a jaw crusher and subsequently sieved to grain sizes < 62 µm. A McCrone micronization mill provided the < 10 µm grain size fraction for the quantitative mineralogical analysis.

X-ray diffraction analysis
The mineralogical compositions of the samples were determined by X-ray diffraction (XRD) analysis on random powder samples loaded from the back side of the sample holders (∅ 16 mm). XRD analyses were performed with a PANalytical Empyrean X-ray diffractometer operating with Bragg-Brentano geometry at 40 mA and 40 kV with CuKα radiation and a PIXel3D detector at a step size of 0.013 • 2 from 4.6 to 85 • 2 and 60 s per step.
Mineralogy was determined qualitatively with the EVA software (version 11.0.0.3) by Bruker. Rietveld refinement for quantitative mineralogy was performed using the program BGMN and the graphical user interface Profex (version 3.10.2, Döbelin and Kleeberg, 2015) calibrated for the used diffractometer. The error of quantitative analyses is expected to be in the range of 3 wt %. Goodness of fit was assessed by visually inspecting the differences between measured and modeled diffractograms and by aiming to obtain R wp values lower than 3, a value arbitrarily chosen (Toby, 2006).
Bulk powder subsamples (45 mg), dispersed in 1.5 mL of de-ionized water and disaggregated in an ultrasonic bath, were placed on round (∅ 3.1 cm) glass slides and air-dried for subsequent smectite and chlorite analyses. Identification of smectite was performed after ethylene glycolation at 45 • C for at least 12 h (2.4-25 • 2 ), and heating to 500 • C for 1 h allowed us to differentiate between chlorite and kaolinite (2.4-25 • 2 , Moore and Reynolds Jr., 1997).
Illite polytype analyses to differentiate between authigenic and detrital (i.e., comminuted muscovite) species were conducted on selected samples with high illite and/or muscovite concentrations (Moore and Reynolds Jr., 1997;Haines and van der Pluijm, 2008). Diffractograms for polytype analyses were acquired between 16 and 44 • 2 at a step size of 0.013 • 2 with 200 s per step. All diffractograms and R wp values are provided in the Supplement.
3.3.2 X-ray fluorescence spectrometry X-ray fluorescence spectrometry (XRF) with a PANalytical AXIOS Advanced spectrometer enabled study of geochemical variations of major elements. For this purpose, after drying of ∼ 1.1 g of bulk rock powder overnight at 105 • C, 1 g of sample material, 6 g of LiBO 2 and 0.5-0.7 g of ammonium nitrate, acting as oxidizing agent, were fusion digested at increasing temperatures from 400 to 1150 • C with a dilution of 1 : 6. The loss on ignition (LOI) providing volatile contents was determined on 20-30 mg of powder samples with an Euro EA Elemental Analyzer.
To evaluate fluid-related element mobilization, XRF data were analyzed based on the equation of Gresens (1967) for composition-volume relationships resulting from metasomatic alterations, employing the isocon method (Grant, 1986(Grant, , 2005. In this, the chemical composition of each investigated fault rock sample (i.e., altered rock) is plotted against the chemical composition of the host rock (i.e., unaltered rock). An isocon, a straight line through the origin, separates species enriched relative to the host rock plotting above from those depleted plotting below. It follows the equation with C A and C H being element concentrations of the altered and the host rock, respectively, and m referring to the isocon's slope. Enrichment or depletion ( C i ) of an element i relative to its host rock equivalent is calculated as follows: The inverse of the isocon's slope (m −1 ) indicates overall mass gain (m −1 > 1) or loss (m −1 < 1), respectively. Mass gains (%) are calculated as (m −1 − 1) × 100 %. A negative mass gain corresponds to a mass loss (%).
To derive isocon slopes, Grant (2005) suggested five approaches: (I) clustering of C A i /C H i , (II) a best fit of data forming a linear array through the origin on an isocon diagram as graphical equivalent to (I), (III) assuming that certain elements are immobile, or assuming (IV) mass or (V) volume to be constant during alteration. However, it has previously been shown that Al and Ti -elements commonly assumed to be not affected by fluid-related alteration (e.g., Grant, 2005;Dolejš and Manning, 2010) -are mobilized by fluid-related alteration at the Alpine Fault. Furthermore, authigenic mineral formation is thought to affect both fault rock mass and rock volume (Schuck et al., 2018). As identification of clusters (I) or linear arrays (II) selects only a subset of all data points, and these selections may be biased, we used a linear best fit through all data points to determine isocon slopes.
Identification of the host rock is just as important as correctly determining the slopes of the isocons (Grant, 1986). We defined a reference protolith composition for fault rocks from Havelock Creek, Waikukupa Thrust, Gaunt Creek and DFDP-1A (Table S1 in the Supplement) based on the average of geochemical compositions of Alpine Schist data reported by Grapes et al. (1982), Roser and  and Toy et al. (2017). These analyses are derived from the entire central segment of the Alpine Fault and individual data points display little variation.
At Martyr River, the geochemical compositions by Spandler et al. (2005, see their Table 1) have been averaged, and used as host rock composition (Table S1 in the Supplement). Their analyses from Bluff Complex and Takitimu Mountains were excluded, because the origin and some geochemical features, respectively, are different from the bulk BST (Spandler et al., 2005). We also note the isocon method has limitations for understanding fluid-assisted alteration related to fault activity at Martyr River. This is because the width of the Alpine Fault deformation zone at this location is larger than the thickness of the BST (Fig. 2b), so faulting-related fluid-rock interaction likely occurred with the DMOB and the geochemical composition of the reference protolith was already affected by fluid-related alteration. Furthermore, the BST is compositionally very heterogeneous due to magmatic differentiation so it is difficult to choose a single representative host rock geochemical composition. Fortunately, these variations only substantially affect the absolute quantification. The method nevertheless allows for relative geochemical changes within the fault zone at this location to be accounted for. The limitations regarding absolute quantification of concentration changes as result of hydrothermal alteration also apply to the Alpine Schist.

Field observations
Individual units of the Alpine Fault rock sequence -(ultra)mylonites, cataclasites and footwall gravels -can be identified clearly at all sites investigated. However, weathering and landslides complicate the unambiguous identification of the PSZ at Gaunt Creek and Havelock Creek. This is exemplified by an approximately 1 cm thick, fine-grained, brown, fractured and discontinuous layer interpreted to represent the PSZ at Gaunt Creek ( Fig. 3c-e) that was identified through sampling to be a footwall feature approximately 8 cm below the PSZ (see Sect. 4.2.3). Furthermore, measurements of fault orientation are only robust in the well-exposed outcrops at Waikukupa Thrust (50 • /44 • SE) and Martyr River (59 • /51 • SE) and those from Gaunt Creek (61 • /32 • SE) and  Thicknesses of fault gouges constituting the PSZ vary among individual locations, including the DFDP-1A core , and decrease in the following order: Havelock Creek (∼ 50 cm; Figs. 4a and 5a), Gaunt Creek (∼ 30-40 cm; Figs. 3c-e and 5b), DFDP-1A (20 cm; Fig. 3b), Waikukupa Thrust (∼ 4 cm; Figs. 4b and 5c) and Martyr River (∼ 1-2 cm; Figs. 4c, d and 5d). Except for Martyr River, these locations have the same protolith ( Fig. 2b; see Sect. 4.3). Furthermore, there is no systematic variation of PSZ thickness in distance along strike in these outcrops.

Microstructures
The fault rocks of Alpine Fault hanging wall, PSZ and footwall contain clasts embedded in a fine-grained matrix. We define matrix as fine particles -detrital and authigenicbelonging to the clay-size fraction, i.e., < 2 µm (Fig. 6a). In fact, the typical matrix grain size is 1 µm. The matrix is composed of abundant phyllosilicates of detrital (see Sect. 4.3) and authigenic (e.g., Fig. 6a) origin and detrital particles, mostly quartz and feldspar. The latter are subangular to rounded in shape with elongate to moderate sphericity and commonly with edges that have been affected by dissolution and mineral alteration processes (Fig. 6a).
Clasts embedded in the fault rock matrix are predominantly comminuted quartz, feldspar and mica grains, as well as fragments of Alpine Schist and mylonite. In addition, there are two more groups of clasts. The first one comprises compact, almost pore-free, fine particles with distinct boundaries but of similar composition compared to the surrounding matrix (Fig. 6b). This group of clasts is termed matrix clasts in the following. The second group comprises small areas, which are microstructurally similar to matrix clasts. These features are characterized by brighter grey values in BSE mode, and typically high porosities so they are easy to differentiate from the surrounding material (Fig. 6c). EDX analyses demonstrate that they are compositionally identical to the surrounding matrix. Because comparison of optical and scanning electron microscopy images reveals that these structures look like matrix clasts under plain polarized light ( Fig. 6d and e), we refer to them as bright matrix clasts in the following.
The distinct lithological and (micro-to-macro) structural characteristics of all units are described in detail in the following sections and summarized in Tables 1 (hanging wall), 2 (PSZ) and 3 (footwall).

Hanging wall cataclasites
Clasts embedded in the fine-grained cataclasite matrix at DFDP-1A and Waikukupa Thrust are mostly quartz and feldspars, with commonly fractured, dissolved and altered edges (Fig. 6a). Furthermore, matrix clasts are encountered close to the contact with the PSZ (< 1.35 m at DFDP-1A and < 10 cm at Waikukupa Thrust, respectively). In contrast, quartz and feldspars only constitute ∼ 5 % of clasts at Martyr River. There, most clasts are polymineralic aggregates consisting of quartz, feldspar, chlorite and rarely biotite and fragments of calcite veins ( Fig. 6f and g). There are subordinately fragments of Alpine Schist (all locations) and mylonites (Waikukupa Thrust). At DFDP-1A and Waikukupa Thrust these fragments display ∼ 150 µm wide, gouge-filled microfaults (Table 1). At Waikukupa Thrust, there are also clay-clast aggregates (CCAs; Boutareaud et al., 2008). At both thin section and outcrop scale, the amount and size of clasts generally decrease towards the PSZ and both vary systematically with PSZ thickness (Table 1; see also Toy et al., 2015): locations with a thinner PSZ contain more clasts in the hanging wall, which tend to be larger, compared to locations with a thicker PSZ.
Open cracks, typically several millimeters long and up to hundreds of micrometer wide, increase in abundance towards the PSZ at DFDP-1A and Waikukupa Thrust. Furthermore, microfaults (∼ 5-130 µm thick) cross-cut cataclasites at these locations. They contain comminuted particles and authigenic phyllosilicates and range in maturity: mature microfaults display a well-developed fine-grained fault gouge with particles up to tens of micrometers in size but mostly smaller than 2 µm (Fig. 6h). In contrast, juvenile microfaults contain poorly sorted particles up to a few hundred micrometers in size with the majority larger than 10 µm. In addition, some microfaults have cores made of coarser-grained matrix and clasts (Fig. 6h) or calcite veins (Fig. 6i) surrounded by fine-grained gouge.
Locally, there are patches almost completely devoid of clastic particles but with abundant authigenic flake-to needle-shaped chlorite crystals cementing pores and fractures at Waikukupa Thrust and DFDP-1A ( Fig. 6j and k). The amount of these chlorite-dominated areas and authigenic phyllosilicates in general increases with decreasing distance to the PSZ and progressively cements the fault rocks. Furthermore, the matrix is cemented by calcite (Table 1). Whereas finely dispersed calcite is abundant at Martyr River, calcite at Waikukupa Thrust is usually present only in direct vicinity (< 3 cm) to the PSZ. There, mutually crosscutting calcite veins (∼ 15 µm thick) display small (< 20 µm) dextral offsets (see Fig. 3f in Schuck et al., 2018). Furthermore, there are < 100 µm large calcite-cemented breccias composed of angular fragments of clastic particles and matrix clasts (see Fig. 3c in Schuck et al., 2018). At DFDP-1A, calcite is encountered in lenses and veins, which frequently cross-cut clasts but rarely the matrix and generally increase in abundance towards the PSZ.
Locally, cataclasites are foliated. At DFDP-1A, clasts together with microfaults define a weak foliation (Fig. 6l) varying non-systematically across the investigated interval, locally displaying SC geometry. At Waikukupa Thrust and close to the PSZ, cataclastic lenses exhibit a very weak foliation parallel to the shear plane.

Principal slip zone fault gouges
The type of contact between hanging wall cataclasites and the fault gouge appears to correlate with PSZ thickness: where the PSZ is thicker, contacts are transitionally manifested by decreasing grain sizes (Table 2), as exemplified at Havelock Creek with its transitional, poorly developed contact over 8 cm (Fig. 5a). In contrast, Waikukupa Thrust has a sharp contact on the outcrop and hand specimen scale (Figs. 4b and 5c), but displays a transition over 4 mm from a porous, clast-rich cataclasite matrix to a very dense PSZ matrix as revealed by microscopic inspection (see Fig. 3a in Schuck et al., 2018). Identification of the sharp but undulating contact at Martyr River is simple on the outcrop and hand specimen scale (Figs. 4c, d and 5d), but difficult on the microscale ( Fig. 7a and b).
Matrix clasts are the prevailing type of clast at Waikukupa Thrust and Havelock Creek ( Fig. 7c and d), but constitute less than 5 % of clasts at Martyr River. Furthermore, the phyllosilicate-cemented fault gouge at Martyr River hosts individual mineral phases, which is distinctly different from its hanging wall, where large, fine-grained chloritized particles dominate (Fig. 7b). Bright matrix clasts are present at all locations, but are abundant only at Havelock Creek and Gaunt Creek, where they, in addition to clast-like shapes, occupy areas of several hundred micrometers. In general, the amount and size of clasts weakly correlate with PSZ thickness (Table 2). Where the PSZ is thinner, clasts tend to be more abundant and sizes tend to be larger. Clast sizes increase in the following order: Gaunt Creek (∼ 75 µm), Havelock Creek (∼ 75 µm), Waikukupa Thrust (∼ 100 µm) and Martyr River (tens of micrometers to ∼ 200 µm).
It is not possible to determine the amount and size of clasts within the DFDP-1A core, because more than 90 % of this unit is cemented by authigenic phyllosilicates, mostly chlorite ( Fig. 7e and f). Authigenic phyllosilicates tend to nucleate along crack and grain boundaries of fractured grains. Additionally, newly formed phyllosilicates replace dissolved larger particles and mimic their original shape. Non-cemented areas are restricted to fractures. Many of these fractures contain up to 350 µm wide calcite cores, locally surrounded by gouge, microstructurally similar to microfaults observed within the DFDP-1A cataclasites (see Fig. 6i). These non-cemented, fracture-related locations host patches of fine-grained, randomly oriented and needle-shaped, authigenic chlorite crystallites. The presence of these patches within the PSZ is restricted to DFDP-1A.   We only saw evidence of pressure solution at one location, DFDP-1A, where there is a slickolite (Fig. 7g), a stylolite with teeth oblique to the stylolite surface (see Passchier and Trouw, 2005). However, enrichment of insoluble material or secondary phases along teeth-crowns is not observed.
PSZs at Waikukupa Thrust and Martyr River are layered and display a hanging wall-proximal and a footwall-proximal layer, respectively ( Fig. 5c and d). At Martyr River, the contact is best recognized on the hand specimen scale (Figs. 4d,  5d, 7a and b), where the fault gouge is blueish-grey proximal to the hanging wall and greenish-grey proximal to the footwall. Minor differences between both gouge layers are manifested by slightly more and larger clasts in the ∼ 0.8-1 cm thick hanging wall-proximal layer than the ∼ 0.3-0.5 cm thick footwall-proximal layer (20 %-25 % vs. 25 %-30 %; Table 2).
The dense, clast-poor (< 1 % clasts) fault gouge at Gaunt Creek is microstructurally identical to the footwall-proximal layer at Waikukupa Thrust ( Fig. 7h and i). Furthermore, both gouges exhibit calcite vein networks close to the contact with the footwall (Fig. 7j-m). At Waikukupa Thrust, veins display mutual cross-cutting relationships with dextral offsets. For a detailed description of these calcite veins, the reader is referred to the preceding work of Schuck et al. (2018). At Gaunt Creek, the 350-800 µm wide vein network ( Fig. 7km) has been subsequently affected by intensive fracturing and calcite dissolution, which hampers the identification of cross-cutting relationships. Furthermore, the contact between calcite veins and the surrounding gouge is poorly developed. Veins are typically 5-10 µm wide and tend to be thicker, less dissolved and deformed close to the footwall. CL colors are yellow-orange but slightly brighter and more yellow towards the footwall (Fig. 7l).
Apart from these vein networks, calcite veins crosscut matrix clasts at Martyr River and the hanging wallproximal gouge layer of Waikukupa Thrust, but are very rare at Havelock Creek. Additional hanging wall-proximal microstructures are clasts forming a weak foliation parallel to displacement at Martyr River and CCAs at Waikukupa Thrust. Furthermore, hanging wall-proximal microstructures at Waikukupa Thrust are similar to those of Havelock Creek  Fig. 5b). The white box indicates location of (e). (e) SEM analysis shows that clasts within a lens ∼ 8-8.5 cm below the PSZ at Gaunt Creek are bright matrix clasts. (f, g) Corresponding PPL (f) and XPL (g) photomicrographs of hanging wall clasts at Martyr River. Identification of constituting mineral phases -mostly quartz, feldspar, chlorite -is only possible by EDX analysis. (h) Mature microfault filled by fine-grained authigenic phyllosilicates with typical grain sizes < 2 µm and a core of coarser-grained material. (i) Microfault displaying a calcite vein within its core surrounded by a fine-grained matrix and some larger clasts. (j) Locally, pores and cracks are cemented by fine-grained, authigenic phyllosilicates (arrows). The orange box indicates location of (k). (k) Matrix-cementing authigenic phyllosilicates are mostly flake-to needle-shaped chlorite crystallites. ( Fig. 7c and d), except that clasts are not homogeneously distributed but cluster randomly at Havelock Creek. Apart from the hanging wall-proximal fault gouge layer of Waikukupa Thrust, where typical clast size decreases while the size of matrix clasts increases towards the footwall, trends regarding microstructures are not observed. Most remarkably, discrete slip planes are not observed within the fault gouges.

Footwall gravels
Where sufficiently exposed, the contact between clast-poor PSZ and clast-bearing footwall is sharp (Figs. 4b-d; 5b, c; 7j-m). Footwall gravels are compositionally similar to overlying hanging wall cataclasites and fault gouges (Table 3). Furthermore, the same correlations regarding amount and size of clasts observed in hanging wall cataclasites and fault gouges are seen within footwall gravels (Table 3): locations with thin PSZ tend to contain more and larger clasts than locations with thicker PSZ. In addition, grain size increases with increasing distance from the PSZ. Whereas there are some bright matrix clasts at Martyr River and Gaunt Creek, matrix clasts are only encountered at Gaunt Creek. At Mar- tyr River, footwall clasts are slightly imbricated and form a weak foliation parallel to displacement.
Calcite is absent at Martyr River and restricted to a ∼ 3 mm thick layer immediately adjacent to the PSZ at Waikukupa Thrust. There, calcite constitutes finely dispersed cement, veinlets and < 5 µm thick rims at grain edges parallel to displacement. In contrast, at Gaunt Creek, calcite is predominantly found in veins and lenses cross-cutting the matrix and detrital clasts.
The ∼ 1 cm thick, fine-grained, brown layer erroneously interpreted as PSZ at Gaunt Creek ( Fig. 3d and e), is mainly composed of fractured, subangular to subrounded matrix clasts (Fig. 6d). The amount of quartzofeldspathic clasts within this lens is < 5 %. Fractures between the individual clasts are filled with a fine-grained matrix. SEM analysis reveals that these clasts are actually bright matrix clasts (Fig. 6e).

Mineralogy
There are only minor variations in the qualitative fault rock composition (Table 4). Quartz, feldspar, calcite and phyllosilicates (chlorite, kaolinite, muscovite, illite and biotite) and traces of apatite, pyrite and rutile are always present. Additionally, Waikukupa Thrust and Havelock Creek contain amphiboles, mainly hornblende, and there is epidote at

Hanging wall cataclasites
Quartz and feldspars are the dominant mineral phases in hanging wall cataclasites. They do not exhibit any clear trend towards the PSZ (Table 4). Phyllosilicates appear to decrease in content towards the PSZ. This probably is an artifact of Rietveld refinement based on the < 10 µm instead of the claysize fraction. In the DFDP-1A drill core material, the concentrations of detrital muscovite/illite decrease towards the PSZ compared to authigenic chlorite/kaolinite, which show an opposite trend. Furthermore, calcite also increases towards the PSZ at DFDP-1A.

Principal slip zone fault gouges
Quartz and feldspar contents vary by more than 50 % among individual locations (21 wt % at Waikukupa Thrust vs. 52 wt % at Havelock Creek; Fig. 8; Table 4). In general, these clastic phases tend to be less abundant in thinner than in thicker PSZs (Fig. 8). Phyllosilicate concentrations exhibit the opposite pattern: thinner PSZs tend to be more phyllosilicate-rich than thicker PSZs. Furthermore, chlorite/kaolinite and muscovite/illite show inverse correlations but do not vary systematically with respect to PSZ thickness. Calcite concentrations of Waikukupa Thrust (29 %) and Gaunt Creek (6 %) reflect mineralization manifested in vein networks. As calcite veins are rare at Havelock Creek, the relatively high (11 %) concentrations potentially reflect calcite cementation. In contrast, calcite concentrations within PSZs of DFDP-1A and Martyr River are lower compared to the hanging wall.

Footwall gravels
At Gaunt Creek, quartz and feldspar concentrations increase with increasing distance from the PSZ, while those of phyllosilicates decrease, and calcite concentrations remain constant (Table 4). Conversely, at Waikukupa Thrust, footwall concentrations of clastic phases and calcite are comparable to those of the hanging wall close to the contact with the PSZ, and phyllosilicates are slightly more abundant in the footwall compared to the hanging wall adjacent to the PSZ. Compared to the PSZ, clastic phases are significantly more abundant and amounts of phyllosilicate are lower.

Geochemistry
Element concentrations commonly scatter and, if at all, increase or decrease only slightly towards the PSZs (Table 5). The isocon analysis reveals a slight but recognizable correlation between the amount of fluid-related alteration and PSZ thickness: using the slope of the isocon as a proxy to assess the overall degree of fluid-related alteration, with a larger deviation from one indicating a higher amount of alteration, the hanging walls and footwalls of locations characterized by a thinner PSZ are generally less affected by metasomatism ( Fig. 9; Tables 6 and S1). The PSZs display the same trend, but less pronounced. This is especially valid for the PSZ at Waikukupa Thrust, which is affected by a fluid-related mass change of up to 68 % (Table 6). However, microstructural observations (see Sect. 4.2.2 and Schuck et al., 2018) suggest this is most likely related to the extraordinarily large amount of calcite veins.
Enrichment of individual elements confirms this correlation: enrichment of Al 2 O 3 , CaO, LOI, MnO, P 2 O 5 and TiO 2 in hanging wall and footwall rocks is more pronounced at locations with thicker PSZ compared to those with thinner PSZ (Tables 6, S2). Element enrichment of PSZ fault gouges exhibits a similar but less pronounced trend. Comparable to the general pattern of isocon slopes described above, the PSZ at Waikukupa Thrust represents an outlier, because the large volume of calcite veins results in extraordinarily enriched CaO.
Looking at cross-fault transects, PSZs are more affected by metasomatism than hanging walls and footwalls. Furthermore, footwalls are generally less affected by metasomatism than hanging walls. However, DFDP-1A constitutes an exception, because the PSZ is less affected by metasomatism than the hanging wall cataclasites: the degree of fluid-related alteration increases towards the PSZ and displays the highest  hanging wall amounts immediately adjacent to it (Table 6). This contrasts with Waikukupa Thrust, where samples display a continuing trend to higher amounts of fluid-related alteration from the hanging wall to the base of the PSZ close to the footwall.

Fluid-related alteration and fluid transport within Alpine Fault rocks
Isocon analyses clearly demonstrate that all investigated fault rocks have been substantially altered by fluids ( Fig. 9; Table 6), which is also shown by the presence of authigenic phyllosilicates ( Fig. 6a and k) and calcite vein networks ( Fig. 7j-m). In general, metasomatic alteration of the fault Table 6.
Slope m of the isocon and derived mass change (%), enrichment (positive values) and depletion (negative values) of elements for samples investigated relative to the host rocks.
Geochemical composition of fault rocks is taken from Table 5 and composition of host rocks from rocks resulted in element enrichment and mass gain, which is in agreement with previous work (Schuck et al., 2018). Analyses of Martyr River fault rocks seem to show an opposite pattern but this may be an artifact related to inappropriate choice of the protolith's geochemical composition (see Sect. 3.3.1). It is recognized that the Alpine Fault's PSZ constitutes an impermeable barrier restricting fluid circulation to its hanging wall (Sutherland et al., 2012;Menzies et al., 2016). Consequently, fluid-related element mobilization can be expected in the hanging wall. This is confirmed by isocon analyses: DFDP-1A samples show increasing amounts of fluid-related alteration towards the PSZ with peak values immediately above it. However, at all other locations the fault gouges are more affected by fluid-related alteration than either hanging wall or footwall rocks (Fig. 9). This suggests, in combination with networks of anastomosing calcite veins (Fig. 7j-m; Schuck et al., 2018), that there is repeated transient fluid flow within the PSZ, most likely along-fault, because the veins are parallel to the shear plane.
Overall, results of isocon analyses confirm outcomes of previous studies, which demonstrated that fluid-related alteration significantly changed the petrophysical, mineralogical, geochemical and geomechanical properties of the Alpine Fault Zone (e.g., Townend et al., 2013;Carpenter et al., 2014). This supports the proposal of Sutherland et al. (2012) that an "alteration zone" within ∼ 50 m of the PSZ should be included as a fundamental and additional part of the Alpine Fault Zone architectural model (Fig. 2c).
Except for the slickolite observed in the DFDP-1A fault core (Fig. 7g), unambiguous microstructural evidence for pressure solution such as dissolution seams or indented and embayed clasts was not found. Annual rainfall along the NW edge of the Southern Alps is very high (e.g., Norris and Cooper, 1997, and references therein), and we suggest the observed dissolution of grain edges is related to weathering rather than deformation processes.
In summary, both the microstructures and the predominant enrichment of elements show that fluids are not responsible for stress-driven dissolution processes or for substantial mass transfer out of the fault zone. However, given that the amount of phyllosilicates and the degree of fluid-related alteration are negatively correlated (Figs. 8 and 9), the results of this study do not identify alteration mechanisms that completely explain the observed element mobilization within the fault rocks, except at Waikukupa Thrust and Gaunt Creek, where Ca enrichment is associated with the formation of calcite networks ( Fig. 7j-m).

Influence of shallow-depth conditions on fault gouge structure
In the following, we consider to what extent the fault gouges formed at shallow depths, thus do not provide information about deeper fault zone architecture.
Fault gouges at Waikukupa Thrust and Martyr River are layered. Continuous transects across the fault core would sample contacts of fault gouge with both hanging wall and footwall rocks. Furthermore, they would allow for microstructures to be described across the entire PSZ. However, these continuous transects are only available from Waikukupa Thrust and Martyr River. Consequently, the absence of layered fault gouges at other locations seems to be more appropriately explained by a sampling bias than by fault gouge layers being site-specific features. Cowan et al. (2003) analyzed low-angle detachment faults in Death Valley. There, layered fault gouges, consisting of a hanging wall-and a footwall-proximal layer, separate footwall cataclasites from weakly consolidated hanging wall sediments. Cowan et al. (2003) compare the footwall-proximal gouge layer to the hanging wall-proximal layer (i.e., the layer juxtaposed to the weakly consolidated sediments), finding the former to be fine-grained, clast-poor and supposed to have accommodated most of the strain. They interpret the formation of these layered fault gouges to reflect changing p-T conditions during exhumation, which in combination with fluid-related alteration, formed authigenic phyllosilicates within the PSZs before fault rocks were brought in contact with weakly consolidated hanging wall sediments. This juxtaposition resulted in mechanical layering, enabling the hanging wall-proximal layer to accommodate most of the strain.
By considering the different fault kinematics of the Alpine Fault compared to these faults in Death Valley (thrust vs. low-angle detachment fault), the Alpine Fault with its layered gouges described in this and previous work (e.g., Schuck et al., 2018) is structurally identical to the faults analyzed in Death Valley ( Fig. 10; for details see Biegel and Sammis, 2004). Furthermore, studied fault gouges bound footwall gravels, which are poorly consolidated and thus assumed to be mechanically weak. In addition, fault gouge layers juxtaposed on these weakly consolidated rocks (footwall-proximal layer at the Alpine Fault; hanging wall-proximal layer in Death Valley) are microstructurally comparable, as they contain fewer and preferentially smaller clasts than the other layers of the fault gouges (see also Table 2). Furthermore, microstructures of the calcite network in the footwall-proximal layer at Waikukupa Thrust suggest localized and cyclic faulting (Schuck et al., 2018) supporting the interpretation that this gouge unit accommodated most of the strain. Consequently, layered fault gouges and associated microstructures demonstrate that strain localization and associated structural complexities within the Alpine Fault core are, at least partially, affected by shallow, hence late, processes, considering that exhumation of investigated rocks to subaerial level is accommodated entirely by the Alpine Fault's PSZ (e.g., Little et al., 2005;Toy et al., 2015).
This interpretation has major consequences for structural investigations of the Alpine Fault's PSZ, because it im-  Table S1. Figure 10. Conceptual model providing a potential explanation of how layered principal slip zones observed at Martyr River and Waikukupa Thrust could have been formed. Resulting from the juxtaposition of weakly consolidated sediments in the footwall and relatively competent fault rocks, strain localized non-homogenously in the PSZ. The photograph example is from Waikukupa Thrust. The stippled line indicates where the contact between cataclasite and gouge is obscured by debris. The figure was modified from Cowan et al. (2003). plies that formation of microstructures observed in outcrops was influenced by boundary conditions not representative of those along most of the fault rocks' exhumation path. This highlights the importance of continuing to attempt to obtain samples for microstructural investigations of Alpine Fault rocks from locations and/or depths where hanging wall and PSZ fault rocks are not juxtaposed on weakly consol-idated footwall sediments. Acknowledging that there is no outcrop where the hard rock hanging wall is in contact with hard rock footwall (Townend et al., 2009), drilling projects aiming to penetrate the Alpine Fault at depths where Quaternary gravels do not constitute the footwall appear to be the only promising approach to achieve this end. This would also overcome the problem of data sparsity due to the limited amount of outcrops, improving the robustness of the results. Norris and Cooper (2007) suggested an average brittle damage zone width of ∼ 100 m, which is based on a compilation of outcrop studies in the fault's central segment. More recent investigations indicate wider damage zones: Alpine Fault outcrops in the segment south of Haast have damage zone widths between 90 and 240 m (Barth et al., 2013). Analyses of downhole geophysical data and recovered cores from the DFDP-2 borehole indicate that the Alpine Fault's hydrologically active damage zone contains an "inner zone" affected by earthquake rupture processes that is 0.5 km wide at up to 8 km depth and about 1 km wide at the surface (Fig. 2c). This inner damage zone is surrounded by an "outer zone" up to a few kilometers wide, which might only be a nearsurface, topographically controlled feature Massiot et al., 2018;Williams et al., 2018).

Alpine Fault damage zone
Outcrop and drill-core investigations demonstrate quite constant fracture densities within at least 500 m of the PSZ in the hanging wall and 30 m in the footwall (Williams et al., 2016, which is different from decreasing fracture densities towards the fault cores observed elsewhere (e.g., Chester et al., 2005;Faulkner et al., 2006;Mitchel and Faulkner, 2009). In fact, DFDP investigations show that apparent fracture density decreases with depth, i.e., towards the PSZ (Townend et al., 2013;Williams et al., 2016). However, this is considered as an artifact resulting from rapid exhumation and associated unloading (Williams et al., 2016). Furthermore, despite its presence across the entire damage zone, gouge-filled fractures within < 160 m of the PSZ tend to be thinner (< 1 cm) and more abundant .

Extent of fault core
As a result of the pervasive alteration zone surrounding the PSZ, the fault core of the Alpine Fault is rather diffuse (Sutherland et al., 2012;Townend et al., 2013;Williams et al., 2018). Commonly, the fault gouge forming the PSZ and part of the surrounding cataclasites are considered to constitute the fault core (e.g., Williams et al., 2017). However, there is some ambiguity when it comes to actually defining the width of the fault core: Toy et al. (2015) consider a 20-30 m wide zone around the PSZ as fault core, and Sutherland et al. (2012) and Townend et al. (2017) narrow it down to 2 m surrounding the PSZ with an additional degree of localization manifested in the < 0.5 m wide PSZ (Sutherland et al., 2012;Williams et al., 2016). However, by defining the damage zone as having elevated fracture densities compared to the host rock (Chester et al., 1993;Faulkner et al., 2010), and noting that there is no observable variation in fracture density within 30 m of the PSZ (Williams et al., 2016), it seems reasonable to exclusively consider the fault gouge forming the PSZ as core of the Alpine Fault's central segment. This view is supported by fault gouge microstructures: microfaults ( Fig. 6h and i), predominantly encountered in hanging wall cataclasites but also present in fault gouges (Table 3), show that strain is localized to discrete structures within the PSZ. Furthermore, the fault breccias encountered in these microfaults exhibit a broad range in degree of sorting (poorly to very well) and grain size (< 2-100 s of micrometer), which demonstrates that they accommodated different amounts of strain. However, microfaults appear to be a local phenomenon and there are no continuous, discrete slip planes within the PSZ, despite the large finite displacement accommodated.
Implications of fault core along-strike variations for fault zone architecture Investigated fault rocks are mostly comparable in terms of microstructures and mineralogy (Tables 1-4). However, individual outcrops display distinct differences, which will be discussed in the context of the Alpine Fault Zone's architecture in the following.
As will be elaborated in the following, these explanations do not appear to be responsible for the observed variations. PSZ thickness does not vary systematically along-strike in the sampled locations (Fig. 2b). Furthermore, all locations, except Martyr River, are in the central segment of the Alpine Fault where the Torlesse Terrane Alpine Schist is the protolith. This demonstrates that variations of protolith mineralogy have no significant influence on the observed differences in PSZ thickness. However, fault gouge thickness correlates with mineralogy: thicker gouges at Havelock Creek, Gaunt Creek and the DFDP-1A core are richer in clastic phases and poorer in phyllosilicates than the thin gouges of Martyr River and Waikukupa Thrust (Fig. 8). This is astonishing, because elsewhere phyllosilicate-rich fault gouges generally tend to be wider (Faulkner et al., 2003) owing to the frictionally weak nature of these phases Moore and Lockner, 2007;Lockner et al., 2011;Carpenter et al., 2015). Conversely, other fault zones with quartzofeldspathic compositions like the Alpine Fault have been found to be fairly narrow (Chester et al., 1993;Faulkner et al., 2003). Consequently, unlike these fault zones, the Alpine Fault displays a negative correlation between phyllosilicate content and fault core thickness. This suggests that strain localization within the fault core might be governed by processes insensitive to rheological variations caused by differing fault rock composition, although observed slight microstructural differences such as grain size (see also Tables 1-3) could influence localization. However, data are too scarce to derive robust conclusions. In addition, orientation of the fault ( Fig. 2a and b), magnitude of stress and the stress field itself (Boese et al., 2012;Warren-Smith et al., 2017) are fairly constant along the central segment of the Alpine Fault, thus it is less likely that these parameters are responsible for the observed variations of fault gouge thickness. However, it is still possible that the extreme variations in topography typical of the Alpine Fault hanging wall generated quite different stress fields at our investigated sites (see Upton et al., 2017). Biegel and Sammis (2004) suggested that along-strike variations of gouge thickness record rupture arrest. However, our observations (and the nature of Alpine Fault outcrops in general) do not easily lend themselves to a systematic analysis of this kind in the Alpine Fault Zone.
Another approach to explain observed thickness variations is to examine the relation between fault core width and displacement. Wear models based on empirical observations describe a positive correlation between gouge thickness and displacement (e.g., Hull, 1988;Scholz, 1987). However, quantitative models of fault core evolution have so far failed to reproduce data sets compiled on natural faults (e.g., Blenkinsop, 1989;Evans, 1990;Sibson, 2003). Nevertheless, faults that accommodated larger displacements do appear to generally have thicker fault cores (Evans, 1990;Ben-Zion and Sammis, 2003;Sagy and Brodsky, 2009;Faulkner et al., 2010, and references therein). In other words: the more earthquakes a fault has seen, the thicker its PSZ may be (Ma et al., 2006;Li et al., 2013). Four of the five fault gouges investigated are located in the Alpine Fault's central segment, which is believed to rupture entirely during an earthquake (Suther-land et al., 2007;Howarth et al., 2018) and to have done so for at least the last ∼ 8000 years (Berryman et al., 2012). If we employ the simple relationship that shear strain in a fault zone is equal to zone width multiplied by boundary displacement, and assume constant strain distribution within the PSZ gouges, the PSZ thickness variations we describe suggest the studied fault gouges accommodated different amounts of coseismic displacement.
Other examples where PSZ thickness varies by 1 order of magnitude within tens of meters along strike are not comparable to the Alpine Fault, because these faults accommodated substantially less displacement than the Alpine Fault: Shervais and Kirkpatrick (2016) and Kirkpatrick et al. (2018) investigate a fault which accommodated displacement of the order of 10-30 km. The fault core contains at least 13 individual slip layers, suggesting that slip accommodated by these individual slip layers is most likely substantially smaller than the overall fault displacement. Consequently, observed PSZ thickness variations potentially reflect fault roughness, which is supposed to smooth as displacement increases (Sagy and Brodsky, 2009;Brodsky et al., 2011). In contrast, the Alpine Fault not only accommodated at least 470 km of displacement (Sutherland et al., 2007), its fault gouge also does not display discrete slip planes (see Sect. 4.2.2). Furthermore, contacts between the fault gouge and the hanging wall and between the fault gouge and the footwall are -where sufficiently exposed -smooth (Figs. 4b-d and 5b-d).
Consequently, Alpine Fault PSZ thickness variations reflecting different amounts of displacement imply that the studied PSZs are not part of the same shear plane, and that the Alpine Fault Zone hosts multiple fault strands, which may not be active simultaneously (Figs. 1b and 11). Figure 11 demonstrates how difficult it will be to identify the fault strand active during the most recent event.
Further evidence for a more complex Alpine Fault Zone architecture This assumption of a more complex Alpine Fault Zone geometry is further supported by microstructural, mineralogical and geophysical observations to be discussed in the following. Most obviously, DFDP results promote a more complex view of the fault zone architecture of the Alpine Fault: DFDP-1B is the only location in the fault's central segment where there are two PSZs, both bound by brittle fault rocks ( Fig. 3a; Sutherland et al., 2012;Toy et al., 2015). This is substantially different from other locations at the Alpine Fault, where multiple fault strands have been identified in Quaternary sediments, which are poorly consolidated and several tens of meters to a few hundred meters thick, covering the hanging wall and the footwall (Kaiser et al., 2009;Carpentier et al., 2012). The presence of these two PSZs in DFDP-1B could reflect that displacement was accommodated by more than one PSZ, which is dissimilar to the paradigm of uplift along a single fault plane (e.g., Little et al., 2005;Toy et al., 2015), and/or that one PSZ was abandoned and displacement transferred to another PSZ. Fault gouge dating could help answer this question, but is beyond the scope of this paper.
Furthermore, although the shallower DFDP-1B PSZ and the PSZ of DFDP-1A have so far been assumed to belong to the same fault plane (Fig. 3a), they exhibit different mineralogical and hence geomechanical properties (Boulton et al., 2014;Schleicher et al., 2015). This is unexpected given the close proximity of the boreholes and thus sample spacing of less than 100 m (Sutherland et al., 2012), because the distinct mineral compositions of these two PSZs represent high (DFDP-1A) and low (DFDP-1B) temperature alteration, respectively. Considering the proximity of the two drill sites, the distinct mineralogies reflecting different temperaturesand hence depths -of formation appear to go beyond what one would expect as common along-strike variation observed in the field. Similarly, although mineralogically comparable (Fig. 8), fault gouge microstructures of DFDP-1A and Gaunt Creek, also within only ∼ 100 m distance (Fig. 3a), differ remarkably, as exemplified by the almost completely cemented DFDP-1A PSZ compared to the PSZ at Gaunt Creek (Table 2; Fig. 7e, f, k-m).
In light of this interpretation, the ∼ 1 cm thick, matrix clast-bearing layer sampled ∼ 8 cm below the PSZ at Gaunt Creek (Figs. 3d, e, 5b and 6d), which is microstructurally similar to the investigated fault gouges at Havelock Creek and Waikukupa Thrust (hanging wall-proximal layer), could be considered to constitute another PSZ of the Alpine Fault, currently starting to accommodate larger amounts of slip, or an inactive one as observed in the fault's southern segment (Barth et al., 2013).
Geometric and spatial relationships, mineralogical composition and geomechanical properties of the various (potential) PSZs at Gaunt Creek/DFDP-1 (Fig. 3a) support our hypothesis that displacement along the Alpine Fault is accommodated by more than one PSZ.
The Alpine Fault's more complex geometry: a shallow-depth phenomenon?
The multidisciplinary approach presented here reveals that the Alpine Fault Zone has a complex geometry (see Figs. 1b and 11), which is different from other quartzofeldspathic faults but typical for carbonate or phyllosilicate-hosted ones (Faulkner et al., 2003). However, it cannot be excluded that the observations we present in favor of a more complex fault zone geometry simply reflect shallow-depth phenomena like layered fault gouges do. Active seismic source investigations at the DFDP-2 site and around Haast also support a more complex fault zone model. These studies imaged a series of subparallel reflectors, interpreted as multiple fault strands (Lay et al., 2016;, which are comparable to those hosted by the Quaternary sediments covering the hanging wall and the footwall (Kaiser et al., 2009;Carpentier et al., 2012). There are also indications that these fault strands at the DFDP-2 site and around Haast extend into the basement. So far, it is debated if these separate fault strands are isolated or anastomosing and if they are/were active contemporaneously or consecutively . Note that this structural complexity differs from the situation encountered at Waikukupa Thrust, as the second fault strand encountered at Hare Mare Creek, 700 m to the NE of this location, resulted from incision of the Waikukupa River and subsequent stress reorientation (Norris and Cooper, 1997).

Conclusions
Layered fault gouges reflect a structural complexity of the Alpine Fault's PSZ, which originates from shallow-depth processes. At depths where Quaternary gravels comprise the footwall, rheological contrasts resulting from the juxtaposition of these poorly consolidated and mechanically weak footwall sediments with more competent hanging wall cataclasites control the distribution of strain within the fault gouge and result in the layered structure at shallow depths. Consequently, rheological interpretations based on these rocks are biased as their microstructures mostly record modifications in just the latest part of the exhumation trajectory. These results also show that careful analyses are necessary in order to study actively exhuming fault zones and that samples should be taken at greater depths to fully investigate these kinds of faults.
Variations of PSZ fault gouge thickness among individual study sites suggest that fault rocks at these locations accommodated different amounts of displacement. Assuming that the Alpine Fault's central segment, where most of the sites are located, is activated entirely during an earthquake, these thickness variations suggest that investigated PSZs are not located on the same fault plane. This implies that multiple fault strands accommodate Alpine Fault displacement. These findings, which are based on field and microstructural analyses, are further supported by the results of the Deep Fault Drilling Project and geophysical investigations.
Consequently, these results show that the Alpine Fault Zone architecture is more appropriately described by the broad and complex conceptual model of Faulkner et al. (2003) rather than the simple single PSZ model of Caine et al. (1996).
Data availability. Underlying diffractograms and supporting data of isocon analysis for this paper can be found in the Supplement.
Author contributions. BS, CJ and VGT conducted field work and obtained the samples. Microstructures were analyzed by BS. BS and AMS performed XRD and XRF analyses. All authors contributed to manuscript preparation but BS had primary responsibility.