Dependence of DuctilityResponse Spectra on the Seismogenic Depth from Finite Element Earthquake Rupture Simulations

Document Type : Seismology and Engineering Seismology


Indian Institute of Technology, Hyderabad


Ongoing research in the development of design philosophies for earthquake resistant structures over the past few decades was initially based on the strength and elastic analysis. Later, design philosophies recognized the deformation to be an important parameter to be considered in design following nonlinear analysis. The maximum design lateral force, for a particular earthquake, acting on structures having multiple natural time periods can be obtained from inelastic response spectrum. Scenario earthquakes characterize the spatio-temporal evolution of fault rupture, which when solved together with the elastodynamic equations can give the acceleration time-history at any point on the surface. For any tectonic regime, based on the past seismicity, a seismogenic depth could be defined based on the depth below which no occurrence of earthquakes was observed in the past. Fixing a certain magnitude, we prescribe the slip on a vertical fault based on statistical relations that exist in literature, and simulate ground motion. The ruptured region is varied, is initially assumed to be closer to the free surface, and is later lowered deeper in intervals of 10 km to emulate larger seismogenic depths. Using the simulated ground motion, we compute the fundamental entity of earthquake engineering: the response spectrum for five depths of hypocenter. Earthquake resistant design of structures is mostly done to allow for large inelastic deformations, giving ductile detailing. Choosing ductility ratios of 2, 4, 6 and 8, this paper describes dependency of elastic and inelastic horizontal spectral acceleration on seismogenic depth, by considering kinematic rupture description of a Mw 6.5 earthquake on a vertical strike slip fault.


  1. Mavroeidis, G.P., Dong, G., and Papageorgiou, A.S. (2004) Near-fault ground motions, and the response of elastic and inelastic single-degreeof-freedom (SDOF) systems. Ear thquake Engineering and Structural Dynamics, 33(9), 1023-1049.
  2. Campbell, K.W., and Bozorgnia, Y. (2008) NGA ground motion model for the geometric mean horizontal component of PGA, PGV, PGD and 5% damped linear elastic response spectra for periods ranging from0.01 to 10 s. Earthquake Spectra, 24(1), 139-171.
  3. Aagaard, B., Kientz, S., Knepley, M., Strand, L., and Williams, C. (2013) PyLith User Manual, Version 2.1.0. , Davis, CA: Computational Infrastructure of Geodynamics.
  4. Balay, S., Gropp, W.D.,McInnes, L.C., and Smith, B.F. (1997) Efficient management of parallelism in object-oriented numerical software libraries. Modern Software Tools for Scientific Computing, 163-202, Birkhauser Boston.
  5. Blacker, T.D., Bohnhoff,W.J., and Edwards, T.L. (1994) CUBIT Mesh Generation Environment. Volume 1: Users'Manual (No. SAND--94-1100). Sandia National Labs.,Albuquerque, NM (United States).
  6. Liu, P. and Archuleta, R.J. (2004) A new nonlinear finite fault inversion with three-dimensional Green's functions: Application to the 1989 Loma Prieta, California, earthquake. Journal of Geophysical Research: Solid Earth, 109(B2).
  7. Somerville, P.G. (2003) Magnitude scaling of the near fault rupture directivity pulse. Physics of the Earth and Planetary Interiors, 137(1), 201-212.
  8. Mai, P.M. and Beroza, G.C. (2000) Source scaling properties from finite-fault-rupture models. Bulletin of the Seismological Society of America, 90(3), 604-615.
  9. Mai, P.M. and Beroza, G.C. (2002) A spatial random field model to characterize complexity in earthquake slip. Journal of Geophysical Research: Solid Earth, 107(B11).
  10. Andrews, D.J. (1976) Rupture velocity of plane strain shear cracks. Journal of Geophysical Research, 81(32), 5679-5687.
  11. Liu, P.,Archuleta, R.J., and Hartzell, S.H. (2006) Prediction of broadband ground-motion time histories: Hybrid low/high-frequency method with correlated random source parameters. Bulletin of the Seismological Society of America, 96(6), 2118-2130.
  12. Iman, R.L. (2008) Latin Hypercube Sampling. JohnWiley & Sons, Ltd.
  13. Acton, C.E., Priestley, K., Mitra, S., and Gaur, V.K. (2011) Crustal structure of the darjeeling-SikkimHimalaya and Southern Tibet. Geophysical Journal International, 184(2), 829-852.
  14. Paul, H., Mitra, S., Bhattacharya, S.N., and Suresh, G. (2015) Active transverse faulting within under thrust Indian crust beneath the Sikkim Himalaya. Geophysical Journal International, 201(2), 1070-1081.