A non-destructive hardness testing method has been developed to investigate the amount of plastic strain demand in steel elements subjected to cyclic loading. The focus of this research is on application to the active links of eccentrically braced frames (EBFs), which are a commonly used seismic-resisting system in modern steel framed buildings. The 2010/2011 Christchurch earthquake series, especially the very intense February 22 shaking, which was the first earthquake worldwide to push complete EBF systems fully into their inelastic state, generating a moderate to high level of plastic strain in EBF active links, for a range of buildings from 3 to 23 storeys in height. This raised two important questions: 1) what was the extent of plastic deformation in active links; and 2) what effect does that have to post-earthquake steel properties? This project comprised determining a robust relationship between hardness and plastic strain in order to be able to answer the first question and provide the necessary input into answering the second question. A non-destructive Leeb (portable) hardness tester (model TH170) has been used to measure the hardness, in order to determine the plastic strain, in hot rolled steel universal sections and steel plates. A bench top Rockwell B was used to compare and validated the hardness measured by the portable hardness tester. Hardness was measured from monotonically strained tensile test specimens to identify the relationship between hardness and plastic strain demand. Test results confirmed a good relationship between hardness and the amount of monotonically induced plastic strain. Surface roughness was identified as an important parameter in obtaining reliable hardness readings from a portable hardness reader. A proper surface preparation method was established by using three different cleaning methods, finished with hand sanding to achieve surface roughness coefficients sufficiently low not to distort the results. This work showed that a test surface roughness (Ra) is not more than 1.6 micron meter (μm) is required for accurate readings from the TH170 tester. A case study on an earthquake affected building was carried out to identify the relationship between hardness and amount of plastic strain demand in cyclically deformed active links. Hardness was carried out from active links shown visually to have been the most affected during one of the major earthquake events. Onsite hardness test results were then compared with laboratory hardness test results. A good relationship between hardness from onsite and laboratory was observed between the test methods; Rockwell B bench top and portable Leeb tester TH170. Manufacturing induced plastic strain in the top and bottom of the webs of hot rolled sections were discovered from this research, an important result which explains why visual effects of earthquake induced active link yielding (eg cracked or flaking paint) was typically more prevalent over the middle half depth of the active link. The extent of this was quantified. It was also evident that the hardness readings from the portable hardness tester are influenced by geometry, mass effects and rigidity of the links. The final experimental stage was application of the method to full scale cyclic inelastic tested nominally identical active links subjected to loading regimes comprising constant and variable plastic strain demands. The links were cyclically loaded to achieve different plastic strain level. A novel Digital Image Correlation (DIC) technique was incorporated during the tests of this scale, to confirm the level of plastic strain achieved. Tensile test specimens were water jet cut from cyclically deformed webs to analyse the level of plastic strain. Test results show clear evidence that cyclically deformed structural steel elements show good correlation between hardness and the amount of plastic strain demand. DIC method was found to be reliable and accurate to check the level of plastic strain within cyclically deformed structural steel elements.
Sewerage systems convey sewage, or wastewater, from residential or commercial buildings through complex reticulation networks to treatment plants. During seismic events both transient ground motion and permanent ground deformation can induce physical damage to sewerage system components, limiting or impeding the operability of the whole system. The malfunction of municipal sewerage systems can result in the pollution of nearby waterways through discharge of untreated sewage, pose a public health threat by preventing the use of appropriate sanitation facilities, and cause serious inconvenience for rescuers and residents. Christchurch, the second largest city in New Zealand, was seriously affected by the Canterbury Earthquake Sequence (CES) in 2010-2011. The CES imposed widespread damage to the Christchurch sewerage system (CSS), causing a significant loss of functionality and serviceability to the system. The Christchurch City Council (CCC) relied heavily on temporary sewerage services for several months following the CES. The temporary services were supported by use of chemical and portable toilets to supplement the damaged wastewater system. The rebuild delivery agency -Stronger Christchurch Infrastructure Rebuild Team (SCIRT) was created to be responsible for repair of 85 % of the damaged horizontal infrastructure (i.e., water, wastewater, stormwater systems, and roads) in Christchurch. Numerous initiatives to create platforms/tools aiming to, on the one hand, support the understanding, management and mitigation of seismic risk for infrastructure prior to disasters, and on the other hand, to support the decision-making for post-disaster reconstruction and recovery, have been promoted worldwide. Despite this, the CES in New Zealand highlighted that none of the existing platforms/tools are either accessible and/or readable or usable by emergency managers and decision makers for restoring the CSS. Furthermore, the majority of existing tools have a sole focus on the engineering perspective, while the holistic process of formulating recovery decisions is based on system-wide approach, where a variety of factors in addition to technical considerations are involved. Lastly, there is a paucity of studies focused on the tools and frameworks for supporting decision-making specifically on sewerage system restoration after earthquakes. This thesis develops a decision support framework for sewerage pipe and system restoration after earthquakes, building on the experience and learning of the organisations involved in recovering the CSS following the CES in 2010-2011. The proposed decision support framework includes three modules: 1) Physical Damage Module (PDM); 2) Functional Impact Module (FIM); 3) Pipeline Restoration Module (PRM). The PDM provides seismic fragility matrices and functions for sewer gravity and pressure pipelines for predicting earthquake-induced physical damage, categorised by pipe materials and liquefaction zones. The FIM demonstrates a set of performance indicators that are categorised in five domains: structural, hydraulic, environmental, social and economic domains. These performance indicators are used to assess loss of wastewater system service and the induced functional impacts in three different phases: emergency response, short-term recovery and long-term restoration. Based on the knowledge of the physical and functional status-quo of the sewerage systems post-earthquake captured through the PDM and FIM, the PRM estimates restoration time of sewer networks by use of restoration models developed using a Random Forest technique and graphically represented in terms of restoration curves. The development of a decision support framework for sewer recovery after earthquakes enables decision makers to assess physical damage, evaluate functional impacts relating to hydraulic, environmental, structural, economic and social contexts, and to predict restoration time of sewerage systems. Furthermore, the decision support framework can be potentially employed to underpin system maintenance and upgrade by guiding system rehabilitation and to monitor system behaviours during business-as-usual time. In conjunction with expert judgement and best practices, this framework can be moreover applied to assist asset managers in targeting the inclusion of system resilience as part of asset maintenance programmes.
The previously unknown Greendale Fault was buried beneath the Canterbury Plains and ruptured in the September 4th 2010 moment magnitude (Mw) 7.1 Darfield Earthquake. The Darfield Earthquake and subsequent Mw 6 or greater events that caused damage to Christchurch highlight the importance of unmapped faults near urban areas. This thesis examines the morphology, age and origin of the Canterbury Plains together with the paleoseismology and surface-rupture displacement distributions of the Greendale Fault. It offers new insights into the surface-rupture characteristics, paleoseismology and recurrence interval of the Greendale Fault and related structures involved in the 2010 Darfield Earthquake. To help constrain the timing of the penultimate event on the Greendale Fault the origin and age of the faulted glacial outwash deposits have been examined using sedimentological analysis of gravels and optically stimulated luminescence (OSL) dating combined with analysis of GPS and LiDAR survey data. OSL ages from this and other studies, and the analysis of surface paleochannel morphology and subsurface gravel deposits indicate distinct episodes of glacial outwash activity across the Canterbury Plains, at ~20 to 24 and ~28 to 33 kyr separated by a hiatus in sedimentation possibly indicating an interstadial period. These data suggest multiple glacial periods between ~18 and 35 kyr which may have occurred throughout the Canterbury region and wider New Zealand. A new model for the Waimakariri Fan is proposed where aggradation is mainly achieved during episodic sheet flooding with the primary river channel location remaining approximately fixed. The timing, recurrence interval and displacements of the penultimate surface-rupturing earthquake on the Greendale Fault have been constrained by trenching the scarp produced in 2010 at two locations. These excavations reveal a doubling of the magnitude of surface displacement at depths of 2-4 m. Aided by OSL ages of sand lenses in the gravel deposits, this factor-of-two increase is interpreted to indicate that in the central section of the Greendale Fault the penultimate surface-rupturing event occurred between ca. 20 and 30 kyr ago. The Greendale Fault remained undetected prior to the Darfield earthquake because the penultimate fault scarp was eroded and buried during Late Pleistocene alluvial activity. The Darfield earthquake rupture terminated against the Hororata Anticline Fault (HAF) in the west and resulted in up to 400 mm of uplift on the Hororata Anticline immediately above the HAF. Folding in 2010 is compared to Quaternary and younger deformation across the anticline recorded by a seismic reflection line, GPS-measured topographic profiles along fluvial surfaces, and river channel sinuosity and morphology. It is concluded that the HAF can rupture during earthquakes dissimilar to the 2010 event that may not be triggered by slip on the Greendale Fault. Like the Greendale Fault geomorphic analyses provide no evidence for rupture of the HAF in the last 18 kyr, with the average recurrence interval for the late Quaternary inferred to be at least ~10 kyr. Surface rupture of the Greendale Fault during the Darfield Earthquake produced one of the most accessible and best documented active fault displacement and geometry datasets in the world. Surface rupture fracture patterns and displacements along the fault were measured with high precision using real time kinematic (RTK) GPS, tape and compass, airborne light detection and ranging (LiDAR), and aerial photos. This allowed for detailed analysis of the cumulative strike-slip displacement across the fault zone, displacement gradient (ground shear strain) and the type of displacement (i.e. faulting or folding). These strain profiles confirm that the rupture zone is generally wide (~30 to ~300 metres) with >50% of displacement (often 70-80%) accommodated by ground flexure rather than discrete fault slip and ground cracking. The greatest fault-zone widths and highest proportions of folding are observed at fault stepovers.
Liquefaction of sandy soil has been observed to cause significant damage to infrastructure during major earthquakes. Historical cases of liquefaction have typically occurred in sands containing some portion of fines particles, which are defined as 75μm or smaller in diameter. The effects of fines on the undrained behaviour of sand are not however fully understood, and this study therefore attempts to quantify these effects through the undrained testing of sand mixed with non-plastic fines sourced from Christchurch, New Zealand. The experimental program carried out during this study consisted of undrained monotonic and cyclic triaxial tests performed on three different mixtures of sand and fines: the Fitzgerald Bridge mixture (FBM), and two Pinnacles Sand mixtures (PSM1 and PSM2). The fines content of each host sand was systematically varied up to a maximum of 30%, with all test specimens being reconstituted using moist tamping deposition. The undrained test results from the FBM soils were interpreted using a range of different measures of initial state. When using void ratio and relative density, the addition of fines to the FBM sand caused more contractive behaviour for both monotonic and cyclic loadings. This resulted in lower strengths at the steady state of deformation, and lower liquefaction resistances. When the intergranular void ratio was used for the interpretation, the effect of additional fines was to cause less contractive response in the sand. The state parameter and state index were also used to interpret the undrained cyclic test results – these measures suggested that additional fines caused less contractive sand behaviour, the opposite to that observed when using the void ratio. This highlighted the dependency on the parameter chosen as a basis for the response comparison when determining the effects of fines, and pointed out a need to identify a measure that normalizes such effects. Based on the FBM undrained test results and interpretations, the equivalent granular void ratio, e*, was identified from the literature as a measure of initial state that normalizes the effects of fines on the undrained behaviour of sand up to a fines content of 30%. This is done through a parameter within the e* definition termed the fines influence factor, b, which quantifies the effects of fines from a value of zero (no effect) to one (same effect as sand particles). The value of b was also determined to be different when interpreting the steady state lines (bSSL) and cyclic resistance curves (bCR) respectively for a given mixture of sand and fines. The steady state lines and cyclic resistance curves of the FBM soils and a number of other sand-fines mixtures sourced from the literature were subsequently interpreted using the equivalent granular void ratio concept, with bSSL and bCR values being back-calculated from the respective test data sets. Based on these interpretations, it was concluded that e* was conceptually a useful parameter for characterizing and quantifying the effects of fines on the undrained behaviour of sand, assuming the fines influence factor value could be derived. To allow prediction of the fines influence factor values, bSSL and bCR were correlated with material and depositional properties of the presented sand-fines mixtures. It was found that as the size of the fines particles relative to the sand particles became smaller, the values of bSSL and bCR reduced, indicating lower effect of fines. The same trend was also observed as the angularity of the sand particles increased. The depositional method was found to influence the value of bCR, due to the sensitivity of cyclic loading to initial soil fabric. This led to bSSL being used as a reference for the effect of fines, with specimens prepared by moist tamping having bCR > bSSL, and specimens prepared by slurry deposition having bCR < bSSL. Finally the correlations of the fines influence factor values with material and depositional properties were used to define the simplified estimation method – a procedure capable of predicting the approximate steady state lines and cyclic resistance curves of a sand as the non-plastic fines content is increased up to 30%. The method was critically reviewed based on the undrained test results of the PSM1 and PSM2 soils. This review suggested the method could accurately predict undrained response curves as the fines content was raised, based on the PSM1 test results. It also however identified some key issues with the method, such as the inability to accurately predict the responses of highly non-uniform soils, a lack of consideration for the entire particle size distribution of a soil, and the fact the errors in the prediction of bSSL carry through into the prediction of bCR. Lastly some areas of further investigation relating to the method were highlighted, including the need to verify the method through testing of sandy soils sourced from outside the Christchurch area, and the need to correlate the value of bCR with additional soil fabrics / depositional methods.
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.
<strong>Sea level rise is one consequence of Earth’s changing climate. Century-long tide gauge records show that global-mean sea-level rise reached 11-16 cm during the twentieth century at a mean rate of 1.2 mm/y. Today, the average rate of global-mean sea-level rise is higher at 3-4 mm/y and is expected to increase in the future. This represents a hazard to low elevation coastal zones worldwide. Yet, before global sea level projections can be used to characterise future coastal flood hazard at a local scale, the effects of tectonics (and other processes) that drive vertical land motion (VLM) must be considered. VLM is defined as the vertical velocity (uplift or subsidence) of the solid surface with respect to the centre of Earth. In this study, new VLM maps are generated over coastal strips in New Zealand, using Sentinel-1 InSAR and GNSS data.</strong>In New Zealand, measuring VLM using InSAR on naturally vegetated or agricultural land is difficult due to signal decorrelation. Along the rural Bay of Plenty coastal strip, I use a persistent-scatterer approach to generate a VLM map from both east-looking ascending and west-looking descending Sentinel-1 data between 2015-2021. Using time-series data over the same time period from a dense network of 20 GNSS sensors, I tie InSAR-derived line-of-sight velocity to the 2014 ITRF reference frame. I test two different methods for measuring VLM and compare the results against GNSS vertical velocity along the Bay of Plenty coast. Best results are achieved by first removing the interpolated horizontal GNSS velocity field from each of the InSAR datasets, before averaging the two VLM estimates. Measured VLM is between -3 and 3 mm/y, with negative values (subsidence) occurring within the low-lying Rangitāiki Plain and Ōpōtiki valley, and uplift across the elevated region west of Matatā.This thesis integrates geomorphological, geological, and historical levelling VLM records with modern satellite datasets to assess VLM across timescales ranging from 10 to 100,000 years at Matatā. Uplift rate has been variable through time, with average uplift over the last 300,000 years of 1 mm/y, 4.5 mm/y since 1720 years, 2 mm/y between 1950-1978, and 10 mm/y between 2004-2011. Previous modelling has shown that the best fit to the 2004-2011 rapid uplift rates is an inflating magmatic source at ~10 km depth beneath Matatā. To reconcile all data, I present a VLM model that consists of short-lived periods (7 years) of rapid uplift (10 mm/y), separated by longer periods (30 years) of lower background uplift (3 mm/y). The episodic nature of VLM at Matatā likely reflects short-lived periods of magmatic intrusion. Episodic VLM characterised by large rates of uplift (10 mm/y) has been seen at Taupō volcano, and other volcanic centers globally. It has been 12 years since the end of the last intrusion episode; this modelling suggest one may expect to observe increased uplift rates at Matatā in the coming decades. Densely populated urban coastal strips are most at risk from the effects of relative sea-level rise. At the same time, anthropogenic activities associated with urbanization, such as groundwater withdrawal, and land reclamation can lead to local land subsidence (LLS), further exacerbating the risk to urban infrastructure. LLS refers to subsidence relative to nearby land area assumed to be stable. In this thesis, I create the first high-resolution (10 m) maps of LLS at six urban coastal strips in New Zealand, with a combined length of 285 km, using Sentinel-1 InSAR data between 2018-2021. This analysis reveals 89% of urban coastal strips are subsiding at rates of -0.5 mm/y or greater, and 11% is subsiding at higher rates of -3.0 mm/y or greater. On average, subsidence is -0.6 to -2.9 mm/y higher at the coastal strip, compared to inland areas occupied by GNSS stations. This analysis also documents highly-localised hotspots of LLS, with subsidence rates of up to -15 mm/y. In Christchurch, rapid and localised subsidence (-8 mm/y) is observed within coastal suburbs New Brighton and Southshore. In most cities, the highest subsidence rates occur on land reclaimed in the early-late twentieth century, and in areas built on Holocene sediment. Time-series analysis of LLS at sites of reclaimed land shows both linear and non-linear rates of deformation over time periods of up to 6-8 years. This thesis highlights the variable exposure to relative sea-level rise of New Zealand coastal strips, and demonstrates that in many cases current rates of VLM should be expected to continue for the next few decades.
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.
