In this paper, we perform hybrid broadband (0-10 Hz) ground motion simulations for the ten most significant events (Mw 4.7-7.1) in the 2010-2011 Canterbury earthquake sequence. Taking advantage of having repeated recordings at same stations, we validate our simulations using both recordings and an empirically-developed ground motion prediction equation (GMPE). The simulation clearly captures the sedimentary basin amplification and the rupture directivity effects. Quantitative comparisons of the simulations with both recordings and the GMPE, as well as analyses of the total residuals (indicating model bias) show that simulations perform better than the empirical GMPE, especially for long period. To scrutinize the ground motion variability, we partitioned the total residuals into different components. The total residual appears to be unbiased, and the use of a 3D velocity structure reduces the long period systematic bias particularly for stations located close to the Banks Peninsula volcanic area.
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.
Seismic isolation is an effective technology for significantly reducing damage to buildings and building contents. However, its application to light-frame wood buildings has so far been unable to overcome cost and technical barriers such as susceptibility of light-weight buildings to movement under high-wind loading. The 1994 Northridge Earthquake (6.7 MW) in the United States, 1995 Kobe Earthquake (6.9 MW) in Japan and 2011 Christchurch Earthquake (6.7 Mw) all highlighted significant loss to light-frame wood buildings with over half of earthquake recovery costs allocated to their repair and reconstruction. This poster presents a value case to highlight the benefits of seismically isolated residential buildings compared to the standard fixed-base dwellings for the Wellington region. Loss data generated by insurance claim information from the 2011 Christchurch Earthquake has been used to determine vulnerability functions for the current light-frame wood building stock. By using a simplified single degree of freedom (SDOF) building model, methods for determining vulnerability functions for seismic isolated buildings are developed. Vulnerability functions are then applied directly in a loss assessment to determine the Expected Annual Loss. Vulnerability was shown to dramatically reduce for isolated buildings compared to an equivalent fixed-base building resulting in significant monetary savings, justifying the value case. A state-of-the-art timber modelling software, Timber3D, is then used to model a typical residential building with and without seismic isolation to assess the performance of a proposed seismic isolation system which addresses the technical and cost issues.
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
At 00:02 on 14th November 2016, a Mw 7.8 earthquake occurred in and offshore of the northeast of the South Island of New Zealand. Fault rupture, ground shaking, liquefaction, and co-seismic landslides caused severe damage to distributed infrastructure, and particularly transportation networks; large segments of the country’s main highway, State Highway 1 (SH1), and the Main North Line (MNL) railway line, were damaged between Picton and Christchurch. The damage caused direct local impacts, including isolation of communities, and wider regional impacts, including disruption of supply chains. Adaptive measures have ensured immediate continued regional transport of goods and people. Air and sea transport increased quickly, both for emergency response and to ensure routine transport of goods. Road diversions have also allowed critical connections to remain operable. This effective response to regional transport challenges allowed Civil Defence Emergency Management to quickly prioritise access to isolated settlements, all of which had road access 23 days after the earthquake. However, 100 days after the earthquake, critical segments of SH1 and the MNL remain closed and their ongoing repairs are a serious national strategic, as well as local, concern. This paper presents the impacts on South Island transport infrastructure, and subsequent management through the emergency response and early recovery phases, during the first 100 days following the initial earthquake, and highlights lessons for transportation system resilience.
At 00:02 on 14th November 2016, a Mw 7.8 earthquake occurred in and offshore of the northeast of the South Island of New Zealand. Fault rupture, ground shaking, liquefaction, and co-seismic landslides caused severe damage to distributed infrastructure, and particularly transportation networks; large segments of the country’s main highway, State Highway 1 (SH1), and the Main North Line (MNL) railway line, were damaged between Picton and Christchurch. The damage caused direct local impacts, including isolation of communities, and wider regional impacts, including disruption of supply chains. Adaptive measures have ensured immediate continued regional transport of goods and people. Air and sea transport increased quickly, both for emergency response and to ensure routine transport of goods. Road diversions have also allowed critical connections to remain operable. This effective response to regional transport challenges allowed Civil Defence Emergency Management to quickly prioritise access to isolated settlements, all of which had road access 23 days after the earthquake. However, 100 days after the earthquake, critical segments of SH1 and the MNL remain closed and their ongoing repairs are a serious national strategic, as well as local, concern. This paper presents the impacts on South Island transport infrastructure, and subsequent management through the emergency response and early recovery phases, during the first 100 days following the initial earthquake, and highlights lessons for transportation system resilience.
Christchurch and Canterbury suffered significant housing losses due to the earthquakes. Estimates from the Earthquake Commission (EQC) (2011) suggest that over 150,000 homes (around three quarters of Christchurch housing stock) sustained damage from the earthquakes. Some areas of Christchurch have been declared not suitable for rebuilding, affecting more than 7,500 residential properties.
This dissertation addresses a diverse range of topics in the physics-based broadband ground motion simulation, with a focus on New Zealand applications. In particular the following topics are addressed: the methodology and computational implementation of a New Zealand Velocity Model for broadband ground motion simulation; generalised parametric functions and spatial correlations for seismic velocities in the Canterbury, New Zealand region from surface-wave-based site characterisation; and ground motion simulations of Hope Fault earthquakes. The paragraphs below outline each contribution in more detail. A necessary component in physics-based ground motion simulation is a 3D model which details the seismic velocities in the region of interest. Here a velocity model construction methodology, its computational implementation, and application in the construction of a New Zealand velocity model for use in physics-based broadband ground motion simulation are presented. The methodology utilises multiple datasets spanning different length scales, which is enabled via the use of modular sub-regions, geologic surfaces, and parametric representations of crustal velocity. A number of efficiency-related workflows to decrease the overall computational construction time are employed, while maintaining the flexibility and extensibility to incorporate additional datasets and re- fined velocity parameterizations as they become available. The model comprises explicit representations of the Canterbury, Wellington, Nelson-Tasman, Kaikoura, Marlborough, Waiau, Hanmer and Cheviot sedimentary basins embedded within a regional travel-time tomography-based velocity model for the shallow crust and provides the means to conduct ground motion simulations throughout New Zealand for the first time. Recently developed deep shear-wave velocity profiles in Canterbury enabled models that better characterise the velocity structure within geologic layers of the Canterbury sedimentary basin to be developed. Here the development of depth- and Vs30-dependent para-metric velocity and spatial correlation models to characterise shear-wave velocities within the geologic layers of the Canterbury sedimentary basin are presented. The models utilise data from 22 shear-wave velocity profiles of up to 2.5km depth (derived from surface wave analysis) juxtaposed with models which detail the three-dimensional structure of the geologic formations in the Canterbury sedimentary basin. Parametric velocity equations are presented for Fine Grained Sediments, Gravels, and Tertiary layer groupings. Spatial correlations were developed and applied to generate three-dimensional stochastic velocity perturbations. Collectively, these models enable seismic velocities to be realistically represented for applications such as 3D ground motion and site response simulations. Lastly the New Zealand velocity model is applied to simulate ground motions for a Mw7.51 rupture of the Hope Fault using a physics-based simulation methodology and a 3D crustal velocity model of New Zealand. The simulation methodology was validated for use in the region through comparison with observations for a suite of historic small magnitude earthquakes located proximal to the Hope Fault. Simulations are compared with conventionally utilised empirical ground motion models, with simulated peak ground velocities being notably higher in regions with modelled sedimentary basins. A sensitivity analysis was undertaken where the source characteristics of magnitude, stress parameter, hypocentre location and kinematic slip distribution were varied and an analysis of their effect on ground motion intensities is presented. It was found that the magnitude and stress parameter strongly influenced long and short period ground motion amplitudes, respectively. Ground motion intensities for the Hope Fault scenario are compared with the 2016 Kaikoura Mw7.8 earthquake, it was found that the Kaikoura earthquake produced stronger motions along the eastern South Island, while the Hope Fault scenario resulted in stronger motions immediately West of the near-fault region. The simulated ground motions for this scenario complement prior empirically-based estimates and are informative for mitigation and emergency planning purposes.
The Mw 7.8 Kaikōura earthquake ruptured ~200 km at the ground surface across the New Zealand plate boundary zone in the northern South Island. This study was conducted in an area of ~600 km2 in the epicentral region where the faults comprise two main non-coplanar sets that strike E-NE and NNE-NW with mainly steep dips (60о-80°). Analysis of the surface rupture using field and LiDAR data provides new information on the dimensions, geometries and kinematics of these faults which was not previously available from pre-earthquake active faults or bedrock structure. The more northerly striking fault set are sub-parallel to basement bedding and accommodated predominantly left-lateral reverse slip with net slips of ~1 and ~5 m for the Stone Jug and Leader faults, respectively. The E-NE striking Conway-Charwell and The Humps faults accrued right-lateral to oblique reverse with net slips of ~2 and ~3 m, respectively. The faults form a hard-linked system dominated by kinematics consistent with the ~260° trend of the relative plate motion vector and the transpressional structures recorded across the plate boundary in the NE South Island. Interaction and intersection of the main fault sets facilitated propagation of the earthquake and transfer of slip northwards across the plate boundary zone.
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 poster presents the on-going development of a new 3D seismic velocity model of Canterbury, New Zealand. The intention of the model is to provide the 3D crustal structure in the region at multiple length scales for seismic wave propagation simulations, both broadband ground motion and more localized shallow site response analyses.
The objective of this study is to examine the influence of near-fault motions on liquefaction triggering in Christchurch and neighboring towns during the 2010-2011 Canterbury earthquake sequence (CES). The CES began with the 4 September 2010, Mw7.1 Darfield earthquake and included up to ten events that triggered liquefaction. However, most notably, widespread liquefaction was induced by the Darfield earthquake and the Mw6.2, 22 February 2011 Christchurch earthquake. Of particular relevance to this study is the forward directivity effects that were prevalent in the motions recorded during the Darfield earthquake, and to a much lesser extent, during the Christchurch earthquake. A 2D variant of the Richart-Newmark fatigue theory was used to compute the equivalent number of cycles (neq) for the ground motions, where volumetric strain was used as the damage metric. This study is unique because it considers the contribution and phasing of both the fault-normal and fault-parallel components of motion on neq and the magnitude scaling factor (MSF). It was found that when the fault-normal and fault-parallel motions were treated individually, the former yielded a lower neq than the latter. Additionally, when the combined effects of fault-normal and fault-parallel components were considered, it was found that the MSF were higher than those commonly used. This implies that motions containing near-fault effects are less demanding on the soil than motions that do not. This may be one of several factors that resulted in less severe liquefaction occurring during the Darfield earthquake than the Christchurch earthquake.
In the aftermath of the 2010 and 2011 earthquakes, Christchurch, New Zealand is framed as a ‘transi- tional’ city, moving from its demolished past to a speculative future. The ADA Mesh Cities project asks what role media art and networks may play in the transitional city, and the practices of remembering, and reimagining space.
Tens of thousands of landslides were generated over 10, 000 km2 of North Canterbury and Marlborough as a consequence of the 14 November 2016, MW7.8 Kaikōura Earthquake. The most intense landslide damage was concentrated in 3500 km2 around the areas of fault rupture. Given the sparsely populated area affected by landslides, only a few homes were impacted and there were no recorded deaths due to landslides. Landslides caused major disruption with all road and rail links with Kaikōura being severed. The landslides affecting State Highway 1 (the main road link in the South Island of New Zealand) and the South Island main trunk railway extended from Ward in Marlborough all the way to the south of Oaro in North Canterbury. The majority of landslides occurred in two geological and geotechnically distinct materials reflective of the dominant rock types in the affected area. In the Neogene sedimentary rocks (sandstones, limestones and siltstones) of the Hurunui District, North Canterbury and around Cape Campbell in Marlborough, first-time and reactivated rock-slides and rock-block slides were the dominant landslide type. These rocks also tend to have rock material strength values in the range of 5-20 MPa. In the Torlesse 'basement' rocks (greywacke sandstones and argillite) of the Kaikōura Ranges, first-time rock and debris avalanches were the dominant landslide type. These rocks tend to have material strength values in the range of 20-50 MPa. A feature of this earthquake is the large number (more than 200) of valley blocking landslides it generated. This was partly due to the steep and confined slopes in the area and the widely distributed strong ground shaking. The largest landslide dam has an approximate volume of 12(±2) M m3 and the debris from this travelled about 2.7 km2 downslope where it formed a dam blocking the Hapuku River. The long-term stability of cracked slopes and landslide dams from future strong earthquakes and large rainstorms are an ongoing concern to central and local government agencies responsible for rebuilding homes and infrastructure. A particular concern is the potential for debris floods to affect downstream assets and infrastructure should some of the landslide dams breach catastrophically. At least twenty-one faults ruptured to the ground surface or sea floor, with these surface ruptures extending from the Emu Plain in North Canterbury to offshore of Cape Campbell in Marlborough. The mapped landslide distribution reflects the complexity of the earthquake rupture. Landslides are distributed across a broad area of intense ground shaking reflective of the elongate area affected by fault rupture, and are not clustered around the earthquake epicentre. The largest landslides triggered by the earthquake are located either on or adjacent to faults that ruptured to the ground surface. Surface faults may provide a plane of weakness or hydrological discontinuity and adversely oriented surface faults may be indicative of the location of future large landslides. Their location appears to have a strong structural geological control. Initial results from our landslide investigations suggest predictive models relying only on ground-shaking estimates underestimate the number and size of the largest landslides that occurred.
Editors: G. Williams, P. Statham, N. Brown & B. Cleland ISBN Proceedings USB: 978-1-86295-644-5 published by the University of Tasmania Blended learning plays an important role in many tertiary institutions but little has been written about the implementation of blended learning in times of adversity, natural disaster or crisis. This paper describes how, in the wake of the 22 February Canterbury earthquake, five teacher educators responded to crisis-driven changing demands and changing directions. Our narratives describe how blended learning provided students in initial teacher education programmes with some certainty and continuity during a time of civil emergency. The professional learning generated from our experiences provides valuable insights for designing and preparing for blended learning in times of crisis, as well as developing resilient blended learning programmes for the future.
The M7.8 Kaikoura Earthquake in 2016 presented a number of challenges to science agencies and institutions throughout New Zealand. The earthquake was complex, with 21 faults rupturing throughout the North Canterbury and Marlborough landscape, generating a localised seven metre tsunami and triggering thousands of landslides. With many areas isolated as a result, it presented science teams with logistical challenges as well as the need to coordinate efforts across institutional and disciplinary boundaries. Many research disciplines, from engineering and geophysics to social science, were heavily involved in the response. Coordinating these disciplines and institutions required significant effort to assist New Zealand during its most complex earthquake yet recorded. This paper explores that effort and acknowledges the successes and lessons learned by the teams involved.
This poster provides a summary of the development of a 3D shallow (z<40m) shear wave velocity (Vs) model for the urban Christchurch, New Zealand region. The model is based on a recently developed Christchurch-specific empirical correlation between Vs and cone penetration test (CPT) data (McGann et al. 2014a,b) and the large high-density database of CPT logs in the greater Christchurch urban area (> 15,000 logs as of 01/01/2014). In particular, the 3D model provides shear wave velocities for the surficial Springston Formation, Christchurch Formation, and Riccarton gravel layers which generally comprise the upper 40m in the Christchurch urban area. Point-estimates are provided on a 200m-by- 200m grid from which interpolation to other locations can be performed. This model has applications for future site characterization and numerical modeling efforts via maps of timeaveraged Vs over specific depths (e.g. Vs30, Vs10) and 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 (Lee et al. 2014) currently being developed for the purpose of broadband ground motion simulation.
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.
The Canterbury Earthquake Sequence (CES) of 2010-2011 caused widespread liquefaction in many parts of Christchurch. Observations from the CES highlight some sites were liquefaction was predicted by the simplified method but did not manifest. There are a number of reasons why the simplified method may over-predict liquefaction, one of these is the dynamic interaction between soil layers within a stratified deposit. Soil layer interaction occurs through two key mechanisms; modification of the ground motion due to seismic waves passing through deep liquefied layers, and the effect of pore water seepage from an area of high excess pore water pressure to the surrounding soil. In this way, soil layer interaction can significantly alter the liquefaction behaviour and surface manifestation of soils subject to seismic loading. This research aimed to develop an understanding of how soil layer interaction, in particular ground motion modification, affects the development of excess pore water pressures and liquefaction manifestation in a soil deposit subject to seismic loading. A 1-D soil column time history Effective Stress Analysis (ESA) was conducted to give an in depth assessment of the development of pore pressures in a number of soil deposits. For this analysis, ground motions, soil profiles and model parameters were required for the ESA. Deconvolution of ground motions recorded at the surface during the CES was used to develop some acceleration time histories to input at the base of the soil-column model. An analysis of 55 sites around Christchurch, where detailed site investigations have been carried out, was then conducted to identify some simplified soil profiles and soil characteristics. From this analysis, four soil profiles representative of different levels of liquefaction manifestation were developed. These were; two thick uniform and vertically continuous sandy deposits that were representative of sites were liquefaction manifested in both the Mw 7.1 September 2010 and the Mw 6.3 February 2011 earthquakes, and two vertically discontinuous profiles with interlayered liquefiable and non-liquefiable layers representative of sites that did not manifest liquefaction in either the September 2010 or the February 2011 events. Model parameters were then developed for these four representative soil profiles through calibration of the constitutive model in element test simulations. Simulations were run for each of the four profiles subject to three levels of loading intensity. The results were analysed for the effect of soil layer interaction. These were then compared to a simplified triggering analysis for the same four profiles to determine where the simplified method was accurate in predicting soil liquefaction (for the continuous sandy deposits) and were it was less accurate (the vertically discontinuous deposits where soil layer interaction was a factor).
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.
A magnitude 6.3 earthquake struck the city of Christchurch at 12:51pm on Tuesday 22 February 2011. The earthquake caused 182 fatalities, a large number of injuries, and resulted in widespread damage to the built environment, including significant disruption to the lifelines. The event created the largest lifeline disruption in a New Zealand city in 80 years, with much of the damage resulting from extensive and severe liquefaction in the Christchurch urban area. The Christchurch earthquake occurred when the Canterbury region and its lifelines systems were at the early stage of recovering from the 4 September 2010 Darfield (Canterbury) magnitude 7.1 earthquake. This paper describes the impact of the Christchurch earthquake on lifelines by briefly summarising the physical damage to the networks, the system performance and the operational response during the emergency management and the recovery phase. Special focus is given to the performance and management of the gas, electric and road networks and to the liquefaction ejecta clean-up operations that contributed to the rapid reinstatement of the functionality of many of the lifelines. The water and wastewater system performances are also summarized. Elements of resilience that contributed to good network performance or to efficient emergency and recovery management are highlighted in the paper.
In most design codes, infill walls are considered as non-structural elements and thus are typically neglected in the design process. The observations made after major earthquakes (Duzce 1999, L’Aquila 2009, Christchurch 2011) have shown that even though infill walls are considered to be non-structural elements, they interact with the structural system during seismic actions. In the case of heavy infill walls (i.e. clay brick infill walls), the whole behaviour of the structure may be affected by this interaction (i.e. local or global structural failures such as soft storey mechanism). In the case of light infill walls (i.e. non-structural drywalls), this may cause significant economical losses. To consider the interaction of the structural system with the ‘non-structural ’infill walls at design stage may not be a practical approach due to the complexity of the infill wall behaviour. Therefore, the purpose of the reported research is to develop innovative technological solutions and design recommendations for low damage non-structural wall systems for seismic actions by making use of alternative approaches. Light (steel/timber framed drywalls) and heavy (unreinforced clay brick) non-structural infill wall systems were studied by following an experimental/numerical research programme. Quasi-static reverse cyclic tests were carried out by utilizing a specially designed full scale reinforced concrete frame, which can be used as a re-usable bare frame. In this frame, two RC beams and two RC columns were connected by two un-bonded post tensioning bars, emulating a jointed ductile frame system (PRESSS technology). Due to the rocking behaviour at the beam-column joint interfaces, this frame was typically a low damage structural solution, with the post-tensioning guaranteeing a linear elastic behaviour. Therefore, this frame could be repeatedly used in all of the tests carried out by changing only the infill walls within this frame. Due to the linear elastic behaviour of this structural bare frame, it was possible to extract the exact behaviour of the infill walls from the global results. In other words, the only parameter that affected the global results was given by the infill walls. For the test specimens, the existing practice of construction (as built) for both light and heavy non-structural walls was implemented. In the light of the observations taken during these tests, modified low damage construction practices were proposed and tested. In total, seven tests were carried out: 1) Bare frame , in order to confirm its linear elastic behaviour. 2) As built steel framed drywall specimen FIF1-STFD (Light) 3) As built timber framed drywall specimen FIF2-TBFD (Light) 4) As built unreinforced clay brick infill wall specimen FIF3-UCBI (Heavy) 5) Low damage steel framed drywall specimen MIF1-STFD (Light) 6) Low damage timber framed drywall specimen MIF2-TBFD (Light) 7) Low damage unreinforced clay brick infill wall specimen MIF5-UCBI (Heavy) The tests of the as built practices showed that both drywalls and unreinforced clay brick infill walls have a low serviceability inter-storey drift limit (0.2-0.3%). Based on the observations, simple modifications and details were proposed for the low damage specimens. The details proved to be working effectively in lowering the damage and increasing the serviceability drift limits. For drywalls, the proposed low damage solutions do not introduce additional cost, material or labour and they are easily applicable in real buildings. For unreinforced clay brick infill walls, a light steel sub-frame system was suggested that divides the infill panel zone into smaller individual panels, which requires additional labour and some cost. However, both systems can be engineered for seismic actions and their behaviour can be controlled by implementing the proposed details. The performance of the developed details were also confirmed by the numerical case study analyses carried out using Ruaumoko 2D on a reinforced concrete building model designed according to the NZ codes/standards. The results have confirmed that the implementation of the proposed low damage solutions is expected to significantly reduce the non-structural infill wall damage throughout a building.
The magnitude Mw7.8 ‘Kaikōura’ earthquake occurred shortly after midnight on 14 November 2016. This paper presents an overview of the geotechnical impacts on the South Island of New Zealand recorded during the postevent reconnaissance. Despite the large moment magnitude of this earthquake, relatively little liquefaction was observed across the South Island, with the only severe manifestation occurring in the young, loose alluvial deposits in the floodplains of the Wairau and Opaoa Rivers near Blenheim. The spatial extent and volume of liquefaction ejecta across South Island is significantly less than that observed in Christchurch during the 2010-2011 Canterbury Earthquake Sequence, and the impact of its occurrence to the built environment was largely negligible on account of the severe manifestations occurring away from the areas of major development. Large localised lateral displacements occurred in Kaikōura around Lyell Creek. The soft fine-grained material in the upper portions of the soil profile and the free face at the creek channel were responsible for the accumulation of displacement during the ground shaking. These movements had severely impacted the houses which were built close (within the zone of large displacement) to Lyell Creek. The wastewater treatment facility located just north of Kaikōura also suffered tears in the liners of the oxidation ponds and distortions in the aeration system due to ground movements. Ground failures on the Amuri and Emu Plains (within the Waiau Valley) were small considering the large peak accelerations (in excess of 1g) experienced in the area. Minor to moderate lateral spreading and ejecta was observed at some bridge crossings in the area. However, most of the structural damage sustained by the bridges was a result of the inertial loading, and the damage resulting from geotechnical issues were secondary.
In the period between September 2010 and December 2011, Christchurch (New Zealand) and its surroundings were hit by a series of strong earthquakes including six significant events, all generated by local faults in proximity to the city: 4 September 2010 (Mw=7.1), 22 February 2011 (Mw=6.2), 13 June 2011 (Mw=5.3 and Mw=6.0) and 23 December 2011 (M=5.8 and (M=5.9) earthquakes. As shown in Figure 1, the causative faults of the earthquakes were very close to or within the city boundaries thus generating very strong ground motions and causing tremendous damage throughout the city. Christchurch is shown as a lighter colour area, and its Central Business District (CBD) is marked with a white square area in the figure. Note that the sequence of earthquakes started to the west of the city and then propagated to the south, south-east and east of the city through a set of separate but apparently interacting faults. Because of their strength and proximity to the city, the earthquakes caused tremendous physical damage and impacts on the people, natural and built environments of Christchurch. The 22 February 2011 earthquake was particularly devastating. The ground motions generated by this earthquake were intense and in many parts of Christchurch substantially above the ground motions used to design the buildings in Christchurch. The earthquake caused 182 fatalities, collapse of two multi-storey reinforced concrete buildings, collapse or partial collapse of many unreinforced masonry structures including the historic Christchurch Cathedral. The Central Business District (CBD) of Christchurch, which is the central heart of the city just east of Hagley Park, was practically lost with majority of its 3,000 buildings being damaged beyond repair. Widespread liquefaction in the suburbs of Christchurch, as well as rock falls and slope/cliff instabilities in the Port Hills affected tens of thousands of residential buildings and properties, and shattered the lifelines and infrastructure over approximately one third of the city area. The total economic loss caused by the 2010-2011 Christchurch earthquakes is currently estimated to be in the range between 25 and 30 billion NZ dollars (or 15% to 18% of New Zealand’s GDP). After each major earthquake, comprehensive field investigations and inspections were conducted to document the liquefaction-induced land damage, lateral spreading displacements and their impacts on buildings and infrastructure. In addition, the ground motions produced by the earthquakes were recorded by approximately 15 strong motion stations within (close to) the city boundaries providing and impressive wealth of data, records and observations of the performance of ground and various types of structures during this unusual sequence of strong local earthquakes affecting a city. This paper discusses the liquefaction in residential areas and focuses on its impacts on dwellings (residential houses) and potable water system in the Christchurch suburbs. The ground conditions of Christchurch including the depositional history of soils, their composition, age and groundwater regime are first discussed. Detailed liquefaction maps illustrating the extent and severity of liquefaction across Christchurch triggered by the sequence of earthquakes including multiple episodes of severe re-liquefaction are next presented. Characteristic liquefaction-induced damage to residential houses is then described focussing on the performance of typical house foundations in areas affected by liquefaction. Liquefaction impacts on the potable water system of Christchurch is also briefly summarized including correlation between the damage to the system, liquefaction severity, and the performance of different pipe materials. Finally, the characteristics of Christchurch liquefaction and its impacts on built environment are discussed in relation to the liquefaction-induced damage in Japan during the 11 March 2011 Great East Japan Earthquake.
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.
