Very-small-aperture 3D infrasonic array for volcanic jet observation at Stromboli Volcano

Summary This study tested a very-small-aperture three-dimensional (VSA-3D) infrasonic array. A 3D array is ideal for resolving the back elevation angle (BEL), which has become important in the analysis of volcanic jet noise or geologic flows on steep mountain slopes. Although a VSA infrasonic array, with an aperture as small as a few tens of meters, has recently been shown to have a sufficient resolution of the back azimuth (BAZ) of incident signals, its BEL resolution is considered to be poor. We performed a four-element 3D array experiment with a 20-m aperture and 2-m height at the summit of Stromboli Volcano. We analyzed the direction of arrival (DOA) with the MUSIC algorithm as a function of frequency and conducted a cluster analysis for the estimated DOA–frequency functions of eruption signals. As a result, individual infrasonic signals were successfully related to eruptive vents. We also calculated the standard deviation (STD) of the DOAs in each cluster. Of the observed BAZ-STDs and BEL-STDs, 80 per cent were <2.0° and <4.6°, respectively. A comparison among the array geometries showed that the installation of a sensor above the ground, even at only 2 m, improved the BEL resolution, indicating that the VSA-3D array provides more detailed information about the wavefield than a planar array. The observed signals had higher BELs (−20° to 0°) than the vent direction (−30° to −25°) at 3–6 Hz, although signals above 20 Hz arrived from the vent direction. Our array verified that such DOA deviations were significant by the STD analysis and some tests with synthetic data. We infer that the DOA deviations do not indicate the source location and are caused by topographical diffraction.

Off-Network Earthquake Location by Earthquake Early Warning Systems: Methodology and Validation

ABSTRACT Regional source-based earthquake early warning systems perform three consecutive tasks: (1) detection and epicenter location, (2) magnitude determination, and (3) ground-motion prediction. The correctness of the magnitude determination is contingent on that of the epicenter location, and the credibility of the ground-motion prediction depends on those of the epicenter location and the magnitude determination. Thus, robust epicenter location scheme is key for regional earthquake early warning systems. Available source-based systems yield acceptably accurate locations when the earthquakes occur inside the real-time seismic network, but they return erroneous results otherwise. In this study, a real-time algorithm that is intended as a supplement to an existing regional earthquake early warning systems is introduced with the sole objective of ameliorating its off-network location capacity. The new algorithm combines measurements from three or more network stations that are analyzed jointly using an array methodology to give the P-wave slowness vector and S-phase arrival time. Prior to the S-phase picking, the nonarrival of the S phase is used for determining a minimum epicentral distance. This estimate is updated repeatedly with elapsed time until the S phase is picked. Thus, the system timeliness is not compromised by waiting for the S-phase arrival. After the S wave is picked, an epicentral location can be determined using a single array by intersecting the back-azimuth beam with the S-minus-P annulus. When several arrays are assembled, the back azimuth and P and S picks from all arrays are combined to constrain the epicenter. The performance of the array processing for back azimuth and S-wave picking is assessed using a large number of accelerograms, recorded by nine strong motion sensors of the KiK-net seismic network in Japan. The nine stations are treated as three distinct seismic arrays, comprising three stations each. Good agreement is found between array-based and catalog-reported parameters. Finally, the advantage of the new array methodology with respect to alternative schemes for back azimuth and distance is demonstrated.

Estimating Rupture Front of Large Earthquakes Using a Novel Multi-Array Back-Projection Method

A ruptured front obtained from high-frequency energy radiation is the key to understand the complex source. It is commonly observed that rupture fronts derived from different arrays often show some variations due to the obvious difference of the positioning accuracy of the far-field array between the azimuth and the epicentral distance. We developed a new multi-array back-projection method based on the classical back-projection method and applied the method to the 2015 MW7.8 Nepal earthquake. The back azimuth information with small error is separated from the classical back-projection results, and the azimuth intersection of multiple arrays is used to obtain more accurate spatial and temporal distribution information of the source rupture fronts.

Benchmarking microbarom radiation and propagation model against infrasound recordings: a vespagram-based approach

This study investigates the use of a vespagram-based approach as a tool for multi-directional comparison between simulated microbarom soundscapes and infrasound data recorded at ground-based array stations. Data recorded at the IS37 station in northern Norway during 2014–2019 have been processed to generate vespagrams (velocity spectral analysis) for five frequency bands between 0.1 and 0.6 Hz. The back azimuth resolution between the vespagram and the microbarom model is harmonized by smoothing the modeled soundscapes along the back azimuth axis with a kernel corresponding to the frequency-dependent array resolution. An estimate of similarity between the output of the microbarom radiation and propagation model and infrasound observations is then generated based on the image-processing approach of the mean square difference. The analysis reveals that vespagrams can monitor seasonal variations in the microbarom azimuthal distribution, amplitude, and frequency, as well as changes during sudden stratospheric warming events. The vespagram-based approach is computationally inexpensive, can uncover microbarom source variability, and has the potential for near-real-time stratospheric diagnostics and atmospheric model assessment.

