Surface wave methods gained in the past decades a primary role in many seismic projects. Specifically, they are often used to retrieve a 1D shear wave velocity model or to estimate the V
S,30
at a ...site. The complexity of the interpretation process and the variety of possible approaches to surface wave analysis make it very hard to set a fixed standard to assure quality and reliability of the results. The present guidelines provide practical information on the acquisition and analysis of surface wave data by giving some basic principles and specific suggestions related to the most common situations. They are primarily targeted to non-expert users approaching surface wave testing, but can be useful to specialists in the field as a general reference. The guidelines are based on the experience gained within the InterPACIFIC project and on the expertise of the participants in acquisition and analysis of surface wave data.
The SEIS (Seismic Experiment for Interior Structure) instrument onboard the InSight mission will be the first seismometer directly deployed on the surface of Mars. From studies on the Earth and the ...Moon, it is well known that site amplification in low-velocity sediments on top of more competent rocks has a strong influence on seismic signals, but can also be used to constrain the subsurface structure. Here we simulate ambient vibration wavefields in a model of the shallow sub-surface at the InSight landing site in Elysium Planitia and demonstrate how the high-frequency Rayleigh wave ellipticity can be extracted from these data and inverted for shallow structure. We find that, depending on model parameters, higher mode ellipticity information can be extracted from single-station data, which significantly reduces uncertainties in inversion. Though the data are most sensitive to properties of the upper-most layer and show a strong trade-off between layer depth and velocity, it is possible to estimate the velocity and thickness of the sub-regolith layer by using reasonable constraints on regolith properties. Model parameters are best constrained if either higher mode data can be used or additional constraints on regolith properties from seismic analysis of the hammer strokes of InSight’s heat flow probe HP
3
are available. In addition, the Rayleigh wave ellipticity can distinguish between models with a constant regolith velocity and models with a velocity increase in the regolith, information which is difficult to obtain otherwise.
Mobility of glaciers such as rapid retreat or disintegration of large ice volumes produces a large variety of different seismic signals. Thus, evaluating cryospheric seismic events (e.g., changes of ...their occurrence in space and time) allows to monitor glacier dynamics. We analyze a 1 year data span recorded at the Neumayer seismic network in Antarctica. Events are automatically recognized using hidden Markov models. In this study we focused on a specific event type occurring close to the grounding line of the Ekström ice shelf. Observed waveform characteristics are consistent with an initial fracturing followed by the resonance of a water‐filled cavity resulting in a so‐called hybrid event. The number of events detected strongly correlates with dominant tide periods. We assume the cracking to be driven by existing glacier stresses trough bending. Voids are then filled by seawater, exciting the observed resonance. In agreement with this model, events occur almost exclusively during rising tides where cavities are opened at the bottom of the glacier, i.e., at the sea/ice interface.
Key Points
Antarctic icequake activity correlates with dominant tide periods
Observed waveform characteristics can be explained by hybrid event model
Automatic detection and classification algorithm performs very well
On 26 December 2004, a moment magnitude Mw = 9.3 earthquake occurred along Northern Sumatra, the Nicobar and Andaman islands, resulting in a devastating tsunami in the Indian Ocean region. The rapid ...and accurate estimation of the rupture length and direction of such tsunami-generating earthquakes is crucial for constraining both tsunami wave-height models as well as the seismic moment of the events. Compressional seismic waves generated at the hypocentre of the Sumatra earthquake arrived after about 12 min at the broadband seismic stations of the German Regional Seismic Network (GRSN), located approximately 9,000 km from the event. Here we present a modification of a standard array-seismological approach and show that it is possible to track the propagating rupture front of the Sumatra earthquake over a total rupture length of 1,150 km. We estimate the average rupture speed to be 2.3-2.7 km s(-1) and the total duration of rupture to be at least 430 s, and probably between 480 and 500 s.
Celotno besedilo
Dostopno za:
DOBA, IJS, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
We apply and evaluate a recent machine learning method for the automatic classification of seismic waveforms. The method relies on Dynamic Bayesian Networks (DBN) and supervised learning to improve ...the detection capabilities at 3C seismic stations. A time-frequency decomposition provides the basis for the required signal characteristics we need in order to derive the features defining typical “signal” and “noise” patterns. Each pattern class is modeled by a DBN, specifying the interrelationships of the derived features in the time-frequency plane. Subsequently, the models are trained using previously labeled segments of seismic data. The DBN models can now be compared against in order to determine the likelihood of new incoming seismic waveform segments to be either signal or noise. As the noise characteristics of seismic stations varies smoothly in time (seasonal variation as well as anthropogenic influence), we accommodate in our approach for a continuous adaptation of the DBN model that is associated with the noise class. Given the difficulty for obtaining a golden standard for real data (ground truth) the proof of concept and evaluation is shown by conducting experiments based on 3C seismic data from the International Monitoring Stations, BOSA and LPAZ.
The H/V spectral ratio has emerged as a single station method within the seismic ambient noise analysis field by its capability to quickly estimate the frequency of resonance at a site and through ...inversion the average profile information. Although it is easy to compute from experimental data, its counter theoretical part is not obvious when building a forward model which can help in reconstructing the derived H/V spectrum. This has led to the simplified assumption that the noise wavefield is mainly composed of Rayleigh waves and the derived H/V often used without further correction. Furthermore, only the right (and left) flank around the H/V peak frequency is considered in the inversion for the subsurface 1-D shear wave velocity profile. A new theoretical approach for the interpretation of the H/V spectral ratio has been presented
by Sánchez-Sesma et al. In this paper, the fundamental idea behind their theory is presented as it applies to receivers at depth. A smooth H/V(z, f) spectral curve on a broad frequency range is obtained by considering a fine integration step which is in turn time consuming. We show that for practical purposes and in the context of inversion, this can be considerably optimized by using a coarse integration step combined with the smoothing of the corresponding directional energy density (DED) spectrum. Further analysis shows that the obtained H/V(z, f) spectrum computed by the mean of the imaginary part of Green's function method could also be recovered using the reflectivity method for a medium well illuminated by seismic sources. Inversion of synthetic H/V(z, f) spectral curve is performed for a single layer over a half space. The striking results allow to potentially use the new theory as a forward computation of the H/V(z, f) to fully invert the experimental H/V spectral ratio at the corresponding depth for the shear velocity profile (Vs) and additionally the compressional velocity profile (Vp) using receivers both at the surface and in depth. We use seismic ambient noise data in the frequency range of 0.2–50 Hz recorded at two selected sites in Germany where borehole information is also available. The obtained 1-D Vs and Vp profiles are correlated with geological log information. Results from shallow geophysical experiment are also used for comparison.
Abstract
Analyzing seismic data in a timely manner is essential for potential eruption forecasting and early warning in volcanology. Here, we demonstrate that unsupervised machine learning methods ...can automatically uncover hidden details from the continuous seismic signals recorded during Iceland’s 2021 Geldingadalir eruption. By pinpointing the eruption’s primary phases, including periods of unrest, ongoing lava extrusion, and varying lava fountaining intensities, we can effectively chart its temporal progress. We detect a volcanic tremor sequence three days before the eruption, which may signify impending eruptive activities. Moreover, the discerned seismicity patterns and their temporal changes offer insights into the shift from vigorous outflows to lava fountaining. Based on the extracted patterns of seismicity and their temporal variations we propose an explanation for this transition. We hypothesize that the emergence of episodic tremors in the seismic data in early May could be related to an increase in the discharge rate in late April.
SUMMARY
Advances in the field of seismic interferometry have provided a basic theoretical interpretation to the full spectrum of the microtremor horizontal-to-vertical spectral ratio H/V(f). The ...interpretation has been applied to ambient seismic noise data recorded both at the surface and at depth. The new algorithm, based on the diffuse wavefield assumption, has been used in inversion schemes to estimate seismic wave velocity profiles that are useful input information for engineering and exploration seismology both for earthquake hazard estimation and to characterize surficial sediments. However, until now, the developed algorithms are only suitable for on land environments with no offshore consideration. Here, the microtremor H/V(z, f) modelling is extended for applications to marine sedimentary environments for a 1-D layered medium. The layer propagator matrix formulation is used for the computation of the required Green’s functions. Therefore, in the presence of a water layer on top, the propagator matrix for the uppermost layer is defined to account for the properties of the water column. As an application example we analyse eight simple canonical layered earth models. Frequencies ranging from 0.2 to 50 Hz are considered as they cover a broad wavelength interval and aid in practice to investigate subsurface structures in the depth range from a few meters to a few hundreds of meters. Results show a marginal variation of 8 per cent at most for the fundamental frequency when a water layer is present. The water layer leads to variations in H/V peak amplitude of up to 50 per cent atop the solid layers.
Permafrost inundated since the last glacial maximum is degrading, potentially releasing trapped or stabilized greenhouse gases, but few observations of the depth of ice‐bonded permafrost (IBP) below ...the seafloor exist for most of the arctic continental shelf. We use spectral ratios of the ambient vibration seismic wavefield, together with estimated shear wave velocity from the dispersion curves of surface waves, for estimating the thickness of the sediment overlying the IBP. Peaks in spectral ratios modeled for three‐layered 1‐D systems correspond with varying thickness of the unfrozen sediment. Seismic receivers were deployed on the seabed around Muostakh Island in the central Laptev Sea, Siberia. We derive depths of the IBP between 3.7 and 20.7 m ± 15%, increasing with distance from the shoreline. Correspondence between expected permafrost distribution, modeled response, and observational data suggests that the method is promising for the determination of the thickness of unfrozen sediment.
Key Points
Ambient noise field used to detect unfrozen sediment thickness on arctic shelf
Numerical modeling used to understand H/V curves observed at the seabed
Innovative, no‐impact technique for detection of permafrost in shallow water
We have used both single-station and array methods to determine shallow shear velocity site profiles in the vicinity of the city of Cologne, Germany from ambient vibration records. Based on ...fk-analysis we assume that fundamental-mode Rayleigh waves dominate the analysed wavefield in the frequency range of 0.7–2.2 Hz. According to this view a close relation exists between H/V spectral ratios and the ellipticity of the contributing Rayleigh waves. The inversion of the shape of H/V spectral ratios then provides quantitative information concerning the local shear wave velocity structure. However, based on tests with synthetic data believed to represent a typical situation in the Lower Rhine Embayment, dispersion curves were found to provide stronger constraints on the absolute values of the velocity–depth model than the ellipticities. The shape of the ellipticities was found to be subject to a strong trade-off between the layer thickness and the average layer velocity. We have made use of this observation by combining the inversion schemes for dispersion curves and ellipticities such that the velocity–depth dependence is essentially constrained by the dispersion curves while the layer thickness is constrained by the ellipticities. In order to test this method in practice, we have used array recordings of ambient vibrations from three sites where the subsurface geology is fairly well known and geotechnical information is at least partially available. In order to keep the parameter space as simple as possible we attempted to fit only a single layer over a half-space model. However, owing to earlier studies from the region, we assume a power-law depth dependence for sediment velocities. For all three sites investigated, the inversion resulted in models for which the shear wave velocity within the sediment layer both in absolute value at the surface and in depth dependence are found to be remarkably similar to the results obtained by Budny from downhole measurements. This is strong support for the interpretation of H/V spectral ratios as Rayleigh wave ellipticities. For all three sites the predicted SH-wave site amplification factors at the fundamental frequency are of the order of 5–6 with a slightly smaller value south of Cologne.