This paper presents a methodology by which both site-specific and spatially distributed ground motion intensity can be obtained immediately following an earthquake event. The methodology makes use of both prediction models for ground motion intensity and its correlation over spatial distances. A key benefit of the methodology is that the ground motion intensity at a given location is not a single value but a distribution of values. The distribution is comprised of both a mean and also standard deviation, with the standard deviation being a function of the distance to nearby strong motion stations. The methodology is illustrated for two applications. Firstly, maps of conditional peak ground acceleration (PGA) have been developed for the major events in the Canterbury earthquake sequence. It is illustrated how these conditional maps can be used for post-event evaluation of liquefaction triggering criteria which have been adopted by the Department of Building and Housing (DBH). Secondly, the conditional distribution of response spectral ordinates is obtained at a specific location for the purposes of determining appropriate ground motion records for use in seismic response analyses of important structures at locations where direct recordings are absent.
This manuscript provides a critical examination of the ground motions recorded in the near-source region resulting from the 22 February 2011 Christchurch earthquake. Particular attention is given to reconciling the observed spatial distribution of ground motions in terms of physical phenomena related to source, path and site effects. The large number of near-source observed strong ground motions show clear evidence of: forward-directivity, basin generated surface waves, liquefaction and other significant nonlinear site response. The pseudo-acceleration response spectra (SA) amplitudes and significant duration of strong motions agree well with empirical prediction models, except at long vibration periods where the influence of basin-generated surface waves and nonlinear site response are significant and not adequately accounted for in empirical SA models. Pseudo-acceleration response spectra are also compared with those observed in the 4 September 2010 Darfield earthquake and routine design response spectra used in order to emphasise the amplitude of ground shaking and elucidate the importance of local geotechnical characteristics on surface ground motions. The characteristics of the observed vertical component accelerations are shown to be strongly dependent on source-to-site distance and are comparable with those from the 4 September 2010 Darfield earthquake, implying the large amplitudes observed are simply a result of many observations at close distances rather than a peculiar source effect.
Located on the edge of two tectonic plates, New Zealand has numerous fault lines and seismic risk across the whole country. The way this risk is communicated affects whether people prepare effectively or at all. Research has shown that perceptions of risk are affected by slight changes in wording, and that probabilities commonly reported by experts and media are often interpreted subjectively based on context. In the context of volcanoes, research has found that given a certain probability of a volcano in a specific time window, people perceive risk as higher in later time intervals within that window. The present study examines this pattern with regard to earthquakes and aftershocks in the New Zealand context. Participants in both Wellington (N = 102) and Christchurch (N = 98) were presented an expert statement of earthquake risk within a given time window in Wellington and aftershock risk in Christchurch, and asked to rate their perception of risk in specific intervals across the time window. For a Wellington earthquake, participants perceived risk as incrementally higher toward the end of the 50 year time window whereas for a Christchurch aftershock, risk perception increased slightly for the first three intervals of the 12 month time window. Likelihood of preparing was constant over the time windows, with Wellington citizens rating themselves more likely than Christchurch citizens to prepare for either an earthquake or aftershock, irrespective of current level of preparedness. These findings suggest that people view earthquakes as more likely later toward the end of a given time window and that they view aftershocks very differently to scientific predictions.
© 2019, Springer-Verlag GmbH Germany, part of Springer Nature. Prediction of building collapse due to significant seismic motion is a principle objective of earthquake engineers, particularly after a major seismic event when the structure is damaged and decisions may need to be made rapidly concerning the safe occupation of a building or surrounding areas. Traditional model-based pushover analyses are effective, but only if the structural properties are well understood, which is not the case after an event when that information is most useful. This paper combines hysteresis loop analysis (HLA) structural health monitoring (SHM) and incremental dynamic analysis (IDA) methods to identify and then analyse collapse capacity and the probability of collapse for a specific structure, at any time, a range of earthquake excitations to ensure robustness. This nonlinear dynamic analysis enables constant updating of building performance predictions following a given and subsequent earthquake events, which can result in difficult to identify deterioration of structural components and their resulting capacity, all of which is far more difficult using static pushover analysis. The combined methods and analysis provide near real-time updating of the collapse fragility curves as events progress, thus quantifying the change of collapse probability or seismic induced losses very soon after an earthquake for decision-making. Thus, this combination of methods enables a novel, higher-resolution analysis of risk that was not previously available. The methods are not computationally expensive and there is no requirement for a validated numerical model, thus providing a relatively simpler means of assessing collapse probability immediately post-event when such speed can provide better information for critical decision-making. Finally, the results also show a clear need to extend the area of SHM toward creating improved predictive models for analysis of subsequent events, where the Christchurch series of 2010–2011 had significant post-event aftershocks.
The rapid classification of building damage states or placards after an earthquake is vital for enabling an efficient emergency response and informed decision-making for rehabilitation and recovery purposes. Traditional methods rely heavily on inspector-led on-site surveys, which are often time-consuming, resource-intensive, and susceptible to human error. This study introduces a machine learning-supported surrogate model designed to streamline the assessment of building damage, focusing on the automated assignment of damage placards within the context of New Zealand's post-earthquake evaluation frameworks. The study evaluates two key safety evaluation protocols—Rapid Building Assessment (RBA) and Detailed Damage Evaluation (DDE)—and integrates corresponding databases derived from the 2010–2011 Canterbury Earthquake Sequence (CES) in Christchurch. Six ML classifiers—Multilayer Perceptron (MLP), Random Forest (RF), Support Vector Machine (SVM), K-Nearest Neighbours (KNN), Gradient Boosting Classifier (GBC), and Gradient Bagging (GBag)—were rigorously tested across both databases. The results indicate that the RF-based surrogate model outperforms the other classifiers across both RBA and DDE protocols. Two distinct sets of critical predictors have been further identified for each protocol, allowing for the rapid retrieval of essential data for future on-site surveys, while retaining the RF model's predictive accuracy. The developed surrogate model provides a pragmatic tool for practising engineers to rapidly assign placards to damaged structures and for policymakers and building owners to make informed recovery decisions for earthquake-affected buildings.
The performance of conventionally designed reinforced concrete (RC) structures during the 2011 Christchurch earthquake has demonstrated that there is greater uncertainty in the seismic performance of RC components than previously understood. RC frame and wall structures in the Christchurch central business district were observed to form undesirable cracks patterns in the plastic hinge region while yield penetration either side of cracks, and into development zones, were less than theoretical predictions. The implications of this unexpected behaviour: (i) significantly less available ductility; (ii) less hysteretic energy dissipation; and (iii) the localization of peak reinforcement strains, results in considerable doubt for the residual capacity of RC structures. The significance of these consequences has prompted a review of potential sources of uncertainty in seismic experimentation with the intention to improve the current confidence level for newly designed conventional RC structures. This paper attempts to revisit the principles of RC mechanics, in particular, to consider the influence of loading history, concrete tensile strength, and reinforcement ratio on the performance of ‘real’ RC structures compared to experimental test specimens.
In practice, several competing liquefaction evaluation procedures (LEPs) are used to compute factors of safety against soil liquefaction, often for use within a liquefaction potential index (LPI) framework to assess liquefaction hazard. At present, the influence of the selected LEP on the accuracy of LPI hazard assessment is unknown, and the need for LEP-specific calibrations of the LPI hazard scale has never been thoroughly investigated. Therefore, the aim of this study is to assess the efficacy of three CPT-based LEPs from the literature, operating within the LPI framework, for predicting the severity of liquefaction manifestation. Utilising more than 7000 liquefaction case studies from the 2010–2011 Canterbury (NZ) earthquake sequence, this study found that: (a) the relationship between liquefaction manifestation severity and computed LPI values is LEP-specific; (b) using a calibrated, LEP-specific hazard scale, the performance of the LPI models is essentially equivalent; and (c) the existing LPI framework has inherent limitations, resulting in inconsistent severity predictions against field observations for certain soil profiles, regardless of which LEP is used. It is unlikely that revisions of the LEPs will completely resolve these erroneous assessments. Rather, a revised index which more adequately accounts for the mechanics of liquefaction manifestation is needed.
In the last two decades, New Zealand (NZ) has experienced significant earthquakes, including the 2010 M 7.2 Darfield, 2011 M 6.2 Christchurch, and 2016 M 7.8 Kaikōura events. Amongst these large events, tens of thousands of smaller earthquakes have occurred. While previous event and ground-motion databases have analyzed these events, many events below M 4 have gone undetected. The goal of this study is to expand on previous databases, particularly for small magnitude (M<4) and low-amplitude ground motions. This new database enables a greater understanding of regional variations within NZ and contributes to the validity of internationally developed ground-motion models. The database includes event locations and magnitude estimates with uncertainty considerations, and tectonic type assessed in a hierarchical manner. Ground motions are extracted from the GeoNet FDSN server and assessed for quality using a neural network classification approach. A deep neural network approach is also utilized for picking P and S phases for determination of event hypocentres. Relative hypocentres are further improved by double-difference relocation and will contribute toward developing shallow (< 50 km) seismic tomography models. Analysis of the resulting database is compared with previous studies for discussion of implications toward national hazard prediction models.
Having a quick but reliable insight into the likelihood of damage to bridges immediately after an earthquake is an important concern especially in the earthquake prone countries such as New Zealand for ensuring emergency transportation network operations. A set of primary indicators necessary to perform damage likelihood assessment are ground motion parameters such as peak ground acceleration (PGA) at each bridge site. Organizations, such as GNS in New Zealand, record these parameters using distributed arrays of sensors. The challenge is that those sensors are not installed at, or close to, bridge sites and so bridge site specific data are not readily available. This study proposes a method to predict ground motion parameters for each bridge site based on remote seismic array recordings. Because of the existing abundant source of data related to two recent strong earthquakes that occurred in 2010 and 2011 and their aftershocks, the city of Christchurch is considered to develop and examine the method. Artificial neural networks have been considered for this research. Accelerations recorded by the GeoNet seismic array were considered to develop a functional relationship enabling the prediction of PGAs. http://www.nzsee.org.nz/db/2013/Posters.htm
This paper concerns the explicit consideration of near-fault directivity in conventional ground motion prediction models, and its implication for probabilistic seismic hazard analysis (PSHA) in New Zealand. The proposed approach utilises recently developed models by Shahi & Baker (2011), which account for both the 'narrowband' nature of the directivity pulse on spectral ordinates, and the probability of pulse occurrence at the site of interest. Furthermore, in order to correctly consider directivity, distributed seismicity sources are considered as finite-faults, as opposed to their (incorrect) conventional treatment as point-sources. The significance of directivity on hazard analysis results is illustrated for various vibration periods at generic sites located in Christchurch and Otira, two locations whose seismic hazard is comprised of notably different seismic sources. When compared to the PSHA results considering directivity and distributed seismicity as finite faults, it is shown that the NZS1170.5:2004 directivity factor is notably unconservative for all vibration periods in Otira (i.e. high seismic hazard region); and unconservative for Christchurch at short-to-moderate vibration periods ( < 3s); but conservative at long periods ( > 4s).
Currently there is a worldwide renaissance in timber building design. At the University of Canterbury, new structural systems for commercial multistorey timber buildings have been under development since 2005. These systems incorporate large timber sections connected by high strength post-tensioning tendons, and timber-concrete composite floor systems, and aim to compete with existing structural systems in terms of cost, constructability, operational and seismic performance. The development of post-tensioned timber systems has created a need for improved lateral force design approaches for timber buildings. Current code provisions for seismic design are based on the strength of the structure, and do not adequately account for its deformation. Because timber buildings are often governed by deflection, rather than strength, this can lead to the exceedence of design displacement limitations imposed by New Zealand codes. Therefore, accurate modeling approaches which define both the strength and deformation of post-tensioned timber buildings are required. Furthermore, experimental testing is required to verify the accuracy of these models. This thesis focuses on the development and experimental verification of modeling approaches for the lateral force design of post-tensioned timber frame and wall buildings. The experimentation consisted of uni-direcitonal and bi-directional quasi-static earthquake simulation on a two-thirds scale, two-storey post-tensioned timber frame and wall building with timber-concrete composite floors. The building was subjected to lateral drifts of up to 3% and demonstrated excellent seismic performance, exhibiting little damage. The building was instrumented and analyzed, providing data for the calibration of analytical and numerical models. Analytical and numerical models were developed for frame, wall and floor systems that account for significant deformation components. The models predicted the strength of the structural systems for a given design performance level. The static responses predicted by the models were compared with both experimental data and finite element models to evaluate their accuracy. The frame, wall and floor models were then incorporated into an existing lateral force design procedure known as displacement-based design and used to design several frame and wall structural systems. Predictions of key engineering demand parameters, such as displacement, drift, interstorey shear, interstorey moment and floor accelerations, were compared with the results of dynamic time-history analysis. It was concluded that the numerical and analytical models, presented in this thesis, are a sound basis for determining the lateral response of post-tensioned timber buildings. However, future research is required to further verify and improve these prediction models.
Liquefaction-induced lateral spreading in large seismic events often results in pervasive and costly damage to engineering structures and lifelines, making it a critical component of engineering design. However, the complex nature of this phenomenon leads to designing for such a hazard extremely challenging and there is a clear for an improved understanding and predicting liquefaction-induced lateral spreading. The 2010-2011 Canterbury (New Zealand) Earthquakes triggered severe liquefaction-induced lateral spreading along the streams and rivers of the Christchurch region, causing extensive damage to roads, bridges, lifelines, and structures in the vicinity. The unfortunate devastation induced from lateral spreading in these events also rendered the rare opportunity to gain an improved understanding of lateral spreading displacements specific to the Christchurch region. As part of this thesis, the method of ground surveying was employed following the 4 September 2010 Darfield (Mw 7.1) and 22 February 2011 Christchurch (Mw 6.2) earthquakes at 126 locations (19 repeated) throughout Christchurch and surrounding suburbs. The method involved measurements and then summation of crack widths along a specific alignment (transect) running approximately perpendicular to the waterway to indicate typically a maximum lateral displacement at the bank and reduction of the magnitude of displacements with distance from the river. Rigorous data processing and comparisons with alternative measurements of lateral spreading were performed to verify results from field observations and validate the method of ground surveying employed, as well as highlight the complex nature of lateral spreading displacements. The welldocumented field data was scrutinized to gain an understanding of typical magnitudes and distribution patterns (distribution of displacement with distance) of lateral spreading observed in the Christchurch area. Maximum displacements ranging from less than 10 cm to over 3.5 m were encountered at the sites surveyed and the area affected by spreading ranged from less than 20 m to over 200 m from the river. Despite the highly non-uniform displacements, four characteristic distribution patterns including large, distributed ground displacements, block-type movements, large and localized ground displacements, and areas of little to no displacements were identified. Available geotechnical, seismic, and topographic data were collated at the ground surveying sites for subsequent analysis of field measurements. Two widely-used empirical models (Zhang et al. (2004), Youd et al. (2002)) were scrutinized and applied to locations in the vicinity of field measurements for comparison with model predictions. The results indicated generally poor correlation (outside a factor of two) with empirical predictions at most locations and further validated the need for an improved, analysis- based method of predicting lateral displacements that considers the many factors involved on a site-specific basis. In addition, the development of appropriate model input parameters for the Youd et al. (2002) model led to a site-specific correlation of soil behavior type index, Ic, and fines content, FC, for sites along the Avon River in Christchurch that matched up well with existing Ic – FC relationships commonly used in current practice. Lastly, a rigorous analysis was performed for 25 selected locations of ground surveying measurements along the Avon River where ground slope conditions are mild (-1 to 2%) and channel heights range from about 2 – 4.5 m. The field data was divided into categories based on the observed distribution pattern of ground displacements including: large and distributed, moderate and distributed, small to negligible, and large and localized. A systematic approach was applied to determine potential critical layers contributing to the observed displacement patterns which led to the development of characteristic profiles for each category considered. The results of these analyses outline an alternative approach to the evaluation of lateral spreading in which a detailed geotechnical analysis is used to identify the potential for large spreading displacements and likely spatial distribution patterns of spreading. Key factors affecting the observed magnitude and distribution of spreading included the thickness of the critical layer, relative density, soil type and layer continuity. It was found that the large and distributed ground displacements were associated with a thick (1.5 – 2.5 m) deposit of loose, fine to silty sand (qc1 ~4-7 MPa, Ic 1.9-2.1, qc1n_cs ~50-70) that was continuous along the bank and with distance from the river. In contrast, small to negligible displacements were characterized by an absence of or relatively thin (< 1 m), discontinuous critical layer. Characteristic features of the moderate and distributed displacements were found to be somewhere between these two extremes. The localized and large displacements showed a characteristic critical layer similar to that observed in the large and distributed sites but that was not continuous and hence leading to the localized zone of displacement. The findings presented in this thesis illustrate the highly complex nature of lateral displacements that cannot be captured in simplified models but require a robust geotechnical analysis similar to that performed for this research.
Research Report No.2010-03Ground motion prediction equations (GMPEs) for geometric-mean pseudo-spectral acceleration amplitudes from New Zealand (NZ) earthquakes are developed. A database of 2437 three-component ground motion records is developed by applying stringent quality criteria to the historically recorded events in NZ. Despite the large number of records, the database is deficient in empirical records from large magnitude events recorded at close distances to the fault rupture plane. As a result, the basis for the NZ-specific GMPE development is to examine the applicability of foreign GMPEs for similar tectonic regions and then modify the most applicable GMPEs based on both theoretical and statistically significant empirically-driven arguments. For active shallow crustal events, five different GMPEs are considered. It was found that the McVerry et al. (2006) model, which is the current model upon which seismic design guidelines and site-specific seismic hazard analyses in NZ are based, provided the worst fit to the NZ database, and that the Chiou et al. (2010) (C10) modification of the Chiou and Youngs (2008) model was the most applicable. Discrepancies between the C10 model and the NZ database that were empirically identified and theoretically justified were used to modify the C10 model for: (i) small magnitude scaling; (ii) scaling of short period ground motion from normal faulting events in volcanic crust; (iii) scaling of ground motions on very hard rock sites; (iv) anelastic attenuation in the NZ crust; and (v) consideration of the increased anelastic attenuation in the Taupo Volcanic Zone (TVZ). For subduction slab events, initially three models were considered. It was found that all of the models had some significant biases with respect to applicability for NZ. The Zhao et al. (2006) (Z06) model was selected because of the rigorous database upon which it was developed and modified by: (i) NZ-specific scaling at small magnitudes; (ii) path scaling at large distances; (iii) consideration of the increased TVZ attenuation; and (iv) revision of the standard deviation model. Based on these modifications the developed model showed no bias of the inter- and intra-event residuals as a function of various predictor variables. The standard deviation of the residuals using the revised standard deviation model also indicated that the model has an adequate precision. Three GMPEs were considered for subduction interface events. The Zhao et al. (2006) (Z06) model was the best performing model with only bias exhibited in the site response model, and possible over-prediction of large magnitude events. The Z06 interface model was modified to account for site response and magnitude scaling using the same functional forms as those of the developed active shallow crustal and subduction slab models. The developed model showed no bias of the inter- and intra-event residuals as a function of various predictor variables. The developed GMPEs include specific features as evident in the NZ database; consistent scaling for parameters not well constrained by the NZ database; and pseudo-spectral amplitudes for vibration periods from 0.01 to 10 seconds. Hence, these models represent a significant advance in the state-of-the art for empirical ground motion prediction in NZ.
Predicting building collapse due to seismic motion is critical in design and more so after a major event. Damaged structures can appear sound, but collapse under following major events. There can thus be significant risk in decision making after a major seismic event concerning the safe occupation of a building or surrounding areas, versus the unknown impact of unknown major aftershocks. Model-based pushover analyses are effective if the structural properties are well understood, which is not valid post-event when this risk information is most useful. This research combines Hysteresis Loop Analysis (HLA) structural health monitoring (SHM) and Incremental Dynamic Analysis (IDA) methods to determine collapse capacity and probability of collapse for a specific structure, at any time, a range of earthquake excitations to ensure robustness. The nonlinear dynamic analysis method presented enables constant updating of building performance predictions using post-event SHM results. The resulting combined methods provide near real-time updating of collapse fragility curves as events progress, quantifying the change of collapse probability or seismic induced losses for decision-making - a novel, higher resolution risk analysis than previously available. The methods are not computationally expensive and there is no requirement for a validated numerical model. Results show significant potential benefits and a clear evolution of risk. They also show clear need for extending SHM toward creating improved predictive models for analysis of subsequent events, where the Christchurch series of 2010-2011 had significant post-event aftershocks after each main event. Finally, the overall method is generalisable to any typical engineering demand parameter.
This article presents a quantitative case study on the site amplification effect observed at Heathcote Valley, New Zealand, during the 2010-2011 Canterbury earthquake sequence for 10 events that produced notable ground acceleration amplitudes up to 1.4g and 2.2g in the horizontal and vertical directions, respectively. We performed finite element analyses of the dynamic response of the valley, accounting for the realistic basin geometry and the soil non-linear response. The site-specific simulations performed significantly better than both empirical ground motion models and physics based regional-scale ground motion simulations (which empirically accounts for the site effects), reducing the spectral acceleration prediction bias by a factor of two in short vibration periods. However, our validation exercise demonstrated that it was necessary to quantify the level of uncertainty in the estimated bedrock motion using multiple recorded events, to understand how much the simplistic model can over- or under-estimate the ground motion intensities. Inferences from the analyses suggest that the Rayleigh waves generated near the basin edge contributed significantly to the observed high frequency (f>3Hz) amplification, in addition to the amplification caused by the strong soil-rock impedance contrast at the site fundamental frequency. Models with and without considering soil non-linear response illustrate, as expected, that the linear elastic assumption severely overestimates ground motions in high frequencies for strong earthquakes, especially when the contribution of basin edge-generated Rayleigh waves becomes significant. Our analyses also demonstrate that the effect of pressure-dependent soil velocities on the high frequency ground motions is as significant as the amplification caused by the basin edge-generated Rayleigh waves.
Liquefaction-induced lateral spreading in Christchurch and surrounding suburbs during the recent Canterbury Earthquake Sequence (2010-2011) caused significant damage to structures and lifelines located in close proximity to streams and rivers. Simplified methods used in current engineering practice for predicting lateral ground displacements exhibit a high degree of epistemic uncertainty, but provide ‘order of magnitude’ estimates to appraise the hazard. We wish to compare model predictions to field measurements in order to assess the model’s capabilities and limitations with respect to Christchurch conditions. The analysis presented focuses on the widely-used empirical model of Youd et al. (2002), developed based on multi-linear regression (MLR) of case history data from lateral spreading occurrence in Japan and the US. Two issues arising from the application of this model to Christchurch were considered: • Small data set of Standard Penetration Test (SPT) and soil gradation indices (fines content FC, and mean grain size, D50) required for input. We attempt to use widely available CPT data with site specific correlations to FC and D50. • Uncertainty associated with the model input parameters and their influence on predicted displacements. This has been investigated for a specific location through a sensitivity analysis.
Prior to the devastating 2010 and 2011 earthquakes, parts of the CBD of Christchurch, New Zealand were undergoing revitalisation incorporating aspects of adaptive reuse and gentrification. Such areas were often characterised by a variety of bars, restaurants, and retail outlets of an “alternative” or “bohemian” style. These early 20th century buildings also exhibited relatively low rents and a somewhat chaotic and loosely planned property development approach by small scale developers. Almost all of these buildings were demolished following the earthquakes and a cordon placed around the CBD for several years. A paper presented at the ERES conference in 2013 presented preliminary results, from observation of post-earthquake public meetings and interviews with displaced CBD retailers. This paper highlighted a strongly held fear that the rebuild of the central city, then about to begin, would result in a very different style and cost structure from that which previously existed. As a result, permanent exclusion from the CBD of the types of businesses that previously characterised the successfully revitalised areas would occur. Five years further on, new CBD retail and office buildings have been constructed, but large areas of land between them remain vacant and the new buildings completed are often having difficulty attracting tenants. This paper reports on the further development of this long-term Christchurch case study and examines if the earlier predictions of the displaced retailers are coming true, in that a new CBD that largely mimics a suburban mall in style and tenancy mix, inherently loses some of its competitive advantage?
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).
The Canterbury Earthquake Sequence 2010-2011 (CES) induced widespread liquefaction in many parts of Christchurch city. Liquefaction was more commonly observed in the eastern suburbs and along the Avon River where the soils were characterised by thick sandy deposits with a shallow water table. On the other hand, suburbs to the north, west and south of the CBD (e.g. Riccarton, Papanui) exhibited less severe to no liquefaction. These soils were more commonly characterised by inter-layered liquefiable and non-liquefiable deposits. As part of a related large-scale study of the performance of Christchurch soils during the CES, detailed borehole data including CPT, Vs and Vp have been collected for 55 sites in Christchurch. For this subset of Christchurch sites, predictions of liquefaction triggering using the simplified method (Boulanger & Idriss, 2014) indicated that liquefaction was over-predicted for 94% of sites that did not manifest liquefaction during the CES, and under-predicted for 50% of sites that did manifest liquefaction. The focus of this study was to investigate these discrepancies between prediction and observation. To assess if these discrepancies were due to soil-layer interaction and to determine the effect that soil stratification has on the develop-ment of liquefaction and the system response of soil deposits.
The objective of the study presented herein is to assess three commonly used CPT-based liquefaction evaluation procedures and three liquefaction severity index frameworks using data from the 2010–2011 Canterbury earthquake sequence. Specifically, post-event field observations, ground motion recordings, and results from a recently completed extensive geotechnical site investigation programme at selected strong motion stations (SMSs) in the city of Christchurch and surrounding towns are used herein. Unlike similar studies that used data from free-field sites, accelerogram characteristics at the SMS locations can be used to assess the performance of liquefaction evaluation procedures prior to their use in the computation of surficial manifestation severity indices. Results from this study indicate that for cases with evidence of liquefaction triggering in the accelerograms, the majority of liquefaction evaluation procedures yielded correct predictions, regardless of whether surficial manifestation of liquefaction was evident or not. For cases with no evidence of liquefaction in the accelerograms (and no observed surficial evidence of liquefaction triggering), the majority of liquefaction evaluation procedures predicted liquefaction was triggered. When all cases are used to assess the performance of liquefaction severity index frameworks, a poor correlation is shown between the observed severity of liquefaction surface manifestation and the calculated severity indices. However, only using those cases where the liquefaction evaluation procedures yielded correct predictions, there is an improvement in the correlation, with the Liquefaction Severity Number (LSN) being the best performing of the frameworks investigated herein. However scatter in the relationship between the observed and calculated surficial manifestation still remains for all liquefaction severity index frameworks.
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.
Data from the 2010-2011 Canterbury earthquake sequence (CES) provides an unprecedented opportunity to assess and advance the current state of practice for evaluating liquefaction triggering. Towards this end, select case histories from the CES are used herein to assess the predictive capabilities of three alternative CPT-based simplified liquefaction evaluation procedures: Robertson and Wride (1998); Moss et al. (2006); and Idriss and Boulanger (2008). Additionally, the Liquefaction Potential Index (LPI) framework for predicting the severity of surficial liquefaction manifestations is also used to assess the predictive capabilities of the liquefaction evaluation procedures. Although it is not without limitations, use of the LPI framework for this purpose circumvents the need for selecting “critical” layers and their representative properties for study sites, which inherently involves subjectivity and thus has been a point of contention among researchers. It was found that while all the assessed liquefaction triggering evaluation procedures performed well for the parameter ranges of the sites analyzed, the procedure proposed by Idriss and Boulanger (2008) yielded predictions that are more consistent with field observations than the other procedures. However, use of the Idriss and Boulanger (2008) procedure in conjunction with a Christchurch-specific correlation to estimate fines content showed a decreased performance relative to using a generic fines content correlation. As a result, the fines correction for the Idriss and Boulanger (2008) procedure needs further study.
Predictive modelling provides an efficient means to analyse the coastal environment and generate knowledge for long term urban planning. In this study, the numerical models SWAN and XBeach were incorporated into the ESRI ArcGIS interface by means of the BeachMMtool. This was applied to the Greater Christchurch coastal environment to simulate geomorphological evolution through hydrodynamic forcing. Simulations were performed using the recent sea level rise predictions by the Intergovernmental Panel on Climate Change (2013) to determine whether the statutory requirements outlined in the New Zealand Coastal Policy Statement 2010 are consistent with central, regional and district designations. Our results indicate that current land use zoning in Greater Christchurch is not consistent with these predictions. This is because coastal hazard risk has not been thoroughly quantified during the process of installing the Canterbury Earthquake Recovery Authority residential red zone. However, the Christchurch City Council’s flood management area does provide an extent to which managed coastal retreat is a real option. The results of this research suggest that progradation will continue to occur along the Christchurch foreshore due to the net sediment flux retaining an onshore direction and the current hydrodynamic activity not being strong enough to move sediment offshore. However, inundation during periods of storm surge poses a risk to human habitation on low lying areas around the Avon-Heathcote Estuary and the Brooklands lagoon.
High-quality ground motion records are required for engineering applications including response history analysis, seismic hazard development, and validation of physics-based ground motion simulations. However, the determination of whether a ground motion record is high-quality is poorly handled by automation with mathematical functions and can become prohibitive if done manually. Machine learning applications are well-suited to this problem, and a previous feed-forward neural network was developed (Bellagamba et al. 2019) to determine high-quality records from small crustal events in the Canterbury and Wellington regions for simulation validation. This prior work was however limited by the omission of moderate-to-large magnitude events and those from other tectonic environments, as well as a lack of explicit determination of the minimum usable frequency of the ground motion. To address these shortcomings, an updated neural network was developed to predict the quality of ground motion records for all magnitudes and all tectonic sources—active shallow crustal, subduction intraslab, and subduction interface—in New Zealand. The predictive performance of the previous feed-forward neural network was matched by the neural network in the domain of small crustal records, and this level of predictive performance is now extended to all source magnitudes and types in New Zealand making the neural network applicable to global ground motion databases. Furthermore, the neural network provides quality and minimum usable frequency predictions for each of the three orthogonal components of a record which may then be mapped into a binary quality decision or otherwise applied as desired. This framework provides flexibility for the end user to predict high-quality records with various acceptability thresholds allowing for this neural network to be used in a range of applications.
The overarching goal of this dissertation is to improve predictive capabilities of geotechnical seismic site response analyses by incorporating additional salient physical phenomena that influence site effects. Specifically, multidimensional wave-propagation effects that are neglected in conventional 1D site response analyses are incorporated by: (1) combining results of 3D regional-scale simulations with 1D nonlinear wave-propagation site response analysis, and (2) modelling soil heterogeneity in 2D site response analyses using spatially-correlated random fields to perturb soil properties. A method to combine results from 3D hybrid physics-based ground motion simulations with site-specific nonlinear site response analyses was developed. The 3D simulations capture 3D ground motion phenomena on a regional scale, while the 1D nonlinear site response, which is informed by detailed site-specific soil characterization data, can capture site effects more rigorously. Simulations of 11 moderate-to-large earthquakes from the 2010-2011 Canterbury Earthquake Sequence (CES) at 20 strong motion stations (SMS) were used to validate simulations with observed ground motions. The predictions were compared to those from an empirically-based ground motion model (GMM), and from 3D simulations with simplified VS30- based site effects modelling. By comparing all predictions to observations at seismic recording stations, it was found that the 3D physics-based simulations can predict ground motions with comparable bias and uncertainty as the GMM, albeit, with significantly lower bias at long periods. Additionally, the explicit modelling of nonlinear site-response improves predictions significantly compared to the simplified VS30-based approach for soft-soil or atypical sites that exhibit exceptionally strong site effects. A method to account for the spatial variability of soils and wave scattering in 2D site response analyses was developed and validated against a database of vertical array sites in California. The inputs required to run the 2D analyses are nominally the same as those required for 1D analyses (except for spatial correlation parameters), enabling easier adoption in practice. The first step was to create the platform and workflow, and to perform a sensitivity study involving 5,400 2D model realizations to investigate the influence of random field input parameters on wave scattering and site response. Boundary conditions were carefully assessed to understand their effect on the modelled response and select appropriate assumptions for use on a 2D model with lateral heterogeneities. Multiple ground-motion intensity measures (IMs) were analyzed to quantify the influence from random field input parameters and boundary conditions. It was found that this method is capable of scattering seismic waves and creating spatially-varying ground motions at the ground surface. The redistribution of ground-motion energy across wider frequency bands, and the scattering attenuation of high-frequency waves in 2D analyses, resemble features observed in empirical transfer functions (ETFs) computed in other studies. The developed 2D method was subsequently extended to more complicated multi-layer soil profiles and applied to a database of 21 vertical array sites in California to test its appropriate- ness for future predictions. Again, different boundary condition and input motion assumptions were explored to extend the method to the in-situ conditions of a vertical array (with a sensor embedded in the soil). ETFs were compared to theoretical transfer functions (TTFs) from conventional 1D analyses and 2D analyses with heterogeneity. Residuals of transfer-function- based IMs, and IMs of surface ground motions, were also used as validation metrics. The spatial variability of transfer-function-based IMs was estimated from 2D models and compared to the event-to-event variability from ETFs. This method was found capable of significantly improving predictions of median ETF amplification factors, especially for sites that display higher event-to-event variability. For sites that are well represented by 1D methods, the 2D approach can underpredict amplification factors at higher modes, suggesting that the level of heterogeneity may be over-represented by the 2D random field models used in this study.
This study uses 44 high quality liquefaction case histories taken from 22 locations affected by the 2010-2011 Canterbury earthquake sequence to evaluate four commonly used CPT-VS correlations (i.e., Robertson, 2009; Hegazy and Mayne, 2006; Andrus et al., 2007; McGann et al., 2015b). Co-located CPT soundings and VS profiles, developed from surface wave testing, were obtained at 22 locations and case histories were developed for the Mw 7.1, 4 September 2010 Darfield and Mw 6.2, 22 February 2011 Christchurch earthquakes. The CPT soundings are used to generate VS profiles using each of four CPT-VS correlations. These correlated VS profiles are used to estimate the factor of safety against liquefaction using the Kayen et al. (2013) VS-based simplified liquefaction evaluation procedure. An error index is used to quantify the predictive capabilities of these correlations in relation to the observations of liquefaction (or the lack thereof). Additionally, the error indices from the CPT-correlated VS profiles are compared to those obtained using: (1) the Kayen et al. (2013) procedure with surface wave-derived VS profiles, and (2) the Idriss and Boulanger (2008) CPT-based liquefaction evaluation procedure. Based on the error indices, the evaluation procedures based on direct measurements of either CPT or VS provided more accurate liquefaction triggering estimates than those obtained from any of the CPT-VS correlations. However, the performance of the CPT-VS correlations varied, with the Robertson (2009) and Hegazy and Mayne (2006) correlations performing relatively poorly for the Christchurch soils and the Andrus et al. (2007) and McGann et al. (2015b) correlations performing better. The McGann et al. (2015b) correlation had the lowest error indices of the CPT-VS correlations tested, however, none of the CPT-VS correlations provided accurate enough VS predictions to be used for the evaluation of liquefaction triggering using the VS-based liquefaction evaluation procedures.
Predictive modelling provides an efficient means to analyse the coastal environment and generate knowledge for long term urban planning. In this study, the numerical models SWAN and XBeach were incorporated into the ESRI ArcGIS interface by means of the BeachMMtool. This was applied to the Greater Christchurch coastal environment to simulate geomorphological evolution through hydrodynamic forcing. Simulations were performed using the recent sea level rise predictions by the Intergovernmental Panel on Climate Change (2013) to determine whether the statutory requirements outlined in the New Zealand Coastal Policy Statement 2010 are consistent with central, regional and district designations. Our results indicate that current land use zoning in Greater Christchurch is not consistent with these predictions. This is because coastal hazard risk has not been thoroughly quantified during the process of installing the Canterbury Earthquake Recovery Authority residential red zone. However, the Christchurch City Council’s flood management area does provide an extent to which managed coastal retreat is a real option. The results of this research suggest that progradation will continue to occur along the Christchurch foreshore due to the net sediment flux retaining an onshore direction and the current hydrodynamic activity not being strong enough to move sediment offshore. However, inundation during periods of storm surge poses a risk to human habitation on low lying areas around the Avon-Heathcote Estuary and the Brooklands lagoon.
The 2015 New Zealand strong-motion database provides a wealth of new strong motion data for engineering applications. An important component of this database is the compilation of new site metadata, describing the soil conditions and site response at GeoNet strong motion stations. We have assessed and compiled four key site parameters for the ~460 GeoNet stations that recorded significant historical ground motions. Parameters include: site classification (NZS1170.5), Vs30, fundamental site period (Tsite) and depth to bedrock (Z1.0, i.e. depth to material with Vs > 1000 m/s). In addition, we have assigned a quality estimate (Quality 1 – 3) to these parameters to provide a qualitative estimate of the uncertainty. New highquality Tsite estimates have largely been obtained from newly available HVSR amplification curves and spectral ratios from inversion of regional strong motion data that has been reconciled with available geological information. Good quality Vs30 estimates, typically in urban centres, have also been incorporated following recent studies. Where site-specific measurements of Vs30 are not available, Vs30 is estimated based on surface geology following national Vs30 maps. New Z1.0 values have been provided from 3D subsurface models for Canterbury and Wellington. This database will be used in efforts to guide development and testing of new and existing ground motion prediction models in New Zealand. In particular, it will allow reexamination of the most important site parameters that control and predict site response in a New Zealand setting. Furthermore, it can be used to provide information about suitable rock reference sites for seismological research, and as a guide to site-specific references in the literature. We discuss compilation of the database, preliminary insights so far, and future directions.
The full scale, in-situ investigations of instrumented buildings present an excellent opportunity to observe their dynamic response in as-built environment, which includes all the real physical properties of a structure under study and its surroundings. The recorded responses can be used for better understanding of behavior of structures by extracting their dynamic characteristics. It is significantly valuable to examine the behavior of buildings under different excitation scenarios. The trends in dynamic characteristics, such as modal frequencies and damping ratios, thus developed can provide quantitative data for the variations in the behavior of buildings. Moreover, such studies provide invaluable information for the development and calibration of realistic models for the prediction of seismic response of structures in model updating and structural health monitoring studies. This thesis comprises two parts. The first part presents an evaluation of seismic responses of two instrumented three storey RC buildings under a selection of 50 earthquakes and behavioral changes after Ms=7.1 Darfield (2010) and Ms=6.3 Christchurch (2011) earthquakes for an instrumented eight story RC building. The dynamic characteristics of the instrumented buildings were identified using state-of-the-art N4SID system identification technique. Seismic response trends were developed for the three storey instrumented buildings in light of the identified frequencies and the peak response accelerations (PRA). Frequencies were observed to decrease with excitation level while no trends are discernible for the damping ratios. Soil-structure interaction (SSI) effects were also determined to ascertain their contribution in the seismic response. For the eight storey building, it was found through system identification that strong nonlinearities in the structural response occurred and manifested themselves in all identified natural frequencies of the building that exhibited a marked decrease during the strong motion duration compared to the pre-Darfield earthquakes. Evidence of foundation rocking was also found that led to a slight decrease in the identified modal frequencies. Permanent stiffness loss was also observed after the strong motion events. The second part constitutes developing and calibrating finite element model (FEM) of the instrumented three storey RC building with a shear core. A three dimensional FEM of the building is developed in stages to analyze the effect of structural, non-structural components (NSCs) and SSI on the building dynamics. Further to accurately replicate the response of the building following the response trends developed in the first part of the thesis, sensitivity based model updating technique was applied. The FEMs were calibrated by tuning the updating parameters which are stiffnesses of concrete, NSCs and soil. The updating parameters were found to generally follow decreasing trends with the excitation level. Finally, the updated FEM was used in time history analyses to assess the building seismic performance at the serviceability limit state shaking. Overall, this research will contribute towards better understanding and prediction of the behavior of structures subjected to ground motion.
This thesis is concerned with the effects of lateral confining reinforcement on the ductile behaviour of reinforced concrete columns. The contents of the chapters are summarized as follows. In Chapter one, the general problems in seismic design are discussed and earthquake design methods based on the ductile design approach are described. Japanese, New Zealand and United States design codes are compared. Finally, the scope of this research project is outlined. In Chapter two, after reviewing previous research on confined concrete, the factors which affect the effectiveness of lateral confinement are discussed. Especially the effects of the yield strength of transverse reinforcement, the compressive strength of plain concrete and the strain gradient in the column section due to bending are discussed based on tests which were conducted by the author et al at Kyoto University and Akashi Technological College, Japan. In the axial compression tests on spirally reinforced concrete cylinders (150 mm in diameter by 300 mm in height), the yield strength of transverse reinforcement and the compressive strength of plain concrete were varied from 161 MPa to 1352 MPa and from 17 MPa to 60 MPa, respectively, as experimental parameters. It is found that, when high strength spirals are used as confining reinforcement, the strength and ductility of the confined core concrete are remarkably enhanced but need to be estimated assuming several failure modes which could occur. These are based on the observations that concrete cylinders with high strength spirals suddenly failed at a concrete compressive strain of 2 to 3.5 % due to explosive crushing of the core concrete between the spiral bars or due to bearing failure of the core concrete immediately beneath the spiral bars, while the concrete cylinders with ordinary strength spirals failed in a gentle manner normally observed. In addition, eccentric loading tests were conducted on concrete columns with 200 mm square section confined by square spirals. It is found that the effectiveness of confining reinforcement is reduced by the presence of the strain gradient along the transverse section of column. In Chapter three, the effectiveness of transverse reinforcement with various types of anchorage details which simplify the fabrication of reinforcing cages are investigated. Eight reinforced concrete columns, with either 400 mm or 550 mm square cross sections, were tested subjected to axial compression loading and cyclic lateral loading which simulated a severe earthquake. The transverse reinforcement consisted of arrangements of square perimeter hoops with 135° end hooks, cross ties with 90° and 135° or 180° end hooks, and 'U' and 'J' shaped cross ties and perimeter hoops with tension splices. Conclusions are reached with regard to the effectiveness of the tested anchorage details in the plastic hinge regions of columns designed for earthquake resistance. In Chapter four, the effectiveness of interlocking spirals as transverse reinforcement is studied. Firstly, the general aspects and the related problems of interlocking spirals to provide adequate ductility in the potential plastic hinge region of columns are discussed, referring to the provisions in the New Zealand code,the CALTRANS (California Transportation Authority) code and other related codes. Secondly, based on those discussions, a design method to securely interlock the spirals is proposed. Thirdly, the effectiveness of interlocking spirals is assessed based on column tests conducted as part of this study. Three columns with interlocking spirals and, for comparison, one rectangular column with rectangular hoopsandcross ties, were tested under cyclic horizontal loading which simulated a severe earthquake. The sections of those columns were 400 mm by 600 mm. In Chapter five, analytical models to investigate the buckling behaviour of longitudinal reinforcement restrained by cross ties with 90° and 135° end hooks and by peripheral hoops are proposed. The analyzed results using the proposed models compare well with the experimental observations described in Chapter three. Using those proposed models, a method to check the effectiveness of cross ties with 90° and 135° end hooks is proposed for practical design purposes. In Chapter six, a theory for the prediction of the ultimate longitudinal compressive concrete strain at the stage of first hoop fracture referred to as the "Energy Balance Theory", which has been developed by Mander, Priestley and Park at University of Canterbury, is introduced. After discussing the problems in the "Energy Balance Theory", a modified theory for the prediction of the ultimate longitudinal compressive concrete strain at the stage of first hoop fracture is proposed. The predictions from the modified theory are found to compare well with previous experimental results.
