Search USGS⌕ Search

Geology topics

Research about New Zealand

Source-linked reports with geographic coverage including New Zealand.

At least 91 records · Page 5Linked to original sources

New Zealand’s deadliest quake sounds alarm for cities on fault lines

The catastrophic Christ Church Earthquake is a strong reminder to engineers and scientists of the hazards pose by fault lines, both mapped and unknown, near major cities. In February 2011, the relatively moderate earthquake that struck the cities of Christchurch and Lyttleton in the Canterbury region of New Zealand's South Island surprised many with its destructive power. The magnitude 6.2 temblor killed 181 people, 118 of whom were killed in the collapse of a single building in the city center. The quake damaged or destroyed more than 100,000 buildings. It was the deadliest quake to strike the nation in 80 years-since the 1931 earthquake that struck the Napier and Hastings area of the North Island. The Christchurch quake was part of the aftershock sequence following the September 2010 magnitude 7.1 earthquake near Darfield, 40 kilometers west of the city. The Darfield earthquake was in a sparsely populated area, causing to loss of life. By contrast, the Christchurch earthquake was generated on a fault in close proximity to the city.

Christchurch↗

Distinguishing between stress-induced and structural anisotropy at Mount Ruapehu volcano, New Zealand

We have created a benchmark of spatial variations in shear wave anisotropy around Mount Ruapehu, New Zealand, against which to measure future temporal changes. Anisotropy in the crust is often assumed to be caused by stress-aligned microcracks, and the polarization of the fast quasi-shear wave ( ϕ ) is thus interpreted to indicate the direction of maximum horizontal stress, but can also be due to aligned minerals or macroscopic fractures. Changes in seismic anisotropy have been observed following a major eruption in 1995/96 and were attributed to changes in stress from the depressurization of the magmatic system. Three-component broadband seismometers have been deployed to complement the permanent stations that surround Ruapehu, creating a combined network of 34 three-component seismometers. This denser observational network improves the resolution with which spatial variations in seismic anisotropy can be examined. Using an automated shear wave splitting analysis, we examine local earthquakes in 2008. We observe a strong azimuthal dependence of ϕ and so introduce a spatial averaging technique and two-dimensional tomography of recorded delay times. The anisotropy can be divided into regions in which ϕ agrees with stress estimations from focal mechanism inversions, suggesting stress-induced anisotropy, and those in which ϕ is aligned with structural features such as faults, suggesting structural anisotropy. The pattern of anisotropy that is inferred to be stress related cannot be modeled adequately using Coulomb modeling with a dike-like inflation source. We suggest that the stress-induced anisotropy is affected by loading of the volcano and a lithospheric discontinuity.

Journal of Geophysical Research B: Solid Earth↗

Preface to the Focused Issue on the 22 February 2011 Magnitude 6.2 Christchurch Earthquake

The 22 February 2011 magnitude 6.2 Christchurch earthquake, centered southeast of Christchurch, was part of the aftershock sequence that has been occurring since the September 2010 magnitude 7.1 quake near Darfield, 40 km west of the city. The Christchurch earthquake killed more than 180 people, damaged or destroyed more than 100,000 buildings, and is New Zealand's most deadly disaster since the earthquake that struck the Napier and Hastings area on 3 February 1931. This special focused issue of Seismological Research Letters , which I had the fortune to edit, contains a selected set of 19 original technical papers. These papers cover different aspects of the 2011 Christchurch earthquake from seismological, geodetic, geological, and engineering perspectives.

Christchurch↗

Coulomb stress change sensitivity due to variability in mainshock source models and receiving fault parameters: A case study of the 2010-2011 Christchurch, New Zealand, earthquakes

Strong aftershocks following major earthquakes present significant challenges for infrastructure recovery as well as for emergency rescue efforts. A tragic instance of this is the 22 February 2011 M w 6.3 Christchurch aftershock in New Zealand, which caused more than 100 deaths while the 2010 M w 7.1 Canterbury mainshock did not cause a single fatality (Figure 1). Therefore, substantial efforts have been directed toward understanding the generation mechanisms of aftershocks as well as mitigating hazards due to aftershocks. Among these efforts are the prediction of strong aftershocks, earthquake early warning, and aftershock probability assessment. Zhang et al. (1999) reported a successful case of strong aftershock prediction with precursory data such as changes in seismicity pattern, variation of b-value, and geomagnetic anomalies. However, official reports of such successful predictions in geophysical journals are extremely rare, implying that deterministic prediction of potentially damaging aftershocks is not necessarily more scientifically feasible than prediction of mainshocks.

