Lateral velocity variations in the near-surface reflect the presence of buried geological or anthropic structures, and their identification is of interest for many fields of application. Surface wave tomography (SWT) is a powerful technique for detecting both smooth and sharp lateral velocity variations at very different scales. A surface-wave inversion scheme derived from SWT is here applied to a 2-D active seismic dataset to characterize the shape of an urban waste deposit in an old landfill, located 15 km South of Vienna (Austria). First, the tomography-derived inverse problem for the 2-D case is defined: under the assumption of straight rays at the surface connecting sources and receivers, the forward problem for one frequency reduces to a linear relationship between observed phase differences at adjacent receivers and wavenumbers (from which phase velocities are straightforwardly derived). A norm damping regularization constraint is applied to ensure a smooth solution in space: the choice of the damping parameter is made through a minimization process, by which only phase variations of the order of the average wavelength are modelled. The inverse problem is solved for each frequency with a weighted least-squares approach, to take into account the data error variances. An independent multi-offset phase analysis (MOPA) is performed using the same dataset, for comparison: pseudo-sections from the tomography-derived linear inversion and MOPA are very consistent, with the former giving a more continuous result both in space and frequency and less artefacts. Local dispersion curves are finally depth inverted and a quasi-2-D shear wave velocity section is retrieved: we identify a well-defined low velocity zone and interpret it as the urban waste deposit body. Results are consistent with both electrical and electromagnetic measurements acquired on the same line.

For decades, bad practices in municipal and industrial waste management have had negative environmental impacts, generating high health risks for people and the environment. The use of badly designed, not engineered, and not well-operated landfills has, around the world, produced a large number of potentially contaminated sites, for which there are urgent needs to assess the actual risk and to proceed, in case, with reclamation activities. One of these sites, an abandoned waste disposal site located near a Site of Community Importance on the central-eastern coast of Sardinia (Italy), is the subject of the case history described in this work. As a part of a multi-method geophysical characterisation, a frequency-domain electromagnetic (FDEM) mapping survey was carried out with the specific aim of detecting the presence of buried materials (waste) and of delineating the lateral extent of the landfill by identifying the electrical conductivity anomalies produced, for the most part, by the conductive waste fill. Using an EM31 device in the vertical-dipole configuration, at a height of 0.9 m above the ground, both quadrature and in-phase electromagnetic responses were collected over a 7-hectare area with elevation varying between 6 m and 2.8 m above sea level. After removing the measurements identified as data coming from any recognisable surface man-made features within the survey area or near its perimeter, the filtered quadrature response (expressed as apparent conductivity) ranged from 5.5 mS/m to about 188.6 mS/m. All values are beyond the low induction number (LIN) condition and valid for the classical EM31 mapping, thus requiring advanced data processing. To obtain undistorted, meaningful, and interpretable high-resolution maps, measured data have been processed to correct the bias, introduced by the nonlinearity of the device, as a function of height above ground and the topography. The comparative analysis of the apparent conductivity map, obtained by the properly processed EM31 data and some aerial photos that clearly documented the site history, has allowed unequivocal delineation of the landfill extent, in good agreement with the results obtained with other geophysical methods (not described in this paper) and with the ground truthing data provided by three boreholes, which were core-drilled at the end of the study at three locations selected on the basis of the apparent conductivity map.

This work presents the results of an advanced geophysical characterization of a contaminated site, where a correct understanding of the dynamics in the unsaturated zone is fundamental to evaluate the effective management of the remediation strategies. Large-scale surface electrical resistivity tomography (ERT) was used to perform a preliminary assessment of the structure in a thick unsaturated zone and to detect the presence of a thin layer of clay supporting an overlying thin perched aquifer. Discontinuities in this clay layer have an enormous impact on the infiltration processes of both water and solutes, including contaminants. In the case here presented, the technical strategy is to interrupt the continuity of the clay layer upstream of the investigated site in order to prevent most of the subsurface water flow from reaching the contaminated area. Therefore, a deep trench was dug upstream of the site and, in order to evaluate the effectiveness of this approach in facilitating water infiltration into the underlying aquifer, a forced infiltration experiment was carried out and monitored using ERT and ground-penetrating radar (GPR) measurements in a cross-hole time-lapse configuration. The results of the forced infiltration experiment are presented here, with a particular emphasis on the contribution of hydro-geophysical methods to the general understanding of the subsurface water dynamics at this complex site.

Surface wave tomography (SWT) is a powerful and well-established technique to retrieve 3D shear-wave (S-wave) velocity models at the regional scale from earthquakes and seismic noise measurements. We have applied SWT to 3D active-source data, in which higher modes and heterogeneous spatial sampling make phase extraction challenging. First, synthetic traveltimes calculated on a dense, regular-spaced station array are used to test the performance of three different tomography algorithms (linearized inversion, Markov chain Monte Carlo [MCMC], and eikonal tomography). The tests suggest that the lowest misfit to the input model is achieved with the MCMC algorithm, at the cost of a much longer computational time. Then, real phases were extracted from a 3D exploration data set at different frequencies. This operation included an automated procedure to isolate the fundamental mode from higher order modes, phase unwrapping in two dimensions, and the estimation of the zero-offset phase. These phases are used to compute traveltimes between each source-receiver couple, which are input into the previously tested tomography algorithms. The resulting phase-velocity maps show good correspondence, highlighting the same geologic structures for all three methods. Finally, individual dispersion curves obtained by the superposition of phase-velocity maps at different frequencies are depth inverted to retrieve a 3D S-wave velocity model.

SUMMARY: In this paper, we present the results of an analysis of passive seismic noise recorded around the Scrovegni Chapel in Padua (Italy), using a dense 2-D network with nearly 1500 autonomous seismic nodes. Surface wave tomography using the active records allowed the imaging of several structures located at a depth of few metres, while this study focuses on the processing of about 22 hr of continuous passive records. First, the ambient noise is characterized in terms of amplitude, frequency content and azimuthal distribution, in order to ensure the applicability of the interferometric method. Second, a cross-correlation analysis is performed to retrieve virtual source gathers. Third, traveltimes are extracted from virtual source gathers using the same processing sequence applied to active gathers. Fourth, Eikonal tomography is run to retrieve isotropic phase velocity maps and azimuthal anisotropy. We compare and discuss the results obtained from the active and the passive methods, and finally propose a strategy for the integration of passive and active information. The new quasi-3-D shear wave velocity model obtained from the joint active and passive analysis is more accurate at depth, due to the addition of the passive low-frequency information.