Spatial variation of magnitude scaling factors during the 2010 Darfield and 2011 Christchurch,New Zealand,earthquakes

The combination of well-documented liquefaction response during the Darfield and Christchurch, New Zealand, earthquakes, densely-recorded ground motions for the events, and detailed subsurface characterization provides an unprecedented opportunity to investigate the significance of the spatial variation of magnitude scaling factors (MSF) on liquefaction triggering. Towards this end, MSF were computed at 15 SMS sites across Christchurch and its surroundings using two established approaches. Trends in the spatial variation of the MSF computed using number of equivalent cycles (n_eq) from both approaches were similar, with the spatial variation being more significant for the Christchurch earthquake than the Darfield earthquake. However, there was no consistent trend for regions with lower computed MSF having experienced more severe or widespread liquefaction. Additionally, there is a general correlation between MSF and a_max, but because a_max ranges more widely than MSF it has a greater influence on the resulting seismic demand imposed on the soil than MSF does. Nevertheless, the spatial variation of the MSF is deemed significant enough that it warrants being considered for incorporation into future variants of simplified liquefaction evaluation procedures.     

Observations and implications of damage from the magnitude Mw 6.3 Christchurch, New Zealand earthquake of 22 February 2011

This paper describes the observations made by a reconnaissance team following the 22nd February 2011, M_w 6.3, Christchurch, New Zealand earthquake (GNS Science, 2011). The team comprised of members of the UK based Earthquake Engineering Field Investigation Team who spent 5 days collecting observations on damage resulting from the earthquake. Although the magnitude of this earthquake was not particularly high (M_w 6.3), the shallow focus and close proximity resulted in locally very high ground motions, widespread damage and 182 fatalities. The earthquake is also particularly notable for the widespread liquefaction it caused, landslides and rockfalls in the hills south of Christchurch, and the significant damage suffered by unreinforced masonry and historic structures. Over wide areas of central Christchurch, recorded accelerations were in excess of those required by the current New Zealand seismic loadings standard (NZS1170.5:2004: Standards New Zealand 2004) and therefore the earthquake presented a valuable opportunity to assess performance of modern buildings under code-level ground acceleration.     

Applicability of existing empirical shear wave velocity correlations to seismic cone penetration test data in Christchurch New Zealand

Seismic piezocone (SCPTu) data compiled from 86 sites in the greater Christchurch, New Zealand area are used to evaluate several existing empirical correlations for predicting shear wave velocity from cone penetration test (CPT) data. It is shown that all the considered prediction models are biased towards overestimation of the shear wave velocity of the Christchurch soil deposits, demonstrating the need for a Christchurch-specific shear wave velocity prediction model (McGann et al., 2014) [1]. It is hypothesized that the unique depositional environment of the considered soils and the potential loss of soil ageing effects brought about by the 2010–2011 Canterbury earthquake sequence are the primary source of the observed prediction bias.     

Entanglement of New Zealand fur seals in man-made debris at Kaikoura, New Zealand

New Zealand fur seals in the Kaikoura region breed near a town with expanding tourist and fishing industries and commonly come ashore entangled in nets and plastic debris. However, the rate at which entanglement occurs was previously unknown. A decade of Department of Conservation seal callout data was analysed to determine the level of entanglement in the region and the most common debris type. Monitoring of adult female fur seals released from entanglement provided information on the potential for serious wounds to heal and survivorship of released individuals. Entanglement rates of pinnipeds in Kaikoura are some of the highest reported world-wide (average range: 0.6-2.8%) with green trawl net (42%), and plastic strapping tape (31%) together contributing the most to debris types. Nearly half of the reported entangled seals are successfully released (43%) and post-release monitoring shows that with appropriate intervention the chance of an individual surviving even with a significant entanglement wound is high. Our study demonstrates that while entanglement in the region is high, a successful intervention protocol may help reduce the potential for entanglement-related mortality in the region.     

Effect of shear strain compatibility and incompatibility approaches in the design of high modulus columns against liquefaction: A case study in Christchurch,New Zealand

