Geology topics

Research about New Zealand

Source-linked reports with geographic coverage including New Zealand.

At least 19 recordsLinked to original sources

USGS2026

Determining Volcanic Risk in Auckland (DEVORA) Research Programme—A transdisciplinary approach to address the challenge of distributed volcanism in an urban environment

The Determining Volcanic Risk in Auckland (DEVORA) Research Programme was launched in 2008 to address the challenges associated with monogenetic volcanism in an urban setting and to enhance volcanic risk management in Tāmaki Makaurau Auckland in Aotearoa New Zealand. It is a multi-agency, increasingly transdisciplinary (defined here as research that transcends traditional disciplinary boundaries by integrating diverse types of knowledge, perspectives, and methods from academic and non-academic participants to create novel solutions to complex problems), and collaborative research program jointly led by Waipapa Taumata Rau University of Auckland and Earth Sciences New Zealand (ESNZ; formerly GNS Science), with core funding from Natural Hazards Commission Toka Tū Ake (NHC; formerly the Earthquake Commission, EQC) and Te Kaunihera o Tāmaki Makaurau Auckland Council (AC). The primary research focus of DEVORA is to investigate the geologic history, volcanic hazards, and risk posed by the basaltic intraplate Auckland Volcanic Field. Disruption from ash fall and gas from other Aotearoa New Zealand volcanoes is also considered. DEVORA’s work to explore exposure and vulnerability in Tāmaki Makaurau Auckland is also useful for assessing risks from other non-volcanic natural hazards, such as seismic and tsunami hazards. The greater Tāmaki Makaurau Auckland region has an ethnically and socio-economically diverse population of approximately 1.7 million, representing about one-third of the Aotearoa New Zealand population, and hosts critical infrastructure of national significance. The size and nature of the populace, consequential economic base, and important infrastructure within Tāmaki Makaurau Auckland mean that the effects of a volcanic eruption would be felt nationally, including through the disruption of air travel to Aotearoa New Zealand. The hazards from such an eruption could potentially affect hundreds of thousands of people, businesses, and lifelines (critical infrastructure). A considerable challenge for emergency and risk managers is the monogenetic nature of the volcanic field. It is not known where or when the next eruption will occur, how much warning we may get before an eruption, nor how an eruption and its effects might unfold. In this contribution, we highlight the concept and collaborative intent of the DEVORA Programme and show how it has evolved over the 16 years since its inception. We describe how DEVORA has unified more than 100 researchers (including more than 50 graduate students) and numerous stakeholders to address key issues facing Tāmaki Makaurau Auckland and describe how research findings are being implemented into policy and communicated to stakeholder agencies and the public. We also illustrate the broader influence of the DEVORA Programme and provide some learnings that might benefit others embarking on similar integrated projects, especially those focused on distributed volcanism in and near populated areas.

Auckland↗

USGS2026

Detecting snow avalanche activity using infrasound: Hooker Valley, New Zealand

Snow avalanches pose considerable hazards to people and infrastructure in alpine environments. Traditional avalanche monitoring relies on meteorological data and visual observations, which can be limited in scope and timeliness. Infrasound offers a promising complementary monitoring tool by detecting the low-frequency sound waves generated by avalanches. Here, we present infrasound and camera observations during a 50-day field campaign in the Hooker Valley of Aoraki/Mount Cook National Park, New Zealand. Our study detected seven avalanches with the cameras, whereas the infrasound system identified only one of these events, which was the largest and occurred under conditions that likely favoured infrasound propagation. The infrasound system recorded numerous other events not captured by the cameras, indicating the benefit of further investigation to determine their sources. These findings highlight the potential of infrasound technology for detecting avalanches and providing broad spatial coverage, capturing events in areas not monitored by cameras, while also showcasing limitations in infrasound capabilities. The limited detection of smaller avalanches underscores the opportunity for further research to enhance detection capabilities and understand environmental influences such as snow cover and wind noise. Overall, this study emphasises the utility of multidisciplinary monitoring techniques to improve avalanche detection in alpine environments.

Hooker Valley↗

USGS2026

From collision to subduction: Thermal-kinematic inversions constrain plate boundary structure and dip-slip activity in southwestern New Zealand

