Seismic site characterization with shear wave (SH) reflection and refraction methods

Reflection and critically refracted seismic methods use traveltime measurements of body waves propagating between a source and a series of receivers on the ground surface to calculate subsurface velocities. Body wave energy is refracted or reflected at boundaries where there is a change in seismic impedance, defined as the product of material density and seismic velocity. This article provides practical guidance on the use of horizontally propagating shear wave (SH-wave) refraction and reflection methods to determine shear wave velocity as a function of depth for near-surface seismic site characterizations. Method principles and the current state of engineering practice are reviewed, along with discussions of limitations and uncertainty assessments. Typical data collection procedures are described using basic survey equipment, along with information on more advanced applications and emerging technologies. Eight case studies provide examples of the techniques in real-world seismic site characterizations performed in a variety of geological settings.

use of horizontally propagating shear wave (SH-wave) refraction and reflection methods to determine shear wave velocity as a function of depth for near-surface seismic site characterizations. Method principles and the current state of engineering practice are reviewed, along with discussions of limitations and uncertainty assessments. Typical data collection procedures are described using basic survey equipment, along with information on more advanced applications and emerging technologies. Eight case studies provide examples of the techniques in real-world seismic site characterizations performed in a variety of geological settings.
Keywords COSMOS guidelines · Seismic reflection · Seismic refraction · Shear wave velocity · Vs30 · Seismic hazard · Seismic site characterization The near-surface geological conditions at a site can have a significant impact on the amplitude and duration of seismic waves as they reach the surface. To reduce the effects of significant ground motion on structures, modern building codes worldwide have adopted seismic site classification schemes with categories that can be defined by the traveltime-weighted average shear (S) wave velocity from ground surface to 30 m depth, referred to as Vs 30 . Vs 30 is an average velocity calculated by dividing 30 m by the traveltime to that depth, as defined in the National Building Code of Canada (NRC 2010) and by the International Code Council (ICC 2015). This is different from the Abstract Reflection and critically refracted seismic methods use traveltime measurements of body waves propagating between a source and a series of receivers on the ground surface to calculate subsurface velocities. Body wave energy is refracted or reflected at boundaries where there is a change in seismic impedance, defined as the product of material density and seismic velocity. This article provides practical guidance on the Article highlights • A review of the principles and state-of-practice for shear wave refraction and reflection techniques for seismic site characterization. • Guidance is provided on data collection and processing, along with information on technique limitations and uncertainty assessment. 1 Shear wave refraction

Introduction
In accordance with Snell's law, an elastic wavefront will be refracted when it encounters a boundary between two materials with a seismic impedance contrast (Z = density*velocity). The amplitude partition between reflected and refracted incident plane waves is predicted by the Zoeppritz (1919) equations. A nonplanar wavefront radiating from a point source will refract at the critical angle of incidence along the boundary and radiate energy back to the surface, yielding "head-wave" refractions (see Heelan 1953;Brekhovskikh 1960;Červeny and Ravindra 1971). Velocity of, and depth to, the refracting surface can be calculated through traveltime measurements of the seismic wave between the seismic source and the receivers. Sample ray paths of refracted and reflected waves in some typical site conditions are shown in Fig. 1. Seismic refraction methodologies for near-surface investigations were developed over 80 years ago (see Nettleton 1940;Jakosky 1950;Dobrin 1976). The American Society for Testing and Materials (ASTM) standard D5777-18 describes the methodology and equipment requirements for refraction surveying (ASTM D5777-18 2018). Most early near-surface, high-resolution refraction applications used compressional (P) waves generated by vertical impact weight drop or explosive sources and vertically polarized geophones. The refraction procedures discussed here employ a methodology similar to those described in ASTM D5777 for P-waves but instead use polarized S-waves radiating from horizontal Fig. 1 Schematic of S-wave travel paths (reflected and refracted) through a typical sediment-over-bedrock setting. Figure modified from  1 3 Vol.: (0123456789) (SH) sources and horizontal geophones (Hunter et al. , 1998a(Hunter et al. , b, 2002. Shallow refraction techniques are subject to four notable limitations. First, as low-amplitude events, refractions may be difficult to identify in noisy environments where signal-to-noise (SN) ratios are low. Second, the presence of dipping subsurface layers may create significant variation in down-dip or updip apparent velocities, even for relatively low angles. Because dipping layers are extremely common, collecting forward and reverse shot positions for each geophone array is very important if site geometry allows. In settings with common sediment-bedrock velocity contrasts, Nettleton (1940) recommends that the calculated up-and down-dip velocities be averaged to estimate refractor velocities for dip angles of less than 20°. Third, refraction cannot reliably detect velocity reversals (or inversions), because the method is based on the assumption that each layer has a constant velocity and that these layer velocities increase with depth. In settings where higher velocity layers overlie lower velocity layers, other techniques (e.g., downhole shear, seismic cone penetrometer, multichannel analysis of surface waves) may be more appropriate. Finally, a layer must have a minimum thickness to be identified (Xia et al. 2002). This effect is referred to as the "blind zone" or "hidden layer" problem, and an example is given in Sect. 1.2. Generally, the ratio of the geophone array length to the refractor depth should be quite large (> 5) to observe the first critically refracted arrival. In settings where softer sediments overlie a higher velocity rock or till (e.g., Z > 20), shorter arrays can be used, but this increases the possibility of the hidden layer problem.

