David J. A. Barrell

We combine previous studies, fieldwork, lidar data, and trenching to examine late Quaternary activity of the SE‐dipping ≥23 km long Settlement Fault in the southeastern South Island. Trenching of a scarp crossing a small alluvial fan exposed a > 3 m thick sequence of folded sandy to silty fan alluvium. The uniformity of folding and four Optically Stimulated Luminescence dates indicate that these sediments were deposited relatively rapidly at ∼20 ka and subsequently accommodated 2.5 m of throw following ≥1 ruptures on an underlying buried reverse fault. Other key constraints on late Quaternary Settlement Fault activity are: (1) stranding of a shore platform at ∼3.7 ka, due to 1.5–2 m of fault uplift, and (2) a flat bedrock erosion surface overlain by loess and correlated to the last interglacial maximum (∼125 ka), indicating 0–4 m of fault uplift. From incorporating this information into a stochastic model of fault displacement accumulation we derive a slip rate of ∼0.02 mm/yr since 125 ka, and ∼0.27 mm/yr since 20 ka. We link the Settlement Fault’s current phase of heightened activity to a similarly timed increase in Akatore Fault activity, despite the tips of these faults being separated by between 9 and 25 km.

The New Zealand Active Faults Database (NZAFD) contains underpinning data to help mitigate the impacts of future surface‐rupturing earthquakes in Aotearoa New Zealand. However, defining the associated hazards and risks must be undertaken at relevant scales and as such, the NZAFD contains two scale‐based datasets each serving complementary, but different, purposes. The high‐resolution (‘NZAFD‐HighRes’) dataset contains enough detail on surface traces for cadastral scale land‐use planning purposes, while the other dataset is generalised to 1:250,000 scale (‘NZAFD‐AF250’). Here we document for the first time the NZAFD‐HighRes dataset (v2.0) and describe the recent efforts that have focused on updating the dataset structure to increase useability, compiling data and improving GIS infrastructure. The NZAFD‐HighRes and NZAFD‐AF250 datasets, along with Fault Avoidance Zones and Fault Awareness Areas – land‐use planning tools used to mitigate surface rupture hazard – are publicly accessible via the active faults web service and a webmap application at https://data.gns.cri.nz/af/ . This upgraded webmap enables active fault data to be discovered and used for informing future surface rupture hazard assessments. The relationship of the NZAFD to other active fault datasets and models is also discussed, along with future directions and challenges.

We use new paleoseismic data and lidar to reassess late Quaternary activity of the NW (northwest) Cardrona Fault, a ∼60 km long range‐bounding fault in Otago. Paleoseismic investigations of the NW Cardrona Fault were conducted in the 1980s, but findings were limited by a paucity of materials suitable for dating. Here, re‐exposure and re‐assessment of a 1980s trench at Macdonalds Creek provide stratigraphic evidence for two surface rupturing earthquakes on the NW Cardrona Fault. Through Optically Stimulated Luminescence dating and OxCal modelling, we estimate these earthquakes occurred at (mean ± 2 standard deviations, σ) 20.5 ± 7.8 ka and 8.7 ± 6.6 ka. The timings of these earthquakes are consistent with those proposed from re‐interpretation of another 1980s trench on the NW Cardrona Fault in the Kawarau River valley, ∼25 km farther south. The slip rate and recurrence interval derived from the Macdonalds Creek trench data are 0.13 ± 0.05 mm/yr and 11,800 ± 5100 years (2σ error) respectively. These estimates imply lower slip rates and longer earthquake recurrence intervals on the NW Cardrona Fault than previously assessed.

The New Zealand Community Fault Model (NZ CFM) is a publicly available representation of New Zealand fault zones that have the potential to produce damaging earthquakes. Compiled through collaborative engagement between New Zealand earthquake‐science experts, this first edition (version 1.0) of the NZ CFM builds upon previous compilations of earthquake‐source active fault models with the addition of new and modified information. Developed primarily to support an update of the New Zealand National Seismic Hazard Model, the NZ CFM comprises two principal components. The first dataset is a two‐dimensional map representation of the surface traces of 880 generalised fault zones. Each fault zone is assigned specific geometric and kinematic attributes, including uncertainties, supplemented with a subjective quality ranking focused primarily on the confidence in assigned slip rates. The second component is a three‐dimensional representation of the fault zones as triangulated mesh surfaces that are projected down‐dip from the two‐dimensional mapped traces to a geophysically‐defined maximum fault rupture depth. This article summarises the compilation and parameterisation of the NZ CFM, along with background on its relation to predecessor datasets, and forward applications to probabilistic seismic hazard assessment and physics‐based earthquake models currently being developed for Aotearoa New Zealand.

