Optimized arrays for 2D cross‐borehole electrical tomography surveys

The use of optimized arrays generated using the ‘Compare R’ method for cross‐borehole resistivity measurements is examined in this paper. We compare the performances of two array optimization algorithms, one that maximizes the model resolution and another that minimizes the point spread value. Although both algorithms give similar results, the model resolution maximization algorithm is several times faster. A study of the point spread function plots for a cross‐borehole survey shows that the model resolution within the central zone surrounded by the borehole electrodes is much higher than near the bottom end of the boreholes. Tests with synthetic and experimental data show that the optimized arrays generated by the ‘Compare R’ method have significantly better resolution than a ‘standard’ measurement sequence used in previous surveys. The resolution of the optimized arrays is less if arrays with both current (or both potential) electrodes in the same borehole are excluded. However, they are still better than the ‘standard’ arrays.     

A numerical comparison of 2D resistivity imaging with 10 electrode arrays

Numerical simulations are used to compare the resolution and efficiency of 2D resistivity imaging surveys for 10 electrode arrays. The arrays analysed include pole‐pole (PP), pole‐dipole (PD), half‐Wenner (HW), Wenner‐α (WN), Schlumberger (SC), dipole‐dipole (DD), Wenner‐β (WB), γ‐array (GM), multiple or moving gradient array (GD) and midpoint‐potential‐referred measurement (MPR) arrays. Five synthetic geological models, simulating a buried channel, a narrow conductive dike, a narrow resistive dike, dipping blocks and covered waste ponds, were used to examine the surveying efficiency (anomaly effects, signal‐to‐noise ratios) and the imaging capabilities of these arrays. The responses to variations in the data density and noise sensitivities of these electrode configurations were also investigated using robust (L₁‐norm) inversion and smoothness‐constrained least‐squares (L₂‐norm) inversion for the five synthetic models. The results show the following. (i) GM and WN are less contaminated by noise than the other electrode arrays. (ii) The relative anomaly effects for the different arrays vary with the geological models. However, the relatively high anomaly effects of PP, GM and WB surveys do not always give a high‐resolution image. PD, DD and GD can yield better resolution images than GM, PP, WN and WB, although they are more susceptible to noise contamination. SC is also a strong candidate but is expected to give more edge effects. (iii) The imaging quality of these arrays is relatively robust with respect to reductions in the data density of a multi‐electrode layout within the tested ranges. (iv) The robust inversion generally gives better imaging results than the L₂‐norm inversion, especially with noisy data, except for the dipping block structure presented here. (v) GD and MPR are well suited to multichannel surveying and GD may produce images that are comparable to those obtained with DD and PD. Accordingly, the GD, PD, DD and SC arrays are strongly recommended for 2D resistivity imaging, where the final choice will be determined by the expected geology, the purpose of the survey and logistical considerations.     

Model Stacking (MOST) technique applied in cross-hole ERT field data for the detection of Thessaloniki ancient walls' depth

A cross-hole Electrical Resistivity Tomography (ERT) study was undertaken near the center of Thessaloniki in order to detect the depth of the existing city walls in the planned route of the new city underground train. This cross-hole setup was used for a study of measurements with various electrode arrays in real urban field conditions to evaluate the resolution of the models which is produced by each array and the reliability of the models which is produced by the newly published “MOST” technique. The pole–tripole array (C₂–C₁P₁P₂) produces high resolution models, even when only borehole electrodes are used. The bipole–bipole C₁C₂–P₁P₂ array, when used for cross-hole measurements only, produces higher resolution models compared to the C₁P₁–C₂P₂ array, even with a lower signal-to-noise ratio, which can result in extremely high RMS error, when noise, systematic or not, must be faced. The models of both arrays are greatly improved by the use of surface electrodes. The pole–bipole array (C₁–P₁P₂) is proved to be less accurate in imaging and quite unstable to the noisy urban environment and to systematic errors. Furthermore, the Model Stacking (MOST) interpretation technique leads to better results with models of greater resolution and fewer artifacts compared even with the combined data inversion. Finally, the ERT cross-hole analysis has been reliable in detecting the city walls.     

Imaging Pathways in Fractured Rock Using Three‐Dimensional Electrical Resistivity Tomography

