Effects of surface topography and subsurface structural irregularities of river embankments on Rayleigh wave propagation

A numerical investigation of surface wave propagation in the surface wave exploration of river embankments is implemented using 3D shallow soil models with different subsurface shapes and S-wave velocities. The focus is on the effects of surface and subsurface irregularities on the propagation and dispersion characteristics of Rayleigh waves generated by a surface source and observed at points along the crest of the embankment. Rayleigh wave dispersion curves observed along the embankment crest are found to be contaminated by reflected waves, with the amount of contamination affected more by the position of the velocity boundary at the bottom edge of the embankment than the surface topography. When the velocity boundary is close to the ground surface, especially at the toe of the slope, the effect of reflected waves is sufficiently large to generate biases in the phase velocity estimates.


Introduction
Surface wave exploration is a method used in various geotechnical and geophysical fields for estimating a subsurface S-wave velocity profile in shallow soil layers using the dispersive property of surface waves.Rayleigh waves are generally used because they are easily generated by vertical impacts (e.g., Socco et al. 2010).More than twothirds of the seismic energy propagates as Rayleigh waves in the case of a vertical impulsive force on the surface, such as with a sledge hammer or a dropped weight (e.g., Richart et al. 1970).Since a phase velocity dispersion of Rayleigh waves is frequency-dependent and also related with a S-wave velocity distribution beneath an observation point (e.g., Okada 2003), a 1D subsurface S-wave velocity profile in the depth direction can usually be estimated from the dispersive features of the phase velocity.Nazarian et al. (1983) introduced a method of spectral analysis of surface waves (SASW) to estimate a 1D S-wave velocity profile to a depth of 100 m from Rayleigh wave phase velocities recorded at two receivers.This method has been used in many geotechnical engineering fields (e.g.Stokoe et al. 1994).Park et al. (1999) proposed the multichannel analysis of surface waves (MASW) to determine phase velocities directly from surface wave data at many points.This method has widely been used in soil surveys of, for example, river embankments, landfill and residential land.In MASW, Rayleigh wave phase velocity dispersion curves are estimated with several methods: the F-K transform (e.g.Yilmaz 1987), the tau-p transform (McMechan and Yedlin 1981), phase shift (Park et al. 1998) and the slant stacking algorithm (Xia et al. 2007).Then, a 1D S-wave profile is obtained from an inversion of the dispersion curve (e.g., Xia et al. 1999).A 2D S-wave velocity structure may be also constructed by synthesizing 1D profiles (Miller et al. 1999;Xia et al. 2004;Hayashi and Suzuki 2004;Luo et al. 2009).
In data analysis in the surface wave method, we usually assume a horizontally layered model with a flat ground surface (Xia et al. 1999).Namely, the surface wave theory for a horizontally layered model is used for analysis even if the soil layers actually have surface and subsurface irregularities.An irregularity near an observation point causes wave-scattering that may contaminate the Rayleigh wave dispersion curve that would be expected in a horizontally layered model.Several studies have investigated the applicability of the surface wave analysis to soil layers with surface irregularities or subsurface structural irregularities.Hayashi and Inazaki (2007) performed a numerical experiment with 3D numerical models as well as a field test to investigate the effects of the slope of a river embankment on surface wave propagation.This work confirmed that the effects of the slope on the surface wave dispersion curves were minor.Kim et al. (2010) performed numerical experiments of surface wave propagation at the crest and on the inclined face of an existing rock-fill dam, and verified the applicability of the surface wave method in this case.In order to improve the dispersion curve estimation in the analysis, Hayashi and Suzuki (2004) proposed using common midpoint cross-correlation analysis (CMPCC) in which mulch-channel and multi-shot surface data are used in a cross-correlation.This method can improve the accuracy and resolution of Rayleigh wave phase velocities in a model with a lateral variation in S-wave velocity along an observation line.In fact, most previous studies have focused on examining the applicability of the surface wave method at sites with surface or subsurface velocity variations in the direction of the observation line.However, their effects in the direction perpendicular to the observation line have not been sufficiently investigated.
In this study, we conduct numerical experiments to investigate the effects of surface topography and subsurface velocity irregularities in river embankment models on the propagation of Rayleigh waves generated by a surface source.In particular, we discuss the effect on the surface wave dispersion of reflected waves resulting from 3D irregularity that is present outside the observation line.We use a 3D finite difference method (FDM) to calculate synthetic data in the experiments.The characteristics of the wave propagation and the Rayleigh wave dispersion features are discussed on the basis of the synthetic data to verify the assumption of a horizontally layered structure as often used in the surface wave exploration.

