Study on seismic response characteristics of living stumps slope based on large-scale shaking table test | Scientific Reports

Scientific Reports volume 14, Article number: 23748 (2024) Cite this article

75 Accesses

Metrics details

A large-scale shaking table test of a living stump slope with a geometric similarity ratio of 1:7 was designed and completed. The peak acceleration, acceleration amplification factor, and displacement response patterns of living stumps slopes under different types of seismic waves and excitation intensities were obtained. The time-frequency and energy variation characteristics were analyzed using the Hilbert-Huang Transform (HHT). The results showed that: (1) Regardless of the type of seismic wave, the peak acceleration and acceleration amplification factor of the living stumps slope surface are positively correlated with relative height and seismic excitation intensity. When the excitation intensity is ≤ 0.4 g, the acceleration amplification effect is more pronounced; when the excitation intensity is > 0.4 g, the acceleration amplification effect weakens. (2) Under the action of different seismic waves, the peak displacement of slope surface shows amplification effect along the elevation, and increases with the increase of excitation intensity. In addition, the incremental displacement gradually decreases from the toe to the top of the slope, which is expressed as D2 > D3 > D1 > D4 > D5. The peak displacement at the top of the slope is the greatest, but the incremental displacement is the smallest; the peak displacement at the toe of the slope is the smallest, but the incremental displacement is relatively large. (3) Regardless of the type of seismic wave, living stumps slope shows the characteristics of filtering the low-frequency components of the seismic waves and amplifying their high-frequency components. At the same time, the seismic Hilbert energy gradually accumulates along the elevation. PSHEA and PMSA significantly increase with elevation and excitation intensity, and they reach the maximum at the top of the slope. (4) The seismic Hilbert energy is positively correlated with the relative height and excitation intensity, and reaches the maximum at the top of the slope. With the accumulation of seismic Hilbert energy increases, the dynamic response parameters such as peak acceleration, acceleration amplification coefficient and displacement also increase synchronously, reaching the maximum at the top of the slope. The research conclusions can provide an experimental basis for the seismic design of living stumps slopes.

In plant ecological slope protection technology, plant root systems have shallow root reinforcement and deep root anchoring effects, which can effectively enhance the stability of the slope1,2, and have a good reinforcement effect on shallow landslides (sliding surface depth of 1–2 m). In particular, they can reduce soil moisture content and pore water pressure through plant absorption and transpiration, thereby increasing soil shear strength. Moreover, this technology is highly sustainable. The longer the time is, the more favorable it is to control the shallow sliding and collapse of the slope. It is an environmentally friendly protection technology with “vitality”3. Because the root systems of most slope protection plants are relatively shallow, the current plant slope protection technologies are mainly used for shallow landslide prevention4. Therefore, the research team proposed the living stumps slope support structure, which involves planting vigorous seedlings (such as willows or elms) into deep soil layers or inserting their branches into deep soil layers. Over time, the roots grow outward along the woody stems, forming a root-soil composite in the deep soil layer, thus playing a role in slope reinforcement, as shown in Fig. 1. Living stumps are an environmentally friendly form of slope support, which has the potential to prevent and control deep landslides. The existing research results5,6,7,8 show that the safety factor of slopes supported by living stumps is 30–50% higher than that of original slopes under static conditions. It is highly effective for supporting slopes with slip surfaces within 5 m depth and has broad engineering application prospects.

Living stumps for slope stabilization.

It is well known that earthquakes are a major trigger for slope instability. The existence of living stumps will change the waveform of site seismic waves and the propagation characteristics within the slope’s geotechnical body, making the seismic interaction mechanism between living stumps and the geotechnical body very complicated. In order to study the seismic dynamic response of living stumps and geotechnical bodies, R. Sonnenberg9 used linden wood to make rigid roots and rubber to make flexible roots, and these model roots were used to reinforce clay slopes. Through a series of centrifuge model tests, they explored the trends of axial strain and bending strain in the model roots under seismic conditions, and the failure mechanism of clay slopes reinforced by model roots was revealed. Liang10 used 3D printing technology to create model roots and conducted centrifuge model tests to study the response of vegetated slopes during earthquakes. It was found that plant roots reduced the settlement of the slope top by 67% during the sliding process under the action of earthquakes. The centrifuge model test system of plant slope has obvious limitations, such as only generating force in one direction, unable to reproduce complex seismic wave morphology, and difficult to capture the transient dynamic response of slope under seismic action. In contrast, the shaking table model test can accurately reproduce the seismic waves and record the transient response of the slope under seismic action through the sensor. Additionally, domestic and foreign scholars often use Fourier transform11 and wavelet transform12 to analyze seismic signals. Fourier transform is suitable for stationary signals, but it can only obtain the frequency domain information of seismic signals and cannot provide time domain information. Although wavelet transform can obtain time-frequency domain information, its analysis effect depends on the selection of wavelet basis functions, making it difficult to accurately analyze seismic signals13. In contrast, Hilbert-Huang transform (HHT) has the advantages of strong adaptability, strong ability to deal with nonlinear and non-stationary signals, and high time-frequency resolution14. In this paper, a large-scale shaking table model test of living stumps slopes was designed and conducted. With the help of Hilbert-Huang transform, the seismic response characteristics of living stumps slopes are studied from the perspective of time-frequency domain and energy domain, aiming to provide a basis for related research.

The shaking table model test was completed on the multifunctional shaking table system at the National Engineering Laboratory for High-Speed Railway Construction Technology of Central South University, as shown in Fig. 2. This shaking table system consists of a 4 m×4 m six-degree-of-freedom fixed platform and three six-degree-of-freedom mobile platforms of the same size. The key performance parameters are shown in Table 1.