Crustal seismic structure and anisotropy of Madagascar and south-eastern Africa using receiver function harmonics: interplay of inherited local heterogeneities and current regional stress

Summary This study investigates the seismic structure and anisotropy in the crust beneath Madagascar and south-eastern Africa, using receiver functions. The understanding of seismic anisotropy is essential for imaging past and present deformation in the lithosphere-asthenosphere system. In the upper mantle, seismic anisotropy mainly results from the orientation of olivine, which deforms under tectonic (fossil anisotropy) or flow processes (in the asthenosphere). In the crust, the crystallographic alignment of amphiboles, feldspars(plagioclase) or micas or the alignment of heterogeneities such as fractures, add to a complex geometry, which results in challenges to understanding the Earth's shallow structure. The decomposition of receiver functions into back-azimuth harmonics allows to characterize orientations of lithospheric structure responsible for azimuthally-varying seismic signals, such as a dipping isotropic velocity contrasts or layers of azimuthal seismic anisotropy. By analysing receiver function harmonics from records of 48 permanent or temporary stations this study reveals significant azimuthally-varying signals within the upper crust of Madagascar and south-eastern Africa. At 30 stations crustal anisotropy dominates the harmonics while the signature of a dipping isotropic contrast is dominant at the remaining 18 stations. However, all stations’ back-azimuth harmonics show complex signals involving both dipping isotropic and shallow anisotropic contrasts or more than one source of anisotropy at shallow depth. Our calculated orientations for the crust are therefore interpreted as reflecting either the average or the interplay of several sources of azimuthally-varying signals depending of their strength. However, comparing information between stations allows us to draw the same conclusions regionally: in both southern Africa and Madagascar our measurements reflect the interplay between local, inherited structural heterogeneities and crustal seismic anisotropy generated by the current extensional stress field imposed by the southward propagation of the East-African Rift System. A final comparison of our crustal orientations with SKS orientations attributed to mantle deformation further probes the interplay of crustal and mantle anisotropy on SKS measurements.

Seismically Derived Ground Tilts Related to the 2010 Chilean Tsunami

Developing new ways to observe tsunami contributes to tsunami research. Tidal and deep-ocean gauges are typically used for coastal and offshore observations. Recently, tsunami-induced ground tilts offer a new possibility. The ground tilt signal accompanied by 2010 Mw 8.8 Chilean earthquake were observed at a tiltmeter network in Japan. However, tiltmeter stations are usually not as widely installed as broadband seismometers in other countries. Here, we studied broadband seismic records from Japan’s F-net and found ground tilt signals consistent with previously published tiltmeter dataset for this particular tsunamic event. Similar waveforms can also be found in broadband seismic networks in other countries, such as Taiwan, as well as an ocean-bottom seismometer. We documented a consistent time sequence of evolving back-azimuth directions of the tsunami waves at different stages of tsunami propagation through beamforming-frequency–wavenumber analysis and particle-motion analysis; the outcomes are consistent with the tsunami propagation model provided by the Pacific Tsunami Warning Center. These results shown that dense broadband seismic networks can provide a useful complementary dataset, in addition to tiltmeter arrays and other networks, to study or even monitor tsunami propagation using arrayed methods.

Seismic Observations of Four Thunderstorms Using an Underground Fiber-Optic Array

Thunderstorms are a common atmospheric phenomenon that cause abundant acoustic disturbances, which can interact with the ground surface, creating a link between atmospheric and solid Earth processes. This article reports seismological observations of four thunderstorms through the spring and summer of 2019, as recorded by the distributed acoustic sensing fiber-optic array (4.9 km) on the Penn State campus in State College, Pennsylvania. With a dense sensor array in the local region, we are able to construct the seismic full waveform response of the thunderstorm events (hereafter referred to as thunderquakes) and track the wave propagation across the array. We use a time-domain grid search to obtain the back azimuth and slowness of the waves, and a modified Geiger’s method to pinpoint source locations of the thunderquakes. Correlated with the time of the recorded signal, this data allows reconstruction of thunderstorm movement as well as offering measurements of the seismic velocity.

Distributed Acoustic Sensing Using Dark Fiber for Array Detection of Regional Earthquakes