Constraints on three-dimensional fault system geometry and kinematics are required to understand both long-term tectonics and contemporary seismic hazard. In many cases, these parameters must be inferred from a combination of different data or analyses. Here, we invert low-temperature thermochronological data to quantitatively parameterize both the crustal-scale geometry and late Cenozoic dip-slip activity of the plate boundary system in southwest New Zealand. We find that the modern transpressive plate boundary system has accommodated the bulk of Australian-Pacific Plate convergence since establishment in the middle Miocene-early Pliocene. Shortening has occurred via (a) dextral-reverse slip on the moderate angle, listric central Alpine Fault décollement, and (b) reverse-slip on either the sharply curved, anti-listric northern Puysegur subduction interface or low-angle offshore thrust fault system that splays from the high-angle southern Alpine Fault at depth. We resolve late Cenozoic shortening rates on these structures that are indistinguishable, and interpret these to reflect the transfer of convergent deformation outboard of the predominantly strike-slip southern Alpine Fault. Our work helps elucidate the transition from strain localized oblique continental collision to strain partitioned subduction. We also provide an estimate of the across-strike width of the central-southern Alpine Fault section boundary at seismogenic depths. This is the first quantitative constraint on the geometry of a structural feature inferred to have arrested several major prehistoric earthquakes. This study furthers understanding of fault structure and dip-slip kinematics along the greatest sources of seismic hazard in and offshore of southwest New Zealand.

South Island↗

USGS2026

Prolonged cooling of andesitic-dacitic lava flows produces optimal groundmass material for 40Ar/39Ar dating

Organizing the eruptive histories and petrogenetic models that underpin hazard assessments at active volcanoes typically requires high-precision ages for late Pleistocene to Holocene lava flows. In this study, features of an existing 40Ar/39Ar dataset for andesitic-dacitic rocks from Ruapehu volcano, Aotearoa New Zealand, were interrogated with new microanalytical data to identify the optimal groundmass material for age determination. To understand the behavior of K and Ar within the cooling melt of lavas during their emplacement, we examined and compared the petrological features of samples from the ice-chilled margin and interior of a lava flow. The groundmass of the rapidly cooled sample comprises non-vesicular rhyolitic glass and microlites of plagioclase, orthopyroxene, and magnetite. The sample from the flow interior has minimal amounts of groundmass glass, and interstices are instead occupied by sanidine and tridymite. Modelling supports the interpretation that sanidine crystallized as a stable groundmass phase during slow cooling of degassed melt in the interior zones of lavas that were insulated behind their glassy margins. Our evaluation of the age data shows that sanidine contributed high yields of radiogenic argon during step-heating experiments for samples collected from lava interiors, which enabled high-precision 40Ar/39Ar ages to be determined. Using Ruapehu as an example of an arc volcano affected by lava-ice interaction, we provide a basic field guide for identifying lava flows that contain groundmass sanidine at the expense of glass. Target outcrops also offer opportunities for reconstructing geomagnetic dynamics and volcano–climate relationships during the Pleistocene and Holocene.

Ruapehu volcano↗

USGS2026

The contribution of a surge event to infilling in a barrier-enclosed estuary: Insights from field observations

Many estuaries worldwide face increasing sediment loading caused by catchment land use change and intensification, creating subsequent adverse effects on estuarine ecosystems. Extreme weather events can disproportionately alter sediment pathways and loading. Although storm-driven sediment exchange has been widely examined at open coasts and inlets, key transport mechanisms within constricted, sheltered estuaries remain understudied. This study presents an observational dataset capturing the impact of a 99th percentile water-level event (based on 20 years of records) on sediment transport pathways in a sheltered, barrier-enclosed estuary. This event, driven by a 3-day storm surge (>0.5 m) combined with a spring tide, was recorded during a 3-week field campaign. Sediment transport pathways and riverine contributions were analysed, and observations revealed substantial changes in suspended sediment concentrations increasing from 18 mg/l to 70 mg/l during the event. The elevated water levels and resulting pressure gradient at the constricted study site entrance caused by the storm surge increased local flood dominance. Combined with higher flow velocities and resuspension, the storm led to a sixfold increase in sediment import at the estuary entrance and a 600-fold increase in sediment flux to the upper estuary. The decoupling of peak suspended sediment concentrations from streamflow indicates that the resuspension of estuarine legacy sediment, rather than catchment inputs, dominated the system's response. These findings challenge assumptions about estuarine sediment budgets and emphasise that incorporating high water-level surge events into models can enhance the prediction of long-term estuarine evolution. Given projected increases in storm frequency under climate change, understanding these episodic but highly consequential sediment pulses can support the assessment of wetland resilience and inform estuarine management strategies.

