This study provides an initial examination of source parameter uncertainty in a New Zealand ground motion simulation model, by simulating multiple event realisations with perturbed source parameters. Small magnitude events in Canterbury have been selected for this study due to the small number of source input parameters, the wealth of recorded data, and the lack of appreciable off-fault non-linear effects. Which provides greater opportunity to identify systematic source, path and site effects, required to robustly investigate the causes of uncertainty.
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.
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.
This paper investigates the effects of variability in source rupture parameters on site-specific physics-based simulated ground motions, ascertained through the systematic analysis of ground motion intensity measures. As a preliminary study, we consider simulations of the 22 February 2011 Christchurch earthquake using the Graves and Pitarka (2015) methodology. The effects of source variability are considered via a sensitivity study in which parameters (hypocentre location, earthquake magnitude, average rupture velocity, fault geometry and the Brune stress parameter) are individually varied by one standard deviation. The sensitivity of simulated ground motion intensity measures are subsequently compared against observational data. The preliminary results from this study indicate that uncertainty in the stress parameter and the rupture velocity have the most significant effect on the high frequency amplitudes. Conversely, magnitude uncertainty was found to be most influential on the spectral acceleration amplitudes at low frequencies. Further work is required to extend this preliminary study to exhaustively consider more events and to include parameter covariance. The ultimate results of this research will assist in the validation of the overall simulation method’s accuracy in capturing various rupture parameters, which is essential for the use of simulated ground motion models in probabilistic seismic hazard analysis.
Research indicates that aside from the disaster itself, the next major source of adverse outcomes during such events, is from errors by either the response leader or organisation. Yet, despite their frequency, challenge, complexity, and the risks involved; situations of extreme context remain one of the least researched areas in the leadership field. This is perhaps surprising. In the 2010 and 2011 (Christchurch) earthquakes alone, 185 people died and rebuild costs are estimated to have been $40b. Add to this the damage and losses annually around the globe arising from natural disasters, major business catastrophes, and military conflict; there is certainly a lot at stake (lives, way of life, and our well-being). While over the years, much has been written on leadership, there is a much smaller subset of articles on leadership in extreme contexts, with the majority of these focusing on the event rather than leadership itself. Where leadership has been the focus, the spotlight has shone on the actions and capabilities of one person - the leader. Leadership, however, is not simply one person, it is a chain or network of people, delivering outcomes with the support of others, guided by a governance structure, contextualised by the environment, and operating on a continuum across time (before, during, and after an event). This particular research is intended to examine the following: • What are the leadership capabilities and systems necessary to deliver more successful outcomes during situations of extreme context; • How does leadership in these circumstances differ from leadership during business as usual conditions; • Lastly, through effective leadership, can we leverage these unfortunate events to thrive, rather than merely survive?
Courage has remained an elusive concept to define despite having been in the English lexicon for hundreds of years. The Canterbury earthquake sequence that began in 2010 provided a unique context in which to undertake research that would contribute to further conceptualisation of courage. This qualitative study was undertaken in Christchurch, New Zealand, with adults over the age of 70 who experienced the Canterbury earthquakes and continued to live in the Canterbury region. The population group was chosen because it is an under researched group in post-disaster environments, and one that offers valuable insights because of members' length and breadth of life experiences, and likely reminiscent and reflective life stage. A constructivist grounded theory approach was utilised, with data collected through semi-structured focus groups and individual key informant interviews. The common adverse experience of the participants initially discussed was the earthquakes, which was followed by exploration of courage in their other lived experiences. Through an inductive process of data analysis, conceptual categories were identified, which when further analysed and integrated, contributed to a definition of courage. The definition was subsequently discussed with social work professionals who had remained working in the Canterbury region after experiencing the earthquakes. From the examples and the actions described within these, a process model was developed to support the application of courage. The model includes five steps: recognising an adverse situation, making a conscious decision to act, accessing sources of motivation, mastering emotion and taking action. Defining and utilising courage can help people to face adversity associated with everyday life and ultimately supports self-actualisation and self-development. Recommendations from the study include teaching about courage within social work education, utilising the process model within supervision, intentionally involving older adults in emergency management planning and developing specific social work tasks in hospital settings following a disaster.
Probabilistic Structural Fire Engineering (PSFE) has been introduced to overcome the limitations of current conventional approaches used for the design of fire-exposed structures. Current structural fire design investigates worst-case fire scenarios and include multiple thermal and structural analyses. PSFE permits buildings to be designed to a level of life safety or economic loss that may occur in future fire events with the help of a probabilistic approach. This thesis presents modifications to the adoption of a Performance-Based Earthquake Engineering (PBEE) framework in Probabilistic Structural Fire Engineering (PSFE). The probabilistic approach runs through a series of interrelationships between different variables, and successive convolution integrals of these interrelationships result in probabilities of different measures. The process starts with the definition of a fire severity measure (FSM), which best relates fire hazard intensity with structural response. It is identified by satisfying efficiency and sufficiency criteria as described by the PBEE framework. The relationship between a fire hazard and corresponding structural response is established by analysis methods. One method that has been used to quantify this relationship in PSFE is Incremental Fire Analysis (IFA). The existing IFA approach produces unrealistic fire scenarios, as fire profiles may be scaled to wide ranges of fire severity levels, which may not physically represent any real fires. Two new techniques are introduced in this thesis to limit extensive scaling. In order to obtain an annual rate of exceedance of fire hazard and structural response for an office building, an occurrence model and an attenuation model for office fires are generated for both Christchurch city and New Zealand. The results show that Christchurch city is 15% less likely to experience fires that have the potential to cause structural failures in comparison to all of New Zealand. In establishing better predictive relationships between fires and structural response, cumulative incident radiation (a fire hazard property) is found to be the most appropriate fire severity measure. This research brings together existing research on various sources of uncertainty in probabilistic structural fire engineering, such as elements affecting post-flashover fire development factors (fuel load, ventilation, surface lining and compartment geometry), fire models, analysis methods and structural reliability. Epistemic uncertainty and aleatory uncertainty are investigated in the thesis by examining the uncertainty associated with modelling and the factors that influence post-flashover development of fires. A survey of 12 buildings in Christchurch in combination with recent surveys in New Zealand produced new statistical data on post-flashover development factors in office buildings in New Zealand. The effects of these parameters on temperature-time profiles are evaluated. The effects of epistemic uncertainty due to fire models in the estimation of structural response is also calculated. Parametric fires are found to have large uncertainty in the prediction of post-flashover fires, while the BFD curves have large uncertainties in prediction of structural response. These uncertainties need to be incorporated into failure probability calculations. Uncertainty in structural modelling shows that the choices that are made during modelling have a large influence on realistic predictions of structural response.
This dissertation addresses a diverse range of topics in the physics-based broadband ground motion simulation, with a focus on New Zealand applications. In particular the following topics are addressed: the methodology and computational implementation of a New Zealand Velocity Model for broadband ground motion simulation; generalised parametric functions and spatial correlations for seismic velocities in the Canterbury, New Zealand region from surface-wave-based site characterisation; and ground motion simulations of Hope Fault earthquakes. The paragraphs below outline each contribution in more detail. A necessary component in physics-based ground motion simulation is a 3D model which details the seismic velocities in the region of interest. Here a velocity model construction methodology, its computational implementation, and application in the construction of a New Zealand velocity model for use in physics-based broadband ground motion simulation are presented. The methodology utilises multiple datasets spanning different length scales, which is enabled via the use of modular sub-regions, geologic surfaces, and parametric representations of crustal velocity. A number of efficiency-related workflows to decrease the overall computational construction time are employed, while maintaining the flexibility and extensibility to incorporate additional datasets and re- fined velocity parameterizations as they become available. The model comprises explicit representations of the Canterbury, Wellington, Nelson-Tasman, Kaikoura, Marlborough, Waiau, Hanmer and Cheviot sedimentary basins embedded within a regional travel-time tomography-based velocity model for the shallow crust and provides the means to conduct ground motion simulations throughout New Zealand for the first time. Recently developed deep shear-wave velocity profiles in Canterbury enabled models that better characterise the velocity structure within geologic layers of the Canterbury sedimentary basin to be developed. Here the development of depth- and Vs30-dependent para-metric velocity and spatial correlation models to characterise shear-wave velocities within the geologic layers of the Canterbury sedimentary basin are presented. The models utilise data from 22 shear-wave velocity profiles of up to 2.5km depth (derived from surface wave analysis) juxtaposed with models which detail the three-dimensional structure of the geologic formations in the Canterbury sedimentary basin. Parametric velocity equations are presented for Fine Grained Sediments, Gravels, and Tertiary layer groupings. Spatial correlations were developed and applied to generate three-dimensional stochastic velocity perturbations. Collectively, these models enable seismic velocities to be realistically represented for applications such as 3D ground motion and site response simulations. Lastly the New Zealand velocity model is applied to simulate ground motions for a Mw7.51 rupture of the Hope Fault using a physics-based simulation methodology and a 3D crustal velocity model of New Zealand. The simulation methodology was validated for use in the region through comparison with observations for a suite of historic small magnitude earthquakes located proximal to the Hope Fault. Simulations are compared with conventionally utilised empirical ground motion models, with simulated peak ground velocities being notably higher in regions with modelled sedimentary basins. A sensitivity analysis was undertaken where the source characteristics of magnitude, stress parameter, hypocentre location and kinematic slip distribution were varied and an analysis of their effect on ground motion intensities is presented. It was found that the magnitude and stress parameter strongly influenced long and short period ground motion amplitudes, respectively. Ground motion intensities for the Hope Fault scenario are compared with the 2016 Kaikoura Mw7.8 earthquake, it was found that the Kaikoura earthquake produced stronger motions along the eastern South Island, while the Hope Fault scenario resulted in stronger motions immediately West of the near-fault region. The simulated ground motions for this scenario complement prior empirically-based estimates and are informative for mitigation and emergency planning purposes.
