Probable contribution of flingstep effect on the response spectra at near source site
A. Nicknam^{1} , M. A. Barkhordari^{2} , H. Hamidi Jamnani^{3} , A. Hosseini^{4}
^{1, 2, 3, 4}School of Civil Engineering, Iran University of Science & Technology, P. O. Box 16765163, Tehran, Iran
^{3}Corresponding author
Journal of Vibroengineering, Vol. 16, Issue 1, 2014, p. 334340.
Received 15 September 2013; received in revised form 18 October 2013; accepted 25 October 2013; published 15 February 2014
JVE Conferences
Flingstep effect (static offset) besides forward directivity is one the major consequence of nearfault earthquakes. This effect may impose unanticipated seismic demands on structures. Comparing to forward directivity effect, few researches are concerned with the effect of flingstep. In this study, the contribution of flingstep effect on the elastic response spectra was examined by means of a probabilistic method as a new approach in investigation of flingstep effect. To this end, seismological parameters of a significant nearfault earthquake were obtained using inversion solution and MOGA technique and then by randomly changing hypocenter of causative fault many ground motions were synthesized. Eventually, the influence and also the probable contribution of flingstep on the response spectra were achieved via probabilistic process.
Keywords: probable contribution, flingstep, nearfault earthquake, multitaper spectrum.
1. Introduction
Distinct characteristics of nearfield (nearfault) earthquakes with respect to those of farfield signify the special attention to such earthquakes. The accelerograms recorded near to ruptured fault contain two unique characteristics, (i). Rupture directivity (large velocity pulse) and (ii). Large coseismic displacement (also known as permanent static displacement or flingstep). LandersUSA (1992), KocaeliTurkey (1999), ChichiTaiwan (1999) are the most famous nearfiled earthquakes having the two mentioned characteristics. Earthquakes coupled with directivity effect impose extreme demands on structures and are recognized to have caused drastic damage in previous nearfault earthquakes [1]. Such nearfault consequents will cause the greater part of the seismic energy from the rupture to reach a single coherent longperiod pulse of motion. Valuable investigations have been performed by some researchers [16]. The majority of investigations are generally focused on rupture directivity and fewer studies are concentrated on flingstep effect. Therefore, the need to further investigations on flingstep effect seems necessary [4].
The effects of rupture directivity are generally considered by modifications to the elastic response spectra (5 % damping ratio) in current design codes [3, 79]. In this study the focus is on investigating the effect of flingstep on the response spectra using probabilistic approach. To this end, seismological parameters of Kocaeli Earthquake (a distinct sample of nearfault earthquake) were achieved using inversion solution, evolutionary approach (MOGA) and Theoreticalbased Green’s Function (TGF) method with regard to desired objectives (Permanent displacement, response spectra and multitaper spectra [7, 10]). This methodology was previously introduced by the authors for compatible seismogram simulation [7]. Afterwards, by changing the hypocenter location, many new waveforms were generated and the distinguishing ones having considerable flingstep were chosen for probabilistic operations. By extracting flingstep signal from the primary waveforms using a special technique, the effect of flingstep on the response spectra was investigated.
2. Nearfield charactersirics and simulation procedure
Ground motions near to ruptured fault are meaningfully different from those observed further away from the active fault. These motions come with longperiod pulses which reveal two major effects; flingstep and rupture directivity (forward directivity). These effects may impose sever structural damages on structures [8, 11, 12]. The specific characteristics of such longperiod ground motions are the interest of many researchers. The limited numbers of nearfault earthquake demonstrate the shortcoming of such earthquakes in the process of concerned studies. Therefore, lack of recorded ground motions at stations near the ruptured fault reveals the need to overcome this problem by simulating synthetic ground motions. Modeling and simulating ground motions are also the interest of other researchers and those can be useful for seismic analyses [13, 14]. The simulation procedure should be as much as similar to the nature of faulting; hence finding the seismological parameters of an existing nearfault earthquake is of high importance. Among several approaches in the literature for retrieving nearfault ground motion e. g. omegasquared $\left({\omega}^{2}\right)$ methods, semi empiricalbased methods, theoretical green’s functionbased methods and hybrid methods [3, 1519], the theoretical green’s function (TGF) method due to its capability to model both flingstep and rupture directivity was adopted to simulate the desired ground motions [20]. Using inversion solution, MultiObjective Genetic Algorithm (MOGA) and TGF method at 01.5 Hz frequency, seismological parameters of Kocaeli (1999) Earthquake were achieved via a computer code written in MATLAB [21]. For more details about GA technique and inverse problem, please see the research done by Jimenez [22]; and Royzman and Goroshko [23]. The upperlower bounds of parameters were adopted from the previous studies [24, 25]. The proposed Fitness Function (FF) used in this code is defined as the inverse of Error Function (EF) expressed as Eq. 1:
where $X$ refers to the objective parameters ($j=1$: response spectra, $j=2$: multitaper spectra, $j=3$: permanent displacement), $a1$ and $a2$ are recorded and synthesized seismograms and $n$ denotes vector dimension. Table 1 shows the achieved seismological parameters of Kocaeli EQ.
Table 1. Achieved seismological parameters of Kocaeli EQ using MOGA
Longitude (deg)

