Anomalous surface waves from Lop Nor nuclear explosions: Observations and numerical modeling

Surface waves from the Chinese test site of Lop Nor are analyzed using long-period and broadband stations located at regional and teleseismic distances and at different azimuths. For most azimuths, strong Love waves between 0.02 and 0.045 Hz are observed with an amplitude of up to 10 times that of the Rayleigh waves. In addition, an anomalous early Rayleigh wave train is observed at some stations in western Europe. Due to a particularly favorable station and source configuration, it is possible to isolate the areas where the anomalies are created. The high-amplitude Love waves must be attributed to either source effects or path effects immediately north of Lop Nor. The early wave train is shown to be due to a partial energy conversion between Love and Rayleigh waves, probably at the Tornquist Zone. To estimate the possible contribution from surface wave conversions to the observed anomalies, numerical simulations are carried out with the indirect boundary element method. The simulations show that a relatively small variation of crustal thickness can induce Rayleigh to Love wave conversions between 0.02 and 0.1 Hz frequency. The calculated amplitudes of the Love waves are significant (up to 35% of the amplitude of the incoming Rayleigh waves), but they are too small to fit the observed amplitude anomaly. The observed converted waves and the numerical results nevertheless indicate that surface wave conversions can be significant across strong lateral crustal heterogeneities. In particular, the conversions due to changes in crustal thickness are located in the period interval which is routinely used for estimation of Ms

Site Amplification Effects in the Ubaye Valley (France): Measurements and Modeling

Local soil conditions can yield to significant variations in ground motion generated by earthquakes. During the last two decades, these amplification effects have been observed by numerous authors for well-documented earthquakes. On the other hand, theoretical models and numerical techniques have been proposed to physically understand site effects. However, few comparisons have been made between observations and theoretical results for well-known 2D structures (sediment-filled valleys). The aim of this paper is to present a study of the response of a valley in the French Alps. The structure (Ubaye valley) was chosen for its moderate dimensions (500 meters wide and 65 meters thick) which allow an accurate determination of the deposit characteristics, and for the relatively high seismicity of the region. The study has included the set-up of a temporary array of five seismological stations, a geophysical survey of the valley to determine the dynamic properties and the geometry of the soft deposits, and numerical modeling (1D and 2D cases) of the response. Comparisons between observed and computed amplifications show a good agreement for particular input motions (corresponding to the SH case). In the other cases, the spectral ratios exhibit a great variability between different groups of similar earthquakes and more developed simulations should be used (2D P-SV, 3D)

Emergence of body waves from cross-correlation of short period seismic noise.

International audienceAmbient noise correlation is now widely used in seismology to obtain the surface waves part of Green's function. More difficult is the extraction of body waves from noise correlations. Using 42 temporary broad-band three components stations located on the northern part of the fennoscandian region, we identify high-frequency (0.5-2 Hz) body waves emerging from noise correlations for inter-station distances up to 550 km. The comparison of the noise correlations with earthquake data confirms that the observed waves can be interpreted as P and S waves reflected from the Moho. Because the crustal model of the area is well known, we also compared the noise correlations with synthetic seismograms, and found an excellent agreement between the travel times of all the observed phases. Polarization analysis provides a further argument to confirm the observation of body waves

Stability of Monitoring Weak Changes in Multiply Scattering Media with Ambient Noise Correlation: Laboratory Experiments

Previous studies have shown that small changes can be monitored in a scattering medium by observing phase shifts in the coda. Passive monitoring of weak changes through ambient noise correlation has already been applied to seismology, acoustics and engineering. Usually, this is done under the assumption that a properly reconstructed Green function as well as stable background noise sources are necessary. In order to further develop this monitoring technique, a laboratory experiment was performed in the 2.5MHz range in a gel with scattering inclusions, comparing an active (pulse-echo) form of monitoring to a passive (correlation) one. Present results show that temperature changes in the medium can be observed even if the Green function (GF) of the medium is not reconstructed. Moreover, this article establishes that the GF reconstruction in the correlations is not a necessary condition: the only condition to monitoring with correlation (passive experiment) is the relative stability of the background noise structure

Lunar subsurface investigated from correlation of seismic noise.

International audienceBy correlating seismic noise recorded by four sensors placed on the Moon during the Apollo 17 mission, we have retrieved a well-defined dispersed Rayleigh wave pulse. Inversion of its group velocity provides new constraints on the lunar subsurface structure. The estimated ”signal-to-noise” ratio (SNR) of the retrieved Rayleigh wavetrain is strongly dependent on solar illumination, effectively making solar heating a source of seismic noise on the Moon. This result suggests that in future planetary missions it is feasible to extract information on the internal structure of extraterrestrial objects by correlating seismic noise even when natural quakes are absent

Analysis of micro-seismicity in sea ice with deep learning and Bayesian inference: application to high-resolution thickness monitoring

