During the 2010 - 2011 Canterbury earthquake sequence, extensive liquefaction was observed in many areas of Christchurch city and its surroundings, causing widespread damage to buildings and infrastructure. While existing simplified methods were found to work well in some areas of the city, there were also large areas where these methods did not perform satisfactorily. In some of these cases, researchers have proposed that layers of fine grained material within the soil profile may be responsible for preventing the manifestation of liquefaction. This paper presents preliminary findings on the mechanisms at play when pressure differentials exist across a clay layer. It is found that if the clay layer is unable to distort, then pore fluid is unable to break-through the layer even with relatively high pressures, resulting in dissipation of excess pore pressures by seepage. If the layers are however able to distort, then it is possible for the pore fluid to break through the clay layer, potentially resulting in adverse effects in terms of the severity of liquefaction.
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.
As a consequence of the 2010 – 2011 Canterbury earthquake sequence, Christchurch experienced widespread liquefaction, vertical settlement and lateral spreading. These geological processes caused extensive damage to both housing and infrastructure, and increased the need for geotechnical investigation substantially. Cone Penetration Testing (CPT) has become the most common method for liquefaction assessment in Christchurch, and issues have been identified with the soil behaviour type, liquefaction potential and vertical settlement estimates, particularly in the north-western suburbs of Christchurch where soils consist mostly of silts, clayey silts and silty clays. The CPT soil behaviour type often appears to over-estimate the fines content within a soil, while the liquefaction potential and vertical settlement are often calculated higher than those measured after the Canterbury earthquake sequence. To investigate these issues, laboratory work was carried out on three adjacent CPT/borehole pairs from the Groynes Park subdivision in northern Christchurch. Boreholes were logged according to NZGS standards, separated into stratigraphic layers, and laboratory tests were conducted on representative samples. Comparison of these results with the CPT soil behaviour types provided valuable information, where 62% of soils on average were specified by the CPT at the Groynes Park subdivision as finer than what was actually present, 20% of soils on average were specified as coarser than what was actually present, and only 18% of soils on average were correctly classified by the CPT. Hence the CPT soil behaviour type is not accurately describing the stratigraphic profile at the Groynes Park subdivision, and it is understood that this is also the case in much of northwest Christchurch where similar soils are found. The computer software CLiq, by GeoLogismiki, uses assessment parameter constants which are able to be adjusted with each CPT file, in an attempt to make each more accurate. These parameter changes can in some cases substantially alter the results for liquefaction analysis. The sensitivity of the overall assessment method, raising and lowering the water table, lowering the soil behaviour type index, Ic, liquefaction cutoff value, the layer detection option, and the weighting factor option, were analysed by comparison with a set of ‘base settings’. The investigation confirmed that liquefaction analysis results can be very sensitive to the parameters selected, and demonstrated the dependency of the soil behaviour type on the soil behaviour type index, as the tested assessment parameters made very little to no changes to the soil behaviour type plots. The soil behaviour type index, Ic, developed by Robertson and Wride (1998) has been used to define a soil’s behaviour type, which is defined according to a set of numerical boundaries. In addition to this, the liquefaction cutoff point is defined as Ic > 2.6, whereby it is assumed that any soils with an Ic value above this will not liquefy due to clay-like tendencies (Robertson and Wride, 1998). The method has been identified in this thesis as being potentially unsuitable for some areas of Christchurch as it was developed for mostly sandy soils. An alternative methodology involving adjustment of the Robertson and Wride (1998) soil behaviour type boundaries is proposed as follows: Ic < 1.31 – Gravelly sand to dense sand 1.31 < Ic < 1.90 – Sands: clean sand to silty sand 1.90 < Ic < 2.50 – Sand mixtures: silty sand to sandy silt 2.50 < Ic < 3.20 – Silt mixtures: clayey silt to silty clay 3.20 < Ic < 3.60 – Clays: silty clay to clay Ic > 3.60 – Organics soils: peats. When the soil behaviour type boundary changes were applied to 15 test sites throughout Christchurch, 67% showed an improved change of soil behaviour type, while the remaining 33% remained unchanged, because they consisted almost entirely of sand. Within these boundary changes, the liquefaction cutoff point was moved from Ic > 2.6 to Ic > 2.5 and altered the liquefaction potential and vertical settlement to more realistic ii values. This confirmed that the overall soil behaviour type boundary changes appear to solve both the soil behaviour type issues and reduce the overestimation of liquefaction potential and vertical settlement. This thesis acts as a starting point towards researching the issues discussed. In particular, future work which would be useful includes investigation of the CLiq assessment parameter adjustments, and those which would be most suitable for use in clay-rich soils such as those in Christchurch. In particular consideration of how the water table can be better assessed when perched layers of water exist, with the limitation that only one elevation can be entered into CLiq. Additionally, a useful investigation would be a comparison of the known liquefaction and settlements from the Canterbury earthquake sequence with the liquefaction and settlement potentials calculated in CLiq for equivalent shaking conditions. This would enable the difference between the two to be accurately defined, and a suitable adjustment applied. Finally, inconsistencies between the Laser-Sizer and Hydrometer should be investigated, as the Laser-Sizer under-estimated the fines content by up to one third of the Hydrometer values.
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.
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.
Following the 22nd February 2011, Mw 6.2 earthquake located along a previously unknown fault beneath the Port Hills of Christchurch, surface cracking was identified in contour parallel locations within fill material at Quarry Road on the lower slopes of Mount Pleasant. GNS Science, in the role of advisor to the Christchurch City Council, concluded that these cracks were a part of a potential rotational mass movement (named zone 11A) within the fill and airfall loess material present. However, a lack of field evidence for slope instability and an absence of laboratory geotechnical data on which slope stability analysis was based, suggested this conclusion is potentially incorrect. It was hypothesised that ground cracking was in fact due to earthquake shaking, and not mass movement within the slope, thus forming the basis of this study. Three soil units were identified during surface and subsurface investigations at Quarry Road: fill derived from quarry operations in the adjacent St. Andrews Quarry (between 1893 and 1913), a buried topsoil, and underlying in-situ airfall loess. The fill material was identified by the presence of organic-rich topsoil “clods” that were irregular in both size (∼10 – 200 mm) and shape, with variable thicknesses of 1 – 10 m. Maximum thickness, as indicated by drill holes and geophysical survey lines, was identified below 6 Quarry Road and 7 The Brae where it is thought to infill a pre-existing gully formed in the underlying airfall loess. Bearing strength of the fill consistently exceeded 300 kPa ultimate below ∼500 mm depth. The buried topsoil was 200 – 300 mm thick, and normally displayed a lower bearing strength when encountered, but not below 300 kPa ultimate (3 – 11 blows per 100mm or ≥100 kPa allowable). In-situ airfall loess stood vertically in outcrop due to its characteristic high dry strength and also showed Scala penetrometer values of 6 – 20+ blows per 100 mm (450 – ≥1000 kPa ultimate). All soils were described as being moist to dry during subsurface investigations, with no groundwater table identified during any investigation into volcanic bedrock. In-situ moisture contents were established using bulk disturbed samples from hand augers and test pitting. Average moisture contents were low at 9% within the fill, 11 % within the buried topsoil, and 8% within the airfall loess: all were below the associated average plastic limit of 17, 15, and 16, respectively, determined during Atterberg limit analysis. Particle size distributions, identified using the sieve and pipette method, were similar between the three soil units with 11 – 20 % clay, 62 – 78 % silt, and 11 – 20 % fine sand. Using these results and the NZGS soil classification, the loess derived fill and in-situ airfall loess are termed SILT with some clay and sand, and the buried topsoil is SILT with minor clay and sand. Dispersivity of the units was found using the Emerson crumb test, which established that the fill can be non- to completely dispersive (score 0 – 4). The buried topsoil was always non-dispersive (score 0), and airfall loess completely dispersive (score 4). Values for cohesion (c) and internal friction angle (φ) of the three soil units were established using the direct shear box at field moisture contents. Results showed all soil units had high shear strengths at the moisture contents tested (c = 18 – 24 kPa and φ = 42 – 50°), with samples behaving in a brittle fashion. Moisture content was artificially increased to 16% within the buried topsoil, which reduced the shear strength (c = 10 kPa, φ = 18°) and allowed it to behave plastically. Observational information indicating stability at Quarry Road included: shallow, discontinuous, cracks that do not display vertical offset; no scarp features or compressional zones typical of landsliding; no tilted or deformed structures; no movement in inclinometers; no basal shear zone identified in logged core to 20 m depth; low field moisture contents; no groundwater table; and high soil strength using Scala penetrometers. Limit equilibrium analysis of the slope was conducted using Rocscience software Slide 5.0 to verify the slope stability identified by observational methods. Friction, cohesion, and density values determined during laboratory were input into the two slope models investigated. Results gave minimum static factor of safety values for translational (along buried topsoil) and rotational (in the fill) slides of 2.4 – 4.2. Sensitivity of the slope to reduced shear strength parameters was analysed using c = 10 kPa and φ = 18° for the translational buried topsoil plane, and a cohesion of 0 kPa within the fill for the rotational plane. The only situation that gave a factor of safety <1.0 was in nonengineered fill at 0.5 m depth. Pseudostatic analysis based on previous peak ground acceleration (PGA) values for the Canterbury Earthquake Sequence, and predicted PGAs for future Alpine Fault and Hope Fault earthquakes established minimum factor of safety values between 1.2 and 3.3. Yield acceleration PGAs were computed to be between 0.8g and 1.6g. Based on all information gathered, the cracking at Quarry Road is considered to be shallow deformation in response to earthquake shaking, and not due to deep-seated landsliding. It is recommended that the currently bare site be managed by smoothing the land, installing contour drainage, and bioremediation of the surface soils to reduce surface water infiltration and runoff. Extensive earthworks, including removal of the fill, are considered unnecessary. Any future replacement of housing would be subject to site-specific investigations, and careful foundation design based on those results.
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.
Fine grained sediment deposition in urban environments during natural hazard events can impact critical infrastructure and properties (urban terrain) leading to reduced social and economic function and potentially adverse public health effects. Therefore, clean-up of the sediments is required to minimise impacts and restore social and economic functionality as soon as possible. The strategies employed to manage and coordinate the clean-up significantly influence the speed, cost and quality of the clean-up operation. Additionally, the physical properties of the fine grained sediment affects the clean-up, transport, storage and future usage of the sediment. The goals of the research are to assess the resources, time and cost required for fine grained sediment clean-up in an urban environment following a disaster and to determine how the geotechnical properties of sediment will affect urban clean-up strategies. The thesis focuses on the impact of fine grained sediment (<1 mm) deposition from three liquefaction events during the Canterbury earthquake sequence (2010-2011) on residential suburbs and transport networks in Christchurch. It also presents how geotechnical properties of the material may affect clean-up strategies and methods by presenting geotechnical analysis of tephra material from the North Island of New Zealand. Finally, lessons for disaster response planning and decision making for clean-up of sediment in urban environments are presented. A series of semi-structured interviews of key stakeholders supported by relevant academic literature and media reports were used to record the clean-up operation coordination and management and to make a preliminary qualification of the Christchurch liquefaction ejecta clean-up (costs breakdown, time, volume, resources, coordination, planning and priorities). Further analysis of the costs and resources involved for better accuracy was required and so the analysis of Christchurch City Council road management database (RAMM) was done. In order to make a transition from general fine sediment clean-up to specific types of fine disaster sediment clean-up, adequate information about the material properties is required as they will define how the material will be handled, transported and stored. Laboratory analysis of young volcanic tephra from the New Zealand’s North Island was performed to identify their geotechnical properties (density, granulometry, plasticity, composition and angle of repose). The major findings of this research were that emergency planning and the use of the coordinated incident management system (CIMS) system during the emergency were important to facilitate rapid clean-up tasking, management of resources and ultimately recovery from widespread and voluminous liquefaction ejecta deposition in eastern Christchurch. A total estimated cost of approximately $NZ 40 million was calculated for the Christchurch City clean-up following the 2010-2011 Canterbury earthquake sequence with a partial cost of $NZ 12 million for the Southern part of the city, where up to 33% (418 km) of the road network was impacted by liquefaction ejecta and required clearing of the material following the 22 February 2011 earthquake. Over 500,000 tonnes of ejecta has been stockpiled at Burwood landfill for all three liquefaction inducing earthquake events. The average cost per kilometre for the event clean-up was $NZ 5,500/km (4 September 2010), $NZ 11,650/km (22 February 2011) and $NZ 11,185/km (13 June 2011). The duration of clean-up time of residential properties and the road network was approximately two to three months for each of the three liquefaction ejecta events; despite events volumes and spatial distribution of ejecta. Interviews and quantitative analysis of RAMM data revealed that the experience and knowledge gained from the Darfield earthquake (4 September 2010) clean-up increased the efficiency of the following Christchurch earthquake induced liquefaction ejecta clean-up events. Density, particle size, particle shape, clay content and moisture content, are the important geotechnical properties that need to be considered when planning for a clean-up method that incorporates collection, transport and disposal or storage. The geotechnical properties for the tephra samples were analysed to increase preparedness and reaction response of potentially affected North Island cities from possible product from the active volcanoes in their region. The geotechnical results from this study show that volcanic tephra could be used in road or construction material but the properties would have to be further investigated for a New Zealand context. Using fresh volcanic material in road, building or flood control construction requires good understanding of the material properties and precaution during design and construction to extra care, but if well planned, it can be economically beneficial.
One of the less understood geotechnical responses to the cyclic loading from the MW6.2 Christchurch Earthquake, on the 22nd of February 2011, is the fissuring in the loessial soil-mantled, footslope positions of the north-facing valleys of the Port Hills. The fissures are characterized by mostly horizontal offset (≤500mm), with minor vertical displacement (≤300mm), and they extend along both sides of valleys for several hundred metres in an approximately contour-parallel orientation. The fissure traces correspond to extensional features mapped in other studies. Previous studies have suggested that the fissures are the headscarps of incipient landslides, but the surface and subsurface features are not typical of landslide movement. Whilst there are some features that correlate with landslide movement, there are many features that contradict the landslide movement hypothesis. Of critical importance to this investigation was the fact that there are no landslide flanks, there has been no basal shear surface found, there is little deformation in the so-called ‘landslide body’, and there have been no recorded zones of low shear strength in the soil deposit that are indicative of a basal shear surface. This thesis is a detailed geotechnical study on the fissures along part of Ramahana Road in the Hillsborough Valley, Christchurch. Shallow and deep investigation methods found that the predominant soil is loess-colluvium, to depths of ~20m, and this soil has variable geotechnical characteristics depending on the layer sampled. The factor that has the most influence on shear strength was found to be the moisture content. Direct shear-box testing of disturbed, recompacted loess-colluvium found that the soil had a cohesion of 35-65kPa and a friction angle of 38-43° when the soil moisture content was at 8-10%. However when the moisture content was at 19-20% the soil’s cohesion decreased to 3-5kPa and its friction angle decreased to 33-38°, this moisture content is at or slightly above the plastic limit. An electrical resistivity geophysical survey was conducted perpendicular to multiple fissure traces and through the compressional zone at 17 Ramahana Road. The electrical resistivity line found that there was an area of high resistivity at the toe of the slope, and an area of high conductivity downslope of this and at greater depths. This area correlated to the compressional zone recorded by previous studies. Moisture content testing of the soil in these locations showed that the soil in the resistive area was relatively dry (9%) compared to the surrounding soil (13%), whilst the soil in the conductive area was relatively wet (22%)compared to the surrounding soil (19%). Density tests of the soil in the compressional zone recorded that the resistive area had a higher dry density than the surrounding soil (~1790 kg/m3 compared to ~1650 kg/m3). New springs arose downslope of the compressional zone contemporaneously with the fissures, and it is interpreted that these have arisen from increased hydraulic head in the Banks Peninsula bedrock aquifer system, and earthquake induced-bedrock fracturing. A test pit was dug across an infilled fissure trace at 17 Ramahana Road to a depth of 3m. The fissure trace had an aperture of 450-470mm at the ground surface, but it gradually lost aperture with depth until 2.0-2.1m where it became a segmented fissure trace with 1-2mm aperture. A mixed-colluvium layer was intercepted by the fissure trace at 2.4m depth, and there was no observable vertical offset of this layer. The fissure trace was at an angle of 78° at the ground surface, but it also flattened with depth, which gave it a slightly curved appearance. The fissure trace was at an assumed angle of 40-50° near the base of the test pit. Rotational slide, translational slide and lateral spread landslide movement types were compared and contrasted as possibilities for landslide movement types, whilst an alternative hypothesis was offered that the fissures are tensile failures with a quasi-toppling motion involving a cohesive block of loessial soil moving outwards from the slope, with an accommodating compressional strain in the lower less cohesive soil. The mechanisms behind this movement are suggested to be the horizontal earthquake inertia forces from the Christchurch Earthquake, the static shear stress of the slope, and bedrock uplift of the Port Hills in relation to the subsidence of the Christchurch city flatlands. Extremely high PGA is considered to be a prerequisite to the fissure trace development, and these can only be induced in the Hillsborough Valley from a Port Hills Fault rupture, which has a recurrence interval of ~10,000 years. The current understanding of how the loess-colluvium soil would behave under cyclic loading is limited, and the mechanisms behind the suggested movement type are not completely understood. Further research is needed to confirm the proposed mechanism of the fissure traces. Laboratory tests such as the cyclic triaxial and cyclic shear test would be beneficial in future research to quantitatively test how the soil behaves under cyclic loading at various moisture contents and clay contents, and centrifuge experiments would be of great use to qualitatively test the suggested mode of movement in the loessial soil.
This project was initiated by ENGEO Limited and KiwiRail Holdings Limited to assess the stability of Slovens Creek Viaduct (specifically its western abutment) and a 3km section of rail corridor between Slovens Creek Viaduct and Avoca on the Midland Line (MDL). Commonly known as the scenic TranzAlpine rail journey (through Arthurs Pass National Park) the MDL connects Greymouth to Christchurch via Rolleston, where the MDL meets the Main South Line into Christchurch. The project area is approximately 40km southeast of Arthurs Pass Township, in the eastern extension of the Castle Hill Basin which is part of the Waimakariri Catchment and Canterbury Foothills. The field area is underlain by Rakaia Terrane, which is part of the Torlesse Composite Terrane forming the basement rock unit for the field area. Cretaceous-Tertiary rocks of the Castle Hill Basin overlie the basement strata and record a transgression-regression sequence, as well as mid-Oligocene submarine volcanism. The stratigraphic sequence in the Castle Hill Basin, and its eastern extension to Avoca, comprises two formations of the Eyre group, the older Broken River Formation and the younger Iron Creek Formation. Deep marine Porter Group limestones, marls, and tuffs of Oligocene age succeed the Iron Creek Formation of the Eyre Group, and probably records the maximum of the transgression. The Enys Formation lies disconformably on the Porter Group and is overlain unconformably by Late Pleistocene glacifluvial and glacial deposits. The Tertiary strata in the Slovens-Avoca rail corridor are weak, and the clay-rich tuff derived from mid-Oligocene volcanism is particularly prone to slaking. Extensive mapping carried out for this project has identified that some 90 percent of the surface along the length of the Slovens-Avoca corridor has been subject to mass movement. The landslides of the Slovens-Avoca rail corridor are clearly younger than the Last Glaciation, and Slovens Creek has been downcutting, with associated faulting and uplift, to form the present day geomorphology of the rail corridor. Deep-seated landslides in the rail corridor extend to Slovens Creek, locally deflecting the stream course, and a generic ground failure model for the rail corridor has been developed. Exploratory geotechnical investigations, including core drilling, installation of an inclinometer and a piezometer, enabled the construction of a simple ground model and cross section for the Slovens Creek Viaduct western abutment. Limit-equilibrium and pseudo-static slope stability analyses using both circular and block critical slip surface search methods were applied to the ground model for the western abutment of Slovens Creek Viaduct. Piezometric and strength data obtained during laboratory testing of core material have been used to constrain the western abutment stability assessment for one representative section line (C-C’). Prior to pseudo-static sensitivity analyses peak ground acceleration (PGA) for various Ultimate Limit State (ULS) design return periods, defined by an equation given in NZS1170.5:2004, were calculated and have been used as a calibration technique to find and compare specific PGA values for pseudo-static analyses in the Slovens Creek Viaduct area. The main purpose has been to provide an indication of how railway infrastructure could be affected by seismic events of various return periods defined by ULS design standards for the area. Limit equilibrium circular slip surface search methods, both grid search and auto refine search, indicated the slope is stable with a FoS greater than 1.0 returned from each, although one particular surface returned the lowest FoS in each. This surface is in the lower portion of the slope, adjacent to Slovens Stream and northeast of the MDL. As expected, pseudo-static analyses returned a lower FoS overall when compared to limit equilibrium analyses. The PGA analyses suggest that partial ground failure at the Slovens Creek Viaduct western abutment could occur in a 1 in 25-year return period event within materials on the slower slope beyond the immediate rail corridor. A ULS (1 in 500-year) event in the Slovens Creek Viaduct area would likely produce a PGA of ~0.9g, and the effects on the western abutment and rail infrastructure would most likely be catastrophic. Observed ground conditions for the western abutment of the Slovens Creek Viaduct suggest there is no movement within the landslide at depth within the monitoring timeframe of this project (22 May 2015 – 4 August 2015). Slope stability monitoring is recommended to be continued in two parts: (1) the inclinometer in BH1 is to be monitored on a six monthly basis for one year following completion of this thesis, and then annually unless ground movements become evident; and (2) surface movement monitoring should be installed using a fixed datum on the stable eastern abutment. Long-term stability management strategies for the Slovens Creek Viaduct western abutment are dependent upon future observed changes and ongoing monitoring. Hazard and risk assessment using the KiwiRail Qualitative Risk Assessment Framework (QRA) is recommended, and if slope stability becomes problematic for operation of the Midland Line consideration should be given to deep slope drainage. In the event of a large magnitude or high PGA earthquake all monitoring should be reviewed.