Shaking table test.

The reliability of shaking table model tests lies in whether the model can accurately reflect the dynamic response of the prototype structure. Therefore, the geometric similarity, kinematic similarity, and material similarity between the prototype and the model15,16 are crucial. The geometric size, acceleration and density similarity ratio are taken as the main control parameters, and the geometric similarity ratio of the model is determined to be 1 : 7, the acceleration similarity ratio is 1 : 1, and the density similarity ratio is 1 : 1. The similarity ratios of the remaining parameters are derived by Meymand’s similarity laws17, as shown in Table 2.

The rigid model box was used in the test, and the internal clearance size was 350 cm\(\:\times\:\)150cm\(\:\times\:\)210 cm (length\(\:\times\:\)width\(\:\times\:\)height). In order to reduce the “model box effect”18, the boundaries of the model box were treated accordingly19: (1) A layer of thick fine stones was laid at the bottom, followed by a layer of fine sand to fill the voids. Then the cement mortar was poured on the top to form the friction boundary, and a layer of 30 cm thick composite material was placed on the top to simulate the bedrock. (2) A 5 cm thick polyethylene foam board was pasted on the inner walls of the model box to create a flexible boundary to weaken the reflection of seismic waves by the model box. (3) Polyvinyl chloride film was pasted on the foam board to reduce the friction between the model and the boundary. The boundary treatment is shown in Fig. 3.

Boundary processing. (a) Fine stone layer at the bottom of the box. (b) Sliding boundary and flexible boundary.

The stability of living stumps slope is closely related to the morphology of the living stumps’ root system. Yen20 and Gray21 classified root system morphology into six categories, among which the “VH” type has a well-developed main root growing vertically downward, with lateral roots concentrated in the upper third of the main root and extending widely. As shown in Fig. 1, for living stumps slope, the length of the main root is a crucial factor affecting the slope’s anti-sliding capacity, while the direction and lateral extension width of the lateral roots are key factors influencing the “anchor support effect” of the lateral roots7. In view of this, the “VH” type elm root system was selected to make a living stump mode in the experiment.

Although the real plant root systems can more accurately reflect their mechanical and hydrological effects, they are not repeatable in model tests. Studies have shown that acrylonitrile-butadiene-styrene (ABS) material can not only effectively simulate the stiffness and strength of roots, but also facilitate the accurate construction of root models in proportion by 3D printing technology10,22,23. The experiment used ABS material and 3D printing technology to create elm stump models at a geometric similarity ratio of 1:7. Epoxy resin was used as an adhesive to attach a layer of quartz sand to the surface of the living stumps model to simulate the real contact surface between roots and soil. In addition, some studies have shown that rectangular distribution is the optimal pattern for enhancing the stability of living stumps slope24. Therefore, in this experiment, the living stumps were arranged in a rectangular grid of 3 rows and 3 columns inside the slope. The center of the living tree pile matrix was located on the central axis of the slope. The horizontal spacing between each living stump was 38 cm, and the vertical spacing was 30 cm. The living stumps model is shown in Fig. 4.

3D printed model of a live stump (Unit: cm).

The soil sample used in the experiment was cohesive soil, and the large stone is filtered out by a sieve with a diameter of 1 cm to simulate the main part of the slope. A composite material with a thickness of 30 cm was used to simulate the bedrock at the bottom of the slope. The composition ratio of the composite material was barite powder: quartz sand: lithium-based lubricant = 10:5:1. The specific material parameters are shown in Table 3. The model box test of the living stumps slope is shown in Fig. 5.

Living stumps slope test model diagram.

Acceleration was measured using uniaxial accelerometers model 1221 L-002, with a range of ± 20 m/s² and a sensitivity of 2000mv/g. Displacement was measured using laser displacement sensors model IL1000, with a range of ± 300 mm. The arrangement of the accelerometer and displacement sensor measurement points is shown in Fig. 6. The accelerometer measurement points were numbered A0-A5 (A0 was located on the shaking table surface to verify seismic wave input). The slope surface displacement measurement points coincide with the accelerometer points and were numbered D1-D5. Considering the bidirectional loading of seismic waves, sensors monitoring dynamic responses in both the horizontal (X-axis) and vertical (Z-axis) directions were arranged at each measurement point. All measurement points were arranged along the central axis of the living stumps slope model.

Layout of accelerometer and displacement measurement points (Unit: cm).

The Kobe wave and the Darui artificial wave with a time compression ratio of 2.65 were used as input waves, both of which were bidirectional excited (horizontal + vertical), hereinafter referred to as K-XZ and DR-XZ. Figures 7 and 8 show the time history curves and corresponding Fourier spectrum when the horizontal acceleration is 0.4 g after time compression of these two kinds of seismic waves. The test follows the specification25, and the horizontal acceleration peaks are set to be 0.1 g, 0.2 g, 0.4 g and 0.6 g respectively. The corresponding vertical acceleration peaks are two-thirds of the horizontal peak, namely 0.067 g, 0.133 g, 0.267 g and 0.4 g. The input wave loading sequence is step-by-step loading. White noise is used to scan before and after the test and before the input wave of the same excitation intensity is loaded. The specific loading scheme is shown in Table 4.

Acceleration time-history curve (0.4 g). (a) K-XZ. (b) DR-XZ.

Fourier spectrum curves (0.4 g). (a) K-XZ. (b) DR-XZ.

