Dynamic Monitoring and Vibration Analysis of Ancient Bridges by Ground-Based Microwave Interferometry and the ESMD Method

In this paper, we propose to conduct a dynamic monitoring and vibration analysis of ancient bridges by means of ground-based microwave interferometry and the extreme-point symmetric mode decomposition (ESMD) method. Ground-based microwave interferometry, a novel non-contact technology with a high accuracy, is used to acquire dynamic time series displacements with environmental excitation factors and a transient load with a car, respectively. The ESMD method, a new alternative to the Hilbert-Huang transform (HHT), is adopted to conduct the instantaneous vibration analysis of Zhaozhou Bridge. Firstly, a series of intrinsic mode functions (IMFs) are obtained together with an optimal adaptive global mean (AGM) curve by using a mode symmetric about the maxima and minima points. Secondly, the instantaneous frequency of each IMF is obtained by the use of a direct interpolation algorithm, which can reconcile the conflict between the period and the frequency for the traditional time-frequency analysis methods. As a representative case, Zhaozhou Bridge, a well-known Chinese ancient bridge constructed more than 1400 years ago, is studied in detail. Four kinds of dynamic time series displacements—two of them acquired by considering only environmental excitation factors for the mid-span and 1/4-span points and the others obtained with the transient load of a car for the mid-span and 1/4-span points—are selected to pursue a comparison of the decomposed IMFs and the instantaneous frequencies to perform the instantaneous vibration analysis of Zhaozhou Bridge. By comparing the results obtained with HHT for the decomposed IMFs and the instantaneous frequencies, the results show that the proposed method has a powerful ability to evaluate the instantaneous dynamic response of ancient bridges.


Introduction
With the inevitable aging of ancient bridges, there is now an urgent problem worldwide for dynamic monitoring and vibration analysis [1].To safeguard the stability of ancient bridges, regular inspection and analysis are essential.At present, dynamic testing under ambient excitation (usually due to traffic or wind) is the main experimental method used to acquire the dynamic responses of bridges in service, which is an approach that has attracted significant research attention [2,3].Most of the typical experimental tests have been carried out using piezoelectric accelerometers and fiber optic sensors [4,5].These transducers can accurately and reliably acquire dynamic time series displacements.However, they need to be fixed in specific positions on the measured bridges and require hardwiring from the transducers to the data acquisition system [6].This is a time-consuming and expensive task and may result in some damages to the ancient bridges.
Ground-based microwave interferometry is an alternative technology for vibration measurement that has been widely applied to determine the dynamic displacement of bridges [7][8][9][10].The great advantages of this technique include non-contact displacement measurement, a wide frequency range of responses, sub-millimetric displacement sensitivity, and quick set-up.By regarding the reflected points of the electromagnetic wave as a series of virtual sensors, the non-contact capability can strongly reduce or even nullify the use of the traditional point sensors, whose positioning involves the use of cumbersome and costly scaffolding, which is occasionally undertaken under hazardous conditions [11].Accuracy comparisons have been made in the laboratory and on-site between ground-based microwave interferometry and the traditional transducers (such as accelerometers and seismometers), and it has been reported that ground-based microwave interferometry can achieve an accuracy of between 1/100 and 1/10 of a millimeter [11,12].
Time-frequency analysis is a commonly used method for structural damage detection by analyzing the vibration-based time series displacements directly.At present, the most widely used time-frequency analysis methods are short-time Fourier transform (STFT), wavelet analysis, and Hilbert-Huang transform (HHT).STFT is the first kind of method of time-frequency analysis, and is a Fourier-related transform used to determine the sinusoidal frequency and phase content of local sections of a signal as it changes over time.Due to the time-bandwidth product theorem, a tradeoff exists between the time and frequency resolution, which makes STFT of limited use [2].Wavelet analysis is a representation of a function by wavelets, which are the scaled and translated copies of a finite-length or fast-decaying oscillating waveform.As wavelets can be localized in both time and frequency, they can be easily applied in the field of structural damage detection [13,14].However, wavelet transform is still an adaptive window Fourier method based on the principle of linear superposition, can only handle non-stationary signals for linear systems, and relies on a priori knowledge [15,16].Furthermore, the results of wavelet analysis are limited by the size of the mother wavelet.
Generally speaking, under the designed load, the dynamic responses of structures should be linear.However, non-stationary and non-linear signals often result from the ambient excitation, and these signals are difficult to process using the traditional Fourier-based methods (e.g., STFT and wavelet analysis).Therefore, a great deal of attention has been placed on analyzing non-stationary and non-linear signals in recent years.In this context, the concept of the instantaneous frequency of a signal, an oscillating change rate during the process of moving back and forth, has become very popular [2,16].Therefore, it is appropriate to analyze the time series displacements from dynamic testing under ambient excitation for bridges.At present, a commonly acceptable definition for instantaneous frequency relies on the analytical signal obtained by the use of HHT [17].HHT is an innovative method of adaptive time-frequency analysis, with the ability to analyze non-linear systems or non-stationary data, which is based on empirical mode decomposition (HHT-EMD) and Hilbert spectral analysis [15].By using the HHT-EMD method, any complicated response signal can be decomposed into a series of finite and small intrinsic mode functions (IMFs), and the meaningful instantaneous frequency of each IMF can be determined and used to build the Hilbert spectrum [17].With the advantage of requiring neither an a priori primary function nor a preset window length, HHT has been successfully applied in the field of vibration-based structural damage detection.Huang et al. [18] noted the frequency downshift as indications of structural yield by examining the Hilbert spectrum and marginal Hilbert spectrum for damage detection of a bridge pier because of piling.Roveri and Carcaterra [2] proposed to adopt a novel HHT-based method for the damage detection of bridge structures under a traveling load, and the damage location is revealed by the inspection of the first instantaneous frequency curve of the corresponding decomposed IMF.Quek et al. [19] illustrated the feasibility of HHT for locating an anomaly, in the form of a crack, delamination, stiffness loss, or boundary in beams and plates, based on physically acquired propagating wave signals.Chen et al. [20] performed a dynamic-based damage detection for large structural systems through the use of HHT, and pointed out that the instantaneous frequency of the decomposed IMFs is more sensitive to structural damage.Kunwar et al. [21] presented the health monitoring of an experimental bridge model based on transient vibration data using HHT.
Extreme-point symmetric mode decomposition (ESMD) is a new alternative to HHT proposed by Wang and Li in 2013 [16].Differing from HHT, the ESMD method permits the residual component to possess a certain number of extreme points, without decomposing a signal to the last trend function with, at the most, one extreme point [16].This residual component can better reflect the evolutionary trend of the whole data.In this context, it can be understood as an optimal adaptive global mean (AGM) curve, which reduces the difficulty of determining the sifting time and improves the accuracy of the decomposed IMFs [16,22].A comparison of the HHT-EMD method and the ESMD method has been performed by decomposing a simulation signal, which consists of three sub-signals with different frequencies and white Gaussian noise, with a signal-to-noise ratio (SNR) of 10 dB.The results showed that the ESMD method performed better than the HHT-EMD method, of which the decomposed IMFs were more consistent with the above three sub-signals, and the mode mixing problems were also evidently decreased [17].Moreover, HHT has the disadvantage of computing a meaningful instantaneous frequency, so this method requires a hypothesis that the signal itself is locally smooth and the derivative of the phase function of the input data exists [16,23].Unfortunately, the dynamic time series displacements of bridges do not always conform to the hypothesis.In order to address this issue, the ESMD method has been adopted to determine the instantaneous frequency and amplitude in a direct interpolation approach, which utilizes the data-based direct interpolation approach to calculate the meaningful instantaneous frequency point by point in time [16,24].Liu et al. [23] conducted a stability analysis of the progressive collapse of a building model using the ESMD method, and the results show that the obtained instantaneous frequency of each IMF can not only detect damages of a building model, but can also analyze load-bearing conditions.
The purpose of this study was to perform the instantaneous vibration analysis of ancient bridges through the use of the ESMD method based on the dynamic time series displacements.In this study, we take the ancient bridge of Zhaozhou Bridge as an experimental object, of which dynamic time series displacements were obtained by ground-based microwave interferometry (the IBIS-S instrument).Two kinds of tests were performed to acquire the dynamic time series displacements for Zhaozhou Bridge; one test was undertaken with only environmental excitation factors and the other test was undertaken with the transient load of a car weighing more than two tons.The ESMD method was adopted to yield a series of IMFs to reflect the overall tendency of the projected displacement according to the magnitude of the frequency.The direct interpolation algorithm was adopted to obtain the meaningful instantaneous frequencies of the corresponding decomposed IMFs to analyze the stability of Zhaozhou Bridge under the transient load of a car.The rest of this paper is organized as follows.Zhaozhou Bridge and the measurement method are introduced in Section 2. Section 3 introduces the ESMD method for non-linear, non-stationary data processing, including the decomposition algorithm of the ESMD method for the time series displacements and the direct interpolation algorithm for instantaneous frequency.Section 4 describes the experimental results and analysis for Zhaozhou Bridge, followed by the conclusion in Section 5.