The Hundalee Fault forms part of the southeastern margin of the 2016 Kaikōura Earthquake rupture zone. Its late Quaternary activity and structural character may offer insights to fault interrelationships associated with the 2016 rupture. Mapping of the Hundalee Fault revealed several pre‐existing fault scarps. Trenching of a scarp across a fluvial terrace together with radiocarbon dating shows the occurrence of at least two, probably three, surface ruptures since 3.5 ka. The scarp records about 2.5 m of reverse‐sense vertical deformation, mostly folding, of which ∼0.6 m occurred in 2016, along with subordinate sinistral slip. Fluvial terrace ages inferred from geomorphological relationships, together with the trenching results, indicate a long‐term vertical slip rate of 0.2–0.4 mm/yr since at least 30 ka. Geological data indicate a relatively recently surface emergence of the Hundalee Fault from beneath a fault propagation fold. The 2016 Hundalee Fault rupture is interpreted as a break‐out from an extensive blind thrust fault, whose motion triggered predominantly contractional rupture of an array of steeper near‐surface faults. The thrust is hypothesised as having provided linkage in the 2016 rupture between the predominantly strike‐slip ruptures of The Humps Fault to the southwest and the Kekerengu Fault to the northeast.

Reverse faulting in Otago, southern New Zealand, accommodates distributed tectonic convergence on the eastern side of the Australian‐Pacific plate boundary. Paleoearthquake records from some of the faults in the region show highly variable earthquake recurrence times, with long periods of quiescence separated by periods of earthquake activity. Here we develop a longer‐term context for these records, using cosmogenic radionuclide dating of faulted alluvial fan surfaces to characterize the late Quaternary slip rates on two significant faults within the system, the Hyde and Dunstan faults. We determine an average slip rate of 0.24 mm/yr (0.19–0.29 mm/yr at 95% confidence) for the Hyde Fault since about 115 ka, and an average slip rate of 0.16 mm/yr (0.12–0.21 mm/yr) for the Dunstan Fault since about 320 ka. Both faults show increases in slip rate of a factor of three to five times the average long‐term rate over timescales of 10 kyr. Increases in slip rate are out of phase on the two faults, supporting a hypothesis that strain is shared within the fault system over timescales on the order of 10 kyr. Over longer timescales, on the order of 100 kyr, slip rates can be well‐approximated by a linear fit, providing an upper limit on the variability of fault slip rates with time.

Determining how the crustal structure developed in previous tectonic regimes influences current deformation is important for understanding seismogenic behavior. Paleozoic‐Mesozoic convergent‐margin tectonics produced the southeastern South Island crust, which is currently involved in late Cenozoic transpressional collision. We examine crustal velocity structure through joint inversion of earthquake travel times and ambient noise derived group velocities to increase the resolution of Vp (P‐wave velocity) and Vp/Vs across the region. This method improves the details of the shallow velocity structure and highlights lithological properties, uplifted regions, and distributed fracture zones. The Otago schist shows thick lower crust while contrasting higher Vp below 15 km depth characterizes the Matai terrane and Median batholith, both backstops during accretion and deformation. The complex velocity structure, integrated with late Cenozoic faults, structural domains, seismicity patterns, and geodynamic modeling, shows how Cenozoic structures relate to deeper properties. Seismicity within the schist is distributed within transition regions rather than along specific faults, and these zones contain low Vp and low Qp (1/attenuation). The most prominent seismicity zone is along the transition from northern Otago ranges to the Southern Alps, where faulting and brittle deformation corresponds to the eastern edge of ductile crustal thickening above the bending of the strong oceanic lower crust. The thick weak crust rheology of Otago, relative to Canterbury, enables broad deformation, with numerous late Cenozoic faults that may not persist to great depth. Along the southern stronger plutonic terranes, there is a band of sparse deeper (15–35 km) seismicity, consistent with strong elastic brittle crust.

Areas that experience permanent ground deformation in earthquakes (e.g., surface fault rupture, slope failure, and/or liquefaction) typically sustain greater damage and loss compared to areas that experience strong ground shaking alone. The 2016 Mw 7.8 Kaikōura earthquake generated ≥220 km of surface fault rupture. The amount and style of surface rupture deformation varied considerably, ranging from centimetre-scale distributed folding to metre-scale discrete rupture. About a dozen buildings – mainly residential (or residential-type) structures comprising single-storey timber-framed houses, barns and wool sheds with lightweight roofing material – were directly impacted by surface fault rupture with the severity of damage correlating with both local discrete fault displacement and local strain. However, none of these buildings collapsed. This included a house built directly atop a discrete rupture that experienced ~10 m of lateral offset. The foundation and flooring system of this structure allowed decoupling of much of the ground deformation from the superstructure thus preventing collapse. Nevertheless, buildings directly impacted by surface faulting suffered greater damage than comparable structures immediately outside the zone of surface rupture deformation. From a life-safety standpoint, all these buildings performed satisfactorily and provide insight into construction styles that could be employed to facilitate non-collapse performance resulting from surface fault rupture and, in certain instances, even post-event functionality.