Before formally loading each level of seismic waves, a segment of white noise with amplitude of 0.08 g and a duration of no less than 60 s26 (denoted as WN-1 to WN-5) was input to observe the dynamic characteristics of the model in real-time. Figure 9 shows the horizontal acceleration time history curve of the measuring point A5 of the living tree pile slope after the white noise WN-1 and WN-5 are input into the model. From the comparison of Fig. 9 (a) and (b), it can be seen that under the action of white noise with the same excitation intensity, the horizontal acceleration time history curve of the measuring point A5 before and after the test begins to change little. In addition, the data of the measuring point A5 after the action of WN-1 and WN-5 are changed by fast Fourier transform, and the acceleration Fourier spectrum is obtained, as shown in Fig. 10. Figure 10 (a) and (b) are the Fourier spectrum of the measuring point A5 after WN-1 and WN-5, respectively. The comparison shows that the predominant frequency bands of the two Fourier spectra are 0 ~ 30 Hz, and their morphological distribution is similar, with multiple peaks. Therefore, the internal structural damage of the model can be ignored.

Presents the time history curve of horizontal acceleration at point A5 after white noise scanning at different stages. (a) WN-1. (b) WN-5.

Fourier spectrum of A5 after white noise scanning at different stages. (a) WN-1. (b) WN-5.

Some scholars27,28 proposed that the data measured in shaking table tests is susceptible to noise and other interferences, resulting in certain errors. Therefore, before analysis, the wavelet transform filter is used as a low-pass filters to process the data. The data collected by various sensors is analyzed to study the peak acceleration, acceleration amplification factor, and displacement of living stumps slopes under different types and intensities of seismic waves.

The seismic waves were bidirectionally loaded, causing acceleration in both the horizontal and vertical directions of the living stumps slope. However, the acceleration response in the horizontal direction is significantly higher than that in the vertical direction, so the analysis of the horizontal peak acceleration is focused. The horizontal peak acceleration of each measuring point of the living tree pile slope is shown in Fig. 11. According to Fig. 11 (a), under the 0.1 g excitation intensity of K-XZ, the horizontal peak acceleration values increase from 1.057 at the toe of the slope to 2.875 at the top of the slope. Under the other three excitation intensities, the horizontal peak acceleration also shows a gradually increasing trend from the toe to the top of the slope. In addition, as the excitation intensity increases from 0.1 g to 0.6 g, the horizontal peak acceleration of each measuring point increases significantly. For example, the horizontal peak acceleration of the slope top measuring point increases from 2.875 to 16.008. This indicates that the horizontal peak acceleration of the living stumps slope increases with elevation, and reaches the maximum at the top of the slope; at the same time, the horizontal peak acceleration is significantly positively correlated with the excitation intensity.

As shown in Fig. 11 (a), the horizontal peak acceleration value increases from 3.638 to 11.347 under the action of K-XZ with 0.4 g excitation intensity. Compared with the excitation intensity of 0.1 g, when the excitation intensity is 0.4 g, the horizontal peak acceleration increases more along the elevation. It is found that when the excitation intensity is 0.2 g, the increase of horizontal peak acceleration is similar to that of 0.1 g; the increase under 0.6 g excitation intensity is similar to that under 0.4 g excitation intensity. The above results show that under the action of seismic waves, the increase of horizontal peak acceleration along the elevation is positively correlated with the excitation intensity. It is noteworthy that when the excitation intensity is low (< 0.4 g), the increase of horizontal peak acceleration along the elevation is small, while at higher excitation intensities (≥ 0.4 g), the increase is significantly larger. This indicates that there is a positive correlation between the amplification effect of horizontal peak acceleration along the elevation and the excitation intensity.

It can be seen from Fig. 11 (b) that under the action of DR-XZ, the horizontal peak acceleration of the living stumps slope increases with the increase of the slope height, and is positively correlated with the excitation intensity, reaching the maximum at the top of the slope. Compared with Fig. 11 (a), the trend of peak acceleration under DR-XZ is similar to that under K-XZ, but the horizontal peak acceleration value under DR-XZ is greater than that those under K-XZ. The above results indicate that the horizontal peak acceleration response patterns of the living stumps slope are mainly influenced by elevation and excitation intensity, and has little correlation with the type of seismic wave.

Horizontal Peak Acceleration of a Living Stumps Slope under Different Earthquake Waves and Excitation Intensities. (a) K-XZ. (b) DR-XZ.

The ratio of the peak acceleration of each measuring point of the living stumps slope to the peak acceleration of the measuring point (A0) in the same direction is defined as the acceleration amplification factor in this direction. Under different types of seismic waves and different excitation intensities, the horizontal acceleration amplification factors of the living stumps slope are shown in Fig. 12. It can be seen from Fig. 12 that with the increase of slope height, the horizontal acceleration amplification coefficient of the living stumps slope shows the characteristics of nonlinear growth, and there is a significant “acceleration amplification effect”, and the peak value appears at the top of the slope. Furthermore, under different excitation intensities of K-XZ and DR-XZ, the variation patterns of the horizontal acceleration amplification factors are consistent, both of which show an upward trend. The horizontal acceleration amplification factors range from 1.075 to 3.440 under K-XZ excitation, and from 1.005 to 3.735 under DR-XZ excitation.