Major challenges exist in delineating bedrock fracture zones because these cause abrupt changes in geological and hydrogeological properties over small distances. Borehole observations cannot sufficiently capture heterogeneity in these systems. Geophysical techniques offer the potential to image properties and processes in between boreholes. We used three‐dimensional cross borehole electrical resistivity tomography (ERT) in a 9 m (diameter) × 15 m well field to capture high‐resolution flow and transport processes in a fractured mudstone contaminated by chlorinated solvents, primarily trichloroethylene. Conductive (sodium bromide) and resistive (deionized water) injections were monitored in seven boreholes. Electrode arrays with isolation packers and fluid sampling ports were designed to enable acquisition of ERT measurements during pulsed tracer injections. Fracture zone locations and hydraulic pathways inferred from hydraulic head drawdown data were compared with electrical conductivity distributions from ERT measurements. Static ERT imaging has limited resolution to decipher individual fractures; however, these images showed alternating conductive and resistive zones, consistent with alternating laminated and massive mudstone units at the site. Tracer evolution and migration was clearly revealed in time‐lapse ERT images and supported by in situ borehole vertical apparent conductivity profiles collected during the pulsed tracer test. While water samples provided important local information at the extraction borehole, ERT delineated tracer migration over spatial scales capturing the primary hydrogeological heterogeneity controlling flow and transport. The fate of these tracer injections at this scale could not have been quantified using borehole logging and/or borehole sampling methods alone.     

Constructing the Apparent Geological Model by Fusing Surface Resistivity Survey and Borehole Records

We constructed an apparent geological model with resistivity data from surface resistivity surveys. We developed a data fusion approach by integrating dense electrical resistivity measurements collected with Schlumberger arrays and wellbore logs. This approach includes an optimization algorithm and a geostatistic interpolation method. We first generated an apparent formation factor model from the surface resistivity measurements and groundwater resistivity records with an inverse distance method. We then converted the model into a geology model with the optimized judgment criteria from the algorithms relating the apparent formation factors to the borehole geology. We also employed a non-parametric bootstrap method to analyze the uncertainty of the predicted sediment types, and the predictive uncertainties of clay, gravel, and sand were less than 5%. Overall, our model is capable of capturing the spatial features of the sediment types. More importantly, this approach can be arranged in a self-updated sequence to enable adjustments to the model to accommodate newly collected core records or geophysical data. This approach yields a more detailed apparent geological model for use in future groundwater simulations, which is of benefit to multi-discipline studies.     

Numerical sensitivity test of three-electrode laterolog borehole tool

Finite element numerical simulation has been carried out to investigate quantitatively the response of the three-electrode laterolog borehole tool (LL3) on radial and vertical heterogeneity of the rock. In order to calculate the apparent resistivity from the electric potential and the current discharge of the measurement electrode the probe coefficient of the LL3 tool with finite electrode extent was determined. Two independent methods, a finite element modeling and a semi-analytical solution, resulted in the probe coefficient of approx. 0.15 m with a relative deviation of 2.4% due to the different geometry, resolution and electronics of the models. It was established that LL3 is only slightly sensitive to the presence of mud when the borehole diameter is d ≤ 30 cm and the ratio of the resistivity of rock and the borehole mud is 1 ≤ R _t/R _m ≤ 1000. Vertical heterogeneity test pointed out that the layer boundaries can be localized exactly even for thin bedded layer (with a thickness of 1 m) and the presence of low-resistive borehole mud. Correction factors were suggested to decrease the biasing effect of the low-resistive borehole mud and the shoulder beds on the apparent resistivity observed by LL3. Finally, it was verified that the probe has large penetration depth with excellent vertical resolution, what explains the enduring popularity of the LL3 tool in well logging.     

Generation of a pseudo-2D shear-wave velocity section by inversion of a series of 1D dispersion curves