Whangateau Estuary↗

USGS2025

Simulating demography, monitoring, and management decisions to evaluate adaptive management strategies for endangered species

Adaptive management (AM) remains underused in conservation, partly because optimization-based approaches require real-world problems to be substantially simplified. We present an approach to AM based in management strategy evaluation, a method used largely in fisheries. Managers define objectives and nominate alternative adaptive strategies, whose future performance is simulated by integrating ecological, learning and decision processes. We applied this approach to conservation of hihi (Notiomystis cincta) across Aotearoa-New Zealand. For multiple extant and prospective hihi populations, we jointly simulated demographic trends, monitoring, estimation, and decisions including translocations and supplementary feeding. Results confirmed that food supplementation assisted recovery, but was more intensive and expensive. Over 20 years, actively pursuing learning, e.g., by removing food from populations, provided little benefit. Recovery group members supported continuing current management or increasing priority on existing populations before reintroducing new populations. Our method can complement formal optimization-based approaches and improve AM uptake, particularly for programs involving many complex and coordinated decisions.

Conservation Letters↗

USGS2025

The communication of volcano information in New Zealand - A narrative review

Communication of volcano information is critical for effective volcanic risk management. A variety of information is communicated to inform decisions and guide actions for planning, preparedness, and response. Such information needs to be reliable, and fit-for-purpose across different stages of volcanic activity (quiescence, unrest, short or long-term eruptive stages, and the post-eruptive stage). However, an understanding of communication across these different stages of volcanic activity remains limited. We undertook a narrative review of New Zealand literature to explore what information is communicated about volcanoes, across which stages of activity and by whom. Results highlight that NZ literature only documents certain aspects of volcano information and communication, specifically regarding certain locations, stages of volcanic activity (i.e. quiescence or unrest), or hazards. Literature gaps exist regarding volcano communication during unrest and post-eruptive stages, as well as how volcano information evolves between these phases, and how decision-makers use such information. Additional work would be useful to document existing examples of volcano information for different stages of activity. Further research could help in understanding the information needs of decision-makers during each of these stages to improve information and communication.

New Zealand Journal of Geology and Geophysics↗

USGS2025

High-precision 40Ar/39Ar dating of Quaternary basalts from Auckland Volcanic Field, New Zealand, with implications for eruption rates and paleomagnetic correlations

The Auckland Volcanic Field (AVF), which last erupted ca. 550 years ago, is a late Quaternary monogenetic basaltic volcanic field (ca. 500 km 2 ) in the northern North Island of New Zealand. Prior to this study only 12 out of the 53 identified eruptive centres of the AVF had been reliably dated. Careful sample preparation and 40 Ar/ 39 Ar analysis has increased the number of well-dated centres in the AVF to 35. The high precision of the results is attributed to selection of fresh, non-vesicular, non-glassy samples from lava flow interiors. Sample selection was coupled with separation techniques that targeted only the groundmass of samples with < 5% glass and with groundmass feldspars > 10 μm wide, coupled with ten-increment furnace step-heating of large quantities (up to 200 mg) of material. The overall AVF age data indicate an onset at 193.2 ± 2.8 ka, an apparent six-eruption flare-up from 30 to 34 ka, and a ≤ 10 kyr hiatus between the latest and second-to-latest eruptions. Such non-uniformity shows that averaging the number of eruptions over the life-span of the AVF to yield a mean eruption rate is overly simplistic. Together with large variations in eruption volumes, and the large sizes and unusual chemistry within the latest eruptions (Rangitoto 1 and Rangitoto 2), our results illuminate a complex episodic eruption history. In particular, the rate of volcanism in AVF has increased since 60 ka, suggesting that the field is still in its infancy. Multiple centres with unusual paleomagnetic inclination and declination orientations are confirmed to fit into a number of geomagnetic excursions, with five identified in the Mono Lake, two within the Laschamp, one within the post-Blake or Blake, and two possibly within the Hilina Pali.

