Seismic Signal-Based Investigation on the Disaster Differences between 2021 Yangbi Ms 6.4 Earthquake and 2014 Ms 6.5 Ludian Earthquake in Yunnan Province, China

An Ms 6.4 earthquake occurred in Yangbi County, Yunnan Province, on May 21, 2021.e magnitude is only 0.1 dierence when compared with the 2014 Ms 6.5 earthquake in Ludian County, Yunnan Province, but the disasters show a big dierence. Seismic signal includes the comprehensive information of source, propagation path, and site eect, which directly reect the eect of the earthquake on the ground.e study mainly investigated the dierence between seismic signals of two events in time domain and time-frequency-energy space based on variational mode decomposition and Hilbert transform to clarify the reasons for disaster dierences. Most damaged buildings are column-and-tie timber structures with Earth walls with local collapse or vertical crack at the corner of the Earth wall in the Yangbi earthquake. Multistorey buildings are intact, and a few small landslides are triggered, whereas, in the Ludian earthquake, multistorey buildings with the frame structure and masonry structure occurred with serious X-type cracking, large interlayer displacement, and even complete collapse, and large landslides with length about 165m are triggered. 53LLT and 53YBX are two strong motion stations in the Ludian earthquake and Yangbi earthquake, respectively, and are near the disaster investigation points. Peak ground acceleration in 53LLT is up to 949.2 cm/s, which is larger than 720.3 cm/s in 53YBX, and the duration in 53LLT is longer than 53YBX. e instantaneous energy of 53YBX is concentrated within 0–10Hz and 10–20Hz, respectively; however, the instantaneous energy of 53LLT is concentrated at 5Hz. Cumulative energy of 53YBX is mainly distributed at 0–5Hz and 10–20Hz, which is around the natural frequency of column-and-tie timber structure, while the energy of 53LLT is concentrated at 1.6Hz and is much greater than other stations.us, large energy in high frequency caused the instantaneous and cumulative damage to short-period structures in Yangbi County and may cause the resonance of column-andtie timber structure, and large energy in low frequency caused the instantaneous and cumulative damage to long-period structures of Longtou Village in Ludian earthquake. e result is meaningful to understanding the seismic disasters from the characteristics of seismic signals and can provide a reference for the seismic design of structures.


