The abundance of cone penetration test (CPT) data from subsurface explorations in Christchurch and the surrounding areas provides a useful source of information for a characterization of the near surface shear wave velocity ( ) profile for the region. A portion of the investigations were conducted using seismic CPT, enabling the comparison of measured shear wave velocity with CPT data, and subsequently the evaluation of existing CPT- correlations for applicability to Canterbury-specific soils. The existing correlations are shown to be biased, generally over-predicting the observed with depth, thus demonstrating the need for a Canterbury-specific CPT- correlation.
Turning parts of Christchurch's red zoned land into a flat water sports lake, a community garden or a wave garden are some of the ideas being put forward to revitalise the earthquake-ravaged eastern part of the city.
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.
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.
Deformational properties of soil, in terms of modulus and damping, exert a great influence on seismic response of soil sites. However, these properties for sands containing some portion of fines particles have not been systematically addressed. In addition, simultaneous modelling of the modulus and damping behaviour of soils during cyclic loading is desirable. This study presents an experimental and computational investigation into the deformational properties of sands containing fines content in the context of site response analysis. The experimental investigation is carried on sandy soils sourced from Christchurch, New Zealand using a dynamic triaxial apparatus while the computational aspect is based on the framework of total-stress one-dimensional (1D) cyclic behaviour of soil. The experimental investigation focused on a systematic study on the deformational behaviour of sand with different amounts of fines content (particle diameter ≤ 75µm) under drained conditions. The silty sands were prepared by mixing clean sand with three different percentages of fines content. A series of bender element tests at small-strain range and stress-controlled dynamic triaxial tests at medium to high-strain ranges were conducted on samples of clean sand and silty sand. This allowed measurements of linear and nonlinear deformational properties of the same specimen for a wide strain range. The testing program was designed to quantify the effects of void ratio and fines content on the low-strain stiffness of the silty sand as well as on the nonlinear stress-strain relationship and corresponding shear modulus and damping properties as a function of cyclic shear strains. Shear wave velocity, Vs, and maximum shear modulus, Gmax, of silty sand was shown to be significantly smaller than the respective values for clean sands measured at the same void ratio, e, or same relative density, Dr. However, the test results showed that the difference in the level of nonlinearity between clean sand and silty sands was small. For loose samples prepared at an identical relative density, the behaviour of clean sand was slightly less nonlinear as compared to sandy soils with higher fines content. This difference in the nonlinear behaviour of clean sand and sandy soils was negligible for dense soils. Furthermore, no systematic influence of fines content on the material damping curve was observed for sands with fines content FC = 0 to 30%. In order to normalize the effects of fines on moduli of sands, equivalent granular void ratio, e*, was employed. This was done through quantifying the participation of fines content in the force transfer chain of the sand matrix. As such, a unified framework for modelling of the variability of shear wave velocity, Vs, (or shear modulus, Gmax) with void ratio was achieved for clean sands and sands with fines, irrespective of their fines content. Furthermore, modelling of the cyclic stress-strain behaviour based on this experimental program was investigated. The modelling effort focused on developing a simple constitutive model which simultaneously models the soil modulus and damping relationships with shear strains observed in laboratory tests. The backbone curve of the cyclic model was adopted based on a modified version of Kondner and Zelasko (MKZ) hyperbolic function, with a curvature coefficient, a. In order to simulate the hysteretic cycles, the conventional Masing rules (Pyke 1979) were revised. The parameter n, in the Masing’s criteria was assumed to be a function of material damping, h, measured in the laboratory. As such the modulus and damping produced by the numerical model could match the stress-strain behaviour observed in the laboratory over the course of this study. It was shown that the Masing parameter n, is strain-dependent and generally takes values of n ≤ 2. The model was then verified through element test simulations under different cyclic loadings. It was shown that the model could accurately simulate the modulus and the damping simultaneously. The model was then incorporated within the OpenSees computational platform and was used to scrutinize the effects of damping on one-dimensional seismic site response analysis. For this purpose, several strong motion stations which recorded the Canterbury earthquake sequence were selected. The soil profiles were modelled as semi-infinite horizontally layered deposits overlying a uniform half-space subjected to vertically propagating shear waves. The advantages and limitations of the nonlinear model in terms of simulating soil nonlinearity and associated material damping were further scrutinized. It was shown that generally, the conventional Masing criteria unconservatively may underestimate some response parameters such as spectral accelerations. This was shown to be due to larger hysteretic damping modelled by using conventional Masing criteria. In addition, maximum shear strains within the soil profiles were also computed smaller in comparison to the values calculated by the proposed model. Further analyses were performed to study the simulation of backbone curve beyond the strain ranges addressed in the experimental phase of this study. A key issue that was identified was that relying only on the modulus reduction curves to simulate the stress-strain behaviour of soil may not capture the actual soil strength at larger strains. Hence, strength properties of the soil layer should also be incorporated to accurately simulate the backbone curve.
Prognostic modelling provides an efficient means to analyse the coastal environment and provide effective knowledge for long term urban planning. This paper outlines how the use of SWAN and Xbeach numerical models within the ESRI ArcGIS interface can simulate geomorphological evolution through hydrodynamic forcing for the Greater Christchurch coastal environment. This research followed the data integration techniques of Silva and Taborda (2012) and utilises their beach morphological modelling tool (BeachMM tool). The statutory requirements outlined in the New Zealand Coastal Policy Statement 2010 were examined to determine whether these requirements are currently being complied with when applying the recent sea level rise predictions by the Intergovernmental Panel on Climate Change (2013), and it would appear that it does not meet those requirements. This is because coastal hazard risk has not been thoroughly quantified by the installation of the Canterbury Earthquake Recovery Authority (CERA) residential red zone. However, the Christchurch City Council’s (CCC) flood management area does provide an extent to which managed coastal retreat is a real option. This research assessed the effectiveness of the prognostic models, forecasted a coastline for 100 years from now, and simulated the physical effects of extreme events such as storm surge given these future predictions. 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 similar to the CCC’s flood management area. There are complex interactions at the Waimakariri River mouth with very high rates of accretion and erosion within a small spatial scale due to the river discharge. There is domination of the marine environment over the river system determined by the lack of generation of a distinct river delta, and river channel has not formed within the intertidal zone clearly. The Avon-Heathcote ebb tidal delta aggrades on the innner fan and erodes on the outer fan due to wave domination. The BeachMM tool facilitates the role of spatial and temporal analysis effectively and the efficiency of that performance is determined by the computational operating system.