Zhaozhou Bridge
Zhaozhou Bridge, which is also known as Anji Bridge or Dashi Bridge, is the oldest open-spandrel stone segmental-arch bridge in the world.Credited to the design of a craftsman named Li Chun, the bridge was constructed in the Sui Dynasty, more than 1400 years ago, and is located on the Xiaohe River in Zhaozhou County, Hebei province, China.Figure 1

Ground-Based Microwave Interferometry
In this study, Imaging by Interferometric Survey (IBIS-S), a typical system based on ground-based microwave interferometry for the remote dynamic monitoring of structures, was adopted to acquire accurate dynamic time series displacements.The IBIS-S system was developed by the Ingegneria Dei Sistemi (IDS) Company, Pisa, Italy, in collaboration with the Department of Electronics and Telecommunications of Florence University.The IBIS-S instrument consists of a radar unit, a control PC, a power supply unit, and a tripod, with the ability of displacement monitoring in any weather conditions, independent of daylight [10].The radar unit, a coherent sensor module, can generate, transmit, and receive the electromagnetic signals, and further process the preserved phase information of the received signals to calculate the displacement of the monitored objects.The control PC is connected to the radar unit by a standard USB 2.0 interface (Figure 2), and it can be used to configure the parameters and manage and store the measurements.It also allows the user to view the initial results of the key locations in real time [10].
The stepped-frequency continuous wave (SF-CW) and ground-based microwave interferometry are the two key technologies of the IBIS-S system.
The purpose of the SF-CW technique is to provide the capability of long-distance transmission and a higher range resolution, without the need to install multiple units on the monitored objects [28].By using the SF-CW technique, the radar unit of the IBIS-S system transmits a set of continuous electromagnetic waves at discrete frequency values, with a sampling bandwidth B and a constant interval f  .Therefore, the range resolution R  can be calculated using the corresponding continuous-wave signals by , where c denotes the speed of light in a vacuum.For the IBIS-S system, the bandwidth range B is 300 MHz, and thus a range resolution of up to 0.5 m can be obtained [10].This means that each object point can be observed separately by not less than 0.5 m from another point along the radial direction [29].Moreover, since the significant frequency content of time series displacements is in the range of 0-20 Hz for a civil engineering structure, the system is highly suited to analyzing the instantaneous dynamic response of Zhaozhou Bridge, with a sampling rate of up to 200 Hz to sample the scenario and a bandwidth of 300 MHz.In other words,

Ground-Based Microwave Interferometry
In this study, Imaging by Interferometric Survey (IBIS-S), a typical system based on ground-based microwave interferometry for the remote dynamic monitoring of structures, was adopted to acquire accurate dynamic time series displacements.The IBIS-S system was developed by the Ingegneria Dei Sistemi (IDS) Company, Pisa, Italy, in collaboration with the Department of Electronics and Telecommunications of Florence University.The IBIS-S instrument consists of a radar unit, a control PC, a power supply unit, and a tripod, with the ability of displacement monitoring in any weather conditions, independent of daylight [10].The radar unit, a coherent sensor module, can generate, transmit, and receive the electromagnetic signals, and further process the preserved phase information of the received signals to calculate the displacement of the monitored objects.The control PC is connected to the radar unit by a standard USB 2.0 interface (Figure 2), and it can be used to configure the parameters and manage and store the measurements.It also allows the user to view the initial results of the key locations in real time [10].
The stepped-frequency continuous wave (SF-CW) and ground-based microwave interferometry are the two key technologies of the IBIS-S system.
The purpose of the SF-CW technique is to provide the capability of long-distance transmission and a higher range resolution, without the need to install multiple units on the monitored objects [28].By using the SF-CW technique, the radar unit of the IBIS-S system transmits a set of continuous electromagnetic waves at discrete frequency values, with a sampling bandwidth B and a constant interval ∆ f .Therefore, the range resolution ∆R can be calculated using the corresponding continuous-wave signals by ∆R = c 2•B , where c denotes the speed of light in a vacuum.For the IBIS-S system, the bandwidth range B is 300 MHz, and thus a range resolution of up to 0.5 m can be obtained [10].This means that each object point can be observed separately by not less than 0.5 m from another point along the radial direction [29].Moreover, since the significant frequency content of time series displacements is in the range of 0-20 Hz for a civil engineering structure, the system is highly suited to analyzing the instantaneous dynamic response of Zhaozhou Bridge, with a sampling rate of up to 200 Hz to sample the scenario and a bandwidth of 300 MHz.In other words, the sampling interval of 0.005 s is well suited to providing a good wave definition for the acquired time series displacements [6].Ground-based microwave interferometry is adopted to ensure a higher displacement accuracy at a high sampling rate [30].In general, amplitude and phase information can be acquired for each measurement using ground-based microwave interferometry.Therefore, considering a single target, a tiny displacement d along the radar line-of-sight direction can be calculated through comparing the phase shift at different times by , where  is the length of the electromagnetic wave and   is the phase shift [31].For IBIS-S, the Ku-band of the working Radar is between 16.6-16.9 GHz, the approximate value of f c   is 18.07 mm, the maximum value of   is  , and therefore, the maximum range of displacement d along the radar line-of-sight direction is about 4.5175 mm for the adjacent sampling interval.However, since ground-based microwave interferometry provides a measurement of the line-of-sight displacement, it requires prior knowledge of the direction of motion to evaluate the actual displacement.For the monitored Zhaozhou Bridge in this study, the displacement under environmental excitation factors and traffic loads can be regarded as the projected displacement, which can be calculated by making straightforward geometric projections [10,28].
In summary, in typical measurement conditions, the radar unit has the following characteristics: the maximum detection distance is up to 1 km, the range resolution is up to 0.50 m, the sampling rate is up to 200 Hz, and the displacement measurement accuracy is up to 0.01 mm.However, in the practical application, the displacement measurement accuracy depends on the power of the reflected target.The general measurement accuracy is between 0.01 mm and 0.1 mm.

1/4-span point
Mid-span point Ground-based microwave interferometry is adopted to ensure a higher displacement accuracy at a high sampling rate [30].In general, amplitude and phase information can be acquired for each measurement using ground-based microwave interferometry.Therefore, considering a single target, a tiny displacement d along the radar line-of-sight direction can be calculated through comparing the phase shift at different times by d = λ∆ϕ 4π , where λ is the length of the electromagnetic wave and ∆ϕ is the phase shift [31].For IBIS-S, the Ku-band of the working Radar is between 16.6-16.9GHz, the approximate value of λ = c/ f is 18.07 mm, the maximum value of ∆ϕ is π, and therefore, the maximum range of displacement d along the radar line-of-sight direction is about 4.5175 mm for the adjacent sampling interval.However, since ground-based microwave interferometry provides a measurement of the line-of-sight displacement, it requires prior knowledge of the direction of motion to evaluate the actual displacement.For the monitored Zhaozhou Bridge in this study, the displacement under environmental excitation factors and traffic loads can be regarded as the projected displacement, which can be calculated by making straightforward geometric projections [10,28].
In summary, in typical measurement conditions, the radar unit has the following characteristics: the maximum detection distance is up to 1 km, the range resolution is up to 0.50 m, the sampling rate is up to 200 Hz, and the displacement measurement accuracy is up to 0.01 mm.However, in the practical application, the displacement measurement accuracy depends on the power of the reflected target.The general measurement accuracy is between 0.01 mm and 0.1 mm.