Introduction
Yunnan Province is located on the eastern edge of the collision zone between the Indian Plate and the Eurasian Plate and the southern section of the North-South Seismic Belt of China. ere are 49 active faults developed in the Yunnan Province [1], with strong tectonic movement and signi cant seismic activity. It is one of the areas with earthquakes most active in the world [2,3]. By 2021, 189 earthquakes of Ms 4.6 or higher than Ms 4.6 had occurred in the area, including 49 earthquakes of Ms 4.6-5.0, 107 earthquakes of Ms 5.0-6.0, and 33 earthquakes of Ms 6.0-7.8. According to the earthquake records since the 20th century, earthquakes with M 5.0∼5.9 in Yunnan occur twice a year, earthquakes with M 6.0∼6.9 occur twice every three years, and earthquakes with M 7 and above occur about once every eight years [4]. Large casualties and losses are caused by the frequent occurrence of earthquakes. Particularly in 2014, there were 8 destructive earthquakes with M ≥ 5.0 in Yunnan (including 5 with M 5.0∼5.9 and 3 with M 6.0∼6.9), and the economic losses, deaths, disappearances, and injuries accounted for 84.03%, 99.2%, 100%, and 97.86%, respectively, in the same period at China [5], in which most of casualties and losses are induced by the Ms 6.5 Ludian earthquake. On May 21, 2021, another destructive earthquake with a magnitude Ms 6.4 occurred in Yangbi County, Yunnan Province. e disasters present a big difference when compared with the Ludian earthquake.
Seismic signal is a vibration phenomenon that occurs when energy is suddenly released and propagated to the ground after long-term accumulation, which is the comprehensive reflection of source, propagation path, and site effect [6]. It is the starting point and essential basis of scientific research on seismology and earthquake engineering. At present, we mainly investigate the seismic signal through amplitude, duration, and spectrum. e peak ground acceleration (PGA) and the peak ground velocity (PGV) are usually used as a sign of earthquake intensity [7]. Two factors were also widely used in the structural seismic vulnerability analysis [8,9]. However, PGA value cannot reflect the earthquake damage degree reasonably. For example, in the 1985 Mexico earthquake, the recorded PGA was only 0.17 g, which was much smaller than 0.6 g recorded in the 1986 San Salvador earthquake, but caused more serious earthquake damage [10]. Reference [11] pointed out that PGA represents the amplitude of the high-frequency component of seismic signal and has great discreteness. e experimental and case investigation show that the types of structural damage mainly include instantaneous damage and cumulative damage, and many motions with small amplitude but long duration are very destructive to the structure [12]. erefore, seismic duration has a great impact on the cumulative damage of structures. Amplitude and duration of seismic signals are the characteristics of the time domain. Under the effect of earthquakes, the dynamic response of the structure is closely related to the spectrum characteristics of seismic signal. When the dominant frequency band is close to the natural frequency of the structure, the structural response will be significantly amplified, resulting in serious damage to the structure [13]. If the frequency of the seismic signal is concentrated in the high-frequency band, the short-period structure will cause great harm, and when the predominant frequency of the seismic signal is in the low-frequency band, it will lead to great damage to the long-period structure. In the Peru earthquake on October 17, 1966, the seismic damage in Lima City was mainly concentrated in singlestorey houses, and there was little damage to multistorey houses. According to the seismic signal records of Lima City, the predominant period is nearby 0.1 s, and the natural vibration period of single-storey houses is just close to 0.1 s.
Fourier transform is widely used in the dynamic response of the structure under seismic ground motion [14,15]. However, previous study works show that the Fourier transform is only applicable to linear and stationary data analysis [16,17]. Compared with the Fourier transform, the Hilbert-Huang transform (HHT) can be used for nonlinear and nonstationary signal and has been widely used in analyzing seismic signals in recent years [17][18][19][20]. However, empirical mode decomposition (EMD) of HHT essentially belongs to recursive modal decomposition. e envelope estimation error is amplified by multiple recursive decomposition, which is prone to modal aliasing, and the decomposition result is greatly affected by the sampling rate [21]. Aiming at the problems of EMD, [22] proposed variational mode decomposition (VMD), which is a nonrecursive, selfadaptive, and multiresolution signal decomposition method. VMD overcomes the problems of mode aliasing and error accumulation of EMD [23,24]. VMD combined with Hilbert transform (HT) was used in the study.
We investigated the disasters of the epicenter of the 2021 Yangbi Ms 6.4 earthquake and the 2014 Ms 6.5 Ludian earthquake at first. e disasters show a big difference, even the magnitude with only 0.1 difference. en, we analyzed the characteristics of seismic signals in the two events from the time domain and time-frequency-energy space which is based on VMD and HT. e result of the study can help us know better about the relationship between seismic disasters and the seismic signals and provide a reference for the seismic design of structures.