Method Models
We created 3D shallow soil models with shapes similar to a river embankment.Although embankments may have inside layers with lower S-wave velocity than the surface layer (Hayashi and Inazaki 2007), we consider only two-layer models in which the S-wave velocity of the first layer is smaller than that of the second layer.Figure 1 shows the cross-sections of the two models.The physical parameters are constant along the longitudinal axis of the models.Embankment height and crest width are fixed at 4 m in both models, and the ratio of the height to the horizontal slope length is 1:2.The transverse direction of the embankment models is defined as the X axis, the longitudinal direction is regarded as the Y axis.The vertical direction is the Z axis.In model A, shown in Fig. 1a, the boundary between the first and second layers is horizontal at the central part of the embankment, and terminated at the bottom edges in the left and right sides.On the other hand, in model B, as shown in Fig. 1b, the boundary is not horizontal and the bottom edges are not included.The thickness of the first layer just below the crest of the embankment is 4 m in the two models.In both models, the P-wave and S-wave velocities of the first layer are 260 m/s and 150 m/s, respectively, assuming that P-wave velocity is 1.73 times the S-wave velocity.We, furthermore, assume the density of the first layer is 1600 kg/m 3 .Table 1 shows the physical parameters of the models.We do not consider any attenuation in the model.

Finite difference calculation
We employed a 3D staggered grid FDM (Graves 1996) of the second-order accuracy in time and the 8th-order accuracy in space to generate synthetic waves.The volume of the model is 75 m (L) × 75 m (W) × 37.5 m (D).
The finite difference (FD) grid spacing is 0.125 m in all directions and the number of grid points per shortest wavelength is 24.The time step, Δt, was determined for each velocity model using the following formula.
where v pmax is the maximum P-wave velocity of the model and h is the grid spacing.To avoid artificial reflections from the edges of the numerical models, we defined an absorbing boundary (Clayton and Engquist 1977) and a buffer region (Cerjan et al. 1985) on the four sides and the bottom of each model.A vacuum formulation (Zeng et al. 2012) was employed to satisfy the free-stress condition on the model surface; that is, the P-and S-wave velocity and density on the model surface were set to zero.
In the numerical experiments, we analyzed the synthetic vertical waves in a frequency range of 5-50 Hz.In order to reduce the effects of the numerical dispersion of the synthetic waves in the finite difference (FD) calculation, a sufficient number of the grid points per minimum wavelength is required (e.g., Mittet 2002).Hence, to validate the relationship between the accuracy of the FD approximation and the grid spacing, we conducted a numerical test as detailed in the Appendix.The synthetic waves were generated on the surface of a homogeneous (1) �t < 0.669 h model with several choices of the FD approximation in terms of the order of the spatial finite difference and the grid intervals.This confirmed that the synthetic wave amplitudes remain stable if the FD accuracy of the spatial partial derivatives is more than 4th-order and the number of the grid points per the minimum wavelength exceeds 15.In consideration of reasonable computing time, the 8th-order spatial derivative accuracy was selected with a grid spacing of 0.125 m.As already noted, the number of the grid points per the shortest wavelength (at 50 Hz) was 24.

Source and observation points
The source was placed on the embankment crest as shown in Fig. 2a.A Ricker wavelet with a central frequency of 20 Hz was applied as the source-time function for the vertical source on the surface in order to simulate an impulsive impact such as with a hammer.The Ricker wavelet is defined by where f is the peak frequency, t is the time variable, and t 0 is the time delay.Synthetic waveforms were calculated for a duration of 0.7 s.One set of observation points was placed linearly along the center of the embankment crest at a spacing of 1 m, as marked by (1) in Fig. 2a; the vertical component of the synthetic ground velocity was calculated at each point.We also placed two circular arrays on the crest, as marked by ( 2) and (3) in Fig. 2a, at distances of 12.5 m and 22.5 m from the source, respectively.The observation points in these two arrays are indicated (2) in Fig. 2b.The purpose of these two arrays is to understand the propagation direction of scattered and reflected waves that may be generated by the surface and subsurface irregularities.

