Classification of long-term very long period (VLP) volcanic earthquakes at Whakaari/White Island volcano, New Zealand

We have observed very long period earthquakes (VLPs) over the period 2007 to the end of 2019 at Whakaari/White Island volcano, New Zealand. The earthquakes exhibit similitude between waveforms which suggests repeating source locations and processes. VLPs recorded at two permanent stations were detected using waveform semblance and were then classified into two main families (F1 and F2) using a clustering analysis. The two families are characterized by ‘mirror image’ reverse waveform polarity suggesting that they are genetically related, but occur during different evolutionary phases of volcanic activity. F1 events occurred throughout the observation period, while F2 events mainly occurred as swarms that mark the onset of volcanic unrest. A detailed cluster analysis reveals possible sub-families implying slight temporal evolutions within a family. Our results add to our understanding of the volcanic magma–hydrothermal system at Whakaari/White Island indicating that relatively stable VLP sources may be exploited to improve monitoring for future unrest.


Introduction
Whakaari/White Island (Fig. 1a) is one of New Zealand's most active volcanoes and has been a frequent tourist destination. Since 2007, this volcano has produced seven hazardous eruptions, and a recent eruption in December 2019 tragically resulted in 21 fatalities. Hence, there is a need for improved monitoring and research of Whakaari/ White Island dynamics (e.g., Chardot et al. 2015). Seismic analysis remains a highly valuable tool for volcano monitoring although recent failures to forecast hydrothermal or phreatic eruptions have highlighted the need for additional research worldwide (Yamaoka et al. 2016). Whakaari/White Island has produced various types of volcanic earthquakes with different frequency contents including very long period (VLP; f < 0.1 Hz), long period (LP; f = 0.5-2 Hz), high frequency (HF; f = 2-10 Hz) and volcano tectonic (f > 5 Hz) events (Nishi et al. 1996;Jolly et al. 2017). The long wavelengths of VLP are conducive to source characterization because their waveforms are less impacted by topography and subsurface structural heterogeneity. VLPs have been observed during eruptions at Whakaari/White Island (Caudron et al. 2018;Jolly et al. 2018) as well as many other volcanoes: e.g., Stromboli in 2007 (Giudicepietro et al. 2009), Redoubt in 2009 (Haney et al. 2013), and Merapi in 2010 (Jousset et al. 2013). In contrast, Zuccarello et al. (2013) observed VLP events during degassing without eruptive activity at Mount Etna.
VLP events are often studied using waveform similarity. At Erebus volcano, three VLP families associated with strombolian eruptions were classified based on initial polarity, and waveform and spectral similarity (Mah 2003). Esposito et al. (2008) clustered VLP waveforms according to different eruptive vents of Stromboli volcano. Dawson and Chouet (2014) found that three VLP types are linked to different surficial and subsurface activity at Kilauea volcano (~ 7 years). These studies revealed stationary sources from several days to years. Relatively long-term VLP activity is documented at Kilauea (Dawson and Chouet 2014) and Popocatepetl (Arciniega-Ceballos et al. 2003) volcanoes, and other volcanoes such as Erebus (Mah 2003;Aster et al. 2008), Etna (Saccorotti et al. 2007), and Stromboli Esposito et al. 2008) also revealed short-or long-term stable VLP activity. Most VLP studies evoke pressure variations related to gas and magma intrusion (e.g., Arciniega-Ceballos et al. 1999;Aster et al. 2003;Molina et al. 2008;Nakamichi et al. 2009;Dawson et al. 2010;Lyons and Waite 2011;Maeda and Takeo 2011). Previous VLP studies at Whakaari/White Island have mainly dealt with their source properties. On 19 to 21 August 2011, a swarm of mixed frequency events including VLP, LP, and HF components was recorded on a temporary seismic array . The authors reported a depth of 0.8-1.5 km for VLP, while contemporaneous LP and HF components occurred at shallower depths (< 1 km). They interpreted the coupled events as part of the migration of fluids into the shallower hydrothermal system. Caudron et al. (2018) studied two VLP events recorded during the 2013 eruption episode and found a similar source depth to the 2011 swarm . VLP signals were also observed before and during the 2016 eruption sequence, and their source locations of 0.8-1.0 km depth and mechanisms of volumetric changes were investigated by waveform semblance and source inversion, respectively . For all cases, VLP originated beneath the active vent at ~ 1 km depth and were inferred to result from fluid migration at the interface between the shallow magmatic conduit and the overlying hydrothermal system.
Although some studies have assessed long-term trends of volcanic activity including volcanic tremor (Chardot et al. 2015), chemical and lake data (Christenson et al. 2017), and seismic ambient noise (Yates et al. 2019), there is generally a lack of research on long-term observation of VLP earthquakes at Whakaari/White Island. As the assessment of long-term volcano hazards becomes increasingly important, it is clear that a rigorous assessment of VLP occurrence and their relation to volcanic activity is required. However, the detection and discrimination of these important event types can be challenging.
In this paper, we observe enduring VLP earthquake families recorded at Whakaari/White Island for the period of 2007-2019. The VLP events are distinguished by using both waveform semblance and a templatematching technique. For the resulting catalogue, we examine their temporal variations and compare them to heightened unrest and eruption occurring during the period 2011-2016 and 2019, in order to assess possible links between VLP activity and eruptions. We find that VLP source excitations are pervasive during periods of both unrest and repose, but that VLP families relate to different periods of quiescence and unrest/eruption. The results have implications for the assessment of unrest at Whakaari/White Island volcano. Fig. 1 a Location map of White Island volcano (triangle) including two permanent seismic stations (dots), WIZ and WSRZ, and the active vent (asterisk); b the station availability and methods used in this study. Note that 2sta and 1sta indicate the use of two stations and one-station analysis, respectively. c One hour sampled RSAM of vertical component data (μm/s) recorded at WIZ station from bandpass-filtered data (2-5 Hz). The onset times of seven eruptions and unrest phases are marked by reversed triangles (red) and shaded parts (green), respectively