Earthquake
Events. An Ms 6.4 earthquake occurred in Yangbi County, Dali City, Yunnan Province (Figure 1), on May 21, 2021, with a focal depth of 8 km. It was another destructive earthquake with a magnitude greater than 6 in Yunnan Province after the Ludian Ms 6.5 earthquake, Jinggu Ms 6.6 earthquake, and Yingjiang Ms 6.1 earthquake in 2014. e seismogenic structure of the earthquake is a secondary fault on the west side of the Weixi-Qiaohou fault, which is dominated by dextral strike-slip movement [25]. Before Ms 6.4 earthquake, four earthquakes with a magnitude of Ms 4.2, Ms 5.6, Ms 4.5, and Ms 2.8 occurred successively on the same day. By June 15, 2021, there were more than 40 aftershocks with a magnitude beyond Ms 3.0 in Yangbi County, including 2 earthquakes above Ms 5.0. e largest aftershock was the Ms 5.2 earthquake (25.59°N, 99.97°E) on May 21, 2021 (43 minutes after the mainshock) [26]. erefore, the earthquake sequence is a typical foreshockmainshock-aftershock sequence. e earthquake caused 3 deaths, 32 injuries, 13000 damaged buildings and triggered a few small landslides. e epicenter intensity is up to VIII (China seismic intensity scale). According to the intensity map made by the Yunnan Province earthquake administration (Figure 1), the area of the VI region is about 5500 km 2 , the area of the VII region is 930 km 2 , and the area of the VIII region is 170 km 2 . A total of 46 landslides, 6 collapses, and 7 debris flows were triggered in the earthquake. All of them have small or medium volume and are mainly distributed in the VIII and VII regions [27].
Compared with other earthquakes (Ms > 6) that occurred in Yunnan Province, especially for Ludian Ms 6.5 earthquake, the building damage degree, landslides number, casualties, and economic loss caused by the Yangbi earthquake were much less. Ludian earthquake (Figure 1) occurred on August 3, 2014; the epicenter intensity was up to IX, and the focal depth was about 12 km, which caused 617 deaths, 2400 injuries, 86000 destroyed buildings and triggered many large landslides [28,29]. According to the survey data of the Department of Natural Resources of Yunnan Province, China, 986 landslides and 462 collapses were triggered. Most of them are large and huge landslides. Among them, the Hongshiyan landslide and Ganjiazhai landslide are two largest landslides with a total volume of more than 10 million m 3 . Fifty-five people were buried in the Ganjiazhai landslide.

Topographic Characteristics.
e elevation and slope angle shown in Figures 2 and 3 present the topography characteristics of the area within the intensity map of the Yangbi earthquake and Ludian earthquake, respectively. e elevation within the intensity map of two earthquakes changes from 1070 m to 4180 m and from 434 m to 4003 m, respectively. Figure 4(a) shows that elevation within 1500 m-2500 m in different intensity regions of two earthquakes takes the largest ratio. It takes 72.8%, 67.6%, and 99.3% in the VI region, VII region, and VIII region, respectively, for the Yangbi earthquake and takes 61.6%, 70.0%, 74.7%, and 78.9% in the VI region, VII region, VIII region, and IX region, respectively, for Ludian earthquake. e slope angle within the intensity map of two earthquakes changes from 0°to 71.98°and from 0°to 85.99°, respectively. e slope angle (Figure 4(b)) within 10-20°t akes the largest ratio. Areas with a slope angle larger than 30°t ake 11.5%, 17.9%, and 10.6% in the VI region, VII region, and VIII region, respectively, for the Yangbi earthquake and takes 23.7%, 19.8%, 27.2%, and 30.4% in the VI region, VII region, VIII region, and IX region, respectively, for Ludian earthquake. It shows that the elevation distribution characteristics in different intensity regions of two earthquakes are close, but the slope angle larger than 30°in different intensity regions of the Ludian earthquake takes a higher ratio than that of the Yangbi earthquake.

