A flow chart depicting potential hazards earthquakes pose to coastal geography.
A digitally manipulated image of printer's type, spelling out "Safe Dust". The photographer comments, "After the September 2010 Christchurch earthquake liquefaction poured out of the ground mostly in the East of Christchurch. This silt, which was a form of sand was declared safe and would not harm gardens if it was spread around in moderation. After the February 2011 quake as a result of even more liquefaction and the sewers being ruptured, the liquefaction was declared as toxic. People clearing it up should wear a mask, boots and gloves especially when it had dried up and become dusty. This just so happens to be the words found on an old printing press".
When the 2010 and 2011 earthquakes created a city-wide outdoor research laboratory, UC Civil Engineering Professor Misko Cubrinovski gathered as much information as possible. This work has been recognised by the American Society of Civil Engineers (ASCE), which is presenting him with the 2019 Ralph B. Peck Award for "outstanding contributions to the geotechnical engineering profession through the publication of several insightful field case histories"
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.
Taken soon after the 6.3 magnatude quake hit Christchurch 22 February 2011. A friend tried to get home when the road opened up and her 4 wheel drive went down. We could not get it out! The grey liquefaction has stuck the other car fast. It took less than an hour to grow this thick.
After the liquefaction has been cleaned up we have piles of the stuff on our streets waiting for the council to pick it up. This little Gnome is obviously doing its bit as Christchurch cleans up after the 6.3 magnitude quake 22 February 2011.
None
Photograph captioned by Fairfax, "Damage from the February 22nd earthquake in Christchurch. Liquefaction at the St Martins Shopping Centre and New World supermarket in St Martins, Christchurch".
Photograph captioned by Fairfax, "Damage from the February 22nd earthquake in Christchurch. Liquefaction at the St Martins Shopping Centre and New World supermarket in St Martins, Christchurch".
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.
This report describes the earthquake hazard in Waimate and Mackenzie districts and the part of Waitaki district within Canterbury, and gives details of historic earthquakes. It includes district-scale (1:500,000) active fault, ground shaking zone, liquefaction and landslide susceptibility maps. The report describes earthquake scenarios for a magnitude 7.2-7.4 Ostler Fault earthquake near Twizel, a magnitude 8 Alpine Fault earthquake, and a magnitude 6.9 Hunters Hills Fault Zone earthquake near Waimate. See Object Overview for background and usage information.
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.
Dried liquefaction silt in North New Brighton. The photographer comments, "This is the result of liquefaction which spewed out after the double earthquake in Christchurch. Having flowed into a shallow depression that was deep enough for a fair quantity of the silty liquid to settle and separate: the heavy sand below and a talcum powder like substance on top. Some of these are so delicate that a mouse crossing them would probably crack them. Here the sun has dried them out and they have contracted and curled up towards their centres".
A pothole in a road surface, showing tyre marks where a vehicle has driven through the hole. The photographer comments, "After the earthquake in Christchurch in February 2011 burst underground pipes and liquefaction caused unseen hollows under the road surfaces. Occasionally after all the rest have been exposed by traffic someone would find 'discover' a new one".
A contractor with a roller stands in an enormous hole, and shouts up to a family standing at the top that it's now ready to be built on. Refers to efforts to compact soil in Christchurch to prevent further liquefaction following the Canterbury earthquakes. Quantity: 1 digital cartoon(s).
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.
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).
In a Christchurch street still covered in liquefaction man weeps over his four-wheel drive car, which has a number plate showing the word 'macho'. Two women who are working at clearing the road of silt watch and one of them observes that 'he's not handling the quake well at all... Keeps getting silt smears on the 4 x 4!' Context - The Christchurch earthquake of 22 February 2011. Liquefaction is a particular problem. There is a point being made here about the 'macho' man who sobs over his car and the two staunch women who get on with the cleaning-up effort. Quantity: 1 digital cartoon(s).
A geotechnical expert from the United States has told the Royal Commission into the Canterbury Earthquakes, developers should be required to submit soil reports before building on land prone 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.
A residential street in New Brighton. Liquefaction still lines the street, and lampposts are leaning in different directions. The photographer comments, "This is the New Brighton red zone, which is parallel to the Avon River. The area suffered serious liquefaction during the numerous earthquakes/aftershocks and the land is being bought by the government. Although the houses do not look too bad in the background they have suffered badly. On the day I took this picture the council had just hours before cut the grass, which made the area look less abandoned".
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.
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.
Workers operate a drilling rig, sampling soil as part of EQC's geotechnical investigation of TC3 land. The photographer comments, "The work of getting 'soil' samples from all the areas marked as green/blue zones in Christchurch. These areas may be susceptible to liquefaction if a major earthquake occurs. The soil samples were a failure as all they found was sand".
A very large woman stands wedged between two rows of concrete pillars eating a huge cream bun. She says 'Christmas fare protection... it helps prevent lateral spread!' Context - overeating at Christmas and lateral spreading, which is associated with liquefaction and tends to occur near streams and waterways as the soil mass moves towards them. Reference to the Christchurch earthquake of 4th September 2010. Quantity: 1 digital cartoon(s).
The last pile of liquefied silt to be cleaned up, from our lawn.
Workers inspect a broken sewerage line in New Brighton. The photographer comments, "After the Christchurch earthquake on 23 December 2011 the sewer pipe got badly damaged at New Brighton and was leaking into the Avon River. I think the guy was worried about the fast flow causing him to fill his boots rather than the depth".
Misko Cubrinovski, Civil Engineer, photographed with liquefaction and lateral spreading on Oxford Terrace. Misko's area is geotechnical earthquake engineering and foundation engineering, and he will feature in a UC in the News pull out supplement inThe Press.
Page 7 of Section A of the Christchurch Press, published on Wednesday 15 June 2011.