Hybrid broadband simulation methods typically compute high-frequency portion of ground-motions using a simplified-physics approach (commonly known as “stochastic method”) using the same 1D velocity profile, anelastic attenuation profile and site-attenuation (κ0) value for all sites. However, these parameters relating to Earth structure are known to vary spatially. In this study we modify this conventional approach for high-frequency ground-shaking by using site-specific input parameters (referred to as “site-specific”) and analyze improvements over using same parameters for all sites (referred to as “generic”). First, we theoretically understand how different 1D velocity profiles, anelastic attenuation profiles and site-attenuation (κ0) values affects the Fourier Acceleration Spectrum (FAS). Then, we apply site-specific method to simulate 10 events from the 2010-2011 Canterbury earthquake sequence to assess performance against the generic approach in predicting recorded ground-motions. Our initial results suggest that the site-specific method yields a lower simulation standard deviation than generic case.
During the 2011 M7.8 Kaikōura earthquake, ground motions recorded near the epicentre showed a significant spatial variation. The Te Mara farm (WTMC) station, the nearest to the epicentre, recorded 1g and 2.7g of horizontal and vertical peak ground accelerations (PGA), respectively. The nearby Waiu Gorge (WIGC) station recorded a horizontal PGA of 0.8g. Interestingly, however, the Culverden Airlie Farm (CULC) station that was very closely located to WIGC recorded a horizontal PGA of only 0.25g. This poster demonstrates how the local geological condition could have contributed to the spatially variable ground motions observed in the North Canterbury, based on the results of recently conducted geophysical investigations. The surficial geology of this area is dominated by alluvial gravel deposits with traces of silt. A borehole log showed that the thickness of the sediments at WTMC is over 76 metres. Interestingly, the shear wave velocity (Vs) profiles obtained from the three strong motion sites suggest unusually high shear wave velocity of the gravelly sediments. The velocity of sediments and the lack of clear peaks in the horizontal-to-vertical (H/V) spectral ratio at WTMC suggest that the large ground motion observed at this station was likely caused by the proximity of the station to the causative fault itself; the site effect was likely insignificant. Comparisons of H/V spectral ratios and Vs profiles suggest that the sediment thickness is much smaller at WIGC compared with CULC; the high PGA at WIGC was likely influenced by the high-frequency amplification caused by the response of shallow sediments.
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.
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.
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.
Heathcote Valley school strong motion station (HVSC) consistently recorded ground motions with higher intensities than nearby stations during the 2010-2011 Canterbury earthquakes. For example, as shown in Figure 1, for the 22 February 2011 Christchurch earthquake, peak ground acceleration at HVSC reached 1.4 g (horizontal) and 2 g (vertical), the largest ever recorded in New Zealand. Strong amplification of ground motions is expected at Heathcote Valley due to: 1) the high impedance contrast at the soil-rock interface, and 2) the interference of incident and surface waves within the valley. However, both conventional empirical ground motion prediction equations (GMPE) and the physics-based large scale ground motions simulations (with empirical site response) are ineffective in predicting such amplification due to their respective inherent limitations.
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.
Overview of SeisFinder SeisFinder is an open-source web service developed by QuakeCoRE and the University of Canterbury, focused on enabling the extraction of output data from computationally intensive earthquake resilience calculations. Currently, SeisFinder allows users to select historical or future events and retrieve ground motion simulation outputs for requested geographical locations. This data can be used as input for other resilience calculations, such as dynamic response history analysis. SeisFinder was developed using Django, a high-level python web framework, and uses a postgreSQL database. Because our large-scale computationally-intensive numerical ground motion simulations produce big data, the actual data is stored in file systems, while the metadata is stored in the database. The basic SeisFinder architecture is shown in Figure 1.
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.
On 15 August 1868, a great earthquake struck off the coast of the Chile-Peru border generating a tsunami that travelled across the Pacific. Wharekauri-Rekohu-Chatham Islands, located 800 km east of Christchurch, Aotearoa-New Zealand (A-NZ) was one of the worst affected locations in A-NZ. Tsunami waves, including three over 6 metres high, injured and killed people, destroyed buildings and infrastructure, and impacted the environment, economy and communities. While experience of disasters, and advancements in disaster risk reduction systems and technology have all significantly advanced A-NZ’s capacity to be ready for and respond to future earthquakes and tsunami, social memory of this event and other tsunamis during our history has diminished. In 2018, a team of scientists, emergency managers and communication specialists collaborated to organise a memorial event on the Chatham Islands and co-ordinate a multi-agency media campaign to commemorate the 150th anniversary of the 1868 Arica tsunami. The purpose was to raise awareness of the disaster and to encourage preparedness for future tsunami. Press releases and science stories were distributed widely by different media outlets and many attended the memorial event indicating public interest for commemorating historical disasters. We highlight the importance of commemorating disaster anniversaries through memorial events, to raise awareness of historical disasters and increase community preparedness for future events – “lest we forget and let us learn.”
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).
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.
The purpose of this research is to investigate men’s experiences of the 2016 7.8 magnitude Kaikōura earthquake and Tsunami. While, research into the impacts of the earthquake has been conducted, few studies have examined how gender shaped people’s experiences of this natural hazard event. Analysing disasters through a gender lens has significantly contributed to disaster scholarship in identifying the resilience and vulnerabilities of individuals and communities pre- and post-disaster (Fordham, 2012; Bradshaw, 2013). This research employs understandings of masculinities (Connell, 2005), to examine men’s strengths and challenges in responding, recovering, and coping following the earthquake. Qualitative inquiry was carried out in Northern Canterbury and Marlborough involving 18 face-to-face interviews with men who were impacted by the Kaikōura earthquake and its aftermath. Interview material is being analysed using thematic and narrative analysis. Some of the preliminary findings have shown that men took on voluntary roles in addition to their fulltime paid work resulting in long hours, poor sleep and little time spent with family. Some men assisted wives and children to high ground then drove into the tsunami zone to check on relatives or to help evacuate people. Although analysis of the findings is currently ongoing, preliminary findings have identified that the men who participated in the study have been negatively impacted by the 2016 Kaikōura earthquake. A theme identified amongst participants was an avoidance to seek support with the challenges they were experiencing due to the earthquake. The research findings align with key characteristics of masculinity, including demonstrating risky behaviours and neglecting self or professional care. This study suggests that these behaviours affect men’s overall resilience, and thus the resilience of the wider community.
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.
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.
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.
The purpose of this research is to investigate men’s experiences of the 2016 7.8 magnitude Kaikōura earthquake and Tsunami. While, research into the impacts of the earthquake has been conducted, few studies have examined how gender shaped people’s experiences of this natural hazard event. Analysing disasters through a gender lens has significantly contributed to disaster scholarship in identifying the resilience and vulnerabilities of individuals and communities pre- and post-disaster (Fordham, 2012; Bradshaw, 2013). This research employs understandings of masculinities (Connell, 2005), to examine men’s strengths and challenges in responding, recovering, and coping following the earthquake. Qualitative inquiry was carried out in Northern Canterbury and Marlborough involving 18 face-to-face interviews with men who were impacted by the Kaikōura earthquake and its aftermath. Interview material is being analysed using thematic and narrative analysis. Some of the preliminary findings have shown that men took on voluntary roles in addition to their fulltime paid work resulting in long hours, poor sleep and little time spent with family. Some men assisted wives and children to high ground then drove into the tsunami zone to check on relatives or to help evacuate people. Although analysis of the findings is currently ongoing, preliminary findings have identified that the men who participated in the study have been negatively impacted by the 2016 Kaikōura earthquake. A theme identified amongst participants was an avoidance to seek support with the challenges they were experiencing due to the earthquake. The research findings align with key characteristics of masculinity, including demonstrating risky behaviours and neglecting self or professional care. This study suggests that these behaviours affect men’s overall resilience, and thus the resilience of the wider community.