Characteristics of Disasters.
After the Yangbi Ms 6.4 earthquake, we investigated the damage degree of buildings and geotechnical disasters at Xiajie Village (in Yangbi County), Huaian Village, and Shahe Village, which are within or close to the epicenter and the 53YBX strong motion station, as it is shown in Figure 2(a). Typical seismic damage phenomenon in the Yangbi earthquake is shown in Figure 5. e structural types of buildings in the areas mainly include adobe and timber, brick and timber, and brick and concrete. Based on the investigation, we found that adobe and timber building takes the largest ratio among the damaged buildings. e damage of this structural type mainly occurred at the Earth wall, as it is shown in Figure 5. At Xiajie Village ( Figure 5(a)), the Earth wall of adobe and timber structure mainly made local collapse at corner, while the brick and concrete with multistorey at Xiajie Village is stable, as it is shown in the last two pictures of Figure 5 e damaged adobe and timber buildings in the Yangbi earthquake belong to traditional Chinese column-and-tie timber structures, as shown in Figure 6. e upper load and weight of the whole roof are transmitted from rafters to purlins, then from purlins to column, and through column to the ground. Chuanfang plays an important role in connecting wooden columns and enhancing the stability of the whole frame and has little contribution to load bearing. e column is the main load-bearing member. e construction sequence of the building is as follows: the wooden column is built at first, and then the Chuanfang is used to connect the columns to establish the timber frame. On this basis, the Earth wall is built. Finally, purlin is installed, and rafters and ties are erected to complete the main construction of the building. e seismic behavior of it is poor because of lacking effective earthquake resistance measures. Its natural frequency is around 10 Hz [30].
Compared with the Yangbi earthquake, the damage phenomenon induced by the Ludian earthquake presents different characteristics. Figure 7 shows the typical seismic Shock and Vibration damage phenomenon in Longtou Village, which is near the 53LLT strong motion station, as it is shown in Figure 3(a). e structural types of buildings in Longtou Village mainly include frame structure, masonry structure, timber structure, and adobe and timber. In Longtou Village, frame structure and masonry structure with multistorey are seriously damaged. Figure 7(a) shows that serious X-type cracking and large interlayer displacement occurred at these buildings, and many of them even collapsed completely. Large landslides with a length of about 165 m also occurred at Longtou Village, as shown in Figure 7(b), which destroyed and buried part of the buildings. It shows that seismic damage caused by the Ludian earthquake is much more serious than the Yangbi earthquake, although the two events all occurred in Yunnan Province, and the magnitude is only 0.1 difference.

Strong Motion Data.
In the Yangbi earthquake, 17 strong motion stations distributed within the intensity map record this event, including 1 in the VIII region, 1 in the VII    Shock and Vibration region, and 15 in the VI region. In the Ludian earthquake, 7 strong motion stations distribute in the intensity map, including 1 in the IX region, 1 in the VII region, and 5 in the VI region. Table 1 shows the basic information of the strong motions with the epicentral distance within 40 km of two earthquakes. e 53YBX and 53LLT with an epicentral distance of around 8 km are near the disaster investigation points of two earthquakes, which give a valuable opportunity to study the disaster differences according to the characteristics of seismic signals. e baseline correction and Butterworth filter with 0.1-30 Hz bandwidth were applied to these strong motion recordings.

Variational Mode Decomposition (VMD).
e Intrinsic Mode Function (IMF) is defined as an amplitudemodulated-frequency-modulated signal, as shown in VMD is a process of solving the variational problem, which can minimize the sum of the estimated bandwidth of each mode. In order to estimate the bandwidth of modes, firstly, the unilateral spectrum is obtained by the Hilbert transform. en, the modal spectrum is modulated to the fundamental frequency band by mixing the central frequency. Finally, the demodulated signal is processed by Gaussian smoothing. e constrained variational problem for estimating the bandwidth of modes is shown in the following equation: where u k � u 1 , u 2 , ..., u k are the IMFs, ω k � ω 1 , ω 2 , ..., ω k are the central frequency of each IMF, and δ(t) is the pulse function.   Shock and Vibration 7 e quadratic penalty term (α) and Lagrange multipliers (λ(t)) are introduced to solve the problem. e augmented Lagrangian function is as follows: In order to solve equation (3), Alternating Direction Method of Multipliers (ADMM) is used by iteratively updating u k , ω k , and λ. e specific steps are as follows: (1) u 1 k , ω 1 k , λ 1 , and n are initialized to 0; (2) u k , ω k , and λ are iteratively updated by equations (4)-(6), respectively; (3) repeat steps (2) to (3) until the iteration termination conditions (equation (7)) are met: where ε is the accuracy, ε > 0.
e number of IMFs is set to 3 in the study. It is taken according to the number of predominant frequencies of the Fourier spectrum. Figures 8 and 9 are the VMD result of 53YBX-EW and its Fourier spectrum, respectively. e Fourier spectrum of the original signal and its three IMFs indicates that the three predominant frequencies are around 2 Hz, 10 Hz, and 15 Hz, respectively.