Multichannel Analysis of Surface Waves utilizes a multichannel recording system to estimate near-surface shear (S)-wave velocities from high-frequency Rayleigh waves. A pseudo-2D S-wave velocity (v_S) section is constructed by aligning 1D models at the midpoint of each receiver spread and using a spatial interpolation scheme. The horizontal resolution of the section is therefore most influenced by the receiver spread length and the source interval. The receiver spread length sets the theoretical lower limit and any v_S structure with its lateral dimension smaller than this length will not be properly resolved in the final v_S section. A source interval smaller than the spread length will not improve the horizontal resolution because spatial smearing has already been introduced by the receiver spread.In this paper, we first analyze the horizontal resolution of a pair of synthetic traces. Resolution analysis shows that (1) a pair of traces with a smaller receiver spacing achieves higher horizontal resolution of inverted S-wave velocities but results in a larger relative error; (2) the relative error of the phase velocity at a high frequency is smaller than at a low frequency; and (3) a relative error of the inverted S-wave velocity is affected by the signal-to-noise ratio of data. These results provide us with a guideline to balance the trade-off between receiver spacing (horizontal resolution) and accuracy of the inverted S-wave velocity. We then present a scheme to generate a pseudo-2D S-wave velocity section with high horizontal resolution using multichannel records by inverting high-frequency surface-wave dispersion curves calculated through cross-correlation combined with a phase-shift scanning method. This method chooses only a pair of consecutive traces within a shot gather to calculate a dispersion curve. We finally invert surface-wave dispersion curves of synthetic and real-world data. Inversion results of both synthetic and real-world data demonstrate that inverting high-frequency surface-wave dispersion curves – by a pair of traces through cross-correlation with phase-shift scanning method and with the damped least-square method and the singular-value decomposition technique – can feasibly achieve a reliable pseudo-2D S-wave velocity section with relatively high horizontal resolution.     

Correlation of vertical electrical sounding and electrical borehole log data for groundwater exploration

We have correlated the longitudinal unit conductance CL obtained from interpreted vertical electrical sounding data with the formation resistivity R_t and the formation resistivity factor F, obtained by carrying out electrical borehole logging. Interpreted geophysical data of eleven soundings and two electrical borehole log records are used for the analysis. The geophysical data used were acquired in a sedimentary basin. The study area is called Lower Maner Basin located in the province of Andhra Pradesh, India. Vertical electrical soundings were carried out using a Schlumberger configuration with half current electrode separation varying from 600–1000 m. For logging the two boreholes, a Widco logger‐model 3200 PLS was used. True formation resistivity R_t was calculated from a resistivity log. Formation resistivity factor F was also calculated at various depths using R_t values. An appreciable inverse relation exists between the correlated parameters. The borehole resistivity R_t and the formation resistivity factor F decrease with the increase in the longitudinal unit conductance CL. We have shown the use of such a relation in computing borehole resistivity R_t and formation resistivity factor F at sites that posses only vertical electrical sounding data, with a fair degree of accuracy. Validation of the correlation is satisfactory. Scope for updating the correlation is discussed. Significance and applications of the relation for exploration of groundwater, namely to update the vertical electrical sounding data interpretation by translating the vertical electrical sounding data into electrical borehole log parameters, to facilitate correlations studies and to estimate the porosity (φ), permeability (K) and water saturation S_w of water bearing zones are discussed.     

InterPACIFIC project: Comparison of invasive and non-invasive methods for seismic site characterization. Part II: Inter-comparison between surface-wave and borehole methods

The InterPACIFIC project was aimed at assessing the reliability, resolution, and variability of geophysical methods in estimating the shear-wave velocity profile for seismic ground response analyses. Three different subsoil conditions, which can be broadly defined as soft-soil, stiff-soil, and hard-rock, were investigated. At each site, several participants performed and interpreted invasive measurements of shear wave velocity (Vs) and compression wave velocity (Vp) in the same boreholes. Additionally, participants in the project analysed a common surface-wave dataset using their preferred strategies for processing and inversion to obtain Vs profiles. The most significant difference between the invasive borehole methods and non-invasive surface wave methods is related to resolution of thin layers and abrupt contrasts, which is inherently better for invasive methods. However, similar variability is observed in the estimated invasive and non-invasive Vs profiles, underscoring the need to account for such uncertainty in site response studies. V_S,30 estimates are comparable between invasive and non-invasive methods, confirming that the higher resolution provided by invasive methods is quite irrelevant for computing this parameter.     

Soil characterization using electrical resistivity tomography and geotechnical investigations

Electrical Resistivity Tomography (ERT) has been used in association with Standard Penetration Test (SPT) and Dynamic Cone Penetration Test (DCPT) for Geotechnical investigations at two sites, proposed for thermal power plants, in Uttar Pradesh (UP), India. SPT and DCPT tests were conducted at 28 points and two ERT profiles, each measuring 355 m long, were recorded using 72 electrodes deployed at 5 m spacing. Electrical characterization of subsurface soil was done using borehole data and grain size analysis of the soil samples collected from boreholes. The concept of electrical resistivity variation with soil strength related to the grain size distribution, cementation, porosity and saturation has been used to correlate the transverse resistance of soil with the number of blow counts (N-values) obtained from SPT and DCPT data. It was thus observed that the transverse resistance of soil column is linearly related with the number of blow counts (N-values) at these sites. The linear relationships are site-specific and the coefficients of linear relation are sensitive to the lithology of subsurface formation, which was verified by borehole data. The study demonstrates the usefulness of the ERT method in geotechnical investigations, which is economic, efficient and less time consuming in comparison to the other geotechnical methods, such as SPT and DCPT, used for the purpose.     

