Source spectrum and source time function of volcanic tremor determined with a dense seismic network near the summit crater of Izu-Oshima Volcano, Japan

1994 ◽  
Vol 99 (B5) ◽  
pp. 9523-9532 ◽  
Author(s):  
Jun Oikawa ◽  
Yoshiaki Ida ◽  
Koshun Yamaoka
Keyword(s):  
Volcanic Tremor ◽  
Summit Crater ◽  
Time Function ◽  
Seismic Network ◽  
Source Spectrum ◽  
Source Time Function

Effect of source spectrum on seismic attenuation measurements using the pulse‐broadening method

Geophysics ◽  
1988 ◽  
Vol 53 (12) ◽  
pp. 1520-1526 ◽  
Author(s):  
Hsi‐Ping Liu
Keyword(s):  
Wave Propagation ◽  
Internal Friction ◽  
Rise Time ◽  
Delta Function ◽  
Seismic Attenuation ◽  
Time Function ◽  
Source Spectrum ◽  
Source Time Function ◽  
Band Limited ◽  
D Wave

Numerical examples of one‐dimensional (1-D) wave propagation, using realistic source time functions in an anelastic material characterized by frequency‐independent internal friction, demonstrate that (a) the source time function strongly influences the dependencies of initial rise time τ and pulse width w on internal friction [Formula: see text] and distance x; (b) in general, τ and w have different functional dependencies on [Formula: see text] and x; and (c) the slope ∂τ/∂x for particle displacement computed for a band‐limited source time function can be either greater than or less than the corresponding value computed for a delta‐function displacement source time function. Result (a) corroborates the result for τ given by Blair and Spathis, which implies that the commonly used linear rise time and distance relation, [Formula: see text], where T is the traveltime and C a source‐independent constant, is an oversimplification of 1-D anelastic wave propagation; the source spectrum must be considered when inferring seismic attenuation from changes in rise time. Result (c) contradicts the assertion by Blair and Spathis that ∂τ/∂x of wave trains generated by a band‐limited source time function is always less than the corresponding value generated by a delta‐function source time function. By matching first arrival times and rise times obtained by modeling 1-D wave propagation in an anelastic medium to those obtained in the field, seismic attenuation can be determined.

scholarly journals Fully probabilistic seismic source inversion – Part 1: Efficient parameterisation

2013 ◽  
Vol 5 (2) ◽  
pp. 1125-1162 ◽  
Author(s):  
S. C. Stähler ◽  
K. Sigloch
Keyword(s):  
Moment Tensor ◽  
Source Parameters ◽  
Practical Importance ◽  
Principal Component ◽  
Seismic Source ◽  
Empirical Orthogonal Functions ◽  
Time Function ◽  
Source Inversion ◽  
Source Time Function ◽  
Earthquake Parameters

Abstract. Seismic source inversion is a non-linear problem in seismology where not just the earthquake parameters themselves, but also estimates of their uncertainties are of great practical importance. Probabilistic source inversion (Bayesian inference) is very adapted to this challenge, provided that the parameter space can be chosen small enough to make Bayesian sampling computationally feasible. We propose a framework for PRobabilistic Inference of Source Mechanisms (PRISM) that parameterises and samples earthquake depth, moment tensor, and source time function efficiently by using information from previous non-Bayesian inversions. The source time function is expressed as a weighted sum of a small number of empirical orthogonal functions, which were derived from a catalogue of >1000 STFs by a principal component analysis. We use a likelihood model based on the cross-correlation misfit between observed and predicted waveforms. The resulting ensemble of solutions provides full uncertainty and covariance information for the source parameters, and permits to propagate these source uncertainties into travel time estimates used for seismic tomography. The computational effort is such that routine, global estimation of earthquake mechanisms and source time functions from teleseismic broadband waveforms is feasible.

The 14 September 1995 (M = 7.3) Copala, Mexico, earthquake: A source study using teleseismic, regional, and local data

1997 ◽  
Vol 87 (4) ◽  
pp. 999-1010
Author(s):  
F. Courboulex ◽  
M. A. Santoyo ◽  
J. F. Pacheco ◽  
S. K. Singh
Keyword(s):  
Time Window ◽  
Time Function ◽  
Source Process ◽  
Source Time Function ◽  
Nonlinear Method ◽  
Local Data ◽  
Plate Interface ◽  
Linear Inversion ◽  
Field Source ◽  
Mexico Earthquake

Abstract We analyze source characteristics of the 14 September 1995, Copala, Mexico, earthquake (M = 7.3) using teleseismic, regional, and local seismograms. In the analysis of the teleseismic and the regional seismic waves, we apply the empirical Green's function (EGF) technique. The recording of an appropriate aftershock is taken as the EGF and is used to deconvolve the mainshock seismogram, thus obtaining an apparent far-field source-time function at each station. The deconvolution has been done using surface waves. For teleseismic data, we apply a spectral deconvolution method that enables us to obtain 37 apparent source-time functions (ASTFs) at 29 stations. In the analysis of the regional broadband seismograms, we use two different aftershocks as EGF, and the deconvolution is performed in the time domain with a nonlinear method, imposing a positivity constraint, and the best azimuth for the directivity vector is obtained through a grid-search approach. We also analyze two near-source accelerograms. The traces are inverted for the slip distribution over the fault plane by applying a linear inversion technique. With the aid of a time-window analysis, we obtain an independent estimation of the source-time function and a more detailed description of the source process. The analysis of the three datasets permits us to deduce the main characteristics of the source process. The rupture initiated at a depth of 16 km and propagated in two directions: updip along the plate interface toward 165° N and toward 70° N. The source duration was between 12 and 14 sec, with the maximum of energy release occurring 8 sec after the initiation of the rupture. The estimated rupture dimension of 35 × 45 km is about one-fourth of the aftershock area. The average dislocation over the fault was 1.4 m (with a maximum dislocation of 4.1 m located 10 km south of the hypocenter), which gives roughly 1 MPa as the average static stress drop.