As shown in Fig. 12 (a), under the action of K-XZ with different excitation intensities, the acceleration amplification factor at the top of the slope is the largest when the excitation intensity is 0.4 g, followed by 0.2 g and 0.6 g, and the acceleration amplification factor under the action of 0.1 g is the smallest. This phenomenon is not only observed at the top of the slope, but also at other measurement points on the slope. It is worth noting that similar results are also observed under the action of DR-XZ with different excitation intensities, which indicates that the acceleration amplification effect has little correlation with the type of seismic wave. Based on the above analysis, the reason for this phenomenon may be that when the seismic wave excitation intensity is low (≤ 0.4 g), the soil of the living stumps slope is re-compaction under vibration. Additionally, due to the low excitation intensity, the internal deformation of the slope is minimal, resulting in less dissipation of seismic wave energy. Consequently, the horizontal acceleration amplification factor of the measurement points increases with the increase of excitation intensity, meaning the amplification effect strengthens with the rise in excitation intensity. However, when the intensity of the seismic wave is high (> 0.4 g), the internal deformation of the slope is large, and the deformation gradually recovers. In this process, the energy dissipation of the seismic wave increases significantly. This leads to a decrease in the horizontal acceleration amplification factor of the measurement points, thereby weakening the horizontal acceleration amplification effect.

Horizontal Acceleration Amplification Factor of a Live Tree Stump Slope under Different Earthquake Waves and Excitation Intensities. (a) K-XZ. (b) DR-XZ.

In the process of seismic wave excitation, the living stumps slope generates horizontal dynamic displacement on the slope surface due to seismic dynamic action. The horizontal peak displacement response of the living stumps slope under different excitation intensities of K-XZ and DR-XZ is shown in Fig. 13. As shown in Fig. 13 (a), under the action of K-XZ with an excitation intensity of 0.4 g, the horizontal peak displacement of the living stumps slope increases from 0.351 mm at the toe to 3.449 mm at the top. Under the other three excitation intensities, the horizontal peak displacement also shows a gradually increasing trend from the toe to the top of the slope. This indicates that the horizontal peak displacement of the living stumps slope shows significant elevation amplification effect. And the dynamic response of the living stumps slope gradually increases from bottom to top, reaching the maximum at the top of the slope. In addition, with the increase of excitation intensity, the horizontal peak displacement increases significantly; and the greater the excitation intensity, the greater the increase of the horizontal peak displacement. For example, when the excitation intensity increases from 0.1 g to 0.6 g, the peak displacement at the top measurement point increases from 0.927 mm to 4.942 mm. As shown in Fig. 13(b), under different excitation intensities of DR-XZ, the horizontal peak displacement response patterns of the living stumps slope are similar to those in Fig. 13(a). This shows that the horizontal peak displacement response of the living stumps slope mainly depends on the elevation and the excitation intensity, and is positively correlated, while it is less correlated with the type of seismic wave.

As shown in Fig. 13, regardless of the type of seismic wave, when the excitation intensity is less than 0.4 g, the horizontal peak displacement of the living stumps slope increases slightly along the elevation; whereas, when the excitation intensity is ≥ 0.4 g, the increase is large. Specifically, under the action of 0.1 g DR-XZ, the horizontal peak displacement increases from 0.347 mm to 0.970 mm along the elevation; while under 0.6 g, it jumps from 0.561 mm to 5.468 mm. This shows that regardless of the type of seismic wave, the increase of horizontal peak displacement along the elevation is mainly related to the excitation intensity.

Horizontal Dynamic Peak Displacement of a Live Tree Stump Slope under Different Earthquake Waves and Excitation Intensities. (a) K-XZ. (b) DR-XZ.

Considering the difference in the increase of seismic wave excitation intensity in the model test and the lag effect of propagation in the model, the incremental displacement29 is adopted, and it is normalized by the slope height, as shown in Fig. 14. It can be seen from Fig. 14 that no matter what kind of seismic wave, when the excitation intensity is 0.6, the incremental displacement of each measuring point is greater than the value under the other three excitation intensities. Further analysis of Fig. 14, it can be found that no matter what kind of seismic wave action, the incremental displacement of each measuring point shows D2 > D3 > D1 > D4 > D5. It is noteworthy that although the horizontal peak displacement of D5 in Fig. 11 is significantly higher than that of the other four measurement points, the residual displacement caused by the vibration is larger, resulting in the smallest incremental displacement. In contrast, D1 and D2 located at the toe of the slope have a large incremental displacement due to the small residual displacement after the end of the vibration, which indicates that the toe of the slope has good stability and self-healing ability, and can recover most of the peak displacement.

Incremental Displacement of a Live Tree Stump Slope under Different Earthquake Waves and Excitation Intensities. (a) K-XZ. (b) DR-XZ.

The Hilbert-Huang Transform30,31 (hereinafter referred to as HHT) is an advanced tool for analyzing nonlinear and non-stationary time series signals. HHT can effectively capture the spectral characteristics of slopes under seismic excitation, so it is often used as a vibration signal analysis tool. Compared with Fourier transform and wavelet transform, HHT has stronger adaptability, higher time-frequency resolution and complete data-driven characteristics when dealing with nonlinear and non-stationary signals. At the same time, it can reduce the error of Fourier transform and effectively avoid the time-frequency relationship in Heisenberg uncertainty principle. Therefore, HHT has a high adaptability to sudden signals such as earthquakes. By applying HHT to the time acceleration data of the living stumps slope, the dynamic response of the living stumps slope can be analyzed in both the time and frequency domains simultaneously. Additionally, the Hilbert marginal spectrum can be obtained by integrating the HHT spectrum over time. The area under the marginal spectrum represents the total seismic Hilbert energy in the total frequency range, and provides the total seismic Hilbert energy contribution scale of each frequency value. Therefore, the marginal spectrum change can be obtained by the acceleration of each measuring point, and the energy change law of different positions in the process of seismic excitation can be compared.

