1. Introduction
There has been a constant focus on the characteristic of near-fault earthquake motion since the Northridge event in the U.S., the Kobe event in Japan, and the Chi-Chi event in Taiwan [
1]. The increased availability of recorded ground motions from these seismic events suggests that the dynamic characteristics of ground shaking can significantly vary as a function of the recording station’s location with respect to the fault and evolution of the rupture process [
2]. The long pulse period and large velocity amplitude of the near-fault earthquake motion may cause severe damage to massive engineering structures [
3]. The empirical mode decomposition (EMD) is used to identify and separate the dominant pulse [
4]. Wen et al. [
5] found that the pulse-like near-fault ground motions can significantly increase the displacement demand of structures under a medium period.
The near-fault earthquake motion is roughly defined by the record whose closest distance-to-fault is within 20–30 km. The effect of rupture direction and the permanent ground displacement, are two significant features of the near-fault earthquake motion [
6]. Larger pulse motions exist in the direction normal to the fault than that parallel to it, due to the radiation pattern of shear dislocation on the fault [
7]. The rupture source propagates into the surrounding medium at a speed close to shear velocity, which causes rupture energy arrival in a few pulse patterns. The horizontal component, perpendicular to the fault, contains the long period pulse-like motion [
8]. This type of phenomenon is very similar to the Doppler Effect concept in acoustics. The peak ground acceleration, velocity, and displacement are the critical parameters that influence the structural seismic response [
9]. Near-fault recordings from recent earthquakes indicate that the pulse is a narrow band waveform, whose period increases with magnitude [
10]. The closest distance-to-fault (R) is just a convenient engineering parameter that roughly identifies the pulse-like motion. A more rigorous identification and representation approach should be considered to quantitatively classify the near-fault pulse-like motions.
No universal accepted criteria exist regarding how to identify the most unfavorable seismic input direction with the largest pulse energy, and the knowledge about the characteristic of the seismic signal in that unfavorable direction is limited. Many researchers have investigated the characteristic of the pulse-like earthquake motion, and its influence on the seismic response of structures: Alavi et al., input the equivalent pulse into the elastic and elastic–plastic frame structure to investigate the influence of pulse motion on frame structure. They found that frame structures behave differently depending on whether their base frequency is smaller or larger than pulse frequency [
11]. Arghya et al. investigated the influence of bidirectional near-fault excitations on reinforced concrete (RC) bridge piers. The influence of bidirectional shaking is accounted for by using a simplified 30% rule. Bidirectional interaction under near-fault motion is observed to substantially amplify damage, particularly for a stiff system [
12]. Sehhati et al. studied the effect of near-fault ground motions on the multi-story structures. They found that pulse-like forward-directivity ground motions impose a larger ductility demand on the structure compared to ordinary ground motions [
13]. Based on the value of maximum fractional signal energy contribution by any half-cycle of the velocity time-history, Mukhopadhyay et al. proposed an objective criterion to differentiate directivity pulse-like motions from the available suite of recorded ground motions [
14]. Ghahari et al. utilized the moving average filtering with appropriate cut-off frequency to decompose the near-fault ground motions into two components: Pulse-Type record and BackGround record. They found that the spectra of near-fault ground motions typically have two distinct local peaks that are representatives of the high- and low-frequency components [
15]. Tang et al. proposed an approach to identify the pulse-like motions in earthquake recordings, based on the congruence relationship between the response spectrum and the dimensionless Π-response spectrum [
16]. However, only one channel of the recorded signals was considered in this approach. Baker proposed a wavelet-based approach to quantitatively classify near-fault earthquake records. Pulse-related parameters can be extracted based on the one pulse waveform with the largest pulse energy [
1]. This approach is quite useful at identifying whether the one set near-fault earthquake motion is pulse-like or not. However, the pulse component extracted by this method is composed of one pulse waveform, which is not an adequate representation of the interested pulse feature of earthquake motion records if the multi-pulse waveform is of concern.
For an earthquake event with highly non-uniform slip distribution, such as the Chi-Chi event in Taiwan (1999-September-20 17:47:16 UTC, Mw = 7.6), the type of pulse sequence observed depends on the instrument’s distance relative to the asperities. These factors contribute to the existence of the multi-pulse near-fault earthquake motion [
17]. The number of dominant pulses in the velocity time-history might be related to the number of finite asperities in a fault. However, it is difficult to estimate the slip distribution pattern in a destructive fault a priori [
18].
Many attempts have been made to investigate the characteristics of the near-fault earthquake motion by identifying or extracting the multi-pulse components in the recorded seismic signal. For example, a sum of two and three velocity pulses was utilized to create pulse representations for two records [
19]. However, only one component from strike-normal or strike-parallel was considered. The strike-normal or strike-parallel direction may not be the one containing the largest pulse energy, the unfavorable horizontal direction with the largest pulse energy should be first obtained before analyzing the corresponding characteristic of multi-pulse components [
1]. M&P [
20] and DB4 [
14], are adopted as the mother wavelet to transform near-fault earthquake motions with variant profiles, respectively. This manually selected mother wavelet cannot fit the variant characteristic of near-fault earthquake motions. Furthermore, the most unfavorable horizontal direction is not considered in Reference [
20].
In order to identify and extract the most unfavorable multi-pulse components embedded in the near-fault earthquake motion, we propose a novel adaptive wavelet transform algorithm. Instead of using one fixed mother wavelet for all seismic signals, an iterative algorithm was proposed to obtain the optimum mother wavelet for each potential pulse. Thus, the resulting multi-pulse component could be represented in higher quality.
3. Adaptive Mother Wavelet Choosing
Wavelet analysis is used as a method of transforming time-sequential data into data on a time-frequency plane [
21]. During the last thirty years, this analysis technique has been under rapid theoretical development and has been used to solve many problems [
22]. The fundamental principle of wavelet transform can be found in numerous textbooks. Thus only the near-fault earthquake-related aspects will be discussed in this section.
It is beneficial for the interpretation of the fundamental principle of wavelet transform to be compared with the Fourier transform. In Fourier transform, seismic signals are approximated by summarizing a series of infinite sinusoid signals with solo frequencies. This method uses the time averaging technique to decompose signals, thus phase information cannot be reserved by the Fourier transform. While in the wavelet transform, recorded seismic signals are represented by summarizing a group of wavelet waveforms, each being a narrow-band signal located at the specific time point, which is best suited for the representation of the nonlinear and non-stationary seismic pulse. The dominant pulse waveforms of interest can be represented briefly by just a few wavelets, with the elaborately selected mother wavelets, time, and scale parameters. Both the amplitude and phase information can be reserved by the wavelet transform.
The wavelet function at time t is defined mathematically by Equation (1), where
is the mother wavelet,
s is the scale factor,
l is time location parameter. The parameter
s is intended to dilate the mother wavelet, which scales the central frequency of the mother wavelet to match the interested pulse waveform embedded in the signal. Parameter
l is intended to translate the wavelet along the time axis. There are many kinds of mother wavelets used in scientific research. Different mother wavelets result in different time–frequency planes [
23]. In order to represent the dominant pulse waveform embedded in the strong ground motion with the best resolution, careful selection of the optimum mother wavelet is always necessary. If the profile of the chosen mother wavelet is close to the interested pulse waveform, then a limited number of wavelets is enough to represent the main profile of the signal, with relatively high resolution. Considering the variant feature of the near-fault earthquake motions, different mother wavelets should be adaptively used to extract each individual dominant pulse embedded in the near-fault earthquake motion, instead of using a single mother wavelet for all seismic signals as in References [
1,
20].
Regarding the typical profile of near-fault earthquake motions, the mother wavelets shown in
Figure 2 were used in this paper as the candidate mother wavelets. The repository was composed of 15 types of mother wavelet, including the Haar wavelet, the Gaussian wavelet family from orders 1 to 8, the Daubechies wavelet family from orders 2 to 6, and the Morlet wavelet. These mother wavelets were adequate for the majority of pulse-like seismic signals.
An adaptive procedure when choosing a mother wavelet is necessary to achieve the best resolution for the time–frequency plane by continuous wavelet transform (CWT). For the specific seismic signal , the waveform of extracted dominant pulses should represent the main profile of the source signal to the maximum extent. In order to evaluate the level of similarity for each mother wavelet in the wavelet repository, an iteration process was carried out to extract the dominant pulse embedded in the source signal by CWT, using mother wavelet .
Many quantitative approaches have been proposed in recent years to evaluate the similarity between the signal and candidate mother wavelets, such as the minimum description length criterion, maximum cross-correlation coefficient criterion, the mean squared error of wavelet coefficients, the evaluation criterion, etc. [
23]. In this paper, the level of similarity was evaluated by determining the minimum cross variance (MCV) between the source signal
and the extracted pulse
. The mother wavelet corresponding to the minimum cross variance was regarded as the optimum mother wavelet
. This adaptive procedure was invoked in two situations within the whole analysis procedure: in
Section 4.1, before identifying the most unfavorable horizontal direction, to obtain the optimum mother wavelet for both East-West (EW) and North-South (NS) components; and in
Section 4.2, the procedure was invoked before extracting the individual pulse
from the residual signal
in every iteration until the termination criterion was reached.
4. Identification and Representation of Multi-Pulse Near-Fault Earthquake Motion
At the seismic station, only the two horizontal recorded components were of concern, effects of the vertical component were not considered in the present work. There are two types of wavelet transforms when decomposing a signal into its time–frequency plane: the continuous wavelet transform (CWT), and the discrete wavelet transform (DWT). The main difference is the arrangement of the number of scales and locations used to calculate the wavelet coefficients. CWT continuously calculates the coefficients at every point of the scale-location map, while DWT calculates the coefficients at selected points but with a relatively efficient algorithm. If the profile of the chosen mother wavelet is close to the feature of interest, the corresponding wavelet coefficient will be larger than that around it. Although redundant coefficients are produced by CWT, it is helpful to locate the exact point in the time–frequency plane with the maximum wavelet coefficient, which can be used to represent the dominant pulses embedded in the signal. Usually, the non-pulse part of the signal corresponds to smaller wavelet coefficients. With well-arranged scales and locations, the CWT is capable of providing a high-resolution time–frequency plane. Thus, in this paper, the CWT was adopted instead of DWT to project the seismic signal into its time–frequency plane.
4.1. Identification of the Most Unfavorable Horizontal Direction
Specific directions in the horizontal plane may contain larger pulse energies compared to other horizontal directions [
1]. The most unfavorable horizontal direction (denoted as
hereinafter) may not be the seismic station direction. The structural response was controlled by the structure’s period ratio and the input seismic signal. However, for one specific structure, whose basic period is a specific value, the seismic signal at the direction of
was responsible for the most unfavorable seismic response. Therefore, it is important to identify this direction and investigate its engineering features. In this paper, the horizontal direction
was defined as zero toward the direction of channel 1 at a seismic station, and positive clockwise. Some seismic stations’ main axis were set according to the local fault direction. While, some other seismic station’s main axis were set perpendicular to NS/EW. This paper rotates the seismic station to the NS/EW direction for convenience. Namely, the EW component corresponded to
, and the NS component corressponded to
.
To identify the most unfavorable horizontal direction
, the two-source horizontal components
and
were projected into their time-frequency planes
CEW and
CNS, respectively. In the time–frequency plane, the x-axis was time and the y-axis was the scale value predefined. The scale value in the y-axis was converted into the instantaneous frequency by Equation (2), where
is the central frequency of the chosen mother wavelet,
is the scale value in the y-axis, and
is the sampling period of the signal. A reasonable period band should be set to cover the whole interested range of the period for the multi-pulse near-fault earthquake motion. In the proposed method, this period band was set from 0.25 s to 15 s, which was adequate to cover all periods of interest.
The wavelet coefficient at location () in the time–frequency plane is defined by Equation (3), where is one of the source signals ( or ). The absolute value of the wavelet coefficient in the time–frequency plane is an indicator of the pulse energy level. The two time–frequency planes obtained above ( and ) can be combined into the resultant time–frequency plane by the square sum, since the energy coefficient is a kind of scalar value. The wavelet coefficient corresponding to the maximum pulse energy (denoted by ) was obtained by finding the largest absolute value of the coefficient in the plane , where and are the time and scale locations for , respectively. Because the scope of the predefined scale gave a time–frequency plane with low resolution, it was difficult to locate the exact scale location of . Thus a “zoom-in” approach was carried out by refining the scale around to get the more precise maximum wavelet coefficient in . Finally, was calculated using Equation (4), and the corresponding signal at the direction of was obtained by Equation (5).
Figure 3 summarizes the procedure discussed above. The procedure of identifying the most unfavorable horizontal direction
was identical to that verified and used by [
1], despite the newly designed procedure of adaptive mother wavelet selection being adopted before projecting the seismic signal into its time–frequency plane. Since a different mother wavelet would surely result in a different pattern of the time–frequency plane [
23], it was necessary to evaluate the level of similarity for different mother wavelets before conducting the CWT to achieve better resolution in the time–frequency plane.
4.2. Representation of Dominant Pulses
The strong ground motion may contain single or multiple pulse-like components featured by the long period and large velocity amplitude. For example, the recorded horizontal seismic signals at station TCU075 during the Chi-Chi event shown in
Figure 4a, are representative of the single one pulse waveform. Its closest distance-to-fault was 3.4 km. The closest distance-to-fault for seismic station NSY, TCU060, and TCU128 were 9.1 km, 8.1 km, and 9.1 km, respectively. The multiple pulse waveform can be easily spotted in the recorded horizontal seismic signals as shown in
Figure 4b–d. Thus, the dominant pulse-like components indeed exist in the recorded seismic signals.
Representation of the dominant multi-pulse waveform embedded in the signal
was carried out by applying the CWT to the residual signal
in an iterative manner. The purpose of CWT was to extract dominant pulses, which was different from stage 1 aimed at identifying
. Before extracting the individual dominant pulses, the optimum mother wavelet was determined first by the adaptive mother wavelet selection procedure, considering the specific feature of the signal
. The largest wavelet coefficient
in the time–frequency plane
produced by the CWT, based on the optimum mother wavelet was used to extract the individual wavelet waveform using Equation (6), where
is the i
th individual pulse corresponding to
, and
is the base value of the mother wavelet adaptively selected by the mother wavelet selection procedure. Due to the varied profile of signal
, the actually adopted mother wavelets were different from each other within the iteration process. The residual seismic signal after the
ith iteration of extraction was updated by
, and it was used as the new target signal for the next iteration.
A criterion for terminating the extraction iteration reasonably was required. For this purpose, two types of criteria with different threshold values were verified with the help of some common pulse-like near-fault seismic signals. One type of criterion was designed as the peak ground velocity (PGV) ratio between and . The other type of criterion was designed as the energy ratio between and . The energy signals and was defined by and , respectively. The energy ratio with a threshold 0.4 was found to behave best compared with the criteria for other threshold values or other types.
The resulting multi-pulse waveform
was obtained by superimposing the extracted individual pulse
by Equation (7), where
is the
ith pulse waveform corresponding to the wavelet coefficient
, and
is the total number of dominant pulses embedded in the signal
. This approach can adaptively extract the varied number of dominant pulse components according to the diverse feature of the target signal
. The individual pulse waveform
, with variant frequency
, was located discretely along the time axis, which demonstrates the non-stationary feature of the near-fault earthquake motion.
After obtaining the multi-pulse component
at the direction of
, the three criteria verified and used by Reference [
1] were adopted in this paper to identify whether the extracted signal
should be qualified as pulse-like motion: (1) pulse index (
), (2) early-arrival of pulse velocity, and (3) omitted small PGV. One signal qualified as pulse-like if the following criteria were reached: (1)
, (2)
, and (3)
.
is the time at which the extracted signal
reaches 10% of its cumulative square velocity (CSV),
is the time at which the original signal
reaches 20% of its CSV, and PGV is the peak ground velocity of the original signal
. The cumulative square velocity is defined by Equation (8), where
is the velocity time–history for
or
. Detailed verification for these criteria can be found in [
1], and is not covered in this paper for brevity.
The extracted parameters for the signal
, such as the number of dominant pulses, pulse periods, and mother wavelets adopted, were collected to characterize these extracted pulse components. Specifically, the parameter N defines the number of significant pulses. The pulse period
of each dominant pulse
embedded in the signal
was calculated by Equation (9), where
is the instantaneous frequency of
in Hz,
is the scale location that corresponds to the
,
is the sampling period of
, and
is the central frequency of the specific mother wavelet selected to conduct the CWT in the i
th iteration.
5. Numerical Example
In this section, we present a numerical example of the identification and representation of multi-pulse near-fault seismic signals by the proposed method. The horizontal velocity time–histories at station TCU101, TCU131, CHY029, and TCU070 during the Chi-Chi event (1999-09-20 17:47:16 UTC, Mw = 7.6) are plotted in
Figure 5. The closest distance-to-fault for these four sets of seismic records were 1.9 km, 26.2 km, 16.4 km, and 18.4 km, respectively. A few prominent pulse waveforms with large amplitude exist in the velocity traces. The signal
and
at station TCU101 reached their peak at the same phase, while this pattern was not so clear for the other three stations. The results from
Section 6 indicate that the signal set at station TCU101 and TCU131 qualified as pulse-like earthquake motions, while the other two sets did not. Although some prominent, large amplitude velocity pulses exist in the two orthogonal channels at stations CHY029 and TCU070, the difference of peak phase may result in the offset of pulse energy when the two orthogonal horizontal records were composited at the specific direction in the horizontal plane. Therefore, it is unreliable to identify pulse-like earthquake motions by single, one channel signals recorded at the seismic station.
The signal set at station TCU101 was used as the numerical example to demonstrate the detailed analysis procedure of the proposed approach. Prior to identifying
, the variance values for all candidate mother wavelets in the repository were evaluated by the adaptive mother wavelet selection procedure. Based on the resulting variance shown in
Figure 6, the mother wavelet ‘gaus3’ was the optimum mother wavelet for the signal set at station TCU101, thus it was used as the mother wavelet to project the signal set into the time–frequency planes
and
. Their resultant plane
shown in
Figure 7, was calculated by the square sum of
and
. In
Figure 7, the x-axis is time, the y-axis is the scales (or frequencies), and the color bar illustrates the amplitude of the wavelet coefficients in the time–frequency plane. The column
and row
for the largest wavelet coefficient
in
were 1558 and 318, respectively. The corresponding wavelet coefficient located at (1558, 318) in
and
were 447.0 and 746.5, respectively. Thus, the most unfavorable direction was identified by
, and the resultant signal
was easily obtained by
as shown in
Figure 8a.
The horizontal direction containing the largest pulse energy for the signals at station TCU101 was 59.09° from north to east, while the strike angle of the Chi-Chi event was 5°. The largest pulse energy theoretically lay in the direction normal to the strike line [
7], at 95°. There was an error ~34.35% between the direction from analysis and theory. This error may result from both the unpredicted distribution of the rock medium between the hypocenter and seismic station, and the unevenly distributed slip direction in the rupture plane.
In order to identify the dominant pulse components embedded in , it was projected by CWT into the time–frequency plane by iteration pattern. The signal length of was 49 s. The y-axis was the transient frequency of the signal converted from the central frequency of the mother wavelet, the scale value, and the sampling frequency.
There were two dominant wavelet components embedded in the signal
. Each component comprised only one wavelet waveform, which was the result of translating and dilating the corresponding mother wavelet by the index
and
of the wavelet coefficient
in the time–frequency plane shown in
Figure 9 as dictated by Equation (6). The signal
was obtained by superimposing the two extracted pulse components:
and
by Equation (7) as shown in
Figure 10. These data show that the signal
captured all the dominant pulses in
. The parameters for these extracted pulse components are listed in
Table 1, where
is the central frequency of the pulse component calculated by Equation (9). The first pulse waveform
was represented by the ‘gaus3’ mother wavelet located at 15.58 s with instantaneous frequency 0.1307 Hz (scale = 306), and the value of PGV is 44.9 cm/s. The second pulse waveform
was represented by the ‘db6’ mother wavelet located at 27.6 s with instantaneous frequency 0.1388 Hz (scale = 524), and the value of PGV was 27.0 cm/s. Due to the the lower frequency and the higher PGV value of the extracted pulse component
compared with that of
, its energy was larger than that of
. This was in accordance with the larger wavelet coefficient for
compared to
.
The criterion of early arrival of the pulse energy was required to ensure that the directivity pulse was of primary concern, as predicted by the theoretical seismology. The pulse index PI was 1.0 for signal
, which was in agreement with the criterion (1):
.
Figure 8 shows the resultant signal
and the pulse signal
embedded in
. The time required by signal
to reach 20% of its CSV was 14.59 s. While the time required by signal
to reach 10% of its CSV was 13.45 s. Thus, the signal
qualified as a pulse early-arriving signal as dictated by criterion (2):
. The PGV of signal
was 0.497 m/s, which was in agreement with criterion (3):
. Finally, the recorded horizontal signals at station TCU101 qualified as multi-pulse near-fault earthquake motions.
Figure 11 shows comparison of the extracted signal by both the author’s and Baker’s methods [
1]. The signal extracted by the author’s method captured all the dominant pulses embedded in the signal
. The signal extracted by Baker’s method only captured the first dominant pulse and failed to capture other dominant pulses after 20 s. Meanwhile, the large disagreement of the signal peak by Baker’s method at 17th second in
Figure 11. This indicates that the author’s method made better representation compared with the Baker’s method, in the situation where the neighboring pulse components lay close to each other. This is because when large amplitude pulse components lie very close to each other, only one pulse is extracted by the Baker’s method, but several large wavelet coefficients are extracted by the author’s method.
The variance results in
Figure 12, by the mother wavelet choosing procedure before extracting
, show that the variance for ‘gaus3’ was 70.0796 and the variance for ‘db4’ was 127.0714. Thus, ‘gaus3’ was the best candidate mother wavelet compared with ‘db4’, which was originally used by Baker. The frequency of the extracted pulse by Baker’s method was 0.098 Hz, 25.02% smaller than the corresponding value for
as listed in
Table 1.
6. Characterization of the Pulse-Like Seismic Motion for Chi-Chi Event
In order to investigate the characterization of the multi-pulse near-fault seismic motion, the author’s method was applied to the 221 sets records for the Chi-Chi event collected from the Strong-Motion Virtual Data Center. Fifty-three sets of records qualified as pulse-like motions with a varied number of dominant pulses.
Figure 13 shows the distribution of the closest distance-to-fault for the 53 sets pulse-like strong ground motions. The cumulative frequency for the records whose closest distance-to-fault was within 30 km was 92.45%, which took up the majority of the qualified records. This result agrees well with the commonly adopted engineering assumption that the near-fault earthquake motions are mainly located within 30 km from the epicenter. Nevertheless, it should be noted that there are still many records located within 30 km from the epicenter that do not qualify as pulse-like, such as the records shown in
Figure 5c,d. Thus, the commonly used engineering assumption of 30 km is not suited for the classification of a specific set of strong ground motion.
For near-fault earthquake motions, the largest pulse energy theoretically lies in the direction normal to the strike line [
7]. The strike angle is 5° for the Chi-Chi event, namely the horizontal direction corresponding to the largest pulse energy in the direction of 95°/275° for the Chi-Chi event. The identified directions containing the largest pulse energy are plotted in
Figure 14 to illustrate the relationship between
and the closest distance-to-fault. In the figure, the theta axis is the value of the parameter
, the radius axis for the closest distance-to-fault is limited to 30 km to facilitate the comparison with the theoretical direction. The dashed line represents the direction that is perpendicular to the strike of the earthquake event. These data show that there are some seismic motions existing in the directions that are deviated from the dashed line, which indicate that the pulse-like seismic motions were not limited to the direction normal to the fault, other directions that were not perpendicular to the fault may still contain prominent pulse-like motions.
Figure 14 also shows that the horizontal direction
was not limited to the EW/NS direction. Thirty-three set records were located within the range between
, which took up 62.26% of the total number of the qualified pulse-like earthquake motions. Only 9.43% of the qualified earthquake motions lay in the direction with an error smaller than 10° with respect to the direction normal to the strike. This meant that the majority of horizontal directions containing the largest pulse energy did not coincide with the EW/NS direction at the seismic station for the Chi-Chi event. The use of single one components recorded at the seismic station (EW/NS) may underestimate the energy level of the input seismic signal. The varied relative orientation of the engineering structure with respect to the identified
may impose a different level of seismic energy on the structure.
The number of dominant pulses embedded in
was a key parameter that greatly influenced the seismic response of structures. The number of dominant pulses was closely related to the non-uniform distributed fracture along the fault, which was difficult to estimate a priori. For the Chi-Chi event, the extraction result using the author’s method showed that 29 sets of earthquake motions were represented by one mother wavelet (listed in
Table 2), 18 sets were represented by two mother wavelets (listed in
Table 3), and six sets were represented by three mother wavelets (listed in
Table 4). No clear trend was found for the relationship between the number of pulses and the closest distance-to-fault.
In
Table 2, the mother wavelets chosen for each station were different from each other. Of the frequencies, 79.31% were located within 0.1–0.2 Hz. In
Table 3, all records (except for TCU068) used different types of mother wavelets to represent the profile of the pulse waveform in the seismic signals. The average frequency for
was 0.128 Hz. The average frequency for
was 0.176 Hz, 37.5% higher than that of
. In
Table 4, the average frequencies for
,
, and
were 0.158 Hz, 0.164 Hz, and 0.227 Hz, respectively.
In
Table 3, the wavelet coefficient for the component
was always larger than that for the component
, indicating that the extracted component
was always the one containing larger pulse energy than component
. A similar pattern was also found in
Table 4. Thus, the individual pulse components extracted using the author’s method were sorted by the pulse energy from high to low.
Figure 15 shows the relative frequency distribution for the occurrence of the mother wavelets in the 53 sets qualified earthquake motions. The relative frequency of ‘db6’ for the first and second pulse components were 28.3% and 37.5%, respectively, which were the most frequently used types of mother wavelets for both the first and the second pulse components as shown in
Figure 15a,b. The ‘morl’ mother wavelet was more frequently used compared with other mother wavelets for the third pulse component as shown in
Figure 15c.
Figure 15 also shows that the ‘haar’ mother wavelet was not used for any pulse extraction because its shape was not similar to the common profile of the pulse waveform embedded in the seismic signals.