Auckland Volcanic Field↗

USGS2025

The crystalline silica respiratory hazard from rhyolitic lava dome eruptions in New Zealand's Taupo Volcanic Zone: A case study from the 1315 CE Kaharoa eruption

The rhyolitic Kaharoa 1315 CE eruption was a complex, long-lived event from Tarawera volcano, New Zealand. Explosive phases were followed by around 5 years of lava dome extrusion and collapse which produced block-and-ash flows (BAF). Lava domes generate crystalline silica in the form of cristobalite, and rhyolitic magmas often contain quartz phenocrysts. Fine-grained ash containing crystalline silica can be formed through dome collapses or explosions, generating a respiratory health hazard for communities affected by ashfall. The aims of this study are to: i) determine whether the Kaharoa eruption dome-forming phase generated substantial quantities of crystalline silica and, therefore, to determine the potential for future dome-forming eruptions of Tarawera to do the same; ii) consider the potential hazard of the crystalline silica by studying the crystal habit and chemistry compared to other lava domes, globally; and iii) assess the particle size and crystalline silica content of the Kaharoa ash, to inform a respiratory hazard assessment. Five co-BAF ash samples and one co-ignimbrite (explosive) ash sample from the Kaharoa pyroclastic deposits were analysed for health-pertinent factors: particle size distribution and crystalline silica content. Eight dome-rock samples were collected from the dome complex and associated BAF deposits and groundmass texture (especially forms of crystalline silica) and quantity of crystalline silica were assessed. Cristobalite was present in the 4 ash samples analysed by X-ray diffraction (XRD; 1.3–3.7 wt%) as was quartz (5.7–12.5 wt%). For the 4 dome samples analysed by XRD, all samples contained quartz (4.1–10.4 wt%) and two contained significant quantities of cristobalite (24.7 and 27.3 wt%). Of the two dome samples with minimal cristobalite (visible as individual vapour-phase crystals by SEM but not quantifiable by XRD), one was from the non-devitrified dome carapace and the other was from the compacted interior but had not undergone devitrification. The two dome samples with substantial cristobalite were from dome interiors and were highly devitrified, with well-developed spherulitic textures. Using energy-dispersive X-ray spectroscopy, cristobalite in all samples contained minor aluminium, as has been seen for volcanic cristobalite from other lava domes, which may ameliorate its toxicity. By laser diffraction, the quantities of ash in the health pertinent size fractions varied, with a range of 1.3–8.1 vol% for particles of <4 μm diameter and 1.7–15.6 vol% for particles of <10 μm diameter, which is lower than measured in ash from large-scale dome collapse events at other volcanoes. The findings suggest a potential for substantial crystalline silica to be formed in future Kaharoa-style eruptions, but that cristobalite generation is site-specific, depending on location within the dome and whether the dome remains sufficiently hot for spherulite formation and glass devitrification. Respiratory hazard will therefore vary depending on the collapse of (or explosions through) individual lobes – although all lava is expected to contain quartz phenocrysts – as well as the size and energy of those collapses, which will influence particle size and quantity of ash generated and dispersed.

Taupo Volcanic Zone↗

USGS2024

Shallow storage of the explosive Earthquake Flat Pyroclastics magma body, Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand: Evidence from phase-equilibria experiments