Hilbert Transform (HT)
. IMF is acquired by VMD, and each of them is well behaved in the Hilbert transform.
e Hilbert transform y(t) is shown in the following equation [17]: where x(t) is the time history of the signal, which is equal to u(t) in this study, and P is the Cauchy principal value. An analytic signal z(t) can be defined as e calculation method of a(t) and θ(t) is shown in equations (10) and (11), respectively: e instantaneous frequency (ω(t)) is the first derivative of θ(t) to time equation.
e Hilbert spectrum H(t, ω) is shown in the following equation:

Shock and Vibration
e marginal spectrum h(w) is the integral of H(t, ω) over time, as shown in the following equation:

Time Domain Characteristics.
Major characteristics such as PGA and duration of seismic wave record by the strong motion station can be acquired in the time domain. e acceleration time history of 53YBX and 53LLT is shown in Figure 10. 53YBX station is 8.6 km away from the epicenter, and the PGAs in three directions were 379.8 cm/s 2 (EW), 720.3 cm/s 2 (NS), and 448.4 cm/s 2 (UD), respectively. 53LLT station is 8.3 km away from the epicenter, and the PGAs in three directions were 949.2 cm/s 2 (EW), −705.8 cm/ s 2 (NS), and −504.4 cm/s 2 (UD), respectively. PGA attenuation rules of the two events are shown in Figure 11. In soil site, with the increase of epicentral distance, PGA in three directions has a decreasing tendency. In rock sites, PGA has a certain degree of increase. We should pay attention to the fact that PGA increases obviously at QJX, which may relate to the local landform. e duration of an earthquake can reflect the cumulative damage effect of seismic signals on structures. ere are two kinds of seismic signal duration, which are absolute duration and relative duration. In the earthquake engineering field, we mainly study the relative duration, which focuses on the strong part of seismic signals. Kawashima bracketed duration, which is a kind of relative duration, was used in the study to investigate the duration characteristics. It takes the duration of the first and last times reaching or exceeding 1/ 5∼1/2 of the maximum peak acceleration as the duration, also known as fractional duration or relative peak duration. In this study, 0.3 is selected for calculation and analysis. e duration of strong motions of two earthquakes is shown in Table 2. e distribution of the duration is relatively discrete. In the Yangbi earthquake, the minimum duration is 2.1 s (53YBX-NS), and the maximum duration is 13.4 s (L3003-UD). In the Ludian earthquake, the minimum duration is 4.3 s (53LLT-UD), and the maximum duration is 15.3 s (53LDC-NS). e maximum PGA of two earthquakes occurred at 53YBX-NS and 53LLT-EW, respectively, and the duration of 53YBX-NS is less than 53LLT-EW, which indicates that the seismic signal of 53LLT has much cumulative damage effect on structures.

Time-Frequency-Energy Space Characteristics.
Based on VMD and HT, the time-frequency-energy distribution characteristics of seismic waves can be obtained, which are mainly presented in the Hilbert spectrum and marginal spectrum in the study.
e Hilbert spectrum reflects the distribution of time, frequency, and instantaneous energy. e marginal spectrum is the integral of the Hilbert spectrum over time, which can clarify the cumulative energy of every instantaneous frequency in the whole duration.

