Volcano Ambient Noise Interferometry in the Caribbean (VANIC)

Author(s):  
Reinoud Sleeman

<p><span><span>The hazardous stratovolcanoes in the Lesser Antilles island arc are monitored with sparse seismic networks. The application of ambient noise interferometry to monitor seismic velocity variations (dv/v) on data from such a sparse instrumented volcanic environment often is a challenge. For the purpose of monitoring it is important a) to analyse the applicability of, and differences between, cross- and single-station cross-correlations, b) to estimate the base level of seismic velocity variations during quiet times and c) to understand the characteristics. Within the EUROVOLC instrument “Transnational Access (TA)” a proposal called VANIC was supported to a) use and evaluate different types of ambient noise cross correlations (single stations vs. multiple stations; auto, cross and cross-component correlations) to be applied on seismic recordings from the Guadeloupe seismic network on La Soufriere, b) compare the results with dv/v base level estimates from the sparse Netherlands Caribbean network on The Quill and Mt. Scenery and c) start collaboration between OVSG and KNMI on both monitoring and research levels with a focus on volcano seismology. This presentation will focus is on the results obtained during the TA visit to OVGS.</span></span></p>

2021 ◽  
Author(s):  
Fabian Lindner ◽  
Joachim Wassermann

<p>Permafrost thawing affects mountain slope stability and can trigger hazardous rock falls. As rising temperatures promote permafrost thawing, spatio-temporal monitoring of long-term and seasonal variations in the perennially frozen rock is therefore crucial in regions with high hazard potential. With various infrastructure in the summit area and population in the close vicinity, Mt. Zugspitze in the German/Austrian Alps is such a site and permafrost has been monitored with temperature logging in boreholes and lapse-time electrical resistivity tomography. Yet, these methods are expensive and laborious, and are limited in their spatial and/or temporal resolution.</p><p>Here, we analyze continuous seismic data from a single station deployed at an altitude of 2700 m a.s.l. in a research station, which is separated by roughly 250 m from the permafrost affected ridge of Mt. Zugspitze. Data are available since 2006 (with some gaps) and reveal high-frequency (>1 Hz) anthropogenic noise likely generated by the cable car stations at the summit. We calculate single-station cross-correlations between the different sensor components and investigate temporal coda wave changes by applying the recently introduced wavelet-based cross-spectrum method. This approach provides time series of the travel time relative to the reference stack as a function of frequency and lag time in the correlation functions. In the frequency and lag range of 1-10 Hz and 0.5-5 s respectively, we find various parts in the coda that show clear annual variations and an increasing trend in travel time over the past 15 years of consideration. Converting the travel time variations to seismic velocity variations (assuming homogeneous velocity changes affecting the whole mountain) results in seasonal velocity changes of up to a few percent and on the order of 0.1% decrease per year. Yet, estimated velocity variations do not scale linearly with lag time, which indicates that the medium changes are localized rather than uniform and that the absolute numbers need to be taken with caution. The annual velocity variations are anti-correlated with the temperature record from the summit but delayed by roughly one month.</p><p>The phasing of the annual seismic velocity change (relative to the temperature record) is in agreement with a previous study employing lapse-time electrical resistivity tomography. Furthermore, the decreasing trend in seismic velocity happens concurrently with an increasing trend in temperature. The results therefore suggest that the velocity changes are related to seasonal thaw and refreeze and permafrost degradation and thus highlight the potential of seismology for permafrost monitoring. By adding additional receivers and/or a fiber-optic cable for distributed acoustic sensing, hence increasing the spatial resolution, the presented method holds promise for lapse-time imaging of permafrost bodies with high spatio-temporal resolution from passive measurements.</p>


2020 ◽  
Vol 224 (1) ◽  
pp. 626-636
Author(s):  
René Steinmann ◽  
Céline Hadziioannou ◽  
Eric Larose

SUMMARY About a decade ago, noise-based monitoring became a key tool in seismology. One of the tools is passive image interferometry (PII), which uses noise correlation functions (NCF) to retrieve seismic velocity variations. Most studies apply PII to vertical components recording oceanic low-frequent ambient noise ( < 1 Hz). In this work, PII is applied to high-frequent urban ambient noise ( > 1 Hz) on three three-component sensors. With environmental sensors inside the subsurface and in the air, we are able to connect observed velocity variations with environmental parameters. Temperatures below 0 °C correlate well with strong shear wave velocity increases. The temperature sensors inside the ground suggest that a frozen layer of less than 5 cm thickness causes apparent velocity increases above 2  % , depending on the channel pair. The observations indicate that the different velocity variation retrieved from the different channel pairs are due to different surface wave responses inherent in the channel pairs. With dispersion curve modelling in a 1-D medium we can verify that surfaces waves of several tens of metres wavelength experience a velocity increase of several percent due to a centimetres thick frozen layer. Moreover, the model verifies that Love waves show larger velocity increases than Rayleigh waves. The findings of this study provide new insights for monitoring with PII. A few days with temperature below 0 °C can already mask other potential targets (e.g. faults or storage sites). Here, we suggest to use vertical components, which is less sensitive to the frozen layer at the surface. If the target is the seasonal freezing, like in permafrost studies, we suggest to use three-component sensors in order to retrieve the Love wave response. This opens the possibility to study other small-scale processes at the shallow subsurface with surface wave responses.