In the perspective of upcoming seasonally ice-free Arctic, understanding the dynamics of sea ice in the changing climate is a major challenge in oceanography and climatology. In particular, the new generation of sea ice models will require fine parameterization of sea ice thickness and rheology. With the rapidly evolving state of sea ice, achieving better accuracy, as well as finer temporal and spatial resolutions of its thickness will set new monitoring standards, with major scientific and geopolitical implications. Recent studies have shown the potential of passive seismology to monitor the thickness, density and elastic properties of sea ice with significantly reduced logistical constraints. For example, human intervention is no longer required, except to install and uninstall the geophones. Building up on this approach, we introduce a methodology for estimating sea ice thickness with high spatial and temporal resolutions from the analysis of icequakes waveforms. This methodology is based on a deep convolutional neural network for automatic clustering of the ambient seismicity recorded on sea ice, combined with a Bayesian inversion of the clustered waveforms. By applying this approach to seismic data recorded in March 2019 on fast ice in the Van Mijen fjord (Svalbard), we observe the spatial clustering of icequakes sources along the shore line of the fjord. The ice thickness is shown to follow an increasing trend that is consistent with the evolution of temperatures during the four weeks of data recording. Comparing the energy of the icequakes with that of calibrated seismic sources, we were able to derive a power law of icequake energy, and to relate this energy to the size of the cracks that generate the icequakes.</p

Fluctuations of correlations and Green's function reconstruction: role of scattering

Correlations of ambient seismic or acoustic vibrations are now widely used to reconstruct the impulse response between two passive receivers as if a source was placed at one of them. This provides the opportunity to do imaging without a source, or \textsl{passive imaging}. Applications include terrestrial and solar seismology, underwater acoustics, and structural health monitoring, to cite only a few. Nevertheless, for a given set of data, correlations do not only yield the Green's function between the sensors. They also contain residual fluctuations that result from an imperfect time or source averaging that might eventually blur the images. In this article, we propose a heuristic model to describe the level of fluctuations of the correlations in the case of non-stationary wavefields, and more particularly in the case of scattering media. The work includes theoretical derivations and numerical simulations. The role of multiple scattering is quantitatively evaluated. The level of fluctuations decreases when the duration and intensity of the diffuse waves increase. The role of absorption is also discussed: absorption is properly retrieved by correlation, but the level of fluctuations is greater, thus degrading the Green's function reconstruction. Discrepancies of our simple model in the case of strong multiple scattering ($k\ell^*\leq 18$) are discussed

On the precision of noise correlation interferometry.

International audienceLong duration noisy-looking waveforms such as those obtained in randomly multiple scattering and reverberant media are complex; they resist direct interpretation. Nevertheless, such waveforms are sensitive to small changes in the source of the waves or in the medium in which they propagate. Monitoring such waveforms, whether obtained directly or obtained indirectly by noise correlation, is emerging as a technique for detecting changes in media. Interpretation of changes is in principle problematic; it is not always clear whether a change is due to sources or to the medium. Of particular interest is the detection of small changes in propagation speeds. An expression is derived here for the apparent, but illusory, waveform dilation due to a change of source. The expression permits changes in waveforms due to changes in wave speed to be distinguished with high precision from changes due to other reasons. The theory is successfully compared with analysis of a laboratory ultrasonic data set and a seismic data set from Parkfield California

On the dynamic sliding with friction of a rigid block and of an infinite elastic slab

International audienceWe consider dynamic motions of two elastic systems undergoing frictional slip. The first one is the classical model of the frictional slider loaded through an elastic spring. The second one is an infinite elastic slab bounded by two planes which is in contract with a rigid body and submitted to shearing. Slip weakening and slip rate weakening friction laws are both considered. The two simple systems show very different qualitative behaviors. In the case of the slip dependent friction a slider moves with a single slip event when a critical stress level is reached. Under the same conditions, a series of slip events occur for the infinite slab. This difference between the behavior of the two systems is due to the important part played by inertia in the mass concentrated block slider model. In the case of slip rate weakening, the analysis of the problem for the infinite slab indicates a major difficulty: this problem has no unique solution if the rate of weakening exceeds a limit that is explicitly given. Whatever is the selection rule chosen to discriminate the solution, shocks will occur. The slip history obtained for the slab is very different from the one obtained with a block slider. For the infinite elastic slab the slip velocity exhibits sharp variations (shocks). On the contrary, a block slider does not exhibit this behavior. It is a clear example of the limitation of the use of such a simple analogy to describe the actual properties of the relative motions of two media in contact with friction

Improving Temporal Resolution in Ambient Noise Monitoring of Seismic Wave Speed

10p.International audienceThe use of ambient seismic noise has been intensively investigated to perform passive tomography at various scales. Besides passive tomography, passive monitoring is another application of seismic noise correlation as was shown by the recent observation of postseismic velocity changes around the San Andreas Fault in Parkfield, California. One of the drawbacks of using ambient noise correlation for passive monitoring is the need to average the correlations over a long time period in order to obtain a sufficient signal-to-noise ratio (SNR) for the phase fluctuations to be measured accurately. For the application to passive monitoring, one wants the possibility of following short-term velocity variations (one day or less) using noise correlation functions calculated on short time windows. Another difficulty may then appear when the spatial distribution of noise sources also evolves with time. The aim of this paper is to introduce an adaptive filter to the Parkfield dataset in order to improve the SNR output of the ambient noise correlation functions. When applied to passive monitoring, the temporal resolution can be increased from 30 days up to 1 day. With this improved temporal resolution, the velocity drop observed at Parkfield is shown to be cosesimic with the September 24, 2004 M_w=6.0 event. The relationship between the measured velocity fluctuations and the time-evolution of the spatial distribution of the noise wavefield is also investigated. Finally, the error bar in the amplitudes of the velocity variations are compared with a theoretical expectation