Dynamic Time Series Displacements Acquisition
Two kinds of tests were carried out to acquire the dynamic time series displacements for Zhaozhou Bridge with the IBIS-S instrument, for which the parameter settings are shown in Table 1.One test was undertaken with only environmental excitation factors (Test 1).In other words, no pedestrians or vehicles passed over the bridge during the data acquisition.The other test was undertaken with the transient load of a car weighing more than 2 tons (Test 2), which went back and forth over the bridge four times with a speed of about 20 km/h.For test 2, aiming to avoid any damage to Zhaozhou Bridge, only one set of time series displacements was collected to analyze the instantaneous dynamic response for Zhaozhou Bridge, of which the duration was 1 min and 47 s.For Test 1, five sets of time series displacements were collected, which have similar change tendencies.Moreover, in order to facilitate the subsequent contrast analysis, the experimental duration was between 2 and 3 min for each set.Hence, only one of the five sets of time series displacements was selected to conduct the analysis in this paper, of which the duration was 2 min and 45 s, as shown in Table 1.Obviously, the two tests employed short observation periods, so it was difficult to perform the global stability analysis for Zhaozhou Bridge.However, as the tests were more sensitive to reflecting the instantaneous stability of bridges for the instantaneous frequency of the decomposed IMFs and a high sampling rate of 199.17 Hz for the IBIS-S instrument [20], the duration of the two tests was enough to evaluate the instantaneous dynamic response of Zhaozhou Bridge, especially for Test 2 with four continuous excitations under the transient load of a car.Furthermore, in order to guarantee the uniform initial reference displacement for the two kinds of tests, the IBIS-S instrument was fixed during the two tests.As shown in Figures 1 and 2, the IBIS-S instrument was located at one side of Zhaozhou Bridge.The angle of altitude of the radar unit was set to 12 • , so that the two antennas on the radar unit could be aligned to the mid-span point of the central span.Considering that Zhaozhou Bridge is of great historical heritage, no corner reflectors were installed on the lower surface of the central span.The acquired time series displacements along the radar line-of-sight direction were directly calculated by the reflected microwave signal from the natural surface of the arch, which was then vertically projected to obtain the projected time series displacements [10,28].

Methods
The ESMD method involves two distinct steps: (1) the decomposition step to yield a series of IMFs together with an optimal adaptive global mean (AGM) curve; and (2) the direct interpolation algorithm to yield the instantaneous frequency of each IMF.

Decomposition Algorithm of the ESMD Method for Time Series Displacements
Differing from the scheme of the HHT-EMD method, i.e., making a mode symmetric about its upper and lower envelope interpolated by the local maxima and minima points, respectively, the scheme of the ESMD method is to make a mode symmetric about its own maxima and minima points [16].The sifting process is implemented by the aid of one, two, three, or more inner curves by the midpoints of the line segments connecting the local maxima and minima points [16,23].The purpose of the decomposition of the ESMD method is to yield a series of IMFs together with an optimal AGM curve.The whole decomposition process is shown in Figure 3. purpose of the decomposition of the ESMD method is to yield a series of IMFs together with an optimal AGM curve.The whole decomposition process is shown in Figure 3.

Time series data Y
Step 1: Find all the local extreme points (maxima points plus minima points) of the input data Y and numerate them by E i (1≤i≤n) Step 2: Connect all the adjacent points with line segments and mark their midpoints by F i (1≤i≤n-1) Step 3: Add a left and a right boundary midpoints F 0 and F n through a linear interpolation Step 4: Construct 2 curves interpolated by the odd and even midpoints with all these n + 1 midpoints and calculate their mean value by Step 5: We denote the original input data as and the optimal AGM curve as , where i y is the obtained time series displacements of bridges by ground-based microwave interferometry and i r is the automatically selected time series displacements by the least squares method [16].The variance 0  of the input data relative to the total mean The variance  of the input data relative to the optimal AGM curve is defined as (2):

Direct Interpolation Algorithm for Instantaneous Frequency
There is generally no local meaning for a frequency at a given point, but there is a frequency-modulating phenomenon.Therefore, the definition of instantaneous frequency is a We denote the original input data as Y = {y i } N i=1 and the optimal AGM curve as R = {r i } N i=1 , where y i is the obtained time series displacements of bridges by ground-based microwave interferometry and r i is the automatically selected time series displacements by the least squares method [16].
The variance σ 0 of the input data relative to the total mean Y = N ∑ i=1 y i /N is defined as (1): The variance σ of the input data relative to the optimal AGM curve is defined as (2):

Direct Interpolation Algorithm for Instantaneous Frequency
There is generally no local meaning for a frequency at a given point, but there is a frequency-modulating phenomenon.Therefore, the definition of instantaneous frequency is a controversial issue [23,32].Only when a quantity varies in a periodic oscillating manner can instantaneous frequency be regarded as an oscillating change rate during the process of moving back and forth [16].Therefore, it is appropriate to analyze the time series displacements from dynamic testing under ambient excitation for ancient bridges.Nowadays, the popular method is HHT, which is superior to the Fourier transform, wavelet transform, and other analytical forms [33][34][35].However, HHT has the disadvantage of obtaining a meaningful instantaneous frequency which represents an actual physical phenomenon, and it requires a hypothesis that the derivative of the phase function of the input data exists [23].
In fact, it is actually the uniform running-mean processing for each integral transform of HHT, and the Hilbert spectral analysis is substituted by the data-based one for instantaneous frequency.Moreover, as an instantaneous frequency, it should be capable of reflecting the intermittent case rather than excluding the adjacent equal situation.In this context, the period is defined relative to a segment of time and the frequency needs to be understood point by point.There is a conflict between the period and the frequency during the uniform running-mean processing [16].Therefore, the direct interpolation algorithm is proposed to reconcile the conflict.
We denote the discrete form of each IMF as (t k , y k ), k = 1, 2, • • • N, and the detailed direct interpolation algorithm can be shown as follows.
Step 1: Traverse to find all of the quasi-extreme points of each IMF which satisfies (3), and enumerate them as set E = (t e i , y e i ), Step 2: Define the frequency interpolation coordinates (a i , f i ) by using set E.
1: for 1 to m 2: if y e i == y e i+1 then 3: if i == 1 then 4: if (t e i , y e i ) and t e i+1 , y e i+1 are extreme points then 11: else 13: end for Step 3: Add the boundary points with a linear interpolating method.
for the left boundary point: if Step 4: Obtain a curve f (t) by using cubic spline interpolation with all the discrete points (a i , f i ).
To be meaningful, the instantaneous frequency curve is defined by (4) [16]:

Procedure for the Instantaneous Vibration Analysis
The procedure for the instantaneous vibration analysis of Zhaozhou Bridge based on the ESMD method is designed as follows.
(1) Apply the ESMD method to the dynamic time series displacements (projected displacement).
Determine the optimal variance ratio with the corresponding sifting time to obtain the optimal AGM curve for each set of time series displacements, and yield a series of IMFs.(2) Compare the displacements and sudden variations of the main IMFs between the 1/4 span point and the mid-span point under the environmental excitation factors, and evaluate the instantaneous dynamic response of the 1/4 span point and the mid-span point under the environmental excitation factors.(3) Apply the direct interpolation algorithm to each IMF of the time series displacements to obtain the corresponding instantaneous frequency and amplitude, and perform the instantaneous vibration analysis of Zhaozhou Bridge by analyzing the instantaneous frequencies and amplitudes of the main IMFs related to Zhaozhou Bridge.