The intrinsic array nature of distributed acoustic sensing (DAS) makes it suitable for applying beamforming techniques commonly used in traditional seismometer arrays for enhancing weak and coherent seismic phases from distant seismic events. We test the capacity of a dark-fiber DAS array in the Sacramento basin, northern California, to detect small earthquakes at The Geysers geothermal field, at a distance of ∼100 km from the DAS array, using beamforming. We use a slowness range appropriate for ∼0.5–1.0 Hz surface waves that are well recorded by the DAS array. To take advantage of the large aperture, we divide the ∼20 km DAS cable into eight subarrays of aperture ∼1.5–2.0 km each, and apply beamforming independently to each subarray using phase-weighted stacking. The presence of subarrays of different orientations provides some sensitivity to back azimuth. We apply a short-term average/long-term average detector to the beam at each subarray. Simultaneous detections over multiple subarrays, evaluated using a voting scheme, are inferred to be caused by the same earthquake, whereas false detections caused by anthropogenic noise are expected to be localized to one or two subarrays. Analyzing 45 days of continuous DAS data, we were able to detect all earthquakes with M≥2.4, while missing most of the smaller magnitude earthquakes, with no false detections due to seismic noise. In comparison, a single broadband seismometer co-located with the DAS array was unable to detect any earthquake of M<2.4, many of which were detected successfully by the DAS array. The seismometer also experienced a large number of false detections caused by spatially localized noise. We demonstrate that DAS has significant potential for local and regional detection of small seismic events using beamforming. The ubiquitous presence of dark fiber provides opportunities to extend remote earthquake monitoring to sparsely instrumented and urban areas.

Performance of a rotational sensor at Etna, Italy focusing on back-azimuth estimations of volcano-seismic events

<p>The field of rotational seismology has only recently emerged. Portable 3 component rotational sensors are commercially available since a few years which opens the pathway for a first use in volcano-seismology. The combination of rotational and translational components of the wavefield allows identifying and filtering for specific seismic wave types, estimating the back azimuth of an earthquake, and calculating local seismic phase velocities.</p><p>Our work focuses on back-azimuth calculations of volcano-tectonic and long-period events detected at Etna volcano in Italy. Therefore, a continuous full seismic wavefield of 30 days was recorded by a BlueSeis-3A, the first portable rotational sensor, and a broadband Trillium Compact seismometer located next to each other at Mount Etna in August and September of 2019. In this study, we applied two methods for back-azimuth calculations. The first one is based on the similarity of the vertical rotation rate to the horizontal acceleration and the second one uses a polarization analysis from the two horizontal components of the rotation rate. The estimated back-azimuths for volcano-tectonic events were compared to theoretical back-azimuths based on the INGV event catalog and the long-period event back-azimuths were analyzed for their dominant directions. We discuss the quality of our back azimuths with respect to event locations and evaluate the sensitivity and benefits of the rotational sensor focusing on volcano-seismic events on Etna regarding the signal to noise ratios, locations, distances, and magnitudes.</p>

The use of seismic ground particle motion for snow avalanche release area identification

<p>Knowledge of the snow avalanche release area is key information in snow avalanche studies. However, it is not easy to obtain from a remote location. The study of the seismic vibrations produced in the initial stages of the snow avalanche, makes possible to identify their origin and to link them to the starting area of the snow avalanche. We developed a methodology for this purpose, applied to seismic data acquired from a 3D seismic station (2Hz eigenfrequency) placed at Cavern A in Vallée de la Sionne experimental site (VDLS, WSL-SLF), deployed in 2013 by UB-RISKNAT. This is the closest position to the snow avalanche release areas, at 700 m to the farthest point. We focus on spontaneous triggered snow avalanches to achieve better signal-to-noise ratio and to be more realistic on its application.</p><p>For the isolation of the Signal Onset (SON) section of seismic data, which corresponds to those vibrations produced by the initial stage of the snow avalanche, we use the STA/LTA ratios and seismic signal amplitude, common methodologies in seismology. The STA/LTA is used for the identification of the first vibrations produced by the movement of the snow mass and the seismic signal amplitude thresholds for the identification of the end of the SON section -when the snow avalanche front reaches the seismic sensor position-. The 3D seismic data [ZNE components] of the SON section were processed in time windows. The study of polarization of the particle motion to obtain the direction of the back-azimuth of the signal (Vidale, 1986; Jurckevicks, 1988) was carried out for each time window of the seismic signal. The accumulation of back-azimuth directions for the entire SON section is related to the origin of the vibrations and, by extension, to the snow avalanche release area.</p><p>The entire algorithm has been automated. In its application on all the trigger activations at VDLS since 2015 until 2020, it was achieved a success rate of 78% on snow avalanche release area identification. In addition, we defined an algorithm based on STA/LTA ratio to select the snow avalanches from other seismic events, used with a success rate of 95%.</p><p>We present the application of our method in a case study, a large spontaneous snow avalanche released on 16th February 2018 at VDLS. The snow avalanche had two main release areas, clearly identified in photos of the site. The two developed fronts can be recognized in the seismic data. The directions to the release areas from Cavern A position can be identified using the presented method. Also, more interpretations can be done on the downhill snow avalanche path.</p>