Semi-empirical models based on in-situ geotechnical tests have become the standard of practice for predicting soil liquefaction. Since the inception of the “simplified” cyclic-stress model in 1971, variants based on various in-situ tests have been developed, including the Cone Penetration Test (CPT). More recently, prediction models based soley on remotely-sensed data were developed. Similar to systems that provide automated content on earthquake impacts, these “geospatial” models aim to predict liquefaction for rapid response and loss estimation using readily-available data. This data includes (i) common ground-motion intensity measures (e.g., PGA), which can either be provided in near-real-time following an earthquake, or predicted for a future event; and (ii) geospatial parameters derived from digital elevation models, which are used to infer characteristics of the subsurface relevent to liquefaction. However, the predictive capabilities of geospatial and geotechnical models have not been directly compared, which could elucidate techniques for improving the geospatial models, and which would provide a baseline for measuring improvements. Accordingly, this study assesses the realtive efficacy of liquefaction models based on geospatial vs. CPT data using 9,908 case-studies from the 2010-2016 Canterbury earthquakes. While the top-performing models are CPT-based, the geospatial models perform relatively well given their simplicity and low cost. Although further research is needed (e.g., to improve upon the performance of current models), the findings of this study suggest that geospatial models have the potential to provide valuable first-order predictions of liquefaction occurence and consequence. Towards this end, performance assessments of geospatial vs. geotechnical models are ongoing for more than 20 additional global earthquakes.
This research aims to explore how business models of SMEs revolve in the face of a crisis to be resilient. The business model canvas was used as a tool to analyse business models of SMEs in Greater Christchurch. The purpose was to evaluate the changes SMEs brought in their business models after hit by a series of earthquake in 2010 and 2011. The idea was to conduct interviews of business owners and analyse using grounded theory methods. Because this method is iterative, a tentative theoretical framework was proposed, half way through the data collection. It was realised that owner specific characteristics were more prominent in the data than the elements business model. Although, SMEs in this study experienced several operational changes in their business models such as change of location and modification of payment terms. However, the suggested framework highlights how owner specific attributes influence the survival of a small business. Small businesses and their owners are extremely interrelated that the business models personify the owner specific characteristics. In other words, the adaptation of the business model reflects the extent to which the owner possess these attributes. These attributes are (a) Mindsets – the attitude and optimism of business owner; (b) Adaptive coping – the ability of business owner to take corrective actions; and (c) Social capital – the network of a business owner, including family, friends, neighbours and business partners.
Geospatial liquefaction models aim to predict liquefaction using data that is free and readily-available. This data includes (i) common ground-motion intensity measures; and (ii) geospatial parameters (e.g., among many, distance to rivers, distance to coast, and Vs30 estimated from topography) which are used to infer characteristics of the subsurface without in-situ testing. Since their recent inception, such models have been used to predict geohazard impacts throughout New Zealand (e.g., in conjunction with regional ground-motion simulations). While past studies have demonstrated that geospatial liquefaction-models show great promise, the resolution and accuracy of the geospatial data underlying these models is notably poor. As an example, mapped rivers and coastlines often plot hundreds of meters from their actual locations. This stems from the fact that geospatial models aim to rapidly predict liquefaction anywhere in the world and thus utilize the lowest common denominator of available geospatial data, even though higher quality data is often available (e.g., in New Zealand). Accordingly, this study investigates whether the performance of geospatial models can be improved using higher-quality input data. This analysis is performed using (i) 15,101 liquefaction case studies compiled from the 2010-2016 Canterbury Earthquakes; and (ii) geospatial data readily available in New Zealand. In particular, we utilize alternative, higher-quality data to estimate: locations of rivers and streams; location of coastline; depth to ground water; Vs30; and PGV. Most notably, a region-specific Vs30 model improves performance (Figs. 3-4), while other data variants generally have little-to-no effect, even when the “standard” and “high-quality” values differ significantly (Fig. 2). This finding is consistent with the greater sensitivity of geospatial models to Vs30, relative to any other input (Fig. 5), and has implications for modeling in locales worldwide where high quality geospatial data is available.
Background This study examines the performance of site response analysis via nonlinear total-stress 1D wave-propagation for modelling site effects in physics-based ground motion simulations of the 2010-2011 Canterbury, New Zealand earthquake sequence. This approach allows for explicit modeling of 3D ground motion phenomena at the regional scale, as well as detailed nonlinear site effects at the local scale. The approach is compared to a more commonly used empirical VS30 (30 m time-averaged shear wave velocity)-based method for computing site amplification as proposed by Graves and Pitarka (2010, 2015), and to empirical ground motion prediction via a ground motion model (GMM).
In this paper we apply Full waveform tomography (FWT) based on the Adjoint-Wavefield (AW) method to iteratively invert a 3-D geophysical velocity model for the Canterbury region (Lee, 2017) from a simple initial model. The seismic wavefields was generated using numerical solution of the 3-D elastodynamic/ visco- elastodynamic equations (EMOD3D was adopted (Graves, 1996)), and through the AW method, gradients of model parameters (compression and shear wave velocity) were computed by implementing the cross-adjoint of forward and backward wavefields. The reversed-in-time displacement residual was utilized as the adjoint source. For inversion, we also account for the near source/ station effects, gradient precondition, smoothening (Gaussian filter in spatial domain) and optimal step length. Simulation-to-observation misfit measurements based on 191 sources at 78 seismic stations in the Canterbury region (Figure 1) were used into our inversion. The inversion process includes multiple frequency bands, starting from 0-0.05Hz, and advancing to higher frequency bands (0-0.1Hz and 0-0.2Hz). Each frequency band was used for up to 10 iterations or no optimal step length found. After 3 FWT inversion runs, the simulated seismograms computed using our final model show a good matching with the observed seismograms at frequencies from 0 - 0.2 Hz and the normalized least-squared misfit error has been significantly reduced. Over all, the synthetic study of FWT shows a good application to improve the crustal velocity models from the existed geological models and the seismic data of the different earthquake events happened in the Canterbury region.
Objectives • To develop a system dynamics model of Christchurch post-quake reconstruction process that captures all the critical dynamics influencing its pathway • To investigate the implications of current rebuild pathway • To build a reconstruction module to be integrated in MERIT (Measuring the Economics of Resilient Infrastructure Tool)
Well-validated liquefaction constitutive models are increasingly important as non-linear time history analyses become relatively more common in industry for key projects. Previous validation efforts of PM4Sand, a plasticity model specifically for liquefaction, have generally focused on centrifuge tests; however, pore pressure transducers installed at several free-field sites during the Canterbury Earthquake Sequence (CES) in Christchurch, New Zealand provide a relatively unique dataset to validate against. This study presents effective stress site response analyses performed in the finite difference software FLAC to examine the capability of PM4Sand to capture the generation of excess pore pressures during earthquakes. The characterization of the subsurface is primarily based on extensive cone penetration tests (CPT) carried out in Christchurch. Correlations based on penetration resistances are used to estimate soil parameters, such as relative density and shear wave velocity, which affect liquefaction behaviour. The resulting free-field FLAC model is used to estimate time histories of excess pore pressure, which are compared with records during several earthquakes in the CES to assess the suitability of PM4Sand.
Natural catastrophes are increasing worldwide. They are becoming more frequent but also more severe and impactful on our built environment leading to extensive damage and losses. Earthquake events account for the smallest part of natural events; nevertheless seismic damage led to the most fatalities and significant losses over the period 1981-2016 (Munich Re). Damage prediction is helpful for emergency management and the development of earthquake risk mitigation projects. Recent design efforts focused on the application of performance-based design engineering where damage estimation methodologies use fragility and vulnerability functions. However, the approach does not explicitly specify the essential criteria leading to economic losses. There is thus a need for an improved methodology that finds the critical building elements related to significant losses. The here presented methodology uses data science techniques to identify key building features that contribute to the bulk of losses. It uses empirical data collected on site during earthquake reconnaissance mission to train a machine learning model that can further be used for the estimation of building damage post-earthquake. The first model is developed for Christchurch. Empirical building damage data from the 2010-2011 earthquake events is analysed to find the building features that contributed the most to damage. Once processed, the data is used to train a machine-learning model that can be applied to estimate losses in future earthquake events.
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.
1. Background and Objectives This poster presents results from ground motion simulations of small-to-moderate magnitude (3.5≤Mw≤5.0) earthquake events in the Canterbury, New Zealand region using the Graves and Pitarka (2010,2015) methodology. Subsequent investigation of systematic ground motion effects highlights the prediction bias in the simulations which are also benchmarked against empirical ground motion models (e.g. Bradley (2013)). In this study, 144 earthquake ruptures, modelled as point sources, are considered with 1924 quality-assured ground motions recorded across 45 strong motion stations throughout the Canterbury region, as shown in Figure 1. The majority of sources are Mw≥4.0 and have centroid depth (CD) 10km or shallower. Earthquake source descriptions were obtained from the GeoNet New Zealand earthquake catalogue. The ground motion simulations were performed within a computational domain of 140km x 120km x 46km with a finite difference grid spacing of 0.1km. The low-frequency (LF) simulations utilize the 3D Canterbury Velocity Model while the high-frequency (HF) simulations utilize a generic regional 1D velocity model. In the LF simulations, a minimum shear wave velocity of 500m/s is enforced, yielding a maximum frequency of 1.0Hz.
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.
Tsunami events including the 2004 Indian Ocean Tsunami and the 2011 Tohoku Earthquake and Tsunami confirmed the need for Pacific-wide comprehensive risk mitigation and effective tsunami evacuation planning. New Zealand is highly exposed to tsunamis and continues to invest in tsunami risk awareness, readiness and response across the emergency management and science sectors. Evacuation is a vital risk reduction strategy for preventing tsunami casualties. Understanding how people respond to warnings and natural cues is an important element to improving evacuation modelling techniques. The relative rarity of tsunami events locally in Canterbury and also globally, means there is limited knowledge on tsunami evacuation behaviour, and tsunami evacuation planning has been largely informed by hurricane evacuations. This research aims to address this gap by analysing evacuation behaviour and movements of Kaikōura and Southshore/New Brighton (coastal suburb of Christchurch) residents following the 2016 Kaikōura earthquake. Stage 1 of the research is engaging with both these communities and relevant hazard management agencies, using a survey and community workshops to understand real-event evacuation behaviour during the 2016 Kaikōura earthquake and subsequent tsunami evacuations. The second stage is using the findings from stage 1 to inform an agent-based tsunami evacuation model, which is an approach that simulates of the movement of people during an evacuation response. This method improves on other evacuation modelling approaches to estimate evacuation times due to better representation of local population characteristics. The information provided by the communities will inform rules and interactions such as traffic congestion, evacuation delay times and routes taken to develop realistic tsunami evacuation models. This will allow emergency managers to more effectively prepare communities for future tsunami events, and will highlight recommended actions to increase the safety and efficiency of future tsunami evacuations.
Many buildings with relatively low damage from the 2010-2011 Canterbury were deemed uneconomic to repair and were replaced [1,2]. Factors that affected commercial building owners’ decisions to replace rather than repair, included capital availability, uncertainty with regards to regional recovery, local market conditions and ability to generate cash flow, and repair delays due to limited property access (cordon). This poster provides a framework for modeling decision-making in a case where repair is feasible but replacement might offer greater economic value – a situation not currently modeled in engineering risk analysis.
Our poster will present on-going QuakeCoRE-founded work on strong motion seismology for Dunedin-Mosgiel area, focusing on ground motion simulations for Dunedin Central Business District (CBD). Source modelling and ground motion simulations are being carried out using the SCEC (Southern California Earthquakes Center) Broad Band simulation Platform (BBP). The platform computes broadband (0-10 Hz) seismograms for earthquakes and was first implemented at the University of Otago in 2016. As large earthquakes has not been experienced in Dunedin in the time of period of instrumental recording, user-specified scenario simulations are of great value. The Akatore Fault, the most active fault in Otago and closest major fault to Dunedin, is the source focused on in the present study. Simulations for various Akatore Fault source scenarios are run and presented. Path and site effects are key components considered in the simulation process. A 1D shear wave velocity profile is required by SCEC BBP, and this is being generated to represent the Akatore-to-CBD path and site within the BBP. A 3D shear velocity model, with high resolution within Dunedin CBD, is being developed in parallel with this study (see Sangster et al. poster). This model will be the basis for developing a 3D shear wave velocity model for greater Dunedin-Mosgiel area for future ground motion simulations, using Canterbury software (currently under development).
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.
Motivation This poster aims to present fragility functions for pipelines buried in liquefaction-prone soils. Existing fragility models used to quantify losses can be based on old data or use complex metrics. Addressing these issues, the proposed functions are based on the Christchurch network and soil and utilizes the Canterbury earthquake sequence (CES) data, partially represented in Figure 1. Figure 1 (a) presents the pipe failure dataset, which describes the date, location and pipe on which failures occurred. Figure 1 (b) shows the simulated ground motion intensity median of the 22nd February 2011 earthquake. To develop the model, the network and soil characteristics have also been utilized.
This poster discusses several possible approaches by which the nonlinear response of surficial soils can be explicitly modelled in physics-based ground motion simulations, focusing on the relative advantages and limitations of the various methodologies. These methods include fully-coupled 3D simulation models that directly allow soil nonlinearity in surficial soils, the domain reduction method for decomposing the physical domain into multiple subdomains for separate simulation, conventional site response analysis uncoupled from the simulations, and finally, the use of simple empirically based site amplification factors We provide the methodology for an ongoing study to explicitly incorporate soil nonlinearity into hybrid broadband simulations of the 2010-2011 Canterbury, New Zealand earthquakes.
Liquefaction-induced lateral spreading during the 2011 Christchurch earthquake in New Zealand was severe and extensive, and data regarding the displacements associated with the lateral spreading provides an excellent opportunity to better understand the factors that influence these movements. Horizontal displacements measured from optical satellite imagery and subsurface data from the New Zealand Geotechnical Database (NZGD) were used to investigate four distinct lateral spread areas along the Avon River in Christchurch. These areas experienced displacements between 0.5 and 2 m, with the inland extent of displacement ranging from 100 m to over 600 m. Existing empirical and semi-empirical displacement models tend to under estimate displacements at some sites and over estimate at others. The integrated datasets indicate that the areas with more severe and spatially extensive displacements are associated with thicker and more laterally continuous deposits of liquefiable soil. In some areas, the inland extent of displacements is constrained by geologic boundaries and geomorphic features, as expressed by distinct topographic breaks. In other areas the extent of displacement is influenced by the continuity of liquefiable strata or by the presence of layers that may act as vertical seepage barriers. These observations demonstrate the need to integrate geologic/geomorphic analyses with geotechnical analyses when assessing the potential for lateral spreading movements.
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.
Despite the relatively low seismicity, a large earthquake in the Waikato region is expected to have a high impact, when the fourth-largest regional population and economy and the high density critical infrastructure systems in this region are considered. Furthermore, Waikato has a deep soft sedimentary basin, which increases the regional seismic hazard due to trapping and amplification of seismic waves and generation of localized surface waves within the basin. This phenomenon is known as the “Basin Effect”, and has been attributed to the increased damage in several historic earthquakes, including the 2010-2011 Canterbury earthquakes. In order to quantitatively model the basin response and improve the understanding of regional seismic hazard, geophysical methods will be used to develop shear wave velocity profiles across the Waikato basin. Active surface wave methods involve the deployment of linear arrays of geophones to record the surface waves generated by a sledge hammer. Passive surface wave methods involve the deployment of two-dimensional seismometer arrays to record ambient vibrations. At each site, the planned testing includes one active test and two to four passive arrays. The obtained data are processed to develop dispersion curves, which describe surface wave propagation velocity as a function of frequency (or wavelength). Dispersion curves are then inverted using the Geopsy software package to develop a suite of shear wave velocity profiles. Currently, more than ten sites in Waikato are under consideration for this project. This poster presents the preliminary results from the two sites that have been tested. The shear wave velocity profiles from all sites will be used to produce a 3D velocity model for the Waikato basin, a part of QuakeCoRE flagship programme 1.
Liquefaction-induced lateral spreading during earthquakes poses a significant hazard to the built environment, as observed in Christchurch during the 2010 to 2011 Canterbury Earthquake Sequence (CES). It is critical that geotechnical earthquake engineers are able to adequately predict both the spatial extent of lateral spreads and magnitudes of associated ground movements for design purposes. Published empirical and semi-empirical models for predicting lateral spread displacements have been shown to vary by a factor of <0.5 to >2 from those measured in parts of Christchurch during CES. Comprehensive post- CES lateral spreading studies have clearly indicated that the spatial distribution of the horizontal displacements and extent of lateral spreading along the Avon River in eastern Christchurch were strongly influenced by geologic, stratigraphic and topographic features.
This paper presents on-going challenges in the present paradigm shift of earthquakeinduced ground motion prediction from empirical to physics-based simulation methods. The 2010-2011 Canterbury and 2016 Kaikoura earthquakes are used to illustrate the predictive potential of the different methods. On-going efforts on simulation validation and theoretical developments are then presented, as well as the demands associated with the need for explicit consideration of modelling uncertainties. Finally, discussion is also given to the tools and databases needed for the efficient utilization of simulated ground motions both in specific engineering projects as well as for near-real-time impact assessment.
This is a joint Resilience Framework undertaken by the Electrical, Computer and Software Engineering Department of the University of Auckland in association with West Power and Orion networks and partially funded by the New Zealand National Science Challenge and QuakeCoRE. The Energy- Communication research group nearly accomplished two different researches focusing on both asset resilience and system resilience. Asset resilience research which covers underground cables system in Christchurch region is entitled “2010-2011 Canterbury Earthquake Sequence Impact on 11KV Underground Cables” and system resilience research which covers electricity distribution and communication system in West Coast region is entitled “NZ Electricity Distribution Network Resilience Assessment and Restoration Models following Major Natural Disturbance“. As the fourth milestone of the aforementioned research project, the latest outcome of both projects has been socialised with the stakeholders during the Cigre NZ 2019 Forum.
Knowing how to rapidly rebuild disaster-damaged infrastructure, while deciding appropriate recovery strategies and catering for future investment is a matter of core interest to government decision makers, utility providers, and business sectors. The purpose of this research is to explore the effects of decisions and outcomes for physical reconstruction on the overall recovery process of horizontal infrastructure in New Zealand using the Canterbury and Kaikoura earthquakes as cases. A mixed approach including a systematic review, questionnaire survey and semi-structured interviews is used to capture perspectives of those involved in reconstruction process and gain insights into the effect of critical elements on infrastructure downtime. Findings from this research will contribute towards advancements of a systems dynamics model considering critical decision-making variables across phases of the reconstruction process to assess how these variables affect the rebuild process and the corresponding downtime. This project will improve the ability to explore alternative resilience improvement pathways and test the efficacy of alternative means for facilitating a faster and better reconstruction process.
We present ground motion simulations of the Porters Pass (PP) fault in the Canterbury region of New Zealand; a major active source near Christchurch city. The active segment of the PP fault has an inferred length of 82 km and a mostly strike-slip sense of movement. The PP fault slip makes up approximately 10% of the total 37 mm/yr margin-parallel plate motion and also comprises a significant proportion of the total strain budget in regional tectonics. Given that the closest segment of the fault is less than 45 km from Christchurch city, the PP fault is crucial for accurate earthquake hazard assessment for this major population centre. We have employed the hybrid simulation methodology of Graves and Pitarka (2010, 2015), which combines low (f<1 Hz) and high (f>1 Hz) frequencies into a broadband spectrum. We have used validations from three moderate magnitude events (𝑀𝑤4.6 Sept 04, 2010; 𝑀𝑤4.6 Nov 06, 2010; 𝑀𝑤4.9 Apr 29, 2011) to build confidence for the 𝑀𝑤 > 7 PP simulations. Thus far, our simulations include multiple rupture scenarios which test the impacts of hypocentre location and the finite-fault stochastic rupture representation of the source itself. In particular, we have identified the need to use location-specific 1D 𝑉𝑠/𝑉𝑝 models for the high frequency part of the simulations to better match observations.
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.
Validation is an essential step to assess the applicability of simulated ground motions for utilization in engineering practice, and a comprehensive analysis should include both simple intensity measures (PGA, SA, etc), as well as the seismic response of a range of complex systems obtained by response history analysis. In order to enable a spectrum of complex structural systems to be considered in systematic validation of ground motion simulations in a routine fashion, an automated workflow was developed. Such a workflow enables validation of simulated ground motions in terms of different complex model responses by considering various ground motion sets and different ground motion simulation methods. The automated workflow converts the complex validation process into a routine one by providing a platform to perform the validation process promptly as a built-in process of simulation post-processing. As a case study, validation of simulated ground motions was investigated via the automated workflow by comparing the dynamic responses of three steel special moment frame (SMRF) subjected to the 40 observed and 40 simulated ground motions of 22 February 2011 Christchurch earthquake. The seismic responses of the structures are principally quantified via the peak floor acceleration and maximum inter-storey drift ratio. Overall, the results indicate a general agreement in seismic demands obtained using the recorded and simulated ensembles of ground motions and provide further evidence that simulated ground motions can be used in code-based structural performance assessments in-place of, or in combination with, ensembles of recorded ground motions.
Seismic isolation is an effective technology for significantly reducing damage to buildings and building contents. However, its application to light-frame wood buildings has so far been unable to overcome cost and technical barriers such as susceptibility to movement during high-wind loading. The precursor to research in the field of isolation of residential buildings was the 1994 Northridge Earthquake (6.7 MW) in the United States and the 1995 Kobe Earthquake (6.9 MW) in Japan. While only a small number of lives were lost in residential buildings in these events, the economic impact was significant with over half of earthquake recovery costs given to repair and reconstruction of residential building damage. A value case has been explored to highlight the benefits of seismically isolated residential buildings compared to a standard fixed-base dwellings for the Wellington region. Loss data generated by insurance claim information from the 2011 Christchurch Earthquake has been used by researchers to determine vulnerability functions for the current light-frame wood building stock. By further considering the loss attributed to drift and acceleration sensitive components, and a simplified single degree of freedom (SDOF) building model, a method for determining vulnerability functions for seismic isolated buildings was developed. Vulnerability functions were then applied directly in a loss assessment using the GNS developed software, RiskScape. Vulnerability was shown to dramatically reduce for isolated buildings compared to an equivalent fixed-base building and as a result, the monetary savings in a given earthquake scenario were significant. This work is expected to drive further interest for development of solutions for the seismic isolation of residential dwellings, of which one option is further considered and presented herein.