Bulletin of Earthquake Engineering - Nowadays, investigating the effectiveness of high modulus columns in liquefaction mitigation is one of the important tasks in earthquake geotechnical...     

Magnetic anomalies across Campbell Plateau,New Zealand

Magnetic measurements over the Campbell Plateau, New Zealand, show the existence of a linear positive anomaly belt which extends for 900 km along an approximately east-west trend. This anomaly belt is considered to be a continuation of the Stokes Magnetic Anomaly occurring through New Zealand and associated with rocks of the New Zealand Geosyncline. If this is so, a transcurrent fault with a dextral displacement of about 330 km and aligned in a northeast-southwest direction must occur between the eastern end of the Stokes Magnetic Anomaly and the western end of the anomaly belt across Campbell Plateau.     

Observations on white island volcano,New Zealand

White Island is a complex of two overlapping cones constructed of lava flows, agglomerates and unconsolidated and unsorted ash and tuff beds. Remnants of a welded-tuff flow have been found on the north-east flank of the volcano. Since the extrusion of the youngest lava flow the young cone has been breached to the south-east and deeply eroded. White Island lavas are porphyritic augite-hypersthene-labradorite andesites. One young lava flow is unusually rich in Na₂O and contains groundmass sodian ferroaugite instead of the normal augite and hypersthene. The unusual groundmass features of this andesite are believed to be the result of contamination. Volcanic, plutonic and gneissic xenoliths have been found in the White Island lavas. Three new analyses of White Island andesites are given together with an electron microprobe analysis of a groundmass glass from one of the andesites. The White Island andesites are believed to have formed from the hybridisation of a primary mantle-derived andesitic magma with crustal material below the base of the Mesozoic New Zealand Geosyncline.     

Thermal modeling of the Southern Alps,New Zealand

Finite-element modeling of the thermal regime across the Southern Alps of New Zealand has been carried out along two profiles situated near the Franz Josef and Haast valleys. The modeling involves viscous deformation beneath the Southern Alps, including both uplift and erosion, and crustal/lithospheric thickening, as a result of crustal shortening extending to 20 mm/y of a 25-km thick crust. Published uplift rates and crustal thickness variations along the two profiles are used to constrain the modeled advection of crustal material, and results are compared with the recent heat flow determinations, 190±50 mW/m² in the Franz Josef valley and 90±25 mW/m² in the Haast valley. Comparisons of the model with published K–Ar and fission track ages, show that the observed heat flow in the Franz Josef valley is consistent with observed zircon fission track ages of around 1 Ma, if the present-day uplift rate is close to 10 mm/y. Major thermal differences between the Franz Josef and Haast profiles appear to be due to different uplift and erosion rates. There is weak evidence that frictional heating close to the Alpine fault zone is not significant. The modeling provides explanations for the distribution of seismicity beneath the Southern Alps, and predicts a low surface heat flow over the eastern foothills due to the dominant thermal effect of crustal thickening beneath this region. Predicted temperatures at mid-crustal depth beneath the zone of maximum uplift rate are 50–100°C cooler than those indicated in previously published models, which implies that thermal weakening of the crust may not be the main factor causing the aseismicity of the central Southern Alps. The results of the modeling demonstrate that the different types of reset age data in the region within 25 km of the Alpine fault are critical for constraining models of the deformation and the thermal regime beneath the Southern Alps.     

Geomorphic evolution of the Southern Alps,New Zealand