Results of numerical tests with model A Synthetic waveforms
Firstly, we explain the results of the numerical experiments using model A. Figure 3a shows snapshots of the synthetic velocities on the ground surface.Red arrows in the figure represent the arrival of the Rayleigh wave that propagates along the surface from the source.We describe these as "direct Rayleigh waves".Green arrows show that the direct Rayleigh wave is scattered or reflected at the bottom edges of the embankment and propagates back to the crest again.Thus, in this paper, we call these "reflected Rayleigh waves".Figure 3b shows snapshots of the synthetic velocities of the embankment cross-section including the source point.Green arrows represent waves that are reflected at the bottom edges of the embankment and propagate back to the crest.
Figure 4 shows the vertical components of the synthetic velocities as recorded by the linear array.The red arrow marks the initial rise point of the direct Rayleigh wave propagating along the embankment crest from the source, while the later phase indicated by the green arrow represents the arrival of the reflected Rayleigh waves.The direct and reflected Rayleigh waves can be easily defined by their arrival times at the points near the source.However, as the distance from the source increases, the phases of these waves overlap because of similar their arrival times.

Rayleigh wave dispersion curves
Dispersion curves of the phase velocities were estimated from the synthetic vertical motions for model A using Figure 5 shows the dispersion images of the synthetic velocities shown in Fig. 4 in the frequency-velocity domain.The gray shading in the figure represents the normalized power at each frequency and phase velocity, while the open red dots indicate the derived dispersion curve of the phase velocity with the maximum power.The solid line represents the theoretical dispersion curve of the fundamental mode Rayleigh waves for the 1D profile beneath the center of the arrays in model A. The two dotted lines indicate the phase velocity from the capable maximum wavelength of Rayleigh waves which is same as the array length, and the capable minimum one which is twice the array spacing (Foti et al. 2017).We focus only on the frequency band between the two dotted lines which is the capable wavelength.The observed dispersion curve for the waves in the range of 5-25 m (Fig. 5a) fits well with the theoretical dispersion curve at frequencies below 45 Hz.However, the phase velocities are larger than the theoretical values at frequencies higher than 45 Hz.Thus, the effect of the reflected Rayleigh wave may be large at high frequency.The observed dispersion curves in a distance range of 5-15 m fit well with the theoretical ones at any frequencies, as shown in Fig. 5b.Indicating minor effect of reflected Rayleigh waves.However, the differences between the observed and theoretical dispersion curves in Fig. 5c (10-20 m) are large at frequencies more than 40 Hz.In particular, the discrepancy between the two values is more than 80 m/s at frequencies of 40 to 50 Hz.The theoretical wavelength of the fundamental mode Rayleigh waves in this frequency band is about 3 m.This means that the contribution of the reflected Rayleigh waves to the wave field is larger than that of the direct Rayleigh waves when the distance from the source is more than about five times the wavelength of the Rayleigh waves.
A large discrepancy is also seen in the phase velocity at a distance range of 15-25 m (Fig. 5d).The difference is about 30 m/s at frequencies of about 24 Hz, and 30 to 130 m/s at frequencies of 35 to 50 Hz.The wavelength of the Rayleigh wave at a frequency of 24 Hz is about 6 m, and the distance from the source is 2.5 times the wavelength for the Rayleigh wave.When the wavelength of the observed Rayleigh wave is not very short as compared with the distance from the source, the direct Rayleigh wave cannot be contaminated by the reflected Rayleigh waves.