Rhyolitic tuffs range widely in their crystal contents from nearly aphyric to crystal-rich, and their crystal cargoes inform concepts of upper crustal magma reservoirs. The Earthquake Flat pyroclastics (Okataina Volcanic Center, Taupo Volcanic Zone, New Zealand) are 10 km 3 of rhyolitic tuffs with abundant (~ 40 vol.%) plagioclase and quartz, minor biotite, hornblende, and orthopyroxene, and accessory Fe-Ti oxides, apatite, and zircon, set in high-silica rhyolitic glass. Major minerals form large, euhedral phenocrysts and abundant glomerocrysts with few disequilibrium textures excepting some faintly resorbed quartz. Plagioclase phenocrysts have thick rims of nearly constant composition near An 30 , and hornblende is weakly zoned or unzoned. The abundant and texturally complex mineral assemblage contrasts with the nearby (~ 25 km), nearly synchronous, but more voluminous and crystal-moderate rhyolite tuffs from Rotoiti caldera. New H 2 O-saturated phase-equilibria results on the erupted Earthquake Flat melt (glass) determine its co-saturation with the partial phenocryst assemblage of plagioclase, quartz, biotite, and Fe-Ti oxides at: 140 MPa, 755 ºC. These closely approximate the conditions of the pre-eruptive magma body assuming it was saturated with nearly pure H 2 O and at an f O 2 of ~ Ni–NiO. Absence of hornblende and orthopyroxene from the synthesized assemblages may result from those minerals being in a peritectic reaction relation with melt to produce biotite, so they would not grow from the liquid used as starting material. Experimental results on Rotoiti rhyolite (Nicholls et al. 1992) show that the two bodies resided at similar pressures, temperatures, and f O 2 s. Lower crystal abundance of the Rotoiti tuffs may result from slight compositional differences. We interpret that the Earthquake Flat pyroclastics were sourced from the crystal-rich periphery of a mushy reservoir system with the Rotoiti occupying a more melt-rich central location. Uncertain is whether this was a single intrusion zoned continuously in crystallinity, or discrete adjacent intrusions, but our results illustrate and quantify complexities of magma storage across relatively short distances.

North Island, Taupo Volcanic Zone↗

USGS2024

Time varying crustal anisotropy at Whakaari/White Island volcano

Whakaari/White Island has been the most active New Zealand volcano in the 21st century, producing small phreatic and phreatomagmatic eruptions, which are hard to predict. The most recent eruption occurred in 2019, tragically claiming the lives of 22 individuals and causing numerous injuries. We employed shear-wave splitting analyses to investigate variations in anisotropy during quiescence, unrest, and an eruption between 2018 and 2020. We examined spatial and temporal variations in 7972 shear-wave splitting measurements and derived 1402 V p /IV s ratio measurements. Comparing shear-wave splitting parameters from similar earthquake paths across different times indicates that the observed temporal changes are unlikely to result from variations in earthquake paths through media with spatial variability. Instead, these changes may stem from variations in anisotropy over time, likely caused by changes in crack alignment due to stress or varying fluid content.

Whakaari/White Island volcano↗

USGS2024

Long-term communication of aftershock forecasts: The Canterbury earthquake sequence in New Zealand

On 14 February 2016, a magnitude (M)5.7 earthquake struck in Christchurch New Zealand (Aotearoa in the Maori language). The shaking caused damage to historic facades, power outages, cliff collapses, rock falls, and liquefaction but no reported injuries or fatalities. This Valentine’s Day earthquake was an aftershock in the Canterbury earthquake sequence (CES), which began on 4 September 2010 with the M7.1 Darfield Earthquake and included the destructive and fatal M6.2 Christchurch aftershock on 22 February 2011. This study, eight months after the Valentine’s Day earthquake and six years after the initiation of the CES, is the first to explore long-term aftershock forecast information and communication needs. The exploratory study also aimed to gather feedback on aftershock scenarios, an alternative form for communicating the forecast. The qualitative study involved workshops with emergency managers, public health officials, and members of the public in Christchurch. Key findings for long-term communication throughout an earthquake sequence include: 1. divergent earthquake experiences affect aftershock communication response and information needs; 2. understanding aftershock sequence behavior is foundational to sense-making when large aftershocks occur; 3. strategic earthquake sequence updates from the trusted science agency and local agencies could serve as important reminders for earthquake preparedness; 4. communication of aftershock forecast uncertainty could aid with both the credibility of the information and living with uncertainty, and 5. inclusion of impact information and preparedness advice into aftershock forecast scenarios could provide links to actionable information. The paper derives implications for research and practice of long-term communications during an aftershock sequence.

Christchurch↗

USGS2024

A decade of shaking in the Garden City: The dynamics of preparedness, perceptions, and beliefs in Canterbury, New Zealand, and implications for earthquake information