2021 ◽  
Author(s):  
◽  
Alexander Yates

<p>Seismic velocity changes at volcanoes carry information about stresses present within hydrothermal and magmatic systems. In this thesis, temporal velocity changes are measured at White Island volcano using ambient noise interferometry between 2007–2017. This period contains multiple well-documented eruptions starting in 2012, following an inactive period that extends back over a decade. Three primary objectives are identified: (1) investigate what seismic velocity changes can tell us about dynamic changes beneath the volcano, (2) investigate non-volcanic sources and their possible influence on interpretations, and (3) consider the potential for real-time monitoring using ambient-noise. These objectives extend beyond White Island volcano, with implications for ambient noise monitoring of volcanoes globally.  Two different approaches are used to measure velocity changes at White Island. The first involves cross-correlating noise recorded by pairs of seismic stations. Velocity changes are sought by averaging changes recorded across ten station-pairs that consist of an onshore station and a station on the volcano. The second approach involves cross-correlating the different components of individual seismic stations. This represents a less traditional approach to monitoring volcanoes, but is well-suited to White Island which has one permanent station active throughout eruptive activity. Single seismic stations located onshore are also processed to investigate background regional changes.  Two periods of long-term velocity increases are detected at the volcano. The first occurs during a highly active period in 2012–2013 and the second occurs in the months preceding an explosive eruption in April 2016. Comparison with velocities recorded by onshore stations suggest a meteorological source for these changes is unlikely. Velocity increases are therefore interpreted to reflect cracks closing under increased pressures beneath the volcano. Similarly, a rapid decline in the velocity within 2–3 months of the April 2016 eruption is interpreted to reflect depressurization of the system.  In addition to volcanic sources, we also find clear evidence of non-volcanic processes influencing velocity changes at the volcano. Two clear co-seismic velocity decreases of approximately 0.05–0.1% are associated with a Mw 5.2 earthquake in 2008 — within 10 km of the volcano — and the Mw 7.1 East Cape earthquake in 2016. The East Cape earthquake — located 200 km away from the volcano — produces significant velocity decreases over a large region, as detected by stations onshore and on White Island. This likely reflects dynamic stress changes as a result of passing seismic waves, with an eruption two weeks later interpreted here to have been triggered by this event. Finally, we identify similarities between annual variations recorded by onshore stations and changes at the volcano, suggesting an environmental influence. Velocity changes at White Island therefore represent a complex interaction of volcanic and non-volcanic processes, highlighting the need for improved understanding of external sources of change to accurately detect short-term eruptive precursors.</p>


2022 ◽  
Author(s):  
◽  
Weiwei Wang