Hilbert Spectrum.
e acceleration time history recorded by strong motion station has three components, and each of them contains meaningful information on seismic signals. In the time domain, it is not easy to synthesize the three components for the amplitude depending on the direction. Most studies mainly investigated the single component and ignored the others. However, in time-frequency-energy space, three components all present the distribution of time, frequency, and energy which are scalar. e synthesized Hilbert spectrum can reflect the instantaneous energy distribution of the seismic wave of a strong motion station completely. e Hilbert spectrums of 53YBX in three directions are shown in Figure 12. In the EW direction (Figure 12(a)), the energy is mainly distributed at 0-10 Hz and 10-20 Hz, respectively. However, in NS and UD directions (Figures 12(b)-12(c)), energy is mainly distributed at 10-20 Hz. e synthesized Hilbert spectrum of   Figure 13(a); two predominant frequencies are clearly clarified. Hilbert spectrums of seismic waves recorded by different stations in the Yangbi earthquake are shown in Figures 13(a)-13(e). e shape of 53YBX presents doublepeak, which are within 0-10 Hz and 10-20 Hz, respectively. However, the shape of others presents a single peak, and the instantaneous energy is mainly concentrated within 20 Hz. With the increase of epicentral distance, the energy shows a downward trend. However, the peak instantaneous energy has a slight increase in L3003 (from 52 cm/s 2 in L2201 to 73 cm/s 2 in L3003), which indicates that the soil layer can absorb peak instantaneous energy. Hilbert spectrums of seismic waves recorded by different stations in the Ludian earthquake are shown in Figures 13(f)-13(j). e shape of them shows a single peak, and the energy is mainly concentrated within 20 Hz. e peak instantaneous energy of them occurred within 10 Hz. e energy also shows a downward trend with the increases of epicentral distance, and the peak energy at the rock site (53HYC) is greater than that at the soil site when the epicentral distance is close. It should be noted that the energy at 53QJX has a significant increase in 53QJX (from 50 cm/s 2 in LDC to 130 cm/s 2 in 53QJX), which may be caused by its special site conditions and need to be investigated furtherly in our later works. e energy of 53LLT is concentrated at 5 Hz and the peak up to 1280 cm/s 2 , which has a significant destructive instantaneous effect on long-period buildings. Compared with 53LLT, 53YBX contains relatively large energy in highfrequency (>10 Hz) parts, which has a destructive instantaneous effect on the short-period structures.

Marginal Spectrum.
Marginal spectrums of seismic waves in the Yangbi earthquake are shown in Figure 14(a). In 53YBX, cumulative energy is mainly distributed in two frequency ranges, which are 0-5 Hz and 10-20 Hz, respectively. In other stations, cumulative energy is mainly distributed within 10 Hz. Peak cumulative energy shows a decreasing trend with the increases of epicentral distance, which is independent of site conditions. Marginal spectrums in the Ludian earthquake are shown in Figure 14(b). In 53LLT, the predominant frequency is 1.6 Hz, and the corresponding amplitude is much greater than it in other stations. e predominant frequency of other stations is within 3.8 Hz-4.3 Hz. Peak cumulative energy also shows a decreasing trend with the increases of epicentral distance. However, at 53QJX, cumulative energy shows an abnormal increase. is phenomenon also shows the Hilbert spectrum. Comparing the marginal spectrum of 53YBX with 53LLT, the difference is obvious. For 53LLT, large cumulative energy is concentrated at low frequency, which can induce the cumulative damage of longperiod structures. For 53YBX, a large part of cumulative energy   is concentrated at a high frequency, which can cause cumulative damage to the short-period structures.
According to the distribution characteristics of cumulative energy shown in Figure 14, we mainly investigate the attenuation rules of cumulative energy in three frequency ranges which are 0-5 Hz, 5-10 Hz, and 10-15 Hz, as it is shown in Figure 15. In three frequency ranges, with the increase of epicentral distance, the cumulative energy of the Yangbi earthquake decays faster than the energy of the Ludian earthquake. e energy of the Ludian earthquake is greater than the Yangbi earthquake when the epicentral distance is close. e result reflects that the seismic wave of the Ludian earthquake contains larger cumulative energy than the Yangbi earthquake and attenuates more slowly, which can cause more serious damage. It should be noticed that the energy in 10-15 Hz of 53YBX is even larger than 53LLT, although there is a 0.1 magnitude of difference between the two earthquakes, and the epicentral distance of the two stations is close. As mentioned above, the natural frequency of damaged buildings (column-and-tie timber structure with Earth wall) in the Yangbi earthquake is around 10 Hz. Large cumulative energy in high frequency (10-15 Hz) may cause the resonance of the structure. So, short-period structures in the Yangbi earthquake suffered serious cumulative damage.