The Southern Alps are the topographic expression of late Cenozoic (<8 Ma ago) uplift of the crust of the leading edge of the Pacific plate in South Island, New Zealand. New fission track data on the basement exposed in the Southern Alps quantify the age, amount, and rate of rock uplift, and in combination with geomorphic parameters permit the construction of a new model of the geomorphic evolution of the Southern Alps. The model emphasizes the development over time and space of rock uplift, mean surface elevation, exhumation of crustal section, and relief. The earliest indications of mean surface uplift are between 4 and 5 Ma ago at the Alpine Fault. Mean surface uplift, which lagged the start of rock uplift, propagated southeastward from the Alpine Fault at a rate of 30 km/Ma. By about 4 Ma ago, exhumation had exposed greywacke basement adjacent to and east of the entire 300 km long central section of the Alpine Fault. At 3 Ma ago, greenschist was exposed in the southern parts of the Southern Alps near Lake Wanaka, and since then has become exhumed along a narrow strip east of the Alpine Fault. The model infers that amphibolite grade schist has been exhumed adjacent to the Alpine Fault only in the last 0·3 Ma. The age of the start of rock uplift and the amount and rate of rock uplift, all of which vary spatially, are considered to be the dominant influences on the development of the landscape in the Southern Alps. The Southern Alps have been studied in terms of domains of different rock uplift rate. At present the rate of rock uplift varies from up to 8–10 mm/a adjacent to the Alpine Fault to 0·8–1·0 mm/a along the southeastern margin of the Southern Alps. This spectrum can be divided into two domains, one northwest of the Main Divide where the present rock uplift rates are very high (up to 8–10 mm/a) and exceed the long-term value of 0·8–1·0 mm/a, and another to the southeast of the Main Divide where the long-term rate is 0·8–1·0 mm/a. A domain of no uplift lies immediately to the east of the Southern Alps, and is separated from them by a 1·0–1·5 km step in the basement topography. We argue that this spatial sequence of uplift rate domains represents a temporal one. The existing models of the geomorphic development of the Southern Alps—the dynamic cuesta model of J. Adams and the numerical model of P. Koons—are compared with the new data and evolutionary model. Particular constraints unrealized by these two earlier models include the following: the earlier timing of the start of rock uplift of the Southern Alps (8 Ma ago); the spatial variation in the timing of the start of rock uplift (8 Ma ago to 3 Ma ago); the lower long-term rock uplift rate (0·8–1·0 mm/a) of the Southern Alps for most of the late Cenozoic; the lag between the start of rock uplift and the start of mean surface uplift; and the patterns of the amounts of late Cenozoic rock uplift and erosion across the Southern Alps.     

Origin of dacites of Taupo Volcanic Zone, New Zealand

Dacites form a relatively small proportion of lavas in the Taupo Volcanic Zone, New Zealand (5km³), and occur mainly on the eastern side. In this paper their origin is considered in terms of three processes: (a) partial melting of crustal rocks; (b) fractional crystallisation of basalt and andesite; and (c) sub-surface mixing of basic and acid magma. Modelling techniques are used to calculate the most acceptable degree of fractional crystallisation and magma mixing to fit major-element data, and these values are used to compare calculated and observed trace-element values. The success or failure of the model is determined by the closeness of the two sets of values. For partial-melt models, trace-element values alone are calculated by the batch-melting equation.Results indicate that White Island dacite can best be modelled by fractional crystallisation; Manawahe by fractional crystallisation plus limited crustal contamination; Maungaongaonga by partial melting of Western Basement greywacke, and Tauhara by partial melting of this greywacke together with minor mixing with a more basic magma. Results from Parekauau and Horohoro indicate that these lavas are unlikely to have formed by any of the processes examined.     

Elliptical boundary of an intraplate volcanic field, Auckland, New Zealand

The outer vents in the Auckland Volcanic Field lie within 19 to 559 m (mean 224.75 m) of a ‘best-fit' ellipse with a 28.9-km-long major axis trending almost north–south, and a minor axis 16.5 km long. The ellipse has formed the outer boundary of the field since the inception of volcanism 140,000 years ago. We present the following testable hypotheses as an explanation of this pattern: The boundary is the expression of a corresponding elliptical source area at depth in the lithospheric mantle (possibly asthenosphere material trapped at this level). The ellipse may represent a depth contour on a very small upper mantle dome or of a lens intruding into a neck of an extensional structure. Alternatively it could be the boundary of a flat elliptical area where tensional stresses allow decompressional melting. The elliptical tensional region may either have developed in a releasing bend during strike-slip faulting along a fundamental lithosphere structure inherited from Mesozoic tectonics, associated with the NNW-trending Dun Mountain ophiolite belt, or may represent the tip of a fracture along which the Auckland Volcanic Province is propagating northward.     