The Hilbert-Huang transform (HHT) is used to process the acceleration data of each measuring point at the foot, middle and top of the slope of the living stumps slope under the action of two kinds of seismic waves with different excitation intensities, and the corresponding HHT spectrum is obtained. Figures 15 and 16 show the HHT spectra under the action of K-XZ and DR-XZ when the excitation intensity is 0.2 g and 0.4 g. It can be seen from the figure that the main frequency of the HHT spectrum of the living tree pile slope is distributed between 5 and 15 Hz regardless of the type of seismic wave. Figures 15 and 16 show that regardless of the type of seismic wave, the low-frequency component changes to the high-frequency component during the propagation of the seismic wave to the top of the slope, and the seismic Hilbert energy gradually converges at the top of the slope. This indicates that the living stumps slope has the characteristic of filtering the low-frequency components and amplifying the high-frequency components, and this behavior is less correlated with the type of seismic wave.

Figure 15 (b) shows the HHT spectrum of the living stumps slope under the action of K-XZ with the excitation intensity of 0.4 g, where the PSHEA at the measurement points from the toe to the top of the slope are 2.9, 8.3, and 12.1, respectively. Figure 15 (a) also reveals the gradual increase of PSHEA from bottom to top. This indicates that the PSHEA of the living stumps slope increases continuously along the elevation, and this increasing trend may be related to the filtering of low-frequency components and the amplification of high-frequency components of the living tree pile slope. Comparing Fig. 15 (a) and Fig. 15 (b), when the excitation intensity increases from 0.2 g to 0.4 g, the PSHEA at the top of the slope increases from 6.3 to 12.1. At the same time, the PSHEA at the toe and middle of the slope also increases. This result strongly proves that there is a significant positive correlation between PSHEA and excitation intensity.

HHT spectrum of live tree stump slopes under K-XZ effects at different excitation intensities. (a) 0.2 g. (b) 0.4 g.

Figure 16 shows the HHT spectrum of the living stumps slope under the action of DR-XZ with excitation intensity of 0.2 g and 0.4 g. The analysis results show that the change trend of PSHEA under the action of DR-XZ with elevation and excitation intensity is consistent with that under the action of K-XZ, that is, it is positively correlated with elevation, but there are numerical differences between the two. This shows that PSHEA is mainly affected by elevation and excitation intensity, and has little correlation with seismic wave types.

HHT spectrum of live tree stump slopes under DR-XZ effects at different excitation intensities. (a) 0.2 g. (b) 0.4 g.

Through the integration of HHT spectrum with time, the marginal spectrum is obtained and analyzed to explore the variation patterns of seismic Hilbert energy in the living stumps slope. Due to the limitation of space, this paper only presents the marginal spectrum of the slope toe, middle and top of the living stumps slope under the action of K-XZ and DR-XZ with excitation intensity of 0.2 g and 0.4 g (see Figs. 17 and 18).

Figure 17 is the marginal spectrum of K-XZ with excitation intensities of 0.2 g and 0.4 g. According to Fig. 17 (a), when the excitation intensity is 0.2 g, the marginal spectral peak (PMSA) at the toe of the slope is mainly concentrated in the range of 0 ~ 5 Hz.As the elevation increases, the frequency corresponding to PMSA migrates to the vicinity of 10 Hz. Meanwhile, with the increase of slope height, the PMSA from slope toe to slope top is 0.139, 0.574 and 0.814 respectively. In Fig. 17 (b), the value of PMSA also increases from the foot of the slope to the top of the slope. In addition, when the excitation intensity increases from 0.2 g to 0.4 g, the PMSA of the living stumps slope increases significantly, for example, the PMSA at the top of the slope increases from 0.814 to 1.399. This indicates that the PMSA of the living stumps slope is positively correlated with the elevation and increases with the increase of the excitation intensity.

Figure 18 shows the marginal spectrum under the action of DR-XZ with excitation intensity of 0.2 g and 0.4 g. By comparing with Fig. 17, it can be found that the PMSA variation patterns of the living stumps slope under the action of DR-XZ are consistent with that under the action of K-XZ. This indicates that the variation pattern of PMSA is less correlated with the type of seismic wave, which is mainly determined by elevation and excitation intensity. It is noteworthy that the area under the marginal spectrum reflects the total seismic Hilbert energy over the entire frequency range32. As shown in Figs. 17 and 18, with the increase of slope height and excitation intensity, the seismic Hilbert energy gradually accumulates and reaches the maximum at the top of the slope. At the same time, the increase of elevation significantly enhances the filtering effect of the living stumps slope on the low-frequency component of the earthquake and the amplification effect of the high-frequency component, so that the seismic Hilbert energy is transferred from low frequency to high frequency with elevation, resulting in the most significant dynamic response of the slope top. Comparing Figs. 17 and 18, the type of seismic wave has a minor impact on the variation trend of seismic Hilbert energy, only altering its numerical magnitude. In addition, the seismic Hilbert energy is mainly affected by elevation and excitation intensity, and with the increase of elevation and excitation intensity, the seismic Hilbert energy also increases.

Marginal Spectrum of Live Stump Slope under K-XZ Effect at Different Excitation Intensities. (a) 0.2 g. (b) 0.4 g.

Marginal Spectrum of Live Stump Slope under DR-XZ Effect at Different Excitation Intensities. (a) 0.2 g. (b) 0.4 g.