基于HJ-CCD和MODIS的吉林省中西部湖泊透明度反演对比

水体透明度能够反映光在水体中的穿透程度,影响水生植被及以光为依赖条件的水生生物的分布,获取透明度的传统方法是采用透明度盘进行测量,但也可以通过遥感方法获得.环境减灾卫星是专门用于环境与灾害监测预报的小卫星星座,影像覆盖范围广,空间、时相分辨率较高,可以为水环境遥感提供较好的数据源.MODIS数据在近岸水体和内陆大型湖泊水环境监测中也有广泛应用,它的时相分辨率也很高,但空间分辨率低.利用HJ-1A卫星CCD数据和MODIS日反射率产品(MOD09GA),以2012年9月吉林省石头口门水库、二龙湖、查干湖、月亮泡等地的实测透明度为基础(实测点数74个,最小值为0.134 m,最大值为1.410 m,平均值为0.488 m),根据灰色关联度选取构建模型的波段组合,建立水体透明度反演模型.HJ1A-CCD数据与MOD09GA数据建立的模型R2分别为0.639和0.894,均方根误差(RMSE)分别为0.248和0.135,模型验证的平均相对误差(MRE)分别为17.1%和9.5%,RMSE分别为0.207和0.089.MODIS数据以其较高的辐射分辨率使模型精度较高,但是HJ数据在应用于透明度小于1 m的水体时精度也较高(MRE=13.5%,RMSE=0.066).HJ-CCD数据在空间分辨率上的优势使其能够获得透明度空间分布的细节信息.比较两者反演得到的湖泊平均透明度,结果较为一致.     

Predicting spatio-temporal concentrations of PM<Subscript>2.5</Subscript> using land use and meteorological data in Yangtze River Delta,China

The prediction of PM_2.5 concentrations with high spatiotemporal resolution has been suggested as a potential method for data collection to assess the health effects of exposure. This work predicted the weekly average PM_2.5 concentrations in the Yangtze River Delta, China, by using a spatio-temporal model. Integrating land use data, including the areas of cultivated land, construction land, and forest land, and meteorological data, including precipitation, air pressure, relative humidity, temperature, and wind speed, we used the model to estimate the weekly average PM_2.5 concentrations. We validated the estimated effects by using the cross-validated R² and Root mean square error (RMSE); the results showed that the model performed well in capturing the spatiotemporal variability of PM_2.5 concentration, with a reasonably large R² of 0.86 and a small RMSE of 8.15 (μg/m³). In addition, the predicted values covered 94% of the observed data at the 95% confidence interval. This work provided a dataset of PM_2.5 concentration predictions with a spatiotemporal resolution of 3 km × week, which would contribute to accurately assessing the potential health effects of air pollution.     

Tomographic inversion of near-surface <Emphasis Type="Italic">Q</Emphasis> factor by combining surface and cross-hole seismic surveys

The estimation of the quality factor Q plays a fundamental role in enhancing seismic resolution via absorption compensation in the near-surface layer. We present a new geometry that can be used to acquire field data by combining surface and cross-hole surveys to decrease the effect of geophone coupling on Q estimation. In this study, we drilled number of receiver holes around the source hole, each hole has different depth and each geophone is placed geophones into the bottom of each receiver hole to avoid the effect of geophone coupling with the borehole wall on Q estimation in conventional cross-hole seismic surveys. We also propose a novel tomographic inversion of the Q factor without the effect of the source signature, and examine its stability and reliability using synthetic data. We estimate the Q factors of the near-surface layer in two different frequency bands using field data acquired in the Dagang Oilfield. The results show that seismic absorption in the near-surface layer is much greater than that in the subsurface strata. Thus, it is of critical practical importance to enhance the seismic solution by compensating for near-surface absorption. In addition, we derive different Q factors from two frequency bands, which can be treated, to some extent, as evidence of a frequency-dependent Q.     

EGS hydraulic stimulation monitoring by surface arrays - location accuracy and completeness magnitude: the Basel Deep Heat Mining Project case study