Latitude (deg)

Depth (km)

Focal mechanism
(strike, dip and rake)

Rupture length (km)

Rupture width (km)

Rupture velocity (m/s)

Slip (m)

29.97

40.7

16.5

91.04, 89.88, 170

100

20

2500

var

Having the seismological parameters in hand and changing the location of hypocenter as a random parameter, many waveforms were generated at YPT station (see Figure 1).
Fig. 1. Changing the location of hypocenter at causative fault for generating new waveforms
3. Flingstep recognition and extraction
Flingstep which is related to permanent tectonic displacement near the causative faults, leads to onesided velocity pulse and consequently to twosided acceleration pulse (integration of velocity pulse). The concerned equation of twosided acceleration pulse is written in the form of Eq. 2 [1]. By forming the mentioned diagrams and equation, the flingsignal can be extracted from the initial acceleration (waveform). Figure 2 displays the flingstep signal and corresponding velocity and displacement diagrams.
Fig. 2. The flingstep pulse as an acceleration and corresponding velocity and displacement
where $t$ is variable of time, $D$ is the permanent displacement (maximum amplitude of the displacement), ${T}_{p}$ is the period of sinusoidal pulse and ${T}_{i}$ is the pulse arrival time. To display the aforementioned technique, the displacement, the acceleration and the corresponding response spectra for a waveform before and after the removal of flingstep signal are shown in Figure 3.
Fig. 3. Effect of flingstep signal on response spectra (red: withfling; black: withoutfling)
The figure demonstrates that removing the fling signal (flingstep effect) from the primary acceleration may increase or decrease the value of response spectra ($Sa$) regarding the amount of period ($T$). This implies that with respect to fundamental period of structure, flingstep effect may decrease or increase the imposed seismic demands.
4. Data analysis and results
The abovementioned procedures and techniques were implemented to the selected accelerograms to gain the goal of investigating the effect of flingstep on the response spectra. The methodology is described concisely in the following flowchart (Figure 4).
Fig. 4. Flowchart of applied procedure
The total number of 32 ground motions (waveforms) having significant flingstep were selected among the 567 synthesized ones and the expressed procedure was applied to them. The flingstep signals were recognized visually via the displacement diagrams and they were extracted from the primary waveforms using supplied computer code. The waveforms, the flingstep signals, the displacement diagram and the concerned response spectra are demonstrated in Figures 57.
Fig. 5. Selected ground motions: (a) with flingstep and (b) without flingstep signal
(a)
(b)
Fig. 6. The extracted flingstep signals: (a), and the corresponding displacements (b)
(a)
(b)
Fig. 7. Response spectra: (a) with flingstep and (b) without flingstep signal
(a)
(b)
In the following, a very important class of statistical distributions named “normal distribution” is used to perform the probabilistic process. Normal distribution curve is symmetric and has “bellshaped” density curve having a single peak. Since it approximates many natural phenomena so well, this distribution has developed into a typical reference for many probability problems [26]. The two major quantities of this distribution $\mu $ (mean) and $\sigma $ (standard deviation) were calculated for each set of data (response spectra; with and without flingstep effect) via a computer code written in MATLAB. The closest existing spectra (having minimum error) to the curves corresponding to probability of 50 % ($\mu $) and 84 % ($\mu +\sigma $) were achieved. Figure 8 demonstrates the comparison of the response spectra in two statuses, with and without flingstep effect. The left diagram is related to the probability of 50 % and the right one is related to 84 %. The difference between the response spectra with and without flingstep is visible for Periods ($T$) more than 1.7 sec.
Fig. 8. Response spectra (wf: with fling; wof: without fling)
(a)
(b)
The Fast Fourier Transform (FFT) of the accelerations (concerning the response spectra shown in Figure 8) is carried out and demonstrated in Figure 9. The difference between the FFT Amplitudes in low frequency range (long period) can justify the difference between the response spectra at Figure 8.
Fig. 9. FFT comparison (wf: with fling; wof: without fling)
The displacements and the corresponding accelerations for each case (50 % and 84 %) are shown in Figure 10. They show the flingstep (permanent displacement) occurring at site for expected probability and the corresponding accelerations could be used as the input of dynamic analysis. Also, one can develop his own risk and probability to achieve the desired curves.
Fig. 10. Acceleration (ground motion) and displacement corresponding to achieved response spectra
5. Conclusion
By incorporating the inversion solution technique, evolutionary approach, engineering seismology, probabilistic method and earthquake engineering, the probable contribution of flingstep was studied. The results imply that the existence of flingstep in ground motions can mostly decrease the value of response spectra, while it can also act adversely (please see Figure 8, left one, $T=$2.4 sec). That is why the word “unanticipated” was used in the beginning of this study for describing the effect of flingstep. A new approach was also introduced that provides a technique taking the account of probability in flingstep investigations. Accordingly, one can determine the amount of the desired safety and choose the probability index ($\mu $, $\mu +\sigma $, $\mu +2\sigma $, etc.) and consequently find the corresponding accelerograms for analysis. The results and techniques of this study could be useful for designing and retrofitting of structures at near source sites with respect to the desired risk or safety.
References
 Kalkan E., Kunnath S. K. Effects of fling step and forward directivity on seismic response of buildings. Earthquake Spectra, Vol. 22, Issue 2, 2006, p. 367390. [Search CrossRef]
 Alavi B., Krawinkler H. Consideration of nearfault ground motion effects in seismic design. Proceedings 12th World Conference on Earthquake Engineering, New Zealand, 2000. [Search CrossRef]
 Somerville P. G. Magnitude scaling of the near fault rupture directivity pulse. Physics of the earth and planetary interiors, Vol. 137, 2003, p. 201212. [Search CrossRef]
 Lestuzzi P., Schwab P., Koller M., Lacave C. How to choose earthquake recordings for nonlinear seismic analysis of structures. 13WCEE, Vol. 1241, 2004. [Search CrossRef]
 Gillie J. J., RodriguezMarek A., McDaniel C. Strength reduction factors for nearfault forwarddirectivity ground motions. Engineering Structures, Vol. 32, Issue 1, 2010, p. 273285. [Search CrossRef]
 Nicknam A., Hosseini A., Hamidi Jamnani H., Barkhordari M. A. Investigation of flingstep effect on the selected structures subjected to the simulated Fandoqa nearsource strong motion. Journal of Vibroengineering, Vol. 16, Issue 1, 2013. [Search CrossRef]
 Nicknam A., Barkhordari M. A., Hamidi Jamnani H., Hosseini A. Compatible seismogram simulation at near source site using multitaper spectral analysis approach. Journal of Vibroengineering, Vol. 15, Issue 2, 2013, p. 626638. [Search CrossRef]
 Somerville P. G., Smith N. F., Graves R. W., Abrahamson N. A. Modification of empirical strong ground motion attenuation relations to include the amplitude and duration effects of rupture directivity. Seism. Res. Lett. 68, 1997, p. 199222. [Search CrossRef]
 ASCE710. Minimum Design Loads for Buildings and Other Structures. American Society of Civil Engineers, 2010. [Search CrossRef]
 Hamidi Jamnani H., Hosseini A., Issa M. A., Nicknam A. Seismic structural performance using theoreticalbased synthetic waveform with respect to the site soil conditions. 2nd International Conference on PerformanceBased Design in Earthquake Geotechnical Engineering, Italy, 2012. [Search CrossRef]
 Lili X., Longjun X., RodriguezMarek A. Representation of nearfault pulsetype ground motions. Earthquake Engineering and Engineering Vibration, Vol. 4, Issue 2, 2005, p. 191199. [Search CrossRef]
 Bray J. D., RodriguezMarek A. Characterization of forwarddirectivity ground motions in the nearfault region. Soil Dynamics & Earthquake Eng., Vol. 24, 2004, p. 815828. [Search CrossRef]
 Šliaupa S. Modeling of the ground motion of the maximum probable earthquake and its impact on buildings. Journal of Vibroengineering, Vol. 15, Issue 2, 2013, p. 532543. [Search CrossRef]
 Wei B., Cui R., Dai G. Seismic performance of a rollingdamper isolation system. Journal of Vibroengineering, Vol. 15, Issue 3, 2013, p. 15041512. [Search CrossRef]
 Nicknam A., Eslamian Y. A hybrid method for simulating nearsource, broadband seismograms. Application to the 2003 Bam earthquake. Tectonophysics, Vol. 487, 2010, p. 4658. [Search CrossRef]
 Bouchon M. A simple method to calculate green’s functions for elastic layered media. Bulletin of the Seismology Society of America, Vol. 71, Issue 4, 1981, p. 959971. [Search CrossRef]
 Koketsu K. Miyake H. A seismological overview of longperiod ground motion. J. Seismol., Vol. 12, 2008, p. 133143. [Search CrossRef]
 Pitarka A., K. Irikura, T. Iwata, H. Sekiguchi Threedimensional simulation of the nearfault ground motion for the 1995 Hyogoken Nanbu (Kobe), Japan earthquake. Bull. Seism. Soc. Am., Vol. 88, 1998, p. 428440. [Search CrossRef]
 Mai P. M., Beroza G. C. A hybrid method for calculating nearsource, broadband seismograms: Application to strong motion prediction. Phys. Earth Planet. Int., Vol. 137, Issue 14, 2003, p. 183199. [Search CrossRef]
 Hisada Y., Bielak J. A Theoretical method for computing nearfault ground motions in layered halfspaces considering static offset due to surface faulting, with a physical interpretation of fling step and rupture directivity. Bulletin of the Seismological Society of America, Vol. 90, Issue 2, 2003, p. 387400. [Search CrossRef]
 Hamidi Jamnani H., Nicknam A., Hosseini A., Barkhordari M. A. Simulation of longperiod TGF ground motions having flingstep for structural performance study. KSCE Journal of Civil Engineering, Paper in review, 2012. [Search CrossRef]
 Jimenez A., Garcia M. J., Romacho M. D. Simultaneous inversion of source parameters and attenuation factor using genetic algorithm. Bull. Seism. Soc. Am., Vol. 95, Issue 4, 2005, p. 14011411. [Search CrossRef]
 Royzman V., Goroshko A. Multiple inverse problem. Journal of Vibroengineering, Vol. 14, Issue 3, 2012, p. 14171424. [Search CrossRef]
 Yagi Y., Kikuchi M. Source rupture process of the Kocaeli, Turkey, earthquake of August 17, 1999, obtained by joint inversion of nearfield data and teleseismic data. Geophysical Research Letters, Vol. 27, Issue 13, 2000, p. 19691972. [Search CrossRef]
 Sekiguchi H., Iwata T. Rupture process of the 1999 Kocaeli, Turkey, earthquake estimated from strongmotion waveforms. Bulletin of the Seismological Society of America, Vol. 92, Issue 1, 2002, p. 300311. [Search CrossRef]
 Bryc Wlodzimierz. The normal distribution: characterizations with applications. SpringerVerlag, ISBN 0387979905, 1995. [Search CrossRef]