<p><b>This thesis uses continuous ambient noise data recorded by Ocean Bottom Seismometers (OBSs) to study seismic velocities in the upper crust of the overriding plate. The first and second projects (Chapters 3 and 4) focus on temporal seismic velocity variations in the northern Hikurangi subduction zone offshore the North Island, New Zealand, while the third project (Chapter 5) focuses on shear wave velocities in the southwestern Okinawa Trough offshore northeastern Taiwan. In the first project (Chapters 3), we investigate a region of frequent slow slip events (SSEs) offshore Gisborne, North Island, New Zealand. From September to October 2014, an SSE occurred with a slip over 250 mm and was recorded successfully by the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip deployment II (HOBITSS II). We apply coda wave interferometry on the ambient noise data acquired by nine OBSs deployed by the HOBITSS II to study the seismic velocity variations related to the SSE. The average velocity variations display a decrease on the order of 0.05% during the SSE, followed by an increase of similar magnitude afterwards. Two hypotheses are proposed to explain our observation. The first hypothesis, which has been suggested by previous studies, considers that the velocity decrease during the SSE is caused by more fluids migrating into the upper plate as the SSE breaks a low-permeability seal on the plate boundary. After the SSE, the fluids in the upper plate diffuse gradually and the velocity increases; The second hypothesis is that before the SSE, elastic strain accumulates causing contraction and reduction of porosity and therefore increase of velocity (the velocity increase between SSEs). During the SSE, the velocity decrease is caused by increased porosity as the SSE relieves the accumulated elastic strain on the plate interface, which results in dilation. After the SSE, stress and strain accumulate again, causing a porosity decrease and a velocity increase back to the original value. This study demonstrates that the velocity variations related to SSEs are observable and provides evidence for slow slip mechanism hypotheses.</b></p> <p>The second project (Chapter 4) focuses on the temporal seismic velocity variations associated with an SSE in 2019 offshore Gisborne, North Island, New Zealand. This is a later SSE in the same area as the first project (Chapters 3). Based on the success of the HOBITSS II, more ocean bottom instruments were deployed in the northern Hikurangi subduction zone from 2018 to 2019 (HOBITSS V). An SSE lasting approximately one month from the end of March to the beginning of May 2019 occurred during the deployment and was recorded by the network. The main slip was south of the deployment and the slip beneath the deployment was up to 150 mm. This study applies coda wave interferometry on the ambient noise data acquired by five OBSs and computes seismic velocity variations to investigate their relation to the SSE. A velocity decrease on the order of 0.015% during the SSE and an increase back to the original velocity value are observed at 1–2.5 s. This supports the two hypotheses proposed in Chapters 3: fluid migration and/or stain changes through the SSE cycle. In addition, velocity variations computed from individual stations show velocity increases before the SSE, which are destructively interfered in their average. Such a situation could occur if the SSE migrated across the network. If the velocity increases before the SSE from individual stations are real, they can be only explained by the hypothesis of crustal strain changes (the second hypothesis in project 1). However, fluid migration (the first hypothesis in project 1) may still happen concomitantly.</p> <p>The third project focuses on the tectonics in southwestern Okinawa Trough offshore northeastern Taiwan. The southwestern Okinawa Trough is an active back-arc basin, extending and rifting within the continental lithosphere. The tectonic development of the back-arc basin is still not well-understood. This study uses continuous ambient noise data recorded by 34 OBSs deployed by Academia Sinica at various periods from 2010 to 2018. Cross-correlations on vertical seismic components and pressure gauges are computed to construct Rayleigh/Scholte waves to study the shear wave velocity structure in the southwestern Okinawa Trough. Phase velocities are measured from the Rayleigh/Scholte waves. Shear velocities are inverted from the phase velocities. Results show the velocity in the south of the back-arc rifting axis near the axis is slower than the velocity in the north of the rifting axis, suggesting the velocity structure in the southwestern Okinawa Trough is asymmetric along the rifting axis. Previous studies have shown high heat flows (about 110mW/m 2 on average) in the south of the rifting axis. The low velocity in the south could be caused by the high heat flow that may be related to asymmetric back-arc extension and/or rifting. This study presents the shear wave velocity structure in the southwest Okinawa Trough is asymmetric along the rifting axis, which implies the back-arc extending/rifting is asymmetric in the study region. This study also suggests effective techniques for OBS noise corrections and unwrapping the cycle skipping of phase velocity measurements.</p> <p>In summary, this thesis represents three projects focusing on seismic velocities in two subduction zones using ambient noise data collected by OBSs. The first and second projects study the temporal velocity variations and the relation to SSEs. Both studies observe velocity decreases during the SSEs and increases after the SSEs, supporting two hypotheses of fluid migration and/or stain changes through the SSE cycle. The third project finds the shear velocity structure in the southwestern Okinawa Trough is asymmetric along the rifting center, which may imply the back-arc extension is asymmetric.</p>


2016 ◽  
Vol 206 (3) ◽  
pp. 1795-1813 ◽  
Author(s):  
Zack Spica ◽  
Mathieu Perton ◽  
Marco Calò ◽  
Denis Legrand ◽  
Francisco Córdoba-Montiel ◽  
...  

2019 ◽  
Vol 24 (1) ◽  
pp. 121-131 ◽  
Author(s):  
Sang-Jun Lee ◽  
Junkee Rhie ◽  
Seongryong Kim ◽  
Tae-Seob Kang ◽  
Chang Soo Cho

AbstractMonitoring seismic activity in the north Korean Peninsula (NKP) is important not only for understanding the characteristics of tectonic earthquakes but also for monitoring anthropogenic seismic events. To more effectively investigate seismic properties, reliable seismic velocity models are essential. However, the seismic velocity structures of the region have not been well constrained due to a lack of available seismic data. This study presents 1-D velocity models for both the inland and offshore (western East Sea) of the NKP. We constrained the models based on the results of a Bayesian inversion process using Rayleigh wave dispersion data, which were measured from ambient noise cross-correlations between stations in the southern Korean Peninsula and northeast China. The proposed models were evaluated by performing full moment tensor inversion for the 2013 Democratic People’s Republic of Korea (DPRK) nuclear test. Using the composite model consisting of both inland and offshore models resulted in consistently higher goodness of fit to observed waveforms than previous models. This indicates that seismic monitoring can be improved by using the proposed models, which resolve propagation effects along different paths in the NKP region.