Simulation Experiment
In this study, aiming to demonstrate the validity of the ESMD method compared to the HHT-EMD method [15], a nonlinear and non-stationary signal s(t) = s 1 (t) + s 2 (t) + s 3 (t) + w(t) is simulated by adding white Gaussian noise to three useful sub-signals with an SNR of 20 dB.Where s 1 (t) = cos(20πt + 2 sin(2t)), s 2 (t) = sin(4πt), s 3 (t) = cos(2πt), and w(t) is white Gaussian noise.The waveforms of the simulation signal s(t) and its component signals s 1 (t), s 2 (t), s 3 (t) and w(t) are shown in Figure 4.
Through the use of the ESMD method and the HHT-EMD method, respectively, the simulation signal is decomposed into a series of IMFs, which are the red curves shown in Figure 5.Moreover, the blue curves of the component signals are exhibited along with their corresponding IMFs.As shown in Figure 5a, through the use of the HHT-EMD method, the simulation signal   Through the use of the ESMD method and the HHT-EMD method, respectively, the simulation signal s(t) is decomposed into a series of IMFs, which are the red curves shown in Figure 5.Moreover, the blue curves of the component signals are exhibited along with their corresponding IMFs.As shown in Figure 5a, through the use of the HHT-EMD method, the simulation signal s(t) is decomposed into five IMFs.The decomposed IMF2 is close to s 1 (t) even though some local errors exist due to mode mixing problems.However, in the posterior halves of IMF3 and IMF4, they are very different from component signals s 2 (t) and s 3 (t) on account of some overshoot and undershoot problems.As shown in Figure 5b, through the use of the ESMD method, the simulation signal s(t) is decomposed into four IMFs, which are equal to the number of component signals.Furthermore, the mode mixing problems and overshoot and undershoot problems almost disappear, and IMF2 to IMF4 accord with component signals s 1 (t), s 2 (t), and s 3 (t), indicating that the ESMD method performs better than the HHT-EMD method in decomposing a signal into a series of physically meaningful representations.Through the use of the ESMD method and the HHT-EMD method, respectively, the simulation signal is decomposed into a series of IMFs, which are the red curves shown in Figure 5.Moreover, the blue curves of the component signals are exhibited along with their corresponding IMFs.As shown in Figure 5a, through the use of the HHT-EMD method, the simulation signal

Results of Time Series Displacements from Ground-Based Microwave Interferometry
Through using IBIS-S in this study, the whole scene of the lower surface of Zhaozhou Bridge was measured inside the radar beam.Therefore, according to the range resolution of 0.5 m of IBIS-S and the central span of 37.47 m of Zhaozhou Bridge, the dynamic responses of about 30 points can be obtained.However, in the field of civil engineering, the mid-span point and the 1/4-span point are the two key locations used to reflect the dynamic response of a bridge.Therefore, in this study, four kinds of dynamic time series displacements were acquired to perform the instantaneous vibration analysis of Zhaozhou Bridge at the location of the mid-span point and the 1/4-span point, as shown in Figure 2. Two of them were acquired only with environmental excitation factors for the mid-span point with a line-of-sight distance of 23 m to the radar unit (named Rbin46_N) and the 1/4-span point with a line-of-sight distance of 32 m to the radar unit (named Rbin64_N).The other two kinds of dynamic time series displacements were obtained with the transient load of a car for the mid-span point (named Rbin46_L) and the 1/4-span point (named Rbin64_L).In order to ensure the reflected intensity of the electromagnetic wave from the natural surface of Zhaozhou Bridge, the thermal SNR-a ratio between the power of the received signal of a pixel and the thermal noise power of the sensor-was calculated, which can determine the measurement accuracy of the IBIS-S instrument [36].Generally, when the thermal SNR is better than 30 dB, the standard deviation of measured displacement can reach 0.09 mm for the IBIS-S instrument [37].Therefore, the empirical thermal SNR of 35 dB was set as the critical value to ensure that the measurement accuracy was better than 0.05 mm for the IBIS-S instrument.As shown in Figure 6, the thermal SNR was 53.8 dB for point Rbin46 and 52.5 dB for point Rbin64, which could ensure that accurate displacements could be obtained by ground-based microwave interferometry.Moreover, another interesting peak occurred at about 45 m (point Rbin90, as shown in Figure 6) along the line-of-sight direction to the radar unit, whose thermal SNR was 53.3 dB.However, as the central span is 37.47 m for Zhaozhou Bridge, point Rbin90 was beyond the scope of the monitored Zhaozhou Bridge.As the area within a red rectangle shown in Figure 2, there is a vertical retaining wall with a smooth surface at the end side of Zhaozhou Bridge in front of the IBIS-S instrument, which can be regarded as a strong backscatter.Furthermore, it has a distance of about 45 m from the vertical retaining wall to the IBIS-S instrument.Therefore, the interesting peak of point Rbin90 may be caused by the vertical retaining wall.

Results of Time Series Displacements from Ground-Based Microwave Interferometry
Through using IBIS-S in this study, the whole scene of the lower surface of Zhaozhou Bridge was measured inside the radar beam.Therefore, according to the range resolution of 0.5 m of IBIS-S and the central span of 37.47 m of Zhaozhou Bridge, the dynamic responses of about 30 points can be obtained.However, in the field of civil engineering, the mid-span point and the 1/4-span point are the two key locations used to reflect the dynamic response of a bridge.Therefore, in this study, four kinds of dynamic time series displacements were acquired to perform the instantaneous vibration analysis of Zhaozhou Bridge at the location of the mid-span point and the 1/4-span point, as shown in Figure 2. Two of them were acquired only with environmental excitation factors for the mid-span point with a line-of-sight distance of 23 m to the radar unit (named Rbin46_N) and the 1/4-span point with a line-of-sight distance of 32 m to the radar unit (named Rbin64_N).The other two kinds of dynamic time series displacements were obtained with the transient load of a car for the mid-span point (named Rbin46_L) and the 1/4-span point (named Rbin64_L).In order to ensure the reflected intensity of the electromagnetic wave from the natural surface of Zhaozhou Bridge, the thermal SNR-a ratio between the power of the received signal of a pixel and the thermal noise power of the sensor-was calculated, which can determine the measurement accuracy of the IBIS-S instrument [36].Generally, when the thermal SNR is better than 30 dB, the standard deviation of measured displacement can reach 0.09 mm for the IBIS-S instrument [37].Therefore, the empirical thermal SNR of 35 dB was set as the critical value to ensure that the measurement accuracy was better than 0.05 mm for the IBIS-S instrument.As shown in Figure 6, the thermal SNR was 53.8 dB for point Rbin46 and 52.5 dB for point Rbin64, which could ensure that accurate displacements could be obtained by ground-based microwave interferometry.Moreover, another interesting peak occurred at about 45 m (point Rbin90, as shown in Figure 6) along the line-of-sight direction to the radar unit, whose thermal SNR was 53.3 dB.However, as the central span is 37.47 m for Zhaozhou Bridge, point Rbin90 was beyond the scope of the monitored Zhaozhou Bridge.As the area within a red rectangle shown in Figure 2, there is a vertical retaining wall with a smooth surface at the end side of Zhaozhou Bridge in front of the IBIS-S instrument, which can be regarded as a strong backscatter.Furthermore, it has a distance of about 45 m from the vertical retaining wall to the IBIS-S instrument.Therefore, the interesting peak of point Rbin90 may be caused by the vertical retaining wall.In this study, in order to facilitate the contrast analysis between the time series displacements acquired only with environmental excitation factors and with the transient load of a car, the length of the time series displacements acquired only with environmental excitation factors was shortened to about 107 s, which was the same as the length of the time series displacements acquired with the transient load of a car.For each set of time series displacements, the vertical component of the corresponding projected displacement-a linear displacement perpendicular to the horizontal plane-was calculated as the input time series displacements for further analysis [10].The curves of the projected displacements of points Rbin46_N and Rbin64_N are shown in Figure 7a,c.As a In this study, in order to facilitate the contrast analysis between the time series displacements acquired only with environmental excitation factors and with the transient load of a car, the length of the time series displacements acquired only with environmental excitation factors was shortened to about 107 s, which was the same as the length of the time series displacements acquired with the transient load of a car.For each set of time series displacements, the vertical component of the corresponding projected displacement-a linear displacement perpendicular to the horizontal plane-was calculated as the input time series displacements for further analysis [10].The curves of the projected displacements of points Rbin46_N and Rbin64_N are shown in Figure 7a,c.As a relative displacement of the acquired projected displacement, the initial reference displacement was defined as 0 mm.The inspection of the projected displacement curves clearly highlights that: (1) the trend of the projected displacement curve of point Rbin46_N (Figure 7a) is similar to that of point Rbin64_N (Figure 7c); and (2) most of the displacement differences are less than 0.3 mm between the adjacent wave peak and wave bottom for the two points, which can be considered as an acceptable displacement difference for Zhaozhou Bridge with the environmental excitation factors.
The curves of the projected displacements of points Rbin46_L and Rbin64_L are shown in Figure 7b,d.The inspection of the projected displacement curves clearly highlights that: (1) compared with the projected displacement curve of points Rbin46_N and Rbin64_N acquired only with the environmental excitation factors, the displacement differences are greatly enlarged for the acquired projected displacements with the transient load of a car, especially for point Rbin64_L with a maximum displacement difference of up to 1 mm between the adjacent wave peak and wave bottom at the time of 35 s (Figure 7d); and (2) compared with the projected displacement curve of point Rbin46_L, the projected displacement curves are more complex and the displacement differences are much bigger than those of point Rbin64_L.These results indicate that the main stress point is located at the 1/4 span point of Zhaozhou Bridge.relative displacement of the acquired projected displacement, the initial reference displacement was defined as 0 mm.The inspection of the projected displacement curves clearly highlights that: (1) the trend of the projected displacement curve of point Rbin46_N (Figure 7a) is similar to that of point Rbin64_N (Figure 7c); and (2) most of the displacement differences are less than 0.3 mm between the adjacent wave peak and wave bottom for the two points, which can be considered as an acceptable displacement difference for Zhaozhou Bridge with the environmental excitation factors.
The curves of the projected displacements of points Rbin46_L and Rbin64_L are shown in Figure 7b,d.The inspection of the projected displacement curves clearly highlights that: (1) compared with the projected displacement curve of points Rbin46_N and Rbin64_N acquired only with the environmental excitation factors, the displacement differences are greatly enlarged for the acquired projected displacements with the transient load of a car, especially for point Rbin64_L with a maximum displacement difference of up to 1 mm between the adjacent wave peak and wave bottom at the time of 35 s (Figure 7d); and (2) compared with the projected displacement curve of point Rbin46_L, the projected displacement curves are more complex and the displacement differences are much bigger than those of point Rbin64_L.These results indicate that the main stress point is located at the 1/4 span point of Zhaozhou Bridge.

Results of IMF Decomposition and Analysis
In order to acquire the optimal AGM curves, and to yield a series of IMFs, the optimal variance ratio and the corresponding sifting time needed to be calculated according to (2).In this study, considering the complexity of the time series displacements, the number of extreme points of the last residual was set as 12 (step 6 in Figure 3) and the maximum number of iterations was set as 30 (step 7 in Figure 3).Therefore, the resulting variance ratio figures were acquired as shown in Figure 8.The minimum variance ratio and the corresponding sifting times of the four kinds of time series displacements are shown in Table 2.

Results of IMF Decomposition and Analysis
In order to acquire the optimal AGM curves, and to yield a series of IMFs, the optimal variance ratio and the corresponding sifting time needed to be calculated according to (2).In this study, considering the complexity of the time series displacements, the number of extreme points of the last residual was set as 12 (step 6 in Figure 3) and the maximum number of iterations was set as 30 (step 7 in Figure 3).Therefore, the resulting variance ratio figures were acquired as shown in Figure 8.The minimum variance ratio and the corresponding sifting times of the four kinds of time series displacements are shown in Table 2.According to the acquired minimum variance ratio and the sifting time, a series of IMFs were yielded, together with the optimal AGM curves, for the four kinds of time series displacements.The four optimal AGM curves for the projected displacement curves are shown in Figure 7.
As shown in Figure 9, each set of time series displacements was decomposed into eight IMF components (Mode 1 to Mode 8) via the ESMD method.The inspection of the IMFs in this figure clearly highlights that: (1) For the IMFs of point Rbin46_N (Figure 9a) and point Rbin64_N (Figure 9c), the projected displacement variations gradually decrease from Mode 1 to Mode 8. Unlike the IMFs of point Rbin46_N and point Rbin64_N, the variation of the projected displacements gradually decreases from Mode 1 to Mode 6 for point Rbin46_L (Figure 9b) and point Rbin64_L (Figure 9d), and the projected displacement variation of Mode 7 and Mode 8 is in the interval of −1 mm to 1 mm, which is the same as the corresponding original time series displacements.This reflects the overall vibration trend of Zhaozhou Bridge with the transient load of a car.(2) Many obvious sudden variations appear on the projected displacement curves of Mode 1 and Mode 2 of point Rbin46_N (Figure 9a), but only some small sudden variations occur on the projected displacement curves of Mode 1 and Mode 2 of point Rbin64_N (Figure 9c).This phenomenon indicates that the mid-span point is much more likely to be affected by the environmental excitation factors than the 1/4 span point of Zhaozhou Bridge.According to the acquired minimum variance ratio and the sifting time, a series of IMFs were yielded, together with the optimal AGM curves, for the four kinds of time series displacements.The four optimal AGM curves for the projected displacement curves are shown in Figure 7.
As shown in Figure 9, each set of time series displacements was decomposed into eight IMF components (Mode 1 to Mode 8) via the ESMD method.The inspection of the IMFs in this figure clearly highlights that: (1) For the IMFs of point Rbin46_N (Figure 9a) and point Rbin64_N (Figure 9c), the projected displacement variations gradually decrease from Mode 1 to Mode 8. Unlike the IMFs of point Rbin46_N and point Rbin64_N, the variation of the projected displacements gradually decreases from Mode 1 to Mode 6 for point Rbin46_L (Figure 9b) and point Rbin64_L (Figure 9d), and the projected displacement variation of Mode 7 and Mode 8 is in the interval of −1 mm to 1 mm, which is the same as the corresponding original time series displacements.This reflects the overall vibration trend of Zhaozhou Bridge with the transient load of a car.(2) Many obvious sudden variations appear on the projected displacement curves of Mode 1 and Mode 2 of point Rbin46_N (Figure 9a), but only some small sudden variations occur on the projected displacement curves of Mode 1 and Mode 2 of point Rbin64_N (Figure 9c).This phenomenon indicates that the mid-span point is much more likely to be affected by the environmental excitation factors than the 1/4 span point of Zhaozhou Bridge.Figure 10 shows the decomposed IMF components of point Rbin64_L obtained via the HHT-EMD method [18].Compared with the decomposed IMFs of point Rbin64_L obtained via the ESMD method, as shown in Figure 9d, the displacement curves of the main IMFs of Mode 1 to Mode 4 have the same trend, which validates the correctness of the decomposed IMFs obtained via the ESMD method.Nevertheless, using the HHT-EMD method, 12 IMFs are decomposed, which is more Figure 10 shows the decomposed IMF components of point Rbin64_L obtained via the HHT-EMD method [18].Compared with the decomposed IMFs of point Rbin64_L obtained via the ESMD method, as shown in Figure 9d, the displacement curves of the main IMFs of Mode 1 to Mode 4 have the same trend, which validates the correctness of the decomposed IMFs obtained via the ESMD method.
Nevertheless, using the HHT-EMD method, 12 IMFs are decomposed, which is more than the eight IMFs obtained via the ESMD method, and the displacement curve trends of Mode 5 to Mode 12 are coarse compared with those of Mode 5 to Mode 8 decomposed by the ESMD method.

Results of Time-Frequency Analysis
Using the proposed direct interpolation algorithm, the instantaneous frequency and amplitude of each IMF for each set of time series displacements could be calculated.Figure 11 shows the curves of acquired instantaneous frequency of Mode 1 to Mode 8 (FM1 to FM8) with different colors for

Results of Time-Frequency Analysis
Using the proposed direct interpolation algorithm, the instantaneous frequency and amplitude of each IMF for each set of time series displacements could be calculated.Figure 11 shows the curves of acquired instantaneous frequency of Mode 1 to Mode 8 (FM1 to FM8) with different colors for each set of time series displacements.The instantaneous frequency of each IMF reduces gradually from Mode 1 to Mode 8 for each set of time series displacements.Thereinto, the red curve indicates the instantaneous frequency of Mode 1 for each set of time series displacements.According to the decomposed IMFs, as shown in Figure 9, the main projected displacement variations focus on Mode 1 and Mode 2, and especially Mode 1. Furthermore, compared with the projected displacement curves of Mode 1 for point Rbin46_N and point Rbin64_N, more sudden variations appear on the projected displacement curves of Mode 1 for point Rbin46_L and point Rbin64_L.These results indicate that the projected displacement curve of Mode 1 is the main projected displacement variation related to the bridge, which is the same as HHT [2,21].Moreover, in order to estimate the efficiency of the instantaneous frequency of Mode 1 from noisy data to evaluate the instantaneous dynamic response of Zhaozhou Bridge, the instantaneous frequencies F1 to F4 with a mean (MN) greater than 0.1 Hz, corresponding to the decomposed Mode 1 to Mode 4 for each set of time series displacement, were selected to make a comparison.The standard deviation (SD) of data must always be understood in the context of the mean of the data, and there is a big difference in the mean of the instantaneous frequencies F1 to F4 for each set of time series displacements, as shown in Table 3.Therefore, in this study, the coefficient of variation (CV) was adopted to evaluate the dispersion of the frequency distribution, which is defined as the ratio of the SD to the MN.As shown in Table 3, the CVs of the instantaneous frequency F1 for each set of time series displacements are much less than those of the corresponding instantaneous frequencies F2 to F4.This further indicates that the projected displacement curve of Mode 1 is the main projected displacement variation related to the bridge.In addition, for the time series displacements with the transient load of a car (points Rbin46_Y and Rbin64_Y), the CVs of instantaneous frequency F1 are somewhat larger than the acquired CVs of instantaneous frequency F1 for the time series displacements with environmental excitation factors (points Rbin46_N and Rbin64_N).The reason for this is that the dispersion of instantaneous frequency distribution can be reduced by the external force of the transient load of a car related to the main projected displacement variation of Zhaozhou Bridge.each set of time series displacements.The instantaneous frequency of each IMF reduces gradually from Mode 1 to Mode 8 for each set of time series displacements.Thereinto, the red curve indicates the instantaneous frequency of Mode 1 for each set of time series displacements.According to the decomposed IMFs, as shown in Figure 9, the main projected displacement variations focus on Mode 1 and Mode 2, and especially Mode 1. Furthermore, compared with the projected displacement curves of Mode 1 for point Rbin46_N and point Rbin64_N, more sudden variations appear on the projected displacement curves of Mode 1 for point Rbin46_L and point Rbin64_L.These results indicate that the projected displacement curve of Mode 1 is the main projected displacement variation related to the bridge, which is the same as HHT [2,21].Moreover, in order to estimate the efficiency of the instantaneous frequency of Mode 1 from noisy data to evaluate the instantaneous dynamic response of Zhaozhou Bridge, the instantaneous frequencies F1 to F4 with a mean (MN) greater than 0.1 Hz, corresponding to the decomposed Mode 1 to Mode 4 for each set of time series displacement, were selected to make a comparison.The standard deviation (SD) of data must always be understood in the context of the mean of the data, and there is a big difference in the mean of the instantaneous frequencies F1 to F4 for each set of time series displacements, as shown in Table 3.Therefore, in this study, the coefficient of variation (CV) was adopted to evaluate the dispersion of the frequency distribution, which is defined as the ratio of the SD to the MN.As shown in Table 3, the CVs of the instantaneous frequency F1 for each set of time series displacements are much less than those of the corresponding instantaneous frequencies F2 to F4.This further indicates that the projected displacement curve of Mode 1 is the main projected displacement variation related to the bridge.In addition, for the time series displacements with the transient load of a car (points Rbin46_Y and Rbin64_Y), the CVs of instantaneous frequency F1 are somewhat larger than the acquired CVs of instantaneous frequency F1 for the time series displacements with environmental excitation factors (points Rbin46_N and Rbin64_N).The reason for this is that the dispersion of instantaneous frequency distribution can be reduced by the external force of the transient load of a car related to the main projected displacement variation of Zhaozhou Bridge.Figure 12 shows the instantaneous frequency and amplitude of Mode 1 for each set of time series displacements.The inspection of the instantaneous frequency and amplitude shown in this figure clearly highlights that: (1) For the instantaneous frequency of Mode 1 for each set of time series displacements, there is an apparent saturation effect in the frequency distribution, which is a neat value for the maximum instantaneous frequency, as shown in Figures 11 and 12.The reason for this is that it is a periodic oscillating manner for the structural vibration responses of Zhaozhou Bridge.Most of the resulting adjacent quasi-extreme points (t e i , y e i ) and t e i+1 , y e i+1 are extreme points which have the same time intervals.Therefore, according to Step 2 of the direct interpolation algorithm, most of the instantaneous frequencies of the discrete points are calculated by f i−1 = 1/ t e i − t e i−2 , and are also the maximum instantaneous frequencies.Therefore, there is an apparent saturation effect in the frequency distribution of the instantaneous frequency of Mode 1 for each set of time series displacements.These results further indicate that the projected displacement curve of Mode 1 is the main projected displacement variation related to the bridge.(2) Instantaneous frequency is generally a transient structural vibration response which depends on the structural natural frequency, the damping, the stiffness, and the excitation conditions [20].Therefore, in this study, for the time series displacements acquired only with environmental excitation factors, the maximum instantaneous frequency is about 2.49 Hz for point Rbin46_N (Figure 12a) and point Rbin64_N (Figure 12c).However, when affected by the transient load of a car which is a stronger excitation condition, the maximum instantaneous frequency of Mode 1 of point Rbin46_L (Figure 12b) and point Rbin64_L (Figure 12d) increases to 3.37 Hz. (3) Generally speaking, if some damage has occurred in a structure, this will cause a reduction in the structural natural frequency, which can be reflected by the reduction of the instantaneous frequency of the structures for HHT [33].In this study, for the direct interpolation algorithm, if some damage had occurred in Zhaozhou Bridge with the transient load of a car, the periodic oscillating manner of the structural vibration response of Zhaozhou Bridge would have been affected.In other words, for the time series displacements acquired after the transient load of a car, the number of adjacent quasi-extreme points being extreme points would have been greatly reduced.Therefore, according to Step 2 of the direct interpolation algorithm, the number of instantaneous frequencies being the maximum instantaneous frequency should also be greatly reduced.However, as shown in Figure 12b,d, there is an approximately uniform distribution for the maximum instantaneous frequency of Mode 1.This indicates that Zhaozhou Bridge was in an instantaneous stable state when the car passed over Zhaozhou Bridge.(4) The maximum instantaneous frequency of Mode 1 is obtained in conjunction with most of the lower amplitudes.However, there are some decreased instantaneous frequencies of Mode 1 corresponding to the higher amplitudes.The reason for this is that the higher amplitudes are caused by the environmental excitation factors or the transient load of a car, which can affect the periodic oscillating manner of Zhaozhou Bridge.Furthermore, the maximum amplitude (A1) of Mode 1 for point Rbin46_N is about 0.34 mm, which is a little bigger than the 0.29 mm for point Rbin64_N.However, the maximum amplitude (A1) of Mode 1 for point Rbin46_L increases to 0.55 mm, which is also bigger than the 0.46 mm for point Rbin64_L.These results further indicate that the main stress point is located at the 1/4 span point rather than the mid-span point.
Figure 13 shows the instantaneous frequency of Mode 1 of point Rbin64_L obtained by HHT [18].Compared with the instantaneous frequency of Mode 1 of point Rbin64_L obtained by the ESMD method, as shown in Figure 12c, most of the peaks of instantaneous frequency are lower than 3.37 Hz (red line in Figure 13), which is the constant maximum instantaneous frequency of Mode 1 of point Rbin64_L obtained by the ESMD method, as shown in Figure 12d.Furthermore, there are no sudden changes in the curve of the obtained instantaneous frequency of Mode 1 of point Rbin64_L obtained by HHT, which also indicates that Zhaozhou Bridge was in a steady state when the car passed over the bridge.However, there are some larger differences among the peaks of the instantaneous frequency for Mode 1 of point Rbin64_L obtained by HHT.In other words, it is difficult to identify the minor damage if some smaller sudden changes exist in the curve of the instantaneous frequency obtained by HHT.However, for the ESMD method, because the acquired maximum instantaneous frequency is constant for the undamaged parts of the structure, it is easy to detect the damage, which is an advantage of the ESMD method.

Discussion
This paper explored a non-contact approach for the joint use of ground-based microwave interferometry and the ESMD method, with the purpose of the instantaneous vibration analysis of ancient bridges.Ground-based microwave interferometry was adopted to obtain the accurate time series displacements of the monitored ancient bridge, and the ESMD method was adopted to perform the instantaneous vibration analysis of the monitored ancient bridge by the decomposed IMFs and the corresponding instantaneous frequencies.Based on the results of the aforementioned experiment for Zhaozhou Bridge, there are three issues that need to be discussed, as follows.
(1) By adopting IBIS-S to obtain the time series displacements of Zhanzhou Bridge, ground-based microwave interferometry has been validated as an effective non-contact technique to monitor ancient bridges.Without passive radar reflectors installed on the surface of ancient bridges, the thermal SNR values of all of the monitoring points were larger than 35 dB, which could ensure that accurate displacements were obtained.Therefore, for the purpose of avoiding any damages on Zhaozhou Bridge, the traditional transducers were not attached on the bridge to evaluate the accuracy of ground-based microwave interferometry.However, as ground-based microwave interferometry makes use of the propagation of electromagnetic waves for displacement measurement, the measurements performed are inevitably influenced by the atmospheric refractive index, which is caused by temperature, humidity, air pressure, and other meteorological factors [38].Different meteorological conditions have different atmospheric refractive index values, which can cause some loss of measurement accuracy.Therefore, if there is a need for the periodic monitoring of bridges in different measuring times using ground-based microwave interferometry, it is necessary to improve the measurement accuracy of the dynamic deflection of bridges by the use of atmospheric parameters correction, such as the Permanent Scatterers technique [38][39][40] and model-based approach [10,41].This can improve the contrast of different period data.However, in this study, the acquisition time of all data only lasts about 30 min, and has almost the same influence of meteorological factors for displacement measurement.Moreover, the purpose of this paper is to perform the instantaneous vibration analysis of ancient bridges.Hence, it is not necessary to compensate for the propagation losses in this study, which will not affect the reliability of the instantaneous vibration analysis of ancient bridges.(2) Generally, any one complicated signal can be regarded as being composed of multiple simple signals which represent different physical meanings [32].In this study, through the use of the ESMD method, the acquired time series displacements was decomposed into eight IMFs together with an optimal AGM curve, which are more accurate compared with the decomposed 12 IMFs by the HHT-EMD method.Furthermore, by analyzing the decomposed IMFs, the overall vibration trend and sudden variations of the projected displacement can be obtained to analyze the structural characteristics of ancient bridges.However, it is difficult to obtain the physical meaning of each IMF, which needs to be further studied.Moreover, the thrust of wind, ground-motion, complicated traffic, etc., will inevitably generate a noise signal that reduces the accuracy of the measured dynamic displacement of the bridge.In addition, the monitored bridge itself also has a periodic vibration, of which the vibration frequency is different from the transient vibration by the instant loading of a car.Therefore, the obtained dynamic displacement of the bridge, caused by the instant loading of a car, should become bigger or smaller.To the authors' best knowledge, there is not a theoretical model or a gold standard method to determine the expected dynamic displacement in the bridge caused by the instant loading of a car.Further study of the denoising algorithm is also needed to improve the accuracy of the expected dynamic displacement caused by the instant loading.Nevertheless, the above factors influence the expected dynamic displacement caused by the instant loading a little, which will not affect the instantaneous vibration analysis in this study.
(3) Instantaneous frequency is a transient structural vibration response, which not only depends on the structural natural frequency, but is also influenced by the damping, the stiffness, and the excitation conditions [20,42].In this study, the obtained maximum instantaneous frequency is 2.49 Hz for Mode 1 only with environmental excitation factors.However, due to the stronger excitation condition of a transient load of a car, the obtained maximum instantaneous frequency increases to 3.37 Hz for Mode 1.Therefore, by using the ESMD method, it can accurately identify the variations of instantaneous frequency by a neat value for the maximum instantaneous frequency.Furthermore, it is generally known that the structural natural frequency will be reduced if some damage occurs in a structure, which can be reflected by the sudden reduction of the instantaneous frequency of the monitored structure [43,44].In this study, according to the direct interpolation algorithm, if some damage had occurred in Zhaozhou Bridge with the transient load of a car, the instantaneous frequencies of Mode 1 should be decreased steadily.Therefore, although the maximum instantaneous frequency changed from 2.49 Hz to 3.37 Hz for Mode 1, there was no sudden steady decrease in the curve of instantaneous frequency, which indicates that it was in an instantaneous stable state when the car passed over Zhaozhou Bridge.However, Zhaozhou Bridge has been operational for 1400 years now, so if we want to detect the significant changes for the purpose of a global stability analysis, a large period of several months may be required.

Conclusions
In order to evaluate the instantaneous dynamic response of ancient bridges, this paper proposed an integrated method consisting of ground-based microwave interferometry and the ESMD method and applied it to the well-known Zhaozhou Bridge.Ground-based microwave interferometry was adopted to acquire the dynamic time series displacements with environmental excitation factors and a transient load with a car, respectively.The ESMD method was used to decompose the time series displacements into a series of IMFs and the corresponding instantaneous frequencies, which were then used to perform the instantaneous vibration analysis of Zhaozhou Bridge.More specifically, the results presented in the paper clearly highlight that: (1) In this study, aiming to avoid damage to the great historical heritage for Zhaozhou Bridge, the IBIS-S instrument was only located on one side of Zhaozhou Bridge without corner reflectors attached on the lower surface of Zhaozhou Bridge.The quick and easy installation of the IBIS-S instrument can greatly improve the efficiency of data collection.Moreover, the resulting thermal SNR of all of monitoring points on the lower surface of the bridge were larger than 35 dB, which could ensure that accurate displacements were obtained.Therefore, these results verify the feasibility and accuracy of the dynamic monitoring of Zhaozhou Bridge by the sensing method of ground-based microwave interferometry in the paper, which further indicates that ground-based microwave interferometry is a viable alternative technique to acquire dynamic time series displacements for the instantaneous vibration analysis of ancient bridges.Meanwhile, it can also reduce the inherent risk involved with the placement of the traditional contact transducers.(2) The ESMD method was performed to yield a series of IMFs together with an optimal AGM curve through the use of a mode symmetric about the maxima and minima points.The decomposed IMFs can reflect the overall tendency of the projected displacement according to the magnitude of the frequency.Furthermore, they can also reflect the instantaneous dynamic response of the different monitored points.(3) The instantaneous frequencies were obtained using the direct interpolation algorithm, which can reconcile the conflict between the period and the frequency, compared with the traditional time-frequency analysis methods.The instantaneous frequencies of the decomposed IMFs of each set of time series displacements showed that Zhaozhou Bridge was in a steady state when the car passed over the bridge.
shows a sketch of Zhaozhou Bridge.The bridge is about 50.82 m long, with a central span of 37.47 m, which was the world's longest arch bridge at the time of construction.The bridge stands 7.3 m tall and has a width of 9 m.The arch covers a circular segment of less than half of a semicircle (84 • ), with a radius of 27.27 m [25-27].Remote Sens. 2018, 10, x FOR PEER REVIEW 4 of 25 Remote Sens. 2018, 10, x; doi: FOR PEER REVIEW www.mdpi.com/journal/remotesensingZhaozhou Bridge.The bridge is about 50.82 m long, with a central span of 37.47 m, which was the world's longest arch bridge at the time of construction.The bridge stands 7.3 m tall and has a width of 9 m.The arch covers a circular segment of less than half of a semicircle (84°), with a radius of 27.27 m [25-27].
Remote Sens. 2018, 10, x FOR PEER REVIEW 5 of 25 Remote Sens. 2018, 10, x; doi: FOR PEER REVIEW www.mdpi.com/journal/remotesensingthesampling interval of 0.005 s is well suited to providing a good wave definition for the acquired time series displacements[6].

Figure 2 .
Figure 2. View of the IBIS-S instrument for dynamic time series displacement acquisition at Zhaozhou Bridge.

Figure 2 .
Figure 2. View of the IBIS-S instrument for dynamic time series displacement acquisition at Zhaozhou Bridge.

Figure 3 .
Figure3.The whole decomposition process of the ESMD method to yield a series of IMFs together with an optimal AGM curve.

Figure 3 .
Figure3.The whole decomposition process of the ESMD method to yield a series of IMFs together with an optimal AGM curve.

Figure 4 .
Figure 4.The simulation signal ) (t s and its component signals

3 ,
is decomposed into five IMFs.The decomposed IMF2 is close to ) (t s 1 even though some local errors exist due to mode mixing problems.However, in the posterior halves of IMF3 and IMF4, they are very different from component signals ) of some overshoot and undershoot problems.As shown in Figure 5b, through the use of the ESMD method, the simulation signal ) (t s is decomposed into four IMFs, which are equal to the number of component signals.Furthermore, the mode mixing problems and overshoot and undershoot problems almost disappear, and IMF2 to IMF4 accord with component signals ) indicating that the ESMD method performs better than the HHT-EMD method in decomposing a signal into a series of physically meaningful representations.

Figure 5 .
Figure 5. IMFs of ) (t s decomposed by the HHT-EMD method and the ESMD method.(a) IMFs of ) (t s decomposed by the HHT-EMD method and the exhibited component signals along with their corresponding IMFs; and (b) IMFs of ) (t s decomposed by the ESMD method and the exhibited component signals along with their corresponding IMFs.

Figure 4 .
Figure 4.The simulation signal ) (t s and its component signals

3 3 ,
is decomposed into five IMFs.The decomposed IMF2 is close to ) (t s 1 even though some local errors exist due to mode mixing problems.However, in the posterior halves of IMF3 and IMF4, they are very different from component signals ) on account of some overshoot and undershoot problems.As shown in Figure 5b, through the use of the ESMD method, the simulation signal ) (t s is decomposed into four IMFs, which are equal to the number of component signals.Furthermore, the mode mixing problems and overshoot and undershoot problems almost disappear, and IMF2 to IMF4 accord with component signals ) indicating that the ESMD method performs better than the HHT-EMD method in decomposing a signal into a series of physically meaningful representations.

Figure 5 .
Figure 5. IMFs of s(t) decomposed by the HHT-EMD method and the ESMD method.(a) IMFs of s(t) decomposed by the HHT-EMD method and the exhibited component signals along with their corresponding IMFs; and (b) IMFs of s(t) decomposed by the ESMD method and the exhibited component signals along with their corresponding IMFs.

Figure 6 .
Figure 6.Curves of the thermal SNR of the received signals.(a) Curve of the thermal SNR of the entire monitoring range; (b) Enlarged curve of the thermal SNR of points Rbin46 and Rbin64.

Figure 6 .
Figure 6.Curves of the thermal SNR of the received signals.(a) Curve of the thermal SNR of the entire monitoring range; (b) Enlarged curve of the thermal SNR of points Rbin46 and Rbin64.

Figure 7 .
Figure 7.The original time series displacements and the corresponding optimal AGM curves.(a) The original time series displacements and the optimal AGM curve of point Rbin46_N; (b) The original time series displacements and the optimal AGM curve of point Rbin46_L; (c) The original time series displacements and the optimal AGM curve of point Rbin64_N; (d) The original time series displacements and the optimal AGM curve of point Rbin64_L.

Figure 7 .
Figure 7.The original time series displacements and the corresponding optimal AGM curves.(a) The original time series displacements and the optimal AGM curve of point Rbin46_N; (b) The original time series displacements and the optimal AGM curve of point Rbin46_L; (c) The original time series displacements and the optimal AGM curve of point Rbin64_N; (d) The original time series displacements and the optimal AGM curve of point Rbin64_L.

Figure 8 .
Figure 8. Variance ratio of the original time series displacements.(a) Variance ratio of point Rbin46_N; (b) Variance ratio of point Rbin46_L; (c) Variance ratio of point Rbin64_N; (d) Variance ratio of point Rbin64_L.

Figure 8 .
Figure 8. Variance ratio of the original time series displacements.(a) Variance ratio of point Rbin46_N; (b) Variance ratio of point Rbin46_L; (c) Variance ratio of point Rbin64_N; (d) Variance ratio of point Rbin64_L.
Remote Sens. 2018, 10, x FOR PEER REVIEW 16 of 25 Remote Sens. 2018, 10, x; doi: FOR PEER REVIEW www.mdpi.com/journal/remotesensingthan the eight IMFs obtained via the ESMD method, and the displacement curve trends of Mode 5 to Mode 12 are coarse compared with those of Mode 5 to Mode 8 decomposed by the ESMD method.

Figure 10 .
Figure 10.IMF components of point Rbin64_L decomposed via the HHT-EMD method.

Figure 10 .
Figure 10.IMF components of point Rbin64_L decomposed via the HHT-EMD method.

Figure 11 .
Figure 11.Frequency distributions of eight modes of each set of time series displacements obtained by the direct interpolation algorithm.(a) Frequency distribution of each mode of point Rbin46_N; (b) Frequency distribution of each mode of point Rbin46_L; (c) Frequency distribution of each mode of point Rbin64_N; (d) Frequency distribution of each mode of point Rbin64_L.

Figure 11 .
Figure 11.Frequency distributions of eight modes of each set of time series displacements obtained by the direct interpolation algorithm.(a) Frequency distribution of each mode of point Rbin46_N; (b) Frequency distribution of each mode of point Rbin46_L; (c) Frequency distribution of each mode of point Rbin64_N; (d) Frequency distribution of each mode of point Rbin64_L.
Remote Sens. 2018, 10, x FOR PEER REVIEW 20 of 25 Remote Sens. 2018, 10, x; doi: FOR PEER REVIEW www.mdpi.com/journal/remotesensingmaximuminstantaneous frequency is constant for the undamaged parts of the structure, it is easy to detect the damage, which is an advantage of the ESMD method.

Figure 12 .
Figure 12.Instantaneous frequency and amplitude of Mode 1 decomposed via the ESMD method for time series displacements by the direct interpolation algorithm.(a) Instantaneous frequency and amplitude of Mode 1 of point Rbin46_N; (b) Instantaneous frequency and amplitude of Mode 1 of point Rbin46_L; (c) Instantaneous frequency and amplitude of Mode 1 of point Rbin64_N; (d) Instantaneous frequency and amplitude of Mode 1 of point Rbin64_L.

Figure 13 .
Figure 13.Instantaneous frequency of Mode 1 of point Rbin64_L obtained by HHT.

Figure 12 .
Figure 12.Instantaneous frequency and amplitude of Mode 1 decomposed via the ESMD method for time series displacements by the direct interpolation algorithm.(a) Instantaneous frequency and amplitude of Mode 1 of point Rbin46_N; (b) Instantaneous frequency and amplitude of Mode 1 of point Rbin46_L; (c) Instantaneous frequency and amplitude of Mode 1 of point Rbin64_N; (d) Instantaneous frequency and amplitude of Mode 1 of point Rbin64_L.

Figure 12 .
Figure 12.Instantaneous frequency and amplitude of Mode 1 decomposed via the ESMD method for time series displacements by the direct interpolation algorithm.(a) Instantaneous frequency and amplitude of Mode 1 of point Rbin46_N; (b) Instantaneous frequency and amplitude of Mode 1 of point Rbin46_L; (c) Instantaneous frequency and amplitude of Mode 1 of point Rbin64_N; (d) Instantaneous frequency and amplitude of Mode 1 of point Rbin64_L.

Figure 13 .
Figure 13.Instantaneous frequency of Mode 1 of point Rbin64_L obtained by HHT.Figure 13.Instantaneous frequency of Mode 1 of point Rbin64_L obtained by HHT.

Table 1 .
Parameters of IBIS-S for the dynamic time series displacements acquisition of Zhaozhou Bridge.

Table 2 .
The minimum variance ratios and the corresponding sifting times of the four kinds of time series displacements.

Table 2 .
The minimum variance ratios and the corresponding sifting times of the four kinds of time series displacements.

Table 3 .
Parameters of the quality assessment for the frequency distribution of each set of time series displacements.

Table 3 .
Parameters of the quality assessment for the frequency distribution of each set of time series displacements.