Open AccessArticle

Integration of Single-Frequency GNSS and Strong-Motion Observations for Real-Time Earthquake Monitoring

Rui Tu

^1,2,3,*,

Rui Zhang

^1,3,

Pengfei Zhang

^1,2,

Jinhai Liu

^1,2 and

Xiaochun Lu

^1,2,3

National Time Service Center, Chinese Academy of Sciences, Shu Yuan Road, Xi’an 710600, China

University of Chinese Academy of Sciences, Yu Quan Road, Beijing 100049, China

Key Laboratory of Precision Navigation Positioning and Timing Technology, Chinese Academy of Sciences, Xi’an 710600, China

Author to whom correspondence should be addressed.

Remote Sens. 2018, 10(6), 886; https://doi.org/10.3390/rs10060886

Submission received: 18 April 2018 / Revised: 26 May 2018 / Accepted: 5 June 2018 / Published: 6 June 2018

(This article belongs to the Special Issue Revisiting the 2008 Mw 7.9 Wenchuan, China Earthquake: New evidence and implications from earth observations)

Download

Browse Figures

Graphical abstract
"> Figure 1
Comparison of average energies for GNSS (bottom) and SM (top) detectors. "> Figure 2
Correction of initial baseline shift for GNSS (top) and strong-motion (bottom). The black line represents raw values; the red line, estimated corrections; and the blue line, corrected values. "> Figure 3
Quality evaluation of the displacement and velocity after an earthquake. The upper graph shows the displacement, which is stable around a fixed offset; the lower graph shows the velocity, which is stable around 0 m/s. "> Figure 4
Real-time integration system of GNSS and SM observations. ‘RT’ indicates real-time, ‘Dis, Vel, and Acc’ represent displacement, velocity and acceleration, respectively, ‘obs’ represents observation, and ‘eph’ represents ephemeris. "> Figure 5
Experimental platform (left) and processing facilities (right) of the combined system GPS and SM system. "> Figure 6
Movement sequence of the GPS and SM experiments. The graph covers the time for eight tests (T1–T8), each indicated by a dashed red rectangle. "> Figure 7
Comparisons of displacement, velocity, and acceleration for the GPS (left); SM (middle); and combined system (COM) (right). "> Figure 8
Comparison of the PSD for the displacements of the first experiment by different sensors. Black, blue and red lines represent the GPS, SM, and the combined system, respectively, the (top) represents no motion and (bottom) represents displacement motion. "> Figure 9
Displacement (left); velocity (middle); and acceleration (right) output of the integration system for the GPS and SM data. Note that vertical scales vary between the tests. "> Figure 10
Experimental equipment used to test the integration of BDS and SM data. Hardware consists of a GNSS sensor, an SM sensor, and a tracked platform. "> Figure 11
Movement sequence of eight BDS and SM experiments. Each green dotted rectangle represents one test. "> Figure 12
Comparison of displacements (top) and PSD (bottom) for different GNSS systems. Black, blue, red, and pink lines represent the GPS, GLONASS, BDS, and their combinations, respectively. "> Figure 13
Displacement (left); velocity (middle); and acceleration (right) output of the integration system for the BDS and SM data. Black curves indicate east–west displacement of the test bed, whereas the red curves indicate north–south motion. "> Figure 14
Displacement (left); velocity (middle); and acceleration (right) output of the integration system for the Baja earthquake data. The upper three plots represent pairs of ‘P496-5058’; the lower three plots represent pairs of ‘P744-5028’. "> Figure 15
Differences of displacement waveforms between GPS and combined solutions. The (left) represents pairs of ‘P496-5058’; the (right) represents pairs of ‘P744-5028’. "> Figure 16
PSD of displacement waveforms by different methods (‘‘COM” represents the GPS and SM combined results, ‘‘DIF” represents the difference between GPS and combined results). The (left) represents pairs of ‘P496-5058’; the (right) represents pairs of ‘P744-5028’. ">

Versions Notes

Abstract

In this study, a real-time earthquake monitoring system based on the integration of single-frequency global navigation satellite system (GNSS) and strong motion (SM) observations was developed. This high-precision integrated system can provide full-frequency monitoring information, and it makes full use of SM data to quickly and accurately determine the vibration window for initial baseline shift correction. High-precision displacement data obtained from GNSS epoch-differenced velocity estimation are used to constrain the SM’s low-frequency baseline shift. Hence, full-frequency monitoring information (displacement, velocity, and acceleration) can be provided in real-time. Three different datasets were used for validation and the results confirm that the proposed system can be used for practical earthquake monitoring.

Keywords:

real time; high-rate GNSS; strong-motion; integration system; earthquake monitoring

Graphical Abstract

1. Introduction

Epoch-differenced velocity estimation technology, which is based on single-frequency global navigation satellite system (GNSS) observations and satellite broadcast ephemeris, has the advantages of no ambiguity parameters, rapid initialization, high precision, low cost, and rapid turnaround; therefore, it is widely used in earthquake engineering applications, such as real-time earthquake monitoring, early warning, and rapid response [1]. In addition, other studies have focused on GNSS positioning from mass-market receivers for monitoring purposes and real-time earthquake monitoring activities [2,3,4,5]. Meng et al. used GNSS positioning technology for structural health monitoring and geo-hazard early warning [2], Bellone et al. and Dabove and Li et al. also used the low-cost receiver for fast deformation monitoring (e.g., earthquakes) [3,4,5]. However, owing to its low sampling rate, the GNSS can only provide low-frequency displacement information with high precision; high-frequency information, such as velocity and acceleration, is significantly affected by noise [6,7,8].

