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.
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.
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.
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.
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.
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.
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.
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.