Survey requirements
Refraction surveys require an array of relatively low frequency horizontal geophones (commonly 4.5 Hz or greater), a seismic cable to connect the geophones, a seismograph, field computer, and a seismic source with a triggering signal pulse (contract closure or pressure trigger). Most engineering seismographs have at least 12 input channels, and many instruments now offer 96 or more. One geophone is needed per channel, oriented transverse to the direction of the array. A 4.5-to 8-kg hammer striking a horizontally imbedded metal I-beam or weighted timber is generally a suitable source for shallow surveys (Fig. 2a, b). To reduce the potential for converted S-to P-wave interference that can occur if the source is struck inline with the receivers, SH-polarized waves should be generated by striking the hammer at right angles to the direction of the array.
When selecting array parameters, the anticipated depth of investigation and velocity-depth distribution should guide the geophone spacing. In near-surface applications, 1-to 5-m separation is usually appropriate. As a minimum, shot locations are recommended at each end of the array and in center of the spread to obtain a pseudo-reversed-refraction record suite. Shots offset from the end of the array can be used in place of moving the geophone spread if the assumption of flat-lying subsurface layers is valid. To minimize unresolved static corrections, a true reversed refraction profile can be collected where the geophone location at each end of the array is replaced by a source location. For SH-refraction surveys, field records should be collected by hammering in one polarizing direction and then turning 180° to record a separate shot(s) in the opposite direction. Reverse polarity data at each shot position can assist with picking refraction arrival times during analysis. To improve the SN ratio and enhance the low-amplitude refractions from distant shot points, records from the same shot position may need to be digitally stacked in the field.