Strong motion (SM) data, which is inexpensive, can provide high-frequency information with high precision and high sensitivity; however, the integrated velocity and displacement are seriously distorted owing to a baseline shift caused by low-frequency noise and tilting or rotation of the instrument; consequently, the recovered displacement has a precision of the order of decimeters or worse [9,10,11,12,13,14].

The effective integration of GNSS and SM data to maximise their complementary advantages is a topic of significant research interest [14,15]. The goal is to develop an early warning system that not only provides high-precision and full-frequency information, but is also inexpensive and efficient. Wang et al. used linear trigonometric and time correlation to express the baseline shift of SM data and to correct the displacement to be consistent with the high-frequency GNSS displacement using a least-squares solution [14]. This method constrains the baseline shift and results in high-precision velocity and displacement measurements; however, the assumptions of the linear function are influenced by subjective factors and, therefore, it cannot be employed for real-time applications. Bock et al. used a Kalman filter to realise the combination of GNSS and SM [15]. The inputs are the GNSS’ displacement and the SM’s acceleration; the outputs are the filtered displacement and velocity. The Kalman filter approach can result in high precision and broadband displacement information. However, GNSS uses a precise positioning mode that requires a dual-frequency receiver and high-precision satellite ephemeris, with associated expense and poor turnaround time; hence, the widespread implementation of this approach is limited. Tu et al. proposed an approach to compare the difference between the global positioning system (GPS) displacement and SM-integrated displacement to derive the baseline shift, and to subsequently recover the final displacement and velocity [16]. However, validation trials on an experimental dataset indicate that it worked only for 50 Hz GPS data, which limits the method’s practical application. Subsequently, many researchers have studied the tight integration of high-quality dual-frequency GNSS and SM observations for deformation monitoring and earthquake early warning systems [16,17,18,19,20,21,22,23,24,25,26], which have the advantages of being complementary systems; however, they have the drawback of significant cost. Only a few researchers have strongly considered the cost issues (e.g., Saunders et al. [26]) by using single-frequency GPS and low-cost MEMS accelerometers for earthquake early warning and rapid response.

In this context, our study proposes a real-time earthquake monitoring system based on the integration of a single-frequency GNSS and low-cost SM observations. The system is composed of three modules, all of which work in real time: data input, processing, and output. The system uses a Kalman filter method to combine the complementary systems and allow real-time retrieval of velocity and acceleration information for applications in earthquake monitoring and early warning.

2. Methodology

2.1. Determination of the Signal Window

High-frequency GNSS has issues such as high levels of noise, time-consuming separation of the true vibration signal from the noise, and poor real-time sensitivity. However, high-sensitivity SM detectors can determine the signal window quickly and accurately, which is useful for initialization and quality inspection.

A seismic event can be divided into three periods: pre-seismic, co-seismic, and post-seismic. For the time series of the co-seismic acceleration, we define an average energy function:

W_{n} = \frac{1}{N} (\sum_{t = n + 1}^{t = n + N} | a (t) | \cdot Δ t^{2}), W_{0}

(1)

where

W_{n}

is the average energy at epoch

n

W_{0}

is the rest energy,

N

is the epoch numbers to calculate the average energy,

a (t)

is the acceleration time series, and

Δ t

is the sample interval. The signal start time

t_{1}

is defined as the time when the average energy

W_{n}

is more than five times the rest energy

W_{0}

; the end time

t_{2}

is defined as the time after the start time when the average energy

W_{n}

is less than five times the rest energy

W_{0}

Figure 1 shows a comparison of average energies for GNSS and SM. The GNSS energy function has high levels of noise and, thus, cannot be used to accurately determine the signal start and end times. Consequently, it is unsuitable for real-time initialization of the GNSS velocity shifts, making implementation for earthquake monitoring and early warning in real-time difficult. However, the shaking window can be readily determined from the SM energy function, and it complements the GNSS system for the initialization shift before the earthquake and quality evaluation after the earthquake.

2.2. Initial Baseline Shift Correction of the Integrated System before an Earthquake

Colosimo et al. and Tu et al. reported a linear shift in the estimated velocity in epoch difference calculations due to uncorrected atmosphere errors; this shift is stable over a period of around 10 min [1,16]. One to two minutes of velocity data prior to the earthquakes was selected to determine the initial shift by assuming that the time series of the GNSS velocity

v (t)

can be expressed as:

v (t) = (v_{1}, v_{2}, v_{3} \dots v_{N}) .

(2)

Given the presence of high-frequency noise, we usually integrate the velocity to obtain low-frequency displacement and use least-squares fitting to estimate the initial shift. Equations (3) and (4) are the solution and correction, respectively, for the displacement:

s (t) = (s_{1}, s_{2}, s_{3} \dots s_{N}) = {At}^{2} + Bt + C,

(3)

v_{c} (t) = v (t) - (2 \cdot \hat{A} t + \hat{B}),

(4)

where

s (t)

is the integrated displacement time series,

A, B, and C

are the polynomial fitting coefficients, and

v_{c} (t)

is the corrected velocity.

The SM acceleration is corrected by a constant initial baseline shift to compensate for environmental noise and instrument effects. The time series of the acceleration can be written as:

a (t) = (a_{1}, a_{2}, a_{3} \dots a_{N})

(5)

Similarly, the acceleration integrated to velocity or displacement used least-squares fitting to solve the acceleration initial shift correction:

a_{c} (t) = a (t) - \hat{a} (t),

(6)

where

a_{c} (t)

is the corrected acceleration and

\hat{a} (t)

is the initial shift correction for acceleration. Figure 2 shows the GNSS and SM system’s initial baseline shift corrections. After correction, the GNSS velocity data fluctuate around the zero line, implying that the initial baseline shift correction model is objective and accurate.

2.3. Integration during an Earthquake

