On 4 September 2010, a magnitude Mw 7.1 earthquake struck the Canterbury region on the South Island of New Zealand. The epicentre of the earthquake was located in the Darfield area about 40 km west of the city of Christchurch. Extensive damage was inflicted to lifelines and residential houses due to widespread liquefaction and lateral spreading in areas close to major streams, rivers and wetlands throughout Christchurch and Kaiapoi. Unreinforced masonry buildings also suffered extensive damage throughout the region. Despite the severe damage to infrastructure and residential houses, fortunately, no deaths occurred and only two injuries were reported in this earthquake. From an engineering viewpoint, one may argue that the most significant aspects of the 2010 Darfield Earthquake were geotechnical in nature, with liquefaction and lateral spreading being the principal culprits for the inflicted damage. Following the earthquake, an intensive geotechnical reconnaissance was conducted to capture evidence and perishable data from this event. This paper summarizes the observations and preliminary findings from this early reconnaissance work.
The 22 February 2011, Mw6.2-6.3 Christchurch earthquake is the most costly earthquake to affect New Zealand, causing 181 fatalities and severely damaging thousands of residential and commercial buildings, and most of the city lifelines and infrastructure. This manuscript presents an overview of observed geotechnical aspects of this earthquake as well as some of the completed and on-going research investigations. A unique aspect, which is particularly emphasized, is the severity and spatial extent of liquefaction occurring in native soils. Overall, both the spatial extent and severity of liquefaction in the city was greater than in the preceding 4th September 2010 Darfield earthquake, including numerous areas that liquefied in both events. Liquefaction and lateral spreading, variable over both large and short spatial scales, affected commercial structures in the Central Business District (CBD) in a variety of ways including: total and differential settlements and tilting; punching settlements of structures with shallow foundations; differential movements of components of complex structures; and interaction of adjacent structures via common foundation soils. Liquefaction was most severe in residential areas located to the east of the CBD as a result of stronger ground shaking due to the proximity to the causative fault, a high water table approximately 1m from the surface, and soils with composition and states of high susceptibility and potential for liquefaction. Total and differential settlements, and lateral movements, due to liquefaction and lateral spreading is estimated to have severely compromised 15,000 residential structures, the majority of which otherwise sustained only minor to moderate damage directly due to inertial loading from ground shaking. Liquefaction also had a profound effect on lifelines and other infrastructure, particularly bridge structures, and underground services. Minor damage was also observed at flood stop banks to the north of the city, which were more severely impacted in the 4th September 2010 Darfield earthquake. Due to the large high-frequency ground motion in the Port hills numerous rock falls and landslides also occurred, resulting in several fatalities and rendering some residential areas uninhabitable.
This study investigates evidence for linkages and fault interactions centred on the Cust Anticline in Northwest Canterbury between Starvation Hill to the southwest and the Ashley and Loburn faults to the northeast. An integrated programme of geologic, geomorphic, paleo-seismic and geophysical analyses was undertaken owing to a lack of surface exposures and difficulty in distinguishing active tectonic features from fluvial and/or aeolian features across the low-relief Canterbury Plains. LiDAR analysis identified surface expression of several previously unrecognised active fault traces across the low-relief aggradation surfaces of the Canterbury Plains. Their presence is consistent with predictions of a fault relay exploiting the structural mesh across the region. This is characterised by interactions of northeast-striking contractional faults and a series of re-activating inherited Late Cretaceous normal faults, the latter now functioning as E–W-striking dextral transpressive faults. LiDAR also allowed for detailed analysis of the surface expression of individual faults and folds across the Cust Anticline contractional restraining bend, which is evolving as a pop-up structure within the newly established dextral shear system that is exploiting the inherited, now re-activated, basement fault zone. Paleo-seismic trenches were located on the crest of the western arm of the Cust Anticline and across a previously unrecognised E–W-striking fault trace, immediately southwest of the steeply plunging Cust Anticline termination. These studies confirmed the location and structural style of north-northeast-striking faults and an E–W-striking fault associated with the development of this structural culmination. A review of available industry seismic reflection lines emphasised the presence of a series of common structural styles having the same underlying structural drivers but with varying degrees of development and expression, both in the seismic profiles and in surface elevations across the study area. Based on LiDAR surface mapping and preliminary re-analysis of industry seismic reflection data, four fault zones are identified across the restraining bend structural culminations, which together form the proposed Oxford–Cust–Ashley Fault System. The 2010–2012 Canterbury Earthquake Sequence showed many similarities to the structural pattern established across the Oxford–Cust–Ashley Fault System, emphasising the importance of identification and characterization of presently hidden fault sources, and the understanding of fault network linkages, in order to improve constraints on earthquake source potential. Improved understanding of potentially-interactive fault sources in Northwest Canterbury, with the potential for combined initial fault rupture and spatial and temporal rupture propagation across this fault system, can be used in probabilistic seismic hazard analysis for the region, which is essential for the suitability and sustainability of future social and economic development.
Since the early 1980s seismic hazard assessment in New Zealand has been based on Probabilistic Seismic Hazard Analysis (PSHA). The most recent version of the New Zealand National Seismic Hazard Model, a PSHA model, was published by Stirling et al, in 2012. This model follows standard PSHA principals and combines a nation-wide model of active faults with a gridded point-source model based on the earthquake catalogue since 1840. These models are coupled with the ground-motion prediction equation of McVerry et al (2006). Additionally, we have developed a time-dependent clustering-based PSHA model for the Canterbury region (Gerstenberger et al, 2014) in response to the Canterbury earthquake sequence. We are now in the process of revising that national model. In this process we are investigating several of the fundamental assumptions in traditional PSHA and in how we modelled hazard in the past. For this project, we have three main focuses: 1) how do we design an optimal combination of multiple sources of information to produce the best forecast of earthquake rates in the next 50 years: can we improve upon a simple hybrid of fault sources and background sources, and can we better handle the uncertainties in the data and models (e.g., fault segmentation, frequency-magnitude distributions, time-dependence & clustering, low strain-rate areas, and subduction zone modelling)? 2) developing revised and new ground-motion predictions models including better capturing of epistemic uncertainty – a key focus in this work is developing a new strong ground motion catalogue for model development; and 3) how can we best quantify if changes we have made in our modelling are truly improvements? Throughout this process we are working toward incorporating numerical modelling results from physics based synthetic seismicity and ground-motion models.
Geomorphic, structural and chronological data are used to establish the late Quaternary paleoseismicity of the active dextral-oblique Northern Esk Fault in North Canterbury, New Zealand. Detailed field mapping of the preserved c. 35 km of surface traces between the Hurunui River and Ashley Head reveals variations in strike ranging from 005° to 057°. Along with kinematic data collected from fault plane striae and offset geomorphic markers along the length of the fault these variations are used to distinguish six structural subsections of the main trace, four dextral-reverse and two dextral-normal. Displacements of geomorphic markers such as minor streams and ridges are measured using differential GPS and rangefinder equipment to reveal lateral offsets ranging from 3.4 to 23.7 m and vertical offsets ranging from < 1 to 13.5 m. Characteristic single event displacements of c. 5 m and c. 2 m have been calculated for strike-slip and reverse sections respectively. The use of fault scaling relationships reveals an anomalously high displacement to surface rupture length ratio when compared to global data sets. Fault scaling relationships based on width limited ruptures and magnitude probabilities from point measurements of displacement imply earthquake magnitudes of Mw 7.0 to 7.5. Optically Stimulated Luminescence (OSL) ages from displaced Holocene alluvial terraces at the northern extent of the active trace along with OSL and radiocarbon samples of the central sections constrain the timing of the last two surface rupturing events (11.15 ±1.65 and 3.5 ± 2.8 ka) and suggest a recurrence interval of c. 5612 ± 445 years and late Quaternary reverse and dextral slip rates of c. 0.31 mm/yr and 0.82 mm/yr respectively. The results of this study show that the Northern Esk Fault accommodates an important component of the c. 0.7 – 2 mm/yr of unresolved strain across the plate boundary within the North Canterbury region and affirm the Esk Fault as a source of potentially damaging ground shaking in the Canterbury region.
The Porter's Pass-Amberley Fault Zone (PPAFZ) is a complex zone of anastomosing faults and folds bounding the south-eastern edge of the transition from subducting Pacific Plate to continental collision on the Australia Plate boundary. This study combines mapping of a 2000 km2 zone from the Southern Alps northeast to the coast near Amberley, 40 km north of metropolitan Christchurch, with an analysis of seismicity and a revision of regional seismic hazard. Three structural styles: 1) a western strike-slip, and 2) a more easterly thrust and reverse domain, pass into 3) a northwest verging fold belt on the northern Canterbury Plains, reflecting the structural levels exposed and the evolving west to east propagation. Basal remnants of a Late Cretaceous-Cenozoic, largely marine sedimentary cover sequence are preserved as outliers that unconformably overlie Mesozoic basement (greywacke and argillite of the Torlesse terrain) in the mountains of the PPAFZ and are underlain by a deeply leached zone which is widely preserved. Structure contouring of the unconformity surface indicates maximum, differential uplift of c.2600 m in the southwest, decreasing to c.1200 m in the coastal fold belt to the northeast. Much lower rates (or reversal) of uplift are evident a few kilometres southeast of the PPAFZ range-front escarpment. The youngest elements of the cover sequence are basement-derived conglomerates of Plio-Pleistocene age preserved on the SE margin. The source is more distant than the intervening mountains of the PPAFZ, probably from the Southern Alps, to the west and northwest. The absence of another regional unconformity on Mesozoic basement, older than Pleistocene, indicates that this uplift is post-Pliocene. Late Pleistocene(<100 kyr) differential uplift rates of c.0.5-2.7 m/kyr from uplifted marine terraces at the east coast, and rates of 2.5-3.3 m/kyr for tectonically-induced river-down cutting further west, suggest that uplift commenced locally during the last 1 Ma, and possibly within the last 0.5 Ma, if average rates are assumed to be uniform over time. Analysis of seismicity, recorded during a 10 week regional survey of micro earthquakes in 1990, identified two seismic zones beneath North Canterbury: 1) a sub-horizontal zone of activity restricted to the upper crust (≤12 km); and 2) a seismic zone in the lower crust (below a ceiling of ≤17 km), that broadens vertically to the north and northwest to a depth of c.40 km, with a bottom edge which dips 10°N and 15°NW, respectively. No events were recorded at depths between 12 km and 17 km, which is interpreted as a relatively aseismic, mid-crustal ductile layer. Marked differences (up to 60°) in the trend of strain axes for events above and below the inferred ductile layer are observed only north of the PPAFZ. A fundamental, north-to-south increase in the Wave-length of major geological structures occurs across the PPAFZ, and is interpreted as evidence that the upper crust beneath the Canterbury Plains is coupled to the lower crust, whereas the upper crust further north is not. Most of the recorded micro earthquakes <12 km deep beneath the PPAFZ have strike-slip mechanisms. It is probable that faults splay upward into the thrusts and folds at the surface as an evolving transpression zone in response to deep shear in basement. There have been no historic surface ruptures of the PPAFZ, but the zone has been characterised historically by frequent small earthquakes. Paleoseismic data (dated landslides and surface ruptures) compiled in this study, indicate a return period of 1500-1900 years between the last two M>7-7.5 earthquakes, and 500-700 years have elapsed since the last. The magnitudes of these events are estimated at c.M7.5, which represents a probable maximum magnitude for the PPAFZ. There are insufficient data to determine whether or not the frequency of large earthquakes conforms to a recognised model of behaviour, but comparison of the paleoseismic data with the historic record of smaller earthquakes, suggests that the magnitudes of the largest earthquakes in this zone are not exponentially distributed. A seismicity model for the PPAFZ (Elder et al., 1991) is reviewed, and a b-value of 1.0 is found to be consistent with the newly acquired paleoseismic data. This b-value reduces the predicted frequency of large earthquakes (M≥7.0) in this zone by a factor of 3.5, while retaining a conservative margin that allows for temporal variations in the frequency of large events and the possibility that the geological database is incomplete, suggesting grounds for revising the hazard model for Christchurch.
Structures of the Lowry Peaks Range - Waikari Valley district are complex. The majority comprise three members of a predominantly WSW -ENE striking major northwards-directed, leading edge imbricate thrust system, with associated angular, asymmetric fault-propagation folds. This system forms anomalously within a large NESW trending belt of structures characterising the entire east coast of north Canterbury, both onshore and offshore and terminates westwards against N-S striking, east facing fold-fault zone. The objectives of this study address the origin, geometry and kinematics of the interaction between these diversely trending systems. Stratigraphy and small-scale structures denote three periods of deformation, namely: i) Middle Cretaceous deformation of the basement rocks, ii) weak Middle Oligocene deformation associated with the inception of the plate boundary through the South Island, and iii) major Pliocene - Recent deformation that formed the majority of the above-mentioned structures. Stress tensor analyses within competent basement and limestone cover rocks suggest two sets of sub-horizontal compression, NE-SW and NW-SE, the former likely to relate to a localised earlier period of deformation, now overprinted by the latter. NW-SE oriented sub-horizontal compression correlates well with results from other parts of north Canterbury. The result of NW-SE compression on the W-E to WSW-ENE striking structures is a large component of oblique motion, which is manifest in four ways: i) movement on two, differently oriented splays rather than a single fault strand, ii) the development of a sinuous trace for a number of the major folds, whereby the ends are oriented normal to the compression direction, the centres parallel to the strike of the faults, iii) the development of a number of cross-folds, striking NNE-SSW and iv) the apparently recent development of a strike-slip component on at least one of the major thrust faults. The origin of the W-E, or WSW-ENE striking structures may be reactivation of Late Cretaceous faults, stratigraphic evidence for the existence of a "structural high" (the Hurunui High) over the majority of the area in the Late Cretaceous to Early Eocene times suggests the formation of a W-E trending horst structure, with a corresponding asymmetric graben to the south. The junction of WSW-ENE trending structures with N-S trending structures to the west centres on an alluvial-filled depression, Waikari Flat, into which the structures of the WSW-ENE trending imbricate thrust system plunge, locally curling to the SW at their ends to link with N-S trending structures to the south. Roof thrusting on two orientations, W-E and N-S, towards to SE is currently occurring above these structures. Currently the area is not highly seismically active, although a magnitude ~6.4 Ms earthquake in historic times has been recorded. The effects of tectonics on the drainage of the area does suggest that the majority of the systems, are still potentially active, albeit moving at a comparatively slow rate. The majority of the recent motion appears to be concentrated on the roof-thrusting occurring in Waikari Flat, and uplift along the Lowry Peaks Fault System. Increasing amounts of secondary movement on back-thrusts and cross fractures is also implied for western ends of the major imbricate thrust system. In contrast, the southern-most fault system appears to be largely sustaining dextral strike-slip motion, with some local folding in central portions.
The south Leader Fault (SLF) is a newly documented active structure that ruptured the surface during the Mw 7.8 Kaikoura earthquake. The Leader Fault is a NNE trending oblique left lateral thrust that links the predominantly right lateral ‘The Humps’ and Conway-Charwell faults. The present research uses LiDAR at 0.5 m resolution and field mapping to determine the factors controlling the surface geometries and kinematics of the south Leader Fault ruptures at the ground surface. The SLF zone is up to 2km wide and comprises a series of echelon NE-striking thrusts linked by near-vertical N-S striking faults. The thrusts are upthrown to the west by up to 1 m and dip 35-45°. Thrust slip surfaces are parallel with Cretaceous-Cenozoic bedding and may reflect flexural slip folding. By contrast, the northerly striking faults dip steeply (65° west- 85° east), and accommodate up to 3m of oblique left lateral displacement at the ground surface and displace Cenozoic bedding. Some of the SLF has been mapped in bedrock, although none were known to be active prior to the earthquake or have a strong topographic expression. The complexity of fault rupture and the width of the fault zone appears to reflect the occurrence of faulting and folding at the ground surface during the earthquake.
The magnitude 7.8 earthquake that struck North Canterbury, on the east coast of New Zealand’s South Island on 14 November 2016 had significant impacts and implications for the community of Kaikōura and surrounding settlements. The magnitude and scope of this event has resulted in extensive and ongoing geological and geophysical research into the event. The current paper complements this research by providing a review of existing social science research and offering new analysis of the impact of the earthquake and its aftermath on community resilience in Kaikōura over the past five years. Results demonstrate the significant economic implications for tourism, and primary industries. Recovery has been slow, and largely dependent on restoring transportation networks, which helped catalyse cooperation among local hospitality providers. Challenges remain, however, and not all sectors or households have benefited equally from post-quake opportunities, and long-term recovery trajectories continue to be hampered by COVID-19 pandemic. The multiple ongoing and future stressors faced by Kaikōura require integrated and equitable approaches in order to build capability and capacity for locally based development pathways to ensure long-term community resilience.
Numerous rockfalls released during the 2010–2011 Canterbury earthquake sequence affected vital road sections for local commuters. We quantified rockfall fatality risk on two main routes by adapting a risk approach for roads originally developed for snow avalanche risk. We present results of the collective and individual fatality risks for traffic flow and waiting traffic. Waiting traffic scenarios particularly address the critical spatial-temporal dynamics of risk, which should be acknowledged in operational risk management. Comparing our results with other risks commonly experienced in New Zealand indicates that local rockfall risk is close to tolerability thresholds and likely exceeds acceptable risk.
The structure and geomorphology of active orogens evolves on time scales ranging from a single earthquake to millions of years of tectonic deformation. Analysis of crustal deformation using new and established remote sensing techniques, and integration of these data with field mapping, geochronology and the sedimentary record, create new opportunities to understand orogenic evolution over these timescales. Timor Leste (East Timor) lies on the northern collisional boundary between continental crust from the Australian Plate and the Banda volcanic arc. GPS studies have indicated that the island of Timor is actively shortening. Field mapping and fault kinematic analysis of an emergent Pliocene marine sequence identifies gentle folding, overprinted by a predominance of NW-SE oriented dextral-normal faults and NE-SW oriented sinistral-normal faults that collectively bound large (5-20km2) bedrock massifs throughout the island. These fault systems intersect at non-Andersonian conjugate angles of approximately 120° and accommodate an estimated 20 km of orogen-parallel extension. Folding of Pliocene rocks in Timor may represent an early episode of contraction but the overall pattern of deformation is one of lateral crustal extrusion sub-parallel to the Banda Arc. Stratigraphic relationships suggest that extrusion began prior to 5.5 Ma, during and after initial uplift of the orogen. Sedimentological, geochemical and Nd isotope data indicate that the island of Timor was emergent and shedding terrigenous sediment into carbonate basins prior to 4.5 Ma. Synorogenic tectonic and sedimentary phases initiated almost synchronously across much of Timor Leste and <2 Myr before similar events in West Timor. An increase in plate coupling along this obliquely converging boundary, due to subduction of an outlying continental plateau at the Banda Trench, is proposed as a mechanism for uplift that accounts for orogen-parallel extension and early uplift of Timor Leste. Rapid bathymetric changes around Timor are likely to have played an important role in evolution of the Indonesian Seaway. The 2010 Mw 7.1 Darfield (Canterbury) earthquake in New Zealand was complex, involving multiple faults with strike-slip, reverse and normal displacements. Multi-temporal cadastral surveying and airborne light detection and ranging (LiDAR) surveys allowed surface deformation at the junction of three faults to be analyzed in this study in unprecedented detail. A nested, localized restraining stepover with contractional bulging was identified in an area with the overall fault structure of a releasing bend, highlighting the surface complexities that may develop in fault interaction zones during a single earthquake sequence. The earthquake also caused river avulsion and flooding in this area. Geomorphic investigations of these rivers prior to the earthquake identify plausible precursory patterns, including channel migration and narrowing. Comparison of the pre and post-earthquake geomorphology of the fault rupture also suggests that a subtle scarp or groove was present along much of the trace prior to the Darfield earthquake. Hydrogeology and well logs support a hypothesis of extended slip history and suggests that that the Selwyn River fan may be infilling a graben that has accumulated late Quaternary vertical slip of <30 m. Investigating fault behavior, geomorphic and sedimentary responses over a multitude of time-scales and at different study sites provides insights into fault interactions and orogenesis during single earthquakes and over millions of years of plate boundary deformation.
