Asset management in power systems is exercised to improve network reliability to provide confidence and security for customers and asset owners. While there are well-established reliability metrics that are used to measure and manage business-as-usual disruptions, an increasing appreciation of the consequences of low-probability high-impact events means that resilience is increasingly being factored into asset management in order to provide robustness and redundancy to components and wider networks. This is particularly important for electricity systems, given that a range of other infrastructure lifelines depend upon their operation. The 2010-2011 Canterbury Earthquake Sequence provides valuable insights into electricity system criticality and resilience in the face of severe earthquake impacts. While above-ground assets are relatively easy to monitor and repair, underground assets such as cables emplaced across wide areas in the distribution network are difficult to monitor, identify faults on, and repair. This study has characterised in detail the impacts to buried electricity cables in Christchurch resulting from seismically-induced ground deformation caused primarily by liquefaction and lateral spread. Primary modes of failure include cable bending, stretching, insulation damage, joint braking and, being pulled off other equipment such as substation connections. Performance and repair data have been compiled into a detailed geospatial database, which in combination with spatial models of peak ground acceleration, peak ground velocity and ground deformation, will be used to establish rigorous relationships between seismicity and performance. These metrics will be used to inform asset owners of network performance in future earthquakes, further assess component criticality, and provide resilience metrics.
We present initial results from a set of three-dimensional (3D) deterministic earthquake ground motion simulations for the northern Canterbury plains, Christchurch and the Banks Peninsula region, which explicitly incorporate the effects of the surface topography. The simu-lations are done using Hercules, an octree-based finite-element parallel software for solving 3D seismic wave propagation problems in heterogeneous media under kinematic faulting. We describe the efforts undertaken to couple Hercules with the South Island Velocity Model (SIVM), which included changes to the SIVM code in order to allow for single repetitive que-ries and thus achieve a seamless finite-element meshing process within the end-to-end ap-proach adopted in Hercules. We present our selection of the region of interest, which corre-sponds to an area of about 120 km × 120 km, with the 3D model reaching a depth of 60 km. Initial simulation parameters are set for relatively high minimum shear wave velocity and a low maximum frequency, which we are progressively scaling up as computing resources permit. While the effects of topography are typically more important at higher frequencies and low seismic velocities, even at this initial stage of our efforts (with a maximum of 2 Hz and a mini-mum of 500 m/s), it is possible to observe the importance of the topography in the response of some key locations within our model. To highlight these effects we compare the results of the 3D topographic model with respect to those of a flat (squashed) 3D model. We draw rele-vant conclusions from the study of topographic effects during earthquakes for this region and describe our plans for future work.
Results from a series of 1D seismic effective stress analyses of natural soil deposits from Christchurch are summarized. The analysed soil columns include sites whose performance during the 2010-2011 Canterbury earthquakes varied significantly, from no liquefaction manifestation at the ground surface to very severe liquefaction, in which case a large area of the site was covered by thick soil ejecta. Key soil profile characteristics and response mechanisms affecting the severity of surface liquefaction manifestation and subsequent damage are explored. The influence of shaking intensity on the triggering and contribution of these mechanisms is also discussed. Careful examination of the results highlights the importance of considering the deposit as a whole, i.e. a system of layers, including interactions between layers in the dynamic response and through pore water pressure redistribution and water flow.
High demolition rates were observed in New Zealand after the 2010-2011 Canterbury Earthquake Sequence despite the success of modern seismic design standards to achieve required performance objectives such as life safety and collapse prevention. Approximately 60% of the multi-storey reinforced concrete (RC) buildings in the Christchurch Central Business District were demolished after these earthquakes, even when only minor structural damage was present. Several factors influenced the decision of demolition instead of repair, one of them being the uncertainty of the seismic capacity of a damaged structure. To provide more insight into this topic, the investigation conducted in this thesis evaluated the residual capacity of moderately damaged RC walls and the effectiveness of repair techniques to restore the seismic performance of heavily damaged RC walls. The research outcome provided insights for developing guidelines for post-earthquake assessment of earthquake-damaged RC structures. The methodology used to conduct the investigation was through an experimental program divided into two phases. During the first phase, two walls were subjected to different types of pre-cyclic loading to represent the damaged condition from a prior earthquake, and a third wall represented a repair scenario with the damaged wall being repaired using epoxy injection and repair mortar after the pre-cyclic loading. Comparisons of these test walls to a control undamaged wall identified significant reductions in the stiffness of the damaged walls and a partial recovery in the wall stiffness achieved following epoxy injection. Visual damage that included distributed horizontal and diagonal cracks and spalling of the cover concrete did not affect the residual strength or displacement capacity of the walls. However, evidence of buckling of the longitudinal reinforcement during the pre-cyclic loading resulted in a slight reduction in strength recovery and a significant reduction in the displacement capacity of the damaged walls. Additional experimental programs from the literature were used to provide recommendations for modelling the response of moderately damaged RC walls and to identify a threshold that represented a potential reduction in the residual strength and displacement capacity of damaged RC walls in future earthquakes. The second phase of the experimental program conducted in this thesis addressed the replacement of concrete and reinforcing steel as repair techniques for heavily damaged RC walls. Two walls were repaired by replacing the damaged concrete and using welded connections to connect new reinforcing bars with existing bars. Different locations of the welded connections were investigated in the repaired walls to study the impact of these discontinuities at the critical section. No significant changes were observed in the stiffness, strength, and displacement capacity of the repaired walls compared to the benchmark undamaged wall. Differences in the local behaviour at the critical section were observed in one of the walls but did not impact the global response. The results of these two repaired walls were combined with other experimental programs found in the literature to assemble a database of repaired RC walls. Qualitative and quantitative analyses identified trends across various parameters, including wall types, damage before repair, and repair techniques implemented. The primary outcome of the database analysis was recommendations for concrete and reinforcing steel replacement to restore the strength and displacement capacity of heavily damaged RC walls.
Damage distribution maps from strong earthquakes and recorded data from field experiments have repeatedly shown that the ground surface topography and subsurface stratigraphy play a decisive role in shaping the ground motion characteristics at a site. Published theoretical studies qualitatively agree with observations from past seismic events and experiments; quantitatively, however, they systematically underestimate the absolute level of topographic amplification up to an order of magnitude or more in some cases. We have hypothesized in previous work that this discrepancy stems from idealizations of the geometry, material properties, and incident motion characteristics that most theoretical studies make. In this study, we perform numerical simulations of seismic wave propagation in heterogeneous media with arbitrary ground surface geometry, and compare results with high quality field recordings from a site with strong surface topography. Our goal is to explore whether high-fidelity simulations and realistic numerical models can – contrary to theoretical models – capture quantitatively the frequency and amplitude characteristics of topographic effects. For validation, we use field data from a linear array of nine portable seismometers that we deployed on Mount Pleasant and Heathcote Valley, Christchurch, New Zealand, and we compute empirical standard spectral ratios (SSR) and single-station horizontal-to-vertical spectral ratios (HVSR). The instruments recorded ambient vibrations and remote earthquakes for a period of two months (March-April 2017). We next perform two-dimensional wave propagation simulations using the explicit finite difference code FLAC. We construct our numerical model using a high-resolution (8m) Digital Elevation Map (DEM) available for the site, an estimated subsurface stratigraphy consistent with the geomorphology of the site, and soil properties estimated from in-situ and non-destructive tests. We subject the model to in-plane and out-of-plane incident motions that span a broadband frequency range (0.1-20Hz). Numerical and empirical spectral ratios from our blind prediction are found in very good quantitative agreement for stations on the slope of Mount Pleasant and on the surface of Heathcote Valley, across a wide range of frequencies that reveal the role of topography, soil amplification and basin edge focusing on the distribution of ground surface motion.
Natural hazards continue to have adverse effects on communities and households worldwide, accelerating research on proactively identifying and enhancing characteristics associated with resilience. Although resilience is often characterized as a return to normal, recent studies of postdisaster recovery have highlighted the ways in which new opportunities can emerge following disruption, challenging the status quo. Conversely, recovery and reconstruction may serve to reinforce preexisting social, institutional, and development pathways. Our understanding of these dynamics is limited however by the small number of practice examples, particularly for rural communities in developed nations. This study uses a social–ecological inventory to document the drivers, pathways, and mechanisms of resilience following a large-magnitude earthquake in Kaikōura, a coastal community in Aotearoa New Zealand. As part of the planning and implementation phase of a multiyear project, we used the tool as the basis for indepth and contextually sensitive analysis of rural resilience. Moreover, the deliberate application of social–ecological inventory was the first step in the research team reengaging with the community following the event. The inventory process provided an opportunity for research partners to share their stories and experiences and develop a shared understanding of changes that had taken place in the community. Results provide empirical insight into reactions to disruptive change associated with disasters. The inventory also informed the design of targeted research collaborations, established a platform for longer-term community engagement, and provides a baseline for assessing longitudinal changes in key resilience-related characteristics and community capacities. Findings suggest the utility of social–ecological inventory goes beyond natural resource management, and that it may be appropriate in a range of contexts where institutional, social, and economic restructuring have developed out of necessity in response to felt or anticipated external stressors.
The 2010-2011 Canterbury earthquake sequence, and the resulting extensive data sets on damaged buildings that have been collected, provide a unique opportunity to exercise and evaluate previously published seismic performance assessment procedures. This poster provides an overview of the authors’ methodology to perform evaluations with two such assessment procedures, namely the P-58 guidelines and the REDi Rating System. P-58, produced by the Federal Emergency Management Agency (FEMA) in the United States, aims to facilitate risk assessment and decision-making by quantifying earthquake ground shaking, structural demands, component damage and resulting consequences in a logical framework. The REDi framework, developed by the engineering firm ARUP, aids stakeholders in implementing resilience-based earthquake design. Preliminary results from the evaluations are presented. These have the potential to provide insights on the ability of the assessment procedures to predict impacts using “real-world” data. However, further work remains to critically analyse these results and to broaden the scope of buildings studied and of impacts predicted.
We examined the stratigraphy of alluvial fans formed at the steep range front of the Southern Alps at Te Taho, on the north bank of the Whataroa River in central West Coast, South Island, New Zealand. The range front coincides with the Alpine Fault, an Australian-Pacific plate boundary fault, which produces regular earthquakes. Our study of range front fans revealed aggradation at 100- to 300-year intervals. Radiocarbon ages and soil residence times (SRTs) estimated by a quantitative profile development index allowed us to elucidate the characteristics of four episodes of aggradation since 1000 CE. We postulate a repeating mode of fan behaviour (fan response cycle [FRC]) linked to earthquake cycles via earthquake-triggered landslides. FRCs are characterised by short response time (aggradation followed by incision) and a long phase when channels are entrenched and fan surfaces are stable (persistence time). Currently, the Te Taho and Whataroa River fans are in the latter phase. The four episodes of fan building we determined from an OxCal sequence model correlate to Alpine Fault earthquakes (or other subsidiary events) and support prior landscape evolution studies indicating ≥M7.5 earthquakes as the main driver of episodic sedimentation. Our findings are consistent with other historic non-earthquake events on the West Coast but indicate faster responses than other earthquake sites in New Zealand and elsewhere where rainfall and stream gradients (the basis for stream power) are lower. Judging from the thickness of fan deposits and the short response times, we conclude that pastoral farming (current land-use) on the fans and probably across much of the Whataroa River fan would be impossible for several decades after a major earthquake. The sustainability of regional tourism and agriculture is at risk, more so because of the vulnerability of the single through road in the region (State Highway 6).
Buildings subject to earthquake shaking will tend to move not only horizontally but also rotate in plan. In-plan rotation is known as “building torsion” and it may occur for a variety of reasons, including stiffness and strength eccentricity and/or torsional effects from ground motions. Methods to consider torsion in structural design standards generally involve analysis of the structure in its elastic state. This is despite the fact that the structural elements can yield, thereby significantly altering the building response and the structural element demands. If demands become too large, the structure may collapse. While a number of studies have been conducted into the behavior of structures considering inelastic building torsion, there appears to be no consensus that one method is better than another and as a result, provisions within current design standards have not adopted recent proposals in the literature. However, the Canterbury Earthquakes Royal Commission recently made the recommendation that provisions to account for inelastic torsional response of buildings be introduced within New Zealand building standards. Consequently, this study examines how and to what extent the torsional response due to system eccentricity may affect the seismic performance of a building and considers what a simple design method should account for. It is concluded that new methods should be simple, be applicable to both the elastic and inelastic range of response, consider bidirectional excitation and include guidance for multi-story systems.
This study explicitly investigates uncertainties in physics-based ground motion simulation validation for earthquakes in the Canterbury region. The simulations utilise the Graves and Pitarka (2015) hybrid methodology, with separately quantified parametric uncertainties in the comprehensive physics and simplified physics components of the model. The study is limited to the simulation of 148 small magnitude (Mw 3.5 – 5) earthquakes, with a point source approximation for the source rupture representations, which also enables a focus on a small number of relevant uncertainties. The parametric uncertainties under consideration were selected through sensitivity analysis, and specifically include: magnitude, Brune stress parameter and high frequency rupture velocity. Twenty Monte Carlo realisations were used to sample parameter uncertainties for each of the 148 events. Residuals associated with the following intensity measures: spectral acceleration, peak ground velocity, arias intensity and significant duration, were ascertained. Using these residuals, validation was performed through assessment of systematic biases in site and source terms from mixed-effects regression. Based on the results to date, initial standard deviation recommendations for parameter uncertainties, based on the Canterbury simulations have been obtained. This work ultimately provides an initial step toward explicit incorporation of modelling uncertainty in simulated ground motion predictions for future events, which will improve the use of simulation models in seismic hazard analysis. We plan to subsequently assess uncertainties for larger magnitude events with more complex ruptures, and events across a larger geographic region, as well as uncertainties due to path attenuation, site effects, and more general model epistemic uncertainties.
A review of the literature showed the lack of a truly effective damage avoidance solution for timber or hybrid timber moment resisting frames (MRFs). Full system damage avoidance selfcentring behaviour is difficult to achieve with existing systems due to damage to the floor slab caused by beam-elongation. A novel gravity rocking, self-centring beam-column joint with inherent and supplemental friction energy dissipation is proposed for low-medium rise buildings in all seismic zones where earthquake actions are greater than wind. Steel columns and timber beams are used in the hybrid MRF such that both the beam and column are continuous thus avoiding beam-elongation altogether. Corbels on the columns support the beams and generate resistance and self-centring through rocking under the influence of gravity. Supplemental friction sliders at the top of the beams resist sliding of the floor whilst dissipating energy as the floor lifts on the corbels and returns. 1:20 scale tests of 3-storey one-by-two bay building based on an earlier iteration of the proposed concept served as proof-of-concept and highlighted areas for improvement. A 1:5 scale 3-storey one-by-one bay building was subsequently designed. Sub-assembly tests of the beam-top asymmetric friction sliders demonstrated repeatable hysteresis. Quasi-static tests of the full building demonstrated a ‘flat bottomed’ flag-shaped hysteresis. Shake table tests to a suite of seven earthquakes scaled for Wellington with site soil type D to the serviceability limit state (SLS), ultimate limit state (ULS) and maximum credible event (MCE) intensity corresponding to an average return period of 25, 500 and 2500 years respectively were conducted. Additional earthquake records from the 22 February 2011 Christchurch earthquakes we included. A peak drift of 0.6%, 2.5% and 3.8% was reached for the worst SLS, ULS and MCE earthquake respectively whereas a peak drift of 4.5% was reached for the worst Christchurch record for tests in the plane of the MRF. Bi-directional tests were also conducted with the building oriented at 45 degrees on the shake table and the excitation factored by 1.41 to maintain the component in the direction of the MRF. Shear walls with friction slider hold-downs which reached similar drifts to the MRF were provided in the orthogonal direction. Similar peak drifts were reached by the MRF in the bi-directional tests, when the excitation was amplified as intended. The building self-centred with a maximum residual drift of 0.06% in the dynamic tests and demonstrated no significant damage. The member actions were magnified by up to 100% due to impact upon return of the floor after uplift when the peak drift reached 4.5%. Nonetheless, all of the members and connections remained essentially linearelastic. The shake table was able to produce a limited peak velocity of 0.275 m/s and this limited the severity of several of the ULS, MCE and Christchurch earthquakes, especially the near-field records with a large velocity pulse. The full earthquakes with uncapped velocity were simulated in a numerical model developed in SAP2000. The corbel supports were modelled with the friction isolator link element and the top sliders were modelled with a multi-linear plastic link element in parallel with a friction spring damper. The friction spring damper simulated the increase in resistance with increasing joint rotation and a near zero return stiffness, as exhibited by the 1:5 scale test building. A good match was achieved between the test quasi-static global force-displacement response and the numerical model, except a less flat unloading curve in the numerical model. The peak drift from the shake table tests also matched well. Simulations were also run for the full velocity earthquakes, including vertical ground acceleration and different floor imposed load scenarios. Excessive drift was predicted by the numerical model for the full velocity near-field earthquakes at the MCE intensity and a rubber stiffener for increasing the post joint-opening stiffness was found to limit the drift to 4.8%. Vertical ground acceleration had little effect on the global response. The system generates most of its lateral resistance from the floor weight, therefore increasing the floor imposed load increased the peak drift, but less than it would if the resistance of the system did not increase due to the additional floor load. A seismic design procedure was discussed under the framework of the existing direct displacement-based design method. An expression for calculating the area-based equivalent viscous damping (EVD) was derived and a conservative correction factor of 0.8 was suggested. A high EVD of up to about 15% can be achieved with the proposed system at high displacement ductility levels if the resistance of the top friction sliders is maximised without compromising reliable return of the floor after uplift. Uniform strength joints with an equal corbel length up the height of the building and similar inter-storey drifts result in minimal relative inter-floor uplift, except between the first floor and ground. Guidelines for detailing the joint for damage avoidance including bi-directional movement were also developed.