The potential and limits of monitoring induced seismicity by surface-based mini arrays was evaluated for the hydraulic stimulation of the Basel Deep Heat Mining Project. This project aimed at the exploitation of geothermal heat from a depth of about 4,630?m. As reference for our results, a network of borehole stations by Geothermal Explorers Ltd. provided ground truth information. We utilized array processing, sonogram event detection and outlier-resistant, graphical jackknife location procedures to compensate for the decrease in signal-to-noise ratio at the surface. We could correctly resolve the NNW?CSSE striking fault plane by relative master event locations. Statistical analysis of our catalog data resulted in M _L 0.36 as completeness magnitude, but with significant day-to-night dependency. To compare to the performance of borehole data with M _W 0.9 as completeness magnitude, we applied two methods for converting M _L to M _W which raised our M _C to M _W in the range of 0.99?C1.13. Further, the b value for the duration of our measurement was calculated to 1.14 (related to M _L), respectively 1.66 (related to M _W), but changes over time could not be resolved from the error bars.     

Characterization of seasonal local recharge using electrical resistivity tomography and magnetic resonance sounding

A groundwater recharge process of heterogeneous hard rock aquifer in the Moole Hole experimental watershed, south India, is being studied to understand the groundwater flow behaviour. Significant seasonal variations in groundwater level are observed in boreholes located at the outlet area indicating that the recharge process is probably taking place below intermittent streams. In order to localize groundwater recharge zones and to optimize implementation of boreholes, a geophysical survey was carried out during and after the 2004 monsoon across the outlet zone. Magnetic resonance soundings (MRS) have been performed to characterize the aquifer and measure groundwater level depletion. The results of MRS are consistent with the observation in boreholes, but it suffers from degraded lateral resolution. A better resolution of the regolith/bedrock interface is achieved using electrical resistivity tomography (ERT). ERT results are confirmed by resistivity logging in the boreholes. ERT surveys have been carried out twice—before and during the monsoon—across the stream area. The major feature of recharge is revealed below the stream with a decrease by 80% of the calculated resistivity. The time‐lapse ERT also shows unexpected variations at a depth of 20 m below the slopes that could have been interpreted as a consequence of a deep seasonal water flow. However, in this area time‐lapse ERT does not match with borehole data. Numerical modelling shows that in the presence of a shallow water infiltration, an inversion artefact may take place thus limiting the reliability of time‐lapse ERT. A combination of ERT with MRS provides valuable information on structure and aquifer properties respectively, giving a clue for a conceptual model of the recharge process: infiltration takes place in the conductive fractured‐fissured part of the bedrock underlying the stream and clayey material present on both sides slows down its lateral dissipation. Copyright © 2007 John Wiley & Sons, Ltd.     

Super-resolution least-squares prestack Kirchhoff depth migration using the <Emphasis Type="Italic">L</Emphasis><Subscript>0</Subscript>-norm

Least-squares migration (LSM) is applied to image subsurface structures and lithology by minimizing the objective function of the observed seismic and reverse-time migration residual data of various underground reflectivity models. LSM reduces the migration artifacts, enhances the spatial resolution of the migrated images, and yields a more accurate subsurface reflectivity distribution than that of standard migration. The introduction of regularization constraints effectively improves the stability of the least-squares offset. The commonly used regularization terms are based on the L₂-norm, which smooths the migration results, e.g., by smearing the reflectivities, while providing stability. However, in exploration geophysics, reflection structures based on velocity and density are generally observed to be discontinuous in depth, illustrating sparse reflectance. To obtain a sparse migration profile, we propose the super-resolution least-squares Kirchhoff prestack depth migration by solving the L₀-norm-constrained optimization problem. Additionally, we introduce a two-stage iterative soft and hard thresholding algorithm to retrieve the super-resolution reflectivity distribution. Further, the proposed algorithm is applied to complex synthetic data. Furthermore, the sensitivity of the proposed algorithm to noise and the dominant frequency of the source wavelet was evaluated. Finally, we conclude that the proposed method improves the spatial resolution and achieves impulse-like reflectivity distribution and can be applied to structural interpretations and complex subsurface imaging.     

Tilt observations from shallow and deep borehole arrays