Data processing
The arrival time of shear wave energy can be analyzed using processing software available from many third-party vendors and open-source websites. 'Timedistance' plots are created by graphing traveltime against the distance between the source and the geophones (Fig. 3). These plots are used to calculate an interval velocity for each of the refracting layers from the average of the forward and reverse plots using the inverse of the slope for each layer (see Case Study 1). The "intercept-time" method allows for the calculation of layer thicknesses from the velocities and intercept times of the slopes. The "critical distance" method can also be used to make layered interpretations. These and other methods are documented by Nettleton (1940), Musgrave (1967), Mooney (1984), Palmer (1988), Telford et al. (1990, and others. Software packages (see Sect. 2.2) are available to perform basic filtering, gain adjustments, pick arrival times, and interpret layer velocities and thicknesses. Refraction tomography techniques are becoming more common (e.g., Sheehan et al. 2005). Tomographic methods do not require an assumption of continuous layers of constant velocity and are therefore useful in settings with varying velocity gradients and changing layer thicknesses.
When considering uncertainty in refraction processing, Williams et al. (2003a) estimated that a ± 10% error on refraction velocity measurements is possible based on sensitivity analysis of refraction first break picking. Dipping or irregular layers, velocity reversals, velocity gradients, and lateral heterogeneity in shallow materials are all sources of uncertainty that influence the results of a velocity analysis using refraction methods. Site factors such as background noise levels, a sloping ground surface, practical limitations in array length, and coupling of source to ground surface also contribute to uncertainty. An interpreter also must be aware of interpretation complexities such as possible hidden layers, low-amplitude refractions, first arrival picking errors, and subtleties in selecting segments along the slope. To reduce the impact of these factors on data quality, forward and reverse shots in the field can help identify inclined or irregular layering. Stacking of recordings in the field can improve the identification of low-amplitude refractions. Data should be acquired at times when noise levels are acceptably low, defined here as a timeframe when the refraction signal is visibly as high as reasonable on test field records. Increasing the number of shots and receivers improves the definition of subsurface structure, allowing for more detailed interpretation routines that can better define lateral subsurface variability. Figure 3 presents the potential impact of the hidden layer problem, where a layer of intermediate velocity may be missed from first arrival time picking. In this example, a 17.5-m-thick layer of material with a Vs of 550 m/s would not produce a first arrival. An experienced interpreter may be able to identify the event from the presence of later arriving high-amplitude reflections.
When reporting on results for a shallow seismic site investigation, the equipment and survey parameters (e.g., geophone coordinates, frequency, and orientation) must be described along with any survey limitations (e.g., noise levels, line length restrictions, or array topography). Any information from nearby boreholes suggesting potential for velocity reversals or hidden layers should be described. Sample seismic field records should be included to indicate data quality. Interpretations of picked first arrivals of forward and reverse shots and the resulting time-distance plots with an error analysis of the slopes should be available for review in the report. Finally, the calculation of the average shear wave velocity to the depth of interest (Vs 30 or other) should be demonstrated in a table of interpreted layer thicknesses (z), velocities (Vs), and the calculated traveltime within each layer. Sample calculations are shown in Table 1.

Ottawa, Ontario, Canada: S-wave refraction in a seismic microzonation study
Between 2004 and 2011, the GSC and Carleton University collaborated in a microzonation project investigating the regional variation of soft sediment thickness and shear wave velocity across Canada's National Capital Region (cities of Ottawa, Ontario, and Gatineau, Quebec). During the field program, seismic reflection and refraction data were collected at 750 test sites, and refraction events were visible on records at 527 of these sites. Key project outcomes were published of microzonation maps showing the distribution of seismic site classes and fundamental site periods across the region Motazedian et al. 2020). The large impedance contrasts between the finegrained soft sediment overlying the compact glacial till or bedrock produced high-quality refraction (and reflection) records at many of the test sites (Hunter et al., 2010a). The stratigraphic sequence in the study area consisted of limestone, dolostone, shale, or granitic bedrock, overlain by thin glacial till deposits (average thickness ~ 7 m), capped by Holocene glaciomarine sediments composed of silt-and clay-size grains. Prior to identifying survey locations, surficial geology maps were consulted to select sites representing a range of near-surface materials. A sediment thickness map was also developed from a regional borehole database that provided guidance on survey design. A 3-m geophone separation was typically used, occasionally increased to 5 m in areas where soft sediments exceeded thicknesses of 100 m. Shot points were located at the center of the array and at off-end positions at 1, 1.5, and 10 times the geophone spacing. The selected source was a loaded metal I-beam struck with a 7.3-kg sledgehammer (Fig. 2a). Where bedrock was within 25-30 m of the surface, the 3-m array configuration with off-end shot points captured clear bedrock refractions (Fig. 4).
During data analyses, it was common to observe a high-velocity surface crust of either overconsolidated glaciomarine sediments or fill materials (Vs = 250-400 m/s), shown as the earliest arrivals of the near traces in Fig. 4. This high-velocity layer was approximately 1 to 5 m thick, and neglecting this layer had a limited effect on the interpretation of a traveltime-weighted average Vs interpretation to 30 m depth (Vs 30 ) (Hunter et al. 2010a). The interpreted arrival times were exported from the picking software to produce traveltime-distance plots for forward and reverse shots (Fig. 5). The inverse of the slope was used to calculate the interval Vs for each material. Conventional "layered earth model" refraction interpretation methods (e.g., Telford et al. 1990) were used to calculate the thicknesses of the unit(s) above the bedrock. Differences between the forward and reverse bedrock velocities and time intercepts were interpreted as a dip along the bedrock surface. Forward and reverse bedrock velocities were averaged at sites where dip angles were interpreted to be less than 20° (Fig. 5). At rare sites where the dip was interpreted to be greater than 20°, the traveltime to the bedrock at the center of the array and the apparent velocities calculated at either end of the array were used to calculate a corrected sediment thickness at the array endpoints. This has the effect of deepening the up-dip sediment thicknesses and shallowing the down-dip thicknesses. In cases where site classes were found to be different at either end of the array, the lower of the two site classes was selected.

Fraser River Delta, British Columbia, Canada: S-wave refraction for earthquake hazard studies
The GSC collected and compiled shear wave surface (reflection, refraction) and borehole seismic data in the mid-1980s in the Fraser River Delta, British Columbia, to support regional earthquake hazard studies carried out by government, academic, and industry partners (Hunter et al. 1998a(Hunter et al. , 1998b. Refraction surveys were conducted at 112 sites using an array of 8-Hz geophones oriented in SH mode, with shot positions selected at either end of the array (Fig. 6a). A loaded I-beam was used as a source at sites where firm ground was within 40 m of the surface (e.g., Fig. 2a), and an 8-gauge in-hole shotgun was needed where sediment thickness exceeded 100 m. Interpreted shear wave first arrivals from forward and reverse shots were exported to produce traveltime-distance plots (Fig. 6b). Two interpretation Fig. 4 a A single off-end shot record showing a thin, high-velocity surface crust, the refraction in the soft sediment, and the bedrock refraction. b The time-distance plot created from first arrival picks was used to interpret the interval velocity of each material approaches were used to analyze the refraction datasets: the traditional "layer-case" method and the "velocity-depth" routine developed after Hunter (1971). The latter approach produced a running least-squares fit centered at each of the points on the traveltime-distance curve. This was seen as an advantage over the "layer-case" method which requires an interpreter to select straight line segments through the slope of the traveltime-distance points. The two techniques were shown to produce similar results (Figs. 6c, d), but the "velocity-depth" approach offers a more gradational increase in velocity with depth which is commonly associated with thick, near-surface, soft sediments.

San Fernando Valley, California, USA: S-wave refraction for seismic site effect investigations
The objective of the San Fernando Valley geophysical study was to investigate differences in shallow S-wave velocities at sites with varying amounts of building damage caused by the 1994 M6.7 Northridge earthquake ( Fig. 7a; Williams et al. 1996). Due to generally limited space in this urbanized environment, array lengths ranged between 70 and 100 m for 11 sites, with example sites POT and MCK (shown in Fig. 7b) limited to 74-m offset. These data were collected with 24 14-Hz single-component geophones separated by 3 m. A 4.0-kg sledgehammer striking a 15.2-cm-wide shear timber with metal endplates was used as the source (Fig. 2b). Source locations were at the array midpoint and endpoints with reversed polarity records collected at each location. The maximum depth of investigation at each site was estimated to be approximately 30 m based on the recommendation of Mooney (1984). As described in Williams et al. (1996), the S-wave field records were first read into processing/analysis software and sorted by offset (sourceto-receiver distance) so that the reversed polarity traces are overlain. The S-wave refraction arrivals in the seismic records were identified as the significant and laterally continuous polarity-reversed arrivals (Fig. 7b), which were then picked from the overlain records for each site. Williams et al. (1996) interpreted the refraction data using the slope-intercept method as described by Mooney (1984), which assumes continuous and possibly dipping layers across the length of the spread. Williams et al. (1996) averaged the resultant reflection and refraction depths to produce a onedimensional depth-velocity profile that represents the site in the middle of the geophone spread. The seismic arrays were located on a generally flat surface; thus, no elevation corrections were needed (Williams et al., 1996).

Introduction
Snell's law and the Zoeppritz equations predict the partition of reflected and refracted energy when seismic waves encounter a boundary between two materials with a seismic impedance contrast (Z = density*velocity) (e.g., Telford et al. 1990). Seismic reflection methods measure the time for seismic energy to travel down from a surface source to an impedance boundary and back up to receivers on the ground surface (Fig. 8). The shape of a traveltime curve from a reflection signal on a multi-receiver record is hyperbolic (Fig. 9) and can be used to calculate an average velocity from the ground surface to  (Hunter et al. 1998a(Hunter et al. , 1998b. a Refraction data were acquired at 112 refraction sites throughout the study area. The location of the site shown in this figure is circled. b Traveltime-distance plot of the interpreted first arrivals of the forward and reverse shots. c and d Resulting S-wave velocity-depth profiles, presented as layered interpretations and as "velocity-depth" fits using a routine developed by Hunter (1971). A 5-pt running least-squares fit is applied here along the traveltime-distance curve, yielding a continually varying velocity profile rather than the traditional single-velocity "layers." Figure modified from  1 3 Vol.: (0123456789) the reflecting horizon using the X 2 -T 2 method. For a flat-lying reflector, graphing the source-receiver distance squared (X 2 ) against the reflection arrival time squared (T 2 ) produces a plot with a slope equal to the inverse of the velocity squared (e.g., Telford et al. 1990). Reflection data can be collected using a single array of geophones or continuously collected along a longer survey alignment creating a two-way traveltime cross section of the subsurface. In this survey mode, the signal-to-noise (SN) ratio is improved by stacking data obtained with different source-receiver locations at the same common midpoint (CMP) (Fig. 10). The two-way traveltime sections can then be converted into depth sections using either downhole velocity logs from nearby borehole(s) or velocity-depth functions.
Deep seismic reflection has been a well-established technique in the oil exploration industry for many decades. Shallow reflection techniques were subsequently adapted in the 1980s when technological advances allowed for portable engineering seismographs and field computers (Hunter et al. 1982Doornenbal and Helbig 1983;Knapp and Steeples 1986b). The technique has been successfully applied in numerous shallow subsurface studies (e.g., Steeples and Miller 1990Steeples 1998Steeples , 2005Brouwer and Helbig 1998;Pullan and Hunter 1999;Brabham et al. 2005;Rabbel 2006) and developed into an ASTM standard (ASTM D7128-18 2018). Although seismic reflection surveys conventionally use compressional (P) waves, shallow S-wave reflection methods are also in practice. In unconsolidated materials, the velocity of P-waves is highly dependent on the degree of water saturation or the presence of gas, whereas S-waves are only sensitive to solid grains and are not transmitted through liquid. Because S-waves travel with lower velocities than P-waves, S-wavelengths are relatively shorter than P-wavelengths, often resulting in increased subsurface resolution (Helbig and Mesdag 1982;Stumpel et al. Williams et al. (1996) 1984; Carr et al. 1998;Pugin et al. 2006). Therefore, shallow Vs reflection methods are highly useful for earthquake hazard studies (e.g., Woolery et al. 1993;Harris and Street 1997;Benjumea et al. 2003;Motazedian and Hunter 2008;Harris 2009Harris , 2010Hunter et al. 2010b;Odum et al. 2010). Shallow multicomponent reflection surveying recently is showing great potential (Pugin et al. 2009(Pugin et al. , 2010Pugin and Yilmaz 2019) and is now being adopted into numerous projects with complex near-surface stratigraphy (e.g., Maries et al. 2017;. As with refraction techniques, shallow reflection techniques are subject to limitations. First, to successfully resolve a survey target, it must be large relative to the wavelength (where wavelength is equal to the ratio of velocity over frequency) of the seismic signal. Second, near-surface materials, especially dry, highporosity unconsolidated sediments, can strongly attenuate high-frequency energy. Thus, a major factor in determining survey quality and resolution is the ability to transmit high-frequency energy into the ground, which in turn depends on the near-surface materials, the frequency and energy of the seismic source, the source-receiver spacings, and the coupling of the receivers and source with the ground. Noise levels also affect survey outcomes by reducing the SN ratio. Common sources include wind, local traffic, machinery operating nearby, and even interference from other types of seismic energy, like ground-coupled airwaves or surface waves (Miller et al. 1986). Finally, to calculate accurate velocities, the change in the arrival times of reflection events across the array (defined as the "moveout") must be observed over a large range of source-receiver offsets. For a fixed array length, deeper reflections have a decreased moveout, resulting in reduced accuracy of velocity estimates with increasing depth.

Survey requirements
Shallow seismic reflection surveying requires a seismograph, geophones, multichannel cables to connect the array of receivers, and a seismic source connected to a highly accurate triggering unit. In Fig. 8 a Seismic waves generated on the ground surface travel from the source to an acoustic impedance boundary where they are partially refracted into the lower layer, and partially reflected back toward the surface. b Subsurface travel paths of reflections from a 12-channel field record. Figure  modified from Pullan et al. (2015) traditional applications, geophones are manually pushed into the ground, and the survey progresses by picking up and advancing ("rolling") the cables and geophones down the line. Data are recorded at a series of shot positions along the array, and in this configuration, reflection and refraction data can often be collected together using a single array of receivers. The development of landstreamers has simplified this acquisition process by towing geophone arrays on sleds, an approach proven effective on roadways with motorized vehicles (e.g., Eiken et al. 1989;van der Veen and Green 1998;van der Veen et al. 2001;Inazaki 2004;Pugin et al. 2004Pugin et al. , 2013. The seismic source also influences data quality, collection time, and thus, survey costs. Past comparisons between controlled seismic sources were carried out to study source attributes (Miller et al. 1992(Miller et al. , 1994Doll et al. 1998;van der Veen et al. 2000). Impulsive sources (e.g., sledgehammers, weight drops, explosives) are more traditional sources for shallow seismic surveys (Fig. 11a). Non-destructive vibratory sources (Fig. 11b) are heavy and relatively expensive but produce highly repeatable signals over a broad range of frequencies resulting in higher SN ratios in noisy settings (e.g., Pugin et al. 2013). In recent decades, small, lower cost vibratory sources have been designed for shallow seismic surveys (Ghose et al. 1998;Matsubara et al. 2002;Truskowski et al. 2004;Haines 2006;Pugin et al. 2020).
For a small survey with a single array, 12 to 96 geophones are typically set out, and records are collected using the same procedures described for refraction acquisition. In landstreamer applications, highefficiency data collection is achieved as the source vehicle tows the receiver array and stops at a regular spacing to record data. For shallow surveys, the source-receiver spacing is generally on the order of a few meters, and receiver spacings are from sub-meter to meters.

Data processing
Several low-cost software programs are available to compute simple velocity and depth estimates from seismic reflection records using the X 2 -T 2 method, which assumes a hyperbolic reflection event from a flat-lying reflector. Considerably, more processing  . 10 Subsurface travel paths in a common midpoint (CMP) gather from six different field records. The traces will be processed together to produce one trace on the final seismic reflection section. The number of traces is referred to as the "fold" of the data. Figure modified from Pullan et al. (2015) effort is required when using longer arrays of data to produce a seismic reflection section in two-way traveltime. Data are sorted and grouped into common midpoint (CMP) gathers, a processing approach that improves the SN ratio by stacking a number of traces together with a range of source-receiver spacings. CMP processing requires a series of steps including trace editing, static corrections, filtering, gain scaling, velocity analyses, and normal moveout (NMO) corrections where each trace is corrected for offset according to a velocity-depth function determined from the data or from nearby borehole information. The NMO-corrected traces in each CMP gather are stacked together, producing one reflection trace for the final section. Once a section is created, surface topography adjustments can be applied, and the seismic time section can be converted to a depth section using velocity information. Specialized software is required to process seismic reflection profiles and can range from open-source code to programs costing thousands of dollars.
Two factors are critical when considering uncertainty in processed seismic reflection data: seismic wavelength and the SN ratio. Wavelength is the fundamental property affecting resolution, and thus the uncertainty in velocity and depth estimates. Under optimum conditions, S-wavelengths in near-surface soft sediments may be less than 1 m resulting in highresolution sections, but typically, resolution is lower than this. Uncertainty will similarly increase when interpreting noisy data where the SN ratio is low. Noise can result from a combination of geological factors (e.g., near-surface material property variations creating "statics," rough/steep reflecting horizons), cultural (e.g., traffic, vibrations from machinery), or environmental effects (e.g., rain, wind).
The velocities and depths to reflectors interpreted from reflection data are also affected by uncertainties in arrival time picking of the wavelet, which may experience phase shifts with increasing offset and the assumption of hyperbolic moveout across the receiver spread. Slight overestimations of arrival times commonly occur when the peak or trough of the wavelet is picked rather than the exact shear onset. Errors in picking a consistent reflection signal across the record can also arise when changes in the source-receiver offset create phase shifts, resulting in velocity calculation errors. With careful picking or by estimating the picking delay, these effects can be identified. Finally, the accuracy of a velocity estimate using the X 2 -T 2 method depends on the moveout of arrival times across the array. Uncertainty in velocity calculations can be reduced by increasing both the array length and number of geophones used in the array to improve the definition of the moveout. Using CMP gathers during analysis aids in the assumption of a hyperbolic moveout, but this assumption may be violated in some situations (e.g., refractive effects of a high-velocity surface layer resulting in non-hyperbolic traveltime-distance plots), and velocity estimates generally become less accurate with increasing depth.
When reporting on survey results, digital records, commonly in SEGY format (SEG Technical Standards Committee 2017), must be provided with field notes detailing instrumentation, recording parameters, geodetic coordinates of survey source and receiver locations if collected, and site conditions. A sample field record(s) should be included in a report to indicate data quality and consideration of uncertainties in velocity estimates discussed above. Information about the software and the processing technique or steps followed must be provided with the final datasets.

Ottawa, Ontario, Canada: S-wave reflection in a seismic microzonation study
Between 2004 and 2011, the GSC and Carleton University collaborated in a microzonation project investigating the regional variation of soft sediment thickness and shear wave velocity across Canada's National Capital Region. Project deliverables included regional Vs 30 and fundamental site period maps of the cities of Ottawa and Gatineau, Canada (Motazedian and Hunter 2008;Hunter et al. 2010aHunter et al. , 2010bHunter et al. , 2012Motazedian et al. 2020). Fault-controlled bedrock basins produce rapid lateral changes in the thicknesses of fine-grained (silt and clay) postglacial sediments that overlie thin glacial deposits and bedrock. Seismic reflection/refraction records collected at 750 sites around the cities were analyzed to produce average traveltime-weighted shear wave velocity-depth profiles, from which a seismic site class could be assigned. A 24-channel array of 4.5-Hz horizontal geophones oriented in horizontal shear (SH) mode was laid out at each site with either a 3-or 5-m spacing. A 4.5-kg hammer striking a steel I-beam with one edge dug into the ground in SH orientation was used as the seismic source (see Fig. 2a). Shot locations were at one and one-and-a-half geophone spacings off each end of the array, and one in the center of the array. This survey geometry allowed for trace-to-trace recognition of wide-angle reflections from the glacial or bedrock surface and provided enough sourcegeophone spacing to see bedrock refractions from depths up to 30 m. Where the postglacial sediments were thick, a bedrock refraction first arrival could not be detected. At these sites, only the reflections from the bedrock, glacial, and intra-sedimentary horizons were used to produce the site's average velocity-depth profile (Fig. 12). Other researchers provide additional examples of the application of these techniques (Harris and Street 1997;Williams et al. 1999Williams et al. , 2003aWilliams et al. , 2003b.

Fraser River Delta, British Columbia, Canada: S-wave profiling for seismic hazard modeling
The Fraser River delta in southwestern British Columbia is known to experience seismic wave amplification (Jackson et al. 2017), leading to greater levels of shaking during earthquakes. The delta's location within the most seismically active region in Canada and the rapid population growth on the delta have raised concerns about the area's seismic vulnerability. To aid in the assessment of earthquake hazards in the Fraser delta, seismic reflection surveys were carried out to image the structure of the delta Fig. 12 Reflection records used to generate an average velocity-depth profile for Vs 30 assessment. The sample record shown in a is a forward shot from a deep (81 m) bedrock site, with several intra-sedimentary reflectors. b Time-averaged shear wave velocities were derived from hyperbolic curves fit to the reflections using seismic processing software 1 3 Vol:. (1234567890) sediments (Pullan and Hunter 1987;Pullan et al. 1989). The depth to, shape of, and velocity contrast at the Holocene-Pleistocene boundary, defining the interface between younger, softer postglacial sediments and older, stiffer glacial till, were needed to model ground motion amplification in the delta. Reflection profile alignments were selected near boreholes that intersected the Pleistocene boundary. Four S-wave reflection profiles were collected beneath the Arthur Laing Bridge, connecting Vancouver to the Vancouver International Airport on Sea Island. An array of 24 4.5-Hz horizontal geophones at 3-m spacings was used to record the data. The source was a 4.5-kg hammer striking an I-beam at an offset of 3 m from the end of the array. Profile Laing 1 is shown in Fig. 13. A strong reflector can be seen in the reflection profile, identified as the top of the Pleistocene till in borehole 90-1 located just beyond the southeast end of the profile (Hunter 1995). Other researchers have used S-wave reflection sections in earthquake hazard studies for the interpretation of near-surface faulting and sediment disturbance (e.g., Woolery et al. 1993Woolery et al. , 1996Woolery and Street 2002;Wang et al. 2003;Pullan et al. 2011;Brooks and Pugin 2019).

Ottawa, Ontario, Canada: S-wave profiling using landstreamer arrays in microzonation studies
The GSC routinely uses vibrating sources coupled with a multicomponent landstreamer receiver array to collect both P-and S-wave data at the same time (Pugin et al. 2009(Pugin et al. , 2013(Pugin et al. , 2020Fig. 14a, b). With this acquisition system, it is possible to collect ~ 1000 multichannel records per day, an improvement many times beyond the hand-planting of geophones. In urban settings, the use of a vibratory source also provides greater stacking capability, resulting in less interference from wind or traffic noise. Figure 15 shows a 1-km segment from an S-wave reflection profile that delineated a significant buried bedrock valley underneath a suburb of Ottawa, Ontario ). The section shows the shallow bedrock surface at the south end of the line, increasing in depth across two bedrock benches to a depth of nearly 90 m at the base of a buried escarpment. A nearby borehole confirmed the depths to the top of bedrock and a thin till layer that underlies the thick, soft, postglacial sediments. These soft sediments are associated with ground motion amplification that occurs during weak earthquake events recorded in the Ottawa area Khaheshi Banab et al. 2012;Motazedian et al. 2020).

Combined shear wave refraction and reflection analyses
As geophysical methods continue to advance, the state of practice in shear wave site characterization will evolve toward full wavefield analysis. Numerous investigations have explored joint analysis of reflection and refraction traveltimes, but such methods are not used consistently in the state of practice today. In the following examples, simple approaches are used to compare and combine reflection and refraction traveltime analyses as a crosscheck for shear wave site characterization.  Fig. 16. These data were recorded on 4.5-Hz horizontal geophones spaced at 1.5 m with total source-receiver offsets from 0 to 120 m. A weighted shear timber (Fig. 2b) and 4.0-kg sledgehammer provided the source energy. Williams et al. (2007) used the slope-intercept method as described by Williams et al. (1996) for interpreting these refraction data (Fig. 16). They interpreted the reflections using a reflection curve fitting application (in this case, a graphical X 2 -T 2 analysis approach) commonly available in commercial engineering seismic data processing software. From the reflection curve fitting, the moveout velocities and zero-offset traveltimes were picked and reflector depths were estimated using the Dix equation as described by Dobrin (1976). To get an estimate of the error bounds for both the reflection and refraction data, the interpreter can intentionally mis-position the line fit from the preferred slopes to a tolerable but extreme amount (Williams et al. 2003a).
In this example, the maximum possible velocity variation is about 15% (see Williams et al. (2007)). Given the recorded dominant S-wave bandwidth of Fig. 16  about 20-60 Hz and the interpreted seismic velocity structure, the resulting minimum vertical resolution is about 2 m, based on the quarter wavelength rule (Widess 1973).

Salt Lake Valley, Utah, USA: S-wave data analyses for seismic site characterization
In some geologic settings, S-wave site characterization can resolve velocity structure to depths significantly greater than 30 m. An example of such an investigation is shown in Fig. 17  . Records were acquired at each end of the array and at 100-m intervals out to 300 m from the array ends. For reflection processing, additional records were acquired at every third station within the array to improve CMP coverage. As detailed in Stephenson et al. (2007), the refraction arrivals were first picked and assigned to subsurface seismic layers in an initial refraction model based on traveltime-offset analysis. A robust leastsquares linear regression algorithm was then used to estimate intercept time and apparent refraction velocity for each picked layer. Finally, these intercept times and apparent velocities were modeled using the slope-intercept method to develop a "best" 1D refraction velocity model (Fig. 17a).
The reflection processing was conventional and included amplitude correction, filtering, and velocity analysis prior to the final stacked section. The stacking velocity field was converted to average-interval velocities based on prominent reflection locations in the stack. The resulting profile with the average-interval velocity field overprinted is shown in Fig. 17b, with the refraction model inset in the middle of the stacked section for comparison. Note that many of the refraction layers are imaged in the reflection data as well, and the average velocities of the models are in reasonable agreement. The maximum usable depth of the velocity model is roughly 225 m, based on the deepest prominent reflector depth.

Summary
A key option for seismic site category definition in modern building codes worldwide continues to be the traveltime-weighted measurement of shear wave velocities to a depth of 30 m (Vs 30 ). Non-invasive reflection and refraction techniques are used to directly characterize body wave velocities based on the measurement and analysis of traveltimes recorded at the ground surface. This article has demonstrated that the techniques are quite adaptable for shallow seismic site investigations, ranging from single-array surveys to more complex multi-array towed landstreamer surveys for high-impact projects or complex lithologies. Basic theory, equipment, survey configuration, data processing, limitations, and reporting requirements were discussed to support practitioners in the standardization of the techniques for shallow seismic site classification surveys.
While surface wave methods have gained prominence in site characterization over the past several decades, body waves provide important constraints on subsurface wave speed and wave propagation effects unique from those inferred from surface waves, which often require inverse modeling to obtain shear wave velocity from phase velocity. As the state-of-practice of seismic site characterization continues to evolve, including implementation of full waveform inversion, body wave methods will need to play a critical role in that evolution.

Conflict of interest
The authors from the GSC and the U.S. Geological Survey have no conflicts of interest to declare. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. or Canadian Governments.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.