Visualisation of the 2010 Christchurch Boxing Day aftershock
Research papers, University of Canterbury Library
None
None
Improving community resilience requires a way of thinking about the nature of a community. Two complementary aspects are proposed: the flows connecting the community with its surrounding environment and the resources the community needs for its ongoing life. The body of necessary resources is complex, with many interactions between its elements. A systems approach is required to understand the issues adequately. Community resilience is discussed in general terms together with strategies for improving it. The ideas are then illustrated and amplified by an extended case study addressing means of improving the resilience of a community on the West Coast of New Zealand to natural disasters. The case study is in two phases. The first relies on a mix of on-the-ground observations and constructed scenarios to provide recommendations for enhancing community resilience, while the second complements the first by developing a set of general lessons and issues to be addressed from observations of the Christchurch earthquakes of 2010 and 2011.
This paper summarizes the development of a high-resolution surficial shear wave velocity model based on the combination of the large high-spatial-density database of cone penetration test (CPT) logs in and around Christchurch, New Zealand and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. This near-surface shear wave velocity model has applications for site characterization efforts via the development of maps of time-averaged shear wave velocities over specific depths, as well as use in site response analysis and ground motion simulation.
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.
Recent field investigations were carried out to define the shear wave velocity (VS) profile and site periods across the Canterbury region, supplementing earlier efforts in urban Christchurch. Active source surface wave testing, ambient wave field (passive) and H/V spectral ratio methods were used to characterise the soil profile in the region. H/V spectral ratio peaks indicate site periods in the range of 5-7 seconds across much of the Canterbury Plains, broadly consistent with those based on a 1D velocity model for the region. Site periods decrease rapidly in the vicinity of the Canterbury foothills and the Banks Peninsula outcrops. In Christchurch, the Riccarton Gravels result in a significant mode of vibration that has a much shorter period than the site period of the entire soil column down to basement rock.
This presentation summarizes the development of high-resolution surficial soil velocity models in the Canterbury, New Zealand basin. Shallow (<30m) shear wave velocities were primarily computed based on a combination of a large database of over 15,000 cone penetration test (CPT) logs in and around Christchurch, and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. Large active-source testing at 22 locations and ambient-wavefield surface wave and H/V testing at over 80 locations were utilized in combination with 1700 water well logs to constrain the inter-bedded stratigraphy and velocity of Quaternary sediments up to depths of several hundred meters. Finally, seismic reflection profiles and the ambient-wavefield surface wave data provide constraint on velocities from several hundred meters to several kilometres. At all depths, the high resolution data illustrates the complexity of the soil conditions in the region, and the developed 3D models are presently being used in broadband ground motion simulations to further interpret the observed strong ground motions in the 2010-2011 Canterbury earthquake sequence.
Children are often overlooked in the aftermath of a natural disaster, and children’s use of coping strategies plays an important part in their post-disaster adaptation (Vernberg, La Greca, Silverman, & Prinstein, 1996). The aim of this qualitative study was to explore the coping strategies of children with adequate self-regulation skills and minimal behaviour problems, living in Christchurch following the major 2010 and 2011 earthquakes. This aim was achieved through the use of semi-structured interviews with five seven-year-old children, their parents, and their teachers. These interviews were analysed using Directed Content Analysis and results showed that children most often reported using active and adaptive coping strategies, followed by avoidant strategies. Results in the current literature regarding children’s coping suggest that children exposed to natural disasters are able to utilise strategies that involve some personal control over their environment and emotions, through the use of active and adaptive coping strategies. Findings from this study contribute to the current understanding of children’s use of coping strategies when faced with commonly occurring childhood upsets. Further research is required regarding the outcomes associated with the use of effective coping strategies following traumatic events.
This paper summarizes the development of a region-wide surficial shear wave velocity model based on the combination of the large high-spatial-density database of cone penetration test (CPT) logs in and around Christchurch, New Zealand and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. The ongoing development of this near-surface shear wave velocity model has applications for site characterization efforts via the development of maps of time-averaged shear wave velocities over specific depths, and the identification of regional similarities and differences in soil shear stiffness.
None
Current research in geotechnical engineering at the University of Canterbury includes a number of laboratory testing programmes focussed on understanding the behaviour of natural soil deposits in Christchurch during the 2010-2011 Canterbury Earthquake Sequence. Many soils found in Christchurch are sands or silty sands with little to no plasticity, making them very difficult to sample using established methods. The gel-push sampling methodology, developed by Kiso-Jiban Consultants in Japan, was developed to address some of the deficiencies of existing sampling techniques and has been deployed on two projects in Christchurch. Gel push sampling is carried out with a range of samplers which are modified versions of existing technology, and the University of Canterbury has acquired three versions of the tools (GP-S, GP-Tr, GP-D). Soil samples are extracted from the bottom of a freshly drilled borehole and are captured within a liner barrel, close to 1m in length. A lubricating polymer gel coats the outside of the soil sample as it enters the liner barrel. The frictional rubbing which normally occurs on the sides of the soil samples using existing techniques is eliminated by the presence of the polymer gel. The operation of the gel-push samplers is significantly more complicated than conventional push-tube samplers, and in the initial trials a number of operational difficulties were encountered, requiring changes to the sampling procedures. Despite these issues, a number of high quality soil samples were obtained on both projects using the GP-S sampler to capture silty soil. Attempts were made to obtain clean sands using a different gel-push sampler (GP-TR) in the Red Zone. The laboratory testing of these sands indicated that they were being significantly disturbed during the sampling and/or transportation procedures. While it remains too early to draw definitive conclusions regarding the performance of the gel-push samplers, the methodology has provided some promising results. Further trialling of the tools are required to refine operating procedures understand the full range of soil conditions which can be successfully sampled using the tools. In parallel with the gel-push trials, a Dames and Moore fixed-piston sampler has been used by our research partners from Berkeley to obtain soil samples at a number of sites within Christchurch. This sampler features relatively short (50cm), thin-walled liner barrels which is advanced into the ground under the action of hydraulic pressure. By reducing the overall length of the soil being captured, the disturbance to the soil as it enters the liner barrel is significantly reduced. The Dames and Moore sampler is significantly easier to operate than the gel-push sampler, and past experience has shown it to be successful in soft, plastic materials (i.e. clays and silty clays). The cyclic resistance of one silty clay obtained using both the gel-push and Dames & Moore samplers has been found to be very similar, and ongoing research aims to establish whether similar results are obtained for different soil types, including silty materials and clean sands.
In the aftermath of the 2010-2011 Canterbury Earthquake Sequence (CES), the location of Christchurch-City on the coast of the Canterbury Region (New Zealand) has proven crucial in determining the types of- and chains of hazards that impact the city. Very rapidly, the land subsidence of up to 1 m (vertical), and the modifications of city’s waterways – bank sliding, longitudinal profile change, sedimentation and erosion, engineered stop-banks… - turned rainfall and high-tides into unprecedented floods, which spread across the eastern side of the city. Within this context, this contribution presents two modeling results of potential floods: (1) results of flood models and (2) the effects of further subsidence-linked flooding – indeed if another similar earthquake was to strike the city, what could be the scenarios of further subsidence and then flooding. The present research uses the pre- and post-CES LiDAR datasets, which have been used as the boundary layer for the modeling. On top of simple bathtub model of inundation, the river flood model was conducted using the 2-D hydrodynamic code NAYS-2D developed at the University of Hokkaido (Japan), using a depth-averaged resolution of the hydrodynamic equations. The results have shown that the area the most at risk of flooding are the recent Holocene sedimentary deposits, and especially the swamplands near the sea and in the proximity of waterways. As the CES drove horizontal and vertical displacement of the land-surface, the surface hydrology of the city has been deeply modified, increasing flood risks. However, it seems that scientists and managers haven’t fully learned from the CES, and no research has been looking at the potential future subsidence in further worsening subsidence-related floods. Consequently, the term “coastal quake”, coined by D. Hart is highly topical, and most especially because most of our modern cities and mega-cities are built on estuarine Holocene sediments.
The standard way in which disaster damages are measured involves examining separately the number of fatalities, of injuries, of people otherwise affected, and the financial damage that natural disasters cause. Here, we implement a novel way to aggregate these separate measures of disaster impact and apply it to two catastrophic events from 2011: the Christchurch (New Zealand) earthquakes and the Greater Bangkok (Thailand) flood. This new measure, which is similar to the World Health Organization's calculation of Disability Adjusted Life Years (DALYs) lost due to the burden of diseases and injuries, is described in detail in Noy [7]. It allows us to conclude that New Zealand lost 180 thousand lifeyears as a result of the 2011 events, and Thailand lost 2644 thousand lifeyears. In per capita terms, the loss is similar, with both countries losing about 15 days per person due to the 2011 catastrophic events in these two countries. © This manuscript version is made available under the CC-BY-NC-ND 4.0 license https://creativecommons.org/licenses/by-nc-nd/4.0/
This paper presents a seismic velocity model of Canterbury, New Zealand based on 3D geologic surfaces and velocities from a range of data sources. The model provides the 3D crustal structure in the region at multiple length scales for seismic wave propagation simulations, such as broadband ground motion and shallow site response analyses related to understanding the ground motions and site responses during the 2010- 2011 Canterbury earthquakes. Pre-Quaternary geologic horizons are calculated based on the reinterpretation of a comprehensive network of seismic reflection surveys from seven different campaigns over the past 50 years, as well as point constraints across an array of petroleum industry drill holes. Particular attention is given to a detailed representation of Quaternary stratigraphy, representing shallow (z<250m) near-surface layers in the model. Seismic velocities are obtained from seismic reflection processing (for Vp) and also recently performed active and passive surface wave analyses (for Vs). Over 1,700 water wells in the region are used to constrain the complex inter-bedded Quaternary stratigraphy (gravels, sands, silts, organics etc.) near the coastline, including beneath urban Christchurch, which has resulted from fluvial deposition and marine regression and transgression. For the near-surface Springston and Christchurch Formations in the Christchurch urban area (z<50m), high-spatial resolution seismic velocities (including Vs30 ) were obtained from over 13,000 cone penetration tests combined with a recently developed CPT-Vs correlation.
This paper presents an examination of ground motion observations from 20 near-source strong motion stations during the most significant 10 events in the 2010-2011 Canterbury earthquake to examine region-specific systematic effects based on relaxing the conventional ergodic assumption. On the basis of similar site-to-site residuals, surfical geology, and geographical proximity, 15 of the 20 stations are grouped into four sub-regions: the Central Business District; and Western, Eastern, and Northern suburbs. Mean site-to-site residuals for these sub-regions then allows for the possibility of non-ergodic ground motion prediction over these sub-regions of Canterbury, rather than only at strong motion station locations. The ratio of the total non-ergodic vs. ergodic standard deviation is found to be, on average, consistent with previous studies, however it is emphasized that on a site-by-site basis the non-ergodic standard deviation can easily vary by ±20%.
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.
Recurrent liquefaction in Christchurch during the 2010-2011 Canterbury earthquake sequence created a wealth of shallow subsurface intrusions with geometries and orientations governed by (1) strong ground motion severity and duration, and (2) intrinsic site characteristics including liquefaction susceptibility, lateral spreading severity, geomorphic setting, host sediment heterogeneity, and anthropogenic soil modifications. We present a suite of case studies that demonstrate how each of these characteristics influenced the geologic expressions of contemporary liquefaction in the shallow subsurface. We compare contemporary features with paleo-features to show how geologic investigations of recurrent liquefaction can provide novel insights into the shaking characteristics of modern and paleo-earthquakes, the influence of geomorphology on liquefaction vulnerability, and the possible controls of anthropogenic activity on the geologic record. We conclude that (a) sites of paleo-liquefaction in the last 1000-2000 years corresponded with most severe liquefaction during the Canterbury earthquake sequence, (b) less vulnerable sites that only liquefied in the strongest and most proximal contemporary earthquakes are unlikely to have liquefied in the last 1000-2000 years or more, (c) proximal strong earthquakes with large vertical accelerations favoured sill formation at some locations, (d) contemporary liquefaction was more severe than paleoliquefaction at all study sites, and (e) stratigraphic records of successive dike formation were more complete at sites with severe lateral spreading, (f) anthropogenic fill suppressed surface liquefaction features and altered subsurface liquefaction architecture.
The 2010-2011 Canterbury earthquakes were recorded over a dense strong motion network in the near-source region, yielding significant observational evidence of seismic complexities, and a basis for interpretation of multi-disciplinary datasets and induced damage to the natural and built environment. This paper provides an overview of observed strong motions from these events and retrospective comparisons with both empirical and physics-based ground motion models. Both empirical and physics-based methods provide good predictions of observations at short vibration periods in an average sense. However, observed ground motion amplitudes at specific locations, such as Heathcote Valley, are seen to systematically depart from ‘average’ empirical predictions as a result of near surface stratigraphic and topographic features which are well modelled via sitespecific response analyses. Significant insight into the long period bias in empirical predictions is obtained from the use of hybrid broadband ground motion simulation. The comparison of both empirical and physics-based simulations against a set of 10 events in the sequence clearly illustrates the potential for simulations to improve ground motion and site response prediction, both at present, and further in the future.
None
Following the devastating 1931 Hawke's Bay earthquake, buildings in Napier and surrounding areas in the Hawke's Bay region were rebuilt in a comparatively homogenous structural and architectural style comprising the region's famous Art Deco stock. These interwar buildings are most often composed of reinforced concrete two-way space frames, and although they have comparatively ductile detailing for their date of construction, are often expected to be brittle, earthquake-prone buildings in preliminary seismic assessments. Furthermore, the likelihood of global collapse of an RC building during a design-level earthquake became an issue warranting particular attention following the collapse of multiple RC buildings in the February 22, 2011 Christchurch earthquake. Those who value the architectural heritage and future use of these iconic Art Deco buildings - including building owners, tenants, and city officials, among others - must consider how they can be best preserved and utilized functionally given the especially pressing implications of relevant safety, regulatory, and economic factors. This study was intended to provide information on the seismic hazard, geometric weaknesses, collapse hazards, material properties, structural detailing, empirically based vulnerability, and recommended analysis approaches particular to Art Deco buildings in Hawke's Bay as a resource for professional structural engineers tasked with seismic assessments and retrofit designs for these buildings. The observed satisfactory performance of similar low-rise, ostensibly brittle RC buildings in other earthquakes and the examination of the structural redundancy and expected column drift capacities in these buildings, led to the conclusion that the seismic capacity of these buildings is generally underrated in simple, force-based assessments.
The Canterbury Earthquake Sequence (CES) of 2010-2011 produced large seismic moments up to Mw 7.1. These large, near-to-surface (<15 km) ruptures triggered >6,000 rockfall boulders on the Port Hills of Christchurch, many of which impacted houses and affected the livelihoods of people within the impacted area. From these disastrous and unpredicted natural events a need arose to be able to assess the areas affected by rockfall events in the future, where it is known that a rockfall is possible from a specific source outcrop but the potential boulder runout and dynamics are not understood. The distribution of rockfall deposits is largely constrained by the physical properties and processes of the boulder and its motion such as block density, shape and size, block velocity, bounce height, impact and rebound angle, as well as the properties of the substrate. Numerical rockfall models go some way to accounting for all the complex factors in an algorithm, commonly parameterised in a user interface where site-specific effects can be calibrated. Calibration of these algorithms requires thorough field checks and often experimental practises. The purpose of this project, which began immediately following the most destructive rupture of the CES (February 22, 2011), is to collate data to characterise boulder falls, and to use this information, supplemented by a set of anthropogenic boulder fall data, to perform an in-depth calibration of the three-dimensional numerical rockfall model RAMMS::Rockfall. The thesis covers the following topics: • Use of field data to calibrate RAMMS. Boulder impact trails in the loess-colluvium soils at Rapaki Bay have been used to estimate ranges of boulder velocities and bounce heights. RAMMS results replicate field data closely; it is concluded that the model is appropriate for analysing the earthquake-triggered boulder trails at Rapaki Bay, and that it can be usefully applied to rockfall trajectory and hazard assessment at this and similar sites elsewhere. • Detailed analysis of dynamic rockfall processes, interpreted from recorded boulder rolling experiments, and compared to RAMMS simulated results at the same site. Recorded rotational and translational velocities of a particular boulder show that the boulder behaves logically and dynamically on impact with different substrate types. Simulations show that seasonal changes in soil moisture alter rockfall dynamics and runout predictions within RAMMS, and adjustments are made to the calibration to reflect this; suggesting that in hazard analysis a rockfall model should be calibrated to dry rather than wet soil conditions to anticipate the most serious outcome. • Verifying the model calibration for a separate site on the Port Hills. The results of the RAMMS simulations show the effectiveness of calibration against a real data set, as well as the effectiveness of vegetation as a rockfall barrier/retardant. The results of simulations are compared using hazard maps, where the maximum runouts match well the mapped CES fallen boulder maximum runouts. The results of the simulations in terms of frequency distribution of deposit locations on the slope are also compared with those of the CES data, using the shadow angle tool to apportion slope zones. These results also replicate real field data well. Results show that a maximum runout envelope can be mapped, as well as frequency distribution of deposited boulders for hazard (and thus risk) analysis purposes. The accuracy of the rockfall runout envelope and frequency distribution can be improved by comprehensive vegetation and substrate mapping. The topics above define the scope of the project, limiting the focus to rockfall processes on the Port Hills, and implications for model calibration for the wider scientific community. The results provide a useful rockfall analysis methodology with a defensible and replicable calibration process, that has the potential to be applied to other lithologies and substrates. Its applications include a method of analysis for the selection and positioning of rockfall countermeasure design; site safety assessment for scaling and demolition works; and risk analysis and land planning for future construction in Christchurch.
This article argues that active coordination of research engagement after disasters has the potential to maximize research opportunities, improve research quality, increase end-user engagement, and manage escalating research activity to mitigate ethical risks posed to impacted populations. The focus is on the coordination of research activity after the 22nd February 2011 Mw6.2 Christchurch earthquake by the then newly-formed national research consortium, the Natural Hazards Research Platform, which included a social science research moratorium during the declared state of national emergency. Decisions defining this organisation’s functional and structural parameters are analyzed to identify lessons concerning the need for systematic approaches to the management of post disaster research, in collaboration with the response effort. Other lessons include the importance of involving an existing, broadly-based research consortium, ensuring that this consortium's coordination role is fully integrated into emergency management structures, and ensuring that all aspects of decision-making processes are transparent and easily accessed.
The Sendai Framework for Disaster Risk Reduction 2015-2030 finds that, despite progress in disaster risk reduction over the last decade “evidence indicates that exposure of persons and assets in all countries has increased faster than vulnerability has decreased, thus generating new risk and a steady rise in disaster losses” (p.4, UNISDR 2015). Fostering cooperation among relevant stakeholders and policy makers to “facilitate a science-policy interface for effective decisionmaking in disaster risk management” is required to achieve two priority areas for action, understanding disaster risk and enhancing disaster preparedness (p. 13, p. 23, UNISDR 2015). In other topic areas, the term science-policy interface is used interchangeably with the term boundary organisation. Both terms are usually used refer to systematic collaborative arrangements used to manage the intersection, or boundary, between science and policy domains, with the aim of facilitating the joint construction of knowledge to inform decision-making. Informed by complexity theory, and a constructivist focus on the functions and processes that minimize inevitable tensions between domains, this conceptual framework has become well established in fields where large complex issues have significant economic and political consequences, including environmental management, biodiversity, sustainable development, climate change and public health. To date, however, there has been little application of this framework in the disaster risk reduction field. In this doctoral project the boundary management framework informs an analysis of the research response to the 2010-2011 Canterbury Earthquake Sequence, focusing on the coordination role of New Zealand’s national Natural Hazards Research Platform. The project has two aims. It uses this framework to tell the nuanced story of the way this research coordination role evolved in response to both the complexity of the unfolding post-disaster environment, and to national policy and research developments. Lessons are drawn from this analysis for those planning and implementing arrangements across the science-policy boundary to manage research support for disaster risk reduction decision-making, particularly after disasters. The second aim is to use this case study to test the utility of the boundary management framework in the disaster risk reduction context. This requires that terminology and concepts are explained and translated in terms that make this analysis as accessible as possible across the disciplines, domains and sectors involved in disaster risk reduction. Key findings are that the focus on balance, both within organisations, and between organisations and domains, and the emphasis on systemic effects, patterns and trends, offer an effective and productive alternative to the more traditional focus on individual or organisational performance. Lessons are drawn concerning the application of this framework when planning and implementing boundary organisations in the hazard and disaster risk management context.