This study explored earthquake preparedness over time - before, during, and 10 years after the Canterbury Earthquake Sequence (CES) in Aotearoa New Zealand (NZ). Surveys of Canterbury residents were conducted in 2009, 2013 and again in 2021, using variables derived from Community Engagement Theory (CET). The surveys measured earthquake perceptions and beliefs, participation and engagement, and preparedness actions. Results were compared across the three samples. Findings indicate that perceptions and beliefs (e.g. risk perception, outcome expectancy beliefs), and types of preparedness actions taken (e.g. collection of survival items, structural preparedness, community and agency relationships), differed over time, depending on people’s experiences before, during, and after the CES. For example, during and after the CES people were more likely to believe that preparing provided a benefit to daily life, but less likely to think it could reduce property damage, perhaps due to people’s experiences of disruption and damage during the earthquakes. An understanding of such dynamics can assist with the provision and timing of risk and preparedness information. This study highlights the importance of providing applicable and actionable preparedness information, that is relevant to people’s experiences, throughout an earthquake sequence. Such information might evolve and change in focus over time depending on risks and needs. Focus could also be given to information that builds peoples beliefs and capacities to undertake preparedness in evolving situations. Understanding preparedness in the context of different experiences and timeframes is useful in helping update models such as the CET, where the dynamics of time might be better incorporated.

Canterbury↗

USGS2023

Cost-benefit analysis for evacuation decision-support: Challenges and possible solutions for applications in areas of distributed volcanism

During a volcanic crisis, evacuation is the most effective mitigation measure to preserve life. However, the decision to call an evacuation is typically complex and challenging, in part due to uncertainties related to the behaviour of the volcano. Cost-benefit analysis (CBA) can support decision-makers: this approach compares the cost of evacuating versus the expected loss from not evacuating, expressed as a ‘break-even’ probability of fatality. Here we combine CBA with a Bayesian Event Tree for Short-term Volcanic Hazard (BET_VHst) to create an evacuation decision-support tool to identify locations that are cost-beneficial to evacuate in the event of volcanic unrest within a distributed volcanic field. We test this approach with the monogenetic Auckland Volcanic Field (AVF), situated beneath the city of Auckland, New Zealand. We develop a BET_VHst for the AVF, extending a recently revised Bayesian Event Tree for Eruption Forecasting (BET_EF) to consider the eruptive style, phenomena produced, and the impact exceedance probability as a function of distance. The output of the BET_VHst is a probability of volcanic hazard impact at a given location. Furthermore, we propose amending the weight of the monitoring component within the BET_VHst framework to a transitional parameter, addressing limitations identified in a previous study. We examine how three possible transitional monitoring component weights affect the spatial vent likelihood and subsequent BET_VHst outputs, compared to the current default weight. For the CBA, we investigate four thresholds, based on two evacuation durations and two different estimates for the value of life that determine the cost of not evacuating. The combinations of CBA and BET_VHst are tested using a synthetic unrest dataset to define an evacuation area for each day. While suitable evacuation areas were identified, there are further considerations required before such an approach can be applied operationally to support crisis management.

Auckland Volcanic Field vent↗

USGS2023

Substantial upper plate faulting above a shallow subduction megathrust earthquake: Mechanics and implications of the surface faulting during the 2016 Kaikoura, New Zealand, earthquake

The 2016 moment magnitude 7.8 Kaikoura, New Zealand, earthquake occurred at the southern end of the Hikurangi subduction zone where the upper plate above the shallow megathrust is exposed sub-aerially. As a result, the substantial co-seismic deformation in the upper plate above the megathrust rupture was observed geologically and geodetically. We explore the relationship between this surface faulting and the subduction megathrust rupture and find that the greatest upper plate fault slip occurred coincident (in time and location) with the megathrust rupture. Models of Coulomb stress change demonstrate that these surface faults become positively loaded as the upper plate rebounds during the megathrust event, favoring fault slip. In addition, during the megathrust rupture these faults terminate against an uncoupled subduction plate interface. We simulate the effects of decoupling at the base of these faults and find that very large fault slip is an expected consequence of this decoupling, allowing near-complete strain release. In contrast, typical strike-slip faults, pinned at their base, would have lower amounts of fault slip. These two conditions—increased Coulomb stress and basal decoupling—combine to produce the extreme co-seismic upper plate faulting observed above the shallow Kaikoura megathrust earthquake. Similar conditions occur in other global subduction zones, but in most subduction zones the region above the coupled megathrust is underwater and poorly observed. Our analysis of the Kaikoura earthquake indicates a need to reevaluate patterns of strain accumulation and release in these regions, rather than assuming simple models of elastic rebound.

Tectonics↗

USGS2023

Identification of seismo-volcanic regimes at Whakaari/White Island (New Zealand) via systematic tuning of an unsupervised classifier