Background at Whakaari/White Island
GeoNet has operated two permanent stations, WIZ and WSRZ, each consisting of three-component Guralp CMG-3ESP 60-s broadband sensors and 24-bit Quanterra Q330 digitizers which have operated since 20 April 2013; between 2007 and 2013 a single broadband sensor at WIZ was operated (Fig. 1a, b). Prior to this, the WIZ station comprised a short period sensor which was incapable of detecting VLP events . In this study, we use the continuous seismic data recorded from 19 April 2007 to 31 December 2019 to identify VLP signals. All seismic records have 100 Hz sampling frequency.
Whakaari/White Island experienced two unrest episodes, including eruptions, since broadband seismic monitoring was established. The onset of volcanic unrest is defined by integration analysis of a full range of geological and geophysical observations including seismicity, gas discharge and deformation; hence a study of past unrest may provide evidence of potential upcoming unrest (Potter et al. 2015). One of the crucial elements for identifying unrest at Whakaari/White Island is the realtime seismic amplitude measurement (RSAM; Fig. 1c). The RSAM is computed from the RMS amplitude of 1-h sampled velocity data in the frequency range 2-5 Hz and allows an assessment of volcanic tremor energy (Chardot et al. 2015). The RSAM increased from early 2011 (Christenson et al. 2017), and the onset of the first unrest period was marked in August 2011 at Whakaari/ White Island after integrated analysis with visual observation (Chardot et al. 2015;Edwards et al. 2017). The amplitude rapidly decreased to background levels in May 2016. During the first unrest period in 2011-2016, three eruptive sequences occurred in 2012, 2013 and 2016 containing six distinct eruptions. In early August 2012, the Volcanic Alert Level (VAL) was raised to 2 indicating moderate volcanic unrest due to an increase in tremor and degassing activity. Following the level increased, the volcano erupted ash and blocks, and minor less distinct mud eruptions occurred later in August (Chardot et al. 2015). During the period of January to July 2013, mud and dry ash venting frequently occurred depending on the condition of the crater lake Edwards et al. 2017). The second eruptive episode occurred in August and October 2013 with four short lived eruptions which included steam, mud and explosive activity (Chardot et al. 2015;Caudron et al. 2018). On 27 April 2016, a phreatic eruption was observed Kilgour et al. 2019;Walsh et al. 2019). From April 2019, the RSAM slightly increased, followed by an increasing trend in June. GeoNet raised the VAL to 2 in the middle of June and November 2019. Following the VAL change, the tremor levels increased further and culminated in an explosive eruption on 9 December 2019 at 14:11 NZDT (UTC 01:11).

