Very-Long-Period (VLP) Seismic Artifacts during the 2018 Caldera Collapse at Kīlauea, Hawai‘i

2020 ◽  
Vol 91 (6) ◽  
pp. 3417-3432
Author(s):  
Ashton F. Flinders ◽  
Ingrid A. Johanson ◽  
Phillip B. Dawson ◽  
Kyle R. Anderson ◽  
Matthew M. Haney ◽  
...  

Abstract Throughout the 2018 eruption of Kīlauea volcano (Hawai‘i), episodic collapses of a portion of the volcano’s summit caldera produced repeated Mw 4.9–5.3 earthquakes. Each of these 62 events was characterized by a very-long-period (VLP) seismic signal (>40  s). Although collapses in the later stage of the eruption produced earthquakes with significant amplitude clipping on near-summit broadband seismometers, the first 12 were accurately recorded. For these initial collapse events, we compare average VLP seismograms at six near-summit locations to synthetic seismograms derived from displacements at collocated Global Positioning System stations. We show that the VLP seismic signal was generated by a radially outward and upward ramp function in displacement. We propose that at local distances the period of the VLP seismic signal is solely dependent on the duration of this ramp function and the instrument transfer function, that is, the seismic VLP is an artifact of the bandlimited instrument response and not representative of real ground motion. The displacement ramp function imposes a sinc-function velocity amplitude spectrum that cannot be fully recovered through standard seismic instrument deconvolution. Any near-summit VLP signals in instrument-response-corrected velocity or displacement seismograms from these collapse events are subject to severe band limitation. Similarly, the seismic amplitude response is not flat through the low-frequency corner, for example, instrument-response-uncorrected seismograms scaled by instrument sensitivity are equally prone to band limitation. This observation is crucial when attempting to clarify the different contributions to the VLP source signature. Not accounting for this effect could lead to misunderstanding of the magmatic processes involved.

2012 ◽  
Vol 707 ◽  
pp. 482-495 ◽  
Author(s):  
Ofer Manor ◽  
Leslie Y. Yeo ◽  
James R. Friend

AbstractThe classical Schlichting boundary layer theory is extended to account for the excitation of generalized surface waves in the frequency and velocity amplitude range commonly used in microfluidic applications, including Rayleigh and Sezawa surface waves and Lamb, flexural and surface-skimming bulk waves. These waves possess longitudinal and transverse displacements of similar magnitude along the boundary, often spatiotemporally out of phase, giving rise to a periodic flow shown to consist of a superposition of classical Schlichting streaming and uniaxial flow that have no net influence on the flow over a long period of time. Correcting the velocity field for weak but significant inertial effects results in a non-vanishing steady component, a drift flow, itself sensitive to both the amplitude and phase (prograde or retrograde) of the surface acoustic wave propagating along the boundary. We validate the proposed theory with experimental observations of colloidal pattern assembly in microchannels filled with dilute particle suspensions to show the complexity of the boundary layer, and suggest an asymptotic slip boundary condition for bulk flow in microfluidic applications that are actuated by surface waves.


1983 ◽  
Vol 73 (6A) ◽  
pp. 1499-1511
Author(s):  
Paul Silver

Abstract A method is proposed for retrieving source-extent parameters from far-field body-wave data. At low frequency, the normalized P- or S-wave displacement amplitude spectrum can be approximated by |Ω^(r^,ω)| = 1 − τ2(r^)ω2/2 where r^ specifies a point on the focal sphere. For planar dislocation sources, τ2(r^) is linearly related to statistical measures of source dimension, source duration, and directivity. τ2(r^) can be measured as the curvature of |Ω^(r^,ω)| at ω = 0 or the variance of the pulse Ω^(r^,t). The quantity ωc=2τ−1(r^) is contrasted with the traditional corner frequency ω0, defined as the frequency at the intersection of the low- and high-frequency trends of |Ω^(r^,ω)|. For dislocation models without directivity, ωc(P) ≧ ωc(S) for any r^. A mean corner frequency defined by averaging τ2(r^) over the focal sphere, ω¯c=2<τ2(r^)>−1/2, satisfies ωc(P) > ωc(S) for any dislocation source. This behavior is not shared by ω0. It is shown that ω0 is most sensitive to critical times in the rupture history of the source, whereas ωc is determined by the basic parameters of source extent. Evidence is presented that ωc is the corner frequency measured on actual seismograms. Thus, the commonly observed corner frequency shift (P-wave corner greater than the S-wave corner), now viewed as a shift in ωc is simply a result of spatial finiteness and is expected to be a property of any dislocation source. As a result, the shift cannot be used as a criterion for rejecting particular dislocation models.