Wave propagation characteristics
F-K power spectra (Capon 1969) were calculated from the synthetic waves recorded in the circular arrays to investigate the propagation direction of the waves.Figure 6 shows the F-K spectra at frequencies of 15.9, 26.8, 36.6 and 47.6 Hz in the two circular arrays marked as ( 2) and (3) in Fig. 2a.The horizontal and vertical axes represent the wavenumbers in the perpendicular and longitudinal axes of the embankment, respectively.The gray shading in the figure indicates the normalized power of the synthetic velocity at each wavenumber.The F-K spectra from the circular array (2) show only a single peak related to the direct Rayleigh waves propagating in the longitude direction at each frequency, as shown by the red rings in Fig. 6a-d.Similarly, for the circular array (3), only single peaks of the direct Rayleigh waves are seen at frequencies of 15.9 Hz, 26.8 Hz and 36.6 Hz (Fig. 6e-g).On the other hand, three peaks can be identified in the F-K spectrum at a frequency of 47.6 Hz in Fig. 6h.The two peaks marked with green rings in the figure correspond to arrivals of the reflected surface waves from the velocity boundary at the bottom edges of the embankment, which propagate at an angle to the longitudinal axis.The white dotted line of Fig. 6h represents the wavenumber of the fundamental mode Rayleigh wave.These results indicate that the waves observed at points far from the source are more contaminated by the fundamental mode of Rayleigh waves generated at the bottom edges of the embankment, especially at a high frequency.

Dispersion curves at high frequency
Because the significant effect of reflected Rayleigh waves on the dispersion curves was identified in the wave fields at the high frequency (Fig. 5), we calculated dispersion curves for every 5 m ranges along the survey line to focus on the high-frequency components.Figure 7 shows the dispersion curves of the synthetic waves recorded in the ranges of 5-10 m, 8-13 m, 11-16 m, 14-19 m, 17-22 m and 20-25 m from the source.The plotted points and solid line represent the observed and theoretical phase velocities, respectively.The two dotted lines indicate the phase velocity estimated from the capable maximum and minimum wavelengths of Rayleigh wave.The dispersion curves of the phase velocity obtained from the synthetic waves at relatively small distances from the source (in the ranges 5-10 m and 8-13 m) fit well with the theoretical ones at any frequency.In the distance range of 11-16 m, the differences in the phase velocity are more than 100 m/s at frequencies above 45 Hz.At greater distances (in the ranges 14-19 m, 17-22 m and 20-25 m), these differences in the phase velocity are larger at frequencies above 30 Hz, reaching 100 m/s or more at some frequencies.The wavelength of the Rayleigh wave at a frequency of 30 Hz is 4.5 m.When the distance from the source is three times this wavelength, the reflected Rayleigh waves become    of the direct Rayleigh wave.Green arrows represent those of the reflected Rayleigh waves.The amplitude of the reflected Rayleigh waves is considered to be smaller in model B than A. Figure 8b shows snapshots of the synthetic velocities of the embankment cross-section including the source.The reflected Rayleigh wave generated at the bottom edge of the embankment has a small amplitude and difficult to see. Figure 9 shows the vertical component of the synthetic velocities calculated at the points along the linear array in model B. A red arrow marks the initial rise point of the arrival of the direct Rayleigh wave propagating along the linear array from the source, while a green arrow represents the arrival of the reflected Rayleigh waves.The amplitudes of the reflected Rayleigh waves in this figure are smaller than those in model A (c.f. Figure 4).Models A and B have the same surface topography, but the positions of the edge of the velocity boundaries are different.Therefore, it is clear that the generation of the reflected Rayleigh waves depends more on the position and shape of the velocity boundary than on the shape of the embankment surface.

Rayleigh waves dispersion curves
Figure 10 presents dispersion images of the synthetic velocities shown in Fig. 9 in a distance range of 20 or 10 m from the source along the linear array.The difference between the observed dispersion curves and the theoretical ones is small at any frequency.The influence of the reflected Rayleigh waves on the phase velocity is less than in model A.
Figure 11 shows the dispersion curves of the phase velocity obtained from the synthetic waves recorded in a range of 5 m along the survey line.The dispersion curves obtained at relatively small distances from the source (5-10 m and 8-13 m) fit well with the theoretical values at any frequencies.The phase velocity differences at distances of 11-25 m are about 20 m/s at frequencies above 40 Hz.The increase in the discrepancy between the observed and theoretical phase velocities with distances from the source is much less than model A.

Realistic embankment model
In this section, we discuss the effects of subsurface irregularities in the case of a realistic model of a river embankment.Figure 12a shows an S-wave velocity structure of a real river embankment obtained by an S-wave tomography (Suzuki and Kanazawa 2017).A model of embankment, shown in Fig. 12b, was developed from characterization with homogeneous layers in Fig. 12a.We refer to this as model C. The surface topography of this model is the same as models A and B, while the velocity boundary is set at different depths on the river (right) side and on the landward (left) side of the embankment.The physical parameters are the same as in the previous models.

