scholarly journals Empirical models for predicting lateral spreading and evaluation using New Zealand data

Author(s):  
Jian Zhang ◽  
Dick Beetham ◽  
Grant Dellow ◽  
John X. Zhao ◽  
Graeme H. McVerry

A New empirical model has been developed for predicting liquefaction-induced lateral spreading displacement and is a function of response spectral displacements and geotechnical parameters. Different from the earlier model of Zhang and Zhao (2005), the application of which was limited to Japan and California, the new model can potentially be applied anywhere if ground shaking can be estimated (by using local strong-motion attenuation relations). The new model is applied in New Zealand where the response spectral displacement is estimated using New Zealand strong-motion attenuation relations (McVerry et al. 2006). The accuracy of the new model is evaluated by comparing predicted lateral displacements with those which have been measured from aerial photos or the width of ground cracks at the Landing Road bridge, the James Street loop, the Whakatane Pony Club and the Edgecumbe road and rail bridges sites after the 1987 Edgecumbe earthquake. Results show that most predicted errors (defined as the ratio of the difference between the measured and predicted lateral displacements to the measured one) from the new model are less than 40%. When compared with earlier models (Youd et al. 2002, Zhang and Zhao 2005), the new model provides the lowest mean errors.

Author(s):  
G. H. McVerry

Probabilistic techniques for seismic hazard analysis have
come into vogue in New Zealand for both the assessment of major projects and the development and review of seismic design codes. However, there are considerable uncertainties in the modelling
 of the strong-motion attenuation, which is necessarily based largely on overseas data. An excellent agreement is obtained between an average 5% damped response spectrum for New Zealand alluvial sites in the 20 to 59 km distance range and 5.4 to 6.0 magnitude class and that given by a Japanese model. Unfortunately, this corresponds to only about half the amplitude levels of 150 year spectra relevant to code design. The much more rapid decay
of ground shaking with distance in New Zealand has led to a considerable modification based on maximum ground acceleration
data from the Inangahua earthquake of the distance-dependence
of the Japanese response spectra model. Less scatter in New Zealand data has resulted in adopting a lower standard deviation for the attenuation model, which is important in reducing the considerable "probabilistic enhancement" of the hazard estimates. Regional differences in attenuation shown by intensities are difficult to resolve from the strong-motion acceleration data, apart from lower accelerations in Fiordland.


Author(s):  
Mark E. Stringer ◽  
Sarah Bastin ◽  
Christopher R. McGann ◽  
Claudio Cappellaro ◽  
Maya El Kortbawi ◽  
...  

The magnitude Mw7.8 ‘Kaikōura’ earthquake occurred shortly after midnight on 14 November 2016. This paper presents an overview of the geotechnical impacts on the South Island of New Zealand recorded during the post-event reconnaissance. Despite the large moment magnitude of this earthquake, relatively little liquefaction was observed across the South Island, with the only severe manifestation occurring in the young, loose alluvial deposits in the floodplains of the Wairau and Opaoa Rivers near Blenheim. The spatial extent and volume of liquefaction ejecta across South Island is significantly less than that observed in Christchurch during the 2010-2011 Canterbury Earthquake Sequence, and the impact of its occurrence to the built environment was largely negligible on account of the severe manifestations occurring away from the areas of major development. Large localised lateral displacements occurred in Kaikōura around Lyell Creek. The soft fine-grained material in the upper portions of the soil profile and the free face at the creek channel were responsible for the accumulation of displacement during the ground shaking. These movements had severely impacted the houses which were built close (within the zone of large displacement) to Lyell Creek. The wastewater treatment facility located just north of Kaikōura also suffered tears in the liners of the oxidation ponds and distortions in the aeration system due to ground movements. Ground failures on the Amuri and Emu Plains (within the Waiau Valley) were small considering the large peak accelerations (in excess of 1g) experienced in the area. Minor to moderate lateral spreading and ejecta was observed at some bridge crossings in the area. However, most of the structural damage sustained by the bridges was a result of the inertial loading, and the damage resulting from geotechnical issues were secondary.


Author(s):  
L. S. Hogan ◽  
L. M. Wotherspoon ◽  
J. M. Ingham

During seismic assessments of bridges where there is a lack of construction documentation, one method of determining likely structural detailing is to use historic design standards. An overview of the New Zealand bridge seismic standards and the agencies that have historically controlled bridge design and construction is presented. Standards are grouped into design era based upon similar design and loading characteristics. Major changes in base shear demand, ductility, foundation design, and linkage systems are discussed for each design era, and loadings and detailing requirements from different eras were compared to current design practices. Bridges constructed using early seismic standards were designed to a significantly lower base shear than is currently used but the majority of these bridges are unlikely to collapse due to their geometry and a preference for monolithic construction. Bridges constructed after the late 1970s are expected to perform well if subjected to ground shaking, but unless bridges were constructed recently their performance when subjected to liquefaction and liquefaction-induced lateral spreading is expected to be poor.


