Analysis of the Seismic Properties for Engineering Purposes of the Shallow Subsurface: Two Case Studies from Italy and Croatia

: We present two case studies of the application of seismic surveys to estimate the elastic properties of soil and rock in the shallow subsurface. The two sites present very different geological characteristics. The ﬁrst test site is a town on the Croatian coast, not far from the city of Split, built on hard rock, where we acquired three seismic lines. The second site is located in the outskirts of the city of Ferrara, in Italy, in an alluvial plain, where two lines were acquired. In both sites, for detailed characterization, we acquired surface-, compressional- and shear-waves, further distinguishing the latter between horizontally (SH) and vertically (SV) polarized waveﬁelds. We processed the data by performing a Multichannel Analysis of Surface Waves to compute a preliminary one-dimensional shear wave velocity proﬁle. Then, we performed ﬁrst-break tomography to compute P-, SH- and SV-velocity proﬁles. Such unusual acquisition allowed us to compute not only basic engineering parameters such as the equivalent shear-wave velocity of the ﬁrst 30 m of subsurface (V S30 ) from the SH proﬁles but also other useful parameters such as the V P /V S and estimate the anisotropy of the medium thanks to the V SV /V SH . Given the level of detail of the results and their engineering value, we conclude that the method of investigation we applied in the two test sites is a valuable tool for characterizing the shallow subsurface.


Introduction
Geotechnical characterization of the shallow geological structures of a site is an essential tool to evaluate the response of the terrain to a macroseismic event. Typically, the geotechnical properties can be estimated by analyzing samples in a laboratory [1] or with in situ measurements such as the core penetration test and the load-bearing test.
However, these have the main disadvantage that they only provide a punctual estimation of the parameters. For this reason, geophysical surveys are becoming increasingly popular in seismic engineering thanks to their ability to accurately compute the main geotechnical parameters of the terrain and therefore its ability to amplify the seismic waves over a relatively wide area [2]. Specifically, the seismic method is able to evaluate the seismic velocities of the medium, which are a proxy for its elastic moduli [3]. More precisely, the seismic velocities are related to the bulk and shear moduli of the terrain, which describe the stiffness of the medium and therefore its ability to amplify seismic waves [4]. Based on this, in the past several decades, seismologists have established the vital importance of near-surface shear-wave velocity characterization both from a theoretical [5] and observational point of view [6]. Legislation followed these studies, and currently, the most commonly used geotechnical parameter is the equivalent shear-wave velocity of the first 30 m of subsoil (V S30 ) [7][8][9]. The most widely used geophysical method to compute V S30 is the Multichannel Analysis of Surface Waves (MASW) [10]. Technically, MASW relies on the dispersive behavior of the surface waves and consists in extracting a dispersion curve from a recorded seismogram and inverting it, obtaining a 1D shear-wave velocity profile. Yust et al. (2018) [11] show that the method is able to estimate accurately V s profiles and therefore compute the V S30 with a relatively small error. Another method that relies on the surface waves is the Horizontal to Vertical Spectral Ratio (HVSR, [12]). The method consists in processing the data from 3C single-station data by computing the ratio between the horizontal and vertical spectra and thereby extracting a 1D shear-wave velocity profile together with the resonant frequency of the terrain.
A method that is gaining popularity is the travel-time tomography of first breaks on an SH active seismic survey [13,14]. The main advantage of this method is that it provides a 2D (potentially even 3D) S-wave velocity profile and therefore a computation of the V S30 along the entire seismic line.
A more accurate characterization of the shallow subsurface, beyond the computation of the V S30 , can be of high interest for engineers. For instance, from the computation of Poisson's ratio or V P /V S , information regarding the fracturing of the rocks [1] and its fluid saturation [15] can be inferred.
In this work, we present the application of an innovative survey technique in two case studies from two very different sites. The first was recorded along the roads of a small town lying on hard rock in the area of Split in Croatia. The second case study consists of a dataset recorded in an alluvial plain in an intensively cultivated countryside close to the city of Ferrara. The novelty of the method proposed lies in the fact that we recorded separately compressional-and shear-wavefields, distinguishing in the latter between the components orthogonal (SH) and parallel (SV) to the seismic line. In fact, while integrating P-and SH-waves is increasingly popular in the near-surface geophysics community [16,17], the distinction between SH-and SV-wavefields is much rarer. Seismologists performed several numerical simulations [18] and analyzed field data [19] aiming at the study of shear-wave splitting and anisotropy, but hardly any field example using controlled sources can be found in the literature.
In the same surveys, data specifically for surface-wave analysis was also recorded. The final aim of our analyses is to provide detailed engineering information of the shallow subsurface without invasive tests at a relatively large scale. Successful testing of the method in two such different locations should allow us to confirm whether the method is replicable in most situations.