Detection of VLP volcanic signals
To identify VLP signals associated with volcanic activity, we applied a three-step process depending on station availability (Fig. 1b). We (1) detected VLP volcanic signals based on waveform semblance using the two-station data recorded 2013-2019 and (2) determined template events using the 2013-2019 semblance triggered signals by a clustering analysis for VLP identification prior to 2013. We then (3) applied a template-matching method over the 2007-2013 period when only one station was recording, followed by the removal of spurious signals. For detailed explanations refer to the following sub-sections. All seismic recordings are bandpass filtered at 0.03-0.1 Hz to isolate the most energetic part of the recorded VLP signals to be consistent with previous studies (Caudron et al. 2018;Jolly et al. 2017Jolly et al. , 2018. This frequency band, however, can be affected by ocean microseism at approximately 0.07 Hz (Yates 2018).

Triggering VLPs using waveform semblance
A Whakaari/White Island VLP detection approach was outlined in Jolly et al. (2018) where several hours of preeruption VLPs were identified prior to the April 27, 2016 eruption by application of a moving window grid search for the waveform semblance (Kawakatsu et al. 2000). This approach exploits radial energy propagation among seismograms at the selected VLP passband. The detection algorithm computed the 3-min waveform semblance within a coarse 3D grid centred on the main vent area in which specific VLP activity may occur. The grid was 2 × 2 × 1.5 km having lateral and vertical spacing of 1 and 0.5 km, respectively. The grid geometry allows discrimination between VLP excitations originating at a range of depths beneath the volcanic vent (coherence from the centre grid elements) and non-volcanic excitations (coherence on margin elements of the grid which relate to teleseismic excitations beyond the volcano) ). This sparse array algorithm was initially tested against the 2013 and 2016 eruption sequences with various semblance thresholds beneath the vent area regardless of depth to assure that it could successfully capture the target VLP seismicity. After testing, we extended the algorithm to the period April 2013 to December 2019 using the semblance threshold of 0.65, when two permanent GeoNet stations were operated. Note that the semblance is stable at high signal-to-noise ratio (SNR > 3) in Whakaari/White Island dataset . We thus selected VLP events, which have high median SNR of 3-component bandpass-filtered records from both stations, giving 92 events ( Fig. 2a and Additional file 1).

Determining template events
The 92 selected events were used to calculate the correlation coefficient (r) between all pairs of 3-component data at both stations. We then completed a hierarchical clustering analysis based on dissimilarities or Pearson correlation distance defined as 1 − r (Fig. 2b). The distance ranges from 0 (r = 1) to 2 (r = − 1). For clarity, we use similarity given by r for the rest of the paper. The concept of hierarchical clustering and corresponding algorithms are well documented in Duda and Hart (1973) and Aminzadeh and Chatterjee (1984). In this study, the cluster tree or dendrogram was created using an average linkage algorithm. By considering correlation of multi components, we can determine families, potentially produced by similar sources/processes. In the dendrogram (Fig. 2b), two main clusters (F1 and F2) are determined at a similarity threshold of r = 0.6, and the remaining signals (U) belong to neither F1 nor F2 clusters possibly due to noise contamination or variations associated with a nonrepeating source. We also confirmed that U is irrelevant to the clusters by a clustering analysis based on the absolute value of correlation coefficient (1 − |r|; Additional file 2). Several sub-groups from the main clusters can be divided as the threshold increases. However, they are dissimilar to each other at only high correlation threshold, hence the two main families with lower threshold are Fig. 2 Triggered VLP signals based on two-station waveform semblance. a Peak-to-peak (P2P) amplitude (μm/s) of VLP families (F1 and F2) and unrelated events (U; gray) determined by b a cluster tree. The amplitude is computed using the vertical component of WIZ station. This dendrogram is created based on the correlation distance in all pairs of 92 VLP events and its x axis is event index regardless of a time order (can check VLP triggered time and index information in Additional file 1). c Stacked waveforms (thick coloured lines) of individual events (thin gray lines) for the two identified families recorded at WIZ and WSRZ stations. Frequency spectra are computed using each unfiltered example waveform recorded at WIZ (black) and WSRZ (green). All amplitudes are normalized. The stacking procedure is based on actual amplitude, and the stacked events are used as templates (see in "Template-matched VLPs" section) regarded as the most robust grouping. Template events were determined by stacking each family to minimize noise effect (Fig. 2c). VLP signals within each family display similar waveforms event-to-event (Fig. 2c) which may indicate a stable source process (e.g., Chouet 1996;Stephens and Chouet 2001;Green and Neuberg 2006;Petersen et al. 2006;Park et al. 2019), despite the variation in surface activity (Edwards et al. 2017). They also show temporal variations in amplitude implying a source intensity change.