Effect of Frequency of Seismic Wave on Slope and Building.
Disasters caused by two events show a big difference. In the Yangbi earthquake, the casualties and economic loss are minor. e most damaged building type is adobe and timber with single-storey and with local collapse or vertical crack occurring at the corner of retaining wall. ere is no obvious damage to multistorey buildings. However, the epicenter of the Ludian earthquake suffered great damage. Frame structure and masonry structure with multistorey caused serious X-type cracking and large interlayer displacement and even complete collapse. Meanwhile, large landslides with a length of 165 m were triggered. e effect of the frequency of a seismic wave on the movement of the slope is shown in Figure 16. In theory, a high-frequency seismic signal with a shorter wavelength than the dimensions of slopes is unlikely to trigger landslides because high-frequency seismic waves drive different parts of the slope in the opposite direction simultaneously, so the whole slope would remain stable. On the other hand, lowfrequency seismic wave with a longer wavelength than dimensions of slopes can lead slope to move in the same direction, which is more likely to induce landslides [31]. Similarly, the frequency of seismic wave has the same effect on the building (Figure 17). e low-frequency seismic wave would lead multistorey building to move in the same direction, and the building easily collapsed. On the contrary, the building would remain stable under the effect of high frequency of the seismic wave. It is in accord with the damage phenomenon of slope and buildings in Yangbi County and Longtou Village. is study shows that 53LLT contains large instantaneous energy and cumulative energy in low-frequency parts; however, 53YBX contains large energy in high-frequency parts. us, the large energy in high frequency caused instantaneous and cumulative damage to short-period structures in Yangbi County. Accordingly, the large energy in low frequency caused instantaneous and cumulative damage to long-period structures (multistorey buildings) and large landslides in Longtou Village.

Site Effect Characteristics.
In order to investigate the site effect in two events, we use the horizontal to vertical Fourier spectral ratio (HVSR) of seismic signals to evaluate it. e method is proposed by Nakamura [32][33][34] and has developed into a standard technique to study site effects [35,36]. e expression is shown in equation (15). In the study, the Fourier spectrum of the horizontal seismic signal is the synthetic seismic signal of EW and NS. Similarly, we calculated the horizontal to vertical marginal spectral based on the VMD and Hilbert transform (VHT) ratio of seismic signals in order to compare with the ratio based on the fast Fourier transform (FFT).
HVSR results based on FFT and VHT, respectively, are shown in Figure 18. e change rules of the FFT ratio and VHT ratio of different stations in two earthquakes are close. Two methods can recognize the predominant frequency of the sites, and the predominant frequency value is close. However, the amplitude has a certain difference, and the predominant frequency recognized by the VHT ratio is clearer than that by the FFT ratio. Peak amplitude (with FFT ratio 26 and VHT ratio 49) occurred at 53QJX, and the predominant frequency is about 5 Hz, which indicates that relatively large instantaneous energy and cumulative energy here are related to the local site conditions. e amplification effects (amplitude) of 53YBX and 53LLT are not obvious compared with other stations, which reflects that the distribution of instantaneous energy and cumulative

Multi-storey building
Multi-storey building Single-storey building Single-storey building where H(f) is the Fourier spectrum of the horizontal seismic signal and V(f) is the Fourier spectrum of the vertical seismic signal.