Semi-determinate hydraulic geometry of river channels,South Island,New Zealand

Data to describe the morphologic, hydrologic and sedimentologic characteristics of 72 South Island, New Zealand, rivers were collected and analysed. Nearly 70 per cent of variation in channel morphology is accounted for by differences in cross-sectional area, slope, and cross-section shape; only 53 per cent of the morphologic variability could be statistically ‘explained’ by the hydrologic and sediment variables used. The level of explanation varied for different morphologic variables; nearly 90 per cent of the variability in cross-sectional area could be explained, but aspect ratio (maximum depth divided by hydraulic radius) was completely independent. Apart from the inadequacy of the measured variables as indices of the true underlying controlling factors, and the imperfect measurement and sampling procedures, the low level of explanation is probably due to the influence of factors such as floodplain vegetation, high quasi-random variability in bark sediment character, boundary effects imposed by bedrock bluffs, and the precise sequence of flood events, none of which are easily quantified. In addition, observations indicate that there is a large random variation in channel form which cannot be related to any factor. An attempt to relate channel morphology to flow variability, using simple indices of the latter, was unsuccessful.     

Magnitude reassessment of New Zealand earthquakes

The surface-wave magnitudes of a selection of New Zealand earthquakes have been determined on a consistent basis using the ‘Prague formula’ and station corrections. The earthquakes range in magnitude from about 5 to 7.8, covering the instrumental period 1901–1988. Magnitudes for many of the earlier events had not been properly determined previously; and some significant discrepancies from the traditional magnitudes were found. The use of European station data (160° < D < 175°) is important to New Zealand because of its geographical isolation. These distant data were found to give consistently slightly higher M_s than closer stations, but could be used without bias through the station correction procedure. The relationship between M_s and M_L was found for 31 ‘shallow’ New Zealand events and much of the scatter was explained as a function of depth. Significant differences in M_s/M_L expressions from Europe and California were also found. The limited New Zealand data for M_w and M₀ related well to Californian and global relationships with M_s.     

Forecasting tsunamis in Poverty Bay, New Zealand, with deep-ocean gauges

The response/transfer function of a coastal site to a remote open-ocean point is introduced, with the intent to directly convert open-ocean measurements into the wave time history at the site. We show that the tsunami wave at the site can be predicted as the wave is measured in the open ocean as far as 1,000+ km away from the site, with a straightforward computation which can be performed almost instantaneously. The suggested formalism is demonstrated for the purpose of tsunami forecasting in Poverty Bay, in the Gisborne region of New Zealand. Directional sensitivity of the site response due to different conditions for the excitation of the shelf and the bay’s normal modes is investigated and used to explain tsunami observations. The suggested response function formalism is validated with available records of the 2010 Chilean tsunami at Gisborne tide gauge and at the nearby deep-ocean assessment and reporting of tsunamis (DART) station 54401. The suggested technique is also demonstrated by hindcasting the 2011 Tohoku tsunami and 2012 Haida Gwaii tsunami at Monterey Bay, CA, using an offshore record of each tsunami at DART station 46411.     

Timing of volcanic, plutonic and geothermal activity at Ngatamariki, New Zealand

Four ⁴⁰Ar/³⁹Ar dates on mineral separates from fresh and hydrothermally altered volcanic and plutonic rocks from the Ngatamariki geothermal field indicate that andesitic volcanism took place in the eastern portion of the Taupo Volcanic Zone (TVZ) prior to 1.2 Ma and probably considerably earlier. These data significantly extend the onset and duration of andesitic volcanism in the east-central TVZ over previous estimates. Intrusive activity is represented at Ngatamariki by a dioritic pluton, the only such pluton yet recognized in the entire TVZ. Hornblende from the pluton yields a crystallization age of near 550 ka. Hydrothermal alteration spatially associated with the pluton produced sericite of a similar age. Overlying and postdating the most intense hydrothermal alteration zone is the Whakamaru Ignimbrite (or its equivalent) which was emplaced at 330 ka. Two distinct geothermal systems may have been active at nearly the same site from 550 ka to present. The most intense activity occurred before 330 ka and was associated with emplacement of the Ngatamariki diorite. This was followed by the less intense system that is currently active. The geothermal regime at Ngatamariki has, therefore, probably been active intermittently for at least 550 ka.     