Christchurch↗

Ground motion attenuation during M 7.1 Darfield and M 6.2 Christchurch, New Zealand, earthquakes and performance of global predictive models

The M 7.1 Darfield earthquake occurred 40 km west of Christchurch (New Zealand) on 4 September 2010. Six months after, the city was struck again with an M 6.2 event on 22 February local time (21 February UTC). These events resulted in significant damage to infrastructure in the city and its suburbs. The purpose of this study is to evaluate the performance of global predictive models (GMPEs) using the strong motion data obtained from these two events to improve future seismic hazard assessment and building code provisions for the Canterbury region. The Canterbury region is located on the boundary between the Pacific and Australian plates; its surface expression is the active right lateral Alpine fault (Berryman et al. 1993). Beneath the North Island and the north South Island, the Pacific plate subducts obliquely under the Australian plate, while at the southwestern part of the South Island, a reverse process takes place. Although New Zealand has experienced several major earthquakes in the past as a result of its complex seismotectonic environment ( e.g. , M 7.1 1888 North Canterbury, M 7.0 1929 Arthur's Pass, and M 6.2 1995 Cass), there was no evidence of prior seismic activity in Christchurch and its surroundings before the September event. The Darfield and Christchurch earthquakes occurred along the previously unmapped Greendale fault in the Canterbury basin, which is covered by Quaternary alluvial deposits (Forsyth et al. 2008). In Figure 1, site conditions of the Canterbury epicentral area are depicted on a V S 30 map. This map was determined on the basis of topographic slope calculated from a 1-km grid using the method of Allen and Wald (2007). Also shown are the locations of strong motion stations. The Darfield event was generated as a result of a complex rupture mechanism; the recordings and geodetic data reveal that earthquake consists of three sub-events (Barnhart et al. 201 1 , page 815 of this issue). The first event was due to rupturing of a blind reverse fault with M 6.2, followed by a second event (M 6.9), releasing the largest portion of the energy on the right-lateral Greendale fault. The third sub-event (M 5.7) is due to a reverse fault with a right-lateral component (Holden et al. 2011). The Christchurch earthquake occurred on an oblique thrust fault. The comparison of spectral acceleration values at stations near Christchurch reveals that the second event produced much larger amplitudes of shaking than the Darfield event due to its proximity to the epicenter. Both events resulted in noticeably large amplitudes of the vertical motion, often exceeding horizontal motion in the near-fault area. The vertical motions, showing asymmetric acceleration traces and pulses, reached 1.26 g during the Darfield earthquake and 2.2 g during the Christchurch event. These events were recorded by more than 100 strong motion stations operated by the Institute of Geological and Nuclear Sciences (http://www.geonet.org.nz/). Using the processed data from these stations, peak ground acceleration (PGA) and 5%-damped spectral acceleration values at 0.3, 1, and 3 s are used for performance evaluation of the global ground motion predictive equations (GMPEs). The selected GMPEs are the Next Generation Attenuation (NGA) models of Abrahamson and Silva (2008), Boore and Atkinson (2008), Campbell and Bozorgnia (2008), and Chiou and Youngs (2008). The Graizer and Kalkan (2007, 2009) model, which is based on the NGA project database, is also included. These GMPEs are abbreviated respectively as AS08, BA08, CB08, CY08, and GK07. Because they have been used widely for seismic hazard analysis for crustal earthquakes, their performance assessment becomes a critical issue especially for immediate response and recovery planning after major events. The occurrence of aftershocks similar to the Christchurch event will most probably control seismic hazard in the broader area, as confirmed by the recent M 6.0 event on June 13, 2011.

Christchurch, Darfield↗

Seismologists meet “Down Under”

Planned presentations on hazard assessment and research weren't the only topics for discussion at the 27th General Assembly of the International Association of Seismology and Physics of the Earth's Interior (IASPEI), held from January 10–21 at the University of Victoria in Wellington, New Zealand. There was a more timely one unfolding: California's Northridge earthquake. When conference participants awoke to news of the earthquake on January 18, they scrambled to action to gather the latest data available. Just 12 hours after the event, they held a general discussion about the earthquake.

Wellington↗

Don't forget about the Christchurch earthquake: Lessons learned from this disaster

In the aftermath of the devastating magnitude-9.0 earthquake and tsunami that struck the Tohoku region of Japan on March 11, attention quickly turned away from a much smaller, but also highly destructive earthquake that struck the city of Christchurch, New Zealand, just a few weeks earlier, on Feb. 22. Both events are stark reminders of human vulnerability to natural disasters and provide a harsh reality check: Even technologically advanced countries with modern building codes are not immune from earthquake disasters. The Christchurch earthquake carried an additional message: Urban devastation can be triggered even by moderate-sized earthquakes.

Christchurch↗

Estimating seismic site response in Christchurch City (New Zealand) from dense low-cost aftershock arrays

The Mw 7.1 September 2010 Darfield earthquake, New Zealand, produced widespread damage and liquefaction ~40 km from the epicentre in Christchurch city. It was followed by the even more destructive Mw 6.2 February 2011 Christchurch aftershock directly beneath the city’s southern suburbs. Seismic data recorded during the two large events suggest that site effects contributed to the variations in ground motion observed throughout Christchurch city. We use densely-spaced aftershock recordings of the Darfield earthquake to investigate variations in local seismic site response within the Christchurch urban area. Following the Darfield main shock we deployed a temporary array of ~180 low-cost 14-bit MEMS accelerometers linked to the global Quake-Catcher Network (QCN). These instruments provided dense station coverage (spacing ~2 km) to complement existing New Zealand national network strong motion stations (GeoNet) within Christchurch city. Well-constrained standard spectral ratios were derived for GeoNet stations using a reference station on Miocene basalt rock in the south of the city. For noisier QCN stations, the method was adapted to find a maximum likelihood estimate of spectral ratio amplitude taking into account the variance of noise at the respective stations. Spectral ratios for QCN stations are similar to nearby GeoNet stations when the maximum likelihood method is used. Our study suggests dense low-cost accelerometer aftershock arrays can provide useful information on local-scale ground motion properties for use in microzonation. Preliminary results indicate higher amplifications north of the city centre and strong high-frequency amplification in the small, shallower basin of Heathcote Valley.

Christchurch↗

The Rakiura Titi Islands Restoration Project: Community action to eradicate Rattus rattus and Rattus exulans for ecological restoration and cultural wellbeing

In 2003, a non-profit group, Ka Mate Nga Kiore, was set up to oversee the restoration of four Maori-owned islands off the south coast of Stewart Island, New Zealand. The first step in the restoration was to eradicate ship rats ( Rattus rattus ) from three islands and Pacific rats ( R. exulans ) from another. The eradication was funded by the Command Oil Spill Trustee Council which managed the mitigation money from an oil spill off the Californian coast in 1998. The funding was coordinated via Oikonos Ecosystem Knowledge, a non-profit USA group primarily involved in seabird research and restoration. The project was primarily to benefit sooty shearwater ( Puffinus griseus ) and to sustain a culturally important customary harvest of their chicks by Rakiura Maori. However, like all island eradications, a wide range of other species also benefited from the removal of rats. The New Zealand Department of Conservation provided technical advice and assistance for the planning and implementation of the eradication programme. This paper describes how, with appropriate funding, community and technical support, rodent eradications can be achieved on private islands. In this case, a range of institutions and individuals joined to achieve a common goal that highlighted a significant international conservation action. We urge that more international and local-community-led restoration projects be initiated in the future.

Rakiura Titi Islands↗

Displaying seismic deaggregation: The importance of the various sources

Seismic hazard deaggregation has become a standard part of probabilistic seismic hazard assessment (PSHA). The first product of PSHA is calculation of the likely severity of ground motion at a given range of annual probability levels, and this is extremely important for seismic design of structures to be built at the site under examination. However, for full analysis of proposed structural designs, engineers also need to examine scenario events to produce detailed time histories. To select such scenarios, a deaggregation of the hazard is performed, whereby the details of sources that contribute to the annual frequency of exceeding specified levels of ground motion, or P exc , are identified. A common format for such a deaggregation is shown in Figure 1 . This relates to the 475-year peak ground acceleration (pga) at Wellington, New Zealand (41.28°S 174.77°E), and shows the distribution in magnitude and distance of sources that contribute to P exc . Return period is approximately the reciprocal of P exc . Stiff soil site conditions ( Standards New Zealand 2004 ) were assumed. The analysis in Figure 1 used the interim version of the updated seismic hazard model for New Zealand ( Stirling et al. 2007 ), with the attenuation function developed by McVerry et al. (2007). Based on a Poisson time dependence model, a return period of 475 years corresponds to a 10% probability of exceedance in 50 years. From Figure 1 , it is apparent that for this site the main contribution to ground motion of this severity is from earthquakes of magnitude about 7.6 less than 10 km from the site (blue), and there is another strong contribution from larger events in the distance range 10 to 20 km (red). These correspond to the Wellington and Wairarapa faults, respectively (see Table 1 ). There are other events less than 10 km from the site and small contributions from other sources. At this site the major contributions are from specific faults nearby, which are readily identified. At sites where there is significant background seismicity, however, the plot will be much more complicated and not so easy to interpret. Figure 1 deaggregates probabilistic pga at the site; other parameters are also commonly deaggregated in the same way, in particular response spectral acceleration at a variety of natural periods. But the figure has a major shortcoming in that it represents only one return period; to obtain a full appreciation of the various contributing sources it is necessary to perform a succession of analyses to cover the full range of return periods.

Seismological Research Letters↗

A multiple-approach radiometric age estimate for the Rotoiti and Earthquake Flat eruptions, New Zealand, with implications for the MIS 4/3 boundary

Pyroclastic fall deposits of the paired Rotoiti and Earthquake Flat eruptions from the Taupo Volcanic Zone (New Zealand) combine to form a widespread isochronous horizon over much of northern New Zealand and the southwest Pacific. This horizon is important for correlating climatic and environmental changes during the Last Glacial period, but has been the subject of numerous disparate age estimates between 35.1±2.8 and 71±6 ka (all errors are 1 s.d.), obtained by a variety of techniques. A potassium–argon (K–Ar) age of 64±4 ka was previously determined on bracketing lavas at Mayor Island volcano, offshore from the Taupo Volcanic Zone. We present a new, more-precise 40 Ar/ 39 Ar age determination on a lava flow on Mayor Island, that shortly post-dates the Rotoiti/Earthquake Flat fall deposits, of 58.5±1.1 ka. This value, coupled with existing ages from underlying lavas, yield a new estimate for the age of the combined eruptions of 61.0±1.4 ka, which is consistent with U–Th disequilibrium model-age data for zircons from the Rotoiti deposits. Direct 40 Ar/ 39 Ar age determinations of plagioclase and biotite from the Rotoiti and Earthquake Flat eruption products yield variable values between 49.6±2.8 and 125.3±10.0 ka, with the scatter attributed to low radiogenic Ar yields, and/or alteration, and/or inheritance of xenocrystic material with inherited Ar. Rotoiti/Earthquake Flat fall deposits occur in New Zealand in association with palynological indicators of mild climate, attributed to Marine Isotope Stage (MIS) 3 and thus used to suggest an age that is post-59 ka. The natures of the criteria used to define the MIS 4/3 boundary in the Northern and Southern hemispheres, however, imply that the new 61 ka age for the Rotoiti/Earthquake Flat eruption deposits will provide the inverse, namely, a more accurate isochronous marker for correlating diverse changes across the MIS 4/3 boundary in the southwest Pacific.

Quaternary Science Reviews↗

Quantifying production of salmon fry in an unscreened irrigation system: A case study on the Rangitata River, New Zealand

Diversion of out-migrant juvenile salmon into unscreened irrigation and hydroelectric canals is thought to have contributed significantly to declining populations of anadromous salmonids in the Pacific Northwest but is seldom studied in detail. Here we describe a program to study the fate of Chinook salmon Oncorhynchus tshawytscha fry diverted into the unscreened Rangitata Diversion Race (RDR) on the Rangitata River, New Zealand, by trapping fish in a random sample of on-farm canals in irrigation schemes (systems) served by the RDR. The catch rate at a site 9 km below the intake was strongly related to Rangitata River flow, but catches further downstream were unrelated to flow. Most fish entering the RDR were fry or early postfry (<70 mm fork length (FL)), but up to 92% of the fish trapped in on-farm canals were fingerlings (>70 mm FL), suggesting that many such fish became resident in the RDR for up to 3 months. Consequently, our estimate of the total number of fish leaving the RDR via on-farm canals (204,200 fish; 95% confidence limits = 127,100 and 326,700) is a conservative measure of the number lost from the Rangitata River because it does not allow for mortality within the RDR. We did not quantify the proportion of Rangitata River out-migrants that entered the RDR, but our results suggest that this figure was at least 5% and that it may have been as high as 25%, depending on mortality rates within the Rangitata River main stem and the RDR itself.

Rangitata River↗

Magma generation at a large, hyperactive silicic volcano (Taupo, New Zealand) revealed by U-Th and U-Pb systematics in zircons

Young (<∼65 ka) explosive silicic volcanism at Taupo volcano, New Zealand, has involved the development and evacuation of several crustal magmatic systems. Up to and including the 26·5 ka 530 km 3 Oruanui eruption, magmatic systems were contemporaneous but geographically separated. Subsequently they have been separated in time and have vented from geographically overlapping areas. Single-crystal (secondary ionization mass spectrometry) and multiple-crystal (thermal ionization mass spectrometry) zircon model-age data are presented from nine representative eruption deposits from ∼45 to ∼3·5 ka. Zircon yields vary by three orders of magnitude, correlating with the degrees of zircon saturation in the magmas, and influencing the spectra of model ages. Two adjacent magma systems active up to 26·5 ka show wholly contrasting model-age spectra. The smaller system shows a simple unimodal distribution. The larger system, using data from three eruptions, shows bimodal model-age spectra. An older ∼100 ka peak is interpreted to represent zircons (antecrysts) derived from older silicic mush or plutonic rocks, and a younger peak to represent zircons (phenocrysts) that grew in the magma body immediately prior to eruption. Post-26·5 ka magma batches show contrasting age spectra, consistent with a mixture of antecrysts, phenocrysts and, in two examples, xenocrysts from Quaternary plutonic and Mesozoic–Palaeozoic metasedimentary rocks. The model-age spectra, coupled with zircon-dissolution modelling, highlight contrasts between short-term silicic magma generation at Taupo, by bulk remobilization of crystal mush and assimilation of metasediment and/or silicic plutonic basement rocks, and the longer-term processes of fractionation from crustally contaminated mafic melts. Contrasts between adjacent or successive magma systems are attributed to differences in positions of the source and root zones within contrasting domains in the quartzo-feldspathic (<15 km deep) crust below the volcano.

Taupo volcano↗

Shear wave velocity variation across the Taupo Volcanic Zone, New Zealand, from receiver function inversion

The Taupo Volcanic Zone (TVZ), New Zealand is a region characterized by very high magma eruption rates and extremely high heat flow, which is manifest in high-temperature geothermal waters. The shear wave velocity structure across the region is inferred using non-linear inversion of receiver functions, which were derived from teleseismic earthquake data. Results from the non-linear inversion, and from forward synthetic modelling, indicate low S velocities at ∼6–16 km depth near the Rotorua and Reporoa calderas. We infer these low-velocity layers to represent the presence of high-level bodies of partial melt associated with the volcanism. Receiver functions at other stations are complicated by reverberations associated with near-surface sedimentary layers. The receiver function data also indicate that the Moho lies between 25 and 30 km, deeper than the 15 ± 2 km depth previously inferred for the crust–mantle boundary beneath the TVZ.

Taupo volcanic zone↗

Element content of Xanthoparmelia scabrosa growing on asphalt in urban and rural New Zealand

Xanthoparmelia scabrosa is a foliose lichen that grows abundantly on pedestrian and automobile asphalt in New Zealand, which are considered inhospitable habitats for lichens. Samples were collected at eight localities ranging from urban streets to very rural roads and analyzed for 28 chemical elements in order to determine elemental chemistry and to test hypotheses about tolerance mechanisms. Anthropogenic elements (Cu, Pb, and Zn) decreased significantly from urban to rural areas, while nutritional elements (K, P, and S) increased. Samples from urban areas contained 10% calcium. Sulfur was elevated at both urban and rural sites, possibly due to pollution in the former site and higher levels of sulfur-containing scabrosin esters at the rural sites. The ability of this lichen to accumulate high levels of Cu, Pb and Zn may make it useful as a remediation tool.

North Island↗

Estimating the sources and transport of nutrients in the Waikato River Basin, New Zealand

We calibrated SPARROW (Spatially Referenced Regression on Watershed Attributes) surface water‐quality models using measurements of total nitrogen and total phosphorus from 37 sites in the 13,900‐km 2 Waikato River Basin, the largest watershed on the North Island of New Zealand. This first application of SPARROW outside of the United States included watersheds representative of a wide range of natural and cultural conditions and water‐resources data that were well suited for calibrating and validating the models. We applied the spatially distributed model to a drainage network of nearly 5000 stream reaches and 75 lakes and reservoirs to empirically estimate the rates of nutrient delivery (and their levels of uncertainty) from point and diffuse sources to streams, lakes, and watershed outlets. The resulting models displayed relatively small errors; predictions of stream yield (kg ha −1 yr −1 ) were typically within 30% or less of the observed values at the monitoring sites. There was strong evidence of the accuracy of the model estimates of nutrient sources and the natural rates of nutrient attenuation in surface waters. Estimated loss rates for streams, lakes, and reservoirs agreed closely with experimental measurements and empirical models from New Zealand, North America, and Europe as well as with previous U.S. SPARROW models. The results indicate that the SPARROW modeling technique provides a reliable method for relating experimental data and observations from small catchments to the transport of nutrients in the surface waters of large river basins.

Waikato River↗

The enigma of the Arthur's Pass, New Zealand, earthquake: 1. Reconciling a variety of data for an unusual earthquake sequence

The 1994 Arthur's Pass earthquake ( M W 6.7) is the largest in a recent sequence of earthquakes in the central South Island, New Zealand. No surface rupture was observed, the aftershock distribution was complex, and routine methods of obtaining the faulting orientation of this earthquake proved contradictory. We use a range of data and techniques to obtain our preferred solution, which has a centroid depth of 5 km, M 0 =1.3 × 10 19 N m, and a strike, dip, and rake of 221°, 47°, 112°, respectively. Discrepancies between this solution and the Harvard centroid moment tensor, together with the Global Positioning System (GPS) observations and unusual aftershock distribution, suggest that the rupture may not have occurred on a planar fault. A second, strike slip, subevent on a more northerly striking plane is suggested by these data but neither the body wave modeling nor regional broadband recordings show any complexity or late subevents. We relocate the aftershocks using both one-dimensional and three-dimensional velocity inversions. The depth range of the aftershocks (1–10 km) agrees well with the preferred mainshock centroid depth. The aftershocks near the hypocenter suggest a structure dipping toward the NW, which we interpret to be the mainshock fault plane. This structure and the Harper fault, ∼15 km to the south, appear to have acted as boundaries to the extensive aftershock zone trending NNW-SSE. Most of the M L ≥ 5 aftershocks, including the two largest ( M L 6.1 and M L 5.7), clustered near the Harper fault and have strike slip mechanisms consistent with motion on this fault and its conjugates. Forward modeling of the GPS data suggests that a reverse slip mainshock, combined with strike slip aftershock faulting in the south, is able to match the observed displacements. The occurrence of this earthquake sequence implies that the level of seismic hazard in the central South Island is greater than previous estimates.

Arthur's Pass↗