Synthetic waveforms
Figure 13 shows snapshots of the surface velocities.

Rayleigh wave dispersion curves
Figure 15 shows the dispersion curves of the Rayleigh wave phase velocities for the synthetic waveforms in Fig. 14.When the distance range is 5-25 m (Fig. 15a), the discrepancy between the observed and theoretical phase velocities is about 40 m/s at frequencies below 20 Hz.The observed dispersion curves are similar to the theoretical ones when the distance range from the source is 5-15 m (Fig. 15b).This result is similar to those obtained in the experiments using models A and B. When the distance range is 10-20 m (Fig. 15c), a small discrepancy (< 30 m/s) between the observed and theoretical phase velocities is seen at frequencies of 20 Hz.The theoretical wavelengths of the fundamental mode Rayleigh wave are about 10 m at 20 Hz.A discrepancy was found when the phase velocity was obtained from the synthetic waveform at offset of more than 15 m.Therefore, the reflected Rayleigh waves may be dominate over the direct Rayleigh wave when the source distance is more than 1.5 times the wavelength.In the distance range of 15-25 m (Fig. 15d), we can also see a small discrepancy of about 20 m/s at frequencies of 20 Hz.The discrepancy can be found in models A and C when the offset distance is long.However, the frequency at which the discrepancy between the observed and theoretical phase velocities for models A differs from that of model C. It is considered that the shape, depth and thickness of the edge part is so important for the frequency of the reflected Rayleigh waves.

Dispersion curve at high frequency
Figure 16 shows the dispersion curves obtained from the synthetic waves recorded in range of 5 m along the survey line.The observed phase velocities fit well with the theoretical ones at all frequencies and the ranges of array.The effects of the reflected Rayleigh waves are not significant at the high frequency in this model (Table 1).

Discussion
In models A and B, reflected Rayleigh waves were found in addition to the direct Rayleigh waves.The reflected Rayleigh waves have higher amplitude in model A than The main difference between the two models is the position of the velocity boundary; when the velocity boundary is located on the ground surface and the edge has a role to trap the waves inside the surface layers, the reflected surface waves generally have higher amplitude (e.g., Kawase and Sato 1992).In model A, the velocity boundary is at the ground surface at the bottom of the embankment.It can easily and effectively generate and trap the Rayleigh wave within the embankment.However, the velocity boundary in model B is in the subsurface and the amplitude of the reflected wave is small.The reflected Rayleigh wave from the left slope of the embankment is only dominant in model C, where the velocity boundary is close to the bottom of the embankment.These results indicate that reflected Rayleigh waves are likely to arise when the velocity boundary is located near the bottom edge of the embankment and the layer is narrow.The actual topographical shape of the embankment itself does not have a large effect on the generation of reflected waves.Table 2 The ratio of the observed and the theoretical phase velocity difference

Table 1 Physical parameters of embankment models
The unit is %

Source distance range [m]
Frequency Table 2 shows the differences value between observed and theoretical phase velocities for each model, source offset distance and frequency.The incorrect estimations of the direct Rayleigh wave phase velocity occurred in model A and C because of the contamination of the direct Rayleigh wave by the reflected Rayleigh waves.In both models, the difference between the observed phase velocity and the theoretical (estimated) one increases with distances from the source.This tendency was especially significant at high frequencies (30 Hz <) in model A and at low frequencies (< 30 Hz) in model C. Also, in model A, when the source offset distance is more than 15 m, a large difference was seen even at low frequency.In model B, we found little effect of the reflected Rayleigh wave on the phase velocity estimation.These demonstrate that the effects of reflected Rayleigh waves on observed phase velocity depend on the location and shape of the velocity boundary.
However, when the source offset distance was within 1.5 times the wavelength in model A, and 2.5 times in model C, the reflected Rayleigh wave does not contaminate the direct wave, and an appropriate phase velocity was obtained from the synthetic waveforms.The minimum source offset distance must be considered as well as the maximum one.Rayleigh waves are not well developed if the distance from the source is less than half of the wavelength (Stokoe et al. 1994).Therefore, it is necessary to use waveforms obtained where the source offset distance is at least 0.5 to 1.5 times of the wavelengths of Rayleigh waves to calculate the phase velocity.

