Seismic reflection imaging is a well-established tool to determine subsurface boundaries of geological, hydrogeological, or other nature. Applying compressional waves (P-waves) is the common method, e.g. in exploration seismology, for deep targets up to several km. For environmental and engineering geophysics (Steeples and Miller, 1988), it provides images of shallow targets up to several 100 m in high resolution. However, some challenges occur during P-wave application on land that either inflate the acquisition costs or limit the structural resolution. Refracted and guided waves as well as surface waves and source-generated noise diminish the commonly-used part of seismic reflections in the seismograms (Steeples and Miller, 1998). Typically, the different propagation velocities of seismic waves separate the wavefield and allow discriminating different parts of a seismogram for shallow reflection imaging. Hunter et al. (1984) defined one of this parts as an “optimum window” for seismic reflection imaging. In this window, reflections do not interfere with P-wave first arrivals or surface wave ground-roll and can be used for reflection imaging without filtering these signals. Its utilization requires large source-receiver distances (offsets) relative to the depth of the reflectors. On the other hand, dense source and receiver spacings are mandatory to avoid spatial aliasing of the recorded wavefields. The combination of large offsets and dense spacing raises the acquisition costs for near-surface P-wave reflection seismic imaging. In addition, structural resolution is limited by the maximum usable frequency of the source, e.g. explosives or seismic vibrators that possess a maximum frequency range for excitation of the seismic wave (Buness, 2007). An extensive seismic processing can improve the structural resolution of near-surface data (Bradford et al., 2006; Burschil et al., 2018), but only in the mentioned limits of P-waves.
The usage of shear waves (S-waves) in the seismic reflection method (Pugin et al., 2013; Krawczyk et al., 2013) is a tool to overcome the limited resolution of P-waves and to decrease the seismic imaging costs. S-waves propagate slower than P-waves by nature so that they offer higher resolution at a similar frequency bandwidth range (Burschil et al., 2015). Due to the lower propagation velocity, the usable part of S-wave reflections occurs in the seismograms at smaller offsets compared to P-waves. They often show interference with high-amplitude surface waves, where the latter often hide the S-wave reflections due to their significantly stronger amplitudes. Unless the site characteristics involve a high-velocity top layer that suppresses the generation of interfering surface waves during acquisition (e.g. Inazaki, 2004; Polom et al., 2010), surface waves need to be removed during seismic processing workflow (e.g. Wadas et al., 2016).
S-wave seismic reflection imaging has been applied to near-surface either in the horizontally-polarized S-wave domain (Inazaki, 2004; Polom et al., 2010) or in the vertically-polarized S-wave domain (Pugin et al., 2009). Several studies successfully investigated shallow targets of <30 m by structural interpretation of reflection patterns, e.g. of the bedrock topography (Polom et al., 2017), fractured environments in context to sinkhole formation (Wadas et al., 2017; Polom et al., 2018), or pathways for aquifer contamination (Malehmir et al., 2018). Furthermore, seismic velocity can be used for interpretation, e.g. a reduced interval velocity may indicate dissolved rock (Krawczyk et al., 2012; Tschache et al., 2018), or a change in interval velocities within layers may indicate different delta deposits (Winsemann et al., 2018). In addition, ground-engineering parameters, like the dynamic shear modulus, can be derived from the propagation velocity (e.g. Mavko et al., 2009; Uhlemann et al., 2016).
Data quality varies between the above-mentioned S-wave studies, even in the same study area (e.g. Carvalho et al., 2016). For all of these studies, one specifically adjusted seismic processing scheme was applied to process the data to an optimized seismic section. The aim of this paper is to highlight different processing strategies for horizontally-polarized S-wave imaging. The paper reflects our experience in S-wave seismic imaging and does not claim to be comprehensive. The geological target of investigation is an overdeepened Quaternary basin in the Rhine Glacier controlled Alpine foreland, i.e. the Tannwald Basin (TB) (Fig. 1), previously studied using P-wave reflection seismic method (Burschil et al., 2018).
The Alpine Orogeny caused the formation of the Molasse Basin in the northern Alpine foreland (Pfiffner, 1986). During the Pleistocene and the glaciations of the Alps, several ice advances of the Alpine Rhine Glacier moved over the Alpine Molasse Basin (Seguinot et al., 2018) and controlled the landscape evolution. Glacial erosion formed a system of overdeepened valleys and basins (Ellwanger et al., 2011). Glacial, lacustrine, and fluvial deposits sequentially filled these basins. Topographic highs remain in areas between the basins in the Rhine glacier affected area. Classification in terms of erosion and accumulation cycles determines three main formations (Fm.), based on the lithostratigraphy (from oldest to youngest): Dietmanns Fm., Illmensee Fm., and Hassenweiler Fm. Chronostratigraphically, these correspond to periods that comprise subsequent glaciation stages: Hosskirchian-Rissian, Rissian-Wurmian, and Wurmian-Present, respectively.
The studied area of the Tannwald Basin (TB) is located near the margin of the last glacial maximum (Fig. 1). The depth of the north-south-elongated basin varies between 80 m and 240 m. Burschil et al. (2018) interpreted the seismic facies on five seismic P-wave profiles, supported by a research borehole (locations in Fig. 1b). The seismic reflections image the substratum of the molasse units and reveal the sedimentary succession: glacial deposits of basal till, lacustrine deposits, and till sequences of the Dietmanns Fm. (Hosskirchian-Rissian stage), till and till sequences of the Illmensee Fm. (Rissian-Wurmian stage) as well as fluvial deposits of the outwash plain (Wurmian-Present stage). Thus, the TB forms an overdeepened Dietmanns basin and an area between Ilmensee basins, where no overdeepened erosion occurred during this stage.
Reflection seismic method
Since 2014, Leibniz Institute of Applied Geophysics (LIAG) conducted several reflection seismic surveys in the TB. Data of two P-wave surveys using vertically-oriented sources and receivers constitute the base for the seismic interpretation (Burschil et al., 2018). In 2015, we recorded two S-wave profiles along previously measured P-wave profiles. In our acquisition setup, both sources and receivers are oriented perpendicular to the profile direction to enable generation and analysis of horizontally polarized S-waves (SH-waves). The acquisition parameters for P-wave and S-wave surveys are summarized in Table 1. Weather conditions were sunny and calm during all surveys. For the S-wave survey, we acquired 2.1 km of profile in total, in five days of operation, including half a day of source tests. The crew consist of 4 members. In both P-wave surveys, we acquired 7.9 km of profile in total in 16 days of operation with 5 crew members.
In the following, we focus on profile TB-1 of the S-wave survey (Fig. 1b), because it contains various surface conditions (Fig. 2a). LIAG’s horizontal hydraulic vibrator MHV4-S was used as seismic source and geophones mounted on two landstreamers of 120 m length each (1 m receiver spacing), recorded the earth response (Fig. 2a). The ground condition varied along the profile from asphalted road surface via a consolidated gravel path to a field path (Fig. 2b). At each shot location, we recorded four sweeps and changed the sweep polarization after two excitations (++–). The shot gathers show reproducible signals at the same shot point (SP). We chose a source-receiver layout that enables a balanced updip-downdip layer detection sensitivity. During production, SPs were located at 4 m intervals in the first quarter of the 240 m-long spread (split-spread). After 60 m, both landstreamers were moved one quarter of the spread ahead (roll-along). At the SW end of the profile, shots were vibrated over the entire receiver spread plus 60 m end-on SPs. To prevent traffic obstruction while crossing the K7529 road with the landstreamers, we measured two individual profile sections: Both sections start at the road, operating away from the road. To close the coverage gap and to gain the same offset range as within the profile segments, end-on SPs on the opposite sides of the road allowed the fusion of the profiles. The operating and SP arrangement of the profile is shown for the first, two consecutive, and the last spreads in Fig. 2c.
The field data show low environmental noise (Fig. 3) and clear reflections in unprocessed shot gathers from shallow parts up to 0.7 s recording time (blue arrows in Fig. 3). Love surface waves with higher amplitudes occur in areas of unconsolidated ground surface (red arrows in Fig. 3), while they are not obvious in parts with asphalted or consolidated ground conditions (yellow arrows in Fig. 3). At a profile distance of 1420 m, the surface waves in the shot gathers show a significant variation that is not evident in the ground conditions, the landscape, or the geological map (green arrow in Fig. 3). The affected region is approximately 20 m wide. Seismic velocities of the surface waves (~500 m/s) and refracted waves (~700 m/s) are alike on both side of this zone (not shown), whereas the velocity of the direct wave is significantly lower (~250 m/s).
For this comparison study of different processing strategies, we adapted an S-wave processing scheme that was used for various studies (Krawczyk et al., 2012; Wadas et al., 2016; Polom et al., 2018). Field data were processed using the Seisspace® software package. The major processing steps comprised signal enhancement (Fig. 4), normal move-out (NMO) correction with stacking velocities received by an iterative velocity analysis (Fig. 5), common midpoint (CMP) stacking, and poststack time migration. The application of different filters and the adaption of the processing parameters were the main adjustments to highlight different details in the final seismic section (Fig. 6; Table 2; see also Supplementary material).
All processing flows comprised adding geometry to the field data, vertical stacking, eliminating noisy traces, and static correction for the topography (Fig. 4a). After several tests with amplitude scaling, including true amplitude recovery, we decided to use automatic gain control (AGC) with different window lengths, which showed best results (Fig. 4b). Spectral whitening flattened the power spectrum range and increased resolution (Figs. 4c, e). It also reduced high-frequency, source-generated noise, which was undesirably amplified by AGC (green arrow in Fig. 4). A crucial step was the elimination of surface waves, tested by different approaches. Filtering in the f-k domain reduced the surface waves best (Fig. 4d). A selective application of the f-k-filter only to shot gathers, in which surface waves are present, generated post-stack artefacts at the edge of the region where we applied the filter. It resulted in less coherent reflectors in the stacked section.
Several iterations of interactive velocity analysis of the reflection hyperbolas in the CMP gathers were required to gain sufficiently accurate stacking velocities along the whole profile. Subsequently, we derived interval velocities in depth from stacking velocities by the smoothed gradient method (Sheriff and Geldart, 1995). The velocity distribution can be divided into two areas (black dashed line in Fig. 5): (1) an upper unit of 40 m thickness with a mean velocity of 517 m/s and a slightly increasing velocity gradient in depth (466 m/s at surface to 580 m/s at 30 m depth), and (2) a lower unit of 638 m/s with small variation (standard deviation of 25 m/s). However, small-scale, shallow velocity variations within the upper unit were not resolved in the stacking velocities. A tomography study of the S-wave first arrivals did not show a strongly varying near-surface velocity field either. Nonetheless, reflections exhibit deviations up to 10 ms from a hyperbolic shape in the shot gather (blue arrows in Fig. 4) as well as in the CMP gather, which we assign to near-surface velocity anomalies. Residual static corrections after NMO correction in a maximum shift range of ±5 ms reduced these effects without risking cycle shifts. 100% stretch mute applied during NMO correction muted the S-wave first arrivals. It includes a significant wider range for constructive wavelet stacking than is common for P-wave processing (e.g. 30–50%), as a result of the stronger hyperbola curvatures and velocity gradients at depth compared to P-waves.
After NMO correction and CMP stacking, processing comprised f-k-filtering (for one strategy only), finite difference (FD) time migration, data enhancement using a dip scan stack that reduces high-frequency, steep-dipping artefacts in the seismic section, and time-to-depth conversion using a smoothed version of the stacking velocity field.
Processing towards specific imaging purposes
Four different strategies pursue the goals: to highlight small-scale structures (1), locate the origin and incidence of reflections in the seismograms (2), generate a coherent image for the interpretation of large-scale structures (3), and sequentially compare relative amplitudes in P-wave and S-wave images (4).
The seismic section for small-scale structures images the structure of the TB in detail. The basin base is observed between 450 m and 350 m elevation (red arrow in Fig. 6a), i.e. 130 m to 200 m below the ground surface. Short coherent reflectors show internal structures of the basin fines (green arrow and black box in Fig. 6a). The section also reveals the till layers (purple arrow in Fig. 6a) and very shallow reflections inside fluvial deposits just below the surface (blue arrow in Fig. 6a).
To focus on steep angles of incidence and reflection, we restricted the offset range to 60 m (red line in Fig. 4) and processed the section using the same parameters as in the strategy to highlight small-scale structures (1). The seismic section with offset restriction (Fig. 6b) reveals the main structural image as in the entire offset range section (Fig. 6a) except for minor variations. The very shallow part shows the same coherent reflectors (e.g. blue arrows in Figs. 6a, b). Below 560 m elevation, differences in the reflection pattern occur, e.g., the prominent reflectors are less resolved (purple arrows in Figs. 6a, b). S-wave first arrivals cross the 60 m offset mark at ~0.2 s (Fig. 4d). By restricting the offset, parts of earlier reflection hyperbolas (e.g. purple arrows in Fig. 4) are excluded from the stacking process. Significant differences in the imaging quality of the basin base are not obvious (red arrows in Fig. 6a, b), even if prominent parts of the basin base reflection occur outside the 60 m offset range (blue arrows in Fig. 4). In the data beyond 60 m offset and between first breaks and basin base reflection (Fig. 4d), no significant signal contributes to the overall seismic image (cf. Fig. 6a, b).
The processing for a large-scale structural interpretation omits spectral whitening but applies additionally an f-k filtering after CMP stacking that smooth the stacked section prior to the poststack migration. This procedure enhances the reflector coherency of basin base (red arrow in Fig. 6c) as well as the till layers (purple arrows in Fig. 6c). The reflectors appear smoother and of lower central frequency, compared to the small-scale processed section (red/purple arrows in Fig. 6a). The reflection pattern of the basin fines (green arrow in Fig. 6c) and other coherent patterns (cyan arrows in Fig. 6c) show low frequency horizontal reflectors, which are not obvious in the small-scale processed section (green arrow in Fig. 6a). Shallow reflections are imaged as well, but their structural information is less clear compared to the small-scale processed section (blue arrows in Figs. 6a, c).
To enable a structural comparison of P-waves and S-waves, it is required to process the data with a comparable amplitude scaling in depth. Major imaged structures are now the basin base reflector (red arrow in Fig. 6d) and the reflectors of the till layers and sequences (purple arrows in Fig. 6d). Shallow structures are still visible, as for the smaller AGC window (blue arrows in Figs. 6a, d). The internal structure of the basin fines is not visible anymore due to low amplitudes (green arrow in Fig. 6d). The adjusted scaling is therefore important: either to highlight small-scale structures or to structurally compare the seismic section with P-wave data.
Comparison of P-wave and S-wave imaging
The P-wave imaging (Burschil et al., 2018) reveals the sedimentary succession of the TB in detail as coherent reflectors (Fig. 7a). The interpretation is based on a seismic facies analysis aided by a cored research borehole and is described in Burschil et al. (2018). In summary, it shows molasse units below the basin base as well as the basin fill, from bottom to top: allochthon molasse, basal till, basin fines, till layers and sequences, and uppermost gravel/coarse sand. All S-wave sections, processed with different strategies, image the same overall structures (Fig. 6), especially in the large-scale version (Fig. 6c). The relative amplitudes in the P-wave and the S-wave sections show a comparable sequence (Fig. 7a vs. Fig. 6d). However, the benefit of the S-wave imaging becomes visible in the small-scale structure version. The shorter wavelength offers improved resolution and improved imaging of shallow depths up to 30 m below surface compared to the P-wave section (Fig. 7, box 1). In areas of conflicting dips or wedging units, S-waves resolve more structural details and lead to a different interpretation because the smooth reflectors of the P-wave section do not image the wedging of the till layer in detail (Fig. 7, box 2). The detailed vertical offset of reflectors enables the interpretation of faults (blue dashed line in Fig. 7, box 2). In addition, the P-wave section shows some significant areas of low amplitude response, which appear as transparent zones in the image, most likely caused by near-surface scattering (discussed in Burschil et al., 2018). This effect does not occur in the S-wave images, e.g. at the reflector of the basin base (Fig. 7, box 3). Furthermore, internal structures of the basin fines are resolved in the small-scale processing of the S-wave data.
The different processing strategies comprise diverse processing steps and result in different seismic sections. Therefore, one has to keep in mind the applied processing steps when examining or interpreting the seismic sections.
A crucial part of the processing flow was the reduction of surface waves. Their stronger amplitudes compared to the reflections dominate during CMP-stacking and make recognition of reflections difficult. Several studies report a natural suppression of surface waves by very dry and compacted surfaces (Ghose et al., 2013) or by a paved surface layer (e.g. Inazaki, 2004; Krawczyk et al., 2013), as we observe for the asphalted ground condition of the profile (Fig. 3). Several tools are reported to reduce surface waves during processing by applying, e.g. a 45 Hz low-cut filter (Pugin et al., 2013), an unspecified dip filter (Carvalho et al., 2016), or f-k filtering (Beilecke et al., 2016; Polom et al., 2010, Polom et al., 2018). Schmelzbach et al. (2007) used a different approach to reduce linear onsets of Rayleigh surface waves, since f-k-filtering led to artefacts of the spatially irregular sampled data. They modeled Rayleigh surface waves in the τ-p domain and subtracted the synthetic signals from the original data in the time domain. The adaption of this approach did not succeed on real horizontally-polarized S-wave data (Van Zanen, 2004). For our data, this subtraction approach was not successful as well, probably due to the differences between Love and Rayleigh surface waves. Finally, f-k-filtering led to the best results for our (regularly sampled) data. Since an selective application led to artefacts, we applied this filter on all shot gathers over the entire profile, even if surface waves are not obvious in data recorded on consolidated ground (Fig. 4f).
Beside surface waves, the applied f-k filer also removed portions of direct and refracted S-wave events. Top muting by NMO stretch mute ensured that the final section does not contain stacked direct and refracted S-waves. We also tested a manually picked top mute, but it decreased the upper detection limit for shallow reflections significantly. However, a change of shallow seismic reflection signature coincides with the transition of the ground condition from asphalt road to gravel road, and, thus, could be related to stacked S-waves first arrivals. By carefully reviewing the shot gathers and considering all available data, we address this change to the geological structure on top of the terminal moraine of the LGM.
Amplitude scaling is a major issue in seismic processing. The majority of studies that uses S-waves skip true amplitude recovery for different purposes. Krawczyk et al. (2012) applied AGC scaling, because the high noise level in an urban setting impedes a statistically-based surface consistent processing. Brodic et al. (2018) scaled the amplitudes first using spherical divergence correction, but also applied AGC after NMO correction. Pugin et al. (2013) reported amplitude scaling by trace amplitude normalization, but applied AGC of a 1 s window prior to correlation. Ghose et al. (2013) and Carvalho et al. (2016) applied an unspecified gain correction. Beilecke et al. (2016) scaled using spherical divergence and subsequent trace normalization. In our case, true amplitude recovery was not successful, probably due to highly spatial-varying near-surface velocity variations and high amplitudes of source-generated noise. We decided to apply AGC for scaling, since it led to most applicable results. However, if P-wave and S-wave seismic sections are to be compared, one has to take the amplitude scaling into account. For the TB, the basin reflections are at ~0.2 s for P-waves and ~0.6 s for S-waves. The P-wave interval velocity is approximately three times higher than the S-wave velocity. To compare the relative sequential amplitudes of reflections, we set the AGC window to 1.0 s with respect to a 0.3 s AGC window length for the corresponding P-wave processing.
Elevation static correction served the changing topography along the profile. Considering near-surface velocity anomalies that we assume to be present in static corrections, was not possible, since we could not identify first arrival times of S-waves clearly due to interference with P-waves. This is probably the reason, why we do not detect the anomalies in a previous tomography study. However, Pugin et al. (2013) state that near-surface S-wave velocities are more homogenous than P-wave velocities, so that static corrections are not necessary. This is true considering the groundwater boundary, when the P-wave velocity often increase by more than a factor of two. It contrasts with the older work of Edelmann et al. (1983). They found S-wave statics to be 6 times larger than P-wave statics on their profile and argue that the groundwater table tends to equalize P-wave static corrections and does not influence S-wave propagation. Beilecke et al. (2016) required refraction static correction, since the overburden over karstic limestone shows a significant velocity increase along the limestone topography. Carvahlo et al. (2016) reports the successful application of residual static. In the case of our data, the application of residual static corrections did improved the coherency of the reflection signals.
A seismic processing with NMO correction and CMP stacking worked quite well for our site, as for other studies (e.g. Inazaki, 2004; Pugin et al., 2013; Wadas et al., 2016). Neither dip moveout correction nor prestack time migration improved the imaging quality of our data, as reported for study sites with conflicting and crosscutting dips (Carvalho et al., 2016; Beilecke et al., 2016). A test of prestack depth migration (PSDM), which we successfully applied on the P-wave data in this area, failed as well, probably because the offsets used for the tomographic inversion step are considerably smaller than those used in the P-wave data. Since we did not observe conflicting dips in our data and NMO correction with CMP stacking worked satisfactorily, a poststack migration seemed to be adequate.
The processing strategy for interpretation of large-scale structures generates an image of coherent reflections. Since we cannot trace all reflection segments back to the individual shot gathers, the seismic section may contain artefacts that are amplified by the smoothing effect of the second applied f-k filter. To avoid any geological interpretation of artefacts, we recommend to consider the seismic sections of other strategies as well for the final interpretation. Nonetheless, geological large-scale structures produced by this strategy are better interpretable in the seismic images.
We calculate a vertical resolution (VR) for our study, on base of an average velocity in the range of 500 m/s and a dominant frequency around 70 Hz, to be ~1.75 m using the quarter wavelength criterion. However, our shallow reflections show wavelengths as short as 4 m, which corresponds to a VR of ~1 m. This is comparable to the VR reported in other studies. Krawczyk et al. (2012) calculates the VR to be 0.5 m theoretically, but estimates VR to be ~1 m due to other effects. Beilecke et al. (2016) state that VR was ~2 m for their survey, caused by the higher S-wave velocity of ~1100 m/s. Pugin et al. (2013) report a wavelength of 1.2 m to 3 m, which corresponds to a VR of 0.3 m to 0.75 m, based on the quarter wavelength criterion. This is approximately half the wavelength estimated for the corresponding P-wave data of 2–3.5 m maximum VR (Burschil et al., 2018), and, thus, S-waves provide a higher VR resolution than P-waves.
Between distances 1320 m and 1340 m in the profile, the surface waves show a significant lower velocity in an approximately 20 m wide region (Fig. 3). The origin for these anomaly that indicate a shallow low velocity anomaly, remains unclear. A pipeline from the pumping station to an elevated water tank (PW and Wbh in Fig. 1b, respectively) crosses the profile in that area. Nonetheless, the dimension of the affected region argues for a geological structure in the fluvial deposits of the outwash plain. A change of the seismic facies at shallow depth indicates a geological origin as well (cf. Fig. 6).
S-wave reflection imaging using near-offset data only produces a seismic section with a similar structural information for the entire offset range and it motivated us to conduct a 3-D S-wave survey (Burschil et al., 2019). In addition, near-offset S-wave acquisition also reduces costs compared to P-wave acquisition with a layout for the use of the “optimum window”. This is supported by the costs of acquisition in terms of duration and man-power (Table 1). The operation progress of profile meters per day per person is higher for the S-wave survey than for the P-wave surveys for all profiles (117 m/day for S-wave and 99/day for P-wave) and for TB-1 in particular (162 m/day for S-wave and 96 m/day for P-wave). However, to assess precise costs for P-wave and S-wave seismic reflection imaging, other acquisition parameters have to be considered as well, since we used planted geophones for P-wave acquisition and landstreamer technique for S-wave application. However, the use of landstreamers is more feasible for short offset of S-wave recording. Furthermore, S-wave operation can be optimized by using only near-offsets, which we demonstrate to be sufficient for S-wave imaging.
Processing strategies of S-wave reflection data need to be adapted for the specific requirements of seismic imaging and interpretation purposes. Tailored processing strategies serve this issue and lead to seismic images for specific interpretations. For an improved resolution compared to P-waves, processing enhancing small-scale features with short window automatic gain control is able to unfold the significant advantages of S-waves. For interpretation of large-scale structures, we recommend the post-stack application of spatial filters to improve the coherency of reflectors. This strategy is a good choice, if structural interpretation based on P-wave data are not available. Adapted amplitude scaling during processing supports a classification of the reflection sequences for interpretation: display reflection strength vs. emphasize internal small-scale structures.
P-wave reflection seismic imaging is a robust method for near-surface imaging. S-waves complement P-wave imaging by contributing shallower reflectors, improved resolution, and structure imaging, where P-wave imaging fails. The parallel application of both methods is mutually beneficial and, thus, it is our favorable choice. However, using the adequate processing strategy, S-wave imaging is a cost-efficient alternative to high-resolution P-wave reflection imaging, since relatively small offsets produce highly-resolved images. Most information in the S-wave records was found within small offsets <60 m for an imaging depth up to 250 m, after we reduced surface waves by f-k-filtering. In the TB, acquisition costs of S-wave surveying were lower than for P-wave acquisition and could be reduced further by limiting the maximum offset. Nonetheless, S-wave imaging is highly site specific but it is an alternative to P-waves, at least for the presented site.
Source: S-wave seismic imaging of near-surface sediments using tailored processing strategies
Authors: Thomas Burschil, Hermann Buness