Author(s):  
Graeme H. McVerry ◽  
John X. Zhao ◽  
Norman A. Abrahamson ◽  
Paul G. Somerville

Attenuation relations are presented for peak ground accelerations (pga) and 5% damped acceleration response spectra in New Zealand earthquakes. Expressions are given for both the larger and the geometric mean of two randomly-oriented but orthogonal horizontal components of motion. The relations take account of the different tectonic types of earthquakes in New Zealand, i.e., crustal, subduction interface and dipping slab, and of the different source mechanisms for crustal earthquakes. They also model the faster attenuation of high-frequency earthquake ground motions in the volcanic region than elsewhere. Both the crustal and subduction zone attenuation expressions have been obtained by modifying overseas models for each of these tectonic environments to better match New Zealand data, and to cover site classes that relate directly to those used for seismic design in New Zealand codes. The study used all available data from the New Zealand strong-motion earthquake accelerograph network up to the end of 1995 that satisfied various selection criteria, supplemented by selected data from digital seismographs. The seismographs provided additional records from rock sites, and of motions involving propagation paths through the volcanic region, classes of data that are sparse in records produced by the accelerograph network. The New Zealand strong-motion dataset lacks records in the nearsource region, with only one record from a distance of less than 10 km from the source, and at magnitudes greater than Mw 7.23. The New Zealand data used in the regression analyses ranged in source distance from 6 km to 400 km (the selected cutoff) and in moment magnitude from 5.08 to 7.23 for pga, with the maximum magnitude reducing to 7.09 for response spectra data. The required near-source constraint has been obtained by supplementing the New Zealand dataset with overseas peak ground acceleration data (but not response spectra) recorded at distances less than 10 km from the source. Further near-source constraints were obtained from the overseas attenuation models, in terms of relationships that had to be maintained between various coefficients that control the estimated motions at short distances. Other coefficients were fitted from regression analyses to better match the New Zealand data. The need for different treatment of crustal and subduction zone earthquakes is most apparent when the effects or source mechanism are taken into account. For crustal earthquakes, reverse mechanism events produce the strongest motions, followed by strike-slip and normal events. For subduction zone events, the reverse mechanism interface events have the lowest motions, at least in the period range up to about ls, while the slab events, usually with normal mechanisms, are generally strongest. The attenuation relations presented in this paper have been used in many hazard studies in New Zealand over the last five years. In particular, they have been used in the derivation of the elastic site spectra in the new Standard for earthquake loads in New Zealand, NZS 1170.5:2004.


Author(s):  
Aasha Pancha ◽  
John Taber

Attenuation relations using weak ground motion recordings have been determined using data from the New Zealand National Seismograph Network and several temporary seismograph deployments. Models have been developed for earthquake sources in four regions: the Eastern North Island deep and shallow regions and the Central North Island (CNI) deep and shallow regions. Deep events were those with hypocenters below 33 km. Regression coefficients have been determined using the attenuation models of Joyner and Boore (1981) and Molas and Yamazaki (1995). The anelastic attenuation rates in the Eastern North Island expressions are comparable to that of Joyner and Boore (1981), suggesting that weak motion attenuation can be used to estimate variations in strong motion attenuation. However, the absolute level of the strong-motion attenuation curves greatly differs from those of the weak-motion. The anelastic attenuation rate for the shallow CNI is of the order of two to three times that observed for the Eastern North Island. The lowest attenuation rate was found for events within the deep CNI, whose ray paths did not cross the shallow Central North Island region. This is consistent with a low rate of attenuation in the subducting Pacific plate. Azimuthal dependence of PGA is evident within each of the regions. Within the Eastern North Island, the attenuation rate is lowest in the direction of 35-55° from North, which is roughly along the strike of the subducting Pacific plate. A similar azimuthal dependence was also noted within the deep CNI region, while a slightly different minimum direction (5°) was determined for the shallow CNI region.


Author(s):  
J. B. Berrill

The principal aim of the present network of strong motion accelerographs is to record the response of structures to earthquakes, and instruments are concentrated in the larger cities where modern, tall buildings are found. However, the behaviour of structures during earthquakes is now comparatively well understood. At the present time, estimating design ground motions is the weakest part in the process of designing structures to resist earthquakes. There is a strong need for more recordings of ground shaking, particularly
sets of several accelerograms from single earthquakes. It is not certain that the present accelerograph network would capture any significant record of strong motion during a major earthquake in New Zealand; and the chance of a set of three or more strong accelerograms being recorded is quite small. It is recommended that 25 additional instruments be installed promptly, to fill the main gaps in the present network, and to extend the capacity of the existing local network in the Wellington area.


Author(s):  
D. J. Dowrick ◽  
D. A. Rhoades