Conclusions
Numerical experiments were conducted to investigate the assumptions of the horizontally layering in the application of the surface wave method to river embankments.The characteristics of the synthetic wave propagation and the Rayleigh wave phase velocity were investigated for several models with different subsurface irregularities.All models have exactly same surface topography.The results show that the Rayleigh wave propagating along the observation line on the embankment crest can be contaminated by the reflected Rayleigh waves generated at the bottom edges of the embankment.The contamination is pronounced when the velocity boundary is located on or near the bottom edges of the embankment, while the reflected Rayleigh waves are less dominant in models with deeper velocity boundaries.Therefore, the role of subsurface structural irregularities is more significant than that of surface topographical irregularities in the generation of the reflected Rayleigh waves in embankments.
We also found that biases in the phase velocity estimates of Rayleigh waves could depend on the source offset distance.In our numerical models, even if a large reflected Rayleigh wave was generated, it is effective to use the waveform record at a source offset distance of at most 1.5 times the wavelength of the Rayleigh wave to estimate the phase velocity.In addition, it is necessary to excluded waveforms obtained where the source offset distance is less than 0.5 times of the wavelengths of Rayleigh waves from the phase velocity calculation, because Rayleigh wave does not occur at that distance.

Appendix
In this study, we employ a 3D finite difference method (FDM) with 8th-order accuracy in space and secondorder accuracy in time to solve the equation of motion.The grid spacing is set to 0.125 m in all three directions.Prior to starting the numerical experiments, the accuracy of the finite difference (FD) calculations was confirmed using several FD schemes.The candidates for the spatial finite difference schemes in this confirmation process were 4th-order, 8th-order (Furumura 2009) and 16thorder approximations.We refer the above existing FD equations and only explain the 16-order approximation in the following.
In the FDM, a 3D region is discretized with a fine grid spacing (�x, �y, �z) .The 16th-order spatial differentials based on the 3D staggered FDM are defined by (3)

Fig. 1
Fig. 1 Cross-sections of embankment models used in numerical experiments.Embankment height and crest width are 4 m.The gradient of the embankment slope is 1:2.The first layer immediately below the embankment crest has the same depth in both models.a Model A. The boundary between the first and second layers is horizontal and at the bottom of the embankment.b Model B. The boundary is not horizontal and is partially in the subsurface

Fig. 3
Fig. 3 Snapshots of synthetic vertical velocities (a) on the surface and (b) in cross-section including the source point of model A. Colors represent the amplitude of the vertical component.Red and green arrows indicate the direct wave and reflected Rayleigh waves, respectively

Fig. 4
Fig. 4 Synthetic vertical velocities along the linear array on the crest of model A. Red and green arrows indicate the direct wave and reflected Rayleigh waves, respectively

Fig. 5
Fig. 5 Rayleigh wave phase velocities from vertical motions in Fig. 4. Phase velocity analysis of synthetic waveforms in 20 or 10 m ranges from the source: a 5-25 m, b 5-10 m, c 10-20 m; and d 15-25 m.Gray shading represents normalized power of synthetic velocity at each frequency.Open red dots represent observed phase velocity.Solid lines are the theoretical dispersion curves of the Rayleigh wave fundamental mode for model A. The two dotted lines indicate the maximum and minimum wavelengths resolved by the linear array.The frequency bands of capable wavelength are focused

Figure
Figure 8a shows snapshots of the synthetic velocities on the surface of model B. Red arrows represent the arrivals

Fig. 6 F
Fig.6F-K spectra from synthetic velocities recorded at frequencies of (a) 15.9, (b) 26.8, (c) 36.6 and (d) 47.6 Hz in array (2) and at frequencies of (e) 15.9, (f) 26.8, (g) 36.6 and (h) 47.6 Hz in array (3).Horizontal and vertical axes are wavenumbers in the perpendicular and longitudinal axes of the embankment.Gray shading represents the normalized power of the synthetic waveform at each wavenumber.The white dotted line in h) represents the wavenumber of the fundamental mode Rayleigh wave

