Quasi Two-dimensional Transfer of Elastic Waves

A theory for multiple scattering of elastic waves is presented in a random medium bounded by two ideal free surfaces, whose horizontal size is infinite and whose transverse size is smaller than the mean free path of the waves. This geometry is relevant for seismic wave propagation in the Earth crust. We derive a time-dependent, quasi-2D radiative transfer equation, that describes the coupling of the eigenmodes of the layer (surface Rayleigh waves, SH waves, and Lamb waves). Expressions are found that relate the small-scale fluctuations to the life time of the modes and to their coupling rates. We discuss a diffusion approximation that simplifies the mathematics of this model significantly, and which should apply at large lapse times. Finally, coherent backscattering is studied within the quasi-2D radiative transfer equation for different source and detection configurations.Comment: REVTeX, 36 pages with 10 figures. Submitted to Phys. Rev.

Coda-Q in the 2.5s -20s period band from seismic noise - Application to the greater Alpine area

International audienceCoda-Q is used to estimate the attenuation and scattering properties of the Earth (Aki & Chouet 1975). So far focus has been on earthquake data at frequencies above 1 Hz, as the high noise level in the first and second microseismic peak, and possibly lower scattering coefficient, hinder stable measurements at lower frequencies. In this work, we measure and map coda-Q in the period bands 2.5 s-5 s, 5 s-10 s and 10 s-20 s in the greater Alpine region using noise cross-correlations between station pairs, based on data from permanent seismic stations and from the temporary AlpArray experiment. The observed coda-Q for short interstation distances is independent of azimuth so there is no indication of influence of the directivity of the incoming noise field on our measurements. In the 2.5 s-5 s and 5 s-10 s period bands, our measurements are self-consistent, and we observe stable geographic patterns of low and high coda-Q in the period bands 2.5 s-5 s and 5 s-10 s. In the period band 10 s-20 s, the dispersion of our measurements increases and geographic patterns become speculative. The coda-Q maps show that major features are observed with high resolution, with a very good geographical resolution of for example low coda-Q in the Po Plain. There is a sharp contrast between the Po Plain and the Alps and Apennines where coda-Q is high, with the exception a small area in the Swiss Alps which may be contaminated by the low coda-Q of the Po Plain. The coda of the correlations is too short to make independent measurements at different times within the coda, so we cannot distinguish between intrinsic and scattering Q. Measurements on more severely selected datasets and longer timeseries result in identical geographical patterns but lower numerical values. Therefore, high coda-Q values may be overestimated, but the geographic distribution between high and low coda-Q areas is respected. Our results demonstrate that noise correlations are a promising tool for extending coda-Q measurements to frequencies lower than those analysed with earthquake data

Phase statistics of seismic coda waves

We report the analysis of the statistics of the phase fluctuations in the coda of earthquakes recorded during a temporary experiment deployed at Pinyon Flats Observatory, California. The practical measurement of the phase is discussed and the main pitfalls are underlined. For large values, the experimental distributions of the phase first, second and third derivatives obey universal power-law decays whose exponents are remarkably well predicted by circular Gaussian statistics. For small values, these distributions are flat. The details of the transition between the plateau and the power-law behavior are governed by the wavelength. The correlation function of the first phase derivative along the array shows a simple algebro-exponential decay with the mean free path as the only length scale. Although only loose bounds are provided in this study, our work suggests a new method to estimate the degree of heterogeneity of the crComment: 4 figures, submitted to Physical Review Letter

Seismic energy envelopes in volcanic media : in need of boundary conditions

Peer reviewedPublisher PD

Energy partition of seismic coda waves in layered media: theory and application to Pinyon Flats Observatory

We have studied the partition of shear, compressional and kinetic energies in the coda of ten earthquakes recorded on a dense array, located at Pinyon Flats Observatory (PFO), California. We observe a clear stabilization of the shear to compressional ($W^s/W^p$) energy ratio in the coda, with an average value of about 2.8. The ratio between the vertical and horizontal kinetic energies ($V^2/H^2$) can be measured from 5 to 25Hz and shows an abrupt transitionfrom 0.1 in the 5-10Hz band, to about 0.8 in the 15-25Hz band. These measured values are in sharp contrast with the theoretical prediction for equipartitioned elastic waves in a homogeneous half-space. To explain these observations, we have developed a theory of equipartition in a layered elastic half-space. Using a rigorous spectral decomposition of the elastic wave equation, we define equipartition as a white noise distributed over the complete set of eigenfunctions. The theory predicts that close to the resonance frequency of a low-velocity layer, the ratio between shear and compressional energies strongly decreases. Using a detailed model of the subsurface at PFO, this conterintuitive result is found to be in good qualitative and quantitative agreement with the observations

Controlling Light Through Optical Disordered Media : Transmission Matrix Approach

We experimentally measure the monochromatic transmission matrix (TM) of an optical multiple scattering medium using a spatial light modulator together with a phase-shifting interferometry measurement method. The TM contains all information needed to shape the scattered output field at will or to detect an image through the medium. We confront theory and experiment for these applications and we study the effect of noise on the reconstruction method. We also extracted from the TM informations about the statistical properties of the medium and the light transport whitin it. In particular, we are able to isolate the contributions of the Memory Effect (ME) and measure its attenuation length

Autocorrelation of the Ground Vibrations Recorded by the SEIS-InSight Seismometer on Mars

Since early February 2019, the SEIS (Seismic Experiment for Interior Structure) seismometer deployed at the surface of Mars in the framework of the InSight mission has been continuously recording the ground motion at Elysium Planitia. In this study, we take advantage of this exceptional data set to put constraints on the crustal properties of Mars using seismic interferometry (SI). To carry out this task, we first examine the continuous records from the very broadband seismometer. Several deterministic sources of environmental noise are identified and specific preprocessing strategies are presented to mitigate their influence. Applying the principles of SI to the single-station configuration of InSight, we compute, for each Sol and each hour of the martian day, the diagonal elements of the time-domain correlation tensor of random ambient vibrations recorded by SEIS. A similar computation is performed on the diffuse waveforms generated by more than a hundred Marsquakes. A careful signal-to-noise ratio analysis and an inter-comparison between the two datasets suggest that the results from SI are most reliable in a narrow frequency band around 2.4 Hz, where an amplification of both ambient vibrations and seismic events is observed. The average autocorrelation functions (ACFs) contain well identifiable seismic arrivals, that are very consistent between the two datasets. Interpreting the vertical and horizontal ACFs as, respectively, the P- and S- seismic reflectivity below InSight, we propose a simple stratified velocity model of the crust, which is mostly compatible with previous results from receiver function analysis. Our results are discussed and compared to recent works from the literature.This study is InSight contribution number 164. The authors acknowledge both “Université Fédérale de Toulouse Midi Pyrénées” and the “Région Occitanie” for funding the PhD grant of Nicolas Compaire. The French authors acknowledge the French Space Agency CNES and ANR (ANR-14-CE36-0012-02 and ANR-19-CE31-0008-08) for funding the InSight Science analysis

The Polarization of Ambient Noise on Mars

Seismic noise recorded at the surface of Mars has been monitored since February 2019, using the InSight seismometers. This noise can reach −200 dB. It is 500 times lower than on Earth at night and it increases of 30 dB during the day. We analyze its polarization as a function of time and frequency in the band 0.03–1 Hz. We use the degree of polarization to extract signals with stable polarization independent of their amplitude and type of polarization. We detect polarized signals at all frequencies and all times. Glitches correspond to linear polarized signals which are more abundant during the night. For signals with elliptical polarization, the ellipse is in the horizontal plane below 0.3 Hz. In the 0.3-1Hz high frequency band (HF) and except in the evening, the ellipse is in the vertical plane and the major axis is tilted. While polarization azimuths are different in the two frequency bands, they both vary as a function of local hour and season. They are also correlated with wind direction, particularly during the daytime. We investigate possible aseismic and seismic origins of the polarized signals. Lander or tether noise can be discarded. Pressure fluctuations transported by wind may explain part of the HF polarization but not the tilt of the ellipse. This tilt can be obtained if the source is an acoustic emission coming from high altitude at critical angle. Finally, in the evening when the wind is low, the measured polarized signals may correspond to the seismic wavefield of the Mars background noise

Measurement of 19Ne spectroscopic properties via a new method of inelastic scattering to study novae

The accuracy of the predictions of the γ flux produced by a classical nova during the first hours after the outburst is limited by the uncertainties on several reaction rates, including the 18F(p,α)15O one. Better constraints on this reaction rate can be obtained by determining the spectroscopic properties of the compound nucleus 19Ne. This was achieved in a new inelastic scattering method using a 19Ne radioactive beam (produced by the GANIL-SPIRAL 1 facility) impinging onto a proton target. The experiment was performed at the VAMOS spectrometer. In this article the performances (excitation energy range covered and excitation energy resolution) and limitations of the new technique are discussed. Excitation energy resolution of σ = 33 keV and low background were obtained with this inverse kinematics method, which will allow extracting the spectroscopic properties of 19Ne