Narrowbandness of seismic ground displacement on a broader area of the lithosphere and importance on base motion in isolated structures
César A. Morales^{1}
^{1}Facultad de Ingeniería, Universidad Peruana de Ciencias Aplicadas, Lima, Perú
Journal of Vibroengineering, Vol. 23, Issue 2, 2021, p. 400406.
https://doi.org/10.21595/jve.2020.21682
Received 31 August 2020; received in revised form 8 October 2020; accepted 22 October 2020; published 13 January 2021
JVE Conferences
Power spectral density of horizontal ground displacement of large and recent events in earthquakeprone Latin America is analyzed. The results confirm in a larger region of the world that strong motion horizontal displacement is a narrowband process, which was previously demonstrated solely on a very limited area: the State of California. Nonetheless, those limited and previous results proved to be important in seismic base isolation; particularly, in the solution of the problem of large displacements at the structure base. This is a current problem for which more expensive techniques than passive control are presently being implemented, as active or hybrid control; therefore, it is emphasized that a solution exists within simple base isolation, and it is based on ground displacement narrowbandness.
 Narrowbandness of seismic ground displacement is confirmed
 A study of truly major earthquakes in a full continent is now completed
 A possible important application in Seismic Base Isolation is presented
Keywords: ground displacement, spectral characterization, Latin American earthquakes, base isolation, frequency content, base displacement.
1. Introduction
Records of ground seismic displacement have not attracted much awareness because in dynamic analysis of conventional fixed structures, or the vast majority of buildings, ground acceleration records are the important signals [13]. However, in dynamics of baseisolated structures, ground displacement can be a more important seismic record than the acceleration one [47].
Base isolation is a very successful antiseismic technique [1, 4, 8] that has proved great performance indeed during real earthquakes [9, 10]; nevertheless, it has a problem, that one of large base displacements which is a side effect of the technique. The large base displacement problem (LBD problem) can be explained and tackled considering the relative displacement of the base [11, 12] or studying the absolute base displacement [4, 13]; when considering relative displacement the problem is normally or has been more lately called isolator displacement demand [14, 15].
The LBD problem has been solved, or it has been shown that the absolute base displacement can be reduced to well under ground displacement levels [4]. It has been suggested that the reduction in absolute base displacement implies a reduction in relative base displacement because a) the latter is simply equal to the first minus the ground displacement record (${u}_{b}\left(t\right)s\left(t\right)$) which is set or fixed (at least in a strong motion by strong motion basis) and b) the beforereduction absolute displacement levels are larger than the ground displacement ones $\text{(}{u}_{b}>s\text{)}$ [4]. Nevertheless, this idea or suggestion has not been demonstrated yet; thus, actual relative or isolator displacement mitigation will be discussed in this work as an extension of the application of the results.
Now, this solution of the LBD problem was based on an important hypothesis regarding the modeling of the seismic input. This hypothesis in turn was based on a frequencydomain characterization of actual strong motion ground displacement, which was established through the analysis of the power spectral density (PSD) of displacement records of four (4) large California earthquakes (moment magnitude above 7.0) [4]. Subsequent research or broader results were necessary in the context of Northern America, and next ten (10) California strongmotion events were analyzed [16]; these later spectral results corroborated the previous and limitedinnumber ones; nevertheless, there is still the geographic limitation, or the fact that the earthquakes analyzed were confined or associated solely to the Californian section of the Pacific and North American plates boundary. Thus, results on a broader region of the lithosphere are still necessary.
In this work, broader results are presented, which are from strong motion events associated with the Nazca, South American, Cocos, Caribbean, North American and Pacific plates; that is, the spectral content of ground displacement of 8 very large Latin American earthquakes is studied (2 events with magnitude above 8); more exactly, the PSD of the critical horizontalcomponent record of the closest station to epicenter is analyzed. The results would let us state –on a much wider geographical base– that horizontal ground displacement is a narrowband process. One important point and also an enhancement relative to previous research, is that seismic stations in this larger region of the world are newer or have broadband sensors; therefore, concerns on the low and high frequency noise levels associated with seismic records are reduced, and this point is discussed in the last Section.
The narrowband main result shall be important by itself, but apart from scientific importance, the results will also corroborate that the engineering problem of large base displacements in isolated structures can be solved [4] within simple passive control. More expensive and complex solutions to the LBD problem or the displacement demand in isolated buildings are presently being proposed: active control [6, 17], semiactive control [18], hybrid control with mangnetorheological dampers [19], tuned mass dampers [12, 20] and the inerter or gyromass [21]. In addition, it is pointed out that the study of seismic ground displacement is important also in engineering of submerged tunnels, pipelines and longspan structures [22, 23]. Finally, it is also noted that the suggested solution of the LBD problem is based on a linear model of the isolation system; nonetheless, there are nonlinear models [24] which require advanced numerical methods for the structural response [25].
2. Method and analysis of seismic base isolation application
The procedure to obtain the main results of this investigation are presented first, which are broader ground displacement power spectral density results; secondly, the extension of their application to mitigating the relative base displacement in isolated structures is explained; this motion is actually the isolator deformation or displacement, which is important in the engineering and construction of real isolation systems.
2.1. Methodology
Eight very large and very recent Latin American events are shown in Table 1; a quarter of which with a moment magnitude above 8.0, and all occurred in this decade. The fact that the earliest or oldest earthquake is just from 10 years ago would not in general imply much regarding the issue of noise in seismic records, but in Latin America the difference is in that most active sensors are newer; thus, the lowfrequency noise problem is limited, as it will be discussed later on. The table includes the name of the closesttoepicenter station in each case, which were also the chosen stations in the previous work [16]; moreover, from the usual two horizontalcomponent records provided by each instrument, the one with the largest peak displacement is selected in the study.
The power spectral density (PSD) function of the selected displacement signals are obtained by means of a periodogram technique which is a nonparametric method to estimate spectral density. All 8 signals are corrected records and available in this form in public seismic databases, as strongmotioncenter.org by the USGS. Formally, the plots are estimates of periodogram PSD of samples of ground displacement after regarding the underlying stochastic process as widesense stationary; a rectangular window has been used in all cases to alleviate the leakage phenomenon.
Table 1. Large and recent Latin American events and analyzed stations
Event

