Reliable seismic assessment of high-speed railway track-bridge systems using validated physics-based spatial wavefield simulations

Abstract

Empirically based seismic - analysis methods that rely on recorded ground motions are routinely used to evaluate structural damage and guide the design of high-speed railway track-bridge systems (HSRTBS). Their accuracy, however, is constrained by the sparse distribution of seismic stations, especially in rugged terrain where ground-motion variability is pronounced. To overcome this limitation, this work utilizes high-precision, physics-based numerical simulations to model complex spatial seismic wave fields. Leveraging the sub-fault input interface provided by SPECFEM3D, we construct a kinematic hybrid source model of the 2016 Meinong M_w 6.4 earthquake and simulate the complete source-path-site process. The resulting wave field is validated against observed records and by comparing RotD50 intensity measures with the NGA-West 2 ground-motion prediction equations. After verification, a detailed dynamic analysis of the HSRTBS is performed under both the synthetic and recorded seismic excitations. The results show excellent agreement, demonstrating that physics-based ground-motion simulations can effectively characterize the spatiotemporal evolution of seismic wave fields in complex terrain. This method addresses the limitations of sparse seismic station records and provides reliable, site-specific, and spatially variable seismic inputs. These findings confirm the feasibility and robustness of physics-based ground-motion simulations as a reliable alternative to conventional empirical approaches for the seismic analysis of the HSRTBS.

Keywords

high-speed railway / track-bridge systems / seismic analysis / spatial seismic wavefields / spectral element method / kinematic hybrid source model

References

[1] Zhou W-b, Peng D-h, Liu L-l, et al. . Transverse seismic analysis of high-speed railway bridge in China based on a simplified calculation model [J]. Journal of Central South University, 2023, 30(1): 351-364

[2] He X-h, Wu T, Zou Y-f, et al. . Recent developments of high-speed railway bridges in China [J]. Structure and Infrastructure Engineering, 2017, 13(12): 1584-1595

[3] Jiang L-z, Zhou W-b, Wei B, et al. . Research progress of train-track-bridge system safety of high-speed railway under earthquake action [J]. China Civil Engineering Journal, 2020, 53(9): 1-13 (in Chinese)

[4] Song T-r, Pu H, Yang T Y, et al. . Optimizing net present values of risk avoidance for mountain railway alignments with seismic performance evaluation [J]. Computer-Aided Civil and Infrastructure Engineering, 2024, 39(7): 944-962

[5] Graves R W. Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences [J]. Bulletin of the Seismological Society of America, 1996, 86(4): 1091-1106

[6] Satoh T. Three-dimensional finite-difference waveform modeling of strong motions observed in the Sendai Basin, Japan [J]. Bulletin of the Seismological Society of America, 2001, 91(4): 812-825

[7] Pitarka A. 3D Elastic finite-difference modeling of seismic motion using staggered grids with nonuniform spacing [J]. Bulletin of the Seismological Society of America, 1999, 89(1): 54-68

[8] Bielak J, Ghattas O, Kim E J. Parallel octree-based finite element method for large-scale earthquake ground motion simulation [J]. Computer Modeling in Engineering & Sciences, 2005, 10(2): 99-112

[9] Furumura T, Kennett B L N, Takenaka H. Parallel 3-D pseudospectral simulation of seismic wave propagation [J]. Geophysics, 1998, 63(1): 279-288

[10] Furumura T, Koketsu K. Specific distribution of ground motion during the 1995 Kobe Earthquake and its generation mechanism [J]. Geophysical Research Letters, 1998, 25(6): 785-788

[11] Sánchez-Sesma F J, Campillo M. Diffraction of P, SV, and Rayleigh waves by topographic features: A boundary integral formulation [J]. Bulletin of the Seismological Society of America, 1991, 81(6): 2234-2253

[12] Ba Z-n, Liu Y, Liang J-w, et al. . Near-fault broadband seismograms synthesis in a stratified transversely isotropic half-space using a semi-analytical frequency-wavenumber method [J]. Engineering Analysis with Boundary Elements, 2023, 146: 1-16

[13] Komatitsch D, Vilotte J P. The spectral element method: An efficient tool to simulate the seismic response of 2D and 3D geological structures [J]. Bulletin of the Seismological Society of America, 1998, 88(2): 368-392

[14] Lee S J, Chan Y C, Komatitsch D, et al. . Effects of realistic surface topography on seismic ground motion in the Yangminshan Region of Taiwan based upon the spectral-element method and LiDAR DTM [J]. Bulletin of the Seismological Society of America, 2009, 99(2A): 681-693