Geological Context
The Croatian site lies in the village of Kaštela Kambelovac, which is a town lying on the shores of the homonymous bay, just northwest of the city of Split. Geologically, the Bay lies in a compressional environment [20], where the rocks are subject to strong thrusting and folding. This caused the thrusting of the Senonian limestone on top of the Eocene flysch [21,22]. The town therefore is built on top of such flysch, the thickness of which has been estimated to be several hundreds of meters and the composition of which is known to be mainly of alternating layers of marls and calcarenites/calcirudites [23]. Previous geological and geotechnical investigations in the area showed that the bedding of the flysch is subvertical because of the tectonic compressional context. More detailed geological information about the Croatian site can be found in [24].
As for the Italian site, the acquisition took place in the southeastern suburbs of the city of Ferrara, in the central Po Valley, an alluvial, subsiding plain. Specifically, the area where the studies were carried out lies at the heart of such a plain, filled with olocenic sediments of marine, lagunal and riverine origin [25]. These are horizontally layered and composed of sand mixed with clay and silt, the proportion of which is variable both depending on the location and on the depth. These sediments are often not consolidated, giving place to phenomena of liquefaction during macroseismic events [26]. At the base of the olocenic sediments, located at approximately 1000 m depth, lie the so-called "Ferrara folds", pre-pleistocene rocks folded by the Appennine orogeny ("Buried Appennine" [27]). In fact, tectonically, the area lies at the margin of the compressional environment that created the Appennine mountain range.

Data Acquisition
In both sites, we acquired three wavefields separately (P, SH and SV). To do this, we used as source a wheelbarrow-mounted vibrator, which is capable of generating both P-and S-waves, reorienting it orthogonal and parallel to the line to generate SH-and SV-waves, respectively. As for the receivers, we used 10 Hz vertical geophones to record the P-waves and 14 Hz horizontal geophones to record the S-waves. The latter were re-oriented to record SH and SV, respectively, orthogonal and parallel to the survey line. We vibrated twice at each shot point, to perform stacking of the seismograms and therefore increase the signal-noise ratio.

Ferrara Site
We acquired two high-resolution seismic lines, the location of which can be seen in Figure 1a, while the UTM coordinates of the two lines are reported in Table 1. For both lines, we deployed 150 active channels, spaced every 2 m, and we shot every 4 m. The length of both lines is therefore approximately 300 m.
Appl. Sci. 2022, 12, 4535 3 of 12 of the olocenic sediments, located at approximately 1000 m depth, lie the so-called "Ferrara folds", pre-pleistocene rocks folded by the Appennine orogeny ("Buried Appennine" [27]). In fact, tectonically, the area lies at the margin of the compressional environment that created the Appennine mountain range.

Data Acquisition
In both sites, we acquired three wavefields separately (P, SH and SV). To do this, we used as source a wheelbarrow-mounted vibrator, which is capable of generating both Pand S-waves, reorienting it orthogonal and parallel to the line to generate SH-and SVwaves, respectively. As for the receivers, we used 10 Hz vertical geophones to record the P-waves and 14 Hz horizontal geophones to record the S-waves. The latter were re-oriented to record SH and SV, respectively, orthogonal and parallel to the survey line. We vibrated twice at each shot point, to perform stacking of the seismograms and therefore increase the signal-noise ratio.

Ferrara Site
We acquired two high-resolution seismic lines, the location of which can be seen in Figure 1a, while the UTM coordinates of the two lines are reported in Table 1. For both lines, we deployed 150 active channels, spaced every 2 m, and we shot every 4 m. The length of both lines is therefore approximately 300 m.   Along Line 1 (highlighted in red in Figure 1a), we acquired surface-wave data by deploying 48 4.5 Hz receivers, which were spaced every 4 m. As a source, we used a 100 kg weight-drop. We shot four times: twice at the beginning and twice at the end of the line.
The data of Line 1 were acquired along a gravel road. The main sources of noise were the agricultural activities ongoing in the nearby fruit farms (tractors, walking, irrigation pumps) and the car traffic.
The data of Line 2 were acquired in a fruit farm along a line of pear trees. Therefore, the geophones were planted in the soft soil. The main sources of noise are, similarly to Line 1, the activities going on in the fruit farm, including tractors and irrigation pumps.