Low-Frequency Pulse Signal Extraction Based on VMD
Method. Many studies have been done on the reasons for the large disaster in Longtou town caused by the Ludian earthquake. One main reason among them is the largeamplitude pulse contained in the time history of 53LLT [37][38][39]. Pulse-type seismic signal has more complex characteristics than ordinary seismic signal and can lead to serious damage to engineering structures [40,41]. e Port Hueneme Mw 4.7 earthquake in 1957, the Parkfield Mw 5.6 earthquake in 1966, and the San Fernando Mw 6.6 earthquake in 1971 all showed the characteristics of small magnitude and serious disasters, and pulse-type seismic signals were recorded in these events [42,43]. So, there is a need to identify, extract, and model the pulse-type seismic signal for further structural seismic analysis. However, the identification method requires visual inspection or repeat using wavelet transform [44,45], and the method is based on the velocity time history of the seismic signal. In this study, an identification or extraction method using acceleration time history is proposed.   Figure 19(a). e result shows that 53LLT (EW) contains a short-term velocity pulse of around 25 s, and IMF1 is the major component of the short-term velocity pulse of the original signal. e method based on VMD can effectively extract the pulse components from the acceleration time history of the seismic signals.

Conclusions
An Ms 6.4 earthquake occurred in Yangbi County, Yunnan Province, on May 21, 2021. Compared with Ms 6.5 Ludian earthquake that occurred in Yunnan Province on 2014, the disasters have a big difference, even the magnitude with 0.1 difference. e study mainly introduced the difference in seismic signals of epicenters of two earthquakes. e main conclusions are shown as follows. e Yangbi Ms 6.4 earthquake caused 3 deaths, 32 injuries, 13000 damaged buildings and triggered a few small landslides. e epicenter intensity is VIII. Column-and-tie timber structure with Earth wall caused local collapse or vertical crack at the corner of Earth wall, which takes the largest ratio among the damaged buildings at the epicenter, and the multistorey buildings do not have obvious damage phenomenon.
e Ludian Ms 6.5 earthquake caused 617 deaths, 2400 injuries, 86000 destroyed buildings and triggered many large landslides. e epicenter intensity is up to IX. Frame structure and masonry structure with multistorey in epicenter caused serious X-type cracking, large interlayer displacement, and even complete collapse.
PGAs of 53YBX in three directions are 379.8 cm/s 2 (EW), 720.3 cm/s 2 (NS), and 448.4 cm/s 2 (UD), and the PGAs of 53LLT in three directions are 949.2 cm/s 2 (EW), −705.8 cm/s 2 (NS), and −504.4 cm/s 2 (UD). PGA in the Ludian earthquake attenuates faster than that in the Yangbi earthquake with the increases of epicentral distance, and the duration in 53LLT is longer than 53YBX. e large PGA and long duration of the seismic signal may cause serious damage to the building in Longtou Village. e Hilbert spectrum shape of 53YBX presents a double peak, and the instantaneous energy is concentrated within 0-10 Hz and 10-20 Hz, respectively. e Hilbert spectrum shape in the Ludian earthquake shows a single peak. e instantaneous energy of 53LLT is concentrated at 5 Hz, and the peak is up to 1280 cm/s 2 , which may have a significant destructive instantaneous effect on long-period buildings. e marginal spectrum shows that, in 53YBX, cumulative energy is mainly distributed in two frequency ranges, which are 0-5 Hz and 10-20 Hz, respectively. In 53LLT, the predominant frequency is 1.6 Hz, and the cumulative energy is much greater than that in other stations. e large energy in low frequency caused instantaneous and cumulative damage to long-period structures in Longtou Village. e large energy in high frequency caused instantaneous and cumulative damage to short-period structures in Yangbi County and may cause the resonance of column-and-tie timber structure.

Data Availability
Basic data can be obtained from the corresponding author when needed. Data for this study are provided by the Institute of Engineering Mechanics, China Earthquake Administration.

Conflicts of Interest
e authors declare that they have no conflicts of interest.