Template-matched VLPs
The 2-station semblance catalogue only encompasses about half of the 2011-2016 unrest period, hence we sought to develop a method to robustly evaluate the earlier period when only one broadband sensor was available (Fig. 1b). In order to detect VLP events recorded prior to April 2013 at WIZ station, we applied a template-matching method which is based on waveform correlation (e.g., Chouet 1996;Stephens and Chouet 2001;Green and Neuberg 2006;Petersen et al. 2006;Park et al. 2019) between the two stacked waveforms (Fig. 2c) and continuous 3-component data. The optimal correlation threshold was found by applying the method to the seismic recording for the 2013-2019 period containing the triggered semblance events; 91% of these events were successfully detected by using the correlation threshold of 0.7. The template-matching catalogue for the period 2007 to 2013 contains VLP signals and irrelevant events which may be produced outside of the vent area including local earthquakes and teleseisms (Additional file 3). To exclude unrelated long period events, we applied a single-station semblance filter. Although this approach is inadequate to completely eliminate teleseisms, it contributes by capturing radiational transients from the vent direction. To select an optimum semblance threshold of one-station data, we compared both semblance values derived from twostation data and one-station data, and the consequent single-station VLP catalogue was limited by semblance of ≥ 0.70. We also restricted the catalogue by the aforementioned SNR condition (≥ 3) and filtered out a small number of extraneous transients from the single-station catalogue using both an empirical SNR filter and visual inspection (Fig. 3). A total of 76 VLP events are selected from 2007-2013 catalogue by template-matching and single-station semblance, and 92 events from 2013-2019 data by semblance of two stations. We then assess both the unrest periods and a long quiescent period using the combined Whakaari/White Island VLP catalogue which consists of 168 events (Fig. 3).

Fig. 3
A combined template-matching catalogue after application to one-station semblance filtering and SNR filters. The yellow filled circles indicate eruptive VLPs. The seven eruptions and swarm periods are denoted by reversed triangles (red) and horizontal bars (black) on the top and bottom of the figure, respectively. Volcanic unrest phases are marked by shaded parts (green), and white arrows indicate VLP events used in previous studies Park et al. Earth, Planets and Space (2020) 72:92 Combined VLP catalogue A total of 168 VLP pulses are extracted from 2007-2019 seismic data using the template-matching and semblance filters described above (Fig. 3). With the two-station semblance and clustering analysis, we classified triggers into two super-families (F1-F2) based on the waveform similarity. A super-family was defined as a group of similar signals implying analogous source process even with slight variations (Park et al. 2019). The F1 event type represents 64 VLPs, whereas F2 type was more active and included 99 events. The remaining events (U) are not defined as either of the families and may be impacted by noise or may relate to a non-repeating source process. F1 type occurred over the entire observation period. By contrast, about 74% of F2 family are distributed in swarm activity in 2010-2011 and 2017-2018 (Fig. 3). The super-families have possible sub-clusters depending on the selected similarity threshold. This can be interpreted in the context of the temporal evolution of volcanic activity. To examine variations in sub-families, we re-conducted the cluster analysis for each family using 3-component data (Fig. 4). Under a similarity threshold of 0.80, the F1 family is divided into seven slightly different family types, and their distribution yields temporal variations without meaningful features (Fig. 4a).
In contrast, the F2 family shows systematic variations between the primary swarms and non-swarm sub-clusters (Fig. 4b). Note that this suggests a more stable source process for swarm activity compared to non-swarm events over the observation period. If we choose a higher threshold of 0.85, one of F1 sub-families (orange circles in Additional file 4) was separated for the first unrest, but has a similar waveform to others which occurred for the repose period. As the correlation threshold increases, the number of clusters also increases, but the waveform characteristics become difficult to distinguish. Due to the potential for noise contamination, caution should be exercised with an interpretation of sub-family occurrence in the case where a higher similarity threshold is applied.
The two super-families show slight waveform evolution within a family group, and their peak-to-peak amplitude varied over the observed time. To check the evolutionary process, we also computed correlation coefficients between individual VLP events and each family stacked waveform as outlined by Park et al. (2019). The results indicate multimodal features implying that the two distinct families occur simultaneously without slow or rapid evolution between different families. Interestingly, a strong negative correlation of − 0.97 was observed between stacked waveforms of F1 and F2 events (Fig. 2c),   Fig. 4 Temporal variations in sub-families within a F1 and b F2 resulted from correlation threshold of r = 0.80 in the cluster tree. Note that a range of correlation distance is originally from 0 to 2, and y-scale is adjusted. All markers on the figure are the same in Fig. 3. The colours in both peak-to-peak amplitude and cluster tree correspond which suggests a 'mirror image' reversed waveform polarity for the two families. We also note that the occurrence of F1 events usually coincided with a dearth of F2 events before 2017 (Fig. 3).