Revised palaeoclimatic significance of Mueller Glacier moraines,Southern Alps,New Zealand

A sequence of Late Holocene moraines on the foreland of the Mueller Glacier, Southern Alps, New Zealand, forms part of a local moraine‐age database used to establish a regional glacier chronology and subsequently to investigate potential intra‐hemispheric and global climate forcing mechanisms. We present new sedimentological and geomorphological evidence that a set of these moraine ridges, previously considered to represent individual advances, constitutes a single moraine complex (the ‘Mueller Memorial Moraine’) formed by supraglacial transport of a large volume of landslide debris to the glacier terminus. Because a moraine formed in this way is not necessarily associated with an advance triggered by a climate event, we question the palaeoclimatic significance of the Mueller Memorial Moraine, as well as that of other moraines in comparable settings. Our findings suggest that the mode of formation and glacio‐dynamical context of moraines whose ages contribute to existing palaeoclimate reconstructions need to be re‐examined in order to assess the reliability of these reconstructions. Copyright © 2015 John Wiley & Sons, Ltd.     

Forecasting tsunamis in Poverty Bay,New Zealand,with deep-ocean gauges

The response/transfer function of a coastal site to a remote open-ocean point is introduced, with the intent to directly convert open-ocean measurements into the wave time history at the site. We show that the tsunami wave at the site can be predicted as the wave is measured in the open ocean as far as 1,000+ km away from the site, with a straightforward computation which can be performed almost instantaneously. The suggested formalism is demonstrated for the purpose of tsunami forecasting in Poverty Bay, in the Gisborne region of New Zealand. Directional sensitivity of the site response due to different conditions for the excitation of the shelf and the bay’s normal modes is investigated and used to explain tsunami observations. The suggested response function formalism is validated with available records of the 2010 Chilean tsunami at Gisborne tide gauge and at the nearby deep-ocean assessment and reporting of tsunamis (DART) station 54401. The suggested technique is also demonstrated by hindcasting the 2011 Tohoku tsunami and 2012 Haida Gwaii tsunami at Monterey Bay, CA, using an offshore record of each tsunami at DART station 46411.

     

Lacustrine shore platforms at Lake Waikaremoana,North Island,New Zealand

This paper examines the morphology and processes governing the development of shore platforms at Lake Waikaremoana, North Island, New Zealand. Shore platforms at Lake Waikaremoana are recent features, and were formed when a new sequence of shoreline development was initiated, due to lowering of the lake by 5 m in 1946 for hydroelectric power development. Three predominant platform morphologies were identified around the lake. These include gently sloping platforms (c.1·5 to 3·9°), ramp platforms (c.6·8 to 9·2°), and concave ramp platforms (c.7·9 to 12°). Platform widths ranged from 11 to 31 m, with the gently sloping platforms characterized by the widest morphologies. Erosion rates were estimated using perched sandstone boulders and were found to range from 3·4 to 12·5 mm a^?1, with a mean erosion rate of 5·9 mm a^?1. Higher rates of erosion were identified at lower platform elevations, due to a greater frequency of wetting and drying cycles coincident with storm waves, while lower erosion rates were identified at higher elevations. Field evidence suggests that shore platforms at Lake Waikaremoana were likely initiated and continue to develop as a result of subaerial wetting and drying cycles. Waves, coincident with fluctuating lake levels, play an important role by removing the weathered material from the platforms, and appear to control the width of the platforms. A conceptual model of platform development is presented, and analogies are drawn between this model, and the formation of shore platforms in oceanic environments. Copyright © 2002 John Wiley & Sons, Ltd.