Next Article in Journal
Tire-Pressure Identification Using Intelligent Tire with Three-Axis Accelerometer
Next Article in Special Issue
Observing the Ocean Submesoscale with Enhanced-Color GOES-ABI Visible Band Data
Previous Article in Journal
Electrocardiogram Classification Based on Faster Regions with Convolutional Neural Network
Previous Article in Special Issue
Rapid Invasion of Spartina Alterniflora in the Coastal Zone of Mainland China: Spatiotemporal Patterns and Human Prevention
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Accurate Evaluation of Vertical Tidal Displacement Determined by GPS Kinematic Precise Point Positioning: A Case Study of Hong Kong

School of Geosciences and Info-physics, Central South University, Changsha 410083, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(11), 2559; https://doi.org/10.3390/s19112559
Submission received: 8 May 2019 / Revised: 29 May 2019 / Accepted: 30 May 2019 / Published: 5 June 2019
(This article belongs to the Special Issue Remote Sensing Applications in Coastal Areas)

Abstract

:
Global Positioning System (GPS) kinematic precise point positioning (KPPP) is an effective approach for estimating the Earth’s tidal deformation. The accuracy of KPPP is usually evaluated by comparing results with tidal models. However, because of the uncertainties of the tidal models, the accuracy of KPPP-estimated tidal displacement is difficult to accurately determine. In this paper, systematic vector differences between GPS estimates and tidal models were estimated by least squares methods in complex domain to analyze the uncertainties of tidal models and determine the accuracy of KPPP-estimated tidal displacements. Through the use of GPS data for 12 GPS reference stations in Hong Kong from 2008 to 2017, vertical ocean tide loading displacements (after removing the body tide effect) for eight semidiurnal and diurnal tidal constituents were obtained by GPS KPPP. By an in-depth analysis of the systematic and residual differences between the GPS estimates and nine tidal models, we demonstrate that the uncertainty of the tidal displacement determined by GPS KPPP for the M2, N2, O1, and Q1 tidal constituents is 0.2 mm, and for the S2 constituent it is 0.5 mm, while the accuracy of the GPS-estimated K1, P1, and K2 tidal constituents is weak because these three tidal constituents are affected by significant common-mode errors. These results suggest that GPS KPPP can be used to precisely constrain the Earth’s vertical tidal displacement in the M2, N2, O1, and Q1 tidal frequencies.

1. Introduction

The tidal force of the Sun and Moon can cause deformation of the solid Earth, which is named “body tide”, and it also causes periodic rise and fall of the ocean mass, which is called “ocean tide”. The response of the solid Earth to this mass distribution is known as “ocean tide loading” (OTL). Modern high-precision geodetic technology can determine large-scale deformation at the sub-millimeter level, but the effect of tidal displacement must be considered [1,2,3,4]. Body tide can now be modeled with an accuracy of usually less than 1% [5]. OTL displacement can be obtained by convolution integral of an ocean tide model and the Green’s function calculated from an Earth model [6], which we call an OTL model. Penna et al. [7] reported that the height error of M2 OTL displacement can reach around 20% (~8 mm), depending on the ocean tide model adopted and the handling strategy for the grid cells in coastal areas. The change of the vertical OTL displacement in some areas when using various Earth models can reach 2–3 mm [8]. Therefore, OTL is still the main uncertain factor in tidal displacement analysis.
Since 2000, the Global Positioning System (GPS) has been widely used to estimate OTL displacement, including eight semidiurnal and diurnal tidal constituents. The methods for estimating OTL displacement can be divided into static and kinematic methods. Adding OTL displacements (amplitude and phase) to the static (usually one-day) GPS solution model and estimating them together with other parameters such as coordinates is called a static method. Schenewerk et al. [9] first determined the vertical OTL displacements by a static method, and reported that the amplitude differences of the M2 height component between the GPS estimates and the OTL model were less than 5 mm at 90% of the GPS reference stations. Over the next 10 years, many scholars [10,11,12,13] used static methods to estimate OTL displacements in different areas, for which the M2 vector difference between the GPS estimates and the OTL model can reach the sub-millimeter level. Yuan et al. [14] estimated global OTL displacements, which indicated that the vertical root-mean-square (RMS) residuals of M2 were 0.6–1.2 mm between the GPS estimates and the different OTL models.
Estimating OTL displacement from a high-sampling GPS coordinate series is called the kinematic method. Khan and Tscherning [15], Khan and Scherneck [16], Yun et al. [17], Melachroinos et al. [18], and Vergnolle et al. [19] used a GPS (1–4 h) relative coordinate series to estimate relative OTL displacements in different regions, for which the M2 vector difference between the GPS estimates and OTL models was usually 1–5 mm. King et al. [20] were the first to estimate absolute OTL displacements in Antarctica using a GPS kinematic precise point positioning (KPPP) coordinate series every 300 s. Since then, Ito et al. [21], Ito and Simons [22], Bos et al. [8], and Tu et al. [23] have used KPPP to estimate the tidal deformation in different areas. These authors reported that the vertical M2 vector difference between the GPS estimates and OTL models can reach the sub-millimeter level. Zhao et al. [24] used an 18.6-year KPPP series to estimate OTL displacements at Shanghai Astronomical Observatory (SHAO) station in China, and reported an M2 RMS agreement of around 0.6 mm in the vertical component compared with OTL models.
In previous studies, the accuracy of GPS-estimated OTL displacement has been evaluated by comparing the results with OTL models. However, this is not accurate enough, because the difference between GPS estimates and OTL models also includes the uncertainties of the OTL models, the unmodeled body tide, the numerical calculation error, and so on. Penna et al. [25] added synthetic signals with an amplitude of 0–6 mm and a period of 13.96 h to a 4-year KPPP coordinate series and demonstrated that the added signals could be determined with an uncertainty of 0.2 mm. However, whether the accuracy of the KPPP estimates of the eight semidiurnal and diurnal tidal constituents can reach 0.2 mm needs to be further validated.
In a small region, the unmodeled body tide and the uncertainties of the OTL model have almost the same effect on all the GPS reference stations, which can be modeled as a systematic component for all the stations. In this study, the systematic and residual components at all the stations were obtained by least squares methods, and used to evaluate the accuracy of the eight KPPP-estimated semidiurnal and diurnal tidal displacements. The Hong Kong area has many uniformly distributed GPS reference stations (Figure 1), so it was chosen as the research area in this study. Through the GPS-observed data from 2008 to 2017, an in-depth study was conducted that how precisely GPS KPPP can provide constraints for the eight semidiurnal and diurnal tidal constituents.