2021 ◽  
Author(s):  
◽  
Adrian Shelley

<p>This thesis is concerned with scrutinising the source, distribution and detectability of seismic velocity phenomena that may be used as proxies to study conditions in the crust. Specifically, we develop modelling techniques in order to analyse the directional variation of seismic wave speed in the crust and test them at Mt. Asama in Japan and Canterbury, New Zealand. We also implement both active source and noise interferometry to identify velocity variations at Mt. Ruapehu, New Zealand.  Observations of temporal variation of anisotropic seismic velocity parameters at Asama volcano in Japan indicate that there is some process (or processes) affecting anisotropy, attributed to closure of microcracks in the rock as it is subjected to volcanic stress in the crust. To test this assertion, a 3D numerical model is created incorporating volcanic stress, ray tracing and estimation of the anisotropy to produce synthetic shear wave splitting results using a dyke stress model. Anisotropy is calculated in two ways; by considering a basic scenario where crack density is uniform and a case where the strength of anisotropy is related to dry crack closure from deviatoric stress. We find that the approach is sensitive to crack density, crack compliance, and the regional stress field. In the case of dry crack closure, modelled stress conditions produce a much smaller degree of anisotropy than indicated by measurements. We propose that the source of anisotropy changes at Asama is tied to more complex processes that may precipitate from stress changes or other volcanic processes, such as the movement of pore fluid.  We develop a generalised anisotropy inversion model based on the linearised, iterative least-squares inversion technique of Abt and Fischer [2008]. The model is streamlined for use with results from the MFAST automatic shear wave splitting software [Savage et al., 2010]. The method iteratively solves for the best fitting magnitude and orientation of anisotropy in each element of the model space using numerically calculated partial derivatives. The inversion is applied to the Canterbury plains in the region surrounding the Greendale fault, using shear-wave splitting data from the 2010 Darfield earthquake sequence. Crustal anisotropy is resolved down to a depth of 20 km at a spatial resolution of 5 km, with good resolution near the Greendale fault. We identify a lateral variation in anisotropy strength across the Greendale fault, possibly associated with post-seismic stress changes.  We perform active source and noise interferometry at Ruapehu in order to investigate potential seismic velocity changes and assess their use as a possible eruption forecasting method. Six co-located 100 kg ammonium nitrate fuel oil explosives were set off serially at Lake Moawhango, situated approximately 20 km south-east of Mount Ruapehu. Two methods of interferometry, using moving window cross correlation in the time and frequency domains, respectively, were applied to the recorded signal from each explosion pair in order to determine velocity changes from the signal coda waves. We identify possible diurnal velocity variations of ~ 0:7% associated with strain caused by the solid Earth tide. Synthetic testing of velocity variation recoverability was also performed using both methods. Interferometry of noise cross-correlations during the period was also performed using moving window cross correlation in the frequency domain. Analysis of velocity variations in the ZZ, RR and TT component pairs show little coherency. This, combined with results from synthetic testing that show that the frequency domain interferometry technique employed is unstable above velocity variations of 0.1%, indicate that the method may not be suitable for determining velocity variations at Ruapehu.</p>


2019 ◽  
Vol 5 (11) ◽  
pp. eaax6642 ◽  
Author(s):  
C. Donaldson ◽  
T. Winder ◽  
C. Caudron ◽  
R. S. White

Seismic noise interferometry is an exciting technique for studying volcanoes, providing a continuous measurement of seismic velocity changes (dv/v), which are sensitive to magmatic processes that affect the surrounding crust. However, understanding the exact mechanisms causing changes in dv/v is often difficult. We present dv/v measurements over 10 years in central Iceland, measured using single-station cross-component correlation functions from 51 instruments across a range of frequency bands. We observe a linear correlation between changes in dv/v and volumetric strain at stations in regions of both compression and dilatation associated with the 2014 Bárðarbunga-Holuhraun dike intrusion. Furthermore, a clear seasonal cycle in dv/v is modeled as resulting from elastic and poroelastic responses to changing snow thickness, atmospheric pressure, and groundwater level. This study comprehensively explains variations in dv/v arising from diverse crustal stresses and highlights the importance of deformation modeling when interpreting dv/v, with implications for volcano and environmental monitoring worldwide.


Sign in / Sign up

Export Citation Format

Share Document