Numerous scholars33,34,35 have revealed a series of general patterns of seismic dynamic response of slopes through large-scale shaking table tests. Specifically, the peak acceleration, its amplification coefficient and the peak displacement are positively correlated with elevation and excitation intensity. The dynamic response law of the living stumps slope observed in this study is consistent with the research conclusions of many scholars, indicating that the existence of the living stumps does not significantly change the dynamic response law of the slope. However, the existing research on slope displacement mostly focuses on peak displacement, ignoring the gradual increase of seismic wave intensity and its hysteresis effect in the model. Therefore, this paper introduces incremental displacement, as shown in Fig. 14. In this paper, the incremental displacement is analyzed, which effectively supplements the research on slope displacement. The incremental displacement and the peak displacement are mutually verified, which further enhances the reliability of the conclusion.

This paper not only explores the seismic dynamic response patterns of the living stumps slope from the time domain, but also obtains the HHT spectra and marginal spectra based on Hilbert-Huang Transform (HHT), and studies time-frequency and energy variation characteristics of the living stumps slope. The research verifies the patterns of the dynamic response of the slope surface acceleration, acceleration amplification coefficient and peak displacement of the living stumps slope observed in the paper, and further proves the influence of the seismic Hilbert energy on the dynamic response of the living stumps slope during the propagation process. The results indicate that as the seismic waves propagate along the slope height, the seismic Hilbert energy continuously accumulates and reaches its peak at the top of the slope. At this time, various parameters representing the seismic dynamic response of the living stumps, such as peak acceleration, show a positive correlation with the elevation, which is most significant at the top of the slope. At the same time, with the increase of seismic excitation intensity, the seismic Hilbert energy also increases, and the dynamic response parameters such as peak acceleration, acceleration amplification factor and displacement also increase synchronously in value. This further illustrates that the dynamic response parameters of the living stumps slope are closely related to the seismic Hilbert energy, indicating that dynamic response parameters such as acceleration are a manifestation of the seismic Hilbert energy.

In this paper, a large-scale shaking table test of the living stumps slope with a geometric similarity ratio of 1:7 was designed and completed. The peak acceleration, acceleration amplification factor and displacement response of the living stumps slope under the action of K-XZ and DR-XZ with different excitation intensities were obtained. Based on the Hilbert-Huang transform, the HHT spectrum and marginal spectrum of the living stumps slope were acquired. The time-frequency and energy variation characteristics were studied, as well as the relationships between seismic Hilbert energy and dynamic response parameters such as peak acceleration, acceleration amplification factor, and displacement. The main conclusions are as follows:

(1) Whether under the action of K-XZ or DR-XZ, the peak acceleration and the acceleration amplification coefficient of the living stumps slope are significantly positively correlated with the relative height and excitation intensity. When the excitation intensity is ≤ 0. 4 g, the acceleration amplification coefficient increases significantly with the increase of excitation intensity, and the acceleration amplification effect is significant. When the excitation intensity is >0.4 g, the acceleration amplification coefficient is less than the value when the excitation intensity is low, and the amplification effect is weakened.

(2) Under the action of different seismic waves, the peak displacement of the soil of the living stumps slope is amplified along the elevation and increases with the increase of excitation intensity. When the excitation intensity is < 0.4 g, the growth rate of the peak displacement along the slope height is relatively gentle. However, when the excitation intensity is ≥ 0.4 g, the growth rate significantly accelerates, and the amplification effect is notably enhanced. In addition, the incremental displacement of the living stumps slope follows the pattern D2 > D3 > D1 > D4 > D5. And the peak displacement is the largest at the top of the slope, but the incremental displacement is the smallest. Conversely, the peak displacement is the smallest at the toe of the slope, but the incremental displacement is relatively larger.

(3) During the propagation of seismic waves, the low-frequency component changes to the high-frequency component, and the living stumps slope shows the characteristics of filtering the low-frequency component of seismic wave and amplifying its high-frequency component.

(4) Under the action of K-XZ or DR-XZ, PSHEA and PMSA increase significantly with elevation and excitation intensity, and reach the maximum at the top of the slope. At the same time, the marginal spectrum is further studied and analyzed. The seismic Hilbert energy shifts from low frequency to high frequency with the increase of height, and gradually accumulates with the slope height, reaching the maximum at the top of the slope. The values of the dynamic response parameters of the living stumps slope are closely related to the seismic Hilbert energy, that is, with the upward accumulation of the seismic Hilbert energy, the values of the dynamic response parameters such as the peak acceleration, the acceleration amplification factor and the displacement of the living stumps slope increase significantly.

The datasets used and analyzed during the current study are available from the corresponding author upon rea-sonable request.