2. Data Processing and Analysis Method

2.1. Method of KPPP-estimated OTL Displacement

The influence of OTL in the horizontal direction is about one-third of that in the vertical direction [8]. Therefore, in this paper, we mainly analyzed the OTL displacement in the vertical direction. According to the International Earth Rotation and Reference Systems Service (IERS) Convention 2010 [26], the vertical OTL displacement of a station can be represented by four semidiurnal (M2, S2, N2 and K2) and four diurnal (K1, O1, P1 and Q1) tidal constituents, depending on the frequencies:
Δ c = j = 1 8 f j A c j cos ( ω j t + χ j ( t 0 ) + u j ϕ c j )
where A c j and ϕ c j respectively represent the amplitude and phase of the j-th tidal constituent at the station; f j and u j denote the nodal modulation corrections; ω j is the angular frequency; and χ j ( t 0 ) is the astronomical argument at reference time t 0 . Table 1 lists the periods of the eight semidiurnal and diurnal tidal constituents. For a given location, except that the amplitude and phase of the j-th tidal constituent need to determined, the remaining parameters can be obtained according to the expression of these arguments in the IERS Convention 2010 [26].
GPS kinematic precise point positioning (KPPP) is an absolute positioning method, which can directly determine the position of a receiver based on satellite constellations [27]. Therefore, an absolute displacement series in vertical direction with a high sampling rate can be obtained by the GPS KPPP solution. OTL displacements can usually be effectively estimated by using a 4-years of GPS observation [25], but simultaneous nodal correction is also needed. The OTL amplitude and phase of the semidiurnal and diurnal tidal constituents can be estimated from the KPPP coordinate series by using the classical harmonic analysis method [28].

2.2. GPS Data Acquisition and Processing