2013 ◽  
Vol 56 (4) ◽  
Author(s):  
Paola Cusano ◽  
Simona Petrosino ◽  
Francesca Bianco ◽  
Edoardo Del Pezzo

<p>The typical earthquakes occurring at Mt. Vesuvius are Volcano-Tectonic. On July 20, 2003, an unusual earthquake with low and narrow frequency content was detected. The seismograms presented an emergent onset and a nearly monochromatic spectrum at all stations of the Osservatorio Vesuviano (Istituto Nazionale di Geofisica e Vulcanologia) seismic network. The event was located at about 4 km b.s.l. close to the crater axis and an equivalent duration magnitude of 0.6 was estimated. The nature of this event was investigated by comparing its features with those of two typical Volcano-Tectonic earthquakes occurred inside the same source volume. We compared the spectral content calculating the spectrograms and the coda patterns using the Hilbert Transform. A Seismic Moment Tensor inversion was performed on the low frequency earthquake. The focal mechanisms for the two Volcano-Tectonic earthquakes were estimated with a classical technique and resulted compatible with the stress field acting on the volcano. Taking into account the clear differences with the typical Volcano-Tectonic events as well as the peculiarities retrieved from our analyses (monochromatic, low frequency spectral content, and sustained coda) and also some geochemical observations, we classify the unusual low frequency seismic event detected at Mt. Vesuvius as Long Period earthquake and propose that its origin could be linked to a pressure drop in the deep hydrothermal system.</p>


2020 ◽  
Vol 11 (2) ◽  
pp. 155-162
Author(s):  
A. F. Sabitov ◽  
I. A. Safina

The spectral method for establishing dynamic response of measuring instruments basically requires determining the amplitude spectrum of the signal in its informative part that includes the amplitude spectrum at zero frequency. The operating frequency range of existing low-frequency spectrum analyzers is above zero frequency that leads to an uncertainty in dynamic response of measuring instruments determined by the spectral method. The purpose of this paper is to develop a program for calculating the signal amplitude spectrum, starting from zero frequency, to implement a spectral method for determining the dynamic response of measuring instruments on computers equipped with the MatLab package.To implement the spectral method for determining the dynamic response of measuring instruments, we developed a program in the MatLab 2013b environment that determines the signal amplitude spectrum from zero Hertz. The program reads the source data from Excel tables and presents the calculated amplitude spectrum as a chart and a report table.It is shown that the developed program calculates the signal amplitude spectrum with a standard deviation of not more than 3.4 % in the frequency range of 0 to 10 rad/s. The calculated amplitude spectrum allows determining the time constant of first-order aperiodic measuring instruments with an uncertainty of not more than 0.166 % at any noise level, if their frequencies are outside the information part of the spectrum.We demonstrated the claimed advantage of the spectral method for determining dynamic response using the developed program by the example of a high-frequency noise in the transient response of some measuring instruments.


2021 ◽  
Author(s):  
Leonard Seydoux ◽  
Michel Campillo ◽  
René Steinmann ◽  
Randall Balestriero ◽  
Maarten de Hoop

&lt;p&gt;Slow slip events are observed in geodetic data, and are occasionally associated with seismic signatures such as slow earthquakes (low-frequency earthquakes, tectonic tremors). In particular, it was shown that swarms of slow earthquake can correlate with slow slip events occurrence, and allowed to reveal the intermittent behavior of several slow slip events. This observation was possible thanks to detailed analysis of slow earthquakes catalogs and continuous geodetic data, but in every case, was limited to particular classes of seismic signatures. In the present study, we propose to infer the classes of seismic signals that best correlate with the observed geodetic data, including the slow slip event. We use a scattering network (a neural network with wavelet filters) in order to find meaningful signal features, and apply a hierarchical clustering algorithm in order to infer classes of seismic signal. We then apply a regression algorithm in order to predict the geodetic data, including slow slip events, from the occurrence of inferred seismic classes. This allow to (1) identify seismic signatures associated with the slow slip events as well as (2) infer the the contribution of each classes to the overall displacement observed in the geodetic data. We illustrate our strategy by revisiting the slow-slip event of 2006 that occurred beneath Guerrero, Mexico.&lt;/p&gt;


2020 ◽  
Vol 39 (7) ◽  
pp. 480-487
Author(s):  
Patrick Smith ◽  
Brandon Mattox