Kaštela Site
We acquired three seismic lines along the roads of the village, planting the sensors in holes drilled in the tarmac. The details of the survey are outlined in [14], and the location of the lines can be seen in Figure 1b, while the UTM coordinates are reported in Table 1. Line 1 and 2 are 300 m long, while Line 3 is 150 m long. The receiver spacing in all lines is 2 m, and the shot spacing is 4 m. It is to be noted that during the acquisition of Line 1 and Line 2, traffic was quite intense in the village, and the weather was windy and even rainy at times. On the other hand, Line 3 was less affected by traffic noise, and the weather was fine.

Processing Methods
Multichannel Analysis of Surface Waves [10] provided information only on the Ferrara site, as no dispersive event could be identified in the Kaštela dataset, which was probably due the presence of tubes or other voids just below the road.
In Figure 2a, we show the frequency-wavenumber (f-k) spectrum of the Ferrara dataset, from which we picked the most energetic linear (i.e., dispersive) event to compute the dispersion curve shown in Figure 2b. The frequency range is from 3.5 to 15 Hz, giving a maximum wavelength of 58 m, i.e., a penetration depth of 29 m (λ/2, where λ is the wavelength). Finally, we inverted such a dispersion curve using the neighborhood algorithm in the GeoPsy software [27]. The initial model consists of three layers + Halfspace, the properties of which are described in Table 2. One million profiles were computed; in Figure 2c, we show the 500 best-that is, those generating a synthetic dispersion curve having the lowest misfit with respect to the picked one.  We performed first break tomography on all the lines for all of the wavefields. Before this, some pre-processing was performed. Specifically, to increase the signal-to-noise ratio, for each shot, we applied a predictive deconvolution to the uncorrelated signals as in [28]. Furthermore, to remove random noise, we stacked in the time domain the seismograms relative to the same shot point.
The first breaks were manually picked on the seismograms. Examples of first-break picking can be seen in Figure 3.  We performed first break tomography on all the lines for all of the wavefields. Before this, some pre-processing was performed. Specifically, to increase the signal-to-noise ratio, for each shot, we applied a predictive deconvolution to the uncorrelated signals as in [28]. Furthermore, to remove random noise, we stacked in the time domain the seismograms relative to the same shot point.
The first breaks were manually picked on the seismograms. Examples of first-break picking can be seen in Figure 3. Then, the travel times were inverted with Cat3D software [29], which is a tomographic package developed at the National Institute of Oceanography and Applied Geophysics-OGS that uses a ray-tracing algorithm based on a minimum time principle and SIRT (Simultaneous Iterative Reconstruction Technique) as an inversion method [30].

Results and Discussion
In Figure 4, we show the velocity profiles from the tomographic inversion of all wavefields relative to the three lines acquired in Kaštela, while in Figure 5, we show those relative to the two lines acquired in Ferrara. The two velocity profiles show very different features. In Kaštela, below a slower weathered layer, the velocity rapidly increases to 4.5 km/s and 3.0 km/s for P-and S-waves, respectively. Furthermore, we observe an overall increase in the velocities in the eastern lines (Lines 2 and 3) compared to the western line 1. In Figure 4d-f, we show the SH-velocity profiles and the corresponding VS30. Since the velocity is >800 m/s along all three profiles, we conclude that the soil can be classified as A-class following the Eurocode-8 provisions [9]. In the bottom-right corner of each velocity profile, the mean RMS error is reported. All profiles present a mean RMS lower than 5%, except for the P-wave profile of line 2. This is probably due to the low signal-to-noise ratio of the data caused by the rainfall during the acquisition. Having said this, all other profiles show low mean RMS and can therefore be considered reliable. Then, the travel times were inverted with Cat3D software [29], which is a tomographic package developed at the National Institute of Oceanography and Applied Geophysics-OGS that uses a ray-tracing algorithm based on a minimum time principle and SIRT (Simultaneous Iterative Reconstruction Technique) as an inversion method [30].

