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

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.


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.

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: 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.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.

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: 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: where s(t) is the integrated displacement time series, A, B, and C are the polynomial fitting coefficients, and v (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: Similarly, the acceleration integrated to velocity or displacement used least-squares fitting to solve the acceleration initial shift correction: where a (t) is the corrected acceleration and 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.

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: ( 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: 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: Similarly, the acceleration integrated to velocity or displacement used least-squares fitting to solve the acceleration initial shift correction: where a c (t) is the corrected acceleration and â(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.

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]: 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 ) and can be ignored.An accurate acceleration can be obtained from: 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:

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]: 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 c1 , a c2 , a c3 . . .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: 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: 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: and the filter estimation can be expressed as: 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.

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.

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.

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.

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.

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.

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.

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.

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 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 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 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 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.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.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.

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.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.

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.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.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.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.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.

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.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.

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 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.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 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.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.

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,

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.

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.

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

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

Figure 2 .
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 2 .
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 .
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 3 .
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.

Remote 17 Figure 3 .
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 5 .
Figure 5. Experimental platform (left) and processing facilities (right) of the combined system GPS and SM system.

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

Figure 6 .
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 6 .
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.

17 Figure 6 .
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 8 .
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 .
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 8 .
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.

17 Figure 8 .
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 .
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 9 .
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 .
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 10 .
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 .
Figure 11.Movement sequence of eight BDS and SM experiments.Each green dotted rectangle represents one test.

Figure 12 .
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.

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

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

Figure 12 .
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 12 .
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 .
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 13 .
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 .
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 14 .
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'.

17 Figure 14 .
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 16 .
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).

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

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).

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).

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

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