Abstract
Local soil conditions play a crucial role in shaping ground surface responses and impacting the intensity of ground shaking. In this study, the influence of different site profiles on computed ground motion was investigated using a 2D equivalent linear analysis approach. Four distinct site profiles: sand, clay, sand-clay-sand-clay, and clay-sand-clay-sand profiles were considered. The results were presented using multiple metrics, including surface acceleration, displacement, modulus decreasing ratio, and coherence analysis. Notably, the clay profile significantly influenced ground motion, while the sand profile exhibited relatively lower seismic activity. This suggests that softer sites significantly influence ground motion, leading to potentially high levels of shaking.
1 Introduction
Seismic site response analysis is a crucial technique for assessing the impact of local soil profiles on ground motion. The characteristics of local soil conditions including shear wave Velocities (Vs) profiles, soil layering, modulus reduction and damping properties have significant effect on amplitude of ground shaking. The impact of local soil nature on site response is often assessed using a 1D equivalent linear analysis approach. In this approach, the analysis results are commonly interpreted in terms of surface spectral acceleration, amplification factors, and surface acceleration time-histories [1]. The analysis assumes that the vertically propagating horizontally polarized seismic waves affect soil deposits. The non-linear dynamic properties of the deposits (shear modulus and damping behavior) are approximated by an iterative procedures [2]. Several factors contribute to variability in seismic site response analysis, including variations in soil types and layering, changes in dynamic characteristics of soils, analysis approaches and differences in the intensity of seismic waves [3]. In the recent years, researchers have made attempts to incorporate these variations into site response analysis and evaluate their impact on computed ground motion [1, 4–9].
The characteristics of soil conditions have a significant impact on the amplitude of ground shaking during seismic events. Performing in-situ tests to investigate soil deposit characteristics may not be practical or economically feasible, except for a few critical facilities. Hence, it becomes imperative to include a wide range of soil types (e.g., soft clay to stiff sand) in site response analysis, to assess their impact on ground motion. Often, the variability of soil properties is modeled using specific probability distribution functions to capture the uncertainty inherent in soil characteristics. A great deal of previous researches have been conducted to examine how soil variability influence ground motion [10–13]. Bazzurro and Cornell [14] studied the influence of soil layering and uncertainty associated with layering on surface motion intensities. In addition to soil variability, the intensity of input rock motion characteristics may have considerable impact on level of ground shaking [15]. Proper characterization of soil layering and non-linear dynamic properties of soil layers is also important.
The magnitudes of ground motion induced by incoming seismic waves are profoundly influenced by the unique properties of soil strata. In stratified soil deposits, the layering may involve the presence of softer layers (e.g., clay) sandwiched between stiffer layers (e.g., sand), or the alternation of sand layers within clay layers. These heterogeneous soil layers play a significant role in determining the amplitude of ground motion and should be carefully considered in site response analyses. Incorporating the complexities of stratified soil profiles is essential for a comprehensive understanding of how seismic waves interact with stratified soil deposits. The aim of this study is to evaluate the impact of stratified soil layers on amplitude of computed ground motion. For this purpose, four distinct types of soil profiles were evaluated using 2D equivalent linear analysis method. The results of the analyses were presented in terms of surface acceleration and displacement time-histories, shear modulus decreasing ratio and wavelet coherence analysis. The accuracy and reliability of the 2D equivalent linear analysis were cross-checked by comparing representative results obtained from the 2D equivalent linear method with the STRATA (1D) equivalent linear analysis.
2 Site profiles and material properties
In this study, 2D equivalent linear site response analysis was performed using the MIDAS GTS NX commercial software [22]. The 2D equivalent linear seismic site response analysis requires soil profile data, material properties, and input ground motion in the form of acceleration time histories. The study considered hypothetical site profiles, including Sa, Cl, Sand-Clay-Sand-Clay (SaClSaCl), and Clay-Sand-Clay-Sand (ClSaClSa) profiles to evaluate influence of soil variability on computed ground response (Fig. 1). Seismic ground response is greatly affected by the dynamic properties of soils, specifically the normalized shear modulus reduction and damping curves. These curves define the variations of modulus reduction and damping of soils with induced shear strains. The MIDAS GTS NX software allowed the incorporation of different dynamic soil properties for each soil profile. The modulus reduction and damping curves for sands proposed by Seed and Idriss [23] were considered. Additionally, the modulus reduction and damping curves for clays (with plasticity indices of 30 and 50) proposed by Vucetic and Dobry [24] were considered. Figure 2 presents dynamic soil properties of site profiles considered for response analysis. Table 1 shows the unit weight, initial shear modulus, and dynamic soil properties (symbolized with the letter M) considered for the response analysis.
Site profiles considered for response analyses (Source: Authors' result)
Citation: Pollack Periodica 19, 2; 10.1556/606.2024.00995
a) Modulus reduction, b) damping ratio with respect to shear strain (Source: on the basis of [23, 24] plotted by the Authors)
Citation: Pollack Periodica 19, 2; 10.1556/606.2024.00995
Material properties for each profile
Profiles | Unit weight, γ (kN m−3) | Gmax (MPa) | Dynamic soil properties (M) |
Cl | 19 | 20 | M2 |
ClSaClSa | 18, 20, 19, 21 | 15, 40, 20, 50 | M1, M3, M2, M4 |
Sa | 20 | 40 | M3 |
SaClSaCl | 20, 18, 21, 19 | 40, 15, 50, 20 | M2, M1, M4, M3 |
3 Rock motion selection
To perform site response analysis, seismic input loads in the form of acceleration time histories are required. These records are essential for assessing how a site's soil and geological conditions influence the behavior of seismic waves. In most cases, geotechnical engineers select real earthquake records based on the seismic hazard level specific to the site [25]. These records represent the ground motion characteristics that the site may experience during an actual seismic event. Several databases compile real seismic records from monitoring stations worldwide [26, 27]. Synthetic ground motion can also be generated to perform seismic site responses [28].
In this study, records from the 1989 Loma Prieta Earthquake, characterized by magnitude of 6.93 recorded at an epicentral distance of 74 km was chosen from Pacific Earthquake Engineering Research (PEER) strong motion database center [26]. Figure 3 illustrates the earthquake records and its spectral acceleration considered for response simulations. The peak acceleration and time interval of the recorded data were 0.276 g and 0.01 s, respectively (Fig. 3a). The record was applied at the base of each site profile and allowed to propagate through the soil to evaluate responses at the ground level.
The 1989 Loma Prieta Earthquake, a) Acceleration time histories, b) spectral acceleration (Source: on the basis of [26] plotted by the Authors)
Citation: Pollack Periodica 19, 2; 10.1556/606.2024.00995
4 2D equivalent linear seismic site response analysis
The 2D equivalent linear site response analysis was conducted using MIDAS GTS NX software [22], as previously indicated. It employs an equivalent linear analysis approach to simulate dynamic properties of soil under cyclic seismic loads. The soil layer is simplified as layers in a horizontal plane with different properties (e.g., shear modulus and damping). Transmitting boundary conditions available in MIDAS GTS NX library are considered to minimize reflection effects on vertical boundaries [29]. The transmitting boundary conditions typically assume that the horizontal properties of each ground layer are equal. Figure 4 depicts the meshed finite element model representing the ClSaClSa profile, created using MIDAS GTS NX software. Each profile has overall dimensions of 200 m by 60 m. The meshing process involved a rigorous evaluation of various element sizes to assess their impact on the analysis outcomes. After thorough assessment, it was determined that an element size of 2 m yielded optimal results. By selecting this element size, the analysis results remain reliable while maintaining computational efficiency.
2 D finite element model mesh grid of ClSaClSa soil profiles (Source: Authors' result)
Citation: Pollack Periodica 19, 2; 10.1556/606.2024.00995
In the equivalent linear analysis approach, the nonlinear dynamic properties of soil are approximated by equivalent linear models. This involves updating the shear modulus and damping as a function of the shear strain during the seismic event. Typically, the values of modulus reduction and damping are estimated at 65% of maximum shear strain, commonly referred to as effective shear strain [30]. The equivalent linear analysis follows an iterative procedure, as outlined in [2, 10]:
Start with the initial assumption of shear modulus and damping values (commonly low-strain values are considered);
Compute the ground response using these initial values;
Compute the effective shear strain for each layer (65% of maximum shear strain induced within each layer);
Select new shear modulus and damping values based on computed effective shear strain;
Repeat steps 2 to 4 until the difference between estimated values in two successive iterations falls below predetermined threshold value, typically ranging from 5 to 10%.
5 Results and discussions
2D seismic site response analyses were conducted for all site profiles under identical rock motion intensity using MIDAS GTS NX. The simulation results are recorded at the ground surface for each soil profile and are presented in terms of time histories for surface acceleration and displacement. Furthermore, modulus decreasing ratio and wavelet coherence analysis were evaluated. The maximum surface acceleration and displacement were recorded from Cl profile while the minimum peak acceleration and displacement were recorded from Sa profile (Fig. 5). The relative difference between the maximum and minimum acceleration and displacement were 80% and 64%, respectively. The maximum acceleration and displacement computed from SaClSaCl profile are very close to those computed from ClSaClSa profile. Figure 6 display samples of simulation results for peak ground acceleration and displacement obtained from the Cl profile. The maximum ground acceleration and displacements are indicated in the legend (left of each image).
a) Acceleration time histories, b) surface displacement time histories of response analyses (Source: Authors' results)
Citation: Pollack Periodica 19, 2; 10.1556/606.2024.00995
Simulated peak ground a) acceleration (m s−2) and b) displacement (m) for Cl profile (Source: Authors' result)
Citation: Pollack Periodica 19, 2; 10.1556/606.2024.00995
Color gradient was used to visualize the shear modulus decreasing ratio. Figure 7 illustrates the variation of shear modulus decreasing ratio across different layers of the site profiles. Sa profile exhibits lower seismic action as indicated by black color gradient across the entire depth of profile. In contrast, the Cl profile experiences higher seismic action, noticeable through dark to light grey gradient throughout the profile depth. As expected, clay soils behave differently in response to seismic loading compared to sands. Clay soils recorded lower shear modulus decreasing ratio than sands, indicating that clay is less stiff and more prone to deformation under seismic forces (e.g., see SaClSaCl and ClSaClSa profiles). This suggests that softer materials may experience amplified ground motion, potentially resulting in higher levels of shaking.
Color gradient of shear modulus decreasing ratio (Source: Authors' result)
Citation: Pollack Periodica 19, 2; 10.1556/606.2024.00995
Figure 8 shows acceleration time histories (Fig. 8a) and displacement (Fig. 8b) of Sa profile obtained from both 2D and 1D equivalent linear analyses. This comparison serves as a validation of the accuracy and reliability of MIDAS GTS NX software. The 1D equivalent analysis was conducted using the STRATA computer program, a widely adopted software for site response analysis. Notably, the results exhibit excellent agreement, underscoring the robust capabilities of the MIDAS GTS NX software.
Comparison between MIDAS GTS NX and STRATA response results, a) acceleration time histories, and b) displacement (Source: Authors' result)
Citation: Pollack Periodica 19, 2; 10.1556/606.2024.00995
6 Ground response wavelet coherence analysis
Figure 9 illustrates the results of wavelet coherence analysis between surface acceleration of different profiles. The cone of influence is shown as shade under broken line. Regions of varying coherency are depicted with a heat map, ranging from low coherency (light regions) to high coherency (dark regions). The y-axis represents scale or frequencies, and x-axis represents the period of vibration. The arrows represent variations between two time-series. Right-pointing arrows indicate movement in the same direction (co-movement). Downward-pointing arrows indicate lagging of the first index. Left-pointing arrows denote variables that are out-of-phase. The results show a clear strong coherence between accelerations computed from SaClSaCl and ClSaClSa profiles. While some signs of low coherence exist between accelerations computed from Cl and Sa, Cl and ClSaClSa, Cl and SaClSaCl, there is strong overall coherence. Arrows mostly point to the right demonstrating co-movement with high coherence effects present throughout the period of vibration.
Wavelet coherence analysis between acceleration time series computed from different site profiles (Source: Authors' result)
Citation: Pollack Periodica 19, 2; 10.1556/606.2024.00995
7 Conclusion
In this study, the impact of soil variability on seismic site response analyses were analyzed using four distinct site profiles (Cl, Sa, ClSaClSa, and SaClSaCl) through a 2D equivalent linear analysis approach. The seismic response of each profile was simulated using the MIDAS GTS NX commercial software. For each site profile, a strong rock motion record was selected and applied at the base of site profiles to calculate ground responses. Dynamic properties of the soils were chosen from the MIDAS library. The results of the analyses were presented in terms of shear modulus decreasing ratio, surface acceleration, and displacement time histories. Furthermore, the complex interaction between seismic waves and local soil deposits was evaluated through seismic signal coherence analysis. Based on the analyses results, the following conclusion can be drawn:
The analysis peak ground acceleration, peak ground displacement, and modulus decreasing ratio computed from Cl profile clearly demonstrated that the Cl profile significantly influenced computed ground motion. This influence can be attributed to the inherent deformability of clays under cyclic seismic loading;
Compared to other profiles, the Sa profile exhibited lower seismic activity. The peak ground acceleration, peak ground displacement, and modulus decreasing ratio computed from the Sa profile indicate minimal impact on seismic amplifications;
The coherence analysis reveals a robust and clear correlation between accelerations computed from SaClSaCl and ClSaClSa profiles. While there are indications of limited coherence between accelerations computed from Cl and Sa, as well as between Cl and ClSaClSa, and Cl and SaClSaCl profiles, the overall coherence remains strong;
Overall, the significance of this study lies in its ability to shed light on the complex interactions between seismic waves and variability of local soil deposits. This contribution was part of the ongoing research on modeling soil variability and quantification of uncertainty in soils.
Acknowledgment
This study was supported by the Stipendium Hungaricum Scholarship.
References
- [1]↑
E. M. Rathje, A. R. Kottke, and W. L. Trent, “Influence of input motion and site property variabilities on seismic site response analysis,” J. Geotech. Geoenvironmental Eng., vol. 136, no. 4, pp. 607–619, 2010.
- [3]↑
I. M. Idriss, “Evolution of the state of the practice,” in Slides at International Workshop on the Uncertainties in Nonlinear Soil Properties and their Impact on Modeling Dynamic Soil Response, Richmond, California, USA, March 18, 2004, pp. 1–15.
- [4]↑
Y. Guzel, M. Rouainia, and G. Elia, “Effect of soil variability on nonlinear site response predictions: Application to the Lotung site,” Comput. Geotech., vol. 121, 2020, Art no. 103444.
- [5]
D. Stanko, Z. Gülerce, S. Markušić, and R. Šalić, “Evaluation of the site amplification factors estimated by equivalent linear site response analysis using time series and random vibration theory based approaches,” Soil Dyn. Earthq. Eng., vol. 117, pp. 16–29, 2019.
- [6]
O. Kegyes-Brassai, Á. Wolf, Z. Szilvágyi, and R. P. Ray, “Effects of local ground conditions on site response analysis results in Hungary,” in 19th Int. Conf. Soil Mech. Geotech. Eng., Seoul, Korea, September 17–22, 2017, pp. 2003–2006.
- [7]
Q. Sun, X. Guo, and D. Dias, “Evaluation of the seismic site response in randomized velocity profiles using a statistical model with Monte Carlo simulations,” Comput. Geotech., vol. 120, 2020, Art no. 103442.
- [8]
T. Katona, “Options for the treatment of uncertainty in seismic probabilistic safety assessment of nuclear power plants,” Pollack Period., vol. 5, no. 1, pp. 121–136, 2010.
- [9]
T. Katona, “Safety assessment of the liquefaction for nuclear power plants,” Pollack Period., vol. 10, no. 1, pp. 39–52, 2015.
- [10]↑
T. T. Tran, S. R. Han, and D. Kim, “Effect of probabilistic variation in soil properties and profile of site response,” Soils Found, vol. 58, no. 6, pp. 1339–1349, 2018.
- [11]
G. R. Toro, “Probabilistic models of site velocity profiles for generic and site-specific ground-motion amplification studies,” Techical report, no. 779574, 1995.
- [12]
T. Bong, Y. Son, S. Noh, and J. Park, “Probabilistic analysis of consolidation that considers spatial variability using the stochastic response surface method,” Soils Found, vol. 54, no. 5, pp. 917–926, 2014.
- [13]
N. Roy, A. Shiuly, R. B. Sahu, and R. S. Jakka, “Effect of uncertainty in VS- N crrelations on seismic site response analysis,” J. Earth Syst. Sci., vol. 127, pp. 103–123, 2018.
- [14]↑
P. Bazzurro and C. A. Cornell, “Ground-motion amplification in nonlinear soil sites with uncertain properties,” Bull. Seismol. Soc. Am., vol. 94, no. 6, pp. 2090–2109, 2004.
- [15]↑
F. Lopez-Caballero, C. Gelis, J. Regnier, and L. F. Bonilla, “Site response analysis including earthquake input ground motion and soil dynamic properties variability,” in 15th World Conference on Earthquake Engineering, Lisbon, Portugal, September 24–28, 2012, Art no. 23796.
- [16]↑
A. R. Kottke and E. M. Rathje, “Technical manual for Strata,” Pacific Earthquake Engineering Research Report, no. 2008/10, University of California, Berkeley.
- [17]↑
Y. M. A. Hashash, M. I. Musgrove, J. A. Harmon, D. R. Groholski, C. A. Phillips, and D. Park, “Nonlinear and equivalent linear seismic site response of one-dimensional soil columns,” DEEPSOIL 6.1, User Manual, 2016.
- [18]↑
I. M. Idriss and J. I. Sun, “User’s manual for SHAKE91: a computer program for conducting equivalent linear seismic response analyses of horizontally layered soil deposits,” Center for Geotechnical Modeling, Dept. of Civil and Environmental Engineering, University of California, 1992.
- [19]↑
G. Zalachoris and E. M. Rathje, “Evaluation of one-dimensional site response techniques using borehole arrays,” J. Geotech. Geoenviron. Eng., vol. 141, no. 12, pp. 1–15, 2015.
- [20]↑
The R project for statistical computing, Vienna, Austria, 2024. [Online]. Available: https://www.r-project.org/. Accessed: Sep. 15, 2022.
- [21]↑
D. Dimitriou, D. Kenourgios, and T. Simos, “Are there any other safe haven assets? Evidence for ‘exotic’ and alternative assets,” Int. Rev. Econ. Financ., vol. 69, pp. 614–628, 2020.
- [22]↑
MIDAS GTS NX, 2D&3D Geotechnical Finite Element Analysis. [Online], Available: https://www.midasgeotech.com/solution/gtsnx. Accessed: Dec. 08, 2023.
- [23]↑
H. B. Seed, R. T. Wong, I. M. Idriss, and K. Tokimatsu, “Moduli and damping factors for dynamic analyses of cohesionless soils,” J. Geotech. Eng., vol. 112, no. 11, pp. 1016–1032, 1986.
- [24]↑
M. Vucetic and R. Dobry, “Effect of soil plasticity on cyclic response,” J. Geotech. Eng., vol. 117, no. 1, pp. 89–107, 1991.
- [25]↑
J. J. Bommer and A. B. Acevedo, “The use of real earthquake accelerograms as input to dynamic analysis,” J. Earthq. Eng., vol. 8, no. sup001, pp. 43–91, 2004.
- [26]↑
T. D. Ancheta, R. B. Darragh, J. P. Stewart, E. Seyhan, W. J. Silva, B. S. J. Chiou, K. E. Wooddell, R. W. Graves, A. R. Kottke, D. M. Boore, T. Kishida, and J. L. Donahue, “NGA-West2 database,” Earthq. Spectra, vol. 30, no. 3, pp. 989–1005, 2014.
- [27]↑
N. N. Ambraseys, P. Smit, J. Douglas, B. N. Margaris, R. Sigbjörnsson, S. Olafsson, P. Suhadolc, and G. Costa, “Internet site for European strong-motion data,” Boll. Di Geofis. Teor. Ed. Appl., vol. 45, pp. 113–129, 2004.
- [28]↑
N. Roy and R. B. Sahu, “Site specific ground motion simulation and seismic response analysis for Micro zonation of Kolkata,” Geomech. Eng., vol. 4, no. 1, pp. 1–18, 2012.
- [29]↑
M. Barla, A. Di Donna, and D. Sterpi, Eds., Challenges and Innovations in Geomechanics, Proceedings of the 16th International Conference of IACMAG, Turin, Italy, August 30–September 2, 2022, Lecture Notes in Civil Engineering, vol. 288, 2021.
- [30]↑
N. Roy, S. Mukherjee, and R. B. Sahu, “Influence of trapped soft/stiff soil layer in seismic site response analysis,” J. Earth Syst. Sci., vol. 129, 2020, Art no. 171.