This paper presents a new approach to modelling the spatial distribution of intensities in crustal earthquakes, using a distributed source. The source is represented by one or two rectangular fault rupture planes of chosen dip, discretised into small rectangles each with its own share of the total seismic moment, and modelling chosen distributions of asperities. The Modified Mercalli (MM) intensity of shaking is represented by isoseismals. Comparisons are made with the actual isoseismals (particularly of intensities MM9 and MM10) of selected large historical crustal New Zealand earthquakes and those predicted by the simpler models of Dowrick & Rhoades (2005). Important differences and insights are found regarding near-source spatial distributions of ground shaking of shallow earthquakes with rupture length greater than about 28 km (Mw > 6.8) with any dip, and for Mw > c. 5.5 with dip < 60º. The influence of asperities relative to that of non-asperities is seen as modest near-fault increases in intensity. The new model can be applied to planar or biplanar fault ruptures of any length, width and dip. In the absence of isoseismal data on large earthquakes with normal focal mechanisms the current model is only verified for use on strike-slip and reverse events. A new concept, seismic-source intensity, is introduced and utilized. The new model can also be applied to earthquakes in other regions of the world with adjustments for local attenuation rates as necessary.


Author(s):  
J. H. Wood ◽  
G. R. Martin

Investigations carried out to evaluate the seismic design parameters, including acceleration response spectra and time-histories, for the design of the Museum of New Zealand, Te Papa Tongarewa, on the Wellington waterfront are described. The procedures used to assess the site stability under strong ground shaking and to determine the maximum likely lateral spreading and settlements are also summarised.


1981 ◽  
Vol 71 (6) ◽  
pp. 2011-2038 ◽  
Author(s):  
William B. Joyner ◽  
David M. Boore

Abstract We have taken advantage of the recent increase in strong-motion data at close distances to derive new attenuation relations for peak horizontal acceleration and velocity. This new analysis uses a magnitude-independent shape, based on geometrical spreading and anelastic attenuation, for the attenuation curve. An innovation in technique is introduced that decouples the determination of the distance dependence of the data from the magnitude dependence. The resulting equations are log A = − 1.02 + 0.249 M − log r − 0.00255 r + 0.26 P r = ( d 2 + 7.3 2 ) 1 / 2 5.0 ≦ M ≦ 7.7 log V = − 0.67 + 0.489 M − log r − 0.00256 r + 0.17 S + 0.22 P r = ( d 2 + 4.0 2 ) 1 / 2 5.3 ≦ M ≦ 7.4 where A is peak horizontal acceleration in g, V is peak horizontal velocity in cm/ sec, M is moment magnitude, d is the closest distance to the surface projection of the fault rupture in km, S takes on the value of zero at rock sites and one at soil sites, and P is zero for 50 percentile values and one for 84 percentile values. We considered a magnitude-dependent shape, but we find no basis for it in the data; we have adopted the magnitude-independent shape because it requires fewer parameters.


1995 ◽  
Vol 85 (6) ◽  
pp. 1821-1834
Author(s):  
Toshimi Satoh ◽  
Toshiaki Sato ◽  
Hiroshi Kawase

Abstract We evaluate the nonlinear behavior of soil sediments during strong ground shaking based on the identification of their S-wave velocities and damping factors for both the weak and strong motions observed on the surface and in a borehole at Kuno in the Ashigara Valley, Japan. First we calculate spectral ratios between the surface station KS2 and the borehole station KD2 at 97.6 m below the surface for the main part of weak and strong motions. The predominant period for the strong motion is apparently longer than those for the weak motions. This fact suggests the nonlinearity of soil during the strong ground shaking. To quantify the nonlinear behavior of soil sediments, we identify their S-wave velocities and damping factors by minimizing the residual between the observed spectral ratio and the theoretical amplification factor calculated from the one-dimensional wave propagation theory. The S-wave velocity and the damping factor h (≈(2Q)−1) of the surface alluvial layer identified from the main part of the strong motion are about 10% smaller and 50% greater, respectively, than those identified from weak motions. The relationships between the effective shear strain (=65% of the maximum shear strain) calculated from the one-dimensional wave propagation theory and the shear modulus reduction ratios or the damping factors estimated by the identification method agree well with the laboratory test results. We also confirm that the soil model identified from a weak motion overestimates the observed strong motion at KS2, while that identified from the strong motion reproduces the observed. Thus, we conclude that the main part of the strong motion, whose maximum acceleration at KS2 is 220 cm/sec2 and whose duration is 3 sec, has the potential of making the surface soil nonlinear at an effective shear strain on the order of 0.1%. The S-wave velocity in the surface alluvial layer identified from the part just after the main part of the strong motion is close to that identified from weak motions. This result suggests that the shear modulus recovers quickly as the shear strain level decreases.


Sign in / Sign up

Export Citation Format

Share Document