Link between VLP seismicity and unrest
This analysis revealed two main VLP types during the period of 2007-2019. Regardless of event types, the occurrence of VLP earthquakes suggests ongoing activity in the volcanic edifice since they originate directly beneath the active vent and are interpreted to represent the response to fluid mass transport . To understand long-term VLP activity, we first attempted to compare our observations with previous studies and highlight particular VLP events associated with volcanic unrest at Whakaari/White Island.
We note that the occurrence of VLP swarm activity was observed before the commencement of both unrest periods (Fig. 3). At the beginning of the first unrest, VLPs occurred as the same event type of F2 family and similar amplitude distribution in June, August and December 2011. The August swarm was interpreted to be from gas slug ascent related to the subsequent 2012-2013 eruptions . Following the swarm activity, the occurrence of F2 events substantially decreased and F1 events were more prevalent. Although repeating F2 VLP activity decreased, the volcano first erupted on 4 August 2012 followed by a sequence of degassing and dome forming eruption along with the strong volcanic tremor (Chardot et al. 2015).
The first unrest period ended after the April 2016 eruption and weak tremor was observed for more than 1 year. Vigorous VLP activity initiated from late 2017 and consisted of F2 type with a small number of F1 events. Unlike the 2011 activity, higher amplitude F2 swarm again occurred instead of eruptions in the following year. However, a small number of F2 events were also detected in August 2010. We also assessed the occurrence of F2 events having lower amplitude in the template-matching results without the SNR and semblance filters (Additional file 3). Following the 2018 swarm F2 events diminished and the seismic activity evolved into elevated RSAM beginning in 2019. The RSAM unrest was then followed by the eruption on 9 December 2019.
Based on the swarm observations, we note that a VLP swarm composed primarily of F2 family events has preceded two unrest and eruption episodes at Whakaari/ White Island. This observation requires additional example unrest episodes to develop this possible evolutionary pathway at Whakaari/White Island. Additionally, we note that the occurrence of VLP swarms and the onset RSAM is indistinct and the 2017 swarm in particular, did not include an associated increase in RSAM.
In contrast, F1 type events were irregularly observed over the entire period. Although the sub-clusters in the similarity threshold of 0.85 show the difference between repose and unrest event types, it is uncertain whether this represents a source change under our data evaluation conditions. Interestingly, three high-amplitude VLP events from the F1 family (yellow filled circles in Fig. 3), contemporaneously accompany the phreatic eruptions on 3 October 2013 (Caudron et al. 2018) and 27 April 2016 . However, because high-amplitude F1 events were also observed during the repose period and the five other eruptions produced other types of volcanic earthquakes, this event type requires additional study and more eruption examples.
Waveforms within a family show modest changes over time, but the changes are too indistinct to resolve features that might relate to changes in source location or mechanism. Their amplitudes also varied over time which might imply change in source intensity. Consequently, we can assume that a family is generated by a similar source process and potential location. According to previous studies, both F1 and F2 families probably related to the upward migration of gas slugs Caudron et al. 2018). However, Jolly et al. (2018) noted strong static deformation during the 2016 eruption. They suggested that this deformation could have resulted from magmatic injection into the overlying magma carapace. Hence, large F1 events might include near-conduit magma injection. Gas-related VLP mechanisms are commonly invoked at volcanoes. For example, at Stromboli volcano, Ripepe et al. (2005) interpreted a close link between rates of gas and VLP rate during the 2003 eruption. Kazahaya et al. (2011) also examined that SO 2 flux was proportional to moment of isolated VLP pulses recorded for a day during gas bursts at Mount Asama. In such short-term analysis, VLP activity is likely to display a positive relationship with gas measurements. However, a complex interplay within a volcano conduit and its overlying crater lake, may lead to inconsistent phenomena leading to ambiguous interpretation of the long-term data. In a conceptual model of Whakaari/ White Island (Cole et al. 2000;Christenson et al. 2017;Jolly et al. 2017), the volcano can be described as a complex vent system in elastic altered rock that is brittlely fragmented by fumarolic pathways with more ductile rock where pressure and porosity allow (Heap et al. 2015). Due to this complicated system, a study on a possible link between surface observations, other monitoring and VLP activity requires additional assessment that is beyond the scope of this study.
The long-term stability of VLP families, at ~ 1 km depth, implies a stable and non-destructive source process regardless of eruptive activity. The aforementioned small-scale deformations for VLP are noted to produce modelled volume changes larger than 10 4 m 3 . The long-term stability of the source process implies that the subsurface plumbing system at the VLP source depth is not modified by explosive activity or the large volumes involved with individual events. Indeed, the observations may reflect the limited range of persistent fluid flux from the underlying magma system and into the shallow hydrothermal system that may be persistent and stable over decades.
It is additionally noted that the opposite polarity for the two primary sources might imply a change in the volumetric or slip direction within the excited source. This change is seen in both vertical and lateral seismometer components (Fig. 2c). Hence, the reverse polarity is robustly observed. The two event types in this study were previously modelled by Jolly et al. (2018;F1 family type) and Jolly et al. (2017;F2 family type). The prior modelling suggests that the mirror image source process relates to an excitation of the same oriented crack but having the opposite sense of slip and volumetric change through time. Future modelling of these events will focus on the nature of inflation/deflation (feed/release) modes to understand plausible source mechanisms.
Although slight waveform evolutions, amplitude variations and the number of events were recognized within a family which might be due to volcanic activity, no distinct new family types are seen over the observation period.
From the above sequence, we highlight three key observations: (1) VLPs show remarkably distinct stable repeating characteristics with little evolution over a 12-year study period; (2) swarms composed of F2 event type were observed before the commencement of both unrest periods, and (3) the 2013 and 2016 eruptions produced F1-type events having high amplitude. F1 events, however, were prevalent over the entire observation time, and an assessment of such features requires further analysis and more eruption cases.
Eruptive activity at Whakaari/White Island is usually observed with an increase in occurrence of volcanic tremor which can be detected by the increasing RSAM trend (e.g., Chardot et al. 2015). Unlike RSAM observation, most VLP events have no clear increase in event size before an eruption (Additional file 5). Whilst a rigorous statistical assessment cannot be completed given the infrequent nature of eruptive activity at Whakaari/ White Island, the observations of the two super-families might be useful for future monitoring assessments at this volcano.

Conclusions
This study provides the first attempt of the comprehensive detection and analysis of long-term VLP earthquakes at Whakaari/White Island. We have detected and discriminated long-term VLP earthquakes at the volcano using template-matching and waveform semblance methods. We note that the volcano produced two distinct types of VLP families which occur throughout the progression from repose to build-up and unrest including eruptions. The occurrence of these distinct event types is stable over the 12 years and implies a relatively stationary source and mechanism throughout this time. We emphasize that VLP swarms of type F2 appear to relate closely to the onset of volcanic unrest.
We propose that automatic detection and classification of these event types would contribute to improved volcano monitoring of unrest at Whakaari/White Island. In addition, we expect that the detection methods introduced in this study may be effective for sparse stations monitoring for VLP activity at volcanoes worldwide.
Although beyond the scope of this study, interesting swarm observation and the waveform polarity switch might indicate changing source locations and processes which relate to the evolution of subsurface conditions during periods of repose, build-up, unrest and eruptions. Further work may include wider data observations such as InSAR/GPS measurements, gas discharge and lake level changes. Future monitoring efforts should emphasize the collection of higher density observation of such data. These data together with new dense observations of the VLP source will improve our understanding of Whakaari/White Island volcanic system.