Zhou, D. & Zhang, J. 1 Vegetation Slope Protection Engineering Technology. (People’s Transportation Publishing House, (2003).

Donn, S., Wheatley, R. E., Mckenzie, B. M., Loades, K. W. & Hallett, P. D. Improved soil fertility from compost amendment increases root growth and reinforcement of surface soil on slopes. Ecol. Eng. 71, 458–465. https://doi.org/10.1016/j.ecoleng.2014.07.066 (2014).

Article Google Scholar

Burri, K., Graf, F. & Böll, A. Revegetation measures improve soil aggregate stability: a case study of a landslide area in Central Switzerland. For. Snow Landsc. Res. 82, 45–60 (2009).

Google Scholar

Prasad, A., Kazemian, S., Kalantari, B., Huat, B. B. K. & Mafian, S. Stability of Tropical residual soil Slope Reinforced by Live Pole: Experimental and Numerical investigations. Arab. J. Sci. Eng. 37, 601–618. https://doi.org/10.1007/s13369-012-0209-2 (2012).

Article Google Scholar

Xueliang, J., Wenjie, L., Hui, Y., Haodong, W. & Zhenyu, L. Study on Mechanical Characteristics of Living Stumps and Reinforcement Mechanisms of Slopes. Sustainability. 16, 4294. https://doi.org/10.3390/su16104294 (2024).

Article Google Scholar

Hui, Y. et al. Study on the supporting effect of Bamboo Anchor and wood Frame Beam reinforcing cohesive soil slope. Indian Geotech. J. 53, 844–858. https://doi.org/10.1007/s40098-022-00702-3 (2023).

Article Google Scholar

Xueliang, J., Wang, Z., Yang, H., Guo, J. & Wang, H. Dynamic Response and Stability of Slopes with Living Stumps under Train Load. Journal of Vibration Engineering 1–10, (2024). https://link.cnki.net/urlid/32.1349.tb.20231222.20231052.20231004

Xueliang, J. et al. A 3D Model Applied to Analyze the Mechanical Characteristic of Living Stump Slope with Different Tap Root Lengths. Applied Sciences 13, (1978). (2023) https://doi.org/10.3390/app13031978

Sonnenberg, R., Bransby, M. F., Bengough, A. G., Hallett, P. D. & Davies, M. C. R. Centrifuge modelling of soil slopes containing model plant roots. Revue Canadienne De Géotechnique. 49, 1–17. https://doi.org/10.1139/T11-081 (2012).

Article Google Scholar

Liang, T., Knappett, J. A., Bengough, A. G. & Ke, Y. X. Small-scale modelling of plant root systems using 3D printing, with applications to investigate the role of vegetation on earthquake-induced landslides. LANDSLIDES 14(5), 1747–1765 (2017). (2017). https://doi.org/10.1007/s10346-017-0802-2

Wang, K. L. & Lin, M. L. Initiation and displacement of landslide induced by earthquake—a study of shaking table model slope test. Eng. Geol. 122, 106–114. https://doi.org/10.1016/j.enggeo.2011.04.008 (2011).

Article Google Scholar

Wen, S., Qingguo, L., Xiangjin, Q., Xiaoping, C. & Lili, W. Study on dynamic response of Loess Slopes with different instability patterns. J. China Railway Soc. 44, 123–130. https://doi.org/10.3969/j.issn.1001-8360.2022.06.015 (2022).

Article Google Scholar

Pai, L. & Wu, H. Experimental study on the dynamic response of lining structures in Tunnels Orthogonally Underpassing landslides under Seismic Action. Chin. J. Rock Mechan. Eng. 41, 979–994. https://doi.org/10.13722/j.cnki.jrme.2021.0542 (2022).

Article Google Scholar

Song, D., Liu, X., Huang, J. & Zhang, J. Energy-based analysis of seismic failure mechanism of a rock slope with discontinuities using Hilbert-Huang transform and marginal spectrum in the time-frequency domain. Landslides. 18, 105–123. https://doi.org/10.1007/s10346-020-01491-7 (2021).

Article Google Scholar

I. A. I. Similitude for shaking table tests on soil-structure-fluid model in 1 g gravitational field. Rep. Port Horbour Res. Inst. 27, 3–24. https://doi.org/10.3208/sandf1972.29.105 (1989).

Article Google Scholar

Xueliang, J., Jiayong, N., Pengyuan, L., Changping, W. & Feifei, W. Large-scale shaking table test study on seismic dynamic characteristics of Rock Slopes with Small Clear-Distance tunnels. Eng. Mech. 34, 11. https://doi.org/10.6052/j.issn.1000-4750.2015.11.0936 (2017).

Article Google Scholar

Meymand, P. J. Shaking Table Scale Model Tests of Nonlinear soil-pile-superstructure Interaction in soft clay (University of California, 1998).

Jianbo, D. & Chengtao, H. Design and testing of a bidirectional layered shear-type Continuum Model Box for shaking table tests. J. Vib. Shock. 42, 164–171. https://doi.org/10.13465/j.cnki.jvs.2023.18.018 (2023).

Article Google Scholar

Xueliang, J. et al. Dynamic response of shallow-buried tunnels under asymmetrical pressure distributions. J. Test. Evaluation. 46, 1574–1590. https://doi.org/10.1520/JTE20170026 (2018).

Article Google Scholar

Yen, C. in Proceedings of the International Workshop on Soil Erosion and its Conuntermeasures. 92–111 (Soil and Water Conservation Society of Thailand).

Gray, D. H. & Sotir, R. B. Biotechnical stabilization of highway cut slope. J. Geotech. Eng. 118, 1395–1409 (1992). https://doi.org/https://doi.org/10.1061/ (ASCE)0733-9410(1992)118:9(1395).

Article Google Scholar

Liang, T., Knappett, J. & Bengough, A. G. in 8th International Conference on Physical Modelling in Geotechnics 2014. (ed Gaudin and White) 361–366 (International Society for Soil Mechanics and Geotechnical Engineering).

Liang, T. & Knappett, J. A. Centrifuge modelling of the influence of slope height on the seismic performance of rooted slopes. Géotechnique. 67, 855–869. https://doi.org/10.13722/j.cnki.jrme.2021.0542 (2017).

Article Google Scholar

Temgoua, A. G. T., Kokutse, N. K. & Kavazovi, Z. Influence of forest stands and root morphologies on hillslope stability. ECOL ENG 95, 622–634 (2016). (2016). https://doi.org/10.1016/j.ecoleng.2016.06.073

A.C.M.A. & S. Vol. GB 55002 – 2021 (China Architecture Publishing and Media Co., 2021).

Haohua, H. Design and Application Technology of Seismic Simulation Shaking Tables (Seismological, 2008).

Xinhao, T., Jing, L. & Liang, Z. Damage evolution mechanism of rock-soil mass of bedrock and overburden layer slopes based on shaking table test. J. Mt. Sci. 19, 3645–3660. https://doi.org/10.1007/s11629-022-7403-9 (2022).

Article Google Scholar

Xinhao, T., Jing, L., Changwei, Y. & Liang, Z. Shaking table test on dynamic damage characteristics of bedrock and overburden layer slopes. J. Test. Eval. 51, 989–1009. https://doi.org/10.1520/JTE20220314 (2023).

Article Google Scholar

Mao, Y. et al. Dynamic response characteristics of shaking table model tests on the gabion reinforced retaining wall slope under seismic action. Geotext. Geomembr. 52, 167–183. https://doi.org/10.1016/j.geotexmem.2023.10.001 (2024).

Article Google Scholar

Yi, Z., Jun, Y., Feng, X., Jian, C. & S. & Studies of filtering effect on internal solitary wave flow field data in the South China Sea using EMD. Adv. Mater. Res. 518, 1422–1425 (2012). https://doi.org/https://doi.org/10.4028/www.scientific.net/AMR.518-523.1422

Google Scholar

Yi, Z. et al. Current field features and propagation characteristics of suspected internal solitary wave packet. Ocean Eng. 72, 448–452. https://doi.org/10.1016/j.oceaneng.2013.07.018 (2013). https://doi.org/https://doi.org/

Article ADS Google Scholar

Gang, F., Limin, Z., Jian-Jing, Z. & Ouyang, F. Energy-BasedAnalysis of mechanisms of Earthquake- InducedL and slide using Hilbert-Huang Transform and marginal spectrum. Rockmechanics Rock. Eng. 50, 2425r2441. https://doi.org/10.1007/s00603-017-1245-8 (2017).

Article Google Scholar

Changping, W., Xueliang, J., Guolin, Y., Hongbin, X. & Zhongqiu, X. Shaking table test study on seismic dynamic characteristics of gravity retaining walls with secondary support slopes. J. Vib. Eng. 27, 426–432. https://doi.org/10.16385/j.cnki.issn.1004-4523.2014.03.014 (2014).

Article Google Scholar

Zhouzhou, X. et al. Study on Seismic Response Differences of Soil-Rock Mixture Slopes with different rock contents based on shaking table tests. Rock. Soil. Mech., 1–14 https://doi.org/10.16285/j.rsm.2023.1379

Zhang, Z., Wang, T., Wu, S., Tang, H. & Liang, C. Seismic performance of loess-mudstone slope in Tianshui – Centrifuge model tests and numerical analysis. Eng. Geol. 222https://doi.org/10.1016/j.enggeo.2017.04.006 (2017).

Download references

This research were financially supported by the following projects: (1) The Hunan Provincial Natural Science Foundation of China (Grant No.: 2022JJ31005); (2) Tertiary Education Scientific research project of Guangzhou Municipal Education Bureau (Grant No.: 2024312379, 2024312426); (3)Key Fields of Higher Education in Guangdong Province (Grant No.: 2023ZDZX4044, 2023ZDZX4045); (4) Research Capacity Enhancement Project of Key Construction Discipline in Guangdong Province (Grant No.: 2022ZDJS092); (5) The Forest Science and Technology Innovation Program of Hunan Province (Grant No.: XLK202105-3); (6)Project supported by National Natural Science Foundation of China (Grant No.: 31971727).

Guangdong Provincial Key Laboratory of Green Construction and Intelligent Operation & Maintenance for Offshore Infrastructure, Guangzhou Maritime University, Guangzhou, 510725, China

Hui Yang & Xueliang Jiang

school of civil engineering, Central South University of Forestry and Technology, Changsha, 410004, China

Hui Yang, Jun Yin, Xueliang Jiang, Bo Shen, Haodong Wang & Zhenzhen Wei

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

H. Y. was responsible for the investigation, conception, and design of the whole study, undertook the task of data processing and analysis, and participated in the writing and revision of the paper. J. Y. was mainly responsible for the design and operation of the experiment and participated in the writing and revision of the paper and was responsible for deepening and integrating the theoretical framework, providing rich theoretical support for the research. X.l. J. is responsible for funding applications, ethical reviews, and other matters to ensure the smooth progress of the research. B. S. was responsible for drawing and modifying the charts and tables on the paper and participated in the modification of the paper format. H.d. W. focused on polishing the language and stylistic unification to ensure the fluency and professionalism of the articles. Z.z. W. provided strategic and theoretical guidance throughout the research process, especially when encountering key problems or challenges, giving valuable solutions and suggestions. All the authors have reviewed the manuscript.

Correspondence to Jun Yin.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. 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://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

Yang, H., Yin, J., Jiang, X. et al. Study on seismic response characteristics of living stumps slope based on large-scale shaking table test. Sci Rep 14, 23748 (2024). https://doi.org/10.1038/s41598-024-75676-8

Download citation

Received: 14 July 2024

Accepted: 07 October 2024

Published: 10 October 2024

DOI: https://doi.org/10.1038/s41598-024-75676-8

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative