Reconstructing the Dynamic Processes of the Taimali Landslide in Taiwan Using the Waveform Inversion Method

: As a landslide occurs, seismic signals generated by the mass sliding on the slope can be recorded by seismometers nearby. Using waveform inversion techniques, we can explore the dynamic processes (e.g., sliding direction, velocity, and runout distance) of a landslide with the inverted force–time function. In this study, the point force history (PFH) inversion method was applied to the Taimali landslide in Taiwan, which was triggered by a heavy rainstorm in 2009. The inverted force–time function for the landslide revealed the complicated dynamic processes. The time series of velocity indicated three di ﬀ erent sliding directions during the landslide. Hence, three propagating stages of the Taimali landslide were determined and were consistent with an investigation using remote sensing images and a digital elevation model of the landslide. In addition, the PFH inversion was implemented using high-quality single-station records and maintained good performance compared with the inversion by multistation records.


Introduction
Landslides caused by either earthquakes or rainfall have caused the destruction of public facilities, economic loss, and loss of life and property in Taiwan over the past few decades [1,2]. One of the severest slope disasters in Taiwan took place during Typhoon Morakot, which affected Taiwan from 8-10 August 2009. That disaster left nearly 700 persons dead or missing and caused about 3.3 billion USD in economic loss [2,3]. The unexpected heavy rainfall caused thousands of landslides and debris flows in southern Taiwan, including the giant Hsiaolin landslide and the Taimali landslide [2]. In particular, significantly more landslides with larger disturbed areas occurred during Typhoon Morakot in 2009 than during previous typhoon events. These giant landslides, known as large-scale landslides (LSLs), are defined as the landslides with areas of larger than 10 5 m 2 , volumes of more than 10 6 m 3 , or depths greater than 10 m [4]. After Typhoon Morakot, the study of LSLs became one of the most important issues related to slope disasters in Taiwan. The features of LSLs include spontaneity, fast movement, and destruction; hence, landslide processes are a hard-to-observe natural phenomenon [5], and it is difficult to mitigate landslide disasters with traditional study methods.
A few studies have reported that block sliding processes on the slope would generate ground motions that can be recorded by seismic stations [5][6][7][8]. The landslide-generated ground motions contain low-frequency (<1 Hz) and high-frequency (>1 Hz) waves. The high-frequency waves are produced by the unloading and reloading cycle of the sliding materials [7], and the high-frequency waves result from the complex impacts between the sliding debris [8]. In the frequency domain, The Taimali landslide is located on the Eocene Bilusan Formation composed of phyllite and slate with sandstone interbeds [20]. Folds and cleavages are also well developed in this region and hence lead to the fractured geologic condition of the Taimali landslide. According to the Central Geological Survey [21], a region-scale anticlinorium lies to the west of the Taimali landslide and is seven kilometers from the landslide. The source area of the Taimali landslide was on a dip slope and prone to collapse during heavy rainfall.
Pre-and post-event digital elevation models (DEMs) with resolutions of five meters were used to map the disturbed area of the landslide and to calculate the variance in topography ( Figure 2). The main source area (Zone 1) was about 9.3 × 10 5 m 2 , and its total volume was about 7.5 × 10 7 m 3 . The maximum depth, roughly at the center of the landslide block, was 220 m. Based on the variance in topography, two deposit areas (Zone 2 and 3) could be recognized. The thickness of the landslide deposit in Zone 2 was thicker than that in Zone 3 (Table 1). This difference indicated that Zone 2 might be the main deposit area for this event. A relatively smaller landslide (3 × 10 5 m 2 ) occurred and provided some amount of debris to the deposit area (Zone 2 and Zone 3). Because the volume of the smaller source was not considered, the sum of the volume of the two major deposit areas is greater than the volume of the only major source area. Based on the average density of 2700 kg/m 3 for phyllite, the mass of the Taimali landslide was estimated at 2 × 10 11 kg, and therefore it generated high-energy seismic signals.

Seismic Data Processing
In this study, the seismic data recorded by broadband seismic stations were collected from the Broadband Array in Taiwan for Seismology (BATS), which was deployed by Academia Sinica and Taiwan's Central Weather Bureau (CWB) since 1992 [22]. BATS provides high-precision continuous seismic records to document the ground motions of landslides. All permanent BATS stations are equipped with broadband sensors with a frequency band of 0.00833-8 Hz. All stations are capable of internet connection for immediate retrieval of real-time data. The high-energy seismic signals generated by the Taimali landslide were detected by 15 three-component broadband seismometers (Figure 3), including the nearest station, MASB, at a distance of 19 km, and the farthest station, KMNB, at a distance of 325 km. Even in the Taiwan Strait, the seismic signals caused by the landslide could be detected by the PHUB and KMNB stations on small islands.

Seismic Data Processing
In this study, the seismic data recorded by broadband seismic stations were collected from the Broadband Array in Taiwan for Seismology (BATS), which was deployed by Academia Sinica and Taiwan's Central Weather Bureau (CWB) since 1992 [22]. BATS provides high-precision continuous seismic records to document the ground motions of landslides. All permanent BATS stations are equipped with broadband sensors with a frequency band of 0.00833-8 Hz. All stations are capable of internet connection for immediate retrieval of real-time data. The high-energy seismic signals generated by the Taimali landslide were detected by 15 three-component broadband seismometers (Figure 3), including the nearest station, MASB, at a distance of 19 km, and the farthest station, KMNB, at a distance of 325 km. Even in the Taiwan Strait, the seismic signals caused by the landslide could be detected by the PHUB and KMNB stations on small islands.
The ground motion induced by the Taimali landslide can be observed in the seismic record at the MASB station ( Figure 4a). However, it is not easy to determine the initiation time and end time of an event by observing complex time-domain waveforms. Therefore, a time-frequency spectrogram and a spectrum were obtained through the fast Fourier transform (Figure 4b,c). In the seismic record, the event initiated after 30 s and terminated at around the 210th second. From the spectrum, the presence of distinct low-frequency and high-frequency components was apparent. The values f = 3.0 Hz and f = 0.5 Hz were taken as the representative frequency of the signal in the high-frequency and low-frequency ranges, respectively. Several wave packets can be observed in the two filtered waveforms ( Figure 4d,e), which indicates the event experienced several stages of movement. Lin et al. [8] proposed an empirical formula to evaluate landslide seismic magnitude (L m ) and reported that the L m of the Taimali landslide was 4.3.
The seismic data were processed by deconvolving the instrument response, integrating the seismic records from ground velocity to displacement, and resampling seismic data from 20 to 1 Hz. Then the horizontal components of the seismograms of each station were rotated from the north and east directions to the radial and transverse directions. The radial direction refers to the direction from the seismic source to the station, and the transverse direction is perpendicular to the radial direction. The processed signals were filtered by a fourth-order Butterworth band-pass filter with a corner frequency of 0.01-0.033 Hz. Finally, the seismic data were weighted according to the signal-to-noise ratio (SNR). SNR was calculated by dividing the maximum of the absolute amplitude by the mean of the absolute amplitude of the seismic waveform.  The ground motion induced by the Taimali landslide can be observed in the seismic record at the MASB station ( Figure 4a). However, it is not easy to determine the initiation time and end time of an event by observing complex time-domain waveforms. Therefore, a time-frequency spectrogram and a spectrum were obtained through the fast Fourier transform (Figure 4b,c). In the seismic record, the event initiated after 30 s and terminated at around the 210th second. From the spectrum, the presence of distinct low-frequency and high-frequency components was apparent. The values f = 3.0 The processed signals were filtered by a fourth-order Butterworth band-pass filter with a corner frequency of 0.01-0.033 Hz. Finally, the seismic data were weighted according to the signal-to-noise ratio (SNR). SNR was calculated by dividing the maximum of the absolute amplitude by the mean of the absolute amplitude of the seismic waveform.

Point Force History Inversion
The point force history (PFH) inversion method was proposed by Allstadt [17]. Even if the landslide material gradually breaks up during the movement, it can still be regarded as a single point under the large terrain scale. Complicated and local dynamic characteristics changes in collapsed materials are not easily reflected in distant seismic records at large spatial scales. Based on a hypothesis that a landslide can be simplified as a block of constant mass m sliding on a slope, the slope parallel force ∥ from the driving force of gravity opposed by the friction force can be expressed as where g is the gravity acceleration, θ is the slope angle, and μ' is the dynamic frictional force. If the sliding force is larger than the frictional force, the landslide block will be unstable and start to slide. Based on Newton's second law, Equation (1) can be modified as

Point Force History Inversion
The point force history (PFH) inversion method was proposed by Allstadt [17]. Even if the landslide material gradually breaks up during the movement, it can still be regarded as a single point under the large terrain scale. Complicated and local dynamic characteristics changes in collapsed materials are not easily reflected in distant seismic records at large spatial scales. Based on a hypothesis that a landslide can be simplified as a block of constant mass m sliding on a slope, the slope parallel force F from the driving force of gravity opposed by the friction force can be expressed as where g is the gravity acceleration, θ is the slope angle, and µ is the dynamic frictional force. If the sliding force is larger than the frictional force, the landslide block will be unstable and start to slide. Based on Newton's second law, Equation (1) can be modified as According to Newton's third law, when body A exerts a force on body B, body B simultaneously exerts a force equal in magnitude and opposite in direction on body A. Therefore, as the landslide block receives a force downslope, the sliding surface receives an equal point force Fe in the opposite direction. This force Fe is a time-dependent three-component vector that can be expressed as The long-period seismic signals generated by landslides are attributed to the single force Fe [10,17,23]. Equation (3) dictates that the direction of the force exerted on the sliding surface will be in the opposite direction to the landslide acceleration. That is, as the landslide accelerates downslope, the direction of the force heads upslope; as the landslide decelerates, the sliding surface receives a force in the downslope direction.
After the instrument responses are removed, the seismograms recorded by seismometers represent the results of the seismic sources that pass through radiation paths. The responses of an impulse force between each source and station are called Green's functions, and they consider all seismic waves and attenuation on the radiation paths [17,24]. If the velocity structure of the material through which the seismic waves are passing is known, Green's functions can be calculated. Due to the long-period band-pass filter (30-100 s) used in this study, the seismic waveforms radiated by the landslide had wavelengths of tens to hundreds of kilometers and had low sensitivity to small-scale heterogeneities in the velocity structure and regional topography. Hence, a generalized one-dimension velocity model is sufficient to calculate Green's functions. In this study, the Taiwan one-dimension velocity model proposed by Chen and Shing [25] was used. At the top 2 km of the velocity structure, the velocity of the P wave reaches 4.64 km/s with a velocity gradient of 0.37 km/s per km, and the velocity of the S wave is 2.64 km/s. The velocity structure was derived from P and S wave travel time and was confirmed that it can indicate the thickness of recent sediment on the Taiwan island.
The Computer Program in Seismology (CPS) software was used to calculate the Green's functions between each station and the landslide location. For the single-force mechanism of the landslide processes suggested by Kanamori and Given [7,8], five Green's functions are required [24]. The five Green's functions are listed as follows: From the waveform perspective, the first four components correspond to the velocity system of P wave and S wave, while THF corresponds to the SH system. Herrmann [22] described the displacement seismograms for any point force time series by where u z , u r , and u t are the synthetic ground displacement in the vertical, radial, and transverse directions, and φ is the source to station azimuth measured clockwise from the north. The single force F = (f N , f E , f Z ) is in the north, east, and vertical direction, respectively. The displacement seismograms are defined as positive in the following directions: up, away from the source, and right angle clockwise in the vertical, radial, and transverse components, respectively. The displacement seismogram d i (t) at each component out of the i station is the result of the force-time function f j (t) in the j direction convoluted (*) with the Green's functions G ij (t), which can be expressed as If the displacement seismogram and Green's functions are available, the landslide force-time function can be extracted by the waveform inversion. To execute the inversion, Equation (7) should be written as a multiplication matrix rather than a convolution. A convolution is defined as the integral of the product of two functions after one is reversed and shifted. To convolve a Green's function with a force-time series through matrix multiplication, the Green's functions were shifted by one sample in each successive row to generate a convolution matrix. An example of the convolution through matrix multiplication between a four-sample Green's function g and a four-sample force-time function f to obtain a seven-sample long seismogram d can be written as Equation (8) where dz, dr, and dt are the column vectors that contain the seismograms for each station at the vertical (z), the radial (r), and the transverse (t) components, respectively. The 3 × 3 matrix (G matrix) consists of the five Green's functions and the azimuth calculation. Each element of the G matrix is a Green's function convolution matrix for each station.
Since large-scale landslides often occur during typhoons, the seismic records can be affected by high-level ambient noise. Thus, different weighting coefficients were given to the seismograms according to their SNR values (Table 2). Once the weighting coefficients are determined, the complete forward equation is established and the inversion procedure can be conducted. A least-square sense is applied to solve the landslide force-time function f(t) by where G w is the weighted Green's function, d w is the weighted seismograms, and superscript T means the transpose. With Equation (10), it is easy to obtain a landslide force-time function in MATLAB on a computer.

Results of PFH Inversion
The inversion results showed that the force induced by the Taimali landslide was initiated at about 09:31:34 and terminated at 09:34:28 ( Figure 5). Among the three components of force, the largest force existed in the E-W direction, reaching 3.7 × 10 10 N. The maximum force in the N-S direction and the maximum vertical forces were 1.6 × 10 10 and 1.2 × 10 10 N, respectively. f = (Gw T Gw) -1 Gw T dw (10) where Gw is the weighted Green's function, dw is the weighted seismograms, and superscript T means the transpose. With Equation (10), it is easy to obtain a landslide force-time function in MATLAB on a computer.

Results of PFH Inversion
The inversion results showed that the force induced by the Taimali landslide was initiated at about 09:31:34 and terminated at 09:34:28 ( Figure 5). Among the three components of force, the largest force existed in the E-W direction, reaching 3.7 × 10 10 N. The maximum force in the N-S direction and the maximum vertical forces were 1.6 × 10 10 and 1.2 × 10 10 N, respectively. After the force-time function was obtained, the inverted force was applied to get the synthetic waveforms of each seismic station ( Figure 6). Most synthetic waveforms were consistent with the station records: 47% of the synthetic waveforms had correlation coefficients (cc) values higher than 0.9, and 70% of the synthetic waveforms had cc values higher than 0.7. They both indicated that the inversion results were robust. The ANPB and CHGB stations both had strong ambient noise, which led to the low SNR and cc values. In addition, the poor results of the KMNB and PHUB stations might be attributable to the velocity model used in this study. Despite the poor results caused by the strong ambient noise and the influence of the velocity model, the high cc values suggested that the inversion had good performance and the results were highly reliable.
Based on the analysis of satellite imagery and the results of waveform inversion, the Taimali landslide slid toward the northeast at the beginning, which implied that the force would first head toward the southwest. However, the first pulse of force from 09:31:21 to 09:31:34 was toward the northeast. Thus, the first pulse of force was considered a precursory signal caused by Gibbs effects [7]. The landslide propagation began at 09:31:34, at which time the force headed toward the southwest. The force after the 250th second was considered fluctuations of ambient noise, not the propagation of the Taimali landslide. The 174-s duration of the landslide was consistent with the waveform observation ( Figure 6); therefore, the force-time function indicated that the Taimali landslide initiated at 09:31:34 and terminated at 09:34:28. After the force-time function was obtained, the inverted force was applied to get the synthetic waveforms of each seismic station ( Figure 6). Most synthetic waveforms were consistent with the station records: 47% of the synthetic waveforms had correlation coefficients (cc) values higher than 0.9, and 70% of the synthetic waveforms had cc values higher than 0.7. They both indicated that the inversion results were robust. The ANPB and CHGB stations both had strong ambient noise, which led to the low SNR and cc values. In addition, the poor results of the KMNB and PHUB stations might be attributable to the velocity model used in this study. Despite the poor results caused by the strong ambient noise and the influence of the velocity model, the high cc values suggested that the inversion had good performance and the results were highly reliable.
Based on the analysis of satellite imagery and the results of waveform inversion, the Taimali landslide slid toward the northeast at the beginning, which implied that the force would first head toward the southwest. However, the first pulse of force from 09:31:21 to 09:31:34 was toward the northeast. Thus, the first pulse of force was considered a precursory signal caused by Gibbs effects [7]. The landslide propagation began at 09:31:34, at which time the force headed toward the southwest. The force after the 250th second was considered fluctuations of ambient noise, not the propagation of the Taimali landslide. The 174-s duration of the landslide was consistent with the waveform observation ( Figure 6); therefore, the force-time function indicated that the Taimali landslide initiated at 09:31:34 and terminated at 09:34:28.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 10 of 16 Figure 6. Three components of the synthetic waveforms and the seismograms of each station. Z denots vertical direction, R denotes radial direction, and T denotes tangential direction. Correlation coefficients are represented by "cc". Light-colored seismograms represent that the lower weighted coefficients were adopted to conduct the PFH inversion.

Dynamic Processes Determined by PFH Inversion
PFH inversion is based on the assumption that a landslide is a block having a constant mass m. With the constant landslide mass, the landslide momentum (m × v) can be calculated by integrating Figure 6. Three components of the synthetic waveforms and the seismograms of each station. Z denots vertical direction, R denotes radial direction, and T denotes tangential direction. Correlation coefficients are represented by "cc". Light-colored seismograms represent that the lower weighted coefficients were adopted to conduct the PFH inversion.

PFH inversion is based on the assumption that a landslide is a block having a constant mass m.
With the constant landslide mass, the landslide momentum (m × v) can be calculated by integrating the force-time function. As the Taimali landslide first slid toward the northeast, the E-component momentum increased and reached the maximum at the 24th second (Figure 7a). After that, both the E-and N-component momentums decreased, and the E-component momentum reached zero; subsequently, the landslide turned to the west at the 37th second (Figure 7b). The time interval between the 1st and the 37th seconds (09:31:34-09:32:11) can be considered the first stage of the landslide propagation (Figure 7c).  The landslide slid toward the west for only a short time, i.e., 17 s, and turned to the east again. The time interval between the 37th and 54th seconds (09:32:11-09:32:28) can be considered the second stage of the landslide propagation (Figure 7c). After the 54th second (09:32:28), the landslide continued toward the east until it stopped. Therefore, the time interval between the 54th and 174th seconds (09:32:28-09:34:28) can be considered the third stage of the landslide propagation. The The landslide slid toward the west for only a short time, i.e., 17 s, and turned to the east again. The time interval between the 37th and 54th seconds (09:32:11-09:32:28) can be considered the second stage of the landslide propagation (Figure 7c). After the 54th second (09:32:28), the landslide continued toward the east until it stopped. Therefore, the time interval between the 54th and 174th seconds (09:32:28-09:34:28) can be considered the third stage of the landslide propagation. The variations of velocity in the E-and N-components did not decrease to zero at the end of landslide propagation. Yamada et al. [14] noted that long-period noise might accumulate with linear trends as the force is integrated in the time domain. Meanwhile, the slight fluctuations in the momentums in the third stage were possibly caused by the small landslides around the Taimali landslide. During the third stage, some small landslides might have generated weak seismic signals that affected the inversion results. Besides, the trajectories of the landslide movement could be measured from the post-event satellite images and DEM. The starting point was defined at the deepest point of the landslide source area. There were two obvious topographic turning points on the route of the landslide movement. The turning points could be considered the starting points of the second and third stages. In addition, the centroid of the deposit area in the river channel could be considered the termination point of the landslide. Eventually, the runout distances of each stage and average velocities could be calculated (Table 3). Chao et al. [16] performed a waveform inversion using the LFH method for the Taimali landslide. To ensure the equilibrium of forces, the LFH method makes strict assumptions about the force changes during the landslide (unloading/reloading cycles). As a result, the block trajectory deduced by the LFH inversion may not completely correspond to the real terrain, but can only show the overall direction of movement. The PFH method does not confine the force changes during a landslide within certain limits, so a more detailed force-time function can be obtained, which can be used to reflect the correlation between the landslide propagation and the actual terrain. Pan et al. [26] and Huang [27] used a 3D discrete element method to simulate the dynamic processes of the Taimali landslide. The studies both displayed that the landslide moved toward the northeast, which coincides with the results from the PFH inversion method. Besides, Huang [27] reported that the landslide turned to the west at the 40th second, and the movement velocity exceeded 20 m/s. Although the two numerical simulations did not present the force history of the landslide, the simulated movement direction and velocity were both consistent with the results of the PFH inversion method.

The Requirement for the Single-Station Inversion
In general, most waveform inversion methods are limited by the quality of the seismic data and the station coverage. The results of waveform inversion will be more robust if seismic records having high SNR values are used and station coverage is complete. It was particularly noteworthy that PFH inversion could be executed with high-quality seismic records from a single station. Accordingly, single-station waveform inversions were performed on each of the 15 station records. The residual of the M LQ value between the multistation inversion and the single-station inversion was calculated to testify to the quality of the single-station inversion ( Table 4). The M LQ value was defined as the value of the landslide mass multiplied by the runout distance [17]. The results revealed that half of the single-station inversions had residuals below 0.2 and had fine consistency with the multistation inversion. It appeared that the seismic records having higher SNR values led to lower residuals and were more suitable for use in single-station inversion. The M LQ value of the TWKB single-station inversion was overestimated and led to a high residual, which might have been caused by the environmental conditions of the seismometer or geological factors. The seismic amplitudes recorded by the TWKB station were abnormally high and unsuitable for single-station waveform inversion. Although the YHNB station had a high SNR value, the M LQ value of the YHNB single-station inversion was underestimated, possibly due to its distance from the seismic source and to waveform attenuation. The rest of the high residuals of M LQ occurred at the stations having low SNR values. In addition, the single-station inverted force can be applied to synthesize waveforms for other stations to acquire high correlation coefficients, which enhances the reliability of the single-station inversion (Figure 8). to acquire high correlation coefficients, which enhances the reliability of the single-station inversion ( Figure 8).  The residuals between 15-station inversion and single-station inversion reveal that the errors caused by single-station inversion increase with station distance to the landslide (Figure 9). Once the residuals are considered as the errors, the MASB single-station inversion would cause an error of 5%, and result in a slight decrease of 0.1-0.15 in the correlation coefficients between the synthetic and observed waveforms. As the error caused by the single-station inversion is larger than 20% (e.g., YHNB, KMNB, etc.), the correlation coefficient between the synthetic and observed waveforms would decrease significantly to lower than 0.3. The noise level at each station is another factor to influence the waveform inversion. In Figure 9, we can observe that when the number of components with SNR > 8.5 is less than two, the errors of single-station inversion are larger than 20%. Moreover, The residuals between 15-station inversion and single-station inversion reveal that the errors caused by single-station inversion increase with station distance to the landslide (Figure 9). Once the residuals are considered as the errors, the MASB single-station inversion would cause an error of 5%, and result in a slight decrease of 0.1-0.15 in the correlation coefficients between the synthetic and observed waveforms. As the error caused by the single-station inversion is larger than 20% (e.g., YHNB, KMNB, etc.), the correlation coefficient between the synthetic and observed waveforms would decrease significantly to lower than 0.3. The noise level at each station is another factor to influence the waveform inversion. In Figure 9, we can observe that when the number of components with SNR > 8.5 is less than two, the errors of single-station inversion are larger than 20%. Moreover, the azimuthal coverage of the seismic network is not very dense; the azimuth gap is large especially in the east, which will influence the selection of seismic records. Therefore, we deduce that the quality of the seismic data and the station coverage are the main factors in determining the results of waveform inversion. In summary, we suggest that stations having SNR values higher than 8.5 and located within 200 km of a landslide are appropriate for use in single-station waveform inversion and may have residuals of M LQ values below 0.2. The single-station PFH inversion provides a chance to analyze small-scale landslides. Small-scale landslides generate weak seismic signals and may not be detected by many seismic stations, making it difficult to constrain the results of the multistation waveform inversion. The dynamic processes of a landslide can still be extracted by executing single-station PFH inversion with high-SNR seismograms.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 14 of 16 the azimuthal coverage of the seismic network is not very dense; the azimuth gap is large especially in the east, which will influence the selection of seismic records. Therefore, we deduce that the quality of the seismic data and the station coverage are the main factors in determining the results of waveform inversion. In summary, we suggest that stations having SNR values higher than 8.5 and located within 200 km of a landslide are appropriate for use in single-station waveform inversion and may have residuals of MLQ values below 0.2. The single-station PFH inversion provides a chance to analyze small-scale landslides. Small-scale landslides generate weak seismic signals and may not be detected by many seismic stations, making it difficult to constrain the results of the multistation waveform inversion. The dynamic processes of a landslide can still be extracted by executing singlestation PFH inversion with high-SNR seismograms. Figure 9. The correlation between distance to the landslide and residual of single-station inversion.

Conclusions
Waveform inversion techniques provide information for the study of the dynamic processes of a landslide. In this study, the Taimali landslide of 2009 was chosen for the implementation of waveform inversion using the PFH method. According to the pre-and post-event DEMs, the source area of the Taimali landslide had an area of 9.3 × 10 5 m 2 and a volume of 7.5 × 10 7 m 3 . PFH inversion presented the complicated force-time functions and two impact timings during the propagation of the landslide. Thus, the 174-s total of the dynamic processes can be divided into three stages for the different sliding directions: northeast for the first stage, northwest for the second stage, and northeast for the third stage. The average sliding velocities of the three stages were 32.5, 29.4, and 7.5 m/s, respectively. The inversion results indicated that the PFH inversion method could extract the dynamic processes of the landslide, even if the directions of landslide movement varied during propagation. The study also proved that the single-station waveform inversion was feasible by using the PFH inversion method with high-quality seismic data.

Conclusions
Waveform inversion techniques provide information for the study of the dynamic processes of a landslide. In this study, the Taimali landslide of 2009 was chosen for the implementation of waveform inversion using the PFH method. According to the pre-and post-event DEMs, the source area of the Taimali landslide had an area of 9.3 × 10 5 m 2 and a volume of 7.5 × 10 7 m 3 . PFH inversion presented the complicated force-time functions and two impact timings during the propagation of the landslide. Thus, the 174-s total of the dynamic processes can be divided into three stages for the different sliding directions: northeast for the first stage, northwest for the second stage, and northeast for the third stage. The average sliding velocities of the three stages were 32.5, 29.4, and 7.5 m/s, respectively. The inversion results indicated that the PFH inversion method could extract the dynamic processes of the landslide, even if the directions of landslide movement varied during propagation. The study also proved that the single-station waveform inversion was feasible by using the PFH inversion method with high-quality seismic data.