The P-Cable high-resolution 3D marine acquisition system tows many short, closely separated streamers behind a small source. It can provide 3D seismic data of very high temporal and spatial resolution. Since the system is containerized and has small dimensions, it can be deployed at short notice and relatively low cost, making it attractive for time-lapse seismic reservoir monitoring. During acquisition of a 3D high-resolution survey in the Gulf of Mexico in 2014, a pair of sail lines were repeated to form a time-lapse seismic test. We processed these in 2019 to evaluate their geometric and seismic repeatability. Geometric repetition accuracy was excellent, with source repositioning errors below 10 m and bin-based receiver positioning errors below 6.25 m. Seismic data comparisons showed normalized root-mean-square difference values below 10% between 40 and 150 Hz. Refinements to the acquisition system since 2014 are expected to further improve repeatability of the low-frequency components. Residual energy on 4D difference seismic data was low, and timing stability was good. We conclude that the acquisition system is well suited to time-lapse seismic surveying in areas where the reservoir and time-lapse seismic signal can be adequately imaged by small-source, short-offset, low-fold data.


2003 ◽  
Author(s):  
Chikahiro Minowa

In this paper, a new method of baseline correction on strong motion acceleration records is presented and the fundamental concept for baseline corrections on the earthquake strong motions is described. Considering the filtering effect, the earthquakes ground motion displacements of 1995 JMA KOBE, 1999 Kocaeli YPT and 1999 Chi-Chi TCU068 are discussed. Also, the linear sloshing responses of large liquid tanks subjected to these motions were discussed. Since liquid storage tanks show the low frequency (long period) sloshing characteristics and the strong motion characteristics of 1999 Kocaeli and Chi-Chi earthquakes are also low frequencies and large permanent displacements, the sloshing responses in large liquid tanks, especially in long natural periods, were significantly affected by the low frequency motions (large permanent displacements) of these devastating earthquakes. It is very important to use suitable ground motion characterized low frequency content for earthquake resistant design of liquid storage tanks. The baseline correction method presented in the paper may be adequately used to correct strong motion records for large liquid storage design.


1974 ◽  
Vol 11 (5) ◽  
pp. 691-697
Author(s):  
A. J. Wickens ◽  
H. S. Hasegawa ◽  
M. N. Bone

The Canadian Seismograph Network long period instrument parameters are used to calculate theoretical output response curves to various input functions. These curves were used as overlays to the actual response of the instruments to a known input pulse sequence. This provides an independent means of determining instrument response as well as daily monitoring ability.


2014 ◽  
Vol 14 (12) ◽  
pp. 3175-3193 ◽  
Author(s):  
P. Bottelin ◽  
D. Jongmans ◽  
D. Daudon ◽  
A. Mathy ◽  
A. Helmstetter ◽  
...  

Abstract. The eastern limestone cliff of Mount Néron (French Alps) was the theater for two medium-size rockfalls between summer and winter 2011. On 14 August 2011, a ~2000 m3 rock compartment detached from the cliff, fell 100 m below and propagated down the slope. Although most of the fallen rocks deposited on the upper part of the slope, some blocks of about 15 m in size were stopped by a ditch and an earthen barrier after a run-out of 800 m. An unstable overhanging ~2600 m3 compartment remained attached to the cliff and was blasted on 13 December 2011. During this artificially triggered event, 7 blocks reached the same ditch, with volumes ranging from 0.8 to 12 m3. A semi-permanent seismic array located about 2.5 km from the site recorded the two events, providing a unique opportunity to understand and to compare the seismic phases generated during natural and artificially triggered rockfalls. Both events have signal duration of ~100 s with comparable maximum amplitudes recorded at large distances (computed local magnitude of 1.14 and 1.05, respectively), most of the energy lying below 20 Hz. Remote sensing techniques (photogrammetry and lidar) were employed before and after the provoked rockfall, allowing the volume and fracturing to be characterized. This event was filmed by two video cameras, and the generated ground motions were recorded using two temporary 3C seismic sensors and three seismic arrays deployed at the slope toe. Videos and seismogram processing provided estimates of the propagation velocity during the successive rockfall phases, which ranges from 12 to 30 m s−1. The main seismic phases were obtained from combined video and seismic signal analyses. The two most energetic phases are related to the ground impact of fallen material after free fall, and to individual rock block impacts into the ditch and the earthen barrier. These two phases are characterized by similar low-frequency content but show very different particle motions. The discrete element technique allowed reproducing the key features of the rockfall dynamics, yielding propagation velocities compatible with experimental observations.


Sign in / Sign up

Export Citation Format

Share Document