This poster discusses several possible approaches by which the nonlinear response of surficial soils can be explicitly modelled in physics-based ground motion simulations, focusing on the relative advantages and limitations of the various methodologies. These methods include fully-coupled 3D simulation models that directly allow soil nonlinearity in surficial soils, the domain reduction method for decomposing the physical domain into multiple subdomains for separate simulation, conventional site response analysis uncoupled from the simulations, and finally, the use of simple empirically based site amplification factors We provide the methodology for an ongoing study to explicitly incorporate soil nonlinearity into hybrid broadband simulations of the 2010-2011 Canterbury, New Zealand earthquakes.
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).
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.
In this paper, we perform hybrid broadband (0-10 Hz) ground motion simulations for the ten most significant events (Mw 4.7-7.1) in the 2010-2011 Canterbury earthquake sequence. Taking advantage of having repeated recordings at same stations, we validate our simulations using both recordings and an empirically-developed ground motion prediction equation (GMPE). The simulation clearly captures the sedimentary basin amplification and the rupture directivity effects. Quantitative comparisons of the simulations with both recordings and the GMPE, as well as analyses of the total residuals (indicating model bias) show that simulations perform better than the empirical GMPE, especially for long period. To scrutinize the ground motion variability, we partitioned the total residuals into different components. The total residual appears to be unbiased, and the use of a 3D velocity structure reduces the long period systematic bias particularly for stations located close to the Banks Peninsula volcanic area.
Tourism is New Zealand’s fourth largest industry, providing jobs for thousands of New
Zealanders and significant foreign capital for the nation’s economy. Of concern to ministry and industry
decision makers is the “spatial yield” of these tourists which takes into account the spatial and temporal
contributions of their movements in terms of economic, cultural and environmental impacts. We have
developed an agent-based model of tourism movements to simulate these impacts and to allow for the
evaluation of different scenarios (such as increases in petrol prices or variations in currency exchange rates)
on the behaviours of those tourists. In order to develop realistic and grounded heuristics for the model,
interview protocols were developed in order to identify the key drivers in tourists’ decision making process.
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.
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).
This dissertation addresses a diverse range of topics in the area of physics-based ground motion simulation with particular focus on the Canterbury, New Zealand region. The objectives achieved provide the means to perform hybrid broadband ground motion simulation and subsequently validates the simulation methodology employed. In particu- lar, the following topics are addressed: the development of a 3D seismic velocity model of the Canterbury region for broadband ground motion simulation; the development of a 3D geologic model of the interbedded Quaternary formations to provide insight on observed ground motions; and the investigation of systematic effects through ground motion sim- ulation of small-to-moderate magnitude earthquakes. The paragraphs below outline each contribution in more detail. As a means to perform hybrid broadband ground motion simulation, a 3D model of the geologic structure and associated seismic velocities in the Canterbury region is devel- oped utilising data from depth-converted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model explicitly characterises five significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age, including the Banks Peninsula volcanics, which are noted to strongly influence seismic wave propagation. The Basement surface represents the base of the Canterbury sedimentary basin, where a large impedance contrast exists re- sulting in basin-generated waves. Seismic velocities for the lithological units between the geologic surfaces are derived from well logs, seismic reflection surveys, root mean square stacking velocities, empirical correlations, and benchmarked against a regional crustal model, thus providing the necessary information for a Canterbury velocity model for use in broadband seismic wave propagation. A 3D high-resolution model of the Quaternary geologic stratigraphic sequence in the Canterbury region is also developed utilising datasets of 527 high-quality water well logs, and 377 near-surface cone penetration test records. The model, developed using geostatistical Kriging, represents the complex interbedded regional Quaternary geology by characterising the boundaries between significant interbedded geologic formations as 3D surfaces including explicit modelling of the formation unconformities resulting from the Banks Peninsula volcanics. The stratigraphic layering present can result in complex wave propagation. The most prevalent trend observed in the surfaces was the downward dip from inland to the eastern coastline as a result of the dominant fluvial depositional environment of the terrestrial gravel formations. The developed model provides a benefi- cial contribution towards developing a comprehensive understanding of recorded ground motions in the region and also providing the necessary information for future site char- acterisation and site response analyses. To highlight the practicality of the model, an example illustrating the role of the model in constraining surface wave analysis-based shear wave velocity profiling is illustrated along with the calculation of transfer functions to quantify the effect of the interbedded geology on wave propagation. Lastly, an investigation of systematic biases in the (Graves and Pitarka, 2010, 2015) ground motion simulation methodology and the specific inputs used for the Canterbury region is presented considering 144 small-to-moderate magnitude earthquakes. In the simulation of these earthquakes, the 3D Canterbury Velocity Model, developed as a part of this dissertation, is used for the low-frequency simulation, and a regional 1D velocity model for the high-frequency simulation. Representative results for individual earthquake sources are first presented to highlight the characteristics of the small-to-moderate mag- nitude earthquake simulations through waveforms, intensity measure scaling with source- to-site distance, and spectral bias of the individual events. Subsequently, a residual de- composition is performed to examine the between- and within-event residuals between observed data, and simulated and empirical predictions. By decomposing the residuals into between- and within-event residuals, the biases in source, path and site effects, and their causes, can be inferred. The residuals are comprehensively examined considering their aggregated characteristics, dependence on predictor variables, spatial distribution, and site-specific effects. The results of the simulation are also benchmarked against empir- ical ground motion models, where their similarities manifest from common components in their prediction. Ultimately, suggestions to improve the predictive capability of the simulations are presented as a result of the analysis.
© 2017 The Royal Society of New Zealand. This paper discusses simulated ground motion intensity, and its underlying modelling assumptions, for great earthquakes on the Alpine Fault. The simulations utilise the latest understanding of wave propagation physics, kinematic earthquake rupture descriptions and the three-dimensional nature of the Earth's crust in the South Island of New Zealand. The effect of hypocentre location is explicitly examined, which is found to lead to significant differences in ground motion intensities (quantified in the form of peak ground velocity, PGV) over the northern half and southwest of the South Island. Comparison with previously adopted empirical ground motion models also illustrates that the simulations, which explicitly model rupture directivity and basin-generated surface waves, lead to notably larger PGV amplitudes than the empirical predictions in the northern half of the South Island and Canterbury. The simulations performed in this paper have been adopted, as one possible ground motion prediction, in the ‘Project AF8’ Civil Defence Emergency Management exercise scenario. The similarity of the modelled ground motion features with those observed in recent worldwide earthquakes as well as similar simulations in other regions, and the notably higher simulated amplitudes than those from empirical predictions, may warrant a re-examination of regional impact assessments for major Alpine Fault earthquakes.
This study examines the performance of nonlinear total-stress wave-propagation site response analysis 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 3-dimensional ground motion phenomena at the regional scale, as well as detailed site effects and soil nonlinearity 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).
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.
his poster presents the ongoing development of a 3D Canterbury seismic velocity model which will be used in physics-based hybrid broadband ground motion simulation of the 2010-2011 Canterbury earthquakes. Velocity models must sufficiently represent critical aspects of the crustal structure over multiple length scales which will influence the results of the simulations. As a result, numerous sources of data are utilized in order to provide adequate resolution where necessary. Figure 2: (a) Seismic reflection line showing P-wave velocities and significant geologic horizons (Barnes et al. 2011), and (b) Shear wave profiles at 10 locations (Stokoe et al. 2013). Figure 4: Cross sections of the current version of the Canterbury velocity model to depths of 10km as shown in Figure 1: (a) at a constant latitude value of -43.6˚, and (b) at a constant longitude value of 172.64˚. 3. Ground Surface and Geologic Horizon Models Figure 3: (a) Ground surface model derived from numerous available digital elevation models, and (b) Base of the Quaternary sediments derived from structural contours and seismic reflection line elevations. The Canterbury region has a unique and complex geology which likely has a significant impact on strong ground motions, in particular the deep and loose deposits of the Canterbury basin. The Canterbury basin has several implications on seismic wave phenomena such as long period ground motion amplification and wave guide effects. Using a realistic 3D seismic velocity model in physics-based ground motion simulation will implicitly account for such effects and the resultant simulated ground motions can be studied to gain a fundamental understanding of the salient ground motion phenomena which occurred during the Canterbury earthquakes, and the potential for repeat occurrences in the Canterbury region. Figure 1 shows the current model domain as a rectangular area between Lat=[-43.2˚,-44.0˚], and Lon=[171.5˚,173.0˚]. This essentially spans the area between the foot of the Southern Alps in the North West to Banks Peninsula in the East. Currently the model extends to a depth of 50km below sea level.
None
The research is funded by Callaghan Innovation (grant number MAIN1901/PROP-69059-FELLOW-MAIN) and the Ministry of Transport New Zealand in partnership with Mainfreight Limited. Need – The freight industry is facing challenges related to climate change, including natural hazards and carbon emissions. These challenges impact the efficiency of freight networks, increase costs, and negatively affect delivery times. To address these challenges, freight logistics modelling should consider multiple variables, such as natural hazards, sustainability, and emission reduction strategies. Freight operations are complex, involving various factors that contribute to randomness, such as the volume of freight being transported, the location of customers, and truck routes. Conventional methods have limitations in simulating a large number of variables. Hence, there is a need to develop a method that can incorporate multiple variables and support freight sustainable development. Method - A minimal viable model (MVM) method was proposed to elicit tacit information from industrial clients for building a minimally sufficient simulation model at the early modelling stages. The discrete-event simulation (DES) method was applied using Arena® software to create simulation models for the Auckland and Christchurch corridor, including regional pick-up and delivery (PUD) models, Christchurch city delivery models, and linehaul models. Stochastic variables in freight operations such as consignment attributes, customer locations, and truck routes were incorporated in the simulation. The geographic information system (GIS) software ArcGIS Pro® was used to identify and analyse industrial data. The results obtained from the GIS software were applied to create DES models. Life cycle assessment (LCA) models were developed for both diesel and battery electric (BE) trucks to compare their life cycle greenhouse gas (GHG) emissions and total cost of ownership (TCO) and support GHG emissions reduction. The line-haul model also included natural hazards in several scenarios, and the simulation was used to forecast the stock level of Auckland and Christchurch depots in response to each corresponding scenario. Results – DES is a powerful technique that can be employed to simulate and evaluate freight operations that exhibit high levels of variability, such as regional pickup and delivery (PUD) and linehaul. Through DES, it becomes possible to analyse multiple factors within freight operations, including transportation modes, routes, scheduling, and processing times, thereby offering valuable insights into the performance, efficiency, and reliability of the system. In addition, GIS is a useful tool for analysing and visualizing spatial data in freight operations. This is exemplified by their ability to simulate the travelling salesman problem (TSP) and conduct cluster analysis. Consequently, the integration of GIS into DES modelling is essential for improving the accuracy and reliability of freight operations analysis. The outcomes of the simulation were utilised to evaluate the ecological impact of freight transport by performing emission calculations and generating low-carbon scenarios to identify approaches for reducing the carbon footprint. LCA models were developed based on simulation results. Results showed that battery-electric trucks (BE) produced more greenhouse gas (GHG) emissions in the cradle phase due to battery manufacturing but substantially less GHG emissions in the use phase because of New Zealand's mostly renewable energy sources. While the transition to BE could significantly reduce emissions, the financial aspect is not compelling, as the total cost of ownership (TCO) for the BE truck was about the same for ten years, despite a higher capital investment for the BE. Moreover, external incentives are necessary to justify a shift to BE trucks. By using simulation methods, the effectiveness of response plans for natural hazards can be evaluated, and the system's vulnerabilities can be identified and mitigated to minimize the risk of disruption. Simulation models can also be utilized to simulate adaptation plans to enhance the system's resilience to natural disasters. Novel contributions – The study employed a combination of DES and GIS methods to incorporate a large number of stochastic variables and driver’s decisions into freight logistics modelling. Various realistic operational scenarios were simulated, including customer clustering and PUD truck allocation. This showed that complex pickup and delivery routes with high daily variability can be represented using a model of roads and intersections. Geographic regions of high customer density, along with high daily variability could be represented by a two-tier architecture. The method could also identify delivery runs for a whole city, which has potential usefulness in market expansion to new territories. In addition, a model was developed to address carbon emissions and total cost of ownership of battery electric trucks. This showed that the transition was not straightforward because the economics were not compelling, and that policy interventions – a variety were suggested - could be necessary to encourage the transition to decarbonised freight transport. A model was developed to represent the effect of natural disasters – such as earthquake and climate change – on road travel and detour times in the line haul freight context for New Zealand. From this it was possible to predict the effects on stock levels for a variety of disruption scenarios (ferry interruption, road detours). Results indicated that some centres rather than others may face higher pressure and longer-term disturbance after the disaster subsided. Remedies including coastal shipping were modelled and shown to have the potential to limit the adverse effects. A philosophical contribution was the development of a methodology to adapt the agile method into the modelling process. This has the potential to improve the clarification of client objectives and the validity of the resulting model.
This paper presents the ongoing development of a new 3D seismic velocity model of Canterbury, New Zealand. The model explicitly represents the Canterbury sedimentary basin, and other significant geologic horizons, which are expected to have important implications on observed ground motions. The model utilizes numerous sources of data, including 3D regional tomography with a variable-depth inferred Moho, seismic reflection survey lines, geotechnical boreholes and well logs, spectral analysis of surface waves, and CPT logs which provide velocity constraints over their respective ranges of application. The model provides P- and S-wave velocity and density (i.e. Vp, Vs and p) over a grid of input points, and is presently being utilized in broadband ground motion simulations of the 2010-2011 Canterbury earthquakes. Comparison of simulated ground motions with those observed in the 2010-2011 Canterbury earthquakes will help provide a better understanding of the salient physical processes which characterized the unique set of strong ground motions recorded in this sequence of earthquake events.
This poster aims to present fragility functions for pipelines buried in liquefaction-prone soils. Existing fragility models used to quantify losses can be based on old data or use complex metrics. Addressing these issues, the proposed functions are based on the Christchurch network and soil and utilizes the Canterbury earthquake sequence (CES) data, partially represented in Figure 1. Figure 1 (a) presents the pipe failure dataset, which describes the date, location and pipe on which failures occurred. Figure 1 (b) shows the simulated ground motion intensity median of the 22nd February 2011 earthquake. To develop the model, the network and soil characteristics have also been utilized
We present initial results from a set of three-dimensional (3D) deterministic earthquake ground motion simulations for the northern Canterbury plains, Christchurch and the Banks Peninsula region, which explicitly incorporate the effects of the surface topography. The simu-lations are done using Hercules, an octree-based finite-element parallel software for solving 3D seismic wave propagation problems in heterogeneous media under kinematic faulting. We describe the efforts undertaken to couple Hercules with the South Island Velocity Model (SIVM), which included changes to the SIVM code in order to allow for single repetitive que-ries and thus achieve a seamless finite-element meshing process within the end-to-end ap-proach adopted in Hercules. We present our selection of the region of interest, which corre-sponds to an area of about 120 km × 120 km, with the 3D model reaching a depth of 60 km. Initial simulation parameters are set for relatively high minimum shear wave velocity and a low maximum frequency, which we are progressively scaling up as computing resources permit. While the effects of topography are typically more important at higher frequencies and low seismic velocities, even at this initial stage of our efforts (with a maximum of 2 Hz and a mini-mum of 500 m/s), it is possible to observe the importance of the topography in the response of some key locations within our model. To highlight these effects we compare the results of the 3D topographic model with respect to those of a flat (squashed) 3D model. We draw rele-vant conclusions from the study of topographic effects during earthquakes for this region and describe our plans for future work.
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.
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.
Damage distribution maps from strong earthquakes and recorded data from field experiments have repeatedly shown that the ground surface topography and subsurface stratigraphy play a decisive role in shaping the ground motion characteristics at a site. Published theoretical studies qualitatively agree with observations from past seismic events and experiments; quantitatively, however, they systematically underestimate the absolute level of topographic amplification up to an order of magnitude or more in some cases. We have hypothesized in previous work that this discrepancy stems from idealizations of the geometry, material properties, and incident motion characteristics that most theoretical studies make. In this study, we perform numerical simulations of seismic wave propagation in heterogeneous media with arbitrary ground surface geometry, and compare results with high quality field recordings from a site with strong surface topography. Our goal is to explore whether high-fidelity simulations and realistic numerical models can – contrary to theoretical models – capture quantitatively the frequency and amplitude characteristics of topographic effects. For validation, we use field data from a linear array of nine portable seismometers that we deployed on Mount Pleasant and Heathcote Valley, Christchurch, New Zealand, and we compute empirical standard spectral ratios (SSR) and single-station horizontal-to-vertical spectral ratios (HVSR). The instruments recorded ambient vibrations and remote earthquakes for a period of two months (March-April 2017). We next perform two-dimensional wave propagation simulations using the explicit finite difference code FLAC. We construct our numerical model using a high-resolution (8m) Digital Elevation Map (DEM) available for the site, an estimated subsurface stratigraphy consistent with the geomorphology of the site, and soil properties estimated from in-situ and non-destructive tests. We subject the model to in-plane and out-of-plane incident motions that span a broadband frequency range (0.1-20Hz). Numerical and empirical spectral ratios from our blind prediction are found in very good quantitative agreement for stations on the slope of Mount Pleasant and on the surface of Heathcote Valley, across a wide range of frequencies that reveal the role of topography, soil amplification and basin edge focusing on the distribution of ground surface motion.
The 2010-2011 Canterbury earthquakes were recorded over a dense strong motion network in the near-source region, yielding significant observational evidence of seismic complexities, and a basis for interpretation of multi-disciplinary datasets and induced damage to the natural and built environment. This paper provides an overview of observed strong motions from these events and retrospective comparisons with both empirical and physics-based ground motion models. Both empirical and physics-based methods provide good predictions of observations at short vibration periods in an average sense. However, observed ground motion amplitudes at specific locations, such as Heathcote Valley, are seen to systematically depart from ‘average’ empirical predictions as a result of near surface stratigraphic and topographic features which are well modelled via sitespecific response analyses. Significant insight into the long period bias in empirical predictions is obtained from the use of hybrid broadband ground motion simulation. The comparison of both empirical and physics-based simulations against a set of 10 events in the sequence clearly illustrates the potential for simulations to improve ground motion and site response prediction, both at present, and further in the future.
SeisFinder is an open-source web service developed by QuakeCoRE and the University of Canterbury, focused on enabling the extraction of output data from computationally intensive earthquake resilience calculations. Currently, SeisFinder allows users to select historical or future events and retrieve ground motion simulation outputs for requested geographical locations. This data can be used as input for other resilience calculations, such as dynamic response history analysis. SeisFinder was developed using Django, a high-level python web framework, and uses a postgreSQL database. Because our large-scale computationally-intensive numerical ground motion simulations produce big data, the actual data is stored in file systems, while the metadata is stored in the database.
We present ground motion simulations of the Porters Pass (PP) fault in the Canterbury region of New Zealand; a major active source near Christchurch city. The active segment of the PP fault has an inferred length of 82 km and a mostly strike-slip sense of movement. The PP fault slip makes up approximately 10% of the total 37 mm/yr margin-parallel plate motion and also comprises a significant proportion of the total strain budget in regional tectonics. Given that the closest segment of the fault is less than 45 km from Christchurch city, the PP fault is crucial for accurate earthquake hazard assessment for this major population centre. We have employed the hybrid simulation methodology of Graves and Pitarka (2010, 2015), which combines low (f<1 Hz) and high (f>1 Hz) frequencies into a broadband spectrum. We have used validations from three moderate magnitude events (𝑀𝑤4.6 Sept 04, 2010; 𝑀𝑤4.6 Nov 06, 2010; 𝑀𝑤4.9 Apr 29, 2011) to build confidence for the 𝑀𝑤 > 7 PP simulations. Thus far, our simulations include multiple rupture scenarios which test the impacts of hypocentre location and the finite-fault stochastic rupture representation of the source itself. In particular, we have identified the need to use location-specific 1D 𝑉𝑠/𝑉𝑝 models for the high frequency part of the simulations to better match observations.
The Canterbury Earthquake Sequence (CES) of 2010-2011 produced large seismic moments up to Mw 7.1. These large, near-to-surface (<15 km) ruptures triggered >6,000 rockfall boulders on the Port Hills of Christchurch, many of which impacted houses and affected the livelihoods of people within the impacted area. From these disastrous and unpredicted natural events a need arose to be able to assess the areas affected by rockfall events in the future, where it is known that a rockfall is possible from a specific source outcrop but the potential boulder runout and dynamics are not understood. The distribution of rockfall deposits is largely constrained by the physical properties and processes of the boulder and its motion such as block density, shape and size, block velocity, bounce height, impact and rebound angle, as well as the properties of the substrate. Numerical rockfall models go some way to accounting for all the complex factors in an algorithm, commonly parameterised in a user interface where site-specific effects can be calibrated. Calibration of these algorithms requires thorough field checks and often experimental practises. The purpose of this project, which began immediately following the most destructive rupture of the CES (February 22, 2011), is to collate data to characterise boulder falls, and to use this information, supplemented by a set of anthropogenic boulder fall data, to perform an in-depth calibration of the three-dimensional numerical rockfall model RAMMS::Rockfall. The thesis covers the following topics: • Use of field data to calibrate RAMMS. Boulder impact trails in the loess-colluvium soils at Rapaki Bay have been used to estimate ranges of boulder velocities and bounce heights. RAMMS results replicate field data closely; it is concluded that the model is appropriate for analysing the earthquake-triggered boulder trails at Rapaki Bay, and that it can be usefully applied to rockfall trajectory and hazard assessment at this and similar sites elsewhere. • Detailed analysis of dynamic rockfall processes, interpreted from recorded boulder rolling experiments, and compared to RAMMS simulated results at the same site. Recorded rotational and translational velocities of a particular boulder show that the boulder behaves logically and dynamically on impact with different substrate types. Simulations show that seasonal changes in soil moisture alter rockfall dynamics and runout predictions within RAMMS, and adjustments are made to the calibration to reflect this; suggesting that in hazard analysis a rockfall model should be calibrated to dry rather than wet soil conditions to anticipate the most serious outcome. • Verifying the model calibration for a separate site on the Port Hills. The results of the RAMMS simulations show the effectiveness of calibration against a real data set, as well as the effectiveness of vegetation as a rockfall barrier/retardant. The results of simulations are compared using hazard maps, where the maximum runouts match well the mapped CES fallen boulder maximum runouts. The results of the simulations in terms of frequency distribution of deposit locations on the slope are also compared with those of the CES data, using the shadow angle tool to apportion slope zones. These results also replicate real field data well. Results show that a maximum runout envelope can be mapped, as well as frequency distribution of deposited boulders for hazard (and thus risk) analysis purposes. The accuracy of the rockfall runout envelope and frequency distribution can be improved by comprehensive vegetation and substrate mapping. The topics above define the scope of the project, limiting the focus to rockfall processes on the Port Hills, and implications for model calibration for the wider scientific community. The results provide a useful rockfall analysis methodology with a defensible and replicable calibration process, that has the potential to be applied to other lithologies and substrates. Its applications include a method of analysis for the selection and positioning of rockfall countermeasure design; site safety assessment for scaling and demolition works; and risk analysis and land planning for future construction in Christchurch.
This paper presents on-going challenges in the present paradigm shift of earthquakeinduced ground motion prediction from empirical to physics-based simulation methods. The 2010-2011 Canterbury and 2016 Kaikoura earthquakes are used to illustrate the predictive potential of the different methods. On-going efforts on simulation validation and theoretical developments are then presented, as well as the demands associated with the need for explicit consideration of modelling uncertainties. Finally, discussion is also given to the tools and databases needed for the efficient utilization of simulated ground motions both in specific engineering projects as well as for near-real-time impact assessment.
he strong motion station at Heathcote Valley School (HVSC) recorded unusually high peak ground accelerations (2.21g vertical and 1.41g horizontal) during the February 2011 Christchurch earthquake. Ground motions recorded at HVSC in numerous other events also exhibited consistently higher intensities compared with nearby strong motion stations. We investigated the underlying causes of such high intensity ground motions at HVSC by means of 2D dynamic finite element analyses, using recorded ground motions during the 2010-2011 Canterbury earthquake sequence. The model takes advantage of a LiDAR-based digital elevation model (DEM) to account for the surface topography, while the geometry and dynamic properties of the surficial soils are characterized by seismic cone penetration tests (sCPT) and Multi-Channel Analyses of Surface Waves (MASW). Comparisons of simulated and recorded ground motions suggests that our model performs well for distant events, while for near-field events, ground motions recorded at the adopted reference station at Lyttelton Port are not reasonable input motions for the simulation. The simulations suggest that Rayleigh waves generated at the inclined interface of the surficial colluvium and underlying volcanic rock strongly affect the ground motions recorded at HVSC, in particular, being the dominant contributor to the recorded vertical motions.
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.
Artificial Neural Networks (ANN) as a tool offers opportunities for modeling the inherent complexity and uncertainty associated with socio-environmental systems. This study draws on New Zealand ski
fields (multiple locations) as socio- environmental systems while considering their perceived resilience to low
probability but potential high consequences catastrophic natural events (specifically earthquakes). We gathered
data at several ski fields using a mixed methodology including: geomorphic assessment, qualitative interviews,
and an adaptation of Ozesmi and Ozesmi’s (2003) multi-step fuzzy cognitive mapping (FCM) approach. The data
gathered from FCM are qualitatively condensed, and aggregated to three different participant social groups. The
social groups include ski fields users, ski industry workers, and ski field managers. Both quantitative and
qualitative indices are used to analyze social cognitive maps to identify critical nodes for ANN simulations. The
simulations experiment with auto-associative neural networks for developing adaptive preparation, response and
recovery strategies. Moreover, simulations attempt to identify key priorities for preparation, response, and
recovery for improving resilience to earthquakes in these complex and dynamic environments. The novel mixed
methodology is presented as a means of linking physical and social sciences in high complexity, high uncertainty
socio-environmental systems. Simulation results indicate that participants perceived that increases in Social
Preparation Action, Social Preparation Resources, Social Response Action and Social Response Resources have
a positive benefit in improving the resilience to earthquakes of ski fields’ stakeholders.
This paper summarizes the development of a high-resolution surficial shear wave velocity model based on the combination of the large high-spatial-density database of cone penetration test (CPT) logs in and around Christchurch, New Zealand and a recently-developed Christchurch-specific empirical correlation between soil shear wave velocity and CPT. This near-surface shear wave velocity model has applications for site characterization efforts via the development of maps of time-averaged shear wave velocities over specific depths, as well as use in site response analysis and ground motion simulation.
This poster provides a summary of the development of a 3D shallow (z<40m) shear wave velocity (Vs) model for the urban Christchurch, New Zealand region. The model is based on a recently developed Christchurch-specific empirical correlation between Vs and cone penetration test (CPT) data (McGann et al. 2014a,b) and the large high-density database of CPT logs in the greater Christchurch urban area (> 15,000 logs as of 01/01/2014). In particular, the 3D model provides shear wave velocities for the surficial Springston Formation, Christchurch Formation, and Riccarton gravel layers which generally comprise the upper 40m in the Christchurch urban area. Point-estimates are provided on a 200m-by- 200m grid from which interpolation to other locations can be performed. This model has applications for future site characterization and numerical modeling efforts via maps of timeaveraged Vs over specific depths (e.g. Vs30, Vs10) and via the identification of typical Vs profiles for different regions and soil behaviour types within Christchurch. In addition, the Vs model can be used to constrain the near-surface velocities for the 3D seismic velocity model of the Canterbury basin (Lee et al. 2014) currently being developed for the purpose of broadband ground motion simulation.