Low-frequency displacement can be accurately recovered from the GNSS single system’s velocity estimation and the initial shift correction [19]:

s_{c} (t) = \sum_{t = 1}^{t = N} v_{c} (t) \cdot Δ t = (s_{c 1}, s_{c 2}, s_{c 3} \dots s_{cN}) .

(7)

For the SM sensors, there is an additional instantaneous shift and new permanent shift due to the instrument’s tilting or rotation during the shaking. Although they can significantly distort the integrated velocity and displacement, the magnitude of these shifts is small (about 0.001

m / s^{2}

) and can be ignored. An accurate acceleration can be obtained from:

a_{c} (t) = (a_{c 1}, a_{c 2}, a_{c 3} \dots a_{cN}) .

(8)

The integrated system used a Kalman filter to determine velocity [15]. The GNSS’ displacement and SM’s acceleration were used for the inputs; the output is the accurate velocity information. The observation equations are:

Y (t) = A \cdot X (t) + ε_{Y},

(9)

Y (t) = [\begin{matrix} Y_{s (t)} \\ Y_{v (t)} \\ Y_{a (t)} \end{matrix}],

(10)

A = [\begin{matrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{matrix}],

(11)

X (t) = [\begin{matrix} s (t) \\ v (t) \\ a (t) \end{matrix}],

(12)

ε_{Y} \sim N (0, Q_{Y}),

(13)

P = Q_{Y}^{- 1} = [\begin{matrix} \frac{1}{q_{s}} \\ \frac{1}{q_{v}} \\ \frac{1}{q_{a}} \end{matrix}],

(14)

where

Y (t)

denote the observation values,

A

is the coefficient matrix,

X (t)

denote the solution parameters,

ε_{Y}

is the observation noise, and

Q_{Y}

is the noise variance matrix,

P

is the observation weight matrix, and

q_{a}

q_{v}

, and

q_{s}

are the noise variance of acceleration, velocity, and displacement, respectively.

The corresponding state equation can be written as:

X (t) = Φ \cdot X (t - 1) + ε_{s},

(15)

Φ = [\begin{matrix} 1 & τ & \frac{τ^{2}}{2} \\ 0 & 1 & τ \\ 0 & 0 & 1 \end{matrix}],

(16)

ε_{s} \sim N (0, Q_{s}),

(17)

Q_{s} = τ q_{a}^{2} [\begin{matrix} \frac{τ^{4}}{20} & \frac{τ^{3}}{8} & \frac{τ^{2}}{6} \\ \frac{τ^{3}}{8} & \frac{τ^{2}}{3} & \frac{τ}{2} \\ \frac{τ^{2}}{6} & \frac{τ}{2} & 1 \end{matrix}],

(18)

where

Φ

is the state transition matrix,

ε_{s}

is the dynamic noise,

Q_{s}

is the dynamic noise matrix, and

τ

is the sample interval of the acceleration.

The prediction estimation can be expressed as:

\bar{X (t)} = Φ \cdot \hat{X} (t - 1),

(19)

\bar{Q (t)} = Φ (t - 1) \cdot Q (t - 1) \cdot Φ^{T} (t - 1) + Q_{s} (t - 1),

(20)

and the filter estimation can be expressed as:

Q_{Y} (t) = P^{- 1} + A \cdot \bar{Q (t)} \cdot A^{T},

(21)

Q (t) = {({\bar{Q (t)}}^{- 1} + A^{T} \cdot Q_{Y}^{- 1} (t) \cdot A)}^{- 1},

(22)

\hat{X} (t) = \bar{X (t)} + Q (t) \cdot A^{T} \cdot Q_{Y}^{- 1} (t) \cdot (Y (t) - A \cdot \bar{X (t)}) .

(23)

In the real-time filtering system, the predict estimation was carried out on the SM sample point, and the filter estimation was performed on the GNSS sample point. Hence, the integrated system can obtain high-precision displacement, velocity, and acceleration information in real-time.

2.4. Quality Assessment of the Integrated System

SM sensors can determine the signal end time quickly and accurately after shaking, perform quality control, and prepare for the next initialization. The inspection was performed separately for the GNSS system, where the time series of velocity were integrated into displacement to determine whether the permanent displacement was stable based on a criterion of standard deviation (STD) of less than 2 cm. On the other hand, for SM, the solved velocity is assessed by whether it returns to the zero level after shaking stops, using the threshold of STD less than 1 cm/s as the criterion.

Figure 3 is an example for the quality evaluation of the integrated system. The permanent displacement is stable and the velocity is on the zero level after the event, both of which indicate that the baseline shift is appropriate.

2.5. Implementation of the Real-Time Integrated System

The system is divided into three modules: real-time data stream input, real-time data processing, and output, it is shown in Figure 4.

(1): The real-time data stream input module collects SM acceleration data, GNSS observation data, satellite broadcast ephemeris products, and associated parameters.
(2): The real-time data processing module include pre-processing, integration, and quality control. During the pre-process period, after the SM determines the signal start, each subsystem carries out an initial shift correction. Then, the GNSS system provides accurate displacement information and the SM system provides accurate acceleration information. In the integration process period, the GNSS displacement and SM acceleration are the system inputs to the Kalman filter process; the accurate and low-frequency displacement can effectively constrain the SM’s baseline shift and retrieve high-precision velocity information. In the quality control module, the integrated system’s velocity and displacement are inspected; the accuracy and reliability are also verified.
(3): The real-time output module controls the output of the displacement, velocity, acceleration, and accuracy information.

3. Validation Results and Analysis

The integrated system was tested and validated using three different datasets: (1) an experimental dataset including GPS and SM observations collected simultaneously on a dedicated experimental platform; (2) an experimental dataset with the BeiDou Satellite System (BDS) and SM observations collected simultaneously in a dedicated experimental setup; and (3) GPS and SM data collected from the E1 Mayor-Cucapah earthquake. All three datasets were processed in a simulated real-time mode by only using data up to the processed epoch and L1 single-frequency GNSS observations; velocity and acceleration were analysed by the combined system. In addition, all the scenarios which are examined in the experiments are not the transient component, but the co-seismic component of the ground response.

3.1. GPS and SM Data Test

Figure 5 shows the GPS and SM experimental platform, composed of a dynamic GPS antenna (Type: JAV_RINGANT_G3T NONE), a low-cost MEMS-type accelerometer [27], and real-time processing facilities. The experiment was carried out in December 2012; the sample rates were 50 Hz and 100 Hz for the single-frequency GPS and SM sensor, respectively. Real permanent displacement was provided by the Vernier calipers and the high-frequency camera to the front of the rig, and it was measured with a precision of better than 5 mm. Data processing used a high-frequency GPS and SM integration solution package developed by our research group.

Figure 6 shows the displacement sequence that was used for the test. We simulated eight events (T1 to T8) by moving the platform from one side to the other. The reference displacements were recorded by a set of Vernier calipers fixed on the table. For each experiment, we initially kept the rig static at the start point for about 5 min; then, we slid the combined instruments from the start point to the end point, for over about 1 min; next, the rig was kept static at the end point for about 5 min; finally, the rig was returned to the start point for the next experiment. For the eighth experiment, it was not returned to the start point, as this was the end of the experiment. Owing to the design of the platform, movement was performed only in one direction along the rail.

Figure 7 compares the displacement, velocity, and acceleration between GPS and SM individually and in combination for the first test. It can be concluded that the single-frequency GPS velocity estimation, which is based on broadcast ephemeris, can recover high-precision low-frequency displacement information; however, the velocity and acceleration are seriously distorted because of high-frequency noise pollution. In comparison, the SM-based acceleration has a much higher signal-to-noise ratio; however, the baseline-shift-corrected velocity and displacement show larger offsets than those for the GPS, while for the combination of GPS and SM measurements, both velocity and displacement are recovered in addition to the combined displacement, which is recovered with much less noise than the GNSS single system.

Figure 8 shows the power spectral densities (PSD) of the displacements obtained using different sensors; the top represents no motion, and the bottom represents the displacement motion. The GPS and SM combined system complements the advantages of each other, resulting in much smaller noise and much higher frequency than the GPS, and the removal of the baseline shift error in the SM.

Figure 9 shows the displacement, velocity, and acceleration outputs of the integration system, conforming that it accurately corrects the initial baseline shift of each system. Hence, the combined system effectively constrains the instantaneous complex baseline shift during an earthquake allowing high-precision broadband information to be retrieved in real-time.

Table 1 presents the final displacement and velocity biases compared to the references and their standard deviation statistics for all eight tests; all STDs are smaller than the thresholds. From these data, the accuracy of the displacement is 1–2

cm

and the velocity is about 1

mm

/s for the GPS and SM combined system.

3.2. BDS and SM Data Test

BDS is a new GNSS system that is operated by China. Currently, it can provide positioning, navigation, and timing services in the Asia-Pacific region, and will expand to provide global coverage in 2020. We tested the performance of the integration of the BDS and SM system. Figure 10 shows another test platform used in the experiment performed in 2016 August in Xi’an, China. The platform, which can slide along a table, includes a dynamic GNSS antenna that can track BDS, GPS, and GLONASS signals (Type: UR380, Antenna: HX-GG486A), and a low-cost MEMS-type accelerometer [27]. The sampling rate was 1 Hz for GNSS and 100 Hz for the accelerometer. The maximum sliding distance of the platform is restricted to about 0.5 m. All data were processed in a simulated real-time mode. Here, we only analysed BDS and SM data.

Figure 11 shows the movement sequence applied to the BDS and SM integration experiment. We simulated eight earthquake events (E1–E8) by moving the platform from one side to the other. As for the previous experiment, reference displacements were recorded by a set of Vernier calipers fixed on the table. For each experiment, we kept the rig static at the start point for several minutes at first; then, we slid the combined instruments from the start point to the endpoint over about 1 or 2 min; and finally, the rig was kept static at the endpoint for several minutes, before starting the next experiment. We carried out eight experiments in the horizontal component: three in a north–south direction, three in an east–west direction, and two in a northeast–southwest direction.

As shown in previous studies, the performance of multi-GNSS systems is better than that of the single GNSS system [28,29,30,31,32]. Figure 12 shows a comparison of displacements and PSD for different GNSS systems for the first experiment. It can be seen that the precision of the GLONASS system was the worst because of fewer observed satellites. The BDS and GPS system nearly show the same accuracy, and their combination results are more stable. The PSD results reflect that BDS has much larger noise than GPS, and that the GLONASS has a small offset; the combination of multi-GNSS systems represents the average of each single system.

In Figure 13, the results of combining BDS and SM data for all eight experiments are illustrated. As expected, when the baseline shifts of the accelerations are corrected by the combined system, the high-resolution acceleration can help constrain the combination solution. Thus, high-precision and broadband deformation information (displacement, velocity, and acceleration) are recovered in real time.

Table 2 summarises the mean bias and STD between the recovered and references values. The RMS is generally less than 2 mm for velocity, and better than 2 cm for displacement.

3.3. Baja Earthquake Data Test

After validation using two experimental datasets, co-located GPS and seismic data collected during the E1 Mayor-Cucapah earthquake (Mw, 7.2, 4 April 2010, 22:40:42 UTC) were also analysed. Two pairs of co-located GPS and SM seismometers were selected and processed. The GPS station P496, which is located 74 km from the earthquake epicentre, is nearly co-located with the seismic station 5058 of the Southern California Seismic Network (SCSN), within a distance of about 140 m. GPS station P744 located 80 km from the epicentre is about 275 m away from the nearest seismic station, 5028 in SCSN. The accelerometer data are 200 Hz, and the GPS data were recorded at 5 Hz.

Figure 14 shows the acceleration, velocity, and displacement information of the combined system in the east, north, and vertical components. The velocity sequence clearly reflects the waveform information, and the velocity fluctuates about the zero line after the earthquake. In the displacement sequence, the displacement fluctuates on the permanent displacement line.

In addition, the difference between GPS and the combined solution (by GPS and SM sensors) are shown in Figure 15. For the displacement, the residuals from the GPS and the integration method are nearly at the same level, the largest difference between the GPS and the combined solution is about 5 cm and most of them are smaller than 2 cm in the horizontal component. Additionally, it is much worse in the vertical component, the reasons may be owing to the GPS atmosphere error and SM baseline errors are not precisely corrected.

Table 3 shows the statistical results for these datasets. The reference displacement is provided by the GPS precise point positioning post solutions, with a precision of better than 5 mm. For the permanent displacement, which is calculated from the GPS post solutions, the mean difference in the horizontal components is 2–3 cm, and the average difference in the vertical components is about 4 cm. Such levels of displacement precision are much better than those of an individual SM seismometer, which has a displacement precision of the order of decimeters. Furthermore, the combined system can recover acceleration and velocity information, which the single GPS system cannot accurately provide. The vertical component of the GPS signal is subject to greater dithering than the horizontal components because the atmosphere error is difficult to process accurately in GPS epoch-differenced data; thus, the precision of the vertical displacement is worse than that of the horizontal components. These GPS errors will be included in the output of the combined system.

Meanwhile, the corresponding PSD are shows in Figure 16, we can see the GPS with much larger noise and narrower band than the combined results. Thus, the combined sensors best complement the advantages of GPS’s high precision and strong-motion’s high frequency.

4. Discussion

GNSS velocity estimation has significant potential benefits as it can be operated with a single station in real-time, and it does not need initialization. Previous studies have shown that 5 Hz GPS data can be used to recover a significant part of PGV values for large earthquakes [33].

For the application of velocity estimation, attention to details is important. As velocity estimation uses broadcast satellite ephemeris and epoch-differences, uncorrected atmosphere errors and ephemeris residuals will lead to shifts in estimated velocity. In this study, we treated these shifts as a linear model; the final precision may not be as high as with a linear correction, but linear corrections are, themselves, limited because they only work for short periods. In a future study, we plan to estimate GNSS shift as unknown parameters. In addition, for the epoch-difference solution, we plan to use the broadcast ephemeris record for the entire event period to avoid possible jumps in the results caused by discontinuities of adjacent ephemeris records.

Using only one kind of sensor has its limitations; the cooperation of different sensors can complement the advantages of each and it can, therefore, provide more reliable and robust information. The integration of multi-GNSS and SM is an interesting area of development for further research.

In this study, we proposed a combination approach by using single-frequency GNSS and low-cost SM observations. However, some of the validation data were not collected from the real single-frequency and/or low-cost sensors; thus, additional datasets are needed for further validation. Our results consider the complementary advantages of each sensor; however, inconsistencies between the two datasets may also depend on the frequency of the excitation, the direction of the excitation signal, and the distance from the excitation source; these need to be analysed in the future [34]. The analysis of noise in GNSS and SM measurements will also represent the focus of future work.

For real-time earthquake monitoring, the trigger mechanism is the most important. As the acceleration variation in SM data is more sensitive than that in GNSS data, we used the energy change from SM to detect the start of the combined system. When an event is detected, we required only a short period of pre-event data for the initialization; thus, when continuous events are detected, this combined system will not work. Furthermore, energy change not induced by earthquakes has no effect on the combined solution, which is processed by each collocated station; stations without influence are used for deformation monitoring, which is not limited to earthquakes.

The data solution of the combined system is based on the detection of an event. As the GNSS’ linear shift by epoch difference is also changing, the data solution is only effective for a short period; thus, it cannot be used for slow deformation hazards.

5. Conclusions

In this study, we developed a high-precision, full-frequency, and real-time earthquake monitoring system based on the integration of single-frequency GNSS and low-cost SM observations. The integrated system makes full use of the high-sensitivity of SM data to quickly and accurately determine the vibration window for the initial baseline shift correction before the earthquake and quality control after an earthquake. Additionally, it uses high-precision displacement recovered from the GNSS epoch-differenced velocity estimation to effectively constrain the SM low-frequency baseline shift. Hence, the integrated system can determine high-precision and full-frequency vibration information in real-time for the application of earthquake monitoring and early warning. Single-frequency GNSS and simple SM sensors are low cost; the acquisition of satellite broadcast ephemeris and observation data is straightforward in real-time; the data processing module can be embedded in the integrated instruments, enabling automatic operation for individual stations or single machines. Consequently, we believe that a real-time earthquake monitoring system based on this approach will be widely used in the future.

Author Contributions

R.T. and P.Z. conceived and designed the experiments; R.T. and R.Z. performed the experiments, analysed the data, and prepared the manuscript; and R.Z., J.L., and X.L. contributed to the discussions and revisions.

Acknowledgments

The work was partly supported by the program of National Key Research and Development Plan of China (Grant No. 2016YFB0501804), National Natural Science Foundation of China (Grant Nos. 41504006, 41674034) and Chinese Academy of Sciences (CAS) programs of ‘Pioneer Hundred Talents’ (Grant No. Y620YC1701), Light of the west joint scholar’ (Grant No. Y507YR1203), and ‘The Frontier Science Research Project’ (Grant No. QYZDB-SSW-DQC028).

Conflicts of Interest

The authors declare no conflict of interest.

References

Colosimo, G.; Crespi, M.; Mazzoni, A. Real-time GPS seismology with a stand-alone receiver: A preliminary feasibility demonstration. J. Geophys. Res. 2011, 116, B11302. [Google Scholar] [CrossRef]
Meng, X. From Structural Health Monitoring to Geo-Hazard Early Warning: An Integrated Approach Using GNSS Positioning Technology. In Earth Observation of Global Changes (EOGC); Springer: Berlin/Heidelberg, Germany, 2013; pp. 285–293. [Google Scholar]
Bellone, T.; Dabove, P.; Manzino, A.M.; Taglioretti, C. Real-time monitoring for fast deformations using GNSS low-cost receivers. Geomat. Nat. Hazards Risk 2016, 7, 458–470. [Google Scholar] [CrossRef]
Dabove, P.; Manzino, A. Kalman Filter as Tool for the Real-time Detection of Fast Displacements by the Use of Low-cost GPS Receivers. In Proceedings of the 2nd International Conference on Geographical Information Systems Theory, Applications and Management, Rome, Italy, 26–27 April 2016; pp. 15–23, ISBN 978-989-758-188-5. [Google Scholar] [CrossRef]
Li, M.; Li, W.; Fang, R.; Shi, C.; Zhao, Q. Real-time high-precision earthquake monitoring using single-frequency GPS receivers. GPS Solut. 2015, 19, 27–35. [Google Scholar] [CrossRef]
Elo’segui, P.; Davis, J.L.; Oberlander, D.; Baena, R.; Ekström, G. Accuracy of high-rate GPS for seismology. Geophys. Res. Lett. 2006, 33, L11308. [Google Scholar] [CrossRef]
Genrich, J.F.; Bock, Y. Instantaneous geodetic positioning with 10–50 Hz GPS measurements: Noise characteristics and implications for monitoring networks. J. Geophys. Res. 2006, 111, B03403. [Google Scholar] [CrossRef]
Larson, K.M.; Bilich, A.; Axelrad, P. Improving the precision of high-rate GPS. J. Geophys. Res. 2007, 112, B05422. [Google Scholar] [CrossRef]
Iwan, W.; Moser, M.; Peng, C. Some observations on strong-motion earthquake measurement using a digital acceleration. Bull. Seismol. Soc. Am. 1985, 75, 1225–1246. [Google Scholar]
Boore, D.M. Effect of baseline corrections on displacement and response spectra for several recordings of the 1999 Chi-Chi, Taiwan, earthquake. Bull. Seismol. Soc. Am. 2001, 91, 1199–1211. [Google Scholar] [CrossRef]
Zhu, L. Recovering permanent displacements from seismic records of the June 9, 1994 Bolivia deep earthquake. Geophys. Res. Lett. 2003, 30, 1740. [Google Scholar] [CrossRef]
Graizer, V. Tilts in strong ground motion. Bull. Seismol. Soc. Am. 2006, 96, 2090–2106. [Google Scholar] [CrossRef]
Wu, Y.; Wu, C. Approximate recovery of coseismic deformation from Taiwan strong-motion records. J. Seismol. 2007, 11, 159–170. [Google Scholar] [CrossRef]
Wang, R.; Schurr, B.; Milkereit, C.; Shao, Z.; Jin, M. An improved automatic scheme for empirical baseline correction of digital strong-motion records. Bull. Seismol. Soc. Am. 2011, 101, 2029–2044. [Google Scholar] [CrossRef]
Bock, Y.; Melgar, D.; Crowell, B.W. Real-time strong-motion broadband displacements from collocated GPS and accelerometers. Bull. Seismol. Soc. Am. 2011, 101, 2904–2925. [Google Scholar] [CrossRef]
Tu, R.; Wang, R.; Ge, M.; Walter, T.R.; Ramatschi, M.; Milkereit, C.; Bindi, D.; Dahm, T. Cost effective monitoring of ground motion related to earthquakes, landslides or volcanic activities by joint use of a single-frequency GPS and a MEMS accelerometer. Geophys. Res. Lett. 2013, 40, 3825–3829. [Google Scholar] [CrossRef]
Tu, R.; Ge, M.; Wang, R.; Walter, T.R. A new algorithm for tight integration of real-time GPS and strong-motion records, demonstrated on simulated, experimental, and real seismic data. J. Seismol. 2014, 18, 151–161. [Google Scholar] [CrossRef]
Tu, R.; Wang, R.; Walter, T.R.; Diao, F.Q. Adaptive recognition and correction of baseline shifts from collocated GPS and accelerometer using two phases Kalman filter. Adv. Space Res. 2014, 54, 1924–1932. [Google Scholar] [CrossRef]
Tu, R.; Zhang, Q.; Wang, L.; Liu, Z.K.; Huang, G.W. An improved method for tight integration of GPS and strong-motion records: Complementary advantages. Adv. Space Res. 2015, 56, 2335–2344. [Google Scholar] [CrossRef]
Tu, R.; Liu, J.H.; Lu, C.X.; Zhang, R.; Zhang, P.F.; Lu, X.C. Cooperating the BDS, GPS, GLONASS and strong-motion observations for real-time deformation monitoring. Geophys. J. Int. 2017, 209, 1408–1417. [Google Scholar] [CrossRef]
Geng, J.; Bock, Y.; Melgar, D.; Crowell, B.W.; Haase, S. A new seismogeodetic approach applied to GPS and accelerometer observations of the 2012 Brawley seismic swarm: Implications for earthquake early warning. Geochem. Geophys. Geosyst. 2013, 14, 2124–2142. [Google Scholar] [CrossRef] [Green Version]
Li, X.; Ge, M.; Zhang, Y.; Wang, R.; Klotz, J.; Wicket, J. High-rate coseismic displacements from tightly-integrated processing of raw GPS and accelerometer data. Geophys. J. Int. 2013, 195, 612–624. [Google Scholar] [CrossRef]
Wang, R.; Parolai, S.; Ge, M.; Jin, M.; Walter, T.R.; Zschau, J. The 2011 Mw 9.0 Tohoku-Oki earthquake: Comparison of GPS and Strong-motion. Bull. Seismol. Soc. Am. 2012. [Google Scholar] [CrossRef]
Geng, J.; Bock, Y.; Melgar, D.; Crowell, B.; Haase, J. A new seismogeodetic approach applied to GPS and accelerometer observations of the 2012 Brawley seismic swarm: Implications for earthquake early warning. Geochem. Geophys. Geosyst. 2013, 14, 2124–2142. [Google Scholar] [CrossRef] [Green Version]
Geng, J.; Melgar, D.; Bock, Y.; Pantoli, E.; Restrepo, J. Recovering coseismic point ground tilts from collocated high-rate GPS and accelerometers. Geophys. Res. Lett. 2013, 40, 5095–5100. [Google Scholar] [CrossRef] [Green Version]
Saunders, J.; Goldberg, E.; Haase, J.; Bock, Y.; Offield, D.; Melgar, D.; Restrepo, J.; Fleischman, R.; Nema, A.; Geng, J.; et al. Seismogeodesy Using GPS and Low-Cost MEMS Accelerometers:Perspectives for Earthquake Early Warning and Rapid Response. Bull. Seismol. Soc. Am. 2016, 106, 2469–2489. [Google Scholar] [CrossRef]
Fleming, K.; Picozzi, M.; Milkereit, C.; Kuehnlenz, F.; Lichtblau, B.; Fischer, J.; Zulfikar, C.; Ozel, O. The SAFER and EDIM Working Groups. The Self-Organising Seismic Early Warning Information System (SOSEWIN). Seismol. Res. Lett. 2009, 80, 755–771. [Google Scholar] [CrossRef]
Msaewe, H.; Hancock, C.; Psimoulis, P.; Roberts, G.; Bonenberg, L. Monitoring Dynamic Deflections at Towers of Severn Suspension Bridge in the UK Using GNSS Technique. In Proceedings of the ISGNSS, Hong Kong, China, 10–13 December 2017. [Google Scholar]
Roberts, G.; Tang, X.; He, X. Accuracy analysis of GPS/BDS relative positioning using zero-baseline measurements. J. Glob. Position. Syst. 2018, 16, 7. [Google Scholar] [CrossRef] [Green Version]
Psimoulis, P.; Pytharouli, S.; Karambalis, D.; Stiros, S. Potential of GPS to measure frequencies of oscillation of engineering structures. J. Sound Vib. 2008, 318, 606–623. [Google Scholar] [CrossRef]
Moschas, F.; Stiros, S. Dynamic Deflections of a Stiff Footbridge Using 100-Hz GNSS and Accelerometer Data. J. Surv. Eng. 2015, 141, 04015003. [Google Scholar] [CrossRef]
Häberling, S.; Rothacher, M.; Zhang, Y.; Clinton, J.; Geiger, A. Assessment of high-rate GPS using a single-axis shake table. J. Geod. 2015, 89, 697–709. [Google Scholar] [CrossRef]
Michel, C.; Kelevitz, K.; Houlié, N.; Edwards, B.; Psimoulis, P.; Su, Z.; Clinton, J.; Giardini, D. The Potential of High-Rate GPS for Strong Ground Motion Assessment. Bull. Seismol. Soc. Am. 2017, 107, 1849–1859. [Google Scholar] [CrossRef]
Psimoulis, P.; Houlié, N.; Meindl, M.; Rothacher, M. Consistency of GPS and strong-motion records: Case study of Mw9.0 Tohoku-Oki 2011 earthquake. Smart Struct. Syst. 2015. [Google Scholar] [CrossRef]

Figure 1. Comparison of average energies for GNSS (bottom) and SM (top) detectors.

Figure 2. Correction of initial baseline shift for GNSS (top) and strong-motion (bottom). The black line represents raw values; the red line, estimated corrections; and the blue line, corrected values.

Figure 3. Quality evaluation of the displacement and velocity after an earthquake. The upper graph shows the displacement, which is stable around a fixed offset; the lower graph shows the velocity, which is stable around 0 m/s.

Figure 4. Real-time integration system of GNSS and SM observations. ‘RT’ indicates real-time, ‘Dis, Vel, and Acc’ represent displacement, velocity and acceleration, respectively, ‘obs’ represents observation, and ‘eph’ represents ephemeris.

Figure 5. Experimental platform (left) and processing facilities (right) of the combined system GPS and SM system.

Figure 6. Movement sequence of the GPS and SM experiments. The graph covers the time for eight tests (T1–T8), each indicated by a dashed red rectangle.

Figure 7. Comparisons of displacement, velocity, and acceleration for the GPS (left); SM (middle); and combined system (COM) (right).

Figure 8. Comparison of the PSD for the displacements of the first experiment by different sensors. Black, blue and red lines represent the GPS, SM, and the combined system, respectively, the (top) represents no motion and (bottom) represents displacement motion.

Figure 9. Displacement (left); velocity (middle); and acceleration (right) output of the integration system for the GPS and SM data. Note that vertical scales vary between the tests.

Figure 10. Experimental equipment used to test the integration of BDS and SM data. Hardware consists of a GNSS sensor, an SM sensor, and a tracked platform.

Figure 11. Movement sequence of eight BDS and SM experiments. Each green dotted rectangle represents one test.

Figure 12. Comparison of displacements (top) and PSD (bottom) for different GNSS systems. Black, blue, red, and pink lines represent the GPS, GLONASS, BDS, and their combinations, respectively.

Figure 13. Displacement (left); velocity (middle); and acceleration (right) output of the integration system for the BDS and SM data. Black curves indicate east–west displacement of the test bed, whereas the red curves indicate north–south motion.

Figure 14. Displacement (left); velocity (middle); and acceleration (right) output of the integration system for the Baja earthquake data. The upper three plots represent pairs of ‘P496-5058’; the lower three plots represent pairs of ‘P744-5028’.

Figure 15. Differences of displacement waveforms between GPS and combined solutions. The (left) represents pairs of ‘P496-5058’; the (right) represents pairs of ‘P744-5028’.

Figure 16. PSD of displacement waveforms by different methods (‘‘COM” represents the GPS and SM combined results, ‘‘DIF” represents the difference between GPS and combined results). The (left) represents pairs of ‘P496-5058’; the (right) represents pairs of ‘P744-5028’.

Table 1. Quality statistics results for the GPS and SM data. D indicates the displacement (m) and V indicates velocity (m/s).

Test	D-Bias	D-STD	V-Bias	V-STD	Test	D-Bias	D-STD	V-Bias	V-STD
T1	−0.006	0.002	1.3 × 10⁻⁴	0.001	T5	0.012	0.003	1.2 × 10⁻⁴	0.001
T2	0.011	0.003	1.1 × 10⁻⁴	0.001	T6	0.008	0.002	−1.0 × 10⁻⁴	0.001
T3	0.013	0.003	−1.1 × 10⁻⁴	0.001	T7	0.016	0.003	1.4 × 10⁻⁴	0.001
T4	−0.014	0.002	1.2 × 10⁻⁴	0.001	T8	−0.013	0.003	1.1 × 10⁻⁴	0.001

Table 2. Quality statistics for the BDS and SM data (D indicates displacement (m), V indicates velocity (m/s), and N and E represent the north and east components, respectively).

Test	D-Bias	D-STD	V-Bias	V-STD	Test	D-Bias	D-STD	V-Bias	V-STD
E1(N)	−0.013	0.002	1.3 × 10⁻⁴	0.001	E4(E)	−0.014	0.002	0.6 × 10⁻⁴	0.001
E2(N)	0.014	0.002	1.0 × 10⁻⁴	0.001	E5(E)	0.012	0.003	1.1 × 10⁻⁴	0.001
E3(N)	0.013	0.002	0.9 × 10⁻⁴	0.001	E6(E)	0.013	0.002	2.0 × 10⁻⁴	0.001
E7(E)	0.016	0.003	1.0 × 10⁻⁴	0.001	E7(N)	0.015	0.002	1.1 × 10⁻⁴	0.001
E8(E)	0.013	0.003	−0.9 × 10⁻⁴	0.001	E8(N)	0.017	0.003	−1.4 × 10⁻⁴	0.001

Table 3. Quality statistics for Baja earthquake data. D indicates displacement (m), V indicates velocity (m/s).

Pairs	Component	D-Bias	D-STD	V-Bias	V-STD
P496-5058	East	0.026	0.003	1.4 × 10⁻⁴	0.002
	North	0.011	0.003	−1.2 × 10⁻⁴	0.002
	Up	0.044	0.005	1.0 × 10⁻⁴	0.001
P744-5028	East	0.021	0.002	1.3 × 10⁻⁴	0.002
	North	0.018	0.003	1.0 × 10⁻⁴	0.002
	Up	0.036	0.004	−0.9 × 10⁻⁴	0.001

© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Tu, R.; Zhang, R.; Zhang, P.; Liu, J.; Lu, X. Integration of Single-Frequency GNSS and Strong-Motion Observations for Real-Time Earthquake Monitoring. Remote Sens. 2018, 10, 886. https://doi.org/10.3390/rs10060886

AMA Style

Tu R, Zhang R, Zhang P, Liu J, Lu X. Integration of Single-Frequency GNSS and Strong-Motion Observations for Real-Time Earthquake Monitoring. Remote Sensing. 2018; 10(6):886. https://doi.org/10.3390/rs10060886

Chicago/Turabian Style

Tu, Rui, Rui Zhang, Pengfei Zhang, Jinhai Liu, and Xiaochun Lu. 2018. "Integration of Single-Frequency GNSS and Strong-Motion Observations for Real-Time Earthquake Monitoring" Remote Sensing 10, no. 6: 886. https://doi.org/10.3390/rs10060886

APA Style

Tu, R., Zhang, R., Zhang, P., Liu, J., & Lu, X. (2018). Integration of Single-Frequency GNSS and Strong-Motion Observations for Real-Time Earthquake Monitoring. Remote Sensing, 10(6), 886. https://doi.org/10.3390/rs10060886

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Integration of Single-Frequency GNSS and Strong-Motion Observations for Real-Time Earthquake Monitoring

Abstract

1. Introduction

2. Methodology

2.1. Determination of the Signal Window

2.2. Initial Baseline Shift Correction of the Integrated System before an Earthquake

2.3. Integration during an Earthquake

2.4. Quality Assessment of the Integrated System

2.5. Implementation of the Real-Time Integrated System

3. Validation Results and Analysis

3.1. GPS and SM Data Test

3.2. BDS and SM Data Test

3.3. Baja Earthquake Data Test

4. Discussion

5. Conclusions

Author Contributions

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI