Research report N. 2016-2 This report presents the simplified seismic assessment of a case study reinforced concrete (RC) building following the newly developed and refined NZSEE/MBIE guidelines on seismic assessment (NZSEE/MBIE, semi-final draft 26 October 2016). After an overview of the step-by-step ‘diagnostic’ process, including an holistic and qualitative description of the expected vulnerabilities and of the assessment strategy/methodology, focus is given, whilst not limited, to the implementation of a Detailed Seismic Assessment (DSA) (NZSEE/MBIE, 2016c). The DSA is intended to provide a more reliable and consistent outcome than what can be provided by an initial seismic assessment (ISA). In fact, while the Initial Seismic Assessment (ISA), of which the Initial Evaluation Procedure is only a part of, is the more natural and still recommended first step in the overall assessment process, it is mostly intended to be a coarse evaluation involving as few resources as reasonably possible. It is thus expected that an ISA will be followed by a Detailed Seismic Assessment (DSA) not only where the threshold of 33%NBS is not achieved but also where important decisions are intended that are reliant on the seismic status of the building. The use of %NBS (% New Building Standard) as a capacity/demand ratio to describe the result of the seismic assessment at all levels of assessment procedure (ISA through to DSA) is deliberate by the NZSEE/MBIE guidelines (Part A) (NZSEE/MBIE 2016a). The rating for the building needs only be based on the lowest level of assessment that is warranted for the particular circumstances. Discussion on how the %NBS rating is to be determined can be found in Section A3.3 (NZSEE/MBIE 2016a), and, more specifically, in Part B for the ISA (NZSEE/MBIE 2016b) and Part C for the DSA (NZSEE/MBIE 2016c). As per other international approaches, the DSA can be based on several analysis procedures to assess the structural behaviour (linear, nonlinear, static or dynamic, force or displacement-based). The significantly revamped NZSEE 2016 Seismic Assessment Guidelines strongly recommend the use of an analytical (basically ‘by hand’) method, referred to the Simple Lateral Mechanism Analysis (SLaMA) as a first phase of any other numerically-based analysis method. Significant effort has thus been dedicated to provide within the NZSEE 2016 guidelines (NZSEE/MBIE 2016c) a step-by-step description of the procedure, either in general terms (Chapter 2) or with specific reference to Reinforced Concrete Buildings (Chapter 5). More specifically, extract from the guidelines, NZSEE “recommend using the Simple Lateral Mechanism Analysis (SLaMA) procedure as a first step in any assessment. While SLaMA is essentially an analysis technique, it enables assessors to investigate (and present in a simple form) the potential contribution and interaction of a number of structural elements and their likely effect on the building’s global capacity. In some cases, the results of a SLaMA will only be indicative. However, it is expected that its use should help assessors achieve a more reliable outcome than if they only carried out a detailed analysis, especially if that analysis is limited to the elastic range For complex structural systems, a 3D dynamic analysis may be necessary to supplement the simplified nonlinear Simple Lateral Mechanism Analysis (SLaMA).” This report presents the development of a full design example for the the implementation of the SLaMA method on a case study buildings and a validation/comparison with a non-linear static (pushover) analysis. The step-by-step-procedure, summarized in Figure 1, will be herein demonstrated from a component level (beams, columns, wall elements) to a subassembly level (hierarchy of strength in a beam-column joint) and to a system level (frame, C-Wall) assuming initially a 2D behaviour of the key structural system, and then incorporating a by-hand 3D behaviour (torsional effects).
This poster presents work to date on ground motion simulation validation and inversion for the Canterbury, New Zealand region. Recent developments have focused on the collection of different earthquake sources and the verification of the SPECFEM3D software package in forward and inverse simulations. SPECFEM3D is an open source software package which simulates seismic wave propagation and performs adjoint tomography based upon the spectral-element method. Figure 2: Fence diagrams of shear wave velocities highlighting the salient features of the (a) 1D Canterbury velocity model, and (b) 3D Canterbury velocity model. Figure 5: Seismic sources and strong motion stations in the South Island of New Zealand, and corresponding ray paths of observed ground motions. Figure 3: Domain used for the 19th October 2010 Mw 4.8 case study event including the location of the seismic source and strong motion stations. By understanding the predictive and inversion capabilities of SPECFEM3D, the current 3D Canterbury Velocity Model can be iteratively improved to better predict the observed ground motions. This is achieved by minimizing the misfit between observed and simulated ground motions using the built-in optimization algorithm. Figure 1 shows the Canterbury Velocity Model domain considered including the locations of small-to-moderate Mw events [3-4.5], strong motion stations, and ray paths of observed ground motions. The area covered by the ray paths essentially indicates the area of the model which will be most affected by the waveform inversion. The seismic sources used in the ground motion simulations are centroid moment tensor solutions obtained from GeoNet. All earthquake ruptures are modelled as point sources with a Gaussian source time function. The minimum Mw limit is enforced to ensure good signal-to-noise ratio and well constrained source parameters. The maximum Mw limit is enforced to ensure the point source approximation is valid and to minimize off-fault nonlinear effects.
On November 14 2016 a magnitude 7.8 earthquake struck the south island of New Zealand. The earthquake lasted for just two minutes with severe seismic shaking and damage in the Hurunui and Kaikōura districts. Although these are predominantly rural areas, with scattered small towns and mountainous topography, they also contain road and rail routes that are essential parts of the national transport infrastructure. This earthquake and the subsequent recovery are of particular significance as they represent a disaster following in close proximity to another similar disaster, with the Canterbury earthquakes occurring in a neighboring district five years earlier. The research used an inductive qualitative case study to explore the nature of the Kaikōura recovery. That recovery process involved a complex interplay between the three parties; (a) the existing local government in the district, (b) central government agencies funding the recovery of the local residents and the national transport infrastructure, and (c) recovery leaders arriving with recent expertise from the earlier Canterbury disaster. It was evident that three groups: locals, government, and experts represented a multi-party governance debate in which the control of the Kaikōura earthquake recovery was shared amongst them. Each party had their own expertise, adgenda and networks that they brought to the Kaikōura recovery, but this created tensions between external expertise and local, community leadership. Recent earthquake research suggests that New Zealand is currently in the midst of an earthquake cluster, with further seismic disasters likely to occur in relatively close succession. This is likely to be compounded by the increasing frequency of other natural disasters with the effects of climate change. The present study investigates a phenomenon that may become increasingly common, with the transfer of disaster expertise from one event to another, and the interface between those experts with local and national government in directing recoveries. The findings of this study have implications for practitioners and policy makers in NZ and other countries where disasters are experienced in close spatial and temporal proximity.
- The Avon-Ōtākaro Redzone is an 11 kilometer stretch of land along the Avon-Ōtākaro River in Christchurch. - This project focused on the creation of a publicly available biodiversity map of the AvonŌtākaro River Corridor, a project undertaken as part of the ecological restoration of the Christchurch redzone. - This project originated from the Christchurch 2010-2011 earthquake sequence which saw liquefaction damage along 11km of the Avon River. Under guidance from The Nature Lab & Ōtākaro Living Laboratory, and various other experts, the primary research objective was to map historical biodiversity, identify hotspots, and assess areas for potential revegetation. - The data collected came from historical black maps, current iNaturalist data, and soil classification information. - The findings show that, pre-colonialism, the area was composed of herbaceous areas, wetlands, native shrubland, and tussock land, with key plants such as river, fern, tutu, and cabbage trees. - The post-earthquake analysis shows a transition from a residential area to patchy grasslands and swampy areas. - The findings also showed a strong relationship between historic sites and soil classifications, providing knowledge for past and future vegetation patterns and spread. - This map will be a valuable resource for conservation efforts and public engagement as the area transitions into a blue-green corridor.
The initial goal of this research was to explore how SME business models change in response to a crisis. Keeping this in mind, the business model canvas (Osterwalder & Pigneur, 2010) was used as a tool to analyse SME business models in the Canterbury region of New Zealand. The purpose was to evaluate the changes SMEs instituted in their business models after being hit by a series of earthquakes in 2010 and 2011. The idea was to conduct interviews with business owners and analyse them using grounded theory methods. As this method is iterative and requires simultaneous data collection and analysis, a tentative model was proposed after first phase of the data collection and analysis. However, as a result of this process, it became apparent that owner-specific characteristics, action orientation and networks were more prominent in the data than business model elements. Although the SMEs in this study experienced several operational changes in their business models, such as a change of location, modifications to their payment terms or expanded/restricted target markets, the suggested framework highlights how owner-specific attributes ensured the recovery of their businesses. After the initial framework was suggested, subsequent interviews were conducted to test, verify, and modify the tentative model. Three aspects of business recovery emerged: (a) cognitive coping – the business owner’s mind-set and motive; (b) adaptive coping – the ability of business owner to take corrective actions; and (c) social capital – the social network of a business owner, including formal and informal connections and their significance. Three distinct groups were identified; self-sufficient SMEs, socially-based SMEs and surviving SMEs. This thesis proposes a grounded theory of business recovery for SMEs following a disaster. Cognitive coping and social capital enabled the owners to take actions, which eventually led to the desired outcomes for the businesses.
One of the current challenges in physics-based ground-motion simulations is to refine the modeling of local site effects. These effects require a finer spatial resolution in the material modeling than that generally considered in regional-scale simulations. Because of this, empirical amplification factors are typically applied to capture these unmodeled phenomena. The ergodic nature of this approach suggests that there is room for improvement. In this study, the predictive capability of simulations is evaluated using alternative methods for capturing local site effects. In addition to the conventional empirical approach, two methods are examined that allow for more site-specific information to be incorporated: the square-root impedance method and the 1D time domain site-response analysis. The three approaches are tested using 1000+ observed ground motions from 150+ small-magnitude events (3.5 ≤ Mw ≤ 5.0), recorded at 20 strong-motion stationsin the Canterbury, New Zealand, region. These 20 well-characterized sites represent a wide range of soil conditions, including stiff gravels with Vs30 values greater than 500 m/s, and sand and silt deposits with Vs30 valuesless than 200 m/s. Multiple intensity measures are computed and prediction residuals are partitioned using mixed-effects regression to rigorously assess the relative performance of the different approaches considered. The results indicate that the benefit of using more sophisticated methods is highly dependent on the characteristics of the site. Key site parameters and trends are identified and discussed in light of the assumptions and limitations of each approach.
This study followed two similarly affected, but socio-economically disparate suburbs as residents responded to and attempted to recover from the devastating 6.3 magnitude earthquake that struck Christchurch, New Zealand, on February 22, 2011. More specifically, it focuses on the role of local churches, community-based organisations (CBOs) and non-governmental organisations (NGOs), here referred to broadly as civil society, in meeting the immediate needs of local residents and assisting with the longer-term recovery of each neighbourhood. Despite considerable socioeconomic differences between the two neighbourhoods, civil society in both suburbs has been vital in addressing the needs of locals in the short and long term following the earthquake. Institutions were able to utilise local knowledge of both residents and the extent of damage in the area to a) provide a swifter local response than government or civil defence and then help direct the relief these agencies provided locally; b) set up central points for distribution of supplies and information where locals would naturally gather; c) take action on what were perceived to be unmet needs; and d) act as a way of bridging locals to a variety of material, informational, and emotional resources. However the findings also support literature which indicates that other factors are also important in understanding neighbourhood recovery and the role of civil society, including: local leadership; a shared, place-based identity; the type and form of civil society organizations; social capital; and neighbourhood- and household-level indicators of relative vulnerability and inequality. The intertwining of these various factors seems to influence how these neighbourhoods have coped with and taken steps in recovering from this disaster. It is recommended that future research be directed at developing a better understanding of how this occurs. It is suggested that a model similar to Yasui’s (2007) Community Vulnerability and Capacity model be developed as a useful way to approach future research in this area.
The objective of this study is to examine the influence of near-fault motions on liquefaction triggering in Christchurch and neighboring towns during the 2010-2011 Canterbury earthquake sequence (CES). The CES began with the 4 September 2010, Mw7.1 Darfield earthquake and included up to ten events that triggered liquefaction. However, most notably, widespread liquefaction was induced by the Darfield earthquake and the Mw6.2, 22 February 2011 Christchurch earthquake. Of particular relevance to this study is the forward directivity effects that were prevalent in the motions recorded during the Darfield earthquake, and to a much lesser extent, during the Christchurch earthquake. A 2D variant of the Richart-Newmark fatigue theory was used to compute the equivalent number of cycles (neq) for the ground motions, where volumetric strain was used as the damage metric. This study is unique because it considers the contribution and phasing of both the fault-normal and fault-parallel components of motion on neq and the magnitude scaling factor (MSF). It was found that when the fault-normal and fault-parallel motions were treated individually, the former yielded a lower neq than the latter. Additionally, when the combined effects of fault-normal and fault-parallel components were considered, it was found that the MSF were higher than those commonly used. This implies that motions containing near-fault effects are less demanding on the soil than motions that do not. This may be one of several factors that resulted in less severe liquefaction occurring during the Darfield earthquake than the Christchurch earthquake.
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.
A team of earthquake geologists, seismologists and engineering seismologists from GNS Science, NIWA, University of Canterbury, and Victoria University of Wellington have collectively produced an update of the 2002 national probabilistic seismic hazard (PSH) model for New Zealand. The new model incorporates over 200 new onshore and offshore fault sources, and utilises newly developed New Zealand-based scaling relationships and methods for the parameterisation of the fault and subduction interface sources. The background seismicity model has also been updated to include new seismicity data, a new seismicity regionalisation, and improved methodology for calculation of the seismicity parameters. Background seismicity models allow for the occurrence of earthquakes away from the known fault sources, and are typically modelled as a grid of earthquake sources with rate parameters assigned from the historical seismicity catalogue. The Greendale Fault, which ruptured during the M7.1, 4 September 2010 Darfield earthquake, was unknown prior to the earthquake. However, the earthquake was to some extent accounted for in the PSH model. The maximum magnitude assumed in the background seismicity model for the area of the earthquake is 7.2 (larger than the Darfield event), but the location and geometry of the fault are not represented. Deaggregations of the PSH model for Christchurch at return periods of 500 years and above show that M7-7.5 fault and background source-derived earthquakes at distances less than 40 km are important contributors to the hazard. Therefore, earthquakes similar to the Darfield event feature prominently in the PSH model, even though the Greendale Fault was not an explicit model input.
Interagency Emergency Response Teams (IERTs) play acrucial role in times of disasters. Therefore it is crucial to understand more thoroughly the communication roles and responsibilities of interagency team members and to examine how individual members communicate within a complex, evolving, and unstable environment. It is also important to understand how different organisational identities and their spatial geographies contribute to the interactional dynamics. Earthquakes hit the Canterbury region on September, 2010 and then on February 2011 a more devastating shallow earthquake struck resulting in severe damage to the Aged Residential Care (ARC) sector. Over 600 ARC beds were lost and 500 elderly and disabled people were displaced. Canterbury District Health Board (CDHB) set up an interagency emergency response team to address the issues of vulnerable people with significant health and disability needs who were unable to access their normal supports due to the effects of the earthquake. The purpose of this qualitative interpretive study is to focus on the case study of the response and evacuation of vulnerable people by interagencies responding to the event. Staff within these agencies were interviewed with a focus on the critical incidents that either stabilised or negatively influenced the outcome of the response. The findings included the complexity of navigating multiple agencies communication channels; understanding the different hierarchies and communication methods within each agency; data communication challenges when infrastructures were severely damaged; the importance of having the right skills, personal attributes and understanding of the organisations in the response; and the significance of having a liaison in situ representing and communicating through to agencies geographically dispersed from Canterbury. It is hoped that this research will assist in determining a future framework for interagency communication best practice and policy.
Following the 2010-2011 Canterbury (New Zealand) earthquake sequence, lightly reinforced wall structures in the Christchurch central business district were observed to form undesirable crack patterns in the plastic hinge region, while yield penetration either side of cracks and into development zones was less than predicted using empirical expressions. To some extent this structural behaviour was unexpected and has therefore demonstrated that there may be less confidence in the seismic performance of conventionally designed reinforced concrete (RC) structures than previously anticipated. This paper provides an observation-based comparison between the behaviour of RC structural components in laboratory testing and the unexpected structural behaviour of some case study buildings in Christchurch that formed concentrated inelastic deformations. The unexpected behaviour and poor overall seismic performance of ‘real’ buildings (compared to the behaviour of laboratory test specimens) was due to the localization of peak inelastic strains, which in some cases has arguably led to: (i) significantly less ductility capacity; (ii) less hysteretic energy dissipation; and (iii) the fracture of the longitudinal reinforcement. These observations have raised concerns about whether lightly reinforced wall structures can satisfy the performance objective of “Life Safety” at the Ultimate Limit State. The significance of these issues and potential consequences has prompted a review of potential problems with the testing conditions and procedures that are commonly used in seismic experimentations on RC structures. This paper attempts to revisit the principles of RC mechanics, in particular, the influence of loading history, concrete tensile strength, and the quantity of longitudinal reinforcement on the performance of real RC structures. Consideration of these issues in future research on the seismic performance of RC might improve the current confidence levels in newly designed conventional RC structures.
The Porter's Pass-Amberley Fault Zone (PPAFZ) is a complex zone of anastomosing faults and folds bounding the south-eastern edge of the transition from subducting Pacific Plate to continental collision on the Australia Plate boundary. This study combines mapping of a 2000 km2 zone from the Southern Alps northeast to the coast near Amberley, 40 km north of metropolitan Christchurch, with an analysis of seismicity and a revision of regional seismic hazard. Three structural styles: 1) a western strike-slip, and 2) a more easterly thrust and reverse domain, pass into 3) a northwest verging fold belt on the northern Canterbury Plains, reflecting the structural levels exposed and the evolving west to east propagation. Basal remnants of a Late Cretaceous-Cenozoic, largely marine sedimentary cover sequence are preserved as outliers that unconformably overlie Mesozoic basement (greywacke and argillite of the Torlesse terrain) in the mountains of the PPAFZ and are underlain by a deeply leached zone which is widely preserved. Structure contouring of the unconformity surface indicates maximum, differential uplift of c.2600 m in the southwest, decreasing to c.1200 m in the coastal fold belt to the northeast. Much lower rates (or reversal) of uplift are evident a few kilometres southeast of the PPAFZ range-front escarpment. The youngest elements of the cover sequence are basement-derived conglomerates of Plio-Pleistocene age preserved on the SE margin. The source is more distant than the intervening mountains of the PPAFZ, probably from the Southern Alps, to the west and northwest. The absence of another regional unconformity on Mesozoic basement, older than Pleistocene, indicates that this uplift is post-Pliocene. Late Pleistocene(<100 kyr) differential uplift rates of c.0.5-2.7 m/kyr from uplifted marine terraces at the east coast, and rates of 2.5-3.3 m/kyr for tectonically-induced river-down cutting further west, suggest that uplift commenced locally during the last 1 Ma, and possibly within the last 0.5 Ma, if average rates are assumed to be uniform over time. Analysis of seismicity, recorded during a 10 week regional survey of micro earthquakes in 1990, identified two seismic zones beneath North Canterbury: 1) a sub-horizontal zone of activity restricted to the upper crust (≤12 km); and 2) a seismic zone in the lower crust (below a ceiling of ≤17 km), that broadens vertically to the north and northwest to a depth of c.40 km, with a bottom edge which dips 10°N and 15°NW, respectively. No events were recorded at depths between 12 km and 17 km, which is interpreted as a relatively aseismic, mid-crustal ductile layer. Marked differences (up to 60°) in the trend of strain axes for events above and below the inferred ductile layer are observed only north of the PPAFZ. A fundamental, north-to-south increase in the Wave-length of major geological structures occurs across the PPAFZ, and is interpreted as evidence that the upper crust beneath the Canterbury Plains is coupled to the lower crust, whereas the upper crust further north is not. Most of the recorded micro earthquakes <12 km deep beneath the PPAFZ have strike-slip mechanisms. It is probable that faults splay upward into the thrusts and folds at the surface as an evolving transpression zone in response to deep shear in basement. There have been no historic surface ruptures of the PPAFZ, but the zone has been characterised historically by frequent small earthquakes. Paleoseismic data (dated landslides and surface ruptures) compiled in this study, indicate a return period of 1500-1900 years between the last two M>7-7.5 earthquakes, and 500-700 years have elapsed since the last. The magnitudes of these events are estimated at c.M7.5, which represents a probable maximum magnitude for the PPAFZ. There are insufficient data to determine whether or not the frequency of large earthquakes conforms to a recognised model of behaviour, but comparison of the paleoseismic data with the historic record of smaller earthquakes, suggests that the magnitudes of the largest earthquakes in this zone are not exponentially distributed. A seismicity model for the PPAFZ (Elder et al., 1991) is reviewed, and a b-value of 1.0 is found to be consistent with the newly acquired paleoseismic data. This b-value reduces the predicted frequency of large earthquakes (M≥7.0) in this zone by a factor of 3.5, while retaining a conservative margin that allows for temporal variations in the frequency of large events and the possibility that the geological database is incomplete, suggesting grounds for revising the hazard model for Christchurch.
The 2015 New Zealand strong-motion database provides a wealth of new strong motion data for engineering applications. An important component of this database is the compilation of new site metadata, describing the soil conditions and site response at GeoNet strong motion stations. We have assessed and compiled four key site parameters for the ~460 GeoNet stations that recorded significant historical ground motions. Parameters include: site classification (NZS1170.5), Vs30, fundamental site period (Tsite) and depth to bedrock (Z1.0, i.e. depth to material with Vs > 1000 m/s). In addition, we have assigned a quality estimate (Quality 1 – 3) to these parameters to provide a qualitative estimate of the uncertainty. New highquality Tsite estimates have largely been obtained from newly available HVSR amplification curves and spectral ratios from inversion of regional strong motion data that has been reconciled with available geological information. Good quality Vs30 estimates, typically in urban centres, have also been incorporated following recent studies. Where site-specific measurements of Vs30 are not available, Vs30 is estimated based on surface geology following national Vs30 maps. New Z1.0 values have been provided from 3D subsurface models for Canterbury and Wellington. This database will be used in efforts to guide development and testing of new and existing ground motion prediction models in New Zealand. In particular, it will allow reexamination of the most important site parameters that control and predict site response in a New Zealand setting. Furthermore, it can be used to provide information about suitable rock reference sites for seismological research, and as a guide to site-specific references in the literature. We discuss compilation of the database, preliminary insights so far, and future directions.
Recent experiences from the Darfield and Canterbury, New Zealand earthquakes have shown that the soft soil condition of saturated liquefiable sand has a profound effect on seismic response of buildings, bridges and other lifeline infrastructure. For detailed evaluation of seismic response three dimensional integrated analysis comprising structure, foundation and soil is required; such an integrated analysis is referred to as Soil Foundation Structure Interaction (SFSI) in literatures. SFSI is a three-dimensional problem because of three primary reasons: first, foundation systems are three-dimensional in form and geometry; second, ground motions are three-dimensional, producing complex multiaxial stresses in soils, foundations and structure; and third, soils in particular are sensitive to complex stress because of heterogeneity of soils leading to a highly anisotropic constitutive behaviour. In literatures the majority of seismic response analyses are limited to plane strain configuration because of lack of adequate constitutive models both for soils and structures, and computational limitation. Such two-dimensional analyses do not represent a complete view of the problem for the three reasons noted above. In this context, the present research aims to develop a three-dimensional mathematical formulation of an existing plane-strain elasto-plastic constitutive model of sand developed by Cubrinovski and Ishihara (1998b). This model has been specially formulated to simulate liquefaction behaviour of sand under ground motion induced earthquake loading, and has been well-validated and widely implemented in verifcation of shake table and centrifuge tests, as well as conventional ground response analysis and evaluation of case histories. The approach adopted herein is based entirely on the mathematical theory of plasticity and utilises some unique features of the bounding surface plasticity formalised by Dafalias (1986). The principal constitutive parameters, equations, assumptions and empiricism of the existing plane-strain model are adopted in their exact form in the three-dimensional version. Therefore, the original two-dimensional model can be considered as a true subset of the three-dimensional form; the original model can be retrieved when the tensorial quantities of the three dimensional version are reduced to that of the plane-strain configuration. Anisotropic Drucker-Prager type failure surface has been adopted for the three-dimensional version to accommodate triaxial stress path. Accordingly, a new mixed hardening rule based on Mroz’s approach of homogeneous surfaces (Mroz, 1967) has been introduced for the virgin loading surface. The three-dimensional version is validated against experimental data for cyclic torsional and triaxial stress paths.
Background: We are in a period of history where natural disasters are increasing in both frequency and severity. They are having widespread impacts on communities, especially on vulnerable communities, those most affected who have the least ability to prepare or respond to a disaster. The ability to assemble and effectively manage Interagency Emergency Response Teams (IERTs) is critical to navigating the complexity and chaos found immediately following disasters. These teams play a crucial role in the multi-sectoral, multi-agency, multi-disciplinary, and inter-organisational response and are vital to ensuring the safety and well-being of vulnerable populations such as the young, aged, and socially and medically disadvantaged in disasters. Communication is key to the smooth operation of these teams. Most studies of the communication in IERTs during a disaster have been focussed at a macro-level of examining larger scale patterns and trends within organisations. Rarely found are micro-level analyses of interpersonal communication at the critical interfaces between collaborating agencies. This study set out to understand the experiences of those working at the interagency interfaces in an IERT set up by the Canterbury District Health Board to respond to the needs of the vulnerable people in the aftermath of the destructive earthquakes that hit Canterbury, New Zealand, in 2010-11. The aim of the study was to gain insights about the complexities of interpersonal communication (micro-level) involved in interagency response coordination and to generate an improved understanding into what stabilises the interagency communication interfaces between those agencies responding to a major disaster. Methods: A qualitative case study research design was employed to investigate how interagency communication interfaces were stabilised at the micro-level (“the case”) in the aftermath of the destructive earthquakes that hit Canterbury in 2010-11 (“the context”). Participant recruitment was undertaken by mapping which agencies were involved within the IERT and approaching representatives from each of these agencies. Data was collected via individual interviews using a semi-structured interview guide and was based on the “Critical Incident Technique”. Subsequently, data was transcribed verbatim and subjected to inductive analysis. This was underpinned theoretically by Weick’s “Interpretive Approach” and supported by Nvivo qualitative data analysis software. Results: 19 participants were interviewed in this study. Out of the inductive analysis emerged two primary themes, each with several sub-factors. The first major theme was destabilising/disruptive factors of interagency communication with five sub-factors, a) conflicting role mandates, b) rigid command structures, c) disruption of established communication structures, d) lack of shared language and understanding, and e) situational awareness disruption. The second major theme stabilising/steadying factors in interagency communication had four sub-factors, a) the establishment of the IERT, b) emergent novel communication strategies, c) establishment of a liaison role and d) pre-existing networks and relationships. Finally, there was a third sub-level identified during inductive analysis, where sub-factors from both primary themes were noted to be uniquely interconnected by emergent “consequences” arising out of the disaster context. Finally, findings were synthesised into a conceptual “Model of Interagency Communication at the Micro-level” based on this case study of the Canterbury earthquake disaster response. Discussion: The three key dimensions of The People, The Connections and The Improvisations served as a framework for the discussion of what stabilises interagency communication interfaces in a major disaster. The People were key to stabilising the interagency interfaces through functioning as a flexible conduit, guiding and navigating communication at the interagency interfaces and improving situational awareness. The Connections provided the collective competence, shared decision-making and prior established relationships that stabilised the micro-level communication at interagency interfaces. And finally, The Improvisations i.e., novel ideas and inventiveness that emerge out of rapidly changing post-disaster environments, also contributed to stabilisation of micro-level communication flows across interagency interfaces in the disaster response. “Command and control” hierarchical structures do provide clear processes and structures for teams working in disasters to follow. However, improvisations and novel solutions are also needed and often emerge from first responders (who are best placed to assess the evolving needs in a disaster where there is a high degree of uncertainty). Conclusion: This study highlights the value of incorporating an interface perspective into any study that seeks to understand the processes of IERTs during disaster responses. It also strengthens the requirement for disaster management frameworks to formally plan for and to allow for the adaptive responsiveness of local teams on the ground, and legitimise and recognise the improvisations of those in the role of emergent boundary spanners in a disaster response. This needs to be in addition to existing formal disaster response mechanisms. This study provides a new conceptual model that can be used to guide future case studies exploring stability at the interfaces of other IERTs and highlights the centrality of communication in the experiences of members of teams in the aftermath of a disaster. Utilising these new perspectives on stabilising communication at the interagency interfaces in disaster responses will have practical implications in the future to better serve the needs of vulnerable people who are at greatest risk of adverse outcomes in a disaster.
The 22 February 2011, Mw6.2-6.3 Christchurch earthquake is the most costly earthquake to affect New Zealand, causing 181 fatalities and severely damaging thousands of residential and commercial buildings, and most of the city lifelines and infrastructure. This manuscript presents an overview of observed geotechnical aspects of this earthquake as well as some of the completed and on-going research investigations. A unique aspect, which is particularly emphasized, is the severity and spatial extent of liquefaction occurring in native soils. Overall, both the spatial extent and severity of liquefaction in the city was greater than in the preceding 4th September 2010 Darfield earthquake, including numerous areas that liquefied in both events. Liquefaction and lateral spreading, variable over both large and short spatial scales, affected commercial structures in the Central Business District (CBD) in a variety of ways including: total and differential settlements and tilting; punching settlements of structures with shallow foundations; differential movements of components of complex structures; and interaction of adjacent structures via common foundation soils. Liquefaction was most severe in residential areas located to the east of the CBD as a result of stronger ground shaking due to the proximity to the causative fault, a high water table approximately 1m from the surface, and soils with composition and states of high susceptibility and potential for liquefaction. Total and differential settlements, and lateral movements, due to liquefaction and lateral spreading is estimated to have severely compromised 15,000 residential structures, the majority of which otherwise sustained only minor to moderate damage directly due to inertial loading from ground shaking. Liquefaction also had a profound effect on lifelines and other infrastructure, particularly bridge structures, and underground services. Minor damage was also observed at flood stop banks to the north of the city, which were more severely impacted in the 4th September 2010 Darfield earthquake. Due to the large high-frequency ground motion in the Port hills numerous rock falls and landslides also occurred, resulting in several fatalities and rendering some residential areas uninhabitable.
The magnitude Mw7.8 ‘Kaikōura’ earthquake occurred shortly after midnight on 14 November 2016. This paper presents an overview of the geotechnical impacts on the South Island of New Zealand recorded during the postevent reconnaissance. Despite the large moment magnitude of this earthquake, relatively little liquefaction was observed across the South Island, with the only severe manifestation occurring in the young, loose alluvial deposits in the floodplains of the Wairau and Opaoa Rivers near Blenheim. The spatial extent and volume of liquefaction ejecta across South Island is significantly less than that observed in Christchurch during the 2010-2011 Canterbury Earthquake Sequence, and the impact of its occurrence to the built environment was largely negligible on account of the severe manifestations occurring away from the areas of major development. Large localised lateral displacements occurred in Kaikōura around Lyell Creek. The soft fine-grained material in the upper portions of the soil profile and the free face at the creek channel were responsible for the accumulation of displacement during the ground shaking. These movements had severely impacted the houses which were built close (within the zone of large displacement) to Lyell Creek. The wastewater treatment facility located just north of Kaikōura also suffered tears in the liners of the oxidation ponds and distortions in the aeration system due to ground movements. Ground failures on the Amuri and Emu Plains (within the Waiau Valley) were small considering the large peak accelerations (in excess of 1g) experienced in the area. Minor to moderate lateral spreading and ejecta was observed at some bridge crossings in the area. However, most of the structural damage sustained by the bridges was a result of the inertial loading, and the damage resulting from geotechnical issues were secondary.
Tsunami have the potential to cause significant disruptions to society, including damage to infrastructure, critical to the every-day operation of society. Effective risk management is required to reduce the potential tsunami impacts to them. Christchurch city, situated on the eastern coast of New Zealand’s South Island, is exposed to a number of far-field tsunami hazards. Although the tsunami hazard has been well identified for Christchurch city infrastructure, the likely impacts have not been well constrained. To support effective risk management a credible and realistic infrastructure impact model is required to inform risk management planning. The objectives of this thesis are to assess the impacts on Christchurch city infrastructure from a credible, hypothetical far-field tsunami scenario. To achieve this an impact assessment process is adopted, using tsunami hazard and exposure measures to determine asset vulnerability and subsequent impacts. However, the thesis identified a number of knowledge gaps in infrastructure vulnerability to tsunami. The thesis addresses this by using two approaches: a tsunami damage matrix; and the development of tsunami fragility functions. The tsunami damage matrix pools together tsunami impacts on infrastructure literature, and post-event field observations. It represents the most comprehensive ‘look-up’ resource for tsunami impacts to infrastructure to date. This damage matrix can inform the assessment of tsunami impacts on Christchurch city infrastructure by providing a measure of damage likelihood at various hazard intensities. A more robust approach to tsunami vulnerability of infrastructure are fragility functions, which are also developed in this thesis. These were based on post-event tsunami surveys of the 2011 ‘Tohoku’ earthquake tsunami in Japan. The fragility functions are limited to road and bridge infrastructure, but represent the highest resolution measure of vulnerability for the given assets. As well as providing a measure of damage likelihood for a given tsunami hazard intensity, these also indicate a level of asset damage. The impact assessment process, and synthesized vulnerability measures, are used to run tsunami impact models for Christchurch infrastructure to determine the probability of asset damage occurring and to determine if impact will reach or exceed a given damage state. The models suggest that infrastructure damage is likely to occur in areas exposed to tsunami inundation in this scenario, with significant damage identified for low elevation roads and bridges. The results are presented and discussed in the context of the risk management framework, with emphasis on using risk assessment to inform risk treatment, monitoring and review. In summary, this thesis A) advances tsunami vulnerability and impact assessment methodologies for infrastructure and B) provides a tsunami impact assessment framework for Christchurch city infrastructure which will inform infrastructure tsunami risk management for planners, emergency managers and lifelines groups.
To this extent, modern buildings generally demonstrated good resistance to collapse during the recent earthquakes in New Zealand. However, damage to non-structural elements (NSE) has been persistent during these events. NSEs include secondary systems or components attached to the floors, roofs, and walls of a building or industrial facility that are not explicitly designed to participate in the main vertical or lateral load-bearing mechanism of the structure. They play a major role in the operational and functional aspects of buildings and contribute a major portion of the building’s overall cost. Therefore, they are expected to accommodate the effects of seismic actions such as drifts and accelerations. Typical examples of NSEs include internal non-loadbearing partitions, suspended ceilings, sprinkler piping systems, architectural claddings, building contents, mechanical/electrical equipment, and furnishings. The main focus of this thesis is the drift sensitive NSEs: precast concrete cladding panels and internal partition walls. Even though most precast concrete cladding panels performed well from a life-safety point of view during recent earthquakes in NZ, some collapsed panels posed a significant threat to life safety. It is, therefore, important that the design and detailing of the panel-to-structure connections ensure that their strength and displacement capacity are adequate to meet the corresponding seismic demands, at least during design level earthquakes. In contrast, the partition wall is likely to get damaged and lose serviceability at a low inter-story drift unless designed to accommodate the relative deformations between them and the structure. Partition walls suffered wide-ranging damage such as screw failures, diagonal cracking, detachments to the gypsum linings, and anchorage failures during the 2011 Canterbury Earthquake Sequence in NZ. Therefore, the thesis is divided into two parts. Part I of the thesis focuses on developing novel low-damage precast concrete cladding panel connections, i.e. “rocking” connection details comprising vertically slotted steel embeds and weld plates. The low-damage seismic performance of novel “rocking” connection details is verified through experimental tests comprising uni-directional, bi-directional, and multi-storey scaled quasi-static cyclic tests. Comparison with the seismic performance of traditional panel connections reported in the literature demonstrated the system’s significantly improved seismic resilience. Furthermore, the finite element models of panel connections and sealants are developed in ABAQUS. The force-drift responses of the “rocking” panel system modelled in SAP2000 is compared with the experimental results to evaluate their accuracy and validity. Part II of the thesis focuses on a) understanding the seismic performance of traditional rigid timber-framed partition wall, b) development and verification of low-damage connections (i.e. “rocking” connection details comprising of dual-slot tracks), and c) seismic evaluation of partition walls with a novel “bracketed and slotted” connections (comprising of innovative fastener and plastic bracket named Flexibracket) under uni-directional and bidirectional quasi-static cyclic loadings. Moreover, parametric investigation of the partition walls was conducted through several experimental tests to understand better the pros and cons of the rocking connection details. The experimental results have confirmed that the implementation of the proposed low damage solutions of precast cladding panels and internal partition walls can significantly reduce their damage in a building.
This is an interim report from the research study performed within the NHRP Research Project “Impacts of soil liquefaction on land, buildings and buried pipe networks: geotechnical evaluation and design, Project 3: Seismic assessment and design of pipe networks in liquefiable soils”. The work presented herein is a continuation of the comprehensive study on the impacts of Christchurch earthquakes on the buried pipe networks presented in Cubrinovski et al. (2011). This report summarises the performance of Christchurch City’s potable water, waste water and road networks through the 2010-2011 Canterbury Earthquake Sequence (CES), and particularly focuses on the potable water network. It combines evidence based on comprehensive and well-documented data on the damage to the water network, detailed observations and interpretation of liquefaction-induced land damage, records and interpretations of ground motion characteristics induced by the Canterbury earthquakes, for a network analysis and pipeline performance evaluation using a GIS platform. The study addresses a range of issues relevant in the assessment of buried networks in areas affected by strong earthquakes and soil liquefaction. It discusses performance of different pipe materials (modern flexible pipelines and older brittle pipelines) including effects of pipe diameters, fittings and pipeline components/details, trench backfill characteristics, and severity of liquefaction. Detailed breakdown of key factors contributing to the damage to buried pipes is given with reference to the above and other relevant parameters. Particular attention is given to the interpretation, analysis and modelling of liquefaction effects on the damage and performance of the buried pipe networks. Clear link between liquefaction severity and damage rate for the pipeline has been observed with an increasing damage rate seen with increasing liquefaction severity. The approach taken here was to correlate the pipeline damage to LRI (Liquefaction Resistance Index, newly developed parameter in Cubrinovski et al., 2011) which represents a direct measure for the soil resistance to liquefaction while accounting for the seismic demand through PGA. Key quality of the adopted approach is that it provides a general methodology that in conjunction with conventional methods for liquefaction evaluation can be applied elsewhere in New Zealand and internationally. Preliminary correlations between pipeline damage (breaks km-1), liquefaction resistance (LRI) and seismic demand (PGA) have been developed for AC pipes, as an example. Such correlations can be directly used in the design and assessment of pipes in seismic areas both in liquefiable and non-liquefiable areas. Preliminary findings on the key factors for the damage to the potable water pipe network and established empirical correlations are presented including an overview of the damage to the waste water and road networks but with substantially less detail. A comprehensive summary of the damage data on the buried pipelines is given in a series of appendices.
Research Report No.2010-03 Ground motion prediction equations (GMPEs) for geometric-mean pseudo-spectral acceleration amplitudes from New Zealand (NZ) earthquakes are developed. A database of 2437 three-component ground motion records is developed by applying stringent quality criteria to the historically recorded events in NZ. Despite the large number of records, the database is deficient in empirical records from large magnitude events recorded at close distances to the fault rupture plane. As a result, the basis for the NZ-specific GMPE development is to examine the applicability of foreign GMPEs for similar tectonic regions and then modify the most applicable GMPEs based on both theoretical and statistically significant empirically-driven arguments. For active shallow crustal events, five different GMPEs are considered. It was found that the McVerry et al. (2006) model, which is the current model upon which seismic design guidelines and site-specific seismic hazard analyses in NZ are based, provided the worst fit to the NZ database, and that the Chiou et al. (2010) (C10) modification of the Chiou and Youngs (2008) model was the most applicable. Discrepancies between the C10 model and the NZ database that were empirically identified and theoretically justified were used to modify the C10 model for: (i) small magnitude scaling; (ii) scaling of short period ground motion from normal faulting events in volcanic crust; (iii) scaling of ground motions on very hard rock sites; (iv) anelastic attenuation in the NZ crust; and (v) consideration of the increased anelastic attenuation in the Taupo Volcanic Zone (TVZ). For subduction slab events, initially three models were considered. It was found that all of the models had some significant biases with respect to applicability for NZ. The Zhao et al. (2006) (Z06) model was selected because of the rigorous database upon which it was developed and modified by: (i) NZ-specific scaling at small magnitudes; (ii) path scaling at large distances; (iii) consideration of the increased TVZ attenuation; and (iv) revision of the standard deviation model. Based on these modifications the developed model showed no bias of the inter- and intra-event residuals as a function of various predictor variables. The standard deviation of the residuals using the revised standard deviation model also indicated that the model has an adequate precision. Three GMPEs were considered for subduction interface events. The Zhao et al. (2006) (Z06) model was the best performing model with only bias exhibited in the site response model, and possible over-prediction of large magnitude events. The Z06 interface model was modified to account for site response and magnitude scaling using the same functional forms as those of the developed active shallow crustal and subduction slab models. The developed model showed no bias of the inter- and intra-event residuals as a function of various predictor variables. The developed GMPEs include specific features as evident in the NZ database; consistent scaling for parameters not well constrained by the NZ database; and pseudo-spectral amplitudes for vibration periods from 0.01 to 10 seconds. Hence, these models represent a significant advance in the state-of-the art for empirical ground motion prediction in NZ.
In the period between September 2010 and December 2011, Christchurch (New Zealand) and its surroundings were hit by a series of strong earthquakes including six significant events, all generated by local faults in proximity to the city: 4 September 2010 (Mw=7.1), 22 February 2011 (Mw=6.2), 13 June 2011 (Mw=5.3 and Mw=6.0) and 23 December 2011 (M=5.8 and (M=5.9) earthquakes. As shown in Figure 1, the causative faults of the earthquakes were very close to or within the city boundaries thus generating very strong ground motions and causing tremendous damage throughout the city. Christchurch is shown as a lighter colour area, and its Central Business District (CBD) is marked with a white square area in the figure. Note that the sequence of earthquakes started to the west of the city and then propagated to the south, south-east and east of the city through a set of separate but apparently interacting faults. Because of their strength and proximity to the city, the earthquakes caused tremendous physical damage and impacts on the people, natural and built environments of Christchurch. The 22 February 2011 earthquake was particularly devastating. The ground motions generated by this earthquake were intense and in many parts of Christchurch substantially above the ground motions used to design the buildings in Christchurch. The earthquake caused 182 fatalities, collapse of two multi-storey reinforced concrete buildings, collapse or partial collapse of many unreinforced masonry structures including the historic Christchurch Cathedral. The Central Business District (CBD) of Christchurch, which is the central heart of the city just east of Hagley Park, was practically lost with majority of its 3,000 buildings being damaged beyond repair. Widespread liquefaction in the suburbs of Christchurch, as well as rock falls and slope/cliff instabilities in the Port Hills affected tens of thousands of residential buildings and properties, and shattered the lifelines and infrastructure over approximately one third of the city area. The total economic loss caused by the 2010-2011 Christchurch earthquakes is currently estimated to be in the range between 25 and 30 billion NZ dollars (or 15% to 18% of New Zealand’s GDP). After each major earthquake, comprehensive field investigations and inspections were conducted to document the liquefaction-induced land damage, lateral spreading displacements and their impacts on buildings and infrastructure. In addition, the ground motions produced by the earthquakes were recorded by approximately 15 strong motion stations within (close to) the city boundaries providing and impressive wealth of data, records and observations of the performance of ground and various types of structures during this unusual sequence of strong local earthquakes affecting a city. This paper discusses the liquefaction in residential areas and focuses on its impacts on dwellings (residential houses) and potable water system in the Christchurch suburbs. The ground conditions of Christchurch including the depositional history of soils, their composition, age and groundwater regime are first discussed. Detailed liquefaction maps illustrating the extent and severity of liquefaction across Christchurch triggered by the sequence of earthquakes including multiple episodes of severe re-liquefaction are next presented. Characteristic liquefaction-induced damage to residential houses is then described focussing on the performance of typical house foundations in areas affected by liquefaction. Liquefaction impacts on the potable water system of Christchurch is also briefly summarized including correlation between the damage to the system, liquefaction severity, and the performance of different pipe materials. Finally, the characteristics of Christchurch liquefaction and its impacts on built environment are discussed in relation to the liquefaction-induced damage in Japan during the 11 March 2011 Great East Japan Earthquake.
Field trips are one of the most critical pieces of learning for students in sciences like geology, biology, and geography. Virtual field trips (VFT) are being increasingly considered as sophisticated and effective forms of teaching, especially with the rise of new technologies and the growing demand for more inclusive classroom environments. This research developed a virtual field trip for Tertiary students in an introductory-level geology course (GEOL 113: Environmental Geohazards) at the University of Canterbury. This initiative was in partnership with LEARNZ – a highly esteemed virtual fieldtrip team run by CORE Education that creates successful VFTs for Primary and Secondary students in New Zealand. Key components of the Tertiary VFT include a student acting as the virtual field trip teacher interviewing experts and leading the field trip, web-based background material, online assessment, and photos. In two successive academic years, students participated in the VFT during lectures and as pre class assignments prior to a one-day earthquake hazards workshop. In 2016, the virtual field trip used the LEARNZ web platform and occurred synchronously with the class; in 2017 the virtual fieldtrip reused the video, images and word documents from the previous year with the addition of a Google Earth component and with no reliance on the LEARNZ web platform. The goals of the trip were designed to prepare students for the earthquake hazards workshop, in which students analysed earthquake impacts over varying timescales and then applied that knowledge to develop strategies for the recovery of three crucial industries (dairy, mining, or tourism) on the West Coast of New Zealand’s South Island. In both years, number of clicks data showed that students interacted with online material far more during this week of the course than any other. Following the synchronous version in 2016, the students who were surveyed reported (1) they enjoyed the trip, (2) they found background material useful for preparation for the trip and the workshop, and (3) the additional work was at the appropriate level. Despite predominantly positive responses from the students, we experienced some negative feedback from participating staff mainly associated with stress and technical difficulties in running the synchronous VFT. With the asynchronous trip in 2017, staff reported a highly positive overall experience, with a perceived enhanced interaction with class during lecture time, and an increased and enhanced engagement with course material outside of class. The student survey again showed that the majority of students surveyed enjoyed the virtual fieldtrip, and that it was useful preparation for the workshop. Additionally, they reported an improved link between earth processes and society, which was a key overarching aim for the course. We propose that the synchronous version poses more excitement and immersion in the field environment, whereas the reuse of the asynchronous version increases the utility (and hence value for money) of the trip, and minimises technical difficulties and lecturer stress. Additionally, re-using the material in the asynchronous version offered opportunities to improve and supplement the past content, such as the incorporation of following an annotated trip path in Google Earth. As recommendations for others interested in developing virtual fieldtrips, we report that the design of a virtual fieldtrip should include (1) Goal-aligned content and assessment for both practice and marks, (2) a student and instructor experience that is authentic and flexible to both the people and the place. We suggest that these aims can be achieved whatever the budget or timeframe and make our material freely available at https://serc.carleton.edu/index.html.
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.
One of the most controversial issues highlighted by the 2010-2011 Christchurch earthquake series and more recently the 2016 Kaikoura earthquake, has been the evident difficulty and lack of knowledge and guidelines for: a) evaluation of the residual capacity damaged buildings to sustain future aftershocks; b) selection and implementation of a series of reliable repairing techniques to bring back the structure to a condition substantially the same as prior to the earthquake; and c) predicting the cost (or cost-effectiveness) of such repair intervention, when compared to fully replacement costs while accounting for potential aftershocks in the near future. As a result of such complexity and uncertainty (i.e., risk), in combination with the possibility (unique in New Zealand when compared to most of the seismic-prone countries) to rely on financial support from the insurance companies, many modern buildings, in a number exceeding typical expectations from past experiences at an international level, have ended up being demolished. This has resulted in additional time and indirect losses prior to the full reconstruction, as well as in an increase in uncertainty on the actual relocation of the investment. This research project provides the main end-users and stakeholders (practitioner engineers, owners, local and government authorities, insurers, and regulatory agencies) with comprehensive evidence-based information to assess the residual capacity of damage reinforced concrete buildings, and to evaluate the feasibility of repairing techniques, in order to support their delicate decision-making process of repair vs. demolition or replacement. Literature review on effectiveness of epoxy injection repairs, as well as experimental tests on full-scale beam-column joints shows that repaired specimens have a reduced initial stiffness compared with the undamaged specimen, with no apparent strength reduction, sometimes exhibiting higher displacement ductility capacities. Although the bond between the steel and concrete is only partially restored, it still allows the repaired specimen to dissipate at least the same amount of hysteretic energy. Experimental tests on buildings subjected to earthquake loading demonstrate that even for severe damage levels, the ability of the epoxy injection to restore the initial stiffness of the structure is significant. Literature review on damage assessment and repair guidelines suggests that there is consensus within the international community that concrete elements with cracks less than 0.2 mm wide only require cosmetic repairs; epoxy injection repairs of cracks less and 2.0 mm wide and concrete patching of spalled cover concrete (i.e., minor to moderate damage) is an appropiate repair strategy; and for severe damaged components (e.g., cracks greater than 2.0 mm wide, crushing of the concrete core, buckling of the longitudinal reinforcement) local replacement of steel and/or concrete in addition to epoxy crack injection is more appropriate. In terms of expected cracking patterns, non-linear finite element investigations on well-designed reinforced concrete beam-to-column joints, have shown that lower number of cracks but with wider openings are expected to occur for larger compressive concrete strength, f’c, and lower reinforcement content, ρs. It was also observed that the tensile concrete strength, ft, strongly affects the expected cracking pattern in the beam-column joints, the latter being more uniformly distributed for lower ft values. Strain rate effects do not seem to play an important role on the cracking pattern. However, small variations in the cracking pattern were observed for low reinforcement content as it approaches to the minimum required as per NZS 3101:2006. Simple equations are proposed in this research project to relate the maximum and residual crack widths with the steel strain at peak displacement, with or without axial load. A literature review on fracture of reinforcing steel due to low-cycle fatigue, including recent research using steel manufactured per New Zealand standards is also presented. Experimental results describing the influence of the cyclic effect on the ultimate strain capacity of the steel are also discussed, and preliminary equations to account for that effect are proposed. A literature review on the current practice to assess the seismic residual capacity of structures is also presented. The various factors affecting the residual fatigue life at a component level (i.e., plastic hinge) of well-designed reinforced concrete frames are discussed, and equations to quantify each of them are proposed, as well as a methodology to incorporate them into a full displacement-based procedure for pre-earthquake and post-earthquake seismic assessment.
This dissertation addresses several fundamental and applied aspects of ground motion selection for seismic response analyses. In particular, the following topics are addressed: the theory and application of ground motion selection for scenario earthquake ruptures; the consideration of causal parameter bounds in ground motion selection; ground motion selection in the near-fault region where directivity effect is significant; and methodologies for epistemic uncertainty consideration and propagation in the context of ground motion selection and seismic performance assessment. The paragraphs below outline each contribution in more detail. A scenario-based ground motion selection method is presented which considers the joint distribution of multiple intensity measure (IM) types based on the generalised conditional intensity measure (GCIM) methodology (Bradley, 2010b, 2012c). The ground motion selection algorithm is based on generating realisations of the considered IM distributions for a specific rupture scenario and then finding the prospective ground motions which best fit the realisations using an optimal amplitude scaling factor. In addition, using different rupture scenarios and site conditions, two important aspects of the GCIM methodology are scrutinised: (i) different weight vectors for the various IMs considered; and (ii) quantifying the importance of replicate selections for ensembles with different numbers of desired ground motions. As an application of the developed scenario-based ground motion selection method, ground motion ensembles are selected to represent several major earthquake scenarios in New Zealand that pose a significant seismic hazard, namely, Alpine, Hope and Porters Pass ruptures for Christchurch city; and Wellington, Ohariu, and Wairarapa ruptures for Wellington city. A rigorous basis is developed, and sensitivity analyses performed, for the consideration of bounds on causal parameters (e.g., magnitude, source-to-site distance, and site condition) for ground motion selection. The effect of causal parameter bound selection on both the number of available prospective ground motions from an initial empirical as-recorded database, and the statistical properties of IMs of selected ground motions are examined. It is also demonstrated that using causal parameter bounds is not a reliable approach to implicitly account for ground motion duration and cumulative effects when selection is based on only spectral acceleration (SA) ordinates. Specific causal parameter bounding criteria are recommended for general use as a ‘default’ bounding criterion with possible adjustments from the analyst based on problem-specific preferences. An approach is presented to consider the forward directivity effects in seismic hazard analysis, which does not separate the hazard calculations for pulse-like and non-pulse-like ground motions. Also, the ability of ground motion selection methods to appropriately select records containing forward directivity pulse motions in the near-fault region is examined. Particular attention is given to ground motion selection which is explicitly based on ground motion IMs, including SA, duration, and cumulative measures; rather than a focus on implicit parameters (i.e., distance, and pulse or non-pulse classifications) that are conventionally used to heuristically distinguish between the near-fault and far-field records. No ad hoc criteria, in terms of the number of directivity ground motions and their pulse periods, are enforced for selecting pulse-like records. Example applications are presented with different rupture characteristics, source-to-site geometry, and site conditions. It is advocated that the selection of ground motions in the near-fault region based on IM properties alone is preferred to that in which the proportion of pulse-like motions and their pulse periods are specified a priori as strict criteria for ground motion selection. Three methods are presented to propagate the effect of seismic hazard and ground motion selection epistemic uncertainties to seismic performance metrics. These methods differ in their level of rigor considered to propagate the epistemic uncertainty in the conditional distribution of IMs utilised in ground motion selection, selected ground motion ensembles, and the number of nonlinear response history analyses performed to obtain the distribution of engineering demand parameters. These methods are compared for an example site where it is observed that, for seismic demand levels below the collapse limit, epistemic uncertainty in ground motion selection is a smaller uncertainty contributor relative to the uncertainty in the seismic hazard itself. In contrast, uncertainty in ground motion selection process increases the uncertainty in the seismic demand hazard for near-collapse demand levels.
The purpose of this thesis is to conduct a detailed examination of the forward-directivity characteristics of near-fault ground motions produced in the 2010-11 Canterbury earthquakes, including evaluating the efficacy of several existing empirical models which form the basis of frameworks for considering directivity in seismic hazard assessment. A wavelet-based pulse classification algorithm developed by Baker (2007) is firstly used to identify and characterise ground motions which demonstrate evidence of forward-directivity effects from significant events in the Canterbury earthquake sequence. The algorithm fails to classify a large number of ground motions which clearly exhibit an early-arriving directivity pulse due to: (i) incorrect pulse extraction resulting from the presence of pulse-like features caused by other physical phenomena; and (ii) inadequacy of the pulse indicator score used to carry out binary pulse-like/non-pulse-like classification. An alternative ‘manual’ approach is proposed to ensure 'correct' pulse extraction and the classification process is also guided by examination of the horizontal velocity trajectory plots and source-to-site geometry. Based on the above analysis, 59 pulse-like ground motions are identified from the Canterbury earthquakes , which in the author's opinion, are caused by forward-directivity effects. The pulses are also characterised in terms of their period and amplitude. A revised version of the B07 algorithm developed by Shahi (2013) is also subsequently utilised but without observing any notable improvement in the pulse classification results. A series of three chapters are dedicated to assess the predictive capabilities of empirical models to predict the: (i) probability of pulse occurrence; (ii) response spectrum amplification caused by the directivity pulse; (iii) period and amplitude (peak ground velocity, PGV) of the directivity pulse using observations from four significant events in the Canterbury earthquakes. Based on the results of logistic regression analysis, it is found that the pulse probability model of Shahi (2013) provides the most improved predictions in comparison to its predecessors. Pulse probability contour maps are developed to scrutinise observations of pulses/non-pulses with predicted probabilities. A direct comparison of the observed and predicted directivity amplification of acceleration response spectra reveals the inadequacy of broadband directivity models, which form the basis of the near-fault factor in the New Zealand loadings standard, NZS1170.5:2004. In contrast, a recently developed narrowband model by Shahi & Baker (2011) provides significantly improved predictions by amplifying the response spectra within a small range of periods. The significant positive bias demonstrated by the residuals associated with all models at longer vibration periods (in the Mw7.1 Darfield and Mw6.2 Christchurch earthquakes) is likely due to the influence of basin-induced surface waves and non-linear soil response. Empirical models for the pulse period notably under-predict observations from the Darfield and Christchurch earthquakes, inferred as being a result of both the effect of nonlinear site response and influence of the Canterbury basin. In contrast, observed pulse periods from the smaller magnitude June (Mw6.0) and December (Mw5.9) 2011 earthquakes are in good agreement with predictions. Models for the pulse amplitude generally provide accurate estimates of the observations at source-to-site distances between 1 km and 10 km. At longer distances, observed PGVs are significantly under-predicted due to their slower apparent attenuation. Mixed-effects regression is employed to develop revised models for both parameters using the latest NGA-West2 pulse-like ground motion database. A pulse period relationship which accounts for the effect of faulting mechanism using rake angle as a continuous predictor variable is developed. The use of a larger database in model development, however does not result in improved predictions of pulse period for the Darfield and Christchurch earthquakes. In contrast, the revised model for PGV provides a more appropriate attenuation of the pulse amplitude with distance, and does not exhibit the bias associated with previous models. Finally, the effects of near-fault directivity are explicitly included in NZ-specific probabilistic seismic hazard analysis (PSHA) using the narrowband directivity model of Shahi & Baker (2011). Seismic hazard analyses are conducted with and without considering directivity for typical sites in Christchurch and Otira. The inadequacy of the near-fault factor in the NZS1170.5: 2004 is apparent based on a comparison with the directivity amplification obtained from PSHA.
Reinforced concrete structures designed in pre-1970s are vulnerable under earthquakes due to lack of seismic detailing to provide adequate ductility. Typical deficiencies of pre-1970s reinforced concrete structures are (a) use of plain bars as longitudinal reinforcement, (b) inadequate anchorage of beam longitudinal reinforcement in the column (particularly exterior column), (c) lack of joint transverse reinforcement if any, (d) lapped splices located just above joint, and (e) low concrete strength. Furthermore, the use of infill walls is a controversial issue because it can help to provide additional stiffness to the structure on the positive side and on the negative side it can increase the possibility of soft-storey mechanisms if it is distributed irregularly. Experimental research to investigate the possible seismic behaviour of pre-1970s reinforced concrete structures have been carried out in the past. However, there is still an absence of experimental tests on the 3-D response of existing beam-column joints under bi-directional cyclic loading, such as corner joints. As part of the research work herein presented, a series of experimental tests on beam-column subassemblies with typical detailing of pre-1970s buildings has been carried out to investigate the behaviour of existing reinforced concrete structures. Six two-third scale plane frame exterior beam-column joint subassemblies were constructed and tested under quasi-static cyclic loading in the Structural Laboratory of the University of Canterbury. The reinforcement detailing and beam dimension were varied to investigate their effect on the seismic behaviour. Four specimens were conventional deep beam-column joint, with two of them using deformed longitudinal bars and beam bars bent in to the joint and the two others using plain round longitudinal bars and beam bars with end hooks. The other two specimens were shallow beam-column joint, one with deformed longitudinal bars and beam bars bent in to the joint, the other with plain round longitudinal bars and beam bars with end hooks. All units had one transverse reinforcement in the joint. The results of the experimental tests indicated that conventional exterior beam-column joint with typical detailing of pre-1970s building would experience serious diagonal tension cracking in the joint panel under earthquake. The use of plain round bars with end hooks for beam longitudinal reinforcement results in more severe damage in the joint core when compared to the use of deformed bars for beam longitudinal reinforcement bent in to the joint, due to the combination of bar slips and concrete crushing. One interesting outcome is that the use of shallow beam in the exterior beam-column joint could avoid the joint cracking due to the beam size although the strength provided lower when compared with the use of deep beam with equal moment capacity. Therefore, taking into account the low strength and stiffness, shallow beam can be reintroduced as an alternative solution in design process. In addition, the presence of single transverse reinforcement in the joint core can provide additional confinement after the first crack occurred, thus delaying the strength degradation of the structure. Three two-third scale space frame corner beam-column joint subassemblies were also constructed to investigate the biaxial loading effect. Two specimens were deep-deep beam-corner column joint specimens and the other one was deep-shallow beam-corner column joint specimen. One deep-deep beam-corner column joint specimen was not using any transverse reinforcement in the joint core while the two other specimens were using one transverse reinforcement in the joint core. Plain round longitudinal bars were used for all units with hook anchorage for the beam bars. Results from the tests confirmed the evidences from earthquake damage observations with the exterior 3-D (corner) beam-column joint subjected to biaxial loading would have less strength and suffer higher damage in the joint area under earthquake. Furthermore, the joint shear relation in the two directions is calibrated from the results to provide better analysis. An analytical model was used to simulate the seismic behaviour of the joints with the help of Ruaumoko software. Alternative strength degradation curves corresponding to different reinforcement detailing of beam-column joint unit were proposed based on the test results.
This thesis is concerned with modelling rockfall parameters associated with cliff collapse debris and the resultant “ramp” that formed following the high peak ground acceleration (PGA) events of 22 February 2011 and 13 June 2011. The Christchurch suburb of Redcliffs, located at the base of the Port Hills on the northern side of Banks Peninsula, New Zealand, is comprised of Miocene-age volcanics with valley-floor infilling marine sediments. The area is dominated by basaltic lava flows of the Mt Pleasant Formation, which is a suite of rocks forming part of the Lyttelton Volcanic Group that were erupted 11.0-10.0Ma. Fresh exposure enabled the identification of a basaltic ignimbrite unit at the study site overlying an orange tuff unit that forms a marker horizon spanning the length of the field area. Prior to this thesis, basaltic ignimbrite on Banks Peninsula has not been recorded, so descriptions and interpretations of this unit are the first presented. Mapping of the cliff face by remote observation, and analysis of hand samples collected from the base of the debris slopes, has identified a very strong (>200MPa), columnar-jointed, welded unit, and a very weak (<5MPa), massive, so-called brecciated unit that together represent the end-member components of the basaltic ignimbrite. Geochemical analysis shows the welded unit is picrite basalt, and the brecciated unit is hawaiite, making both clearly distinguishable from the underlying trachyandesite tuff. RocFall™ 4.0 was used to model future rockfalls at Redcliffs. RocFall™ is a two-dimensional (2D), hybrid, probabilistic modelling programme for which topographical profile data is used to generate slope profiles. GNS Science collected the data used for slope profile input in March 2011. An initial sensitivity analysis proved the Terrestrial Laser Scan (TLS)-derived slope to be too detailed to show any results when the slope roughness parameter was tested. A simplified slope profile enabled slope roughness to be varied, however the resulting model did not correlate with field observations as well. By using slope profile data from March 2011, modelled rockfall behaviour has been calibrated with observed rockfall runout at Redcliffs in the 13 June 2011 event to create a more accurate rockfall model. The rockfall model was developed on a single slope profile (Section E), with the chosen model then applied to four other section lines (A-D) to test the accuracy of the model, and to assess future rockfall runout across a wider area. Results from Section Lines A, B, and E correlate very well with field observations, with <=5% runout exceeding the modelled slope, and maximum bounce height at the toe of the slope <=1m. This is considered to lie within observed limits given the expectation that talus slopes will act as a ramp on which modelled rocks travel further downslope. Section Lines C and D produced higher runout percentage values than the other three section lines (23% and 85% exceeding the base of the slope, respectively). Section D also has a much higher maximum bounce height at the toe of the slope (~8.0m above the slope compared to <=1.0m for the other four sections). Results from modelling of all sections shows the significance of the ratio between total cliff height (H) and horizontal slope distance (x), and of maximum drop height to the top of the talus (H*) and horizontal slope distance (x). H/x can be applied to the horizontal to vertical ratio (H:V) as used commonly to identify potential slope instability. Using the maximum value from modelling at Redcliffs, the future runout limit can be identified by applying a 1.4H:1V ratio to the remainder of the cliff face. Additionally, the H*/x parameter shows that when H*/x >=0.6, the percentage of rock runout passing the toe of the slope will exceed 5%. When H*/x >=0.75, the maximum bounce height at the toe of the slope can be far greater than when H*/x is below this threshold. Both of these parameters can be easily obtained, and can contribute valuable guideline data to inform future land-use planning decisions. This thesis project has demonstrated the applicability of a 2D probabilistic-based model (RocFall™ 4.0) to evaluate rockfall runout on the talus slope (or ramp) at the base of ~35-70m high cliff with a basaltic ignimbrite source. Limitations of the modelling programme have been identified, in particular difficulties with adjusting modelled roughness of the slope profile and the inability to consider fragmentation. The runout profile using RocFall™ has been successfully calibrated against actual profiles and some anomalous results have been identified.
