The 22 February 2011, Mw6.2 Christchurch earthquake is the most costly earthquake to affect New Zealand, causing an estimated 181 fatalities and severely damaging thousands of residential and commercial buildings. This paper presents a summary of some of the observations made by the NSF-sponsored GEER Team regarding the geotechnical/geologic aspects of this earthquake. The Team focused on documenting the occurrence and severity of liquefaction and lateral spreading, performance of building and bridge foundations, buried pipelines and levees, and significant rockfalls and landslides. Liquefaction was pervasive and caused extensive damage to residential properties, water and wastewater networks, high-rise buildings, and bridges. Entire neighborhoods subsided, resulting in flooding that caused further damage. Additionally, liquefaction and lateral spreading resulted in damage to bridges and to stretches of levees along the Waimakariri and Kaiapoi Rivers. Rockfalls and landslides in the Port Hills damaged several homes and caused several fatalities.
Earthquake-triggered soil liquefaction caused extensive damage and heavy economic losses in Christchurch during the 2010-2011 Canterbury earthquakes. The most severe manifestations of liquefaction were associated with the presence of natural deposits of clean sands and silty sands of fluvial origin. However, liquefaction resistance of fines-containing sands is commonly inferred from empirical relationships based on clean sands (i.e. sands with less than 5% fines). Hence, existing evaluation methods have poor accuracy when applied to silty sands. The liquefaction behaviour of Christchurch fines-containing (silty) sands is investigated through a series of Direct Simple Shear (DSS) tests. This type of test better resembles earthquake loading conditions in soil deposits compared to cyclic triaxial tests. Soil specimens are reconstituted in the laboratory with the water sedimentation technique. This preparation method yields soil fabrics similar to those encountered in fluvial soil deposits, which are common in the Christchurch area. Test results provide preliminary indications on how void ratio, relative density, preparation method and fines content influence the cyclic liquefaction behaviour of sand-silt mixtures depending on the properties of host sand and silt.
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.
Coastal and river environments are exposed to a number of natural hazards that have the potential to negatively affect both human and natural environments. The purpose of this research is to explain that significant vulnerabilities to seismic hazards exist within coastal and river environments and that coasts and rivers, past and present, have played as significant a role as seismic, engineering or socio-economic factors in determining the impacts and recovery patterns of a city following a seismic hazard event. An interdisciplinary approach was used to investigate the vulnerability of coastal and river areas in the city of Christchurch, New Zealand, following the Canterbury Earthquake Sequence, which began on the 4th of September 2010. This information was used to identify the characteristics of coasts and rivers that make them more susceptible to earthquake induced hazards including liquefaction, lateral spreading, flooding, landslides and rock falls. The findings of this research are applicable to similar coastal and river environments elsewhere in the world where seismic hazards are also of significant concern. An interdisciplinary approach was used to document and analyse the coastal and river related effects of the Canterbury earthquake sequence on Christchurch city in order to derive transferable lessons that can be used to design less vulnerable urban communities and help to predict seismic vulnerabilities in other New Zealand and international urban coastal and river environments for the future. Methods used to document past and present features and earthquake impacts on coasts and rivers in Christchurch included using maps derived from Geographical Information Systems (GIS), photographs, analysis of interviews from coastal, river and engineering experts, and analysis of secondary data on seismicity, liquefaction potential, geology, and planning statutes. The Canterbury earthquake sequence had a significant effect on Christchurch, particularly around rivers and the coast. This was due to the susceptibility of rivers to lateral spreading and the susceptibility of the eastern Christchurch and estuarine environments to liquefaction. The collapse of river banks and the extensive cracking, tilting and subsidence that accompanied liquefaction, lateral spreading and rock falls caused damage to homes, roads, bridges and lifelines. This consequently blocked transportation routes, interrupted electricity and water lines, and damaged structures built in their path. This study found that there are a number of physical features of coastal and river environments from the past and the present that have induced vulnerabilities to earthquake hazards. The types of sediments found beneath eastern Christchurch are unconsolidated fine sands, silts, peats and gravels. Together with the high water tables located beneath the city, these deposits made the area particularly susceptible to liquefaction and liquefaction-induced lateral spreading, when an earthquake of sufficient size shook the ground. It was both past and present coastal and river processes that deposited the types of sediments that are easily liquefied during an earthquake. Eastern Christchurch was once a coastal and marine environment 6000 years ago when the shoreline reached about 6 km inland of its present day location, which deposited fine sand and silts over this area. The region was also exposed to large braided rivers and smaller spring fed rivers, both of which have laid down further fine sediments over the following thousands of years. A significant finding of this study is the recognition that the Canterbury earthquake sequence has exacerbated existing coastal and river hazards and that assessments and monitoring of these changes will be an important component of Christchurch’s future resilience to natural hazards. In addition, patterns of recovery following the Canterbury earthquakes are highlighted to show that coasts and rivers are again vulnerable to earthquakes through their ability to recovery. This city’s capacity to incorporate resilience into the recovery efforts is also highlighted in this study. Coastal and river areas have underlying physical characteristics that make them increasingly vulnerable to the effects of earthquake hazards, which have not typically been perceived as a ‘coastal’ or ‘river’ hazard. These findings enhance scientific and management understanding of the effects that earthquakes can have on coastal and river environments, an area of research that has had modest consideration to date. This understanding is important from a coastal and river hazard management perspective as concerns for increased human development around coastlines and river margins, with a high seismic risk, continue to grow.
This study uses 44 high quality liquefaction case histories taken from 22 locations affected by the 2010-2011 Canterbury earthquake sequence to evaluate four commonly used CPT-VS correlations (i.e., Robertson, 2009; Hegazy and Mayne, 2006; Andrus et al., 2007; McGann et al., 2015b). Co-located CPT soundings and VS profiles, developed from surface wave testing, were obtained at 22 locations and case histories were developed for the Mw 7.1, 4 September 2010 Darfield and Mw 6.2, 22 February 2011 Christchurch earthquakes. The CPT soundings are used to generate VS profiles using each of four CPT-VS correlations. These correlated VS profiles are used to estimate the factor of safety against liquefaction using the Kayen et al. (2013) VS-based simplified liquefaction evaluation procedure. An error index is used to quantify the predictive capabilities of these correlations in relation to the observations of liquefaction (or the lack thereof). Additionally, the error indices from the CPT-correlated VS profiles are compared to those obtained using: (1) the Kayen et al. (2013) procedure with surface wave-derived VS profiles, and (2) the Idriss and Boulanger (2008) CPT-based liquefaction evaluation procedure. Based on the error indices, the evaluation procedures based on direct measurements of either CPT or VS provided more accurate liquefaction triggering estimates than those obtained from any of the CPT-VS correlations. However, the performance of the CPT-VS correlations varied, with the Robertson (2009) and Hegazy and Mayne (2006) correlations performing relatively poorly for the Christchurch soils and the Andrus et al. (2007) and McGann et al. (2015b) correlations performing better. The McGann et al. (2015b) correlation had the lowest error indices of the CPT-VS correlations tested, however, none of the CPT-VS correlations provided accurate enough VS predictions to be used for the evaluation of liquefaction triggering using the VS-based liquefaction evaluation procedures.
Well-validated liquefaction constitutive models are increasingly important as non-linear time history analyses become relatively more common in industry for key projects. Previous validation efforts of PM4Sand, a plasticity model specifically for liquefaction, have generally focused on centrifuge tests; however, pore pressure transducers installed at several free-field sites during the Canterbury Earthquake Sequence (CES) in Christchurch, New Zealand provide a relatively unique dataset to validate against. This study presents effective stress site response analyses performed in the finite difference software FLAC to examine the capability of PM4Sand to capture the generation of excess pore pressures during earthquakes. The characterization of the subsurface is primarily based on extensive cone penetration tests (CPT) carried out in Christchurch. Correlations based on penetration resistances are used to estimate soil parameters, such as relative density and shear wave velocity, which affect liquefaction behaviour. The resulting free-field FLAC model is used to estimate time histories of excess pore pressure, which are compared with records during several earthquakes in the CES to assess the suitability of PM4Sand.
Semi-empirical models based on in-situ geotechnical tests have become the standard of practice for predicting soil liquefaction. Since the inception of the “simplified” cyclic-stress model in 1971, variants based on various in-situ tests have been developed, including the Cone Penetration Test (CPT). More recently, prediction models based soley on remotely-sensed data were developed. Similar to systems that provide automated content on earthquake impacts, these “geospatial” models aim to predict liquefaction for rapid response and loss estimation using readily-available data. This data includes (i) common ground-motion intensity measures (e.g., PGA), which can either be provided in near-real-time following an earthquake, or predicted for a future event; and (ii) geospatial parameters derived from digital elevation models, which are used to infer characteristics of the subsurface relevent to liquefaction. However, the predictive capabilities of geospatial and geotechnical models have not been directly compared, which could elucidate techniques for improving the geospatial models, and which would provide a baseline for measuring improvements. Accordingly, this study assesses the realtive efficacy of liquefaction models based on geospatial vs. CPT data using 9,908 case-studies from the 2010-2016 Canterbury earthquakes. While the top-performing models are CPT-based, the geospatial models perform relatively well given their simplicity and low cost. Although further research is needed (e.g., to improve upon the performance of current models), the findings of this study suggest that geospatial models have the potential to provide valuable first-order predictions of liquefaction occurence and consequence. Towards this end, performance assessments of geospatial vs. geotechnical models are ongoing for more than 20 additional global earthquakes.
The Canterbury Earthquake Sequence 2010-2011 (CES) induced widespread liquefaction in many parts of Christchurch city. Liquefaction was more commonly observed in the eastern suburbs and along the Avon River where the soils were characterised by thick sandy deposits with a shallow water table. On the other hand, suburbs to the north, west and south of the CBD (e.g. Riccarton, Papanui) exhibited less severe to no liquefaction. These soils were more commonly characterised by inter-layered liquefiable and non-liquefiable deposits. As part of a related large-scale study of the performance of Christchurch soils during the CES, detailed borehole data including CPT, Vs and Vp have been collected for 55 sites in Christchurch. For this subset of Christchurch sites, predictions of liquefaction triggering using the simplified method (Boulanger & Idriss, 2014) indicated that liquefaction was over-predicted for 94% of sites that did not manifest liquefaction during the CES, and under-predicted for 50% of sites that did manifest liquefaction. The focus of this study was to investigate these discrepancies between prediction and observation. To assess if these discrepancies were due to soil-layer interaction and to determine the effect that soil stratification has on the develop-ment of liquefaction and the system response of soil deposits.
None
The paper presents preliminary findings from comprehensive research studies on the liquefaction-induced damage to buildings and infrastructure in Christchurch during the 2010-2011 Canterbury earthquakes. It identifies key factors and mechanisms of damage to road bridges, shallow foundations of CBD buildings and buried pipelines, and highlights the implications of the findings for the seismic analysis and design of these structures.
Earthquake-triggered soil liquefaction caused extensive damage and heavy economic losses in Christchurch during the 2010-2011 Canterbury earthquakes. The most severe manifestations of liquefaction were associated with the presence of natural deposits of clean sands and silty sands of fluvial origin. However, liquefaction resistance of fines-containing sands is commonly inferred from empirical relationships based on clean sands (i.e. sands with less than 5% fines). Hence, existing evaluation methods have poor accuracy when applied to silty sands. Also, existing methods do not quantify appropriately the influence on liquefaction resistance of soil fabric and structure, which are unique to a specific depositional environment. This study looks at the influence of fines content, soil fabric (i.e. arrangement of soil particles) and structure (e.g. layering, segregation) on the undrained cyclic behaviour and liquefaction resistance of fines-containing sandy soils from Christchurch using Direct Simple Shear (DSS) tests on soil specimens reconstituted in the laboratory with the water sedimentation technique. The poster describes experimental procedures and presents early test results on two sands retrieved at two different sites in Christchurch.
Liquefaction-induced lateral spreading in large seismic events often results in pervasive and costly damage to engineering structures and lifelines, making it a critical component of engineering design. However, the complex nature of this phenomenon leads to designing for such a hazard extremely challenging and there is a clear for an improved understanding and predicting liquefaction-induced lateral spreading. The 2010-2011 Canterbury (New Zealand) Earthquakes triggered severe liquefaction-induced lateral spreading along the streams and rivers of the Christchurch region, causing extensive damage to roads, bridges, lifelines, and structures in the vicinity. The unfortunate devastation induced from lateral spreading in these events also rendered the rare opportunity to gain an improved understanding of lateral spreading displacements specific to the Christchurch region. As part of this thesis, the method of ground surveying was employed following the 4 September 2010 Darfield (Mw 7.1) and 22 February 2011 Christchurch (Mw 6.2) earthquakes at 126 locations (19 repeated) throughout Christchurch and surrounding suburbs. The method involved measurements and then summation of crack widths along a specific alignment (transect) running approximately perpendicular to the waterway to indicate typically a maximum lateral displacement at the bank and reduction of the magnitude of displacements with distance from the river. Rigorous data processing and comparisons with alternative measurements of lateral spreading were performed to verify results from field observations and validate the method of ground surveying employed, as well as highlight the complex nature of lateral spreading displacements. The welldocumented field data was scrutinized to gain an understanding of typical magnitudes and distribution patterns (distribution of displacement with distance) of lateral spreading observed in the Christchurch area. Maximum displacements ranging from less than 10 cm to over 3.5 m were encountered at the sites surveyed and the area affected by spreading ranged from less than 20 m to over 200 m from the river. Despite the highly non-uniform displacements, four characteristic distribution patterns including large, distributed ground displacements, block-type movements, large and localized ground displacements, and areas of little to no displacements were identified. Available geotechnical, seismic, and topographic data were collated at the ground surveying sites for subsequent analysis of field measurements. Two widely-used empirical models (Zhang et al. (2004), Youd et al. (2002)) were scrutinized and applied to locations in the vicinity of field measurements for comparison with model predictions. The results indicated generally poor correlation (outside a factor of two) with empirical predictions at most locations and further validated the need for an improved, analysis- based method of predicting lateral displacements that considers the many factors involved on a site-specific basis. In addition, the development of appropriate model input parameters for the Youd et al. (2002) model led to a site-specific correlation of soil behavior type index, Ic, and fines content, FC, for sites along the Avon River in Christchurch that matched up well with existing Ic – FC relationships commonly used in current practice. Lastly, a rigorous analysis was performed for 25 selected locations of ground surveying measurements along the Avon River where ground slope conditions are mild (-1 to 2%) and channel heights range from about 2 – 4.5 m. The field data was divided into categories based on the observed distribution pattern of ground displacements including: large and distributed, moderate and distributed, small to negligible, and large and localized. A systematic approach was applied to determine potential critical layers contributing to the observed displacement patterns which led to the development of characteristic profiles for each category considered. The results of these analyses outline an alternative approach to the evaluation of lateral spreading in which a detailed geotechnical analysis is used to identify the potential for large spreading displacements and likely spatial distribution patterns of spreading. Key factors affecting the observed magnitude and distribution of spreading included the thickness of the critical layer, relative density, soil type and layer continuity. It was found that the large and distributed ground displacements were associated with a thick (1.5 – 2.5 m) deposit of loose, fine to silty sand (qc1 ~4-7 MPa, Ic 1.9-2.1, qc1n_cs ~50-70) that was continuous along the bank and with distance from the river. In contrast, small to negligible displacements were characterized by an absence of or relatively thin (< 1 m), discontinuous critical layer. Characteristic features of the moderate and distributed displacements were found to be somewhere between these two extremes. The localized and large displacements showed a characteristic critical layer similar to that observed in the large and distributed sites but that was not continuous and hence leading to the localized zone of displacement. The findings presented in this thesis illustrate the highly complex nature of lateral displacements that cannot be captured in simplified models but require a robust geotechnical analysis similar to that performed for this research.
The 2010 Darfield and 2011 Christchurch Earthquakes triggered extensive liquefaction-induced lateral spreading proximate to streams and rivers in the Christchurch area, causing significant damage to structures and lifelines. A case study in central Christchurch is presented and compares field observations with predicted displacements from the widely adopted empirical model of Youd et al. (2002). Cone penetration testing (CPT), with measured soil gradation indices (fines content and median grain size) on typical fluvial deposits along the Avon River were used to determine the required geotechnical parameters for the model input. The method presented attempts to enable the adoption of the extensive post-quake CPT test records in place of the lower quality and less available Standard Penetration Test (SPT) data required by the original Youd model. The results indicate some agreement between the Youd model predictions and the field observations, while the majority of computed displacements error on the side of over-prediction by more than a factor of two. A sensitivity analysis was performed with respect to the uncertainties used as model input, illustrating the model’s high sensitivity to the input parameters, with median grain size and fines content among the most influential, and suggesting that the use of CPT data to quantify these parameters may lead to variable results.
Geologic phenomena produced by earthquake shaking, including rockfalls and liquefaction features, provide important information on the intensity and spatiotemporal distribution of earthquake ground motions. The study of rockfall and liquefaction features produced in contemporary well- instrumented earthquakes increases our knowledge of how natural and anthropogenic environments respond to earthquakes and improves our ability to deduce seismologic information from analogous pre-contemporary (paleo-) geologic features. The study of contemporary and paleo- rockfall and liquefaction features enables improved forecasting of environmental responses to future earthquakes. In this thesis I utilize a combination of field and imagery-based mapping, trenching, stratigraphy, and numerical dating techniques to understand the nature and timing of rockfalls (and hillslope sedimentation) and liquefaction in the eastern South Island of New Zealand, and to examine the influence that anthropogenic activity has had on the geologic expressions of earthquake phenomena. At Rapaki (Banks Peninsula, NZ), field and imagery-based mapping, statistical analysis and numerical modeling was conducted on rockfall boulders triggered by the fatal 2011 Christchurch earthquakes (n=285) and compared with newly identified prehistoric (Holocene and Pleistocene) boulders (n=1049) deposited on the same hillslope. A significant population of modern boulders (n=26) travelled farther downslope (>150 m) than their most-travelled prehistoric counterparts, causing extensive damage to residential dwellings at the foot of the hillslope. Replication of prehistoric boulder distributions using 3-dimensional rigid body numerical models requires the application of a drag-coefficient, attributed to moderate to dense slope vegetation, to account for their spatial distribution. Radiocarbon dating provides evidence for 17th to early 20th century deforestation at the study site during Polynesian and European colonization and after emplacement of prehistoric rockfalls. Anthropocene deforestation enabled modern rockfalls to exceed the limits of their prehistoric predecessors, highlighting a shift in the geologic expression of rockfalls due to anthropogenic activity. Optical and radiocarbon dating of loessic hillslope sediments in New Zealand’s South Island is used to constrain the timing of prehistoric rockfalls and associated seismic events, and quantify spatial and temporal patterns of hillslope sedimentation including responses to seismic and anthropogenic forcing. Luminescence ages from loessic sediments constrain timing of boulder emplacement to between ~3.0 and ~12.5 ka, well before the arrival of Polynesians (ca AD 1280) and Europeans (ca AD 1800) in New Zealand, and suggest loess accumulation was continuing at the study site until 12-13 ka. Large (>5 m3) prehistoric rockfall boulders preserve an important record of Holocene hillslope sedimentation by creating local traps for sediment aggradation and upbuilding soil formation. Sediment accumulation rates increased considerably (>~10 factor increase) following human arrival and associated anthropogenic burning of hillslope vegetation. New numerical ages are presented to place the evolution of loess-mantled hillslopes in New Zealand’s South Island into a longer temporal framework and highlight the roles of earthquakes and humans on hillslope surface process. Extensive field mapping and characterization for 1733 individual prehistoric rockfall boulders was conducted at Rapaki and another Banks Peninsula site, Purau, to understand their origin, frequency, and spatial and volumetric distributions. Boulder characteristics and distributions were compared to 421 boulders deposited at the same sites during the 2010-2011 Canterbury earthquake sequence. Prehistoric boulders at Rapaki and Purau are comprised of two dominant lithofacies types: volcanic breccia and massive (coherent) lava basalt. Volcanic breccia boulders are found in greatest abundance (64-73% of total mapped rockfall) and volume (~90-96% of total rockfall) at both locations and exclusively comprise the largest boulders with the longest runout distances that pose the greatest hazard to life and property. This study highlights the primary influence that volcanic lithofacies architecture has on rockfall hazard. The influence of anthropogenic modifications on the surface and subsurface geologic expression of contemporary liquefaction created during the 2010-2011 Canterbury earthquake sequence (CES) in eastern Christchurch is examined. Trench observations indicate that anthropogenic fill layer boundaries and the composition/texture of discretely placed fill layers play an important role in absorbing fluidized sand/silt and controlling the subsurface architecture of preserved liquefaction features. Surface liquefaction morphologies (i.e. sand blows and linear sand blow arrays) display alignment with existing utility lines and utility excavations (and perforated pipes) provided conduits for liquefaction ejecta during the CES. No evidence of pre-CES liquefaction was identified within the anthropogenic fill layers or underlying native sediment. Radiocarbon dating of charcoal within the youngest native sediment suggests liquefaction has not occurred at the study site for at least the past 750-800 years. The importance of systematically examining the impact of buried infrastructure on channelizing and influencing surface and subsurface liquefaction morphologies is demonstrated. This thesis highlights the importance of using a multi-technique approach for understanding prehistoric and contemporary earthquake phenomena and emphasizes the critical role that humans play in shaping the geologic record and Earth’s surface processes.
None
Recent earthquakes have shown that liquefaction and associated ground deformations are major geotechnical hazards to civil engineering infrastructures, such as pipelines. In particular, sewer pipes have been damaged in many areas in Christchurch as a result of liquefaction-induced lateral spreading near waterways and ground oscillation induced by seismic shaking. In this paper, the addition of a flexible AM liner as a potential countermeasure to increase sewer pipe capacity was investigated. Physical testing through 4-point loading test was undertaken to characterise material properties and the response of both unlined pipe and its lined counterpart. Next, numerical models were created using SAP2000 and ABAQUS to analyse buried pipeline response to transverse permanent ground displacement and to quantify, over a range of pipe segment lengths and soil parameters, the effectiveness of the AM liner in increasing displacement capacity. The numerical results suggest that the addition of the AM liner increases the deformation capacity of the unlined sewer pipe by as much as 50 times. The results confirmed that AM liner is an effective countermeasure for sewer pipes in liquefied ground not only in terms of increased deformation capacity but also the fact that AM-Liner can prevent influx of sand and water through broken pipes, making sewer pipes with liner remaining serviceable even under severe liquefaction condition.
Gravelly soils’ liquefaction has been documented since early 19th century with however the focus being sand-silts mixture – coarse documentation of this topic, that gravels do in fact liquefy was only acknowledged as an observation. With time, we have been impacted by earthquakes, EQ causing more damage to our property: life and environment-natural and built. In this realm of EQ related-damage the ground or soils in general act as buffer between the epicentre and the structures at a concerned site. Further, in this area, upon the eventual acknowledgement of liquefaction of soils as a problem, massive efforts were undertaken to understand its mechanics, what causes and thereby how to mitigate its ill-effects. Down that lane, gravelly soils’ liquefaction was another milestone covered in early 20th century, and thus regarded as a problematic subject. This being a fairly recent acknowledgement, efforts have initiated in this direction (or area of particle size under consideration-gravels>2mm), with this research outputs intended to complement that research for the betterment of our understanding of what’s happening and how shall we best address it, given the circumstances: socio (life) - environment (structures) - economic (cost or cost-“effectiveness’) and of course political (our “willingness” to want to address the problem). Case histories from at least 29 earthquakes worldwide have indicated that liquefaction can occur in gravelly soils (both in natural deposits and manmade reclamations) inducing large ground deformation and causing severe damage to civil infrastructures. However, the evaluation of the liquefaction resistance of gravelly soils remains to be a major challenge in geotechnical earthquake engineering. To date, laboratory tests aimed at evaluating the liquefaction resistance of gravelly soils are still very limited, as compared to the large body of investigations carried out on assessing the liquefaction resistance of sandy soils. While there is a general agreement that the liquefaction resistance of gravelly soils can be as low as that of clean sands, previous studies suggested that the liquefaction behaviour of gravelly soils is significantly affected by two key factors, namely relative density (Dr) and gravel content (Gc). While it is clear that the liquefaction resistance of gravels increases with the increasing Dr, there are inconclusive and/or contradictory results regarding the effect of Gc on the liquefaction resistance of gravelly soils. Aimed at addressing this important topic, an investigation is being currently carried out by researchers at the University of Canterbury, UC. As a first step, a series of undrained cyclic triaxial tests were conducted on selected sand-gravel mixtures (SGMs), and inter-grain state framework concepts such as the equivalent and skeleton void ratios were used to describe the joint effects of Gc and Dr on the liquefaction resistance of SGMs. Following such experimental effort, this study is aimed at providing new and useful insights, by developing a critical state-based method combined with the inter-grain state framework to uniquely describe the liquefaction resistance of gravelly soils. To do so, a series of monotonic drained triaxial tests will be carried out on selected SGMs. The outcomes of this study, combined with those obtained to date by UC researchers, will greatly contribute to the expansion of a worldwide assessment database, and also towards the development of a reliable liquefaction triggering procedure for characterising the liquefaction potential of gravelly soils, which is of paramount importance not only for the New Zealand context, but worldwide. This will make it possible for practising engineers to identify liquefiable gravelly soils in advance and make sound recommendations to minimise the impact of such hazards on land, and civil infrastructures.
Using case studies from the 2010-2011 Canterbury, New Zealand earthquake sequence, this study assesses the accuracies of paleoliquefaction back-analysis methods and explores the challenges, techniques, and uncertainties associated with their application. While liquefaction-based back-analyses have been widely used to estimate the magnitudes of paleoearthquakes, their uncertain efficacies continue to significantly affect the computed seismic hazard in regions where they are relied upon. Accordingly, their performance is evaluated herein using liquefaction data from modern earthquakes with known magnitudes. It is shown that when the earthquake source location and mechanism are known, back-analysis methods are capable of accurately deriving seismic parameters from liquefaction evidence. However, because the source location and mechanism are often unknown in paleoseismic studies, and because accurate interpretation is shown to be more difficult in such cases, new analysis techniques are proposed herein. An objective parameter is proposed to geospatially assess the likelihood of any provisional source location, enabling an analyst to more accurately estimate the magnitude of a liquefaction-inducing paleoearthquake. This study demonstrates the application of back-analysis methods, provides insight into their potential accuracies, and provides a framework for performing paleoliquefaction analyses worldwide.
Motivation This poster aims to present fragility functions for pipelines buried in liquefaction-prone soils. Existing fragility models used to quantify losses can be based on old data or use complex metrics. Addressing these issues, the proposed functions are based on the Christchurch network and soil and utilizes the Canterbury earthquake sequence (CES) data, partially represented in Figure 1. Figure 1 (a) presents the pipe failure dataset, which describes the date, location and pipe on which failures occurred. Figure 1 (b) shows the simulated ground motion intensity median of the 22nd February 2011 earthquake. To develop the model, the network and soil characteristics have also been utilized.
This paper presents a methodology by which both site-specific and spatially distributed ground motion intensity can be obtained immediately following an earthquake event. The methodology makes use of both prediction models for ground motion intensity and its correlation over spatial distances. A key benefit of the methodology is that the ground motion intensity at a given location is not a single value but a distribution of values. The distribution is comprised of both a mean and also standard deviation, with the standard deviation being a function of the distance to nearby strong motion stations. The methodology is illustrated for two applications. Firstly, maps of conditional peak ground acceleration (PGA) have been developed for the major events in the Canterbury earthquake sequence. It is illustrated how these conditional maps can be used for post-event evaluation of liquefaction triggering criteria which have been adopted by the Department of Building and Housing (DBH). Secondly, the conditional distribution of response spectral ordinates is obtained at a specific location for the purposes of determining appropriate ground motion records for use in seismic response analyses of important structures at locations where direct recordings are absent.
There is a growing awareness of the need for the earthquake engineering practice to incorporate in addition to empirical approaches in evaluation of liquefaction hazards advanced methods which can more realistically represent soil behaviour during earthquakes. Currently, this implementation is hindered by a number of challenges mainly associated with the amount of data and user-experience required for such advanced methods. In this study, we present key steps of an advanced seismic effective-stress analysis procedure, which on the one hand can be fully automated and, on the other hand, requires no additional input (at least for preliminary applications) compared to simplified cone penetration test (CPT)-based liquefaction procedures. In this way, effective-stress analysis can be routinely applied for quick, yet more robust estimations of liquefaction hazards, in a similar fashion to the simplified procedures. Important insights regarding the dynamic interactions in liquefying soils and the actual system response of a deposit can be gained from such analyses, as illustrated with the application to two sites from Christchurch, New Zealand.
In the period between September 2010 and December 2011, Christchurch was shaken by a series of strong earthquakes including the MW7.1 4 September 2010, Mw 6.2 22 February 2011, MW6.2 13 June 2011 and MW6.0 23 December 2011 earthquakes. These earthquakes produced very strong ground motions throughout the city and surrounding areas that resulted in soil liquefaction and lateral spreading causing substantial damage to buildings, infrastructure and the community. The stopbank network along the Kaiapoi and Avon River suffered extensive damage with repairs projected to take several years to complete. This presented an opportunity to undertake a case-study on a regional scale of the effects of liquefaction on a stopbank system. Ultimately, this information can be used to determine simple performance-based concepts that can be applied in practice to improve the resilience of river protection works. The research presented in this thesis draws from data collected following the 4th September 2010 and 22nd February 2011 earthquakes. The stopbank damage is categorised into seven key deformation modes that were interpreted from aerial photographs, consultant reports, damage photographs and site visits. Each deformation mode provides an assessment of the observed mechanism of failure behind liquefaction-induced stopbank damage and the factors that influence a particular style of deformation. The deformation modes have been used to create a severity classification for the whole stopbank system, being ‘no or low damage’ and ‘major or severe damage’, in order to discriminate the indicators and factors that contribute to ‘major to severe damage’ from the factors that contribute to all levels of damage a number of calculated, land damage, stopbank damage and geomorphological parameters were analysed and compared at 178 locations along the Kaiapoi and Avon River stopbank systems. A critical liquefiable layer was present at every location with relatively consistent geotechnical parameters (cone resistance (qc), soil behaviour type (Ic) and Factor of Safety (FoS)) across the study site. In 95% of the cases the critical layer occurred within two times the Height of the Free Face (HFF,). A statistical analysis of the geotechnical factors relating to the critical layer was undertaken in order to find correlations between specific deformation modes and geotechnical factors. It was found that each individual deformation mode involves a complex interplay of factors that are difficult to represent through correlative analysis. There was, however, sufficient data to derive the key factors that have affected the severity of deformation. It was concluded that stopbank damage is directly related to the presence of liquefaction in the ground materials beneath the stopbanks, but is not critical in determining the type or severity of damage, instead it is merely the triggering mechanism. Once liquefaction is triggered it is the gravity-induced deformation that causes the damage rather than the shaking duration. Lateral spreading and specifically the depositional setting was found to be the key aspect in determining the severity and type of deformation along the stopbank system. The presence or absence of abandoned or old river channels and point bar deposits was found to significantly influence the severity and type of deformation. A review of digital elevation models and old maps along the Kaiapoi River found that all of the ‘major to severe’ damage observed occurred within or directly adjacent to an abandoned river channel. Whilst a review of the geomorphology along the Avon River showed that every location within a point bar deposit suffered some form of damage, due to the depositional environment creating a deposit highly susceptible to liquefaction.
Advanced seismic effective-stress analysis is used to scrutinize the liquefaction performance of 55 well-documented case-history sites from Christchurch. The performance of these sites during the 2010-2011 Canterbury earthquake sequence varied significantly, from no liquefaction manifestation at the ground surface (in any of the major events) to severe liquefaction manifestation in multiple events. For the majority of the 55 sites, the simplified liquefaction evaluation procedures, which are conventionally used in engineering practice, could not explain these dramatic differences in the manifestation. Detailed geotechnical characterization and subsequent examination of the soil profile characteristics of the 55 sites identified some similarities but also important differences between sites that manifested liquefaction in the two major events of the sequence (YY-sites) and sites that did not manifest liquefaction in either event (NN-sites). In particular, while the YY-sites and NN-sites are shown to have practically identical critical layer characteristics, they have significant differences with regard to their deposit characteristics including the thickness and vertical continuity of their critical zones and liquefiable materials. A CPT-based effective stress analysis procedure is developed and implemented for the analyses of the 55 case history sites. Key features of this procedure are that, on the one hand, it can be fully automated in a programming environment and, on the other hand, it is directly equivalent (in the definition of cyclic resistance and required input data) to the CPT-based simplified liquefaction evaluation procedures. These features facilitate significantly the application of effective-stress analysis for simple 1D free-field soil-column problems and also provide a basis for rigorous comparisons of the outcomes of effective-stress analyses and simplified procedures. Input motions for the analyses are derived using selected (reference) recordings from the two major events of the 2010-2011 Canterbury earthquake sequence. A step-by-step procedure for the selection of representative reference motions for each site and their subsequent treatment (i.e. deconvolution and scaling) is presented. The focus of the proposed procedure is to address key aspects of spatial variability of ground motion in the near-source region of an earthquake including extended-source effects, path effects, and variation in the deeper regional geology.
The Canterbury Earthquake Sequence (CES) of 2010-2011 caused widespread liquefaction in many parts of Christchurch. Observations from the CES highlight some sites were liquefaction was predicted by the simplified method but did not manifest. There are a number of reasons why the simplified method may over-predict liquefaction, one of these is the dynamic interaction between soil layers within a stratified deposit. Soil layer interaction occurs through two key mechanisms; modification of the ground motion due to seismic waves passing through deep liquefied layers, and the effect of pore water seepage from an area of high excess pore water pressure to the surrounding soil. In this way, soil layer interaction can significantly alter the liquefaction behaviour and surface manifestation of soils subject to seismic loading. This research aimed to develop an understanding of how soil layer interaction, in particular ground motion modification, affects the development of excess pore water pressures and liquefaction manifestation in a soil deposit subject to seismic loading. A 1-D soil column time history Effective Stress Analysis (ESA) was conducted to give an in depth assessment of the development of pore pressures in a number of soil deposits. For this analysis, ground motions, soil profiles and model parameters were required for the ESA. Deconvolution of ground motions recorded at the surface during the CES was used to develop some acceleration time histories to input at the base of the soil-column model. An analysis of 55 sites around Christchurch, where detailed site investigations have been carried out, was then conducted to identify some simplified soil profiles and soil characteristics. From this analysis, four soil profiles representative of different levels of liquefaction manifestation were developed. These were; two thick uniform and vertically continuous sandy deposits that were representative of sites were liquefaction manifested in both the Mw 7.1 September 2010 and the Mw 6.3 February 2011 earthquakes, and two vertically discontinuous profiles with interlayered liquefiable and non-liquefiable layers representative of sites that did not manifest liquefaction in either the September 2010 or the February 2011 events. Model parameters were then developed for these four representative soil profiles through calibration of the constitutive model in element test simulations. Simulations were run for each of the four profiles subject to three levels of loading intensity. The results were analysed for the effect of soil layer interaction. These were then compared to a simplified triggering analysis for the same four profiles to determine where the simplified method was accurate in predicting soil liquefaction (for the continuous sandy deposits) and were it was less accurate (the vertically discontinuous deposits where soil layer interaction was a factor).
The majority of current procedures used to deduce liquefaction potential of soils rely on empirical methods. These methods have been proven to work in the past, but these methods are known to overestimate the liquefaction potential in certain regions of Christchurch due to a whole range of factors, and the theoretical basis behind these methods cannot be explained scientifically. Critical state soil mechanics theory was chosen to provide an explanation for the soil's behaviour during the undrained shearing. Soils from two sites in Christchurch were characterised at regular intervals for the critical layers and tested for the critical state lines (CSL). Various models and relationships were then used to predict the CSL and compared with the actual CSL. However none of the methods used managed to predict the CSL accurately, and a separate Christchurch exclusive relationship was proposed. The resultant state parameter values could be obtained from shear-wave velocity plots and were then developed into cyclic resistance ratios (CRR). These were subsequently compared with cyclic stress ratios (CSR) from recent Christchurch earthquakes to obtain the factor of safety. This CSL-based approach was compared with other empirical methods and was shown to yield a favourable relationship with visual observations at the sites' locations following the earthquake.
This paper provides a summary of the ground motions observed in the recent Canterbury, New Zealand earthquake sequence. The sequence occurred in a region of relatively moderate seismicity, 130km to the east of the Alpine Fault, the major plate-boundary in the region. From an engineering perspective, the sequence has been primarily comprised of the initial 04/09/2010 Darfield earthquake (Mw7.1) followed by the 22/02/2011 Christchurch earthquake (Mw6.3), and two aftershocks on 13/06/ 2011 (Mw5.3 and 6.0, respectively). The dense spacing of strong motions in the region, and their close proximity to the respective causative faults, has resulted in strong ground motions far exceeding the previous catalogue of strong motion observed in New Zealand. The observed ground motions have exhibited clear evidence of: (i) near-source directivity; (ii) sedimentary basin focusing, amplification and basin effect refraction; (iii) non-linear site response; (iv) cyclic mobility postliquefaction; and (v) extreme vertical ground motions exceeding 2g, among others.
Geospatial liquefaction models aim to predict liquefaction using data that is free and readily-available. This data includes (i) common ground-motion intensity measures; and (ii) geospatial parameters (e.g., among many, distance to rivers, distance to coast, and Vs30 estimated from topography) which are used to infer characteristics of the subsurface without in-situ testing. Since their recent inception, such models have been used to predict geohazard impacts throughout New Zealand (e.g., in conjunction with regional ground-motion simulations). While past studies have demonstrated that geospatial liquefaction-models show great promise, the resolution and accuracy of the geospatial data underlying these models is notably poor. As an example, mapped rivers and coastlines often plot hundreds of meters from their actual locations. This stems from the fact that geospatial models aim to rapidly predict liquefaction anywhere in the world and thus utilize the lowest common denominator of available geospatial data, even though higher quality data is often available (e.g., in New Zealand). Accordingly, this study investigates whether the performance of geospatial models can be improved using higher-quality input data. This analysis is performed using (i) 15,101 liquefaction case studies compiled from the 2010-2016 Canterbury Earthquakes; and (ii) geospatial data readily available in New Zealand. In particular, we utilize alternative, higher-quality data to estimate: locations of rivers and streams; location of coastline; depth to ground water; Vs30; and PGV. Most notably, a region-specific Vs30 model improves performance (Figs. 3-4), while other data variants generally have little-to-no effect, even when the “standard” and “high-quality” values differ significantly (Fig. 2). This finding is consistent with the greater sensitivity of geospatial models to Vs30, relative to any other input (Fig. 5), and has implications for modeling in locales worldwide where high quality geospatial data is available.
This poster provides a comparison between the strong ground motions observed in the 22 February 2011 Mw6.3 Christchurch earthquake with those observed in Tokyo during the 11 March 2011 Mw9.0 Tohoku earthquake. The destuction resulting from both of these events has been well documented, although tsunami was the principal cause of damage in the latter event, and less attention has been devoted to the impact of earthquake-induced ground motions. Despite Tokyo being located over 100km from the nearest part of the causative rupture, the ground motions observed from the Tohoku earthquake were significant enough to cause structural damage and also significant liquefaction to loose reclaimed soils in Tokyo Bay. The author was fortunate enough (from the perspective of an earthquake engineer) to experience first-hand both of these events. Following the Tohoku event, the athor conducted various ground motion analyses and reconniassance of the Urayasu region in Tokyo Bay affected by liquefaction in collaboration with Prof. Kenji Ishihara. This conference is therefore a fitting opportunity in which to discuss some of authors insights obtained as a result of this first hand knowledge. Figure 1 illustrates the ground motions recorded in the Christchurch CBD in the 22 February 2011 and 4 September 2010 earthquakes, with that recorded in Tokyo Bay in the 11 March 2011 Tohoku earthquake. It is evident that these three ground motions vary widely in their amplitude and duration. The CBGS ground motion from the 22 February 2011 event has a very large amplitude (nearly 0.6g) and short duration (approx. 10s of intense shaking), as a result of the causal Mw6.3 rupture at short distance (Rrup=4km). The CBGS ground motion from the 4 September 2010 earthquake has a longer duration (approx. 30s of intense shaking), but reduced acceleration amplitude, as a result of the causal Mw7.1 rupture at a short-to-moderate distance (Rrup=14km). Finally, the Urayasu ground motion in Tokyo bay during the 11 March 2011 Tohoku earthquake exhibits an acceleration amplitude similar to the 4 September 2010 CBGS ground motion, but a significantly larger duration (approx 150s of intense shaking). Clearly, these three different ground motions will affect structures and soils in different ways depending on the vibration characteristics of the structures/soil, and the potential for strength and stiffness degradation due to cumulative effects. Figure 2 provides a comparison between the arias intensities of the several ground motion records from the three different events. It can be seen that the arias intensities of the ground motions in the Christchurch CBD from the 22 February 2011 earthquake (which is on average AI=2.5m/s) is approximately twice that from the 4 September 2010 earthquake (average AI≈1.25). This is consistent with a factor of approximately 1.6 obtained by Cubrinovski et al. (2011) using the stress-based (i.e.PGA-MSF) approach of liquefaction triggering. It can also be seen that the arias intensity of the ground motions recorded in Tokyo during the 2011 Tohoku earthquake are larger than ground motions in the Christchurch CBD from the 4 September 2011 earthquake, but smaller than those of the 22 February 2011 earthquake. Based on the arias intensity liquefaction triggering approach it can therefore be concluded that the ground motion severity, in terms of liquefaction potential, for the Tokyo ground motions is between those ground motions in Christchurch CBD from the 4 September 2010 and 22 February 2011 events.
Results from a series of 1D seismic effective stress analyses of natural soil deposits from Christchurch are summarized. The analysed soil columns include sites whose performance during the 2010-2011 Canterbury earthquakes varied significantly, from no liquefaction manifestation at the ground surface to very severe liquefaction, in which case a large area of the site was covered by thick soil ejecta. Key soil profile characteristics and response mechanisms affecting the severity of surface liquefaction manifestation and subsequent damage are explored. The influence of shaking intensity on the triggering and contribution of these mechanisms is also discussed. Careful examination of the results highlights the importance of considering the deposit as a whole, i.e. a system of layers, including interactions between layers in the dynamic response and through pore water pressure redistribution and water flow.
In 2010 and 2011 a series of earthquakes hit the central region of Canterbury, New Zealand, triggering widespread and damaging liquefaction in the area of Christchurch. Liquefaction occurred in natural clean sand deposits, but also in silty (fines-containing) sand deposits of fluvial origin. Comprehensive research efforts have been subsequently undertaken to identify key factors that influenced liquefaction triggering and severity of its manifestation. This research aims at evaluating the effects of fines content, fabric and layered structure on the cyclic undrained response of silty soils from Christchurch using Direct Simple Shear (DSS) tests. This poster outlines preliminary calibration and verification DSS tests performed on a clean sand to ensure reliability of testing procedures before these are applied to Christchurch soils.
