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.
The Canterbury earthquakes resulted in numerous changes to the waterways of Ōtautahi Christchurch. These included bank destabilisation, liquefaction effects, changes in bed levels, and associated effects on flow regimes and inundation levels. This study set out to determine if these effects had altered the location and pattern of sites utilised by īnanga (Galaxias maculatus) for spawning, which are typically restricted to very specific locations in upper estuarine areas. Extensive surveys were carried out in the Heathcote/Ōpāwaho and Avon/Ōtākaro catchments over the four peak months of the 2015 spawning season. New spawning sites were found in both rivers and analysis against pre-earthquake records identified that other significant changes have occurred. Major changes include the finding of many new spawning sites in the Heathcote/Ōpāwaho catchment. Sites now occur up to 1.5km further downstream than the previously reported limit and include the first records of spawning below the Woolston Cut. Spawning sites in the Avon/Ōtākaro catchment also occur in new locations. In the mainstem, sites now occur both upstream and downstream of all previously reported locations. A concentrated area of spawning was identified in Lake Kate Sheppard at a distinctly different location versus pre-quake records, and no spawning was found on the western shores. Spawning was also recorded for the first time in Anzac Creek, a nearby waterway connected to Lake Kate Sheppard via a series of culverts.
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.