[15] He C-h, Wang J-t, Zhang C-h, et al. . Simulation of broadband seismic ground motions at dam canyons by using a deterministic numerical approach [J]. Soil Dynamics and Earthquake Engineering, 2015, 76: 136-144

[16] Ba Z-n, Zhao J-x, Zhang Y-s, et al. . The 3D physics-based broadband ground-motion simulation, including topography effects causing damage to structures, for the 2022 mW 6.6 Menyuan, China, earthquake [J]. Bulletin of the Seismological Society of America, 2024, 114(3): 1481-1506

[17] Zhao J-x, Ba Z-n, Kuo C-y, et al. . Broadband ground motion simulations applied to the Luding M_S 6.8 earthquake on September 5, 2022 based on spectral element method [J]. Acta Seismologica Sinica, 2023, 45(2): 179-195 (In Chinese)

[18] Ba Z-n, Fu J-s, Zhu Z-h, et al. . An end-to-end 3d seismic simulation of underground structures due to point dislocation source by using an FK-FEM hybrid approach [J]. Frontiers of Structural and Civil Engineering, 2022, 16(12): 1515-1529

[19] Wu M-t, Wang B, Ba Z-n, et al. . Seismic simulation of seismic source - basin - structural integral physical model: the FK-SE-FE approach [J]. Journal of Vibration Engineering, 2026, 39(1): 1-13 (in Chinese)

[20] Wang J-t, Jin F, Zhang C-han. Nonlinear seismic response analysis of high arch dams to spatially-varying ground motions [J]. International Journal of Civil Engineering, 2019, 17(4): 487-493

[21] Soto V, Sáez E, Magna-Verdugo C. Numerical modeling of 3D site-city effects including partially embedded buildings using spectral element methods. Application to the case of Viña del Mar city, Chile [J]. Engineering Structures, 2020, 223: 111188

[22] Andrews D J. A stochastic fault model: 2. Time-dependent case [J]. Journal of Geophysical Research: Solid Earth, 1981, 86(B11): 10821-10834

[23] Irikura K, Miyake H. Recipe for predicting strong ground motion from crustal earthquake scenarios [J]. Pure and Applied Geophysics, 2011, 168(1): 85-104

[24] Graves R, Pitarka A. Broadband time history simulation using a hybrid approach [C]. 13th World Conference on Earthquake Engineering: WCEE, 2004, 8: 1098

[25] Graves R W, Pitarka A. Broadband ground-motion simulation using a hybrid approach [J]. Bulletin of the Seismological Society of America, 2010, 100(5A): 2095-2123

[26] Graves R, Pitarka A. Refinements to the Graves and pitarka (2010) broadband ground-motion simulation method [J]. Seismological Research Letters, 2015, 86(1): 75-80

[27] Ba Z-n, Zhao J-x, Zhang Y-s, et al. . Broadband ground motion spectral element simulation based on GP14.3 kinematic hybrid source model and SPECFEM 3D [J]. Chinese Journal of Geophysics, 2023, 66(3): 1125-1138

[28] Graves R, Pitarka A. Kinematic ground-motion simulations on rough faults including effects of 3D stochastic velocity perturbations [J]. Bulletin of the Seismological Society of America, 2016, 106(5): 2136-2153

[29] Wang X-c, Wang J-t, Zhang C-han. Deterministic full-scenario analysis for maximum credible earthquake hazards [J]. Nature Communications, 2023, 14: 6600

[30] Jiang W, Tao X-x, Zhao K. Scaling laws of the global parameters of source models from NGA data [J]. China Earthquake Engineering Journal, 2017, 39(2): 221-226 247

[31] Jiang W, Tao X-x, Tao Z-r, et al. . Scaling laws of local parameters of finite fault source model [J]. Earthquake Engineering and Engineering Vibration, 2017, 37(6): 23-30 (in Chinese)

[32] Somerville P, Irikura K, Graves R, et al. . Characterizing crustal earthquake slip models for the prediction of strong ground motion [J]. Seismological Research Letters, 1999, 70(1): 59-80

[33] Mai P M, Beroza G C. A spatial random field model to characterize complexity in earthquake slip [J]. Journal of Geophysical Research: Solid Earth, 2002, 107(B11): 2308

[34] Dreger D S, Beroza G C, Day S M, et al. . Validation of the SCEC broadband platform V_14.3 simulation methods using pseudospectral acceleration data [J]. Seismological Research Letters, 2015, 86(1): 39-47

[35] Ba Z-n, Zhao J-x, Zhu Z-h, et al. . 3D physics-based ground motion simulation and topography effects of the 05 September 2022 MW_6.6 Luding earthquake, China [J]. Soil Dynamics and Earthquake Engineering, 2023, 172: 108048

