A 3D high-resolution model of the geologic structure and associated seismic velocities in the Canterbury, New Zealand region is developed utilising data from depthconverted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model, developed using geostatistical Kriging, explicitly represents the significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age. The model is examined in the form of both geologic surface elevation contour maps as well as vertical cross sections of shear wave velocity, with the most prominent features being the Banks Peninsula Miocene-Pliocene volcanic edifice, and the Pegasus and Rakaia late Mesozoic-Neogene sedimentary basins. The adequacy of the modelled geologic surfaces is assessed through a residual analysis of point constraints used in the Kriging and qualitative comparisons with previous geologic models of subsets of the region. Seismic velocities for the lithological units between the geologic surfaces have also been derived, thus providing the necessary information for a Canterbury velocity model (CantVM) for use in physics-based seismic wave propagation. The developed model also has application for the determination of depths to specified shear wave velocities for use in empirical ground motion modelling, which is explicitly discussed via an example.
Active faults capable of generating highly damaging earthquakes may not cause surface rupture (i.e., blind faults) or cause surface ruptures that evade detection due to subsequent burial or erosion by surface processes. Fault populations and earthquake frequency-‐magnitude distributions adhere to power laws, implying that faults too small to cause surface rupture but large enough to cause localized strong ground shaking densely populate continental crust. The rupture of blind, previously undetected faults beneath Christchurch, New Zealand in a suite of earthquakes in 2010 and 2011, including the fatal 22 February 2011 moment magnitude (Mw) 6.2 Christchurch earthquake and other large aftershocks, caused a variety of environmental impacts, including major rockfall, severe liquefaction, and differential surface uplift and subsidence. All of these effects occurred where geologic evidence for penultimate effects of the same nature existed. To what extent could the geologic record have been used to infer the presence of proximal, blind and / or unidentified faults near Christchurch? In this instance, we argue that phenomena induced by high intensity shaking, such as rock fragmentation and rockfall, revealed the presence of proximal active faults in the Christchurch area prior to the recent earthquake sequence. Development of robust earthquake shaking proxy datasets should become a higher scientific priority, particularly in populated regions.
Deep shear wave velocity (Vs) profiles (>400 m) were developed at 14 sites throughout Christchurch, New Zealand using surface wave methods. This paper focuses on the inversion of surface wave data collected at one of these sites, Hagley Park. This site is located on the deep soils of the Canterbury Plains, which consist of alluvial gravels inter-bedded with estuarine and marine sands, silts, clays and peats. Consequently, significant velocity contrasts exist at the interface between geologic formations. In order to develop realistic velocity models in this complex geologic environment, a-priori geotechnical and geologic data were used to identify the boundaries between geologic formations. This information aided in developing the layering for the inversion parameters. Moreover, empirical reference Vs profiles based on material type and confining pressure were used to develop realistic Vs ranges for each layer. Both the a-priori layering information and the reference Vs curves proved to be instrumental in generating realistic velocity models that account for the complex inter-bedded geology in the Canterbury Plains.
A zone of active tectonism occurs in mid and north Canterbury, from the Rakaia to the Waipara Rivers, which coincides with seismicity concentrations and several Quaternary surface anomalies and is here defined as the Porters Pass Tectonic Zone. Although parallel to the Marlborough faults to the north, the lack of regional definition suggests this zone is much younger in its inception reflecting a southward movement of the plate rotation vector. The objectives of this study were to map the structures associated with this zone in the segment between the Rakaia and Waimakariri Rivers with detailed analysis concentrated in the upper Kawai Valley. Quaternary offsets on the main lineament of the Porters Pass Fault were traced through the area and evidence for the rate of movement, probable magnitudes and return periods of related seismic events was sought. The basement was found to be complicated by pre-existing deformation structures in Torlesse Group rocks which have been subsequently been re-activated or rotated by recent fault movement probably beginning in the Pleistocene. This phase is dominantly thrusting and uplift has lead to the erosion of most of the overlying sedimentary cover. Remnants of the Cret-Tertiary sediments still remain as fault-bounded packets. Evidence suggests that a change to development of a regional lateral shear associated with the Porters Pass Tectonic Zone transects the thrust system with complex interaction between the older reverse and new strike-slip faults. Offset rates along the segments of the Porters Pass Fault are not well constrained but are believed to be approximately in the range of 11-13 mm/year for at least the last 130,000 years. This rate is similar to other large faults in the Marlborough region. Two earthquake events have been identified and dated at 600 and 2000 years ago, with a magnitude of greater than 6.5. Evidence suggests characteristic earthquakes along the Porters Pass Fault are greater than Magnitude 7. This result has some major ramifications for the expected seismic hazards for nearby Christchurch.
The 22 February 2011, Mw6.2 Christchurch earthquake is the most costly earthquake to affect New Zealand, causing an estimated 181 fatalities and severely damaging thousands of residential and commercial buildings. This paper presents a summary of some of the observations made by the NSF-sponsored GEER Team regarding the geotechnical/geologic aspects of this earthquake. The Team focused on documenting the occurrence and severity of liquefaction and lateral spreading, performance of building and bridge foundations, buried pipelines and levees, and significant rockfalls and landslides. Liquefaction was pervasive and caused extensive damage to residential properties, water and wastewater networks, high-rise buildings, and bridges. Entire neighborhoods subsided, resulting in flooding that caused further damage. Additionally, liquefaction and lateral spreading resulted in damage to bridges and to stretches of levees along the Waimakariri and Kaiapoi Rivers. Rockfalls and landslides in the Port Hills damaged several homes and caused several fatalities.
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.
Rock mass defect controlled deep-seated landslides are widespread within the deeply incised landscapes formed in Tertiary soft rock terrain in New Zealand. The basal failure surfaces of deep-seated slope failures are defined by thin, comparatively weak and laterally continuous bedding parallel layers termed critical stratigraphic horizons. These horizons have a sedimentary origin and have typically experienced some prior tectonically induced shear displacement at the time of slope failure. The key controls on the occurrence and form of deep-seated landslides are considered in terms of rock mass defect properties and tectonic and climatic forcing. The selection of two representative catchments (in southern Hawke's Bay and North Canterbury) affected by tectonic and climatic forcing has shown that the spatial and temporal initiation of deep-seated bedrock landslides in New Zealand Tertiary soft rock terrain is a predictable rather than a stochastic process; and that deep-seated landslides as a mass wasting process have a controlling role in landscape evolution in many catchments formed in Tertiary soft rock terrain. The Ella Landslide in North Canterbury is a deep-seated (~85 m) translational block slide that has failed on a 5 - 10 mm thick, kaolinite-rich, pre-sheared critical stratigraphic horizon. The residual strength of this sedimentary horizon, (C'R 2.6 - 2.7 kPa, and Ѳ'R = 16 - 21°), compared to the peak strength of the dominant lithology (C' = 176 kPa, and Ѳ' = 37°) defines a high strength contrast in the succession, and therefore a critical location for the basal failure surface of deep-seated slope failures. The (early to mid Holocene) Ella Landslide debris formed a large landslide dam in the Kate Stream catchment and this has significantly retarded rates of mass wasting in the middle catchment. Numerical stability analysis shows that this slope failure would have most likely required the influence of earthquake induced strong ground motion and the event is tentatively correlated to a Holocene event on the Omihi Fault. The influence of this slope failure is likely to affect the geomorphic development of the catchment on a scale of 10⁴ - 10⁵ years. In deeply incised catchments at the southeastern margin of the Maraetotara Plateau, southern Hawke's Bay, numerous widespread deep-seated landslides have basal failure surfaces defined by critical stratigraphic horizons in the form of thin « 20 mm) tuffaceous beds in the Makara Formation flysch (alternating sandstone and mudstone units). The geometry of deep-seated slope failures is controlled by these regularly spaced (~70 m), very weak critical stratigraphic horizons (C'R 3.8 - 14.2 kPa, and Ѳ'R = 2 - 5°), and regularly spaced (~45 m) and steeply dipping (-50°) critical conjugate joint/fault sets, which act as slide block release surfaces. Numerical stability analysis and historical precedent show that the temporal initiation of deep-seated landslides is directly controlled by short term tectonic forcing in the form of periodic large magnitude earthquakes. Published seismic hazard data shows the recurrence interval of earthquakes producing strong ground motions of 0.35g at the study site is every 150 yrs, however, if subduction thrust events are considered the level of strong ground motion may be much higher. Multiple occurrences of deep-seated slope failure are correlated to failure on the same critical stratigraphic horizon, in some cases in three adjacent catchments. Failure on multiple critical stratigraphic horizons leads to the development of a "stepped" landscape morphology. This slope form will be maintained during successive accelerated stream incision events (controlled by long term tectonic and climatic forcing) for as long as catchments are developing in this specific succession. Rock mass defect controlled deep seated landslides are controlling catchment head progression, landscape evolution and hillslope morphology in the Hawke's Bay study area and this has significant implications for the development of numerical landscape evolution models of landscapes formed in similar strata. Whereas the only known numerical model to consider deep seated landslides as an erosion process (ZSCAPE) considers them as stochastic in time and space, this study shows that this could not be applied to a landscape where the widespread spatial occurrence of deep-seated landslides is controlled by rock mass defects. In both of the study areas for this project, and by implication in many catchments in Tertiary soft rock terrain, deep-seated landslides controlled by rock mass defect strength, spacing and orientation, and tectonic and climatic forcing have an underlying control on landscape evolution. This study quantifies parameters for the development of numerical landscape evolution models that would assess the role of specific parameters, such as uplift rates, incision rates and earthquake recurrence in catchment evolution in Tertiary soft rock terrain.
Oblique convergence of the Pacific and Australian Plates is accommodated in the northern South Island by the Marlborough Fault System. The Hope Fault is the southern of four major dextral strike-slip faults of this system. Hanmer Basin is a probable segment boundary between the Hope River and Conway segments of the Hope Fault. The Conway segment is transpressional and shows increasing structural complexity near the segment boundary at Hanmer Basin, with multiple Late Quaternary traces, and fault-parallel folding in response to across-fault shortening. Between Hossack Station and Hanmer Basin a crush zone in excess of one kilometre wide is exposed in incised streams and rivers. The crush zone has an asymmetrical geometry about the active trace of the Hope Fault, being only 100-300 metres wide south of the fault, and more than 500 metres wide north of the fault. The most intense deformation of Torlesse bedrock occurs at the south side of the fault zone, indicating that strain is accommodated against the fault footwall. North of the fault deformation is less intense, but occurs over a wider area. The wide fault zone at Hossack Station may reflect divergence of the Hanmer Fault, a major splay of the Hope Fault. At Hossack Station, the Hope Fault has accommodated at least 260 metres of dextral displacement during the Holocene. Dating of abandoned stream channels, offset by the Hope Fault, indicated a Late Holocene dextral slip-rate of 18±8 mm-¹ for the west end of the Conway segment. Using empirical formulae and inferred fault parameters, the expected magnitude of an earthquake generated by the Conway segment is M6.9 to M7.4; for an exceedence probability of 10%, the magnitude is M7.7 to M7.9. Effects associated with coseismic rupture of the Conway segment include shaking of up to MMIX along the ruptured fault and at Hanmer Basin. Uplift at the east end of Hanmer Basin, in conjunction with subsidence at the southwest margin of the basin, is resulting in the development of onlapping stratigraphy. Seismic reflection profiles support this theory. Possible along-fault migration of the basin is inferred to be a consequence of non-parallelism of the master faults.
The Amuri Earthquake of September 1, 1888 (magnitude M = 6.5 to 6.8) occurred on the Hope River Segment of the Hope Fault west of Hanmer Plains. The earthquake was felt strongly in North Canterbury and North Westland and caused considerable property damage and landsliding in the Lower Hope Valley. However, damage reports and the spatial distribution of felt intensities emphasize extreme variations in seismic effects over short distances, probably due to topographic focusing and local ground conditions. Significant variations in lateral fault displacement occurred at secondary fault segment boundaries (side-steps and bends in the fault trace) during the 1888 earthquake. This historical spatial variation in lateral slip is matched by the Late Quaternary geomorphic distribution of slip on the Hope River Segment of the Hope Fault. Trenching studies at two sites on the Hope Fault have also identified evidence for five pre-historic earthquakes of similar magnitude to the 1888 earthquake and an average recurrence interval of 134 ± 27 years between events. Magnitude estimates for the 1888 earthquake are combined with a. strong ground motion attenuation expression to provide an estimate of potential ground accelerations in Amuri District during-future earthquakes on the Hope River Segment of the Hope Fault. The predicted acceleration response on bedrock sites within 20 km of the epicentral region is between 0.23 g and 0.34 g. The close match between the historic, inferred pre-historic and geomorphic distribution of lateral slip indicates that secondary fault segmentation exerts a strong structural control on rupture propagation and the expression of fault displacement at the surface. In basement rocks at depth the spatial variations in slip are inferred to be distributed within zones of pervasive cataclastic shear, on either side of the fault segment boundaries. The large variations in surface displacement across fault segment boundaries means that one must know the geometry of the fault in order to evaluate slip-rates calculated from individual locations. The average Late Quaternary slip-rate on the Hope Fault at Glynn Wye Station is between 15.5 mm/yr and 18.25 mm/yr and the rate on the subsidiary Kakapo Fault is between 5.0 mm/yr and 7.5 mm/yr. These rates have been determined from sites which are relatively free of structural complication.
This paper provides a photographic tour of the ground-surface rupture features of the Greendale Fault, formed during the 4th September 2010 Darfield Earthquake. The fault, previously unknown, produced at least 29.5 km of strike-slip surface deformation of right-lateral (dextral) sense. Deformation, spread over a zone between 30 and 300 m wide, consisted mostly of horizontal flexure with subsidiary discrete shears, the latter only prominent where overall displacement across the zone exceeded about 1.5 m. A remarkable feature of this event was its location in an intensively farmed landscape, where a multitude of straight markers, such as fences, roads and ditches, allowed precise measurements of offsets, and permitted well-defined limits to be placed on the length and widths of the surface rupture deformation.
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.
Liquefaction-induced lateral spreading in Christchurch and surrounding suburbs during the recent Canterbury Earthquake Sequence (2010-2011) caused significant damage to structures and lifelines located in close proximity to streams and rivers. Simplified methods used in current engineering practice for predicting lateral ground displacements exhibit a high degree of epistemic uncertainty, but provide ‘order of magnitude’ estimates to appraise the hazard. We wish to compare model predictions to field measurements in order to assess the model’s capabilities and limitations with respect to Christchurch conditions. The analysis presented focuses on the widely-used empirical model of Youd et al. (2002), developed based on multi-linear regression (MLR) of case history data from lateral spreading occurrence in Japan and the US. Two issues arising from the application of this model to Christchurch were considered: • Small data set of Standard Penetration Test (SPT) and soil gradation indices (fines content FC, and mean grain size, D50) required for input. We attempt to use widely available CPT data with site specific correlations to FC and D50. • Uncertainty associated with the model input parameters and their influence on predicted displacements. This has been investigated for a specific location through a sensitivity analysis.
Surface rupture of the previously unrecognised Greendale Fault extended west-east for ~30 km across alluvial plains west of Christchurch, New Zealand, during the Mw 7.1 Darfield (Canterbury) earthquake of September 2010. Surface rupture displacement was predominantly dextral strike-slip, averaging ~2.5 m, with maxima of ~5 m. Vertical displacement was generally less than 0.75 m. The surface rupture deformation zone ranged in width from ~30 to 300 m, and comprised discrete shears, localised bulges and, primarily, horizontal dextral flexure. About a dozen buildings, mainly single-storey houses and farm sheds, were affected by surface rupture, but none collapsed, largely because most of the buildings were relatively flexible and resilient timber-framed structures and also because deformation was distributed over a relatively wide zone. There were, however, notable differences in the respective performances of the buildings. Houses with only lightly-reinforced concrete slab foundations suffered moderate to severe structural and non-structural damage. Three other buildings performed more favourably: one had a robust concrete slab foundation, another had a shallow-seated pile foundation that isolated ground deformation from the superstructure, and the third had a structural system that enabled the house to tilt and rotate as a rigid body. Roads, power lines, underground pipes, and fences were also deformed by surface fault rupture and suffered damage commensurate with the type of feature, its orientation to the fault, and the amount, sense and width of surface rupture deformation.
This project was initiated by ENGEO Limited and KiwiRail Holdings Limited to assess the stability of Slovens Creek Viaduct (specifically its western abutment) and a 3km section of rail corridor between Slovens Creek Viaduct and Avoca on the Midland Line (MDL). Commonly known as the scenic TranzAlpine rail journey (through Arthurs Pass National Park) the MDL connects Greymouth to Christchurch via Rolleston, where the MDL meets the Main South Line into Christchurch. The project area is approximately 40km southeast of Arthurs Pass Township, in the eastern extension of the Castle Hill Basin which is part of the Waimakariri Catchment and Canterbury Foothills. The field area is underlain by Rakaia Terrane, which is part of the Torlesse Composite Terrane forming the basement rock unit for the field area. Cretaceous-Tertiary rocks of the Castle Hill Basin overlie the basement strata and record a transgression-regression sequence, as well as mid-Oligocene submarine volcanism. The stratigraphic sequence in the Castle Hill Basin, and its eastern extension to Avoca, comprises two formations of the Eyre group, the older Broken River Formation and the younger Iron Creek Formation. Deep marine Porter Group limestones, marls, and tuffs of Oligocene age succeed the Iron Creek Formation of the Eyre Group, and probably records the maximum of the transgression. The Enys Formation lies disconformably on the Porter Group and is overlain unconformably by Late Pleistocene glacifluvial and glacial deposits. The Tertiary strata in the Slovens-Avoca rail corridor are weak, and the clay-rich tuff derived from mid-Oligocene volcanism is particularly prone to slaking. Extensive mapping carried out for this project has identified that some 90 percent of the surface along the length of the Slovens-Avoca corridor has been subject to mass movement. The landslides of the Slovens-Avoca rail corridor are clearly younger than the Last Glaciation, and Slovens Creek has been downcutting, with associated faulting and uplift, to form the present day geomorphology of the rail corridor. Deep-seated landslides in the rail corridor extend to Slovens Creek, locally deflecting the stream course, and a generic ground failure model for the rail corridor has been developed. Exploratory geotechnical investigations, including core drilling, installation of an inclinometer and a piezometer, enabled the construction of a simple ground model and cross section for the Slovens Creek Viaduct western abutment. Limit-equilibrium and pseudo-static slope stability analyses using both circular and block critical slip surface search methods were applied to the ground model for the western abutment of Slovens Creek Viaduct. Piezometric and strength data obtained during laboratory testing of core material have been used to constrain the western abutment stability assessment for one representative section line (C-C’). Prior to pseudo-static sensitivity analyses peak ground acceleration (PGA) for various Ultimate Limit State (ULS) design return periods, defined by an equation given in NZS1170.5:2004, were calculated and have been used as a calibration technique to find and compare specific PGA values for pseudo-static analyses in the Slovens Creek Viaduct area. The main purpose has been to provide an indication of how railway infrastructure could be affected by seismic events of various return periods defined by ULS design standards for the area. Limit equilibrium circular slip surface search methods, both grid search and auto refine search, indicated the slope is stable with a FoS greater than 1.0 returned from each, although one particular surface returned the lowest FoS in each. This surface is in the lower portion of the slope, adjacent to Slovens Stream and northeast of the MDL. As expected, pseudo-static analyses returned a lower FoS overall when compared to limit equilibrium analyses. The PGA analyses suggest that partial ground failure at the Slovens Creek Viaduct western abutment could occur in a 1 in 25-year return period event within materials on the slower slope beyond the immediate rail corridor. A ULS (1 in 500-year) event in the Slovens Creek Viaduct area would likely produce a PGA of ~0.9g, and the effects on the western abutment and rail infrastructure would most likely be catastrophic. Observed ground conditions for the western abutment of the Slovens Creek Viaduct suggest there is no movement within the landslide at depth within the monitoring timeframe of this project (22 May 2015 – 4 August 2015). Slope stability monitoring is recommended to be continued in two parts: (1) the inclinometer in BH1 is to be monitored on a six monthly basis for one year following completion of this thesis, and then annually unless ground movements become evident; and (2) surface movement monitoring should be installed using a fixed datum on the stable eastern abutment. Long-term stability management strategies for the Slovens Creek Viaduct western abutment are dependent upon future observed changes and ongoing monitoring. Hazard and risk assessment using the KiwiRail Qualitative Risk Assessment Framework (QRA) is recommended, and if slope stability becomes problematic for operation of the Midland Line consideration should be given to deep slope drainage. In the event of a large magnitude or high PGA earthquake all monitoring should be reviewed.
None
None
None
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 Canterbury Region is susceptible to a variety of natural hazards, including earthquakes, landslides and climate hazards. Increasing population and tourism within the region is driving development pressures and as more and more development occurs, the risk from natural hazards increases. In order to avoid development occurring in unacceptably vulnerable locations, natural hazard assessments are required. This study is a reconnaissance natural hazard assessment of Lakes Lyndon, Coleridge and Tekapo. There is restricted potential for development at Lake Lyndon, because the land surrounding the lake is owned by the Crown and has a number of development restrictions. However, there is the potential for conservation or recreation-linked development to occur. There is more potential for development at Lake Coleridge. Most of the land surrounding the lake is privately owned and has less development restrictions. The majority of land surrounding Lake Tekapo is divided into Crown-owned pastoral leases, which are protected from development, such as subdivision. However, there are substantial areas around the lake, which are privately owned and, therefore, have potential for development. Earthquake, landslide and climate hazards are the main natural hazards threatening Lakes Lyndon, Coleridge and Tekapo. The lakes are situated in a zone of active earth deformation in which large and relatively frequent earthquakes are produced. A large number of active faults lie within 15 km of each lake, which are capable of producing M7 or larger earthquakes. Ground shaking, liquefaction, landslides, tsunami and seiches are among the consequences of earthquakes, all of which have the potential to cause severe damage to lives, lifelines and infrastructure. Landslides are also common in the landscape surrounding the lakes. The majority of slopes surrounding the lakes are at significant risk from earthquake-induced failure under moderate to strong earthquake shaking. This level of shaking is expected to occur in any 50 year period around Lakes Lyndon and Coleridge, and in any 150 year period around Lake Tekapo. Injuries, fatalities and property damage can occur directly from landslide impact or from indirect effects such as flooding from landslide-generated tsunami or from landslide dam outbreaks. Lakes Lyndon, Coleridge and Tekapo are also susceptible to climate hazards, such as high winds, drought, heavy snowfall and heavy rainfall, which can lead to landslides and flooding. Future climate change due to global warming is most likely going to affect patterns of frequency and magnitudes of extreme weather events, leading to an increase in climate hazards. Before development is permitted around the lakes, it is essential that each of these hazards is considered so that unacceptably vulnerable areas can be avoided.
Bulk rock strength is greatly dependent on fracture density, so that reductions in rock strength associated with faulting and fracturing should be reflected by reduced shear coupling and hence S-wave velocity. This study is carried out along the Canterbury rangefront and in Otago. Both lie within the broader plate boundary deformation zone in the South Island of New Zealand. Therefore built structures are often, , located in areas where there are undetected or poorly defined faults with associated rock strength reduction. Where structures are sited near to, or across, such faults or fault-zones, they may sustain both shaking and ground deformation damage during an earthquake. Within this zone, management of seismic hazards needs to be based on accurate identification of the potential fault damage zone including the likely width of off-plane deformation. Lateral S-wave velocity variability provides one method of imaging and locating damage zones and off-plane deformation. This research demonstrates the utility of Multi-Channel Analysis of Surface Waves (MASW) to aid land-use planning in such fault-prone settings. Fundamentally, MASW uses surface wave dispersive characteristics to model a near surface profile of S-wave velocity variability as a proxy for bulk rock strength. The technique can aid fault-zone planning not only by locating and defining the extent of fault-zones, but also by defining within-zone variability that is readily correlated with measurable rock properties applicable to both foundation design and the distribution of surface deformation. The calibration sites presented here have well defined field relationships and known fault-zone exposure close to potential MASW survey sites. They were selected to represent a range of progressively softer lithologies from intact and fractured Torlesse Group basement hard rock (Dalethorpe) through softer Tertiary cover sediments (Boby’s Creek) and Quaternary gravels. This facilitated initial calibration of fracture intensity at a high-velocity-contrast site followed by exploration of the limits of shear zone resolution at lower velocity contrasts. Site models were constructed in AutoCAD in order to demonstrate spatial correlations between S-wave velocity and fault zone features. Site geology was incorporated in the models, along with geomorphology, river profiles, scanline locations and crosshole velocity measurement locations. Spatial data were recorded using a total-station survey. The interpreted MASW survey results are presented as two dimensional snapshot cross-sections of the three dimensional calibration-site models. These show strong correlations between MASW survey velocities and site geology, geomorphology, fluvial profiles and geotechnical parameters and observations. Correlations are particularly pronounced where high velocity contrasts exist, whilst weaker correlations are demonstrated in softer lithologies. Geomorphic correlations suggest that off-plane deformation can be imaged and interpreted in the presence of suitable topographic survey data. A promising new approach to in situ and laboratory soft-rock material and mass characterisation is also presented using a Ramset nail gun. Geotechnical investigations typically involve outcrop and laboratory scale determination of rock mass and material properties such as fracture density and unconfined compressive strength (UCS). This multi-scale approach is espoused by this study, with geotechnical and S-wave velocity data presented at multiple scales, from survey scale sonic velocity measurements, through outcrop scale scanline and crosshole sonic velocity measurements to laboratory scale property determination and sonic velocity measurements. S-wave velocities invariably increased with decreasing scale. These scaling relationships and strategies for dealing with them are investigated and presented. Finally, the MASW technique is applied to a concealed fault on the Taieri Ridge in Macraes Flat, Central Otago. Here, high velocity Otago Schist is faulted against low velocity sheared Tertiary and Quaternary sediments. This site highlights the structural sensitivity of the technique by apparently constraining the location of the principal fault, which had been ambiguous after standard processing of the seismic reflection data. Processing of the Taieri Ridge dataset has further led to the proposal of a novel surface wave imaging technique termed Swept Frequency Imaging (SFI). This inchoate technique apparently images the detailed structure of the fault-zone, and is in agreement with the conventionally-determined fault location and an existing partial trench. Overall, the results are promising and are expected to be supported by further trenching in the near future.
Among the deformation features produced in Christchurch by the September 4th Darfield Earthquake were numerous and widespread “sand volcanoes”. Most of these structures occurred in urban settings and “erupted” through a hardened surface of concrete or tarseal, or soil. Sand volcanoes were also widespread in the Avon‐ Heathcote Estuary and offered an excellent opportunity to readily examine shallow subsurface profiles and as such the potential appearance of such structures in the rock record.
None
This study contains an evaluation of the seismic hazard associated with the Springbank Fault, a blind structure discovered in 1998 close to Christchurch. The assessment of the seismic hazard is approached as a deterministic process in which it is necessary to establish: 1) fault characteristics; 2) the maximum earthquake that the fault is capable of producing and 3) ground motions estimations. Due to the blind nature of the fault, conventional techniques used to establish the basic fault characteristics for seismic hazard assessments could not be applied. Alternative methods are used including global positioning system (GPS) surveys, morphometric analyses along rivers, shallow seismic reflection surveys and computer modelling. These were supplemented by using multiple empirical equations relating fault attributes to earthquake magnitude, and attenuation relationships to estimate ground motions in the near-fault zone. The analyses indicated that the Springbank Fault is a reverse structure located approximately 30 km to the northwest of Christchurch, along a strike length of approximately 16 km between the Eyre and Ashley River. The fault does not reach the surface, buy it is associated with a broad anticline whose maximum topographic expression offers close to the mid-length of the fault. Two other reverse faults, the Eyrewell and Sefton Faults, are inferred in the study area. These faults, together with the Springbank and Hororata Faults and interpreted as part of a sys of trust/reverse faults propagating from a decollement located at mid-crustal depths of approximately 14 km beneath the Canterbury Plains Within this fault system, the Springbank Fault is considered to behave in a seismically independent way, with a fault slip rate of ~0.2 mm/yr, and the capacity of producing a reverse-slip earthquake of moment magnitude ~6.4, with an earthquake recurrence of 3,000 years. An earthquake of the above characteristics represents a significant seismic hazard for various urban centres in the near-fault zone including Christchurch, Rangiora, Oxford, Amberley, Kaiapoi, Darfield, Rollestion and Cust. Estimated peak ground accelerations for these towns range between 0.14 g to 0.5 g.
At 4.35am on Saturday 4 September 2010, a magnitude 7.1 earthquake struck near the township of Darfield in Canterbury leading to widespread damage in Christchurch and the wider central Canterbury region. Though it was reported no lives were lost, that was not entirely correct. Over 3,000 animals perished as a result of the earthquake and 99% of these deaths would have been avoidable if appropriate mitigation measures had been in place. Deaths were predominantly due to zoological vulnerability of birds in captive production farms. Other problems included lack of provision of animal welfare at evacuation centres, issues associated with multiple lost and found pet services, evacuation failure due to pet separation and stress impact on dairy herds and associated milk production. The Canterbury Earthquake has highlighted concerns over a lack of animal emergency welfare planning and capacity in New Zealand, an issue that is being progressed by the National Animal Welfare Emergency Management Group. As animal emergency management becomes better understood by emergency management and veterinary professionals, it is more likely that both sectors will have greater demands placed upon them by national guidelines and community expectations to ensure provisions are made to afford protection of animals in times of disaster. A subsequent and more devastating earthquake struck the region on Monday 22 February 2011; this article however is primarily focused on the events pertaining to the September 4 event.
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 Porters Pass fault (PPF) is a prominent element of the Porters Pass-Amberley Fault Zone (PPAFZ) which forms a broad zone of active earth deformation ca 100 km long, 60-90 km west and north of Christchurch. For a distance of ca 40 km the PPF is defined by a series of discontinuous Holocene active traces between the Rakaia and Waimakariri Rivers. The amount of slip/event and the timing of paleoearthquakes are crucial components needed to estimate the earthquake potential of a fault. Movement was assumed to be, coseismic and was quantified by measuring displaced geomorphic features using either tape measure or surveying equipment. Clustering of offset data suggests that four to five earthquakes occurred on the PPF during the Holocene and these range between ca 5-7 m/event. Timing information was obtained from four trenches excavated across the fault and an auger adjacent to the fault. Organic samples from these sites were radiocarbon dated and used in conjunction with data from previous studies to identify the occurrence of at least four earthquakes at 8500 ± 200, 5300 ± 700, 2500 ± 200 and 1000 ± 100 years B.P. Evidence suggests that an additional event is also possible at 6200 ± 500 years B.P. The ~1000, 5300 and 6200 years B.P. paleoearthquakes were previously unrecognised, while the 500 year event previously inferred from rock-avalanche data has been discarded. The present data set produces recurrence intervals of ~2000-2500 years for the Holocene. The identification of only one Holocene PPF rupture to the west of Red Lakes indicates the presence of a segment boundary that prevents the propagation of rupture beyond this point. This is consistent with displacement data and results in slip rates of 0.5-0.7 mm/yr and 2.5-3.4 mm/yr to the west and east of Red Lakes respectively. It is possible that the nearby extensional Red Hill Fault influences PPF rupture propagation. The combination of geometric, slip rate and timing data has enabled the magnitude of prehistoric earthquakes on the PPF to be estimated. These magnitudes range from an average of between 6.9 for a fault rupture from Waimakariri River to Red Lakes, to a maximum of 7.4 that ruptures the entire length of the PPAFZ, including the full length of the PPF. These estimates are approximately consistent with previous magnitude estimates along the full length of the PPAFZ of between 7.0 and 7.5.
Despite over a century of study, the relationship between lunar cycles and earthquakes remains controversial and difficult to quantitatively investigate. Perhaps as a consequence, major earthquakes around the globe are frequently followed by 'prediction' claims, using lunar cycles, that generate media furore and pressure scientists to provide resolute answers. The 2010-2011 Canterbury earthquakes in New Zealand were no exception; significant media attention was given to lunarderived earthquake predictions by non-scientists, even though the predictions were merely 'opinions' and were not based on any statistically robust temporal or causal relationships. This thesis provides a framework for studying lunisolar earthquake temporal relationships by developing replicable statistical methodology based on peer reviewed literature. Notable in the methodology is a high accuracy ephemeris, called ECLPSE, designed specifically by the author for use on earthquake catalogs, and a model for performing phase angle analysis. The statistical tests were carried out on two 'declustered' seismic catalogs, one containing the aftershocks from the Mw7.1 earthquake in Canterbury, and the other containing Australian seismicity from the past two decades. Australia is an intraplate setting far removed from active plate boundaries and Canterbury is proximal to a plate boundary, thus allowing for comparison based on tectonic regime and corresponding tectonic loading rate. No strong, conclusive, statistical correlations were found at any level of the earthquake catalogs, looking at large events, onshore events, offshore events, and the fault type of some events. This was concluded using Schuster's test of significance with α=5% and analysis of standard deviations. A few weak correlations, with p-5-10% of rejecting the null hypothesis, and anomalous standard deviations were found, but these are difficult to interpret. The results invalidate the statistical robustness of 'earthquake predictions' using lunisolar parameters in this instance. An ambitious researcher could improve on the quality of the results and on the range of parameters analyzed. The conclusions of the thesis raise more questions than answers, but the thesis provides an adaptable methodology that can be used to further investigation the problem.
The Avon and Heathcote Rivers, located in the city of Christchurch, New Zealand, are lowland spring-fed rivers linked with the Christchurch Groundwater System. At present, the flow paths and recharge sources to the Christchurch Groundwater System are not fully understood. Study of both the Avon and Heathcote Rivers can provide greater insight into this system. In addition, during the period 2010-2012, Christchurch has experienced large amounts of seismic activity, including a devastating Mw 6.2 aftershock on February 22nd, 2011, which caused widespread damage and loss of life. Associated with these earthquakes was the release of large amounts of water through liquefaction and temporary springs throughout the city. This provided a unique opportunity to study groundwater surface water interactions following a large scale seismic event. Presented herein is the first major geochemical study on the Avon and Heathcote Rivers and the hydrological impact of the February 22, 2011 Christchurch Earthquake. The Avon, Heathcote, and Waimakariri Rivers were sampled in quarterly periods starting in July 2011 and analyzed for stable Isotopes δ¹⁸O, δD, and δ¹³C and major anion composition. In addition, post -earthquake samples were collected over the days immediately following the February 22, 2011 earthquake and analyzed for stable isotopes δ¹⁸O and δD and major anion composition. A variety of analytical methods were used identify the source of the waters in the Avon-Heathcote System and evaluate the effectiveness of stable isotopes as geochemical tracers in the Christchurch Groundwater System. The results of this thesis found that the waters from the Avon and Heathcote Rivers are geochemically the same, originating from groundwater, and exhibit a strong tidal influence within 5km of the Avon-Heathcote Estuary. The surface waters released following the February 22nd, 2011 earthquake were indistinguishable from quarterly samples taken from the Avon and Heathcote Rivers when comparing stable isotopic composition. The anion data suggests the waters released following the February 22nd, 2011 Christchurch Earthquake were sourced primarily from shallow groundwater, and also suggests a presence of urban sewage at some sites. Attempts to estimate recharge sources for the Avon-Heathcote Rivers using published models for the Christchurch Groundwater System yielded results that were not consistent between models. In evaluating the use of geochemical constituents as tracers in the Christchurch Groundwater System, no one isotope could provide a clear resolution, but when used in conjunction, δ¹⁸O, δ¹³C, and DIC, seem to be the most effective tracers. Sample sizes for δ¹³C were too small for a robust evaluation. Variability on the Waimakariri River appears to be greater than previously estimated, which could have significant impacts on geochemical models for the Christchurch Groundwater System. This research demonstrates the value of using multiple geochemical constituents to enrich our understanding of the groundwater surfaces-water interactions and the Christchurch Groundwater System as a whole.
The Mw 7.1 Darfield earthquake generated a ~30 km long surface rupture on the Greendale Fault and significant surface deformation related to related blind faults on a previously unrecognized fault system beneath the Canterbury Plains. This earthquake provided the opportunity for research into the patterns and mechanisms of co-seismic and post-seismic crustal deformation. In this thesis I use multiple across-fault EDM surveys, logic trees, surface investigations and deformation feature mapping, seismic reflection surveying, and survey mark (cadastral) re-occupation using GPS to quantify surface displacements at a variety of temporal and spatial scales. My field mapping investigations identified shaking and crustal displacement-induced surface deformation features south and southwest of Christchurch and in the vicinity of the projected surface traces of the Hororata Blind and Charing Cross Faults. The data are consistent with the high peak ground accelerations and broad surface warping due to underlying reverse faulting on the Hororata Blind Fault and Charing Cross Fault. I measured varying amounts of post-seismic displacement at four of five locations that crossed the Greendale Fault. None of the data showed evidence for localized dextral creep on the Greendale Fault surface trace, consistent with other studies showing only minimal regional post-seismic deformation. Instead, the post-seismic deformation field suggests an apparent westward translation of northern parts of the across-fault surveys relative to the southern parts of the surveys that I attribute to post-mainshock creep on blind thrusts and/or other unidentified structures. The seismic surveys identified a deformation zone in the gravels that we attribute to the Hororata Blind Fault but the Charing Cross fault was not able to be identified on the survey. Cadastral re-surveys indicate a deformation field consistent with previously published geodetic data. We use this deformation with regional strain rates to estimate earthquake recurrence intervals of ~7000 to > 14,000 yrs on the Hororata Blind and Charing Cross Faults.