A major hazard accompanying earthquake shaking in areas of steep topography is the detachment of rocks from bedrock outcrops that subsequently slide, roll, or bounce downslope (i.e. rockfalls). The 2010-2011 Canterbury earthquake sequence caused recurrent and severe rockfall in parts of southern Christchurch. Coseismic rockfall caused five fatalities and significant infrastructural damage during the 2011 Mw 6.2 Christchurch earthquake. Here we examine a rockfall site in southern Christchurch in detail using geomorphic mapping, lidar analysis, geochronology (cosmogenic 3He dating, radiocarbon dating, optically stimulated luminescence (OSL) from quartz, infrared stimulated luminescence from K-feldspar), numerical modeling of rockfall boulder trajectories, and ground motion prediction equations (GMPEs). Rocks fell from the source cliff only in earthquakes with interpolated peak ground velocities exceeding ~10 cm/s; hundreds of smaller earthquakes did not produce rockfall. On the basis of empirical observations, GMPEs and age chronologies we attribute paleo-rockfalls to strong shaking in prehistoric earthquakes. We conclude that earthquake shaking of comparable intensity to the strongest contemporary earthquakes in Christchurch last occurred at this site approximately 5000 to 7000 years ago, and that in some settings, rockfall deposits provide useful proxies for past strong ground motions.
he strong motion station at Heathcote Valley School (HVSC) recorded unusually high peak ground accelerations (2.21g vertical and 1.41g horizontal) during the February 2011 Christchurch earthquake. Ground motions recorded at HVSC in numerous other events also exhibited consistently higher intensities compared with nearby strong motion stations. We investigated the underlying causes of such high intensity ground motions at HVSC by means of 2D dynamic finite element analyses, using recorded ground motions during the 2010-2011 Canterbury earthquake sequence. The model takes advantage of a LiDAR-based digital elevation model (DEM) to account for the surface topography, while the geometry and dynamic properties of the surficial soils are characterized by seismic cone penetration tests (sCPT) and Multi-Channel Analyses of Surface Waves (MASW). Comparisons of simulated and recorded ground motions suggests that our model performs well for distant events, while for near-field events, ground motions recorded at the adopted reference station at Lyttelton Port are not reasonable input motions for the simulation. The simulations suggest that Rayleigh waves generated at the inclined interface of the surficial colluvium and underlying volcanic rock strongly affect the ground motions recorded at HVSC, in particular, being the dominant contributor to the recorded vertical motions.
The empirical liquefaction triggering chart of Idriss and Boulanger (2008) is compared to direct measurements of the cyclic resistance of Christchurch silty sands via undisturbed and reconstituted lab specimens. Comparisons suggest that overall there is a reasonable agreement between the empirical triggering curve and the interpreted test data. However, the influence of fines on cyclic resistance appears to be over-predicted by the empirical method, particularly for non-plastic silty sands that are commonly encountered in flood over-bank deposits in Christchurch and nearby settlements
The Screw Driving Sounding (SDS) method developed in Japan is a relatively new insitu testing technique to characterise soft shallow sites, typically those required for residential house construction. An SDS machine drills a rod into the ground in several loading steps while the rod is continuously rotated. Several parameters, such as torque, load and speed of penetration, are recorded at every rotation of the rod. The SDS method has been introduced in New Zealand, and the results of its application for characterising local sites are discussed in this study. A total of 164 SDS tests were conducted in Christchurch, Wellington and Auckland to validate/adjust the methodologies originally developed based on the Japanese practice. Most of the tests were conducted at sites where cone penetration tests (CPT), standard penetration tests (SPT) and borehole logs were available; the comparison of SDS results with existing information showed that the SDS method has great potential as an in-situ testing method for classifying the soils. By compiling the SDS data from 3 different cities and comparing them with the borehole logs, a soil classification chart was generated for identifying the soil type based on SDS parameters. Also, a correlation between fines content and SDS parameters was developed and a procedure for estimating angle of internal friction of sand using SDS parameters was investigated. Furthermore, a correlation was made between the tip resistance of the CPT and the SDS data for different percentages of fines content. The relationship between the SPT N value and a SDS parameter was also proposed. This thesis also presents a methodology for identifying the liquefiable layers of soil using SDS data. SDS tests were performed in both liquefied and non-liquefied areas in Christchurch to find a representative parameter and relationship for predicting the liquefaction potential of soil. Plots were drawn of the cyclic shear stress ratios (CSR) induced by the earthquakes and the corresponding energy of penetration during SDS tests. By identifying liquefied or unliquefied layers using three different popular CPT-based methods, boundary lines corresponding to the various probabilities of liquefaction happening were developed for different ranges of fines contents using logistic regression analysis, these could then be used for estimating the liquefaction potential of soil directly from the SDS data. Finally, the drilling process involved in screw driving sounding was simulated using Abaqus software. Analysis results proved that the model successfully captured the drilling process of the SDS machine in sand. In addition, a chart to predict peak friction angles of sandy sites based on measured SDS parameters for various vertical effective stresses was formulated. As a simple, fast and economical test, the SDS method can be a reliable alternative insitu test for soil and site characterisation, especially for residential house construction.
A photograph of a fence post that lifted out of the ground during the earthquake, pulling a clump of soil with it.
Scientists are calling for more ground testing to be carried out before reconstruction starts in earthquake devastated Christchurch.
Soil-structure interaction (SSI) has been widely studied during the last decades. The influence of the properties of the ground motion, the structure and the soil have been addressed. However, most of the studies in this field consider a stand-alone structure. This assumption is rarely justifiable in dense urban areas where structures are built close to one another. The dynamic interaction between adjacent structures has been studied since the early 1970s, mainly using numerical and analytical models. Even though the early works in this field have significantly contributed to understanding this problem, they commonly consider important simplifications such as assuming a linear behaviour of the structure and the soil. Some experimental works addressing adjacent structures have recently been conducted using geotechnical centrifuges and 1g shake tables. However, further research is needed to enhance the understanding of this complex phenomenon. A particular case of SSI is that of structures founded in fine loose saturated sandy soil. An iconic example was the devastating effects of liquefaction in Christchurch, New Zealand, during the Canterbury earthquake in 2011. In the case of adjacent structures on liquefiable soil, the experimental evidence is even scarcer. The present work addresses the dynamic interaction between adjacent structures by performing multiple experimental studies. The work starts with two-adjacent structures on a small soil container to expose the basics of the problem. Later, results from tests considering a more significant number of structures on a big laminar box filled with sand are presented. Finally, the response of adjacent structures on saturated sandy soil is addressed using a geotechnical centrifuge and a large 1g shake table. This research shows that the acceleration, lateral displacement, foundation rocking, damping ratio, and fundamental frequency of the structure of focus are considerably affected by the presence of neighbouring buildings. In general, adjacent buildings reduced the dynamic response of the structure of focus on dry sand. However, the acceleration was amplified when the structures had a similar fundamental frequency. In the case of structures on saturated sand, the presence of adjacent structures reduced the liquefaction potential. Neighbouring structures on saturated sand also presented larger rotation of the footing and lateral displacement of the top mass than that of the stand-alone case.
A box containing drilling cores from soil sampling. The photographer comments, "These are the samples from drilling near my home. As you can see they are not so much samples as sand piles. The drill in a nearby street went down 20m and it was sand all the way. This is the box of samples from the ground level to 4.6m deep".
Photograph captioned by Fairfax, "Demolition underway on the Manchester Courts Building. Street sign dwarfed by the large mound of soil alongside the building".
A photograph of contractor Tony Fisher of Fisher Agricultural Ltd. inspecting a dug out soil pit on a farm near River Road in Lincoln.
The Catholic Cathedral is classified as a category 1 listed heritage building constructed largely of unreinforced stone masonry, and was significantly damaged in the recent Canterbury earthquakes of 2010 and 2011. In the 2010 event the building presented slight to moderta damage, meanwhile in the 2011 one experienced ground shaking in excess of its capacity leading to block failures and partial collapse of parts of the building, which left the building standing but still posing a significant hazard. In this paper we discuss the approach to develop the earthquake analysis of the building by 3D numerical simulations, and the results are compared/calibrated with the observed damage of the 2010 earthquake. Very accurate records were obtained during both earthquakes due to a record station located least than 80 m of distance from the building and used in the simulations. Moreover it is included in the model the soil structure interaction because it was observed that the ground and foundation played an important role on the seismic behavior of the structure. A very good agreement was found between the real observed damage and the nonlinear dynamic simulations described trough inelastic deformation (cracking) and building´s performance.
A photograph of Peter Almond (wearing high visibility jacket) and Derrick Moot inspecting a dug-out soil pit on a farm near River Road in Lincoln.
This report summarizes the development of a region-wide surficial soil shear wave velocity (Vs ) model based on the unique combination of a large high-spatial-density database of cone penetration test (CPT) logs in the greater Christchurch urban area (> 15, 000 logs as of 1 February 2014) and the Christchurch-specific empirical correlation between soil Vs and CPT data developed by McGann et al. [1, 2]. This model has applications for site characterization efforts via maps of time-averaged Vs over specific depths (e.g. Vs30, Vs10), and for numerical modeling efforts via the identification of typical Vs profiles for different regions and soil behaviour types within Christchurch. In addition, the Vs model can be used to constrain the near-surface velocities for the 3D seismic velocity model of the Canterbury basin [3] currently being developed for the purpose of broadband ground motion simulation. The general development of these region-wide near-surface Vs models includes the following general phases, with each discussed in separate chapters of this report. • An evaluation of the available CPT dataset for suitability, and the definition of other datasets and assumptions necessary to characterize the surficial sediments of the region to 30 m depth. • The development of time-averaged shear wave velocity (Vsz) surfaces for the Christchurch area from the adopted CPT dataset (and supplementary data/assumptions) using spatial interpolation. The Vsz surfaces are used to explore the characteristics of the near-surface soils in the regions and are shown to correspond well with known features of the local geology, the historical ecosystems of the area, and observations made following the 2010- 2011 Canterbury earthquakes. • A detailed analysis of the Vs profiles in eight subregions of Christchurch is performed to assess the variablity in the soil profiles for regions with similar Vsz values and to assess Vsz as a predictive metric for local site response. It is shown that the distrubution of soil shear wave velocity in the Christchurch regions is highly variable both spatially (horizontally) and with depth (vertically) due to the varied geological histories for different parts of the area, and the highly stratified nature of the nearsurface deposits. This variability is not considered to be greatly significant in terms of current simplified site classification systems; based on computed Vs30 values, all considered regions can be categorized as NEHRP sites class D (180 < Vs < 360 m/s) or E (Vs < 180 m/s), however, detailed analysis of the shear wave velocity profiles in different subregions of Christchurch show that the expected surficial site response can vary quite a bit across the region despite the relative similarity in Vs30
Recent experiences from the Darfield and Canterbury, New Zealand earthquakes have shown that the soft soil condition of saturated liquefiable sand has a profound effect on seismic response of buildings, bridges and other lifeline infrastructure. For detailed evaluation of seismic response three dimensional integrated analysis comprising structure, foundation and soil is required; such an integrated analysis is referred to as Soil Foundation Structure Interaction (SFSI) in literatures. SFSI is a three-dimensional problem because of three primary reasons: first, foundation systems are three-dimensional in form and geometry; second, ground motions are three-dimensional, producing complex multiaxial stresses in soils, foundations and structure; and third, soils in particular are sensitive to complex stress because of heterogeneity of soils leading to a highly anisotropic constitutive behaviour. In literatures the majority of seismic response analyses are limited to plane strain configuration because of lack of adequate constitutive models both for soils and structures, and computational limitation. Such two-dimensional analyses do not represent a complete view of the problem for the three reasons noted above. In this context, the present research aims to develop a three-dimensional mathematical formulation of an existing plane-strain elasto-plastic constitutive model of sand developed by Cubrinovski and Ishihara (1998b). This model has been specially formulated to simulate liquefaction behaviour of sand under ground motion induced earthquake loading, and has been well-validated and widely implemented in verifcation of shake table and centrifuge tests, as well as conventional ground response analysis and evaluation of case histories. The approach adopted herein is based entirely on the mathematical theory of plasticity and utilises some unique features of the bounding surface plasticity formalised by Dafalias (1986). The principal constitutive parameters, equations, assumptions and empiricism of the existing plane-strain model are adopted in their exact form in the three-dimensional version. Therefore, the original two-dimensional model can be considered as a true subset of the three-dimensional form; the original model can be retrieved when the tensorial quantities of the three dimensional version are reduced to that of the plane-strain configuration. Anisotropic Drucker-Prager type failure surface has been adopted for the three-dimensional version to accommodate triaxial stress path. Accordingly, a new mixed hardening rule based on Mroz’s approach of homogeneous surfaces (Mroz, 1967) has been introduced for the virgin loading surface. The three-dimensional version is validated against experimental data for cyclic torsional and triaxial stress paths.
A number of field testing techniques, such as standard penetration test (SPT), cone penetration test (CPT), and Swedish weight sounding (SWS), are popularly used for in-situ characterisation. The screw driving sounding (SDS) method, which has been recently developed in Japan, is an improved version of the SWS technique and measures more parameters, including the required torque, load, speed of penetration and rod friction; these provide more robust way of characterising soil stratigraphy. It is a cost-efficient technique which uses a machine-driven and portable device, making it ideal for testing in small-scale and confined areas. Moreover, with a testing depth of up to 10-15m, it is suitable for liquefaction assessment. Thus, the SDS method has great potential as an in-situ testing method for geotechnical site characterisation, especially for residential house construction. In this paper, the results of SDS tests performed at a variety of sites in New Zealand are presented. The soil database was employed to develop a soil classification chart based on SDS-derived parameters. Moreover, using the data obtained following the 2010-2011 Christchurch Earthquake Se-quence, a methodology was established for liquefaction potential evaluation using SDS data. http://www.isc5.com.au/wp-content/uploads/2016/09/1345-2-ORENSE.pdf
A photograph of a crack in the pavement of a residential street in Christchurch. The crack was caused by liquefaction erupting out of the soil underneath.
A photograph of University of Canterbury Geology staff and students using a dynamic cone penetrometer to measure soil strength of a paddock on the Greendale fault line.
This paper presents insights from recent advanced laboratory testing of undisturbed and reconstituted specimens of Christchurch silty-sands. The purpose of the testing was to establish the cyclic strength of silty-sands from sites in the Central Business District (CBD), where liquefaction was observed in 4 September 2010, 22 February 2011, and 13 June 2011. Similar overall strengths were obtained from undisturbed and reconstituted tests prepared at similar densities, albeit with higher variability for the reconstituted specimens. Reconstituted specimens exhibited distinctly different response in terms of lower compressibility during initial loading cycles, and exhibited a more brittle response when large strains were mobilised, particularly for samples with high fines content. Given the lower variability in natural sample response and the possibility of age-related strength to be significant for sites not subjected to earthquakes, high quality undisturbed samples are recommended over the use of reconstituted specimens to establish the cyclic strength of natural sands.
The influence of nonlinear soil-foundation-structure interaction (SFSI) on the performance of multi-storey buildings during earthquake events has become increasingly important in earthquake resistant design. For buildings on shallow foundations, SFSI refers to nonlinear geometric effects associated with uplift of the foundation from the supporting soil as well as nonlinear soil deformation effects. These effects can potentially be beneficial for structural performance, reducing forces transmitted from ground shaking to the structure. However, there is also the potential consequence of residual settlement and rotation of the foundation. This Thesis investigates the influence of SFSI in the performance of multi-storey buildings on shallow foundations through earthquake observations, experimental testing, and development of spring-bed numerical models that can be incorporated into integrated earthquake resistant design procedures. Observations were made following the 22 February 2011 Christchurch Earthquake in New Zealand of a number of multi-storey buildings on shallow foundations that performed satisfactorily. This was predominantly the case in areas where shallow foundations, typically large raft foundations, were founded on competent gravel and where there was no significant manifestation of liquefaction at the ground surface. The properties of these buildings and the soils they are founded on directed experimental work that was conducted to investigate the mechanisms by which SFSI may have influenced the behaviour of these types of structure-foundation systems. Centrifuge experiments were undertaken at the University of Dundee, Scotland using a range of structure-foundation models and a layer of dense cohesionless soil to simulate the situation in Christchurch where multi-storey buildings on shallow foundations performed well. Three equivalent single degree of freedom (SDOF) models representing 3, 5, and 7 storey buildings with identical large raft foundations were subjected to a range of dynamic Ricker wavelet excitations and Christchurch Earthquake records to investigate the influence of SFSI on the response of the equivalent buildings. The experimental results show that nonlinear SFSI has a significant influence on structural response and overall foundation deformations, even though the large raft foundations on competent soil meant that there was a significant reserve of bearing capacity available and nonlinear deformations may have been considered to have had minimal effect. Uplift of the foundation from the supporting soil was observed across a wide range of input motion amplitudes and was particularly significant as the amplitude of motion increased. Permanent soil deformation represented by foundation settlement and residual rotation was also observed but mainly for the larger input motions. However, the absolute extent of uplift and permanent soil deformation was very small compared to the size of the foundation meaning the serviceability of the building would still likely be maintained during large earthquake events. Even so, the small extent of SFSI resulted in attenuation of the response of the structure as the equivalent period of vibration was lengthened and the equivalent damping in the system increased. The experimental work undertaken was used to validate and enhance numerical modelling techniques that are simple yet sophisticated and promote interaction between geotechnical and structural specialists involved in the design of multi-storey buildings. Spring-bed modelling techniques were utilised as they provide a balance between ease of use, and thus ease of interaction with structural specialists who have these techniques readily available in practice, and theoretically rigorous solutions. Fixed base and elastic spring-bed models showed they were unable to capture the behaviour of the structure-foundation models tested in the centrifuge experiments. SFSI spring-bed models were able to more accurately capture the behaviour but recommendations were proposed for the parameters used to define the springs so that the numerical models closely matched experimental results. From the spring-bed modelling and results of centrifuge experiments, an equivalent linear design procedure was proposed along with a procedure and recommendations for the implementation of nonlinear SFSI spring-bed models in practice. The combination of earthquake observations, experimental testing, and simplified numerical analysis has shown how SFSI is influential in the earthquake performance of multi-storey buildings on shallow foundations and should be incorporated into earthquake resistant design of these structures.
A photograph of a dug-out soil pit on a farm near River Road in Lincoln, showing saturated sand between the topsoil and the subsoil where the grass roots end.
A photograph of a dug-out soil pit on a farm near River Road in Lincoln, showing saturated sand between the topsoil and the subsoil where the grass roots end.
A photograph of a dug-out soil pit on a farm near River Road in Lincoln, showing saturated sand between the topsoil and the subsoil where the grass roots end.
A photograph of a dug-out soil pit on a farm near River Road in Lincoln, showing saturated sand between the topsoil and the subsoil where the grass roots end.
A photograph of a dug-out soil pit on a farm near River Road in Lincoln, showing saturated sand between the topsoil and the subsoil where the grass roots end.
A photograph of a dug-out soil pit on a farm near River Road in Lincoln, showing saturated sand between the topsoil and the subsoil where the grass roots end.
A photograph of a dug-out soil pit on a farm near River Road in Lincoln, showing saturated sand between the topsoil and the subsoil where the grass roots end.
A photograph of a dug-out soil pit on a farm near River Road in Lincoln, showing saturated sand between the topsoil and the subsoil where the grass roots end.
None
In 2010 and 2011 Christchurch, New Zealand experienced a series of earthquakes that caused extensive damage across the city, but primarily to the Central Business District (CBD) and eastern suburbs. A major feature of the observed damage was extensive and severe soil liquefaction and associated ground damage, affecting buildings and infrastructure. The behaviour of soil during earthquake loading is a complex phenomena that can be most comprehensively analysed through advanced numerical simulations to aid engineers in the design of important buildings and critical facilities. These numerical simulations are highly dependent on the capabilities of the constitutive soil model to replicate the salient features of sand behaviour during cyclic loading, including liquefaction and cyclic mobility, such as the Stress-Density model. For robust analyses advanced soil models require extensive testing to derive engineering parameters under varying loading conditions for calibration. Prior to this research project little testing on Christchurch sands had been completed, and none from natural samples containing important features such as fabric and structure of the sand that may be influenced by the unique stress-history of the deposit. This research programme is focussed on the characterisation of Christchurch sands, as typically found in the CBD, to facilitate advanced soil modelling in both res earch and engineering practice - to simulate earthquake loading on proposed foundation design solutions including expensive ground improvement treatments. This has involved the use of a new Gel Push (GP) sampler to obtain undisturbed samples from below the ground-water table. Due to the variable nature of fluvial deposition, samples with a wide range of soil gradations, and accordingly soil index properties, were obtained from the sampling sites. The quality of the samples is comprehensively examined using available data from the ground investigation and laboratory testing. A meta-quality assessment was considered whereby a each method of evaluation contributed to the final quality index assigned to the specimen. The sampling sites were characterised with available geotechnical field-based test data, primarily the Cone Penetrometer Test (CPT), supported by borehole sampling and shear-wave velocity testing. This characterisation provides a geo- logical context to the sampling sites and samples obtained for element testing. It also facilitated the evaluation of sample quality. The sampling sites were evaluated for liquefaction hazard using the industry standard empirical procedures, and showed good correlation to observations made following the 22 February 2011 earthquake. However, the empirical method over-predicted liquefaction occurrence during the preceding 4 September 2010 event, and under-predicted for the subsequent 13 June 2011 event. The reasons for these discrepancies are discussed. The response of the GP samples to monotonic and cyclic loading was measured in the laboratory through triaxial testing at the University of Canterbury geomechanics laboratory. The undisturbed samples were compared to reconstituted specimens formed in the lab in an attempt to quantify the effect of fabric and structure in the Christchurch sands. Further testing of moist tamped re- constituted specimens (MT) was conducted to define important state parameters and state-dependent properties including the Critical State Line (CSL), and the stress-strain curve for varying state index. To account for the wide-ranging soil gradations, selected representative specimens were used to define four distinct CSL. The input parameters for the Stress-Density Model (S-D) were derived from a suite of tests performed on each representative soil, and with reference to available GP sample data. The results of testing were scrutinised by comparing the data against expected trends. The influence of fabric and structure of the GP samples was observed to result in similar cyclic strength curves at 5 % Double Amplitude (DA) strain criteria, however on close inspection of the test data, clear differences emerged. The natural samples exhibited higher compressibility during initial loading cycles, but thereafter typically exhibited steady growth of plastic strain and excess pore water pressure towards and beyond the strain criteria and initial liquefaction, and no flow was observed. By contrast the reconstituted specimens exhibited a stiffer response during initial loading cycles, but exponential growth in strains and associated excess pore water pressure beyond phase-transformation, and particularly after initial liquefaction where large strains were mobilised in subsequent cycles. These behavioural differences were not well characterised by the cyclic strength curve at 5 % DA strain level, which showed a similar strength for both GP samples and MT specimens. A preliminary calibration of the S-D model for a range of soil gradations is derived from the suite of laboratory test data. Issues encountered include the influence of natural structure on the peak-strength–state index relationship, resulting in much higher peak strengths than typically observed for sands in the literature. For the S-D model this resulted in excessive stiffness to be modelled during cyclic mobility, when the state index becomes large momentarily, causing strain development to halt. This behaviour prevented modelling the observed re- sponse of silty sands to large strains, synonymous with “liquefaction”. Efforts to reduce this effect within the current formulation are proposed as well as future research to address this issue.
The sample of water referred to in the present note was collected by the writer on the 21st January, 1889, in the Otira Gorge, from a spring which is stated to have been first discovered shortly after the earthquake of the 1st September, 1888. From the results obtained this water might be termed siliceous and sulphurous. It is essentially different from the water from the Hanmer Springs, and pertains more to the character of the waters of the Rotorua district. It differs, however, from these waters in having only a portion of its carbonic anhydride replaced by silica, and in containing less dissolved matter.