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.
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.
Heathcote Valley school strong motion station (HVSC) consistently recorded ground motions with higher intensities than nearby stations during the 2010-2011 Canterbury earthquakes. For example, as shown in Figure 1, for the 22 February 2011 Christchurch earthquake, peak ground acceleration at HVSC reached 1.4 g (horizontal) and 2 g (vertical), the largest ever recorded in New Zealand. Strong amplification of ground motions is expected at Heathcote Valley due to: 1) the high impedance contrast at the soil-rock interface, and 2) the interference of incident and surface waves within the valley. However, both conventional empirical ground motion prediction equations (GMPE) and the physics-based large scale ground motions simulations (with empirical site response) are ineffective in predicting such amplification due to their respective inherent limitations.
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.
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.