[36] Kagawa T, Irikura K, Somerville P G. Differences in ground motion and fault rupture process between the surface and buried rupture earthquakes [J]. Earth, Planets and Space, 2004, 56(1): 3-14

[37] Day S M. Memory-efficient simulation of anelastic wave propagation [J]. Bulletin of the Seismological Society of America, 2001, 91: 520-531

[38] Komatitsch D, Tromp J. Spectral-element simulations of global seismic wave propagation: II. Three-dimensional models, oceans, rotation and self-gravitation [J]. Geophysical Journal International, 2002, 150(1): 303-318

[39] Yang Y-h, Chen Q, Diao X, et al. . New interpretation of the rupture process of the 2016 Taiwan Meinong Mw 6.4 earthquake based on the InSAR, 1-Hz GPS and strong motion data [J]. Journal of Geodesy, 2021, 95(11): 121

[40] Rau R J, Wen Y-y, Ching K E, et al. . Origin of coseismic anelastic deformation during the 2016 Mw 6.4 Meinong Earthquake, Taiwan [J]. Tectonophysics, 2022, 836: 229428

[41] Mai H A, Lee J C, Chen K H, et al. . Coulomb stress changes triggering surface pop-up during the 2016 Mw 6.4 Meinong earthquake with implications for earthquake-induced mud diapiring in SW Taiwan [J]. Journal of Asian Earth Sciences, 2021, 218: 104847

[42] Wang L-y, Gao H, Feng G-cai. Insar and GPS earthquake source parameter inversion for the 2016 M_w6.4 Meinong, Taiwan earthquake [J]. Chinese Journal of Geophysics, 2017, 60(4): 346-357

[43] Lu S-d, Chen B-c, Zhang Y-sheng. 20160206 earthquake geological survey report [R], 2016, Taibei, Central Geological Survey R.O.C

[44] Pitarka A, Akinci A, De Gori P, et al. . Deterministic 3D ground-motion simulations (0–5 hz) and surface topography effects of the 30 October 2016 mW 6.5 norcia, Italy, earthquake [J]. Bulletin of the Seismological Society of America, 2022, 112(1): 262-286

[45] Graves R W. Stability and accuracy analysis of coarse-grain viscoelastic simulations [J]. Bulletin of the Seismological Society of America, 2003, 93(1): 283-300

[46] Boore D M. Orientation-independent, nongeometric-mean measures of seismic intensity from two horizontal components of motion [J]. Bulletin of the Seismological Society of America, 2010, 100(4): 1830-1835

[47] Boore D M. Orientation-independent measures of ground motion [J]. Bulletin of the Seismological Society of America, 2006, 96(4A): 1502-1511

[48] Ning C-x, Xie Y-zhou. Convolutional variational autoencoder for ground motion classification and generation toward efficient seismic fragility assessment [J]. Computer-Aided Civil and Infrastructure Engineering, 2024, 39(2): 165-185

[49] Campbell K W, Bozorgnia Y. NGA-West2 ground motion model for the average horizontal components of PGA, PGV, and 5% damped linear acceleration response spectra [J]. Earthquake Spectra, 2014, 30: 1087-1115

[50] Chiou B S, Youngs R R. Update of the chiou and youngs NGA model for the average horizontal component of peak ground motion and response spectra [J]. Earthquake Spectra, 2014, 30: 1117-1153

[51] Idriss I M. An NGA-West2 empirical model for estimating the horizontal spectral values generated by shallow crustal earthquakes [J]. Earthquake Spectra, 2014, 30: 1155-1177

[52] Abrahamson N A, Silva W J, Kamai R. Summary of the ASK14 ground motion relation for active crustal regions [J]. Earthquake Spectra, 2014, 30: 1025-1055

[53] Wei B, Hu Z-l, Zuo C-j, et al. . Effects of horizontal ground motion incident angle on the seismic risk assessment of a high-speed railway continuous bridge [J]. Archives of Civil and Mechanical Engineering, 2021, 21(1): 18

[54] Zhu Z-h, Tang Y-j, Ba Z-n, et al. . Seismic analysis of high-speed railway irregular bridge–track system considering V-shaped canyon effect [J]. Railway Engineering Science, 2022, 30(1): 57-70

[55] GB 50111 - 2006. Code for seismic design of railway engineering [S], 2009, Beijing, China Planning Press (in Chinese)

RIGHTS & PERMISSIONS

Central South University

Just Accepted

This article has successfully passed peer review and final editorial review, and will soon enter typesetting, proofreading and other publishing processes. The currently displayed version is the accepted final manuscript. The officially published version will be updated with format, DOI and citation information upon launch. We recommend that you pay attention to subsequent journal notifications and preferentially cite the officially published version. Thank you for your support and cooperation.