Hong Kong is located on the north coast of the South China Sea, and is mainly composed of the Kowloon Peninsula and many islands. The Hong Kong coastline is complex and is affected by OTL. GPS observation data for 2008–2017 provided by the Surveying and Mapping Office of the Hong Kong Lands Department (https://www.geodetic.gov.hk/smo/gsi/programs/tc/index.htm) were used, in which the HKFN station was replaced by the T430 station after 2014. The distribution of the GPS reference stations in the Hong Kong area are shown in Figure 1.
By using 24-h GPS observation data in RINEX format, the kinematic PPP solution for every 600 s was carried out with Bernese GNSS software version 5.2 [29] (http://www.bernese.unibe.ch/docs/DOCU52.pdf), in which the elevation cutoff angle was set to 3°. The main parameter files, i.e., satellite orbit, 30 s satellite clock, Earth rotation, ionosphere parameter, and differential code bias (DCB), were provided by the European Center for Orbit Determination in Europe (COD). The tropospheric delay was modeled using the Global Mapping Function (GMF) [30]. The antenna phase center file used PCV_COD.I08. Body tide and pole tide were modeled according to the IERS Convention 2010 [26]. No corrections were applied to the ocean and atmospheric tide loading.
The estimated KPPP coordinates were in the IGS08 reference frame. The origin of IGS08 is defined as the center of mass (CM), including the solid Earth, the ocean, and the atmosphere, which is usually called the CM frame. Coordinates under IGS08 reference frame were transformed into the station coordinate system, including east, north, and vertical directions. The linear trend of the vertical series was removed. A threshold of 200 mm was set to eliminate gross errors. The % of Gross Errors are shown in Table 2. The eight semidiurnal and diurnal tidal components in the vertical direction at 12 stations were estimated using the method described in Section 2.1. The amplitudes and phases of OTL displacements for eight tidal constituents determined by GPS KPPP are listed in Table A1 and Table A2.

2.3. OTL Model and Tidal Signal Analysis

OTL displacement predicted by a model is usually based on the Green’s function [6], which is determined by the adopted Earth model and represents the response characteristics of the solid Earth to mass loading. The OTL displacement u for a location r is as follows:
u ( r ) = Ω ρ G ( | r r | ) Z ( r ) d Ω
where Z is a vector expression of the tide at r provided by the ocean tide model; Ω represents the global ocean area; ρ is the density of sea water; and G is the Green’s function, characterizing the load deformation caused by the unit mass of sea water.
In this paper, eight recent ocean tide models were considered: FES2014b, GOT4.10c, TPXO.8, DTU10, EOT11a, HAMTIDE, NAO99b, and OSU12 [31,32]. The OTL displacements predicted by the FES2014b and GOT4.10c models based on the STW105 Earth model were provided by the International Mass Loading Service [33], for which the density of sea water was set to 1025 kg/m3. The predictions of the OTL displacements by FES2014b, TPXO.8, DTU10, EOT11a, HAMTIDE, NAO99b, and OSU12 were provided by the Ocean Tide Loading Provider [34], for which the density of sea water was set to 1030 kg/m3. The FES2014b and TPXO.8 tide models based on the Preliminary Reference Earth Model (PREM) were calculated with CARGA [7]. The DTU10, EOT11a, HAMTIDE, NAO99b, and OSU12 tide models based on the Gutenberg-Bullen A Earth model were calculated with OLFG/OLMPP [35]. A total of nine OTL models were obtained, and we refer to them as GOT4.10c + STW, FES2014b + STW, FES2014b + PREM, TPXO.8 + PREM, DTU10 + GB, EOT11a + GB, HAMTIDE + GB, NAO99b + GB, and OSU12 + GB, respectively.
After modeling the body tide in the GPS data processing, the GPS-estimated tidal deformation is mainly made up of the OTL displacement. A tidal constituent can be conveniently expressed in terms of the in-phase and out-of-phase amplitudes of the phasors on the complex plane [13]. The phasor vector relationship between the GPS estimate and OTL model is shown in Figure 2.
The phasor difference between the GPS estimate and OTL model for the j-th tidal constituent at station k is as follows:
Z k , j = Z G P S , k , j Z O T L , k , j
where Z G P S , k , j is the tidal displacement estimated by GPS, and Z O T L , k , j denotes the tidal deformation predicted by the OTL model. The complex form of Z k , j can be written as:
Z k , j = A k , j ( cos Φ k , j + i sin Φ k , j )
where i is the imaginary unit; and A k , j and Φ k , j represent the residual amplitude and phase of the j-th tidal constituent, respectively.
The RMS of the phasor difference for the j-th tidal constituent between the GPS estimate and OTL model for N stations is calculated as follows:
RMS j = ( 1 N k = 1 N | Z k , j | 2 ) 1 2
In a small region, the phasor difference for the j-th tidal constituent between the GPS estimate and OTL model at station k can be divided into two kinds of components. The first kind of component is related to the specific station, i.e., the multipath effect on GPS observation. The second kind of component is applicable to all stations, i.e., the unmodeled body tide, the uncertainties of the OTL model caused by ocean tide model errors and mantle anelasticity, the GPS orbit errors, and the signal propagation delay. Therefore, the phasor difference Z k , j of the j-th tidal constituent at station k can be decomposed into Z r e s i d u a l , k , j , which is related to the specific station, and Z s y s t e m , j , which is applicable to all stations. For the N stations, the phasor difference of the j-th tidal constituent can be expressed as follows:
Z k , j = Z s y s t e m , j + Z r e s i d u a l , k , j       ( k = 1 , 2 , , N )
Next, we explain how to obtain Z r e s i d u a l , k , j and Z s y s t e m , j on the complex plane by the least-squares method. Equation (6) is written in matrix form:
[ Z 1 , j Z k , j Z N , j ] = [ 1 1 1 ] Z s y s t e m , j + [ Z r e s i d u a l , 1 , j Z r e s i d u a l , k , j Z r e s i d u a l , N , j ]
Equation (7) can be abbreviated as:
Z j N × 1 = I N × 1 Z s y s t e m , j + Z r e s i d u a l , j N × 1
where I denotes the unit matrix, and Z s y s t e m , j is the parameter to be estimated. It is assumed that Z r e s i d u a l , k , j ( k = 1 , , N ) obey a Gaussian distribution and are uncorrelated with each other. The expected value and variance of Z r e s i d u a l , k , j ( k = 1 , , N ) are E ( | Z r e s i d u a l , k , j | ) = 0 and D ( | Z r e s i d u a l , k , j | ) = σ 2 , respectively. Therefore, the variance-covariance matrix of Z j can be written as:
D ( Z j N × 1 ) = D ( Z r e s i d u a l , j N × 1 ) = σ 2 I N × N
where D ( ) is the variance-covariance matrix. If we let V donate the correction value of Z r e s i d u a l , j and Z ^ s y s t e m , j represent the estimated value of Z s y s t e m , j , then formula (8) can be transformed as:
V N × 1 = I N × 1 Z ^ s y s t e m , j Z j N × 1
From formula (9), we know that Z k , j ( k = 1 , , N ) are equal-weight observations. Therefore, to obtain Z ^ s y s t e m , j based on the least-squares principle, we only need to satisfy V T V = m i n , in which V T is the transpose matrix of V. In order to obtain the minimum value of V T V , the derivative of V T V to Z ^ s y s t e m , j should be 0. The detailed process is as follows:
( V T V ) Z ^ s y s t e m , j = 2 V T V Z ^ s y s t e m , j = 2 V T I = 0
According to Equations (10) and (11), the formula after eliminating V can be obtained as:
I T 1 × N I N × 1 Z ^ s y s t e m , j I T 1 × N Z j N × 1 = 0
Therefore, the least-squares estimate of complex Z s y s t e m , j is:
Z ^ s y s t e m , j = ( I T 1 × N I N × 1 ) 1 I T 1 × N Z j N × 1
Finally, Z r e s i d u a l , k , j ( k = 1 , , N ) can be obtained by formula (6).

3. Results and Discussion

3.1. Comparison with the OTL Models

To analyze the KPPP-estimated tidal constituent stability, according to Equations (6) and (13), we calculated the residual and systematic phasor differences between the GPS estimates (using the observed data for years 4 to 10 since 2008) and the GOT4.10c + STW OTL model, for each of the eight tidal constituents at 12 stations. The corresponding RMSs of the residual and systematic phasor difference of the eight tidal constituents were obtained by formula (5), as shown in Figure 3.
Figure 3a,b indicate that, for the M2, N2, O1 and Q1 tidal constituents, after the fourth year, the RMS of the residual and systematic phasors for the four tidal constituents is stabilized. For the S2 tidal constituent, its residual RMS is stabilized after the fourth year; however, it is not until the ninth year that the systematic RMS of the S2 constituent begins to stabilize. Although the residual and systematic RMS of the P1 constituent is stabilized after the fourth year, the magnitude of the residual and systematic RMS for the P1 tidal constituent is significant. The residual RMS of the K1 and K2 constituents is stabilized after the fifth year, but their systematic RMS does not convergence until the tenth year. Figure 4a shows the spatial distributions of the phasor vector differences, between the GPS estimates (using the observed data for 2008–2017) and the GOT4.10c + STW OTL model, of the eight tidal constituents at 12 stations. Figure 4b shows the systematic phasor vector differences estimated from Figure 4a by formula (13) for the eight tidal constituents. The residual vector phasor differences were obtained by differencing Figure 4a,b, as shown in Figure 4c. The expression of the phasor vector difference for a given tidal constituent is given in Figure 2.
From Figure 4a, we can see that the spatial distributions of the phasor vector differences for the M2, K2, K1, O1, and P1 tidal constituents have obvious coherence. The systematic phasor differences could be well derived, as shown in Figure 4b. Figure 4a also shows that the spatial distributions of the phasor vector differences for the S2, N2, and Q1 tidal constituents had no obvious coherence; therefore, there were no distinguishing systematic differences for these three constituents, as shown in Figure 4b. Figure 4c indicates that the residual phasor differences of the eight tidal constituents had no spatial correlation, i.e., they are related to the location.

3.2. Accuracy Assessment for the KPPP-Estimated OTL

To validate the accuracy of the GPS KPPP estimates of the eight tidal constituents, according to Equations (5), (6) and (13), we calculated the RMS of the total phasor differences, the systematic phasor differences and the corresponding residual phasor differences between the GPS estimates (using the observed data for 2008–2017) and the nine OTL models, for each of the eight tidal constituents at 12 stations. The RMS values of the eight tidal constituents are shown in Figure 5.
Figure 5 shows that the systematic RMS of the N2 and Q1 tidal constituents between the GPS estimates and the different OTL models varies from 0.02 to 0.40 mm and from 0.12 to 0.40 mm, respectively. The systematic RMS varies significantly with the different OTL models, which shows that the systematic of the N2 and Q1 tidal constituents was mainly caused by the uncertainties of the OTL models. The uncertainties of the OTL models result in the systematic differences. After removing the systematic differences, the RMS of the residual phasors of the N2 and Q1 tidal constituents are all around 0.2 mm. The residual RMS does not change with the different OTL models, so the residual RMS reflects the uncertainty of the GPS estimates.
As shown in Figure 5, the systematic RMS change of the M2 and O1 tidal constituents with the various OTL models is from 0.01 to 1.02 mm and 0.20 to 1.01 mm, respectively. These results indicate that the systematic RMS of the M2 and O1 constituents was mainly caused by the uncertainties of the OTL models. After removing the systematic differences, the residual RMS of the O1 tidal constituent was around 0.2 mm, and no longer changed with the various OTL models. The residual RMS values of the M2 tidal constituent between the GPS estimates and the six OTL models of GOT4.10c + STW, FES2014b + STW, FES2014b + PREM, TPXO.8 + PREM, NAO99b + GB, and OSU12 + GB were all 0.2 mm. However, the residual RMS of the M2 tidal constituent between the GPS estimates and the other three OTL models (DTU10 + GB, EOT11a + GB, HAMTIDE + GB) were around 0.35 mm. This indicates that the residual RMS of the M2 tidal constituent changes slightly with the different OTL models. Therefore, the residual RMS should absorb the errors of the three OTL models, which may not exist entirely as systematic differences, because the M2 tidal resolution of the DTU10, EOT11a, and HAMTIDE ocean tide models is 1/8° × 1/8° smaller than the 1/2° × 1/2° size of Hong Kong. Overall the uncertainty of the GPS-estimated M2 tidal constituent is at the 0.2 mm level.
Figure 5 also shows that there was obvious systematic RMS of the K1, P1, and K2 tidal constituents between the GPS estimates and the nine OTL models, of around 6 mm, 2 mm, and 9 mm, respectively. This indicates that the K1, K2, and P1 tidal constituents are affected by the significant common-mode errors for all the stations. These errors may include: (1) the unmodeled body tide; (2) the OTL model uncertainty caused by the ocean tide model and mantle inelasticity; and (3) GPS-related error, such as satellite orbit error, which has almost the same effect for all stations. Firstly, for the three tidal constituents of K1, K2, and P1, the unmodeled body tide error should be at the sub-millimeter level [12]. Secondly, the OTL displacement differences for the K1, P1, and K2 tidal constituents among the predictions by the nine different OTL models are at the sub-millimeter level, which indicates that that the OTL models are unlikely to have caused such misfits. Therefore, the common-mode errors are likely caused by the GPS-related errors. After removing the systematic differences of the K1, P1, and K2 tidal constituents, the residual RMS was around 0.8 mm, 0.5 mm, and 1 mm, respectively. The residual RMS was still significant, which indicates that these three tidal constituents are probably affected by the station-related error. Previous studies [13,14] inferred that K1 and K2 may be affected by the satellite orbit error and multipath effect, because the period of the K1 coincides with that of the satellite constellation, and the period of the K2 is almost the same as the satellite orbital period, while P1 may be disturbed by residual GPS signal propagation errors or non-tidal daily cycles in the atmosphere because its period is very close to solar day.
As shown in Figure 5, the total RMS of the S2 tidal constituent between the GPS estimates and the nine OTL models were between 0.5 to 0.6 mm. There are some not obvious systematic differences of 0.1 to 0.3 mm between the GPS estimates and the nine different OTL models, which may arise from the effect of the atmospheric S2 tide. The effect of the atmospheric S2 tide was analyzed using the GEOSFPIT atmospheric model provided by the International Mass Loading Service [8], with a resolution of 0.625° × 0.5°. The results show that the amplitude of the atmospheric tide displacement, as predicted by the GEOSFPIT model, is only 0.22 mm at Hong Kong. Figure 6 shows that, after correcting the atmospheric tide by the GEOSFPIT model, the RMS residual of the S2 tidal constituent between the GPS estimates and some OTL models slightly reduced. Overall there was still an obvious RMS residual of about 0.5 mm between the GPS estimates and the OTL models. Figure 5 shows that the residual RMS of the S2 tidal constituent were all nearly 0.5 mm. This indicates that there are some station-related errors of around 0.5 mm, which arise from GPS-related error. GPS signal propagation error or non-tidal daily cycles in the atmosphere may bias the estimate of S2 because the period of S2 is half of the solar day [13,14]. This leads to the uncertainty of GPS-estimated S2 tidal constituent being at 0.5 mm level.

4. Conclusions

In this study, we used GPS kinematic PPP to estimate the vertical OTL displacements of eight semidiurnal and diurnal tidal constituents with 10-year GPS data from 12 GPS reference stations in Hong Kong. The systematic and residual phasor differences between the GPS estimates and nine OTL models were obtained by least squares methods in the complex domain. Our experimental results suggest that:
For the M2, N2, O1, and Q1 tidal constituents, the RMS residual between the GPS KPPP estimates and the different OTL models is mainly due to the uncertainties of the OTL models, which exist as a systematic component. When the systematic differences are removed, the residual RMS is basically within 0.2 mm, which indicates that the uncertainty of the four KPPP-estimated tidal constituents is at the 0.2 mm level. The K1, P1, and K2 tidal constituents are affected by significant GPS-related common-mode error. After removing the systematic differences, the residual RMS was still significant, which indicates that the three tidal constituents are also affected by GPS-related residual error. The RMS residual of the S2 tidal constituent between the GPS estimates and the OTL models is mainly caused by GPS-related error. Overall the uncertainty of the GPS-estimated S2 tidal constituent is at the 0.5 mm level. These results indicate that GPS KPPP can provide constraints for the Earth’s vertical tidal displacement at the 0.2 mm level for the M2, N2, O1, and Q1 tidal components, and the 0.5 mm level for the S2 tidal component.

Author Contributions

G.W. and Q.W. conceived and designed the experiments; G.W. and W.P. processed data; G.W. and Q.W. carried out the experimental analysis.

Funding

This research was funded by the National Natural Science Foundation of China, grant number U1531128 and the Fundamental Research Funds for the Central Universities of Central South University, grant number 2018zzts680.

Acknowledgments

We thank the Lands Department of Hong Kong for providing GPS observations. We are grateful to the providers of the Bernese software, the International Mass Loading Service and Ocean Tide Loading Provider. We thank the two anonymous reviewers for their careful reading and thorough comments.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The amplitudes and phases of vertical ocean tide loading (OTL) displacements for eight semidiurnal and diurnal tidal constituents determined by GPS kinematic precise point positioning (KPPP) were shown in Table A1 and Table A2, in which all phases are positive for the Greenwich lag.
Table A1. KPPP-estimated vertical OTL displacements of four semidiurnal tides.
Table A1. KPPP-estimated vertical OTL displacements of four semidiurnal tides.
M2S2N2K2
StationA (mm) Φ (°)A (mm) Φ (°)A (mm) Φ (°)A (mm) Φ (°)
HKFN4.86189.21.47181.41.04169.48.1414.7
HKKT5.25187.82.03200.01.10175.07.3230.6
HKLT5.41189.71.70193.21.39181.17.7923.1
HKMW5.88192.91.41193.01.44182.78.0414.7
HKNP6.13194.52.34195.61.27192.78.3121.3
HKOH6.17189.82.09207.31.50181.17.8031.9
HKPC5.73191.11.56210.21.54184.97.9515.2
HKSC5.82189.91.94189.71.42182.511.5518.4
HKSL5.80193.62.08230.11.40184.69.0017.5
HKSS5.55184.91.79190.71.33179.97.6931.7
HKST5.44184.82.02178.31.34176.98.0920.1
HKWS5.51185.52.08186.71.22179.98.9031.5
Table A2. KPPP-estimated vertical OTL displacements of four diurnal tides.
Table A2. KPPP-estimated vertical OTL displacements of four diurnal tides.
StationK1O1P1Q1
A (mm) Φ (°)A (mm) Φ (°)A (mm) Φ (°)A (mm) Φ (°)
HKFN7.2132.76.73306.91.85260.21.20297.4
HKKT7.9535.86.96305.51.98262.41.64281.6
HKLT6.4740.67.25304.41.00255.81.62283.9
HKMW6.8539.67.77306.61.22297.51.93284.1
HKNP8.2132.57.76308.81.17260.81.73291.4
HKOH8.8432.97.86306.91.71296.51.66282.2
HKPC8.0332.17.46306.01.38300.21.79281.4
HKSC8.2037.77.50309.31.07311.21.70284.2
HKSL8.0841.77.63308.11.74289.51.53286.3
HKSS7.4627.67.07306.70.96254.11.47284.1
HKST7.2226.87.27307.61.60291.31.37280.9
HKWS7.2032.57.30305.31.50268.31.35276.5

References

  1. Yuan, L.G.; Ding, X.L.; Zhong, P.; Chen, W.; Huang, D.F. Estimates of ocean tide loading displacements and its impact on position time series in Hong Kong using a dense continuous GPS network. J. Geod. 2009, 83, 999–1015. [Google Scholar] [CrossRef]
  2. Li, Z.; Jiang, W.; Ding, W.; Deng, L.; Peng, L. Estimates of minor ocean tide loading displacement and its impact on continuous GPS coordinate time series. Sensors 2014, 14, 5552–5572. [Google Scholar] [CrossRef] [PubMed]
  3. Lei, M.; Wang, Q.; Liu, X.; Xu, B.; Zhang, H. Influence of ocean tidal loading on InSAR offshore areas deformation monitoring. Geod. Geodyn. 2017, 8, 70–76. [Google Scholar] [CrossRef] [Green Version]
  4. Peng, W.; Wang, Q.; Cao, Y. Analysis of Ocean Tide Loading in Differential InSAR Measurements. Remote Sens. 2017, 9, 101. [Google Scholar] [CrossRef]
  5. Krasna, H.; Bohm, J.; Schuh, H. Tidal Love and Shida numbers estimated by geodetic VLBI. J. Geodyn. 2013, 70, 21–27. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Farrell, W.E. Deformation of the Earth by surface loads. Rev. Geophys. 1972, 10, 761–797. [Google Scholar] [CrossRef]
  7. Penna, N.T.; Bos, M.S.; Baker, T.F.; Scherneck, H.G. Assessing the accuracy of predicted ocean tide loading displacement values. J. Geod. 2008, 82, 893–907. [Google Scholar] [CrossRef] [Green Version]
  8. Bos, M.S.; Penna, N.T.; Baker, T.F.; Clarke, P.J. Ocean tide loading displacements in western Europe: 2. GPS-observed anelastic dispersion in the asthenosphere. J. Geophys. Res. Solid Earth 2015, 120, 6540–6557. [Google Scholar] [CrossRef]
  9. Schenewerk, M.S.; Marshall, J.; Dillinger, W. Vertical ocean-loading deformations derived from a global GPS network. J. Geod. Soc. Jpn. 2001, 47, 237–242. [Google Scholar]
  10. Allinson, C.R.; Clarke, P.J.; Edwards, S.J.; King, M.A. Stability of direct GPS estimates of ocean tide loading. Geophys. Res. Lett. 2004, 31, L15603. [Google Scholar] [CrossRef]
  11. King, M.A.; Penna, N.T.; Clarke, P.J.; King, E.C. Validation of ocean tide models around Antarctica using onshore GPS and gravity data. J. Geophys. Res. 2005, 110, B08401. [Google Scholar] [CrossRef]
  12. Thomas, I.D.; King, M.A.; Clarke, P.J. A Comparison of GPS, VLBI and Model Estimates of Ocean Tide Loading Displacements. J. Geod. 2006, 81, 359–368. [Google Scholar] [CrossRef]
  13. Yuan, L.G.; Chao, B.F. Analysis of tidal signals in surface displacement measured by a dense continuous GPS array. Earth Planet Sci. Lett. 2012, 355–356, 255–261. [Google Scholar] [CrossRef]
  14. Yuan, L.G.; Chao, B.F.; Ding, X.L.; Zhong, P. The tidal displacement field at Earth’s surface determined using global GPS observations. J. Geophys. Res. Solid Earth 2013, 118, 2618–2632. [Google Scholar] [CrossRef]
  15. Khan, S.A.; Tscherning, C.C. Determination of semi-diurnal ocean tide loading constituents using GPS in Alaska. Geophys. Res. Lett. 2001, 28, 2249–2252. [Google Scholar] [CrossRef]
  16. Khan, S.A.; Scherneck, H.-G. The M2 ocean tide loading wave in Alaska: Vertical and horizontal displacements, modelled and observed. J. Geod. 2003, 77, 117–127. [Google Scholar] [CrossRef]
  17. Yun, H.-S.; Lee, D.-H.; Song, D.-S. Determination of vertical displacements over the coastal area of Korea due to the ocean tide loading using GPS observations. J. Geodyn. 2007, 43, 528–541. [Google Scholar] [CrossRef]
  18. Melachroinos, S.A.; Biancale, R.; Llubes, M.; Perosanz, F.; Lyard, F.; Vergnolle, M.; Bouin, M.N.; Masson, F.; Nicolas, J.; Morel, L.; et al. Ocean tide loading (OTL) displacements from global and local grids: Comparisons to GPS estimates over the shelf of Brittany, France. J. Geod. 2007, 82, 357–371. [Google Scholar] [CrossRef]
  19. Vergnolle, M.; Bouin, M.N.; Morel, L.; Masson, F.; Durand, S.; Nicolas, J.; Melachroinos, S.A. GPS estimates of ocean tide loading in NW-France: Determination of ocean tide loading constituents and comparison with a recent ocean tide model. Geophys. J. Int. 2008, 173, 444–458. [Google Scholar] [CrossRef]
  20. King, M. Kinematic and static GPS techniques for estimating tidal displacements with application to Antarctica. J. Geodyn. 2006, 41, 77–86. [Google Scholar] [CrossRef]
  21. Ito, T.; Okubo, M.; Sagiya, T. High resolution mapping of Earth tide response based on GPS data in Japan. J. Geodyn. 2009, 48, 253–259. [Google Scholar] [CrossRef]
  22. Ito, T.; Simons, M. Probing Asthenospheric Density, Temperature, and Elastic Moduli Below the Western United States. Science 2011, 332, 947–951. [Google Scholar] [CrossRef] [PubMed]
  23. Tu, R.; Zhao, H.; Zhang, P.; Liu, J.; Lu, X. Improved Method for Estimating the Ocean Tide Loading Displacement Parameters by GNSS Precise Point Positioning and Harmonic Analysis. J. Surv. Eng. 2017, 143, 04017005. [Google Scholar] [CrossRef]
  24. Zhao, H.; Zhang, Q.; Tu, R.; Liu, Z. Analysis of ocean tide loading displacements by GPS kinematic precise point positioning: A case study at the China coastal site SHAO. Surv. Rev. 2017, 51, 172–182. [Google Scholar] [CrossRef]
  25. Penna, N.T.; Clarke, P.J.; Bos, M.S.; Baker, T.F. Ocean tide loading displacements in western Europe: 1. Validation of kinematic GPS estimates. J. Geophys. Res. Solid Earth 2015, 120, 6523–6539. [Google Scholar] [CrossRef]
  26. Petit, G.; Luzum, B. IERS Conventions (2010). In IERS Technical Note No.36; Verlag des Bundesamts für Kartographie und Geodäsie: Frankfurt am Main, Germany, 2010. [Google Scholar]
  27. Geng, J.; Teferle, F.N.; Meng, X.; Dodson, A.H. Kinematic precise point positioning at remote marine platforms. GPS Solut. 2010, 14, 343–350. [Google Scholar] [CrossRef]
  28. Pawlowicz, R.; Beardsley, B.; Lentz, S. Classical tidal harmonic analysis including error estimates in MATLAB using T_TIDE. Comput. Geosci. 2002, 28, 929–937. [Google Scholar] [CrossRef]
  29. Bernese GNSS Software Version 5.2. Available online: http://www.bernese.unibe.ch/ (accessed on 20 October 2018).
  30. Boehm, J.; Niell, A.; Tregoning, P.; Schuh, H. Global Mapping Function (GMF): A new empirical mapping function based on numerical weather model data. Geophys. Res. Lett. 2006, 33, L07304. [Google Scholar] [CrossRef]
  31. Ray, R.D. Precise comparisons of bottom-pressure and altimetric ocean tides. J. Geophys. Res. Oceans 2013, 118, 4570–4584. [Google Scholar] [CrossRef] [Green Version]
  32. Stammer, D.; Ray, R.D.; Andersen, O.B.; Arbic, B.K.; Bosch, W.; Carrère, L.; Cheng, Y.; Chinn, D.S.; Dushaw, B.D.; Egbert, G.D.; et al. Accuracy assessment of global barotropic ocean tide models. Rev. Geophys. 2014, 52, 243–282. [Google Scholar] [CrossRef] [Green Version]
  33. International Mass Loading Service: Tidal Ocean Loading. Available online: http://massloading.net/toc/ (accessed on 20 October 2018).
  34. Ocean Tide Loading Provider. Available online: http://holt.oso.chalmers.se/loading/ (accessed on 20 October 2018).
  35. Scherneck, H.-G. A parametrized solid earth tide model and ocean tide loading effects for global geodetic baseline measurements. Geophys. J. Int. 1991, 106, 677–694. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Distribution of the GPS reference stations in the Hong Kong area.
Figure 1. Distribution of the GPS reference stations in the Hong Kong area.
Sensors 19 02559 g001
Figure 2. Phasor vector diagram of the relationship between the GPS estimate and OTL model.
Figure 2. Phasor vector diagram of the relationship between the GPS estimate and OTL model.
Sensors 19 02559 g002
Figure 3. GPS-estimated tidal stability for the eight tidal constituents: (a) the RMS of the systematic phasor difference; (b) the RMS of the residual phasor difference corresponding to (a).
Figure 3. GPS-estimated tidal stability for the eight tidal constituents: (a) the RMS of the systematic phasor difference; (b) the RMS of the residual phasor difference corresponding to (a).
Sensors 19 02559 g003
Figure 4. Spatial distribution of the phasor vector differences between the GPS estimates and the GOT4.10c + STW OTL model: (a) total phasor vector differences (black); (b) systematic phasor vector differences (red) estimated from (a); (c) residual phasor vector differences (blue) by (a,b).
Figure 4. Spatial distribution of the phasor vector differences between the GPS estimates and the GOT4.10c + STW OTL model: (a) total phasor vector differences (black); (b) systematic phasor vector differences (red) estimated from (a); (c) residual phasor vector differences (blue) by (a,b).
Sensors 19 02559 g004
Figure 5. RMS values of the eight tidal constituents between the GPS estimates and the nine OTL models.
Figure 5. RMS values of the eight tidal constituents between the GPS estimates and the nine OTL models.
Sensors 19 02559 g005
Figure 6. RMS residual of the S2 tidal constituent between the GPS estimates and the OTL models
Figure 6. RMS residual of the S2 tidal constituent between the GPS estimates and the OTL models
Sensors 19 02559 g006
Table 1. The periods of the eight semidiurnal and diurnal tidal constituents.
Table 1. The periods of the eight semidiurnal and diurnal tidal constituents.
NameHarmonic ConstituentPeriod/h
M2Principal lunar semi-diurnal12.4206
S2Principal solar semi-diurnal12.0000
N2Larger lunar elliptic semi-diurnal12.6583
K2Luni-solar declinational semi-diurnal11.9672
K1Luni-solar declinational diurnal23.9345
O1Principal lunar declinational diurnal25.8193
P1Principal solar declinational diurnal24.0659
Q1Lunar elliptic diurnal26.8680
Table 2. WGS84 coordinates and the Gross Error elimination rate for the 12 stations.
Table 2. WGS84 coordinates and the Gross Error elimination rate for the 12 stations.
StationLongitude (E)Latitude (N)Ellipsoidal Height/mData Length/d% of Gross Errors
HKFN114°08′17″22°29′41″41.21236281.57%
HKKT114°04′00″22°26′41″34.57636401.06%
HKLT113°59′48″22°25′05″125.92236501.23%
HKMW114°00′11″22°15′21″194.94635294.58%
HKNP113°53′38″22°14′57″350.67236431.15%
HKOH114°13′43″22°14′52″166.40136461.08%
HKPC114°02′16″22°17′05″18.13036381.23%
HKSC114°08′28″22°19′20″20.23935942.61%
HKSL113°55′41″22°22′19″95.29736520.92%
HKSS114°16′09″22°25′52″38.71336451.11%
HKST114°11′03″22°23′43″258.70436510.93%
HKWS114°20′07″22°26′03″63.79136520.89%

Share and Cite

MDPI and ACS Style

Wei, G.; Wang, Q.; Peng, W. Accurate Evaluation of Vertical Tidal Displacement Determined by GPS Kinematic Precise Point Positioning: A Case Study of Hong Kong. Sensors 2019, 19, 2559. https://doi.org/10.3390/s19112559

AMA Style

Wei G, Wang Q, Peng W. Accurate Evaluation of Vertical Tidal Displacement Determined by GPS Kinematic Precise Point Positioning: A Case Study of Hong Kong. Sensors. 2019; 19(11):2559. https://doi.org/10.3390/s19112559

Chicago/Turabian Style

Wei, Guoguang, Qijie Wang, and Wei Peng. 2019. "Accurate Evaluation of Vertical Tidal Displacement Determined by GPS Kinematic Precise Point Positioning: A Case Study of Hong Kong" Sensors 19, no. 11: 2559. https://doi.org/10.3390/s19112559

APA Style

Wei, G., Wang, Q., & Peng, W. (2019). Accurate Evaluation of Vertical Tidal Displacement Determined by GPS Kinematic Precise Point Positioning: A Case Study of Hong Kong. Sensors, 19(11), 2559. https://doi.org/10.3390/s19112559

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

Article Metrics

Back to TopTop