We present an algorithm based on Self-Organizing Maps (SOM) and k-means clustering to recognize patterns in a continuous 12.5-year tremor time series recorded at Whakaari/White Island volcano, New Zealand (hereafter referred to as Whakaari). The approach is extendable to a variety of volcanic settings through systematic tuning of the classifier. Hyperparameters are evaluated by statistical means, yielding a combination of “ideal” SOM parameters for the given data set. Extending from this, we applied a Kernel Density Estimation approach to automatically detect changes within the observed seismicity. We categorize the Whakaari seismic time series into regimes representing distinct volcano-seismic states during recent unrest episodes at Whakaari (2012/2013, 2016, and 2019). There is a clear separation in classification results between background regimes and those representing elevated levels of unrest. Onset of unrest is detected by the classifier 6 weeks before the August 2012 eruption, and ca. 3.5 months before the December 2019 eruption, respectively. Regime changes are corroborated by changes in commonly monitored tremor proxies as well as with reported volcanic activity. The regimes are hypothesized to represent diverse mechanisms including: system pressurization and depressurization, degassing, and elevated surface activity. Labeling these regimes improves visualization of the 2012/2013 and 2019 unrest and eruptive episodes. The pre-eruptive 2016 unrest showed a contrasting shape and nature of seismic regimes, suggesting differing onset and driving processes. The 2016 episode is proposed to result from rapid destabilization of the shallow hydrothermal system, while rising magmatic gases from new injections of magma better explain the 2012/2013 and 2019 episodes.

Whakaari, White Island↗

USGS2022

Towards real-time probabilistic ash deposition forecasting for New Zealand

Volcanic ashfall forecasts are highly dependent on eruption source parameters (ESPs) and synoptic weather conditions at the time and location of the eruption. In New Zealand, MetService and GNS Science have been jointly developing an ashfall forecast system that incorporates four-dimensional high-resolution numerical weather prediction (NWP) and ESPs into the HYSPLIT model, a state-of-the art hybrid Eulerian and Lagrangian dispersion model widely used for volcanic ash. However, these forecasts are based on discrete ESPs combined with a deterministic weather forecast and thus provide no information on output uncertainty. This shortcoming hinders stakeholder decision making, particularly near the geographical margin of forecasted ashfall and in areas with large gradients in forecasted ash deposition. Our study presents a new approach that incorporates uncertainty from both eruptive and meteorological inputs to deliver uncertainty in the model output. To this end, we developed probability density functions (PDFs) for three key ESPs (plume height, mass eruption rate, eruption duration) tailored to New Zealand’s volcanoes and combine them with NWP ensemble datasets to generate probabilistic ashfall forecasts using the HYSPLIT model. We show that the Latin Hypercube Sampling (LHS) technique can be used to representatively span this four-dimensional parameter space and allow us to add uncertainty quantification to rapid response forecast systems. For a case study of a hypothetical eruption at Tongariro, New Zealand we suggest that large parts of New Zealand’s North Island would not receive adequate warning for potential ashfall if uncertainties were not included in the forecasts. We also propose new probabilistic summary products to support public information and emergency responders decision making.

Journal of Applied Volcanology↗

USGS2022

Hybrid broadband ground-motion simulation validation of small magnitude active shallow crustal earthquakes in New Zealand

This article presents a comprehensive validation of the hybrid broadband ground-motion simulation approach (via the commonly used Graves and Pitarka method) in a New Zealand context with small magnitude point source ruptures using an extensive set of 5218 ground motions recorded at 212 sites from 479 active shallow crustal earthquakes across the country. Modifications to the simulation method inferred from a previous New Zealand validation are implemented, and the improvements are explicitly quantified. Empirical ground-motion models are also considered to provide a benchmark for simulation prediction accuracy and precision. Examination of intensity measure residuals identifies that the simulation method modifications lead to reduced model prediction bias and within-event variability and provides evidence toward the use of spatially varying coefficient models for simulation parameters, such as the high-frequency Brune stress parameter. Additional biases identified include, among others, underprediction of significant durations at soft soil sites and overprediction of short-period pseudo-spectral accelerations at stiff alluvial gravel and rock sites due to low-estimated 30 m time-averaged shear-wave velocity values.

Earthquake Spectra↗