Measurements of the Earth tide using tilt and strainmeter data can be used to infer crustal structure. However, in almost all measurements of crustal motion using shallow tiltmeters and strainmeters, reference is made to the degradation of the data implied to be caused by influence of meteorological parameters on near surface layers. For example, studies by Wood and King and Goulty et al. indicated that a substantial portion of measured crustal motion by shallow tilt and strain data may be attributed to the weather. These studies were limited to data outside the tidal band. However, it is important to determine the effects of weather on these data, within the tidal band, if shallow tiltmeters are to be used to determine crustal structure.Collocated deep and shallow borehole tiltmeter arrays were operated by The Air Force Geophysics Laboratory in eastern Massachusetts and provided an opportunity to investigate weather effects on shallow tilt. The tiltmeters installed in deep boreholes (100 m) yielded internally consistent M₂ amplitudes (2%) and phases (3%) (Cabaniss and McConnell) but the collocated shallow (4 m) borehole tilt array did not produce the same results. Least squares analyses of the shallow tilt data revealed variations of 20% for M₂ tidal amplitude compared to the deep data. A multichannel Wiener filter technique was used to investigate the hypothesis that this discrepancy was due to meteorological influence. Results indicated that approximately one-half of the anomalous M₂ power was correlated with air and borehole temperatures. It is suggested that the remaining anomalous power may be due to ground deformation caused by the action of groundwater. If this is true, shallow tiltmeters might not be measuring the same deformation as the deep instruments and, therefore, not capable of measuring the Earth tide with the accuracy required to determine crustal structure.     

Compiling and verifying 3D models of 2D induced polarization and resistivity data by geostatistical methods

A 3D model of collected time-domain induced polarization (IP) and electrical resistivity tomography (ERT) data is compiled by geostatistical methods as well as studying spatial correlation among the database. Mesgaran copper deposit, located in Birjand eastern Iran, was chosen to compile and verify the model, leading to five parallel surveyed IP and ERT profiles with dipole–dipole arrays. The collected data were inverted, and then 2D models of IP and ER were prepared; also 3D inversion was done. Afterward, the 3D model has been built by geostatistical methods. Correspondingly, the anomalies threshold was detected by fractal methods and the estimation variance and Kriging efficiency were calculated to validate the modeling. The mineralization zones were determined according to the classified anomalies and those with the lowest error. Results indicated a high correlation between anomalies identified from the model and mineralization. The results made it possible to construct 3D models from surveyed 2D data with acceptable error level.     

A surface seismic approach to liquefaction

The liquefaction potential of soils is traditionally assessed through geotechnical approaches based on the calculation of the cyclical stress ratio (CSR) induced by the expected earthquake and the ‘resistance’ provided by the soil, which is quantified through standard penetration (SPT), cone penetration (CPT), or similar tests. In more recent years, attempts to assess the liquefaction potential have also been made through measurement of shear wave velocity (V_S) in boreholes or from the surface. The latter approach has the advantage of being non-invasive and low cost and of surveying lines rather than single points. However, the resolution of seismic surface techniques is lower than that of borehole techniques and it is still debated whether it is sufficient to assess the liquefaction potential.In this paper we focus our attention on surface seismic techniques (specifically the popular passive and active seismic techniques based on the correlation of surface waves such as ReMi^TM, MASW, ESAC, SSAP, etc.) and explore their performance in assessing the liquefaction susceptibility of soils. The experimental dataset is provided by the two main seismic events of M_L=5.9 and 5.8 (M_W=6.1, M_W=6.0) that struck the Emilia-Romagna region (Northern Italy) on May 20 and 29, 2012, after which extensive liquefaction phenomena were documented in an area of 1200 km².The CPT and drillings available in the area allow us to classify the soils into four classes: A) shallow liquefied sandy soils, B) shallow non-liquefied sandy soils, C) deep non-liquefied sandy soils, and D) clayey–silty soils, and to determine that on average class A soils presented a higher sand content at the depth of 5–8 m compared to class B soils, where sand was dominant in the upper 5 m. Surface wave active–passive surveys were performed at 84 sites, and it was found that they were capable of discriminating among only three soil classes, since class A and B soils showed exactly the same V_S distribution, and it is possible to show both experimentally and theoretically that they appear not to have sufficient resolution to address the seismic liquefaction issue.As a last step, we applied the state-of-the art CSR–V_S method to assess the liquefaction potential of sandy deposits and we found that it failed in the studied area. This might be due to the insufficient resolution of the surface wave methods in assessing the Vs of thin layers and to the fact that Vs scales with the square root of the shear modulus, which implies an intrinsic lower sensitivity of Vs to the shear resistance of the soil compared to parameters traditionally measured with the penetration tests. However, it also emerged that the pure observation of the surface wave dispersion curves at their simplest level (i.e. in the frequency domain, with no inversion) is still potentially informative and can be used to identify the sites where more detailed surveys to assess the liquefaction potential are recommended.