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.
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.
1. Background and Objectives This poster presents results from ground motion simulations of small-to-moderate magnitude (3.5≤Mw≤5.0) earthquake events in the Canterbury, New Zealand region using the Graves and Pitarka (2010,2015) methodology. Subsequent investigation of systematic ground motion effects highlights the prediction bias in the simulations which are also benchmarked against empirical ground motion models (e.g. Bradley (2013)). In this study, 144 earthquake ruptures, modelled as point sources, are considered with 1924 quality-assured ground motions recorded across 45 strong motion stations throughout the Canterbury region, as shown in Figure 1. The majority of sources are Mw≥4.0 and have centroid depth (CD) 10km or shallower. Earthquake source descriptions were obtained from the GeoNet New Zealand earthquake catalogue. The ground motion simulations were performed within a computational domain of 140km x 120km x 46km with a finite difference grid spacing of 0.1km. The low-frequency (LF) simulations utilize the 3D Canterbury Velocity Model while the high-frequency (HF) simulations utilize a generic regional 1D velocity model. In the LF simulations, a minimum shear wave velocity of 500m/s is enforced, yielding a maximum frequency of 1.0Hz.
Severe liquefaction was repeatedly observed during the 2010 - 2011 C hristchurch earthquake s , particularly affecting deposits of fine sands and silty sands of recent fluvial or estuarine origin. The effects of liquefaction included major sliding of soil tow ard water bodies ( i.e. lateral spreading ) rang ing from centimetres to several metres. In this paper, a series of undrained cyclic torsional shear tests were conducted to evaluate the liquefaction and extremely large deformation properties of Christchurch b oiled sand . In these tests, the simple shear conditions were reproduced in order to apply realistic stress conditions that soil s experience in the field during horizontal seismic shaking. Several hollow cylindrical medium dense specimens ( D r = 50%) were pr epared by pluviation method, isotropically consolidated at an effective stress of 100 kPa and then cyclically sheared under undrained conditions up to 10 0% double amplitude shear strain (γ DA ) . The cyclic strength at different levels of γ DA of 7.5%, 15%, 3 0 % and 6 0%, development of extremely large post - liquefaction deformation and shear strain locali s ation properties were assessed from the analysis of the effective stress paths and stress - strain responses . To reveal possible distinctiveness, the cyclic undra ined behaviour of CHCH boiled sand was compared with that of Toyoura sand previously examined under similar testing conditions
his poster presents the ongoing development of a 3D Canterbury seismic velocity model which will be used in physics-based hybrid broadband ground motion simulation of the 2010-2011 Canterbury earthquakes. Velocity models must sufficiently represent critical aspects of the crustal structure over multiple length scales which will influence the results of the simulations. As a result, numerous sources of data are utilized in order to provide adequate resolution where necessary. Figure 2: (a) Seismic reflection line showing P-wave velocities and significant geologic horizons (Barnes et al. 2011), and (b) Shear wave profiles at 10 locations (Stokoe et al. 2013). Figure 4: Cross sections of the current version of the Canterbury velocity model to depths of 10km as shown in Figure 1: (a) at a constant latitude value of -43.6˚, and (b) at a constant longitude value of 172.64˚. 3. Ground Surface and Geologic Horizon Models Figure 3: (a) Ground surface model derived from numerous available digital elevation models, and (b) Base of the Quaternary sediments derived from structural contours and seismic reflection line elevations. The Canterbury region has a unique and complex geology which likely has a significant impact on strong ground motions, in particular the deep and loose deposits of the Canterbury basin. The Canterbury basin has several implications on seismic wave phenomena such as long period ground motion amplification and wave guide effects. Using a realistic 3D seismic velocity model in physics-based ground motion simulation will implicitly account for such effects and the resultant simulated ground motions can be studied to gain a fundamental understanding of the salient ground motion phenomena which occurred during the Canterbury earthquakes, and the potential for repeat occurrences in the Canterbury region. Figure 1 shows the current model domain as a rectangular area between Lat=[-43.2˚,-44.0˚], and Lon=[171.5˚,173.0˚]. This essentially spans the area between the foot of the Southern Alps in the North West to Banks Peninsula in the East. Currently the model extends to a depth of 50km below sea level.
Between 2010 and 2011, Canterbury experienced a series of four large earthquake events with associated aftershocks which caused widespread damage to residential and commercial infrastructure. Fine grained and uncompacted alluvial soils, typical to the Canterbury outwash plains, were exposed to high peak ground acceleration (PGA) during these events. This rapid increase in PGA induced cyclic strain softening and liquefaction in the saturated, near surface alluvial soils. Extensive research into understanding the response of soils in Canterbury to dynamic loading has since occurred. The Earthquake Commission (EQC), the Ministry of Business and Employment (MBIE), and the Christchurch City Council (CCC) have quantified the potential hazards associated with future seismic events. Theses bodies have tested numerous ground improvement design methods, and subsequently are at the forefront of the Canterbury recovery and rebuild process. Deep Soil Mixing (DSM) has been proven as a viable ground improvement foundation method used to enhance in situ soils by increasing stiffness and positively altering in situ soil characteristics. However, current industry practice for confirming the effectiveness of the DSM method involves specific laboratory and absolute soil test methods associated with the mixed column element itself. Currently, the response of the soil around the columns to DSM installation is poorly understood. This research aims to understand and quantify the effects of DSM columns on near surface alluvial soils between the DSM columns though the implementation of standardised empirical soil test methods. These soil strength properties and ground improvement changes have been investigated using shear wave velocity (Vs), soil behaviour and density response methods. The results of the three different empirical tests indicated a consistent improvement within the ground around the DSM columns in sandier soils. By contrast, cohesive silty soils portrayed less of a consistent response to DSM, although still recorded increases. Generally, within the tests completed 50 mm from the column edge, the soil response indicated a deterioration to DSM. This is likely to be a result of the destruction of the soil fabric as the stress and strain of DSM is applied to the un‐mixed in situ soils. The results suggest that during the installation of DSM columns, a positive ground effect occurs in a similar way to other methods of ground improvement. However, further research, including additional testing following this empirical method, laboratory testing and finite 2D and 3D modelling, would be useful to quantify, in detail, how in situ soils respond and how practitioners should consider these test results in their designs. This thesis begins to evaluate how alluvial soils tend to respond to DSM. Conducting more testing on the research site, on other sites in Christchurch, and around the world, would provide a more complete data set to confirm the results of this research and enable further evaluation. Completing this additional research could help geotechnical DSM practitioners to use standardised empirical test methods to measure and confirm ground improvement rather than using existing test methods in future DSM projects. Further, demonstrating the effectiveness of empirical test methods in a DSM context is likely to enable more cost effective and efficient testing of DSM columns in future geotechnical projects.
Extended Direct Analysis (EDA), developed at the University of Canterbury, is an advance on the AISC Direct Analysis method for the analysis of frames subjected to static forces. EDA provides a faster, simple and more rational way to properly consider the second-order effects, initial residual stresses (IRS) and the initial imperfections or steel structures under one directional loading than conventional analysis methods. This research applied the EDA method to quantify the effect of member overstrength on frame behaviour for a single storey frame. Also, the effects of IRS, which were included in the EDA static analysis, but which are not considered explicitly in non-linear seismic analysis, were evaluated in two ways. Firstly, they were considered for simple structures subject to increasing cyclic displacement in different directions. Secondly, incremental dynamic analysis with realistic ground motion was used to quantify the likely effect of IRS in earthquakes. It was found that, contrary to traditional wisdom and practice, greater member strengths can result in lower frame strengths for frames under monotonic lateral loading. The structural lateral capacity of the overstrength case was reduced by 6% compared to the case using the dependable member strengths. Also, it resulted significantly different in member demands. Therefore, it is recommended that when either plastic analysis or EDA is used, that both upper and lower bounds on the likely member strength should be considered to determine the total frame strength and the member demands. Results of push-pull analysis under displacement control showed that for IRS ratio, gamma < 0.5 and axial compressive force ratio, N*/Ns, up to 0.5, IRS did affect the structural behaviour in the first half cycle. However, the behavior in the later cycles was not significantly affected. It also showed that the effect of initial residual stresses in the frame was less significant than for the column alone when the column was subjected to similar axial compressive force. The incremental dynamic analysis results from both cantilever column and the three-storey steel frame showed that by increasing gamma = 0 to 0.5, the effect of IRS on seismic responses, based on the 50% confidence level, was less than 3% for N*/Ns, up to 0.5.
Current research in geotechnical engineering at the University of Canterbury includes a number of laboratory testing programmes focussed on understanding the behaviour of natural soil deposits in Christchurch during the 2010-2011 Canterbury Earthquake Sequence. Many soils found in Christchurch are sands or silty sands with little to no plasticity, making them very difficult to sample using established methods. The gel-push sampling methodology, developed by Kiso-Jiban Consultants in Japan, was developed to address some of the deficiencies of existing sampling techniques and has been deployed on two projects in Christchurch. Gel push sampling is carried out with a range of samplers which are modified versions of existing technology, and the University of Canterbury has acquired three versions of the tools (GP-S, GP-Tr, GP-D). Soil samples are extracted from the bottom of a freshly drilled borehole and are captured within a liner barrel, close to 1m in length. A lubricating polymer gel coats the outside of the soil sample as it enters the liner barrel. The frictional rubbing which normally occurs on the sides of the soil samples using existing techniques is eliminated by the presence of the polymer gel. The operation of the gel-push samplers is significantly more complicated than conventional push-tube samplers, and in the initial trials a number of operational difficulties were encountered, requiring changes to the sampling procedures. Despite these issues, a number of high quality soil samples were obtained on both projects using the GP-S sampler to capture silty soil. Attempts were made to obtain clean sands using a different gel-push sampler (GP-TR) in the Red Zone. The laboratory testing of these sands indicated that they were being significantly disturbed during the sampling and/or transportation procedures. While it remains too early to draw definitive conclusions regarding the performance of the gel-push samplers, the methodology has provided some promising results. Further trialling of the tools are required to refine operating procedures understand the full range of soil conditions which can be successfully sampled using the tools. In parallel with the gel-push trials, a Dames and Moore fixed-piston sampler has been used by our research partners from Berkeley to obtain soil samples at a number of sites within Christchurch. This sampler features relatively short (50cm), thin-walled liner barrels which is advanced into the ground under the action of hydraulic pressure. By reducing the overall length of the soil being captured, the disturbance to the soil as it enters the liner barrel is significantly reduced. The Dames and Moore sampler is significantly easier to operate than the gel-push sampler, and past experience has shown it to be successful in soft, plastic materials (i.e. clays and silty clays). The cyclic resistance of one silty clay obtained using both the gel-push and Dames & Moore samplers has been found to be very similar, and ongoing research aims to establish whether similar results are obtained for different soil types, including silty materials and clean sands.
The Acheron rock avalanche is located in the Red Hill valley almost 80 km west of Christchurch and is one of 42 greywacke-derived rock avalanches identified in the central Southern Alps. It overlies the Holocene active Porters Pass Fault; a component of the Porters Pass-Amberley Fault Zone which extends from the Rakaia River to beyond the Waimakariri River. The Porters Pass Fault is a dextral strike-slip fault system viewed as a series of discontinuous fault scarps. The location of the fault trace beneath the deposit suggests it may represent a possible source of seismic shaking resulting in the formation of the Acheron rock avalanche. The rock mass composition of the rock avalanche source scar is Torlesse Supergroup greywacke consisting of massive sandstone and thinly bedded mudstone sequences dipping steeply north into the centre of the source basin. A stability analysis identified potential instability along shallow north dipping planar defects, and steep south dipping toppling failure planes. The interaction of the defects with bedding is considered to have formed conditions for potential instability most likely triggered by a seismic event. The dTositional area of the rock avalanche covers 7.2 x 105 m2 with an estimated volume of 9 x 10 m3 The mobilised rock mass volume was calculated at 7.5 x 106 m3• Run out of the debris from the top of the source scar to the distal limit reached 3500m, descending over a vertical fall of almost 700m with an estimated Fahrboschung of 0.2. The run out of the rock avalanche displayed moderate to high mobility, travelling at an estimated maximum velocity of 140-160 km/hour. The rapid emplacement of the deposit is confirmed by highly fragmented internal composition and burial of forest vegetation New radiocarbon ages from buried wood retrieved from the base of Acheron rock avalanche deposit represents an emplacement age closely post-dating (Wk 12094) 1152 ± 51 years B.P. This differs significantly from a previous radiocarbon age of (NZ547) 500 ± 69 years B.P. and modal lichenometry and weathering-rind thickness ages of approximately 460 ± 10 yrs and 490 ± 50 years B.P. The new age shows no resemblance to an earthquake event around 700- 500 years B.P. on the Porters Pass-Amberley Fault Zone. The DAN run out simulation using a friction model rheology successfully replicated the long run out and velocity of the Acheron rock avalanche using a frictron angle of 27° and high earth pressure coefficients of 5.5, 5.2, and 5.9. The elevated earth pressure coefficients represent dispersive pressures derived from dynamic fragmentation of the debris within the mobile rock avalanche, supporting the hypothesis of Davies and McSaveney (2002). The DAN model has potential applications for areas prone to large-scale instability in the elevated slopes and steep waterways of the Southern Alps. A paleoseismic investigation of a newly identified scarp of the Porters Pass Fault partially buried by the rock avalanche was conducted to identify any evidence of a coseismic relationship to the Acheron rock avalanche. This identified three-four fault traces striking at 078°, and a sag pond displaying a sequence of overbank deposits containing two buried soils representing an earthquake event horizon. A 40cm vertical offset of the ponded sediment and lower buried soil horizqn was recorded, which was dated to (Wk 13112 charcoal in palosol) 653 ± 54 years B.P. and (Wk 13034 palosol) 661 ± 34 years B.P. The evidence indicates a fault rupture occurred along the Porters Pass Fault, west of Porters Pass most likely extending to the Red Lakes terraces, post-dating 700 years B.P., resulting in 40cm of vertical displacement and an unknown component of dextral strike slip movement. This event post dates the event one (1000 ± 100 years B.P) at Porters Pass previously considered to represent the most recent rupture along the fault line. This points to a probable source for resetting of the modal weathering-rind thicknesses and lichen size populations in the Red Hill valley and possibly the Red Lakes terraces. These results suggest careful consideration must be given to the geomorphic and paleoseismic history of a specific site when applying surface dating techniques and furthermore the origin of dates used in literature and their useful range should be verified. An event at 700-500 years B.P did not trigger the Acheron rock avalanche as previously assumed supporting Howard's conclusions. The lack of similar aged rupture evidence in either of the Porters Pass and Coleridge trenches supports Howard's hypothesis of segmentation of the Porters Pass Fault; where rupture occurs along one fault segment but not along another. The new rock avalanche age closely post-dating 1200-1100 years B.P. resembles the poorly constrained event one rupture age of 1700-800 years B.P for the Porters Pass Fault and the tighter constrained Round Top event of 1010 ± 50 years B.P. on the Alpine Fault. Eight other rock avalanche deposits spread across the central Southern Alps also resemble the new ages however are unable to be assigned specific earthquake events due to the large associated error bars of± 270 years. This clustering of ages does represent compelling lines of evidence for large magnitude earthquake events occurring over the central Southern Alps. The presence of a rock avalanche deposit does not signify an earthquake based on the historical evidence in the Southern Alps however clustering of ages does suggest that large Mw >7 earthquakes occurred across the Southern Alps between 1200-900 years BP.
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.