Fig. 7
Fig. 7 Rayleigh wave phase velocity from synthetic vertical motions in Fig. 3 in 5 m ranges of distance from the source (5-10 m, 8-13 m, 11-16 m, 14-19 m, 17-22 m and 20-25 m).Plotted points represent observed phase velocity in each distance range and the solid line represents the theoretical Rayleigh wave phase velocity.The two dotted lines indicate the maximum and minimum wavelengths resolved by the linear array.The frequency bands of capable wavelength are focused

Fig. 8
Fig. 8 Snapshots of synthetic vertical velocities (a) on surface and (b) in cross-section including the source point of model B. Colors represent the amplitude of waves propagating along the surface

Fig. 10
Fig. 10 Dispersion features of Rayleigh wave phase velocities obtained from synthetic vertical velocities in Fig. 9 for model B. Synthetic waveforms are analyzed in 20 or 10 m ranges from the source: a 5-25 m, b 5-15 m, c 10-20 m, d 15-25 m.Gray shading represents normalized power of synthetic velocity at each frequency.Open red dots represent observed phase velocities and the solid lines are the theoretical dispersion curves of the Rayleigh wave fundamental mode for the 1D profile beneath the linear array in the model.The two dotted lines indicate the maximum and minimum wavelengths resolved by the linear array.The frequency bands of capable wavelength are focused Figure 13 shows snapshots of the surface velocities.Figure 14 shows the synthetic waveforms recorded in the linear array on the crest of model C. Travel times of the Rayleigh waves propagating along the left and right slopes of the embankment are different, with amplitudes on the left slope larger than those on the right slope.This is due to the difference in the position of the bottom edges of the velocity boundary.Since the velocity boundary at the bottom edge on the left side of the embankment is shallower and the layer is thinner than that on the right side, the Rayleigh wave generated at the bottom left edge of the embankment propagates on the slope in the direction of the crest with the large amplitude.

Fig. 11
Fig. 11 Rayleigh wave phase velocity from synthetic vertical velocity in Fig. 9 for model B. Synthetic waveforms are analyzed in ranges of 5 m from the source: 5-10 m, 8-13 m, 11-16 m, 14-19 m, 17-22 m, 20-25 m.Plotted points represent observed phase velocities for the synthetic waves recorded at each source distance.The solid line represents the theoretical phase velocity of the fundamental mode Rayleigh waves.The two dotted lines indicate the maximum and minimum wavelengths resolved by the linear array, respectively.The frequency bands of capable wavelength are focused

Fig. 13
Fig. 13 Snapshots of synthetic vertical velocities (a) on surface and (b) in cross-section including the source point of model C. Colors represent the amplitude of the vertical component

Fig. 15
Fig. 15 Rayleigh waves phase velocities from synthetic vertical velocities in Fig. 14 for model C. Synthetic waveforms are analyzed in 20 or 10 m ranges from the source: a 5-25 m, b 5-15 m, c 10-20 m, d 15-25 m).Gray shading represents normalized power of synthetic waveform velocity at each frequency.Open red dots represent observed phase velocities and solid lines are the theoretical dispersion curves of the Rayleigh wave fundamental mode for the 1D profile beneath the linear array in the model.The two dotted lines indicate the maximum and minimum wavelengths resolved by the linear array, respectively.The frequency bands of capable wavelength are focused

Fig. 16
Fig. 16 Rayleigh waves phase velocity from synthetic vertical velocity in Fig. 14 for model C. Synthetic waveforms are analyzed in ranges of 5 m from the source: 5-10 m, 8-13 m, 11-16 m, 14-19 m, 17-22 m and 20-25 m.Plotted points represent the observed phase velocity recorded at each source distance.The solid line represents the theoretical dispersion curve of the fundamental mode Rayleigh wave for the model.The two dotted lines indicate the maximum and minimum wavelengths resolved by the linear array, respectively.The frequency bands of capable wavelength are focused

Fig. 17 a
Fig. 17 a Cross-section of model AP1 with flat surface.b Cross-section of homogeneous embankment model called model AP2.ρ: density of the model.Vp: P-wave velocity of the model.Vs: S-wave velocity structure of the model.Q: Q-value of the model.c Synthetic vertical velocities recorded along linear array (1) in Fig. 2a for model AP1.d Synthetic vertical velocities recorded along linear array (1) in Fig. 2a for model AP2.N: Number of grid points per shortest wavelength