Results and Discussion
In Figure 4, we show the velocity profiles from the tomographic inversion of all wavefields relative to the three lines acquired in Kaštela, while in Figure 5, we show those relative to the two lines acquired in Ferrara. The two velocity profiles show very different features. In Kaštela, below a slower weathered layer, the velocity rapidly increases to 4.5 km/s and 3.0 km/s for P-and S-waves, respectively. Furthermore, we observe an overall increase in the velocities in the eastern lines (Lines 2 and 3) compared to the western line 1. In Figure 4d-f, we show the SH-velocity profiles and the corresponding V S30 . Since the velocity is >800 m/s along all three profiles, we conclude that the soil can be classified as A-class following the Eurocode-8 provisions [9]. In the bottom-right corner of each velocity profile, the mean RMS error is reported. All profiles present a mean RMS lower than 5%, except for the P-wave profile of line 2. This is probably due to the low signal-to-noise ratio of the data caused by the rainfall during the acquisition. Having said this, all other profiles show low mean RMS and can therefore be considered reliable.   In Figure 4, we show the velocity profiles relative to all wavefields for the two lines acquired in Ferrara. We observe a horizontal layering of the sediments, which is typical of alluvial plains. The highest velocities are reached in Line 2, where they exceed 2 km/s and 0.3 km/s for P-and S-waves, respectively. The VS30 is shown in Figure 4c,d, and it shows values for which the soil can be identified as C-class (VS30 >180 m/s) following the Eurocode-8 provisions [9], matching the estimates performed in the area by previous large-scale studies [31]. This value and the velocity profile at the center of the line (where the 1D S-wave velocity profile from the MASW is ideally located) is compatible with the results of the inversion of the surface waves, confirming the reliability of our results. The sharp transition at approximately 3 m depth from blue to green in both P-wave velocity profiles most likely indicates the water table. This is consistent with previous geotechnical investigations in the area [32]. Similarly to Figure 4, the mean RMS error is shown in the bottom-right corner of each velocity profile. The values are quite variable, depending on the signal-to-noise ratio of the datasets, which depended on the amount of activities ongoing in the fields. Having said this, they usually show quite a low value, in the range of 5%, and therefore, the velocity profiles can be considered reliable.
Further information about the sites can be obtained by computing the VP/VS ratio (proxy for Poisson's ratio). In Figure 6a-c, we show the VP/VS profiles relative to the three lines acquired in Croatia, while in Figure 6d,e, we show those relative to the lines acquired in Italy. The profiles acquired in Croatia show an overall increase in the VP/VS ratio in the In Figure 4, we show the velocity profiles relative to all wavefields for the two lines acquired in Ferrara. We observe a horizontal layering of the sediments, which is typical of alluvial plains. The highest velocities are reached in Line 2, where they exceed 2 km/s and 0.3 km/s for P-and S-waves, respectively. The V S30 is shown in Figure 4c,d, and it shows values for which the soil can be identified as C-class (V S30 >180 m/s) following the Eurocode-8 provisions [9], matching the estimates performed in the area by previous large-scale studies [31]. This value and the velocity profile at the center of the line (where the 1D S-wave velocity profile from the MASW is ideally located) is compatible with the results of the inversion of the surface waves, confirming the reliability of our results. The sharp transition at approximately 3 m depth from blue to green in both P-wave velocity profiles most likely indicates the water table. This is consistent with previous geotechnical investigations in the area [32]. Similarly to Figure 4, the mean RMS error is shown in the bottom-right corner of each velocity profile. The values are quite variable, depending on the signal-to-noise ratio of the datasets, which depended on the amount of activities ongoing in the fields. Having said this, they usually show quite a low value, in the range of 5%, and therefore, the velocity profiles can be considered reliable.
Further information about the sites can be obtained by computing the V P /V S ratio (proxy for Poisson's ratio). In Figure 6a-c, we show the V P /V S profiles relative to the three lines acquired in Croatia, while in Figure 6d,e, we show those relative to the lines acquired in Italy. The profiles acquired in Croatia show an overall increase in the V P /V S ratio in the eastern lines (especially Line 3) compared to the western line. This indicates a more fractured medium that is possibly saturated with water. As for the Italian site, the V P /V S ratio is very high in both profiles, reaching a maximum in Line 2. The sharp increase in both lines at a depth of approximately 3 m confirms the position of the water table. Furthermore, such high V P /V S probably indicate the presence of water-saturated unconsolidated clay sediments.
Appl. Sci. 2022, 12, 4535 9 of 12 eastern lines (especially Line 3) compared to the western line. This indicates a more fractured medium that is possibly saturated with water. As for the Italian site, the VP/VS ratio is very high in both profiles, reaching a maximum in Line 2. The sharp increase in both lines at a depth of approximately 3 meters confirms the position of the water table. Furthermore, such high VP/VS probably indicate the presence of water-saturated unconsolidated clay sediments. In Figure 7, we show profiles of the VSV/VSH. In the Croatian site, shown in Figure 7ac, the presence of anisotropy is due to the layering of the Eocene flysch [14]. From the geometrical considerations applied to the raypaths, it is possible to estimate the bedding angles (dip and strike), as shown in [14]. As for the Italian site, shown in Figure 7d,e the presence of anisotropy confirms that the sediments are composed of horizontal layers of clay and silt, which are known for producing anisotropy [33,34]  In Figure 7, we show profiles of the V SV /V SH . In the Croatian site, shown in Figure 7a-c, the presence of anisotropy is due to the layering of the Eocene flysch [14]. From the geometrical considerations applied to the raypaths, it is possible to estimate the bedding angles (dip and strike), as shown in [14]. As for the Italian site, shown in Figure 7d,e the presence of anisotropy confirms that the sediments are composed of horizontal layers of clay and silt, which are known for producing anisotropy [33,34]

Conclusions
We presented two case studies from two very different sites in order to validate the method in such different circumstances. The first consisted of three high-resolution seismic lines acquired along the roads of a village built on hard rock in Croatia. The second consisted of two high-resolution seismic lines acquired in a cultivated area in northern Italy. In both cases, we acquired P-waves as well as S-waves, the latter in two different polarizations: parallel (SV) and orthogonal (SH) to the seismic line. Furthermore, surfacewaves were acquired in the two sites. However, in Croatia,no dispersive event could be identified and therefore no results could be obtained. This highlights the importance of acquiring SH-waves. We then performed first break tomography on all lines, for all wavefields.
The results output by the presented survey method are of high value for engineering purposes, as well as for the overall geological and geotechnical characterization of the investigated areas. Specifically, from an engineering point a view, the sites were classified based on their VS30 as required by the current legislation. The method allows to investigate also other important properties of the subsoil, such as the VP/VS and the anisotropy of the medium. From the former, it is possible to infer the fracturing of the medium and its water saturation (i.e., the position of the water table). From the latter in the Croatian site, it was possible to compute the dip and strike angles of the bedding of the flysch, providing therefore a geologically valuable information, while in the Italian site, it confirmed the silty/clay composition of the sediments and their horizontal layering.
The method proposed is non-invasive, replicable and relatively cheap, as the acquisition required only 2-3 days for each seismic line. Furthermore, it provides an in-depth characterization of the composition of the subsurface up to a depth of several tens of meters without the need for expensive procedures such as drilling and coring. Finally, these parameters are given along entire profiles, a few hundred meters long, making it attractive for the extensive characterization of geologically heterogenous areas.

Conclusions
We presented two case studies from two very different sites in order to validate the method in such different circumstances. The first consisted of three high-resolution seismic lines acquired along the roads of a village built on hard rock in Croatia. The second consisted of two high-resolution seismic lines acquired in a cultivated area in northern Italy. In both cases, we acquired P-waves as well as S-waves, the latter in two different polarizations: parallel (SV) and orthogonal (SH) to the seismic line. Furthermore, surfacewaves were acquired in the two sites. However, in Croatia,no dispersive event could be identified and therefore no results could be obtained. This highlights the importance of acquiring SH-waves. We then performed first break tomography on all lines, for all wavefields.
The results output by the presented survey method are of high value for engineering purposes, as well as for the overall geological and geotechnical characterization of the investigated areas. Specifically, from an engineering point a view, the sites were classified based on their V S30 as required by the current legislation. The method allows to investigate also other important properties of the subsoil, such as the V P /V S and the anisotropy of the medium. From the former, it is possible to infer the fracturing of the medium and its water saturation (i.e., the position of the water table). From the latter in the Croatian site, it was possible to compute the dip and strike angles of the bedding of the flysch, providing therefore a geologically valuable information, while in the Italian site, it confirmed the silty/clay composition of the sediments and their horizontal layering.
The method proposed is non-invasive, replicable and relatively cheap, as the acquisition required only 2-3 days for each seismic line. Furthermore, it provides an in-depth characterization of the composition of the subsurface up to a depth of several tens of meters without the need for expensive procedures such as drilling and coring. Finally, these parameters are given along entire profiles, a few hundred meters long, making it attractive for the extensive characterization of geologically heterogenous areas.
Given the value of the results and all the above considerations, we conclude that the method proposed is an optimal tool to perform an in-depth characterization of the shallow layers of the subsurface.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.