Date

Magnitude, ${M}_{m}$

Station

Haiti

1/12/2010

7.0

Presa Sabenta

Maule (Chile)

2/27/2010

8.8

Constitución

Sierra Cucapá (México)

4/4/2010

7.2

Bond’s Corner

Coquimbo (Chile)

9/16/2015

8.3

Pedregal

Ecuador

4/16/2016

7.8

Otavalo

Chiloé (Chile)

12/25/2016

7.6

Espejo Luna

Valparaíso (Chile)

4/24/2017

6.9

Torpederas

Puebla (México)

9/19/2017

7.1

Unam

2.2. Extension of the applications to reduce isolator displacement
As previously indicated, one application of the initial PSD results was in the reduction of the absolute base displacement in isolated structures; in fact, the mitigation was successful or to well under ground displacement levels [4]. Nonetheless, the physical design and construction of an isolation system is more concerned with the relative displacement than with the absolute one; this is mainly because the width of the moat and the flexibility of service connections depend on that isolator displacement (displacement demand) which is established by the dynamic design of the system or by seismic codes.
This relative displacement problem is a current issue in Earthquake Engineering; in fact and as explained previously, more complex and expensive solutions to this LBD problem have recently been proposed [1721]; thus, the spectral content of seismic ground displacement should be revisited; moreover, the narrowband results must be confirmed in larger regions of the world, and as importantly, the aim or application must be on the control of relative base displacement this time around.
With this last objective in mind, for further work and for expanding the application of the PSD results, we propose as solution to the isolator (relative) displacement problem to make use of the frequencydomain concept of relative transmissibility from the field vibration isolation of mechanical engineering. The hypothesis is that if the absolute transmissibility function made possible quite good reductions in absolute base motions [4], good levels of isolator displacement mitigation can be as well attained if the relative transmissibility is considered. This function is not new; however, the concept is not at all as widespread as its absolute counterpart: whereas this can be found in all undergraduate textbooks on Mechanical Vibrations, the relative transmissibility does not appear in most textbooks on the subject; this is because isolator clearance is not a big concern in mechanical engineering, and also because absolute transmissibility is associated not only with displacement isolation but also with force isolation whereas the relative transmissibility applies only for displacement attenuation [26]. As a matter of fact, if we consider the most popular textbooks (English language or worldwide) on mechanical vibrations, the book by Thomson [27] is the only one that has the notion although indirectly rather than as a direct presentation as found in a more advanced or professional references as the Harris’ Shock and Vibration Handbook [26]. In these very few references in which the concept appears, it is established that if the system of Fig. 1 is excited by a stationary harmonic displacement of the support, $s$, the relative displacement of the structure, $x$, is also harmonic with the same frequency $\omega $, as shown in the figure where $X$ and $S$ are the amplitudes of mass and ground motions, respectively; more importantly, the relative transmissibility is the ratio of these harmonic amplitudes, and is further defined by [26, 27]:
where $r$ is the frequency ratio and $\zeta $ damping factor.
As it is proposed herein, the optimization of ${T}_{r}$ can or must be the basis of a design procedure of baseisolation systems when the relative LBD problem is a concern, or its control is an objective, in the full design process of isolated structures. Now, the idea is based on the seismic input being narrow band, which is validated in this paper, and the methodology will be similar to the one published previously to solve the absolute LBD problem [4]. All this is currently the subject of further research work, and we emphasize that the proposal or idea of making use of relative transmissibility is original. Two additional comments or annotations, the proposal will be valid under both, nearfield and farfield conditions, and the displacement demand problem in the isolation technique is not exclusive of civil buildings; structures as electric substation transformers and nuclear power plants do have the same issue under seismic load [28, 29], which can further expand the applications of the spectral results presented.
Fig. 1. Standard onedegreeoffreedom system under harmonic motion of support
3. Ground displacements power spectral density results
The PSD function of the selected strongmotion displacement signals is shown in Fig. 2; it includes the corresponding horizontal component or angular information. These results confirm the narrowband characteristic previously noticed in the works by Morales [4, 16]; note that the plot range 01.5 Hz is already a very slender spectral range. As in the previous works, the narrow band is similar in all cases, except for the Maule event which is justified later; in other words, the energy below 0.025 Hz is close to zero, the peak value is between 0.05 and 0.18 Hz, and above 0.40 Hz the spectral content is negligible; the 0.025 Hz frequency is marked by the vertical dashed line, in all PSDF figures. In addition, the mean value of the peakpower frequency is 0.10 Hz, which is below but close to the 0.13 of the previous two works [4, 16].
A discussion is in order regarding the problem of low frequency, or long period, noise levels associated with strong motion records [30] because this can be severe for frequencies below 0.1 Hz. However, in this regard this work also presents an enhancement relative to previous spectral results: most seismic stations in Latin America are newer or these have broadband sensors; therefore, the lowfrequency noise problem is very limited. This can be proven, the lowcut frequency in the processing of all records is between 0.02 and 0.06 Hz (mean value equal to 0.03 Hz) which is below the frequencies in which we are observing rich spectral content; in other words, it was assumed in the processing or by the analysts, that the signaltonoise ratio was, in general, problematic for frequencies below 0.03 Hz, or periods above 30 s. As an example, most important records of the 1995 Northridge event (California) had a low processed frequency of 0.1 Hz.
The Maule at Constitución PSD is distinct because liquefaction was widespread on the rupture zone; in particular, several types of liquefaction occurred in the small city of Constitución [31]. It is important to note that the same type of spectral content was noticed in a site of the previous work that also suffered liquefaction; that is, peaks in the PSD functions around 0.3 Hz rather than around 0.1 Hz, which is the case in a normal or most of the records analyzed [16].
Fig. 2. Power spectral density of ground displacement
a)
b)
c)
d)
e)
f)
g)
h)
4. Conclusions
Strongmotion records on a much broader area of the lithosphere have been analyzed. It can be concluded or it is confirmed that horizontal ground displacement is practically a narrowband process; this is based on the analysis of 8 very large earthquakes (2 events with magnitude above 8) in the vast Latin American region. These results corroborate previous ones, and that the absolute LBD problem in isolated structures can be solved by vibrations analysis in the frequency domain.
Moreover, a proposal has been presented to extend the application of the further and broader results to control the isolator displacement or deformation, which from an engineering or construction point of view, is more important than the control of the absolute base displacement; this is because the relative displacement is the one having direct impact on the seismic gap or the unusable and costly space between structure and moat walls, and on the expense of flexible utility connections and size of the isolators. This is a presentday problem for which more expensive techniques than simple base isolation are being implemented, as active and hybrid control, and also complex inertial additions as tuned mass dampers and inerters.
The lowest frequency of the band for this set can be established as 0.025 Hz; it had previously been set as 0.049 Hz for strongmotion records in a limited region or the State of California, USA; this difference or presence of lower dominant frequencies in some events is the result of geological diversity, which in a whole continent is richer than in a single state of a country.
Acknowledgements
Funding by Universidad Peruana de Ciencias Aplicadas, Grant number C0252020.
References
 Chopra A. K. Dynamics of Structures. Pearson, USA, 2012. [Search CrossRef]
 Berg G. V., Housner G. W. Integrated velocity and displacement of strong earthquake ground motion. Bulletin of the Seismological Society of America, Vol. 51, 1961, p. 175189. [Search CrossRef]
 Lai S. P. Statistical characterization of strong ground motions using power spectral density function. Bulletin of the Seismological Society of America, Vol. 72, 1982, p. 259274. [Search CrossRef]
 Morales C. A. Transmissibility concept to control base motion in isolated structures. Engineering Structures, Vol. 25, 2003, p. 13251331. [Publisher]
 Kelly J. M., Leitmann G., Soldatos A. G. Robust control of base isolated structures under earthquake excitation. Journal of Optimization Theory and Applications, Vol. 53, 1987, p. 159180. [Publisher]
 Meirovitch L., Stemple T. J. Nonlinear control of structures in earthquakes. Journal of Engineering Mechanics, Vol. 123, 1997, p. 10901095. [Publisher]
 Tadjbakhsh I. G., Rofooei F. Optimal hybrid control of structures under earthquake excitation. Earthquake Engineering and Structural Dynamics, Vol. 21, 1992, p. 233252. [Publisher]
 Naeim F., Kelly J. M. Design of Seismic Isolated Structures. John Wiley and Sons, USA, 1999. [Publisher]
 Nagarajaiah S., Xiaohong S. Response of baseisolated USC Hospital building in Northridge Earthquake. Journal of Structural Engineering, Vol. 126, 2000, p. 11771186. [Publisher]
 Komodromos P. Seismic Isolation for EarthquakeResistant Structures. Wit Press, UK, 2000. [Search CrossRef]
 Kelly J. M. The role of damping in seismic isolation. Earthquake Engineering and Structural Dynamics, Vol. 28, 1999, p. 320. [Publisher]
 Taniguchi T., Der Kiureghian A., Melkumyan M. Effect of tuned mass damper on displacement demand of baseisolated structures. Engineering Structures, Vol. 30, 2008, p. 34783488. [Publisher]
 Barbat A. H., Rodellar J., Ryan E. P. Active control of nonlinear baseisolated buildings. Journal of Engineering Mechanics, Vol. 121, 1995, p. 676684. [Publisher]
 De Domenico D., Ricciardi G. An enhanced base isolation system equipped with optimal tuned mass damper inerter (TMDI). Earthquake Engineering and Structural Dynamics, Vol. 47, 2018, p. 11691192. [Publisher]
 Zargar H., Ryan K. L., Rawlinson T. A., Marshall J. D. Evaluation of a passive gap damper to control displacements in a shaking test of a seismically isolated threestory frame. Earthquake Engineering and Structural Dynamics, Vol. 46, 2017, p. 5171. [Publisher]
 Morales C. A. Further seismic displacement PSDF results. Structural Engineering and Mechanics, Vol. 34, 2010, p. 663666. [Publisher]
 Chang C., Shia S., Yang C. Use of active control algorithm for optimal design of baseisolated buildings against earthquakes. Structural and Multidisciplinary Optimization, Vol. 58, 2018, p. 613626. [Publisher]
 Vu D. C., Politopoulos I., Diop S. A new semi‐active control based on nonlinear inhomogeneous optimal control for mixed base isolation, Structural control and health monitoring, Vol. 25, 2018, p. e2032. [Publisher]
 Madalonni G., Caterino N., Ochiuzzi A. Shake table investigation of a structure isolated by recycled rubber devices and magnetorheological dampers, Structural Control and Health Monitoring, Vol. 24, 2017, p. e1906. [Publisher]
 Xiang P., Nishitani A. Optimum design for more effective tuned mass damper system and its application to baseisolated buildings. Structural Control and Health Monitoring, Vol. 21, 2014, p. 98114. [Publisher]
 Saitoh M. On the performance of gyromass device for displacement mitigation in base isolation systems. Structural Control and Health Monitoring, Vol. 19, 2012, p. 246259. [Publisher]
 Iai S., Kurata E., Tsuchida H., Hayashi S. Integration of strong motion accelerograms. Proceedings of the 10th Joint UJNR Panel Conference, Gaithersburg, USA, 1978. [Search CrossRef]
 Powell G. H. Seismic response analysis of above‐ground pipelines. Earthquake Engineering and Structural Dynamics, Vol. 6, 1978, p. 157165. [Publisher]
 Vaiana N., Sessa S., Marmo F., Rosati L. An accurate and computationally efficient uniaxial phenomenological model for steel and fiber reinforced elastomeric bearings. Composite Structures, Vol. 211, 2019, p. 196212. [Publisher]
 Greco F., Luciano R., Serino G., Vaiana N. A mixed explicitimplicit time integration approach for nonlinear analysis of baseisolated structures. Annals of Solid and Structural Mechanics, Vol. 10, 2018, p. 1729. [Publisher]
 Harris C. M., Piersol A. G. Harris’ Shock and Vibration Handbook. McGrawHill, USA, 2002. [Search CrossRef]
 Thomson W. T., Dahleh M. D. Theory of Vibration with Applications. PrenticeHall, USA, 1998. [Search CrossRef]
 Ala M., Feizi B., Kempner L., Alston D. On Seismic response of substation equipment and application of base isolation to transformers. IEEE Transactions on Power Delivery, Vol. 25, 2010, p. 177186. [Publisher]
 Sarebanha A., Mosqueda G., Kim M. K., Kim J. H. Seismic response of base isolated nuclear power plants considering impact to moat walls. Nuclear Engineering and Design, Vol. 328, 2018, p. 5872. [Publisher]
 Boore D. M., Bommer J. J. Processing of strongmotion accelerograms: needs, options and consequences. Soil Dynamics and Earthquake Engineering, Vol. 25, 2005, p. 93115. [Publisher]
 Verdugo R., Gonzalez J. Liquefactioninduced ground damages during the 2010 Chile earthquake. Soil Dynamics and Earthquake Engineering, Vol. 79, 2015, p. 280295. [Publisher]