The 2019 eruptive phase of Stromboli volcano through multiparametric geophysical observations.

2020 ◽  
Author(s):  
Massimo Orazi ◽  
Flora Giudicepietro ◽  
Carmen López ◽  
Giovanni Macedonio ◽  
Salvatore Alparone ◽  
...  
Keyword(s):  
Pyroclastic Flow ◽  
Seismic Station ◽  
Volcanic Tremor ◽  
Summit Crater ◽  
Seismic Signal ◽  
Sea Surface ◽  
Effusive Activity ◽  
Eruptive Column ◽  
Borehole Strainmeter ◽  
Sciara Del Fuoco

<p>In summer 2019, two paroxysmal explosions occurred in Stromboli. The first one occurred on July 3, when the Strombolian ordinary eruptive activity did not show a significant intensification. The explosion formed an eruptive column more than 3 km high. A pyroclastic flow ran down the “Sciara del Fuoco” slope causing a victim and some injuries. Moreover, the pyroclastic flow spread over the sea surface for about one kilometer. On August 28, a second paroxysmal explosion occurred, similar to the previous one. Also in this case the eruption formed an eruptive column of more than 3 km and a pyroclastic flow that expanded along the “Sciara del Fuoco” slope and traveled about 1 km on the sea surface. In the period between the two paroxysms, effusive activity occurred from the summit crater area. The eruptive phase of summer 2019, which began with the paroxysm of 3 July, was not preceded by significant changes in the routinely monitored parameters, such as the hourly frequency (daily average) of the VLP events (typical of Stromboli) and the amplitude of the seismic signal (RSAM). For this reason, we have analyzed the seismic and dilatometric data, which were recorded by the INGV geophysical network in the period November 2018 - September 2019, focusing our attention on other parameters that can give indications on the activity state of the volcano. In particular, we analyzed the data of the broadband seismic stations, equipped with the Guralp CMG40T sensors, and the data of one Sacks-Evertson borehole strainmeter. We defined the "VLP size", which takes into account the waveform of the VLP events, in terms of both amplitude and duration. We also applied time varying Fractal Dimension (FD) analysis to the seismograms of a seismic station close to the crater area and we analyzed the polarization of the same signal. We carried out the polarization analysis both without applying a filter and by filtering the seismic signal in the typical frequency bands of the Stromboli volcanic tremor (1-3 Hz) and of the VLPs (0.5-0.05 Hz). We found that the "VLP size", the FD and the polarization parameters showed significant changes about one month before the paroxysm of July 3. In the short term, we applied an appropriately tuned STA/LTA algorithm to the data of the borehole strainmeter, which is installed on the island at about 2km from the craters, and we obtained an automatic detection of the paroxysmal events 10 and 7.5 minutes before the explosion of July 3 and August 28, respectively.</p>

The characteristic of source spectra and stress drop of earthquakes in the Bucaramanga nest

2020 ◽  
Author(s):  
Wenzheng Gong ◽  
Xiaofei Chen
Keyword(s):  
Focal Mechanism ◽  
Stress Drop ◽  
Source Parameters ◽  
Time Function ◽  
Corner Frequency ◽  
Ratio Method ◽  
Source Time Function ◽  
Spectra Analysis ◽  
Source Spectra ◽  
Very High

<p>Spectra analysis is helpful to understand earthquake rupture processes and estimate source parameters like stress drop. Obtaining real source spectra and source time function isn’t easy, because the station recordings contain path effect and we usually can’t get precise path information. Empirical Green’s function (EGF) method is a popular way to cancel out the path effect, main two of which are the stacking spectra method (Prieto et al, 2006) and the spectral ratio method (Viegas et al, 2010; Imanishi et al, 2006). In our study, we apply the latter with multitaper spectral analysis method (Prieto et al, 2009) to calculate relative source spectra and relative source time function. Target event and EGFs must have similar focal mechanism and be collocated, so we combine correlation coefficient of wave at all stations and focal mechanism similarity to select proper EGFs.</p><p>The Bucaramanga nest has very high seismicity, so it’s suitable to calculate source spectra by using EGF method. We calculate the source spectra and source time function of about 1540 earthquakes (3-5.7ml, 135-160km depth) at Bucaramanga nest in Colombia. Simultaneously we also estimate corner frequency by fitting spectral source model (Brune, 1970; Boatwright, 1980) and stress drop using simple model (Eshelby, 1957) of earthquakes with multiple station recordings or EGFs. We obtain about 30000 events data with 12 stations from National Seismological Network of Colombia (RSNC).</p><p>The result show that the source spectra of most earthquakes fitted well by omega-square model are smooth, and the source spectra of some have obvious ‘holes’ near corner frequency, and the source time function of a few earthquakes appear two separate peeks. The first kind of earthquakes are style of self-arresting ruptures (Xu et al. 2015), which can be autonomously arrested by itself without any outside interference. Abercrombie (2014) and Wen et al. (2018) both researched the second kind of earthquakes and Wen think that this kind of earthquakes are style of the runaway ruptures including subshear and supershear ruptures. The last kind of earthquakes maybe be caused by simultaneous slip on two close rupture zone. Stress drop appear to slightly increase with depth and are very high (assuming rupture velocity/s wave velocity is 0.9). We also investigate the high-frequency falloff n, usually 2, of Brune model and Boatwright model by fitting all spectra, and find that the best value of n for Boatwright model is 2 and for Brune model is 3.5.</p>

Sign in / Sign up

Export Citation Format

Share Document