During a quarter of a century, the main characteristics of the horizontal-to-vertical spectral ratio of ambient noise HVSRN have been extensively used for site effect assessment. In spite of the uncertainties about the optimum theoretical model to describe these observations, several schemes for inversion of the full HVSRN curve for near surface surveying have been developed over the last decade. In this work, a computer code for forward calculation of H/V spectra based on the diffuse field assumption (DFA) is presented and tested.It takes advantage of the recently stated connection between the HVSRN and the elastodynamic Green's function which arises from the ambient noise interferometry theory. The algorithm allows for (1) a natural calculation of the Green's functions imaginary parts by using suitable contour integrals in the complex wavenumber plane, and (2) separate calculation of the contributions of Rayleigh, Love, P-SV and SH waves as well. The stability of the algorithm at high frequencies is preserved by means of an adaptation of the Wang's orthonormalization method to the calculation of dispersion curves, surface-waves medium responses and contributions of body waves. This code has been combined with a variety of inversion methods to make up a powerful tool for passive seismic surveying.Comment: Published in Computers & Geosciences 97, 67-7

In order to evaluate the site effects on seismic ground motion and establish preventive measures to mitigate these effects, the dynamic characterization of sites is mandatory. Among the various geophysical tools aimed to this end, the horizontal to vertical spectral ratio (H/V) is a simple way to assess the dominant frequency of a site from seismic ambient noise. The aim of this communication is contributing to enhance the potential of this measurement with a novel method that allows extracting from the H/V the elastic properties of the subsoil, assumed here as a multilayer medium. For that purpose, we adopt the diffuse field assumption from both the experimental and the modelling perspectives. At the experimental end, the idea is to define general criteria that make the data processing closely supported by theory. On the modelling front, the challenge is to compute efficiently the imaginary part of Green's function. The Cauchy's residue theory in the horizontal wavenumber complex plane is the selected approach. This method allows both identifying the contributions of body and surface waves and computing them separately. This permits exploring the theoretical properties of the H/V under different compositions of the seismic ambient noise. This answers some questions that historically aroused and gives new insights into the H/V method. The efficient forward calculation is the prime ingredient of an inversion scheme based on both gradient and heuristic searches. The availability of efficient forward calculation of H/V allows exploring some relevant relationships between the H/V curves and the parameters. This allows generating useful criteria to speed up inversion. As in many inverse problems, the non-uniqueness issues also emerge here. A joint inversion method that considers also the dispersion curves of surface waves extracted from seismic ambient noise is presented and applied to experimental data. This joint scheme mitigates effectively the non-uniqueness.