This paper summarizes the development of a high-resolution surficial shear wave velocity model based on the combination of the large high-spatial-density database of cone penetration test (CPT) logs in and around Christchurch, New Zealand and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. This near-surface shear wave velocity model has applications for site characterization efforts via the development of maps of time-averaged shear wave velocities over specific depths, as well as use in site response analysis and ground motion simulation.
This paper presents a seismic velocity model of Canterbury, New Zealand based on 3D geologic surfaces and velocities from a range of data sources. The model provides the 3D crustal structure in the region at multiple length scales for seismic wave propagation simulations, such as broadband ground motion and shallow site response analyses related to understanding the ground motions and site responses during the 2010- 2011 Canterbury earthquakes. Pre-Quaternary geologic horizons are calculated based on the reinterpretation of a comprehensive network of seismic reflection surveys from seven different campaigns over the past 50 years, as well as point constraints across an array of petroleum industry drill holes. Particular attention is given to a detailed representation of Quaternary stratigraphy, representing shallow (z<250m) near-surface layers in the model. Seismic velocities are obtained from seismic reflection processing (for Vp) and also recently performed active and passive surface wave analyses (for Vs). Over 1,700 water wells in the region are used to constrain the complex inter-bedded Quaternary stratigraphy (gravels, sands, silts, organics etc.) near the coastline, including beneath urban Christchurch, which has resulted from fluvial deposition and marine regression and transgression. For the near-surface Springston and Christchurch Formations in the Christchurch urban area (z<50m), high-spatial resolution seismic velocities (including Vs30 ) were obtained from over 13,000 cone penetration tests combined with a recently developed CPT-Vs correlation.
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.
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.
This paper summarizes the development of a region-wide surficial shear wave velocity model based on the combination of the large high-spatial-density database of cone penetration test (CPT) logs in and around Christchurch, New Zealand and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. The ongoing development of this near-surface shear wave velocity model has applications for site characterization efforts via the development of maps of time-averaged shear wave velocities over specific depths, and the identification of regional similarities and differences in soil shear stiffness.
This report summarizes the development of a region-wide surficial soil shear wave velocity (Vs ) model based on the unique combination of a large high-spatial-density database of cone penetration test (CPT) logs in the greater Christchurch urban area (> 15, 000 logs as of 1 February 2014) and the Christchurch-specific empirical correlation between soil Vs and CPT data developed by McGann et al. [1, 2]. This model has applications for site characterization efforts via maps of time-averaged Vs over specific depths (e.g. Vs30, Vs10), and for numerical modeling efforts via the identification of typical Vs profiles for different regions and soil behaviour types within Christchurch. In addition, the Vs model can be used to constrain the near-surface velocities for the 3D seismic velocity model of the Canterbury basin [3] currently being developed for the purpose of broadband ground motion simulation. The general development of these region-wide near-surface Vs models includes the following general phases, with each discussed in separate chapters of this report. • An evaluation of the available CPT dataset for suitability, and the definition of other datasets and assumptions necessary to characterize the surficial sediments of the region to 30 m depth. • The development of time-averaged shear wave velocity (Vsz) surfaces for the Christchurch area from the adopted CPT dataset (and supplementary data/assumptions) using spatial interpolation. The Vsz surfaces are used to explore the characteristics of the near-surface soils in the regions and are shown to correspond well with known features of the local geology, the historical ecosystems of the area, and observations made following the 2010- 2011 Canterbury earthquakes. • A detailed analysis of the Vs profiles in eight subregions of Christchurch is performed to assess the variablity in the soil profiles for regions with similar Vsz values and to assess Vsz as a predictive metric for local site response. It is shown that the distrubution of soil shear wave velocity in the Christchurch regions is highly variable both spatially (horizontally) and with depth (vertically) due to the varied geological histories for different parts of the area, and the highly stratified nature of the nearsurface deposits. This variability is not considered to be greatly significant in terms of current simplified site classification systems; based on computed Vs30 values, all considered regions can be categorized as NEHRP sites class D (180 < Vs < 360 m/s) or E (Vs < 180 m/s), however, detailed analysis of the shear wave velocity profiles in different subregions of Christchurch show that the expected surficial site response can vary quite a bit across the region despite the relative similarity in Vs30
This presentation summarizes the development of high-resolution surficial soil velocity models in the Canterbury, New Zealand basin. Shallow (<30m) shear wave velocities were primarily computed based on a combination of a large database of over 15,000 cone penetration test (CPT) logs in and around Christchurch, and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. Large active-source testing at 22 locations and ambient-wavefield surface wave and H/V testing at over 80 locations were utilized in combination with 1700 water well logs to constrain the inter-bedded stratigraphy and velocity of Quaternary sediments up to depths of several hundred meters. Finally, seismic reflection profiles and the ambient-wavefield surface wave data provide constraint on velocities from several hundred meters to several kilometres. At all depths, the high resolution data illustrates the complexity of the soil conditions in the region, and the developed 3D models are presently being used in broadband ground motion simulations to further interpret the observed strong ground motions in the 2010-2011 Canterbury earthquake sequence.
Documenting earthquake-induced ground deformation is significant to assess the characteristics of past and contemporary earthquakes and provide insight into seismic hazard. This study uses airborne light detection and ranging (LiDAR) and conducts multi-disciplinary field techniques to document the surface rupture morphology and evaluate the paleoseismicity and seismic hazard parameters of the Hurunui segment of the Hope Fault in the northern South Island of New Zealand. It also documents and evaluates seismically induced features and ground motion characteristics of the 2010 Darfield and 2011 Christchurch earthquakes in the Port Hills, south of Christchurch. These two studies are linked in that they investigate the near-field coseismic features of large (Mw ~7.1) earthquakes in New Zealand and produce data for evaluating seismic hazards of future earthquakes. In the northern South Island of New Zealand, the Australian-Pacific plate boundary is characterised by strike-slip deformation across the Marlborough Fault System (MFS). The ENE-striking Hope Fault (length: ~230 km) is the youngest and southernmost fault in the MFS, and the second fastest slipping fault in New Zealand. The Hope Fault is a major source of seismic hazard in New Zealand and has ruptured (in-part) historically in the Mw 7.1 1888 Amuri earthquake. In the west, the Hurunui segment of the Hope Fault is covered by beech forest. Hence, its seismic hazard parameters and paleoearthquake chronology were poorly constrained and it was unknown whether the 1888 earthquake ruptured this segment or not and if so, to what extent. Utilising LiDAR and field data, a 29 km-long section of the Hurunui segment of the Hope Fault is mapped. LiDAR-mapping clearly reveals the principal slip zone (PSZ) of the fault and a suite of previously unrecognised structures that form the fault deformation zone (FDZ). FDZ width measurements from 415 locations reveal a spatially-variable, active FDZ up to ~500 m wide with an average width of 200 m. Kinematic analysis of the fault structures shows that the Hurunui segment strikes between 070° and 075° and is optimally oriented for dextral strike-slip within the regional stress field. This implies that the wide FDZ observed is unlikely to result from large-scale fault mis-orientation with respect to regional stresses. The analysis of FDZ width indicates that it increases with increased hanging wall topography and increased topographic relief suggesting that along-strike topographic perturbations to fault geometry and stress states increase fault zone complexity and width. FDZ width also increases where the tips of adjacent PSZ strands locally vary in strike, and where the thickness of alluvial deposits overlying bedrock increases. LiDAR- and photogrammetrically-derived topographic mapping indicates that the boundary between the Hurunui and Hope River segments is characterised by a ~850-m-wide right stepover and a 9º-14° fault bend. Paleoseismic trenching at Hope Shelter site reveals that 6 earthquakes occurred at A.D. 1888, 1740-1840, 1479-1623, 819-1092, 439-551, and 373- 419. These rupture events have a mean recurrence interval of ~298 ± 88 yr and inter-event times ranging from 98 to 595 yrs. The variation in the inter-event times is explained by (1) coalescing rupture overlap from the adjacent Hope River segment on to the Hurunui segment at the study site, (2) temporal clustering of large earthquakes on the Hurunui segment, and/or (3) ‘missing’ rupture events. It appears that the first two options are more plausible to explain the earthquake chronologies and rupture behaviour on the Hurunui segment, given the detailed nature of the geologic and chronologic investigations. This study provides first evidence for coseismic multi-segment ruptures on the Hope Fault by identifying a rupture length of 44-70 km for the 1888 earthquake, which was not confined to the Hope River segment (primary source for the 1888 earthquake). LiDAR data is also used to identify and measure dextral displacements and scarp heights from the PSZ and structures within the FDZ along the Hurunui segment. Reconstruction of large dextrally-offset geomorphic features shows that the vertical component of slip accounts for only ~1% of the horizontal displacements and confirms that the fault is predominantly strike-slip. A strong correlation exists between the dextral displacements and elevations of geomorphic features suggesting the possibility of age correlation between the geomorphic features. A mean single event displacement (SED) of 3.6 ± 0.7 m is determined from interpretation of sets of dextral displacements of ≤ 25 m. Using the available surface age data and the cumulative dextral displacements from Matagouri Flat, McKenzie Fan, Macs Knob and Hope River sites, and the mean SED, a mean slip rate of 12.2 ± 2.4 mm/yr, and a mean recurrence interval of ~320 ± 120 yr, and a potential earthquake magnitude of Mw 7.2 are determined for the Hurunui segment. This study suggests that the fault slip rate has been constant over the last ~15000 yr. Strong ground motions from the 2010 Darfield (Canterbury) earthquake displaced boulders and caused ground damage on some ridge crests in the Port Hills. However, the 2011 Christchurch earthquake neither displaced boulders nor caused ground damage at the same ridge crests. Documentation of locations (~400 m a.s.l.), lateral displacements (8-970 cm), displacement direction (250° ± 20°) of displaced boulders, in addition to their hosting socket geometries (< 1 cm to 50 cm depth), the orientation of the ridges (000°-015°) indicate that boulders have been displaced in the direction of instrumentally recorded transient peak ground horizontal displacements nearby and that the seismic waves have been amplified at the study sites. The co-existence of displaced and non-displaced boulders at proximal sites suggests small-scale ground motion variability and/or varying boulder-ground dynamic interactions relating to shallow phenomena such as variability in soil depth, bedrock fracture density and/or microtopography on the bedrock-soil interface. Shorter shaking duration of the 2011 Christchurch event, differing frequency contents and different source characteristics were all factors that may have contributed to generating circumstances less favourable to boulder displacement in this earthquake. Investigating seismically induced features, fault behaviour, site effects on the rupture behaviour, and site response to the seismic waves provides insights into fault rupture hazards.
Spatial variations in river facies exerted a strong influence on the distribution of liquefaction features observed in Christchurch during the 2010-11 Canterbury Earthquake Sequence (CES). Liquefaction and liquefaction-induced ground deformation was primarily concentrated near modern waterways and areas underlain by Holocene fluvial deposits with shallow water tables (< 1 to 2 m). In southern Christchurch, spatial variations of liquefaction and subsidence were documented in the suburbs within inner meander loops of the Heathcote River. Newly acquired geospatial data, geotechnical reports and eye-witness discussions are compiled to provide a detailed account of the surficial effects of CES liquefaction and ground deformation adjacent to the Heathcote River. LiDAR data and aerial photography are used to produce a new series of original figures which reveal the locations of recurrent liquefaction and subsidence. To investigate why variable liquefaction patterns occurred, the distribution of surface ejecta and associated ground damage is compared with near-surface sedimentologic, topographic, and geomorphic variability to seek relationships between the near-surface properties and observed ground damages. The most severe liquefaction was concentrated within a topographic low in the suburb of St Martins, an inner meander loop of the Heathcote River, with liquefaction only minor or absent in the surrounding areas. Subsurface investigations at two sites in St Martins enable documentation of fluvial stratigraphy, the expressions of liquefaction, and identification of pre-CES liquefaction features. Excavation to water table depths (~1.5 m below the surface) across sand boils reveals multiple generations of CES liquefaction dikes and sills that cross-cut Holocene fluvial and anthropogenic stratigraphy. Based on in situ geotechnical tests (CPT) indicating sediment with a factor of safety < 1, the majority of surface ejecta was sourced from well-sorted fine to medium sand at < 5 m depth, with the most damaging liquefaction corresponding with the location of a low-lying sandy paleochannel, a remnant river channel from the Holocene migration of the meander in St Martins. In the adjacent suburb of Beckenham, where migration of the Heathcote River has been laterally confined by topography associated with the volcanic lithologies of Banks Peninsula, severe liquefaction was absent with only minor sand boils occurring closest to the modern river channel. Auger sampling across the suburb revealed thick (>1 m) clay-rich overbank and back swamp sediments that produced a stratigraphy which likely confined the units susceptible to liquefaction and prevented widespread ejection of liquefied material. This analysis suggests river migration promotes the formation and preservation of fluvial deposits prone to liquefaction. Trenching revealed the strongest CES earthquakes with large vertical accelerations favoured sill formation and severe subsidence at highly susceptible locations corresponding with an abandoned channel. Less vulnerable sites containing deeper and thinner sand bodies only liquefied in the strongest and most proximal earthquakes forming minor localised liquefaction features. Liquefaction was less prominent and severe subsidence was absent where lateral confinement of a Heathcote meander has promoted the formation of fluvial stratum resistant to liquefaction. Correlating CES liquefaction with geomorphic interpretations of Christchurch’s Heathcote River highlights methods in which the performance of liquefaction susceptibility models can be improved. These include developing a reliable proxy for estimating soil conditions in meandering fluvial systems by interpreting the geology and geomorphology, derived from LiDAR data and modern river morphology, to improve the methods of accounting for the susceptibility of an area. Combining geomorphic interpretations with geotechnical data can be applied elsewhere to identify regional liquefaction susceptibilities, improve existing liquefaction susceptibility datasets, and predict future earthquake damage.
Oblique-convergent plate collision between the Pacific and Australian plates across the South Island has resulted in shallow, upper crustal earthquake activity and ground surface deformation. In particular the Porters Pass - Amberley Fault Zone displays a complex hybrid zone of anastomosing dextral strike-slip and thrust/reverse faulting which includes the thrust/reverse Lees Valley Fault Zone and associated basin deformation. There is a knowledge gap with respect to the paleoseismicity of many of the faults in this region including the Lees Valley Fault Zone. This study aimed to investigate the earthquake history of the fault at a selected location and the structural and geomorphic development of the Lees Valley Fault Zone and eastern rangefront. This was investigated through extensive structural and geomorphic mapping, GPS field surveying, vertical aerial photo interpretation, analysis of Digital Elevation Models, paleoseismic trenching and optically stimulated luminescence dating. This thesis used a published model for tectonic geomorphology development of mountain rangefronts to understand the development of Lees Valley. Rangefront geomorphology is investigated through analysis of features such as rangefront sinuosity and faceted spurs and indicates the recently active and episodic nature of the uplifted rangefront. Analysis of fault discontinuity, fault splays, distribution of displacement, fault deformation zone and limited exposure of bedrock provided insight into the complex structure of the fault zone. These observations revealed preserved, earlier rangefronts, abandoned and uplifted within the eastern ranges, indicating a basinward shift in focus of faulting and an imbricate thrust wedge development propagating into the footwall of the fault zone and along the eastern ranges of Lees Valley. Fault scarp deformation analysis indicated multiple events have produced the deformation present preserved by the active fault trace in the northern valley. Vertical deformation along this scarp varied with a maximum of 11.5 m and an average of 5 m. Field mapping revealed fan surfaces of various ages have been offset and deformed, likely during the Holocene, based on expected relative surface ages. Geomorphic and structural mapping highlighted the effect of cross-cutting and inherited structures on the Lees Valley Fault, resulting in a step-over development in the centre of the eastern range-bounding trace. Paleoseismic trenching provided evidence of at least two earthquakes, which were constrained to post 21.6 ± 2.3 ka by optically stimulated luminescence dating. Single event displacements (1.48 ± 0.08 m), surface rupture earthquake magnitudes (Mw 6.7 ± 0.1, with potential to produce ≥ 7.0), and a minimum recurrence interval (3.6 ± 0.3 ka) indicated the Lees Valley Fault is an active structure capable of producing significant earthquake events. Results from this study indicate that the Lees Valley Fault Zone accommodates an important component of the Porters Pass - Amberley Fault Zone deformation and confirms the fault as a source of potentially damaging, peak ground accelerations in the Canterbury region. Remnants of previous rangefronts indicate a thrust wedge development of the Lees Valley Fault Zone and associated ranges that can potentially be used as a model of development for other thrust-fault bounded basins.
he strong motion station at Heathcote Valley School (HVSC) recorded unusually high peak ground accelerations (2.21g vertical and 1.41g horizontal) during the February 2011 Christchurch earthquake. Ground motions recorded at HVSC in numerous other events also exhibited consistently higher intensities compared with nearby strong motion stations. We investigated the underlying causes of such high intensity ground motions at HVSC by means of 2D dynamic finite element analyses, using recorded ground motions during the 2010-2011 Canterbury earthquake sequence. The model takes advantage of a LiDAR-based digital elevation model (DEM) to account for the surface topography, while the geometry and dynamic properties of the surficial soils are characterized by seismic cone penetration tests (sCPT) and Multi-Channel Analyses of Surface Waves (MASW). Comparisons of simulated and recorded ground motions suggests that our model performs well for distant events, while for near-field events, ground motions recorded at the adopted reference station at Lyttelton Port are not reasonable input motions for the simulation. The simulations suggest that Rayleigh waves generated at the inclined interface of the surficial colluvium and underlying volcanic rock strongly affect the ground motions recorded at HVSC, in particular, being the dominant contributor to the recorded vertical motions.
In the aftermath of the 2010-2011 Canterbury Earthquake Sequence (CES), the location of Christchurch-City on the coast of the Canterbury Region (New Zealand) has proven crucial in determining the types of- and chains of hazards that impact the city. Very rapidly, the land subsidence of up to 1 m (vertical), and the modifications of city’s waterways – bank sliding, longitudinal profile change, sedimentation and erosion, engineered stop-banks… - turned rainfall and high-tides into unprecedented floods, which spread across the eastern side of the city. Within this context, this contribution presents two modeling results of potential floods: (1) results of flood models and (2) the effects of further subsidence-linked flooding – indeed if another similar earthquake was to strike the city, what could be the scenarios of further subsidence and then flooding. The present research uses the pre- and post-CES LiDAR datasets, which have been used as the boundary layer for the modeling. On top of simple bathtub model of inundation, the river flood model was conducted using the 2-D hydrodynamic code NAYS-2D developed at the University of Hokkaido (Japan), using a depth-averaged resolution of the hydrodynamic equations. The results have shown that the area the most at risk of flooding are the recent Holocene sedimentary deposits, and especially the swamplands near the sea and in the proximity of waterways. As the CES drove horizontal and vertical displacement of the land-surface, the surface hydrology of the city has been deeply modified, increasing flood risks. However, it seems that scientists and managers haven’t fully learned from the CES, and no research has been looking at the potential future subsidence in further worsening subsidence-related floods. Consequently, the term “coastal quake”, coined by D. Hart is highly topical, and most especially because most of our modern cities and mega-cities are built on estuarine Holocene sediments.
Recent field investigations were carried out to define the shear wave velocity (VS) profile and site periods across the Canterbury region, supplementing earlier efforts in urban Christchurch. Active source surface wave testing, ambient wave field (passive) and H/V spectral ratio methods were used to characterise the soil profile in the region. H/V spectral ratio peaks indicate site periods in the range of 5-7 seconds across much of the Canterbury Plains, broadly consistent with those based on a 1D velocity model for the region. Site periods decrease rapidly in the vicinity of the Canterbury foothills and the Banks Peninsula outcrops. In Christchurch, the Riccarton Gravels result in a significant mode of vibration that has a much shorter period than the site period of the entire soil column down to basement rock.
A major lesson from the 2011 Christchurch earthquake was the apparent lack of ductility of some lightly reinforced concrete (RC) wall structures. In particular, the structural behaviour of the critical wall in the Gallery Apartments building demonstrated that the inelastic deformation capacity of a structure, as well as potentially brittle failure of the reinforcement, is dependent on the level of bond deterioration between reinforcement and surrounding concrete that occurs under seismic loading. This paper presents the findings of an experimental study on bond behaviour between deformed reinforcing bars and the surrounding concrete. Bond strength and relative bond slip was evaluated using 75 pull-out tests under monotonic and cyclic loading. Variations of the experiments include the loading rate, loading history, concrete strength (25 to 70 MPa), concrete age, cover thickness, bar diameter (16 and 20 mm), embedded length, and the position of the embedded bond region within the specimen (deep within or close to free surface). Select test results are presented with inferred implications for RC structures.
The 2010-2011 Canterbury earthquakes were recorded over a dense strong motion network in the near-source region, yielding significant observational evidence of seismic complexities, and a basis for interpretation of multi-disciplinary datasets and induced damage to the natural and built environment. This paper provides an overview of observed strong motions from these events and retrospective comparisons with both empirical and physics-based ground motion models. Both empirical and physics-based methods provide good predictions of observations at short vibration periods in an average sense. However, observed ground motion amplitudes at specific locations, such as Heathcote Valley, are seen to systematically depart from ‘average’ empirical predictions as a result of near surface stratigraphic and topographic features which are well modelled via sitespecific response analyses. Significant insight into the long period bias in empirical predictions is obtained from the use of hybrid broadband ground motion simulation. The comparison of both empirical and physics-based simulations against a set of 10 events in the sequence clearly illustrates the potential for simulations to improve ground motion and site response prediction, both at present, and further in the future.
This paper presents an overview of the soil profile characteristics at strong motion station (SMS) locations in the Christchurch Central Business District (CBD) based on recently completed geotechnical site investigations. Given the variability of Christchurch soils, detailed investigations were needed in close vicinity to each SMS. In this regard, CPT, SPT and borehole data, and shear wave velocity (Vs) profiles from surface wave dispersion data in close vicinity to the SMSs have been used to develop detailed representative soil profiles at each site and to determine site classes according to the New Zealand standard NZS1170.5. A disparity between the NZS1170.5 site classes based on Vs and SPT N60 investigation techniques is highlighted, and additional studies are needed to harmonize site classification based on these techniques. The short period mode of vibration of soft deposits above gravels, which are found throughout Christchurch, are compared to the long period mode of vibration of the entire soil profile to bedrock. These two distinct modes of vibration require further investigation to determine their impact on the site response. According to current American and European approaches to seismic site classification, all SMSs were classified as problematic soil sites due to the presence of liquefiable strata, soils which are not directly accounted for by the NZS1170.5 approach.
Recurrent liquefaction in Christchurch during the 2010-2011 Canterbury earthquake sequence created a wealth of shallow subsurface intrusions with geometries and orientations governed by (1) strong ground motion severity and duration, and (2) intrinsic site characteristics including liquefaction susceptibility, lateral spreading severity, geomorphic setting, host sediment heterogeneity, and anthropogenic soil modifications. We present a suite of case studies that demonstrate how each of these characteristics influenced the geologic expressions of contemporary liquefaction in the shallow subsurface. We compare contemporary features with paleo-features to show how geologic investigations of recurrent liquefaction can provide novel insights into the shaking characteristics of modern and paleo-earthquakes, the influence of geomorphology on liquefaction vulnerability, and the possible controls of anthropogenic activity on the geologic record. We conclude that (a) sites of paleo-liquefaction in the last 1000-2000 years corresponded with most severe liquefaction during the Canterbury earthquake sequence, (b) less vulnerable sites that only liquefied in the strongest and most proximal contemporary earthquakes are unlikely to have liquefied in the last 1000-2000 years or more, (c) proximal strong earthquakes with large vertical accelerations favoured sill formation at some locations, (d) contemporary liquefaction was more severe than paleoliquefaction at all study sites, and (e) stratigraphic records of successive dike formation were more complete at sites with severe lateral spreading, (f) anthropogenic fill suppressed surface liquefaction features and altered subsurface liquefaction architecture.
Following the Mw 6.2 Christchurch Earthquake on 22 February 2011, extensive ground cracking in loessial soils was reported in some areas of the Port Hills, southeast of central Christchurch. This study was undertaken to investigate the mechanisms of earthquake-induced ground damage on the eastern side of the Hillsborough Valley. A zone of extensional cracking up to 40m wide and 600m long was identified along the eastern foot-slope, accompanied by compression features and spring formation at the toe of the slope. An engineering geological and geomorphological model was developed for the eastern Hillsborough Valley that incorporates geotechnical investigation data sourced from the Canterbury Geotechnical Database (CGD), the findings of trenching and seismic refraction surveying carried out for this research, and interpretation of historical aerial photographs. The thickness and extent of a buried peat swamp at the base of the slope was mapped, and found to coincide with significant compression features. Ground cracking was found to have occurred entirely within loess-colluvium and to follow the apices of pre-1920s tunnel-gully fan debris at the southern end of the valley. The ground-cracking on the eastern side of the Hillsborough Valley is interpreted to have formed through tensile failure of the loess-colluvium. Testing was carried out to determine the tensile strength of Port Hills loess colluvium as a function of water content and density, in order to better understand the occurrence and distribution of the observed ground cracking. A comprehensive review of the soil tensile strength testing literature was undertaken, from which a test methodology was developed. Results show remoulded loess-colluvium to possess tensile strength of 7 - 28 kPa across the range of tested moisture contents (10-15%) and dry densities (1650-1900kg/m3). A positive linear relationship was observed between tensile strength and dry density, and a negative linear relationship between moisture content and tensile strength. The observed ground damage and available geotechnical information (inclinometer and piezometer records provided by the Earthquake Commission) were together used to interpret the mechanism(s) of slope movement that occurred in the eastern Hillsborough Valley. The observed ground damage is characteristic of translational movement, but without the development of lateral release scarps, or a basal sliding surface - which was not located during drilling. It is hypothesised that shear displacement has been accommodated by multiple slip surfaces of limited extent within the upper 10m of the slope. Movement has likely occurred within near-saturated colluvial units that have lost strength during earthquake shaking. The eastern Hillsborough Valley is considered to be an ‘incipient translational slide’, as both the patterns of damage and shearing are consistent with the early stages of such slide development. Sliding block analysis was utilised to understand how the eastern Hillsborough Valley may perform in a future large magnitude earthquake. Known cumulative displacements of ~0.3m for eastern Hillsborough Valley during the 2010-2011 Canterbury Earthquake Sequence were compared with modelled slope displacements to back-analyse a lower-bound yield acceleration of 0.2 - 0.25g. Synthetic broadband modelling for future Alpine and Hope Fault earthquakes indicates PGAs of approximately 0.08g for soil sites in the Christchurch area, as such, slope movement is unlikely to be reactivated by an Alpine Fault or Hope Fault earthquake. This does not take into account the possible role of strength loss due to excess pore pressure that may occur during these future events.
This thesis is concerned with springs that appeared in the Hillsborough, Christchurch during the 2010-2011 Canterbury Earthquake Sequence, and which have continued to discharge groundwater to the surface to the present time. Investigations have evolved, measurements of discharge at selected sites, limited chemical data on anions and isotope analysis. The springs are associated with earthquake generated fissures (extensional) and compression zones, mostly in loess-colluvium soils of the valley floor and lower slopes. Extensive peat swamps are present in the Hillsborough valley, with a groundwater table at ~1m below ground. The first appearance of the ‘new’ springs took place following the Mw 7.1 Darfield Earthquake on 4 September 2010, and discharges increased both in volume and extent of the Christchurch Mw 6.3 Earthquake of 22 February 2011. Five monitored sites show flow rates in the range of 4.2-14.4L/min, which have remained effectively constant for the duration of the study (2014-2015). Water chemistry analysis shows that the groundwater discharges are sourced primarily from volcanic bedrocks which underlies the valley at depths ≤50m below ground level. Isotope values confirm similarities with bedrock-sourced groundwater, and the short term (hours-days) influence of extreme rainfall events. Cyclone Lusi (2013-2014) affects were monitored and showed recovery of the bedrock derived water signature within 72 hours. Close to the mouth of the valley sediments interfinger with Waimakiriri River derived alluvium bearing a distinct and different isotope signature. Some mixing is evident at certain locations, but it is not clear if there is any influence from the Huntsbury reservoir which failed in the Port Hills Earthquake (22 February 2011) and stored groundwater from the Christchurch artesian aquifer system (Riccarton Gravel).
The Canterbury Earthquake Sequence (CES) of 2010-2011 produced large seismic moments up to Mw 7.1. These large, near-to-surface (<15 km) ruptures triggered >6,000 rockfall boulders on the Port Hills of Christchurch, many of which impacted houses and affected the livelihoods of people within the impacted area. From these disastrous and unpredicted natural events a need arose to be able to assess the areas affected by rockfall events in the future, where it is known that a rockfall is possible from a specific source outcrop but the potential boulder runout and dynamics are not understood. The distribution of rockfall deposits is largely constrained by the physical properties and processes of the boulder and its motion such as block density, shape and size, block velocity, bounce height, impact and rebound angle, as well as the properties of the substrate. Numerical rockfall models go some way to accounting for all the complex factors in an algorithm, commonly parameterised in a user interface where site-specific effects can be calibrated. Calibration of these algorithms requires thorough field checks and often experimental practises. The purpose of this project, which began immediately following the most destructive rupture of the CES (February 22, 2011), is to collate data to characterise boulder falls, and to use this information, supplemented by a set of anthropogenic boulder fall data, to perform an in-depth calibration of the three-dimensional numerical rockfall model RAMMS::Rockfall. The thesis covers the following topics: • Use of field data to calibrate RAMMS. Boulder impact trails in the loess-colluvium soils at Rapaki Bay have been used to estimate ranges of boulder velocities and bounce heights. RAMMS results replicate field data closely; it is concluded that the model is appropriate for analysing the earthquake-triggered boulder trails at Rapaki Bay, and that it can be usefully applied to rockfall trajectory and hazard assessment at this and similar sites elsewhere. • Detailed analysis of dynamic rockfall processes, interpreted from recorded boulder rolling experiments, and compared to RAMMS simulated results at the same site. Recorded rotational and translational velocities of a particular boulder show that the boulder behaves logically and dynamically on impact with different substrate types. Simulations show that seasonal changes in soil moisture alter rockfall dynamics and runout predictions within RAMMS, and adjustments are made to the calibration to reflect this; suggesting that in hazard analysis a rockfall model should be calibrated to dry rather than wet soil conditions to anticipate the most serious outcome. • Verifying the model calibration for a separate site on the Port Hills. The results of the RAMMS simulations show the effectiveness of calibration against a real data set, as well as the effectiveness of vegetation as a rockfall barrier/retardant. The results of simulations are compared using hazard maps, where the maximum runouts match well the mapped CES fallen boulder maximum runouts. The results of the simulations in terms of frequency distribution of deposit locations on the slope are also compared with those of the CES data, using the shadow angle tool to apportion slope zones. These results also replicate real field data well. Results show that a maximum runout envelope can be mapped, as well as frequency distribution of deposited boulders for hazard (and thus risk) analysis purposes. The accuracy of the rockfall runout envelope and frequency distribution can be improved by comprehensive vegetation and substrate mapping. The topics above define the scope of the project, limiting the focus to rockfall processes on the Port Hills, and implications for model calibration for the wider scientific community. The results provide a useful rockfall analysis methodology with a defensible and replicable calibration process, that has the potential to be applied to other lithologies and substrates. Its applications include a method of analysis for the selection and positioning of rockfall countermeasure design; site safety assessment for scaling and demolition works; and risk analysis and land planning for future construction in Christchurch.
This thesis addresses the topic of local bond behaviour in RC structures. The mechanism of bond refers to the composite action between deformed steel reinforcing bars and the surrounding concrete. Bond behaviour is an open research topic with a wide scope, particularly because bond it is such a fundamental concept to structural engineers. However, despite many bond-related research findings having wide applications, the primary contribution of this research is an experimental evaluation of the prominent features of local bond behaviour and the associated implications for the seismic performance of RC structures. The findings presented in this thesis attempt to address some structural engineering recommendations made by the Canterbury Earthquakes Royal Commission following the 2010-2011 Canterbury (New Zealand) earthquake sequence. A chapter of this thesis discusses the structural behaviour of flexure-dominated RC wall structures with an insufficient quantity of longitudinal reinforcement, among other in situ conditions, that causes material damage to predominantly occur at a single crack plane. In this particular case, the extent of concrete damage and bond deterioration adjacent to the crack plane will influence the ductility capacity that is effectively provided by the reinforcing steel. As a consequence of these in situ conditions, some lightly reinforced wall buildings in Christchurch lost their structural integrity due to brittle fracture of the longitudinal reinforcement. With these concerning post-earthquake observations in mind, there is the underlying intention that this thesis presents experimental evidence of bond behaviour that allows structural engineers to re-assess their confidence levels for the ability of lightly reinforced concrete structures to achieve the life-safety seismic performance objective the ultimate limit state. Three chapters of this thesis are devoted to the experimental work that was conducted as the main contribution of this research. Critical details of the experimental design, bond testing method and test programme are reported. The bond stress-slip relationship was studied through 75 bond pull-out tests. In order to measure the maximum local bond strength, all bond tests were carried out on deformed reinforcing bars that did not yield as the embedded bond length was relatively short. Bond test results have been presented in two separate chapters in which 48 monotonic bond tests and 27 cyclic bond tests are presented. Permutations of the experiments include the loading rate, cyclic loading history, concrete strength (25 to 70 MPa), concrete age, cover thickness, bar diameter (16 and 20 mm), embedded length, and position of the embedded bond region within the specimen (close or far away to the free surface). The parametric study showed that the concrete strength significantly influences the maximum bond strength and that it is reasonable to normalise the bond stress by the square-root of the concrete compressive strength, √(f'c). The generalised monotonic bond behaviour is described within. An important outcome of the research is that the measured bond strength and stiffness was higher than stated by the bond stress-slip relationship in the fib Model Code 2010. To account for these observed differences, an alternative model is proposed for the local monotonic bond stress-slip relationship. Cyclic bond tests showed a significant proportion of the total bond degradation occurs after the loading cycle in the peak bond strength range, which is when bond slip has exceeded 0.5 mm. Subsequent loading to constant slip values showed a linear relationship between the amount of bond strength degradation and the log of the number of cycles that were applied. To a greater extent, the cyclic bond deterioration depends on the bond slip range, regardless of whether the applied load cycling is half- or fully-reversed. The observed bond deterioration and hysteretic energy dissipated during cyclic loading was found to agree reasonably well between these cyclic tests with different loading protocols. The cyclic bond deterioration was also found to be reasonably consistent exponential damage models found in the literature. This research concluded that the deformed reinforcing bars used in NZ construction, embedded in moderate to high strength concrete, are able to develop high local bond stresses that are mobilised by a small amount of local bond slip. Although the relative rib geometry was not varied within this experimental programme, a general conclusion of this thesis is that deformed bars currently available in NZ have a relative rib bearing area that is comparatively higher than the test bars used in previous international research. From the parametric study it was found that the maximum monotonic bond strength is significant enhanced by dynamic loading rates. Experimental evidence of high bond strength and initial bond stiffness generally suggests that only a small amount of local bond slip that can occur when the deformed test bar was subjected to large tension forces. Minimal bond slip and bond damage limits the effective yielding length that is available for the reinforcing steel to distribute inelastic material strains. Consequently, the potential for brittle fracture of the reinforcement may be a more problematic and widespread issue than is apparent to structural engineers. This research has provided information that improve the reliability of engineering predictions (with respect to ductility capacity) of maximum crack widths and the extent of bond deterioration that might occur in RC structures during seismic actions.