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.
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.
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.
Background and methodology The Mw 7.8, 14th November 2016 earthquake centred (item b, figure 1) in the Hurunui District of the South Island, New Zealand, damaged critical infrastructure across North Canterbury and Marlborough. We investigate the impacts to infrastructure and adaptations to the resulting service disruption in four small rural towns (figure 1): Culverden (a), Waiau (c), Ward (d) and Seddon (e). This is accomplished though literary research, interviews and geospatial analysis. Illustrating our methods, we have displayed here a Hurunui District hazard map (figure 2b) and select infrastructure inventories (figures 2a, 3).
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).
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.
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.
