Next Article in Journal
Investigating the Sensitivity of Low-Cost Sensors in Measuring Particle Number Concentrations across Diverse Atmospheric Conditions in Greece and Spain
Next Article in Special Issue
Research on Soil Moisture Estimation of Multiple-Track-GNSS Dual-Frequency Combination Observations Considering the Detection and Correction of Phase Outliers
Previous Article in Journal
Classification of Low Earth Orbit (LEO) Resident Space Objects’ (RSO) Light Curves Using a Support Vector Machine (SVM) and Long Short-Term Memory (LSTM)
Previous Article in Special Issue
Multilayer Model in Soil Moisture Content Retrieval Using GNSS Interferometric Reflectometry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sea-Level Estimation from GNSS-IR under Loose Constraints Based on Local Mean Decomposition

College of Geomatics and Geoinformation, Guilin University of Technology, Guilin 541004, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(14), 6540; https://doi.org/10.3390/s23146540
Submission received: 15 May 2023 / Revised: 14 July 2023 / Accepted: 18 July 2023 / Published: 20 July 2023
(This article belongs to the Special Issue GNSS Sensing and Imaging Based on Monitoring Applications)

Abstract

:
The global navigation satellite system–interferometric reflectometry (GNSS-IR) technique has emerged as an effective coastal sea-level monitoring solution. However, the accuracy and stability of GNSS-IR sea-level estimation based on quadratic fitting are limited by the retrieval range of reflector height (RH range) and satellite-elevation range, reducing the flexibility of this technology. This study introduces a new GNSS-IR sea-level estimation model that combines local mean decomposition (LMD) and Lomb–Scargle periodogram (LSP). LMD can decompose the signal-to-noise ratio (SNR) arc into a series of signal components with different frequencies. The signal components containing information from the sea surface are selected to construct the oscillation term, and its frequency is extracted by LSP. To this end, observational data from SC02 sites in the United States are used to evaluate the accuracy level of the model. Then, the performance of LMD and the influence of noise on retrieval results are analyzed from two aspects: RH ranges and satellite-elevation ranges. Finally, the sea-level variation for one consecutive year is estimated to verify the stability of the model in long-term monitoring. The results show that the oscillation term obtained by LMD has a lower noise level than other signal separation methods, effectively improving the accuracy of retrieval results and avoiding abnormal values. Moreover, it still performs well under loose constraints (a wide RH range and a high-elevation range). In one consecutive year of retrieval results, the new model based on LMD has a significant improvement effect over quadratic fitting, and the root mean square error and mean absolute error of retrieval results obtained in each month on average are improved by 8.34% and 8.87%, respectively.

1. Introduction

Glaciers are melting in large quantities as global warming intensifies, leading to sea-level rise and a constant threat to human life [1]. Consequently, real-time tracking of sea-level variations and studying related patterns is of paramount importance [2]. Over the past two decades, traditional tide gauges (TG) have been used as the primary means of sea-level monitoring. However, their measured results are subject to error due to the influence of crustal movements. Although satellite altimetry can achieve high-accuracy and large-scale sea-level monitoring, its monitoring accuracy is low in the near-coastal area [3]. In recent years, with the maturity of global positioning system–interferometric reflectometry (GNSS-IR) technology, the sea-level estimation method based on this technology provides an effective solution for near-coastal sea-level monitoring. GNSS-IR technology enables the cost-effective, uninterrupted monitoring of sea levels, and the resulting observations are automatically anchored in a stable frame [4].
In 1993, Martin-Neira, a scientist at the European Space Agency (ESA), initially introduced global navigation satellite system–reflectometry (GNSS-R) technology, which has become one of the focuses of research in the field of remote sensing [5]. This technique uses the characteristic parameters of the reflected signal to detect the physical characteristics of the reflecting surface and has a high application value [6,7,8,9]. In 2008, Larson et al. further proposed the GNSS-IR technique and successfully retrieved the variation of soil moisture [10,11]. The technique estimates relevant parameters of the reflecting surface through the characteristics of signal-to-noise ratio (SNR) oscillation (such as frequency, amplitude, and initial phase). With the advantages of low cost, high flexibility, and accessible data acquisition, GNSS-IR technology has been applied in many fields, such as snow depth detection [12,13], soil moisture [14,15], storm surge [16,17], and vegetation change [18,19].
In 2013, Larson et al. systematically described the general method of retrieving sea-level height using GNSS-IR technology, therefore establishing the foundation of GNSS-IR sea-level measurement technology [20]. Lofgren et al. used SNR and phase delay analysis to retrieve the water level variation in the rough sea surface, revealing the superior performance of the SNR analysis method [21]. Given the irrationality of the static sea-level assumption in the retrieval principle, Larson et al. proposed a dynamic sea-level correction method. It was applied to Kachemak Bay, where the tidal level variation was greater than 7 m, and the retrieval results were obtained with the root mean square error (RMSE) of 2.3 cm [22]. Williams et al. found that tropospheric delay could bias sea-level estimates. To address this problem, they proposed using the Global Temperature and Pressure (GPT2w) and Vienna mapping functions (VMF1) to eliminate the effect of tropospheric delay on sea-level retrieval [23]. Jia et al. first used the SNR data of L2, L6, and L7 bands in the BeiDou system to retrieve sea-level height and validated its availability in the GNSS-IR technology for sea-level estimation [24]. Wang et al. analyzed the SNR data from three sites, PBAY, SC02, and BRST, determining the optimal azimuth range for sea-level monitoring at each site [25]. Wang et al. used wavelet analysis to process SNR arc and extract instantaneous frequencies, significantly improving the data utilization rate [26]. Wang et al. used a robust regression method to combine SNR data from quad constellations of GNSS to retrieve the sea-level height, and the accuracy was improved by about 40–75% [27]. After years of development, GNSS-IR sea-level measurement technology has matured significantly. In this technique, the accuracy of the retrieval result largely depends on the quality of the oscillatory term [28]. The traditional model uses quadratic fitting to separate the trend and oscillation terms of the SNR arc. However, this method offers inadequate signal separation and introduces complex noise in the oscillation term, which is not conducive to extracting oscillation frequency [29,30]. Although the effect of noise can be weakened to some extent by setting a retrieval range for the reflector height (RH range) and using the SNR series with a low-elevation range, it also reduces the flexibility of the technique. To this end, Wang et al. used wavelet decomposition to process the SNR arc and selected signal components containing information about the sea surface to construct the oscillation term [31]. Zhang et al. and Hu et al. processed SNR arc using empirical modal decomposition (EMD) and variational modal decomposition (VMD), respectively, which effectively improved the accuracy and stability of retrieval results in the high-elevation range [32,33]. This refined model, based on signal processing, enhances retrieval accuracy by extracting the oscillation term with as low noise as possible from the SNR arc. In addition, the B-sample [34], wavelet decomposition [31], EMD [32], VMD [33], and singular spectrum analysis (SSA) [35] have all been applied to the GNSS-IR sea-level estimation model, all of which have improved the accuracy and stability of retrieval results to some extent. Nevertheless, it is still challenging for technology to achieve long-term, accurate, stable sea-level monitoring. In previous studies, the focus has been on exploring the advantages of improved methods based on signal decomposition at high satellite-elevation ranges, overlooking the impact of the RH range on retrieval accuracy. It is important to note that the elevation range is primarily influenced by the visual range of the antenna. As short-term and long-term tide-level fluctuations often exhibit different elevation ranges, it is essential to appropriately expand the RH range to capture the comprehensive sea-level elevation variation in long-term monitoring. Moreover, the mechanism behind the signal-decomposition method remains unclear in existing research.
To further enhance the retrieval accuracy and stability of the GNSS-IR sea-level estimation model, this study proposes a GNSS-IR model based on local mean decomposition (LMD) for sea-level estimation. The model uses LMD instead of quadratic fitting to process the SNR arc to obtain the oscillation term and uses the Lomb–Scargle periodogram (LSP) to extract the oscillation frequency. LMD can decompose the SNR arc into a series of signal components, and those that contain the sea surface information are selected to construct the oscillation term with a low noise level. Observations from the SC02 site are used to analyze the performance of LMD and the effect of noise on retrieval results from RH ranges and satellite-elevation ranges, and the stability of the model is verified by retrieving the sea-level variation at Friday Harbor for one year.
The rest of this paper is organized as follows. Section 2 introduces the experimental site and data. Section 3.1 introduces the basic principle of GNSS-IR technology, and Section 3.2 briefly introduces the decomposition principle of LMD. Section 4 describes the experimental procedure in detail. In Section 5, experimental results are given and discussed. Finally, Section 6 offers the conclusion and final remarks on the paper.

2. Site and Data

The SC02 site (48°32′46.30″ N, 123°00′27.40″ W) is located northeast of Friday Harbor, Washington, USA, and is one of the observation sites of the Plate Boundary Observatory (PBO) operated by UNAVCO for EarthScope. It is equipped with a geodetic receiver TRIMBLE NETR9 and choke antenna with rectifier TRM59800.80. The antenna is erected on the bedrock by the coast, its phase center is about 5.5 m from the sea level, and the azimuth range of 50°~240° is the sea surface [25,36]. The Friday Harbor TG, about 300 m west of SC02, provides 6-min tidal level monitoring data, which is operated by the National Oceanic and Atmospheric Administration (NOAA) of the United States.
This study used the SNR data of the Global Positioning System (GPS) L1 band at the SC02 site in 2015 with a sampling interval of 15 s, and the accuracy of retrieval results was evaluated using TG data at Friday Harbor. The first Fresnel zone (FFZ) is typically used to represent the sensing range of the signal on the reflecting surface. With a fixed RH, its position and size are determined by the azimuth and elevation of the satellite, respectively. As the elevation increases, the length of the reflector zone diminishes, and the center moves closer to the antenna. In this study, the SNR series in the elevation range of 0.5°~15° was used to retrieve the sea-level height, and the number of SNR observations should be greater than 100 in one observation period. The SNR series with elevations up to 35° was also used to explore the performance of the model at high elevations. Figure 1 shows the surrounding environment of the SC02 site and the FFZ with satellite elevations of 5° and 15°, and the radius of the reflection zone is about 100 m [37].

3. Methodology

3.1. GNSS Interferometric Reflection Theory

To maximize the reception of the reflected signals from the sea surface, antennas of GNSS receivers intended for sea-level estimation are typically installed near the coast. The signal received by the receiver can be divided into two categories: one is the direct signal received directly after the satellite launch, and the other is the reflected signal reflected by the sea surface and then enters the receiver. Figure 2 shows the geometry of the GNSS-IR model, where h is the vertical distance from the antenna phase center to the sea surface, i.e., the RH; θ is the angle between the direct signal and sea level, i.e., the satellite-elevation angle; and D is the additional distance of the reflected signal compared to the direct signal.
If only one reflection from a calm reflecting surface is considered, the additional distance D can be deduced from the geometric relationship as:
D = 2 h sin θ
The interference signal composed of the direct and reflected signal is recorded by the receiver in SNR, which can be described as [20]:
S N R 2 = A d 2 + A m 2 + 2 A d A m cos ψ
where A d and A m are the direct and reflected signal amplitudes, respectively, and ψ is the phase difference between the direct and reflected signal. To ensure accurate positioning, the antenna of a standard geodetic receiver is typically designed with specific features to mitigate the impact of the multipath effect, thus allowing the received signal to satisfy the relationship A d A m . Consequently, the SNR series shows a parabolic trend in the whole by the direct signal. In contrast, the reflected signal leads to local periodic oscillation (Figure 3), in which the information on the reflecting surface is hidden [38].
SNR is a quantitative index to evaluate signal quality [39,40]. As shown in the dashed box in Figure 3, the SNR values are smaller at low-elevation angles. They are more affected by multipath, and the periodic oscillations are more significant. Therefore, extracting the characteristic parameters of SNR oscillations for low-elevation angles is more accessible than for high-elevation angles. Although the SNR series with a larger elevation range contains more complete reflector information, the effective information is often submerged in complex noise and difficult to extract accurately. In traditional GNSS-IR models, a quadratic polynomial is used to fit the trend of the SNR arc. Then the fitting term is subtracted from the SNR arc to eliminate direct signals and a small number of reflected signals, namely the term A d 2 + A m 2 in Equation (2). The oscillation term of SNR obtained can be approximated by the cosine model [10]:
S N R m = 2 A d A m cos ψ = A cos 2 π f x + ϕ
where A is the signal amplitude, f is the oscillation frequency, and ϕ is the phase. From the additional distance D , the phase difference ψ between the direct and reflected signal can be calculated as follows:
ψ = 2 π λ D = 4 π h λ sin θ
where λ is the carrier wavelength. According to Equation (4), there is a linear relationship between the phase difference ψ and the sine value sin θ of the satellite-elevation angle, thus:
2 π f = 4 π h λ
After Equation (5) is simplified, the relationship between the reflector height h and frequency f of S N R m can be obtained as follows:
h = λ f 2
The oscillation term varying with sin θ is a non-equidistant series, so it is difficult to use the fast Fourier transform (FFT) for spectral analysis [41,42]. The conventional sea-level estimation model uses LSP to extract the frequency information of the oscillation term. The frequency f corresponding to the maximum peak amplitude in the periodogram is converted into the reflector height h according to Equation (6), and then the sea surface height is obtained by unifying reflector height to the same datum of TG.

3.2. Signal Decomposition Based on LMD

According to the basic principle of GNSS-IR introduced in Section 3.1, the accurate acquisition of reflection signals is a crucial step in retrieving sea-level height. For the LSP, the input signal must have a zero-mean value, and the noise level of the signal has a significant impact on frequency extraction [28,29]. Low-order polynomials cannot accurately depict the trend change in the SNR arc, and the acquired oscillation term contains a lot of noise, which makes the frequency acquired by LSP uncertain [43]. Therefore, the GNSS-IR technique typically uses SNR series with a low-elevation range to reduce the noise source. In addition, the RH range is set according to the actual sea-level variation during retrieval, and the maximum peak amplitude is retrieved only within this range to avoid abnormal values. Although these limitations improve retrieval accuracy to some extent, they also remarkably reduce the adaptability and flexibility of the technology. As an adaptive signal analysis method, LMD can decompose a complex signal into a series of signal components with different frequencies, which has a significant effect on the processing of non-stationary signals. Therefore, this paper proposes using LMD to process the SNR arc and separate the trend item and noise components from the SNR arc to provide a more stationary and purer oscillation term for LSP.
In 2005, Smith et al. proposed a new adaptive signal analysis based on EMD, namely LMD, and successfully applied it to the processing of electroencephalogram (EEG) signals [44]. LMD can decompose a complex amplitude and frequency modulation (AM-FM) signal into a series of product functions (PF), each being the product of a local envelope and a pure frequency-modulated signal. The traditional LMD algorithm uses adjacent extreme points to calculate the local mean and local envelope. It uses the moving average algorithm to process the local mean and local envelope to construct the local mean function and envelope estimation function. Then, the local mean function is separated from the original signal to obtain the zero-mean signal, and the envelope estimation function is used to demodulate the zero-mean signal. This procedure is repeated for the demodulated signal until a purely FM signal is obtained. At this time, the product of a series of envelope estimation functions obtained in the demodulation process is the final envelope signal, and the envelope signal is multiplied by the purely FM function to obtain the first PF component. After subtracting the PF component from the original signal, the entire process is repeated for the residual signal until a monotone residual signal is produced. If the original signal is designated as x t , the decomposition outcome of LMD can be represented as [44]:
x t = i = 1 k P F i t + u k t
where u t is the monotone residual term. The traditional LMD based on the moving average algorithm is inefficient and may not be able to obtain the convergent envelope estimation function. In response, Hu et al. proposed interpolating upper and lower extremum points with cubic spline to obtain the envelope. The experimental results showed that the decomposition effect of the LMD algorithm based on cubic spline interpolation is better than that of the traditional LMD algorithm [45]. However, the envelope estimation method based on the interpolation algorithm still makes it difficult to obtain an accurate envelope. To this end, Jia et al. proposed an Empirical Optimal Envelope (EOE). The method uses tangent points instead of extreme points, approximates the optimal interpolation point position through an iterative greedy algorithm, and optimizes the envelope distance to obtain the optimal envelope. The experimental results showed that the EOE-LMD algorithm could obtain more accurate envelopes and signal components [46]. In this study, all the LMD algorithms used are in the EOE-LMD algorithm.

4. Sea-Level Estimation

According to the theory introduced in Section 3, this study combines LMD and LSP to construct a GNSS-IR sea-level estimation model based on LMD. The experimental process is mainly divided into five steps (Figure 4).
Step 1: Data preprocessing. The GNSS data were processed using RTKLIB ver. 2.42 software to obtain SNR, azimuth, and satellite-elevation data. Due to the site environment, not all SNR data come from the sea surface. Generally, the effectiveness of the SNR arc can be judged according to the azimuth and the satellite elevation. SNR arcs boasting over 100 observations within the range of azimuth 50°~240° and satellite elevation 0.5°~15° were chosen for subsequent processing.
Step 2: Signal decomposition. The SNR arc was decomposed by LMD, and the signal components containing sea surface information were used to construct the oscillation term. Concurrently, this study utilized quadratic fitting, cubic fitting, wavelet decomposition, and EMD for SNR arc processing. Their retrieval results were compared to analyze the effectiveness and accuracy of LMD.
Step 3: Frequency extraction. LSP was used to extract the effective frequency in the SNR oscillation term. The frequency was converted to RH according to Equation (6), and then the RH was converted to the sea surface height under the same datum of TG. This study retrieved RH in the range of 4~7 m, and a wider RH range was also used in exploring the performance of LMD. In addition, the retrieval values of peak amplitude less than 3 and peak-to-noise ratio (the ratio of the maximum peak amplitude to the average peak amplitude of noise in the RH range) less than 2 were removed [40].
Step 4: Accuracy evaluation. The measured sea-level height at the corresponding time was obtained by interpolating the TG data using the cubic spline method. The gross errors in retrieval results were removed based on the principle of double or triple standard deviation. The correlation coefficient (Pearson correlation coefficient, R ), RMSE, and mean absolute error (MAE) between the retrieved and measured values were calculated to evaluate the accuracy of retrieval results.
Step 5: Comparative analysis. This study first used SNR data for 7 consecutive days to verify the effectiveness of LMD. On this basis, the performance of the GNSS-IR sea-level estimation model based on LMD and the influence of noise on retrieval results were explored from two aspects of different RH ranges and satellite-elevation ranges, respectively. Finally, the data for one consecutive year were processed to verify the stability of the model in long-term sea-level monitoring.

5. Results and Discussion

5.1. Decomposition Results of LMD for SNR

The key to using LMD to extract SNR oscillation is to find the component where reflected information from the sea surface is located. For this purpose, SNR data from DOY 173–179 in 2015 at the SC02 site were processed using LMD. From the examined data over seven consecutive days, there were 380 eligible SNR arcs, the majority of which were decomposed into three to five layers of PF components (as shown in Table 1). Figure 5 shows the decomposition results of three different SNR arcs, which were processed by LMD to obtain the three, four, and five layers of PF components, respectively, and the periodogram of each component is presented. As shown in the figure, although the decomposition layers of different SNR arcs are different, the significant periodic fluctuations are concentrated in PF1 and PF2. Furthermore, their periodograms exhibit significant peaks within the RH range of 4~7 m. The RH corresponding to the maximum peak amplitude of PF2 falls within this range. In contrast, other components show scant and discontinuous fluctuations in the low-elevation range, and the RH corresponding to the significant peak in the periodogram is less than 4 m. In addition, the spectral analysis results of the residual term are also given in the figure. It can be seen that the residual term has essentially the same trend as the SNR arc, which shows a significant peak near the zero axis in the periodogram.
The sea-level variations were estimated using the PF1, PF2, PF3, and PF4 of different SNR arcs, respectively, and the results are shown in Table 2. It can be seen that the retrieval results of PF1 and PF2 have better correlation and accuracy with TG data than other PF components. Accordingly, PF1 and PF2 were selected to construct the oscillation term, which was then compared to quadratic fitting. The oscillation terms obtained by LMD are more regular and stationary (Figure 6).
Following the experimental results, the oscillation term comprising P F 1 + P F 2 was selected to retrieve the sea-level height, and the retrieval results are presented in Figure 7 and Table 3. Additionally, the component combination schemes of wavelet decomposition and EMD were obtained through many experiments. On the premise of ensuring a sufficient number of retrieval values, the combination of the components with the best accuracy in the experimental results was selected as the effective component.
As seen from Figure 7, the R between the retrieval results of different methods and measured values all reaches 0.97. Moreover, the number of retrieval values exceeded 350, with daily averages surpassing 50. The overall trend of retrieval results can better reflect the actual sea-level variation. As shown in Table 3, in descending order of accuracy: LMD > wavelet decomposition > EMD > cubic fitting > quadratic fitting. The accuracy of LMD, wavelet decomposition, and EMD is significantly better than that of low-order polynomials. Cubic fitting shows a slight improvement in accuracy compared to quadratic fitting, and the accuracy of the LMD is slightly better than that of the wavelet decomposition and EMD, although the difference is negligible. The RMSE and MAE of the new model based on LMD were 11.74 cm and 9.30 cm, improving 11.99% and 12.51% over the quadratic fitting, respectively.

5.2. Performance Analysis of Different RH Ranges

To compare the retrieval results of different RH ranges, four RH ranges of 4~7 m, 3~8 m, 2~9 m, and 0~11 m were set in the retrieval. Figure 8 shows the retrieval results of different methods that meet the quality control conditions before removing abnormal values.
As seen in Figure 8, the retrieval results of quadratic fitting stay mostly the same when the RH range is expanded from 4~7 m to 3~8 m. When the RH range is expanded to 2~9 m, a small number of abnormal values appear in the retrieval results. These values are approximately 3 m, which is above sea level and approximate to the ground height. The corresponding frequency represents the reflected information from the ground, and it shows in the periodogram that the peak of the reflected signal from the ground replaces the one from the sea surface as the main peak. As the RH range broadens to 0~11 m, more abnormal values appear in the retrieval results. Their height is concentrated around 5 m, close to the antenna height. The corresponding frequency is close to the zero axis in the periodogram, which represents the residual trend term in the oscillation term, and the peak near the zero axis replaces that of the reflected signal from the sea surface as the main peak. For LMD, when the RH range is expanded from 4~7 m to 3~8 m, the abnormal values appearing in the retrieval results are sparse and distributed. This outcome arises because the effective frequencies of some SNR series reside in higher-order components, which are eliminated as noise, thus causing the loss of effective information and an error in frequency retrieval. As the RH range continues to be expanded, the retrieval results of LMD change little, and it can maintain a high retrieval accuracy and stability all the time. With the expansion of the RH range, the change in retrieval results obtained by cubic fitting is consistent with that of the quadratic fitting. However, when the RH range is extended to 0~11 m, the number of abnormal values in the cubic fitting is significantly reduced compared to the quadratic fitting. The wavelet decomposition and EMD greatly reduce the number of abnormal values compared with the low-order polynomials, but a small number of abnormal values still appear when the RH range is extended to 0~11 m. This indicates that the wavelet decomposition and EMD can also effectively remove the influence of noise in the SNR series. However, some SNR oscillation terms still contain ground information, resulting in certain abnormal values in retrieval results. To illustrate the differences in accuracy and the number of points between different methods, Table 4 lists the accuracy evaluation of retrieval results after removing abnormal values.
Table 4 shows that the new model based on LMD always obtains retrieval results with high accuracy and stable quantity in different RH ranges. LMD can effectively remove the interference of low-frequency noise. Compared with wavelet decomposition and EMD, LMD boasts superior signal decomposition, accurately distinguishing noise from the oscillation term and efficiently circumventing abnormal values.
Due to the poor fitting effect of the quadratic polynomial, the obtained oscillation term often contains a residual trend term. Furthermore, when the SNR series contains noise from other reflecting surfaces, it cannot be removed by quadratic fitting. An oscillatory term containing a trend term and noise exhibits multiple significant peaks in the periodogram (Figure 9). Setting a range for retrieving RH can reduce the interference of other frequency components in the SNR oscillation to a certain extent and avoid abnormal values. However, it is essentially just a constrained retrieval condition. The existence of noise will still impact the accuracy and stability of identified frequency, resulting in low-accuracy retrieval results, as confirmed by the experimental outcomes presented in Section 5.2. In addition, it can be seen from the experimental results that the height of most anomalous results is greater than that of sea level, and only a few anomalous values are lower than sea level due to the loss of effective information. This fact shows that the high-frequency noise in the SNR series has no direct effect on frequency extraction.

5.3. Performance Analysis of Different Elevation Ranges

To explore the retrieval performance of the new model based on LMD at high elevation, SNR series of 0.5°~15°, 0.5°~20°, 0.5°~25°, 0.5°~30°, and 0.5°~35° were used to retrieve sea-level height, respectively. The retrieval results are illustrated in Table 5 and Figure 10.
Combining Table 5 and Figure 10, it can be found that the accuracy of the retrieval results obtained by different methods drops as the elevation range widens. Except for quadratic fitting, the number of retrieved values remains stable for other methods. However, in this process, the accuracy of LMD is always better than the other four methods.
For LSP, the uncertainty of the frequency is inversely proportional to the number of samples [47]. As the elevation range widens, the length of the SNR series input to LSP increases, which is advantageous for the solution of the frequencies. Sea-level estimation based on SNR analysis attempts to find the average frequency of the whole series. It uses the sea-level height transformed from this frequency to represent the sea-level height at the middle moment of the entire series. LSP can only identify significant fluctuations in the series, and the height obtained may be the sea-level height at any time in the observation period. When the SNR series with a wider elevation range retrieves sea-level height, the sea-level variation is more significant over a longer observation time, decreasing retrieval accuracy. Figure 11 displays the spectral analysis results for two SNR arcs with varying elevation ranges, processed using quadratic fitting and LMD, respectively. The first one produced an abnormal value in the elevation range of 0.5°~35° by quadratic fitting, while LMD obtained an effective reversal result. The second one obtained effective retrieval results at the high-elevation range after processing by both methods, but the retrieval value of LMD was more accurate. It can be found from Figure 11 that with the expansion of the elevation range, the fitting deviation of the quadratic polynomial increases, and the trend term contained in the oscillation term greatly reduces the significance of the main peak. When the elevation angle is greater than 15°, the interference components contained in the SNR series become more complex, and the number of peaks with similar frequencies increases in the periodogram. For LSP, the significant trend items and complex noise in the signal intensify the uncertainty of the obtained frequency, which may increase the deviation between the obtained and actual frequency and reduce the accuracy of the retrieval results. The noise peak may even replace the peak of RH as the main peak, resulting in abnormal values. In addition, LMD can only remove part of the noise signals far from RH and cannot change the frequency components in the RH range, and the influence of noise components in the range on retrieving RH still exists.

5.4. Stability of Long-Term Monitoring

As seen from the previous experiments, LMD can effectively weaken the influence of noise and still have excellent performance in a wider RH range. The sea-level variation at Friday Harbor is greater than 3 m during a year, and the RH range should be appropriately expanded during retrieval to meet the complex sea-level variation. To further verify the stability of the new model based on LMD for long-term monitoring, one year of SNR data from the SC02 site in 2015 was processed using the model. The accuracy of the monthly retrieval results was assessed individually and is presented in both Table 6 and Figure 12.
As seen from Table 6 and Figure 12, among the retrieval results of 12 consecutive months, there are 8 months in which the accuracy index of LMD is entirely superior to the other four methods. The average accuracy of monthly retrieval results is also better than the other four methods. The number of retrieval values from different methods remained basically at the same level. Moreover, the retrieval result accuracy of LMD, wavelet, and EMD surpasses that of low-order polynomials, barring the month of December. Compared with the retrieval results of quadratic fitting, the RMSE of LMD is minimally improved in February at 3.44% and maximally improved in January at 16.92%, and the MAE is minimally improved in February at 4.12% and maximally improved in January at 14.45%. The RMSE and MAE of monthly retrieval results are improved by 8.34% and 8.87% on average. Over the one-year retrieval results, the GNSS-IR sea-level estimation model based on LMD shows high accuracy and excellent performance.

6. Conclusions

This study developed a GNSS-IR sea-level estimation model based on LMD. In this model, the SNR arc was decomposed by LMD, and the signal components containing reflected information from the sea surface were selected to construct the oscillation term. LMD provides a lower noise oscillation term for LSP compared to quadratic fitting, therefore effectively enhancing retrieval accuracy. The results show that the new model based on LMD still has excellent performance under loose constraints, which is of great significance for sea-level monitoring in the long term and without a priori tide-level height. In the one-year retrieval results, the GNSS-IR model for sea-level estimation based on LMD shows excellent retrieval accuracy and stability.
Changes in RH range and satellite-elevation range affect retrieval results differently. As the RH range is expanded, more noise peaks are included in the retrieval range in the periodogram. An anomalous retrieval result is obtained when the peak amplitude of the noise component is larger than that of the RH. In contrast, the use of SNR arcs with a high satellite-elevation range reduces the significance of the effective peaks and produces more noise peaks with similar frequencies, reducing the accuracy of the retrieval results. The improved method based on signal decomposition can effectively remove the trend term and the noise outside the RH range in the SNR series, which is the reason for its ability to maintain high retrieval accuracy and a sufficiently large number of retrieved values despite the loose constraints (either in the wide RH range or satellite-elevation angle range). LMD shows significant advantages over other signal separation methods with its excellent signal-decomposition performance. Unfortunately, LMD cannot change the frequency components in the RH range, which is a limitation of the LMD method and all signal analysis methods, which is consistent with the experimental findings of Wang et al. [31].
The retrieval process requires other parameter settings to ensure retrieval result accuracy, such as data length, azimuth range, peak-to-noise ratio, and peak amplitude power. In the comparison experiments about different elevation ranges, the numerical advantage of points in our results is not as significant as that of Zhang et al. [32]. This is mainly due to the different criteria for determining the effective arc. Therefore, the effect of these parameters on retrieval accuracy is equally worth investigating. This paper used fixed two-layer PF components to construct the oscillatory term. LMD is an adaptive signal-decomposition method, and the number of decomposition layers varies for different SNR series. The combination of fixed components can easily cause the loss of effective information, which is unreasonable. However, the existing dynamic selection method depends on the oscillation term obtained by quadratic fitting and cannot be carried out independently. Moreover, it cannot accurately judge the boundary between the trend item and the effective component, which fails to exploit the advantages of the signal-decomposition method fully. In the future, the problem will continue to be studied by combining the adaptive layer selection method with LMD to improve the stability of the method further. In addition, the data union of multi-systems can greatly improve the temporal resolution of retrieval results, but the SNR data of different systems with different frequencies have different sensitivities to noise, and the direct combination cannot significantly improve the accuracy level of retrieval results. Therefore, we hope to improve the retrieval accuracy of multi-mode and multi-frequency by way of combining them after denoising.

Author Contributions

Z.W.: Methodology, Software, Writing—original draft. C.R.: Conceptualization, Methodology, Writing—review and editing. X.L.: Data curation, Formal analysis. Y.L.: Conceptualization, Supervision, Writing—review and editing. A.Y.: Visualization, Validation. J.L.: Data curation, Software. W.Y.: Project administration, Writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 42064003, 41901409) and the Guangxi Natural Science Foundation (Grant No. 2021GXNSFBA220046).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The GNSS data were obtained from the Plate Boundary Observatory operated by UNAVCO for EarthScope (http://www.earthscope.org (accessed on 18 March 2023)). The tide gauge data at Friday Harbor were provided by the NOAA (https://tidesandcurrents.noaa.gov/ (accessed on 18 March 2023)).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abbass, K.; Qasim, M.Z.; Song, H.; Murshed, M.; Mahmood, H.; Younis, I. A review of the global climate change impacts, adaptation, and sustainable mitigation measures. Environ. Sci. Pollut. Res. 2022, 29, 42539–42559. [Google Scholar] [CrossRef] [PubMed]
  2. Melet, A.; Meyssignac, B.; Almar, R.; Le Cozannet, G. Under-estimated wave contribution to coastal sea-level rise. Nat. Clim. Chang. 2018, 8, 234–239. [Google Scholar] [CrossRef]
  3. Jia, Y.; Xiao, K.; Lin, M.; Zhang, X. Analysis of Global Sea Level Change Based on Multi-Source Data. Remote Sens. 2022, 14, 4854. [Google Scholar] [CrossRef]
  4. Wang, X.; Niu, Z.; Chen, S.; He, X. A Correction Method of Height Variation Error Based on One SNR Arc Applied in GNSS–IR Sea-Level Retrieval. Remote Sens. 2021, 14, 11. [Google Scholar] [CrossRef]
  5. Martin-Neira, M. A passive reflectometry and interferometry system (PARIS): Application to ocean altimetry. ESA J. 1993, 17, 331–355. [Google Scholar]
  6. Auber, J.C.; Bibaut, A.; Rigal, J.M. Characterization of Multipath on Land and Sea at GPS Frequencies. In Proceedings of the 7th International Technical Meeting of the Satellite Division of the Institute of Navigation (ION GPS 1994), Salt Lake City, UT, USA, 20–23 September 1994; pp. 1155–1171. [Google Scholar]
  7. Soulat, F.; Caparrini, M.; Germain, O.; Lopez-Dekker, P.; Taani, M.; Ruffini, G. Sea state monitoring using coastal GNSS-R. Geophys. Res. Lett. 2004, 31, 1–4. [Google Scholar] [CrossRef] [Green Version]
  8. Rodriguez-Alvarez, N.; Camps, A.; Vall-Llossera, M.; Bosch-Lluis, X.; Monerris, A.; Ramos-Perez, I.; Valencia, E.; Marchan-Hernandez, J.F.; Martinez-Fernandez, J.; Baroncini-Turricchia, G.; et al. Land geophysical parameters retrieval using the interference pattern GNSS-R technique. IEEE Trans. Geosci. Remote Sens. 2010, 49, 71–84. [Google Scholar] [CrossRef]
  9. Rius, A.; Nogués-Correig, O.; Ribó, S.; Cardellach, E.; Oliveras, S.; Valencia, E.; Park, H.; Tarongí, J.M.; Camps, A.; van der Marel, H.; et al. Altimetry with GNSS-R interferometry: First proof of concept experiment. GPS Solut. 2012, 16, 231–241. [Google Scholar] [CrossRef]
  10. Larson, K.M.; Small, E.E.; Gutmann, E.D.; Bilich, A.L.; Braun, J.J.; Zavorotny, V.U. Use of GPS receivers as a soil moisture network for water cycle studies. Geophys. Res. Lett. 2008, 35, 851–854. [Google Scholar] [CrossRef] [Green Version]
  11. Larson, K.M.; Braun, J.J.; Small, E.E.; Zavorotny, V.U.; Gutmann, E.D.; Bilich, A.L. GPS multipath and its relation to near-surface soil moisture content. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2009, 3, 91–99. [Google Scholar] [CrossRef]
  12. Larson, K.M.; Gutmann, E.D.; Zavorotny, V.U.; Braun, J.J.; Williams, M.W.; Nievinski, F.G. Can we measure snow depth with GPS receivers? Geophys. Res. Lett. 2009, 36, L17502. [Google Scholar] [CrossRef] [Green Version]
  13. Larson, K.M.; Nievinski, F.G. GPS snow sensing: Results from the EarthScope Plate Boundary Observatory. GPS Solut. 2013, 17, 41–52. [Google Scholar] [CrossRef]
  14. Lv, J.; Zhang, R.; Tu, J.; Liao, M.; Pang, J.; Yu, B.; Li, K.; Xiang, W.; Fu, Y.; Liu, G. A GNSS-IR Method for Retrieving Soil Moisture Content from Integrated Multi-Satellite Data That Accounts for the Impact of Vegetation Moisture Content. Remote Sens. 2021, 13, 2442. [Google Scholar] [CrossRef]
  15. Liang, Y.; Lai, J.; Ren, C.; Lu, X.; Zhang, Y.; Ding, Q.; Hu, X. GNSS-IR Multisatellite Combination for Soil Moisture Retrieval Based on Wavelet Analysis Considering Detection and Repair of Abnormal Phases. Measurement 2022, 203, 111881. [Google Scholar] [CrossRef]
  16. Peng, D.; Hill, E.M.; Li, L.; Switzer, A.D.; Larson, K.M. Application of GNSS interferometric reflectometry for detecting storm surges. GPS Solut. 2019, 23, 47. [Google Scholar] [CrossRef] [Green Version]
  17. Larson, K.M.; Lay, T.; Yamazaki, Y.; Cheung, K.F.; Ye, L.; Williams, S.D.; Davis, J.L. Dynamic sea level variation from GNSS: 2020 Shumagin earthquake tsunami resonance and Hurricane Laura. Geophys. Res. Lett. 2021, 48, e2020GL091378. [Google Scholar] [CrossRef]
  18. Zhang, S.; Wang, T.; Wang, L.; Zhang, J.; Peng, J.; Liu, Q. Evaluation of GNSS-IR for retrieving soil moisture and vegetation growth characteristics in wheat farmland. J. Surv. Eng. 2021, 147, 04021009. [Google Scholar] [CrossRef]
  19. Sui, M.; Chen, K.; Shen, F. Monitoring of Wheat Height Based on Multi-GNSS Reflected Signals. Remote Sens. 2022, 14, 4955. [Google Scholar] [CrossRef]
  20. Larson, K.M.; Löfgren, J.S.; Haas, R. Coastal sea level measurements using a single geodetic GPS receiver. Adv. Space Res. 2013, 51, 1301–1310. [Google Scholar] [CrossRef] [Green Version]
  21. Löfgren, J.S.; Haas, R. Sea level measurements using multi-frequency GPS and GLONASS observations. EURASIP J. Adv. Signal Process. 2014, 2014, 50. [Google Scholar] [CrossRef] [Green Version]
  22. Larson, K.M.; Ray, R.D.; Nievinski, F.G.; Freymueller, J.T. The accidental tide gauge: A GPS reflection case study from Kachemak Bay, Alaska. IEEE Geosci. Remote Sens. Lett. 2013, 10, 1200–1204. [Google Scholar] [CrossRef] [Green Version]
  23. Williams, S.D.P.; Nievinski, F.G. Tropospheric delays in ground-based GNSS multipath reflectometry—Experimental evidence from coastal sites. J. Geophys. Res. Solid Earth 2017, 122, 2310–2327. [Google Scholar] [CrossRef] [Green Version]
  24. Jin, S.; Qian, X.; Wu, X. Sea level change from BeiDou Navigation Satellite System-Reflectometry (BDS-R): First results and evaluation. Glob. Planet. Chang. 2017, 149, 20–25. [Google Scholar] [CrossRef]
  25. Wang, X.; Zhang, Q.; Zhang, S. Azimuth selection for sea level measurements using geodetic GPS receivers. Adv. Space Res. 2018, 61, 1546–1557. [Google Scholar] [CrossRef]
  26. Wang, X.; Zhang, Q.; Zhang, S. Sea level estimation from SNR data of geodetic receivers using wavelet analysis. GPS Solut. 2019, 23, 6. [Google Scholar] [CrossRef]
  27. Wang, X.; He, X.; Zhang, Q. Evaluation and combination of quad-constellation multi-GNSS multipath reflectometry applied to sea level retrieval. Remote Sens. Environ. 2019, 231, 111229. [Google Scholar] [CrossRef]
  28. Purnell, D.; Gomez, N.; Chan, N.H.; Strandberg, J.; Holland, D.M.; Hobiger, T. Quantifying the uncertainty in ground-based GNSS-reflectometry sea level measurements. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 4419–4428. [Google Scholar] [CrossRef]
  29. VanderPlas, J.T. Understanding the lomb–scargle periodogram. Astrophys. J. Suppl. Ser. 2018, 236, 16. [Google Scholar] [CrossRef]
  30. Bilich, A.; Larson, K.M. Correction published 29 March 2008: Mapping the GPS multipath environment using the signal-to-noise ratio (SNR). Radio Sci. 2007, 42, 1–16. [Google Scholar]
  31. Wang, X.; Zhang, Q.; Zhang, S. Water levels measured with SNR using wavelet decomposition and Lomb–Scargle periodogram. GPS Solut. 2018, 22, 22. [Google Scholar] [CrossRef]
  32. Zhang, S.; Liu, K.; Liu, Q.; Zhang, C.; Zhang, Q.; Nan, Y. Tide variation monitoring based improved GNSS-MR by empirical mode decomposition. Adv. Space Res. 2019, 63, 3333–3345. [Google Scholar] [CrossRef]
  33. Hu, Y.; Yuan, X.; Liu, W.; Wickert, J.; Jiang, Z.; Haas, R. GNSS-IR Model of Sea Level Height Estimation Combining Variational Mode Decomposition. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 14, 10405–10414. [Google Scholar] [CrossRef]
  34. Strandberg, J.; Hobiger, T.; Haas, R. Improving GNSS-R sea level determination through inverse modeling of SNR data. Radio Sci. 2016, 51, 1286–1296. [Google Scholar] [CrossRef] [Green Version]
  35. Vu, P.L.; Ha, M.C.; Frappart, F.; Darrozes, J.; Ramillien, G.; Dufrechou, G.; Gegout, P.; Morichon, D.; Bonneton, P. Identifying 2010 Xynthia storm signature in GNSS-R-based tide records. Remote Sens. 2019, 11, 782. [Google Scholar] [CrossRef] [Green Version]
  36. Larson, K.M.; Ray, R.D.; Williams, S.D.P. A 10-year comparison of water levels measured with a geodetic GPS receiver versus a conventional tide gauge. J. Atmos. Ocean. Technol. 2017, 34, 295–307. [Google Scholar] [CrossRef] [Green Version]
  37. Roesler, C.; Larson, K.M. Software tools for GNSS interferometric reflectometry (GNSS-IR). GPS Solut. 2018, 22, 80. [Google Scholar] [CrossRef] [Green Version]
  38. Song, M.; He, X.; Wang, X.; Zhou, Y.; Xu, X. Study on the quality control for periodogram in the determination of water level using the GNSS-IR technique. Sensors 2019, 19, 4524. [Google Scholar] [CrossRef] [Green Version]
  39. Larson, K.M.; Small, E.E.; Gutmann, E.; Bilich, A.; Axelrad, P.; Braun, J. Using GPS multipath to measure soil moisture fluctuations: Initial results. GPS Solut. 2008, 12, 173–177. [Google Scholar] [CrossRef]
  40. Löfgren, J.S.; Haas, R.; Scherneck, H.G. Sea level time series and ocean tide analysis from multipath signals at five GPS sites in different parts of the world. J. Geodyn. 2014, 80, 66–80. [Google Scholar] [CrossRef] [Green Version]
  41. Lomb Nicholas, R. Least-squares frequency analysis of unequally spaced data. Astrophys. Space Sci. 1976, 39, 447–462. [Google Scholar] [CrossRef]
  42. Scargle, J.D. Studies in astronomical time series analysis. II-Statistical aspects of spectral analysis of unevenly spaced data. Astrophys. J. 1982, 263, 835–853. [Google Scholar] [CrossRef]
  43. Ran, Q.; Zhang, B.; Yao, Y.; Yan, X.; Li, J. Editing arcs to improve the capacity of GNSS-IR for soil moisture retrieval in undulating terrains. GPS Solut. 2022, 26, 19. [Google Scholar] [CrossRef]
  44. Smith, J.S. The local mean decomposition and its application to EEG perception data. J. R. Soc. Interface 2005, 2, 443–454. [Google Scholar] [CrossRef] [PubMed]
  45. Hu, J.S.; Yang, S.X.; Ren, D.Q. Spline-based local mean decomposition method for vibration signal. J. Data Acquis. Process. 2009, 24, 82–86. (In Chinese) [Google Scholar]
  46. Jia, L.; Zhang, Q.; Zheng, X.; Yao, P.; He, X.; Wei, X. The empirical optimal envelope and its application to local mean decomposition. Digit. Signal Process. 2019, 87, 166–177. [Google Scholar] [CrossRef]
  47. Matz, V.; Ramos, P.M.; Brás, N.B.; Serra, A.C.A. Comparative Evaluation Between Frequency Estimation Algorithms for Power Quality Assesment in DSP Implementation. In Proceedings of the XVIII IMEKO World Congress, Metrology for a Sustainable Development, Department of Electrical and Computer Engineering, Instituto Superior Tcnico, Technical University of Lisbon, Rio de Janeiro, Brazil, 17–22 September 2006. [Google Scholar]
Figure 1. (Left) Surroundings of the SC02 site (http://www.unavco.org (accessed on 18 March 2023)). (Right) First Fresnel zone of the SC02 site. Fresnel reflection zones of 5° and 15° are plotted on the Google Earth image, with the yellow pin indicating site location and ellipses of different colors corresponding to reflection zones of different satellites.
Figure 1. (Left) Surroundings of the SC02 site (http://www.unavco.org (accessed on 18 March 2023)). (Right) First Fresnel zone of the SC02 site. Fresnel reflection zones of 5° and 15° are plotted on the Google Earth image, with the yellow pin indicating site location and ellipses of different colors corresponding to reflection zones of different satellites.
Sensors 23 06540 g001
Figure 2. Schematic diagram of GNSS-IR geometry.
Figure 2. Schematic diagram of GNSS-IR geometry.
Sensors 23 06540 g002
Figure 3. Variation of SNR and elevation angle of PRN 22 satellite on day of year (DOY) 173 in 2015. The black and red lines indicate the SNR and satellite-elevation variation with time, respectively, and the dashed boxes indicate SNR values for elevation angles of 0.5°~15°.
Figure 3. Variation of SNR and elevation angle of PRN 22 satellite on day of year (DOY) 173 in 2015. The black and red lines indicate the SNR and satellite-elevation variation with time, respectively, and the dashed boxes indicate SNR values for elevation angles of 0.5°~15°.
Sensors 23 06540 g003
Figure 4. Flowchart of experiments and analysis of GNSS-IR sea-level estimation model based on LMD.
Figure 4. Flowchart of experiments and analysis of GNSS-IR sea-level estimation model based on LMD.
Sensors 23 06540 g004
Figure 5. Decomposition results of three SNR arcs. Rows: (top) SNR arc 1; (middle) SNR arc 2; (bottom) SNR arc 3. Columns: (left) SNR series variation; (middle) Components obtained by decomposition; (right) Spectral analysis of each component. Reference lines (dashed black lines) are added to the sine at 35° elevation and RH positions at 4 m and 7 m.
Figure 5. Decomposition results of three SNR arcs. Rows: (top) SNR arc 1; (middle) SNR arc 2; (bottom) SNR arc 3. Columns: (left) SNR series variation; (middle) Components obtained by decomposition; (right) Spectral analysis of each component. Reference lines (dashed black lines) are added to the sine at 35° elevation and RH positions at 4 m and 7 m.
Sensors 23 06540 g005
Figure 6. The oscillation terms obtained by LMD and quadratic fitting for different arcs. (Top) SNR arc 1; (Middle) SNR arc 2; (Bottom) SNR arc 3.
Figure 6. The oscillation terms obtained by LMD and quadratic fitting for different arcs. (Top) SNR arc 1; (Middle) SNR arc 2; (Bottom) SNR arc 3.
Sensors 23 06540 g006
Figure 7. Sea-level retrieval results for 7 consecutive days by different methods. The SNR arc was decomposed by 6 layers using Daubechies4 wavelet, and 1–4 layers of components were selected to construct the oscillation term. EMD is an adaptive signal-decomposition method. When it is used to process SNR arc, 5–6 layers of intrinsic mode functions (IMF) are generally obtained, and 1–3 layers of IMF components were selected to construct oscillation term.
Figure 7. Sea-level retrieval results for 7 consecutive days by different methods. The SNR arc was decomposed by 6 layers using Daubechies4 wavelet, and 1–4 layers of components were selected to construct the oscillation term. EMD is an adaptive signal-decomposition method. When it is used to process SNR arc, 5–6 layers of intrinsic mode functions (IMF) are generally obtained, and 1–3 layers of IMF components were selected to construct oscillation term.
Sensors 23 06540 g007
Figure 8. Retrieval results of the five methods in different RH ranges. Rows: from top to bottom are the RH ranges of 4~7 m, 3~8 m, 2~9 m, and 0~11 m, respectively. Columns: from left to right are the retrieval results of LMD, quadratic fitting, cubic fitting, wavelet decomposition, and EMD, respectively.
Figure 8. Retrieval results of the five methods in different RH ranges. Rows: from top to bottom are the RH ranges of 4~7 m, 3~8 m, 2~9 m, and 0~11 m, respectively. Columns: from left to right are the retrieval results of LMD, quadratic fitting, cubic fitting, wavelet decomposition, and EMD, respectively.
Sensors 23 06540 g008
Figure 9. Periodogram of the oscillation term. The frequency of the horizontal axis in the periodogram is converted to RH. The blue circle represents the maximum peak amplitude retrieved in the RH ranges of 4~7 m, 2~9 m, and 0~11 m, respectively, and the coordinate position of the peak is marked.
Figure 9. Periodogram of the oscillation term. The frequency of the horizontal axis in the periodogram is converted to RH. The blue circle represents the maximum peak amplitude retrieved in the RH ranges of 4~7 m, 2~9 m, and 0~11 m, respectively, and the coordinate position of the peak is marked.
Sensors 23 06540 g009
Figure 10. Accuracy comparison of retrieval results by five methods in different elevation ranges. (Top left) R between the retrieval results of five methods in different elevation ranges and TG data; (Top right) RMSE; (Bottom left) MAE; (Bottom right) The number of retrieval values by five methods in different elevation ranges.
Figure 10. Accuracy comparison of retrieval results by five methods in different elevation ranges. (Top left) R between the retrieval results of five methods in different elevation ranges and TG data; (Top right) RMSE; (Bottom left) MAE; (Bottom right) The number of retrieval values by five methods in different elevation ranges.
Sensors 23 06540 g010
Figure 11. Periodograms of the oscillation terms at different elevation angles (Columns) obtained by different methods (Rows) for two SNR arcs. (Top) SNR arc 1, (Bottom) SNR arc 2. Each picture shows the variation of power spectral density with RH (red line) and the maximum peak amplitude retrieved in the range of 4~7 m (blue circles). The absolute deviation between the retrieval and measured value at the corresponding time is marked in the figure (unit: cm).
Figure 11. Periodograms of the oscillation terms at different elevation angles (Columns) obtained by different methods (Rows) for two SNR arcs. (Top) SNR arc 1, (Bottom) SNR arc 2. Each picture shows the variation of power spectral density with RH (red line) and the maximum peak amplitude retrieved in the range of 4~7 m (blue circles). The absolute deviation between the retrieval and measured value at the corresponding time is marked in the figure (unit: cm).
Sensors 23 06540 g011
Figure 12. Accuracy comparison of monthly retrieval results by different methods. (Top left) R between the monthly obtained retrieval results by different methods and TG data; (Top right) RMSE; and (Bottom left) MAE; (Bottom right) The number of monthly obtained retrieval values by different methods.
Figure 12. Accuracy comparison of monthly retrieval results by different methods. (Top left) R between the monthly obtained retrieval results by different methods and TG data; (Top right) RMSE; and (Bottom left) MAE; (Bottom right) The number of monthly obtained retrieval values by different methods.
Sensors 23 06540 g012
Table 1. Statistical results of decomposition layers for all effective arcs.
Table 1. Statistical results of decomposition layers for all effective arcs.
Number of PF ComponentsNumber of SNR Series
335
4256
581
67
71
Table 2. The retrieval results of different PF components.
Table 2. The retrieval results of different PF components.
ResultsPF1PF2PF3PF4
R 0.72520.95250.08661.0000
RMSE (cm)47.3917.27109.92108.47
MAE (cm)28.1611.9084.96108.41
Effective Points115348622
Table 3. Accuracy comparison of results obtained by different methods. The RH range was set to 4~7 m, and the gross errors larger than the triple standard deviation in the search results were removed.
Table 3. Accuracy comparison of results obtained by different methods. The RH range was set to 4~7 m, and the gross errors larger than the triple standard deviation in the search results were removed.
Methods R RMSE (cm)MAE (cm)Effective Points
LMD0.978111.749.30363
Quadratic fitting0.973413.3410.63364
Cubic fitting0.974012.9910.28370
Wavelet0.977511.949.47364
EMD0.977811.969.59365
Table 4. Accuracy evaluation of retrieval results in different RH ranges for five methods. Abnormal values retrieved outside the RH range of 4~7 m and gross errors larger than the double standard deviation in the retrieval results were removed.
Table 4. Accuracy evaluation of retrieval results in different RH ranges for five methods. Abnormal values retrieved outside the RH range of 4~7 m and gross errors larger than the double standard deviation in the retrieval results were removed.
ResultsMethods4~7 m3~8 m2~9 m0~11 m
R LMD0.98030.97740.97740.9774
Quadratic fitting0.97710.97700.97640.9844
Cubic fitting0.97650.97690.97600.9791
Wavelet0.98110.97830.97960.9770
EMD0.98030.98030.98030.9805
RMSE (cm)LMD11.0112.0011.9911.99
Quadratic fitting12.4712.4912.5112.39
Cubic fitting12.2912.2212.2811.22
Wavelet10.9311.6911.2912.09
EMD11.2211.2211.0911.04
MAE (cm)LMD8.879.449.449.44
Quadratic fitting10.1610.1810.2510.31
Cubic fitting9.899.879.939.26
Wavelet8.849.339.119.61
EMD9.139.139.069.01
Effective PointsLMD356366367367
Quadratic fitting358359349114
Cubic fitting365367353222
Wavelet353363354361
EMD357359350343
Table 5. Accuracy evaluation of retrieval results of five methods in different elevation ranges. The gross errors in retrieval results were eliminated based on the double standard deviation.
Table 5. Accuracy evaluation of retrieval results of five methods in different elevation ranges. The gross errors in retrieval results were eliminated based on the double standard deviation.
ResultsMethods0.5°~15°0.5°~20°0.5°~25°0.5°~30°0.5°~35°
R LMD0.98030.97030.96690.96380.9618
Quadratic fitting0.97710.96350.95660.95090.9489
Cubic fitting0.97650.96800.95980.95420.9510
Wavelet0.98110.96900.96330.95620.9469
EMD0.98030.96860.96520.9570.9524
RMSE (cm)LMD11.0113.0613.9114.5715.00
Quadratic fitting12.4714.5115.7817.3818.79
Cubic fitting12.2913.5715.3816.7517.79
Wavelet10.9313.5814.8116.2118.36
EMD11.2213.5314.4616.1817.67
MAE (cm)LMD8.8710.2511.0611.7112.18
Quadratic fitting10.1611.4412.4813.6115.10
Cubic fitting9.8910.9012.2213.1814.27
Wavelet8.8410.8211.9012.9814.40
EMD9.1310.7211.7213.0014.24
Effective PointsLMD356371373370375
Quadratic fitting358347327330333
Cubic fitting365373368364375
Wavelet353376378381390
EMD357375374380387
Table 6. Accuracy statistics of monthly retrieval results by different methods. The RH range was set to 3~8 m, the gross errors larger than the double standard deviation were removed from the retrieval results, and the accuracy of retrieval results was evaluated each month.
Table 6. Accuracy statistics of monthly retrieval results by different methods. The RH range was set to 3~8 m, the gross errors larger than the double standard deviation were removed from the retrieval results, and the accuracy of retrieval results was evaluated each month.
MonthsMethods R RMSE (cm)MAE (cm)Effective Points
JanLMD0.990011.889.461478
Quadratic fitting0.986314.3011.061488
Cubic fitting0.985914.2310.771497
Wavelet0.987113.5310.351528
EMD0.986713.8110.511522
FebLMD0.987512.079.651361
Quadratic fitting0.987212.5010.061321
Cubic fitting0.984713.2910.431379
Wavelet0.988011.749.511364
EMD0.988111.649.451360
MarLMD0.983611.609.301577
Quadratic fitting0.981212.5810.021554
Cubic fitting0.980612.589.921582
Wavelet0.981912.209.751593
EMD0.982711.899.571584
AprLMD0.981512.289.851431
Quadratic fitting0.978413.6210.861430
Cubic fitting0.978213.4510.531458
Wavelet0.980412.8610.271451
EMD0.980712.7610.191455
MayLMD0.985412.469.941631
Quadratic fitting0.981914.1311.211637
Cubic fitting0.982113.8510.881651
Wavelet0.984512.8310.281653
EMD0.984412.9710.301656
JunLMD0.986813.1210.371551
Quadratic fitting0.986013.8011.131522
Cubic fitting0.984714.1311.121563
Wavelet0.987412.8710.291554
EMD0.986513.4010.631568
JulLMD0.987312.9510.071586
Quadratic fitting0.986713.7510.971551
Cubic fitting0.985014.2210.911607
Wavelet0.987612.8810.191590
EMD0.988012.7510.101587
AugLMD0.977414.7211.091598
Quadratic fitting0.975115.7712.241583
Cubic fitting0.975015.4611.651599
Wavelet0.976414.9911.421617
EMD0.976515.0911.441618
SepLMD0.974613.5110.391545
Quadratic fitting0.971414.7611.541540
Cubic fitting0.973413.7610.641531
Wavelet0.975213.2510.341558
EMD0.975213.3810.461556
OctLMD0.982712.359.871558
Quadratic fitting0.981113.1410.661516
Cubic fitting0.979713.2710.551543
Wavelet0.981812.7110.191565
EMD0.981512.8010.221565
NovLMD0.987412.209.611494
Quadratic fitting0.986213.1410.381468
Cubic fitting0.985712.9610.091488
Wavelet0.986712.599.841504
EMD0.987312.319.691493
DecLMD0.988113.2510.441584
Quadratic fitting0.985914.7611.591568
Cubic fitting0.987413.6110.731580
Wavelet0.987613.6110.721597
EMD0.987213.8510.881598
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wei, Z.; Ren, C.; Liang, X.; Liang, Y.; Yin, A.; Liang, J.; Yue, W. Sea-Level Estimation from GNSS-IR under Loose Constraints Based on Local Mean Decomposition. Sensors 2023, 23, 6540. https://doi.org/10.3390/s23146540

AMA Style

Wei Z, Ren C, Liang X, Liang Y, Yin A, Liang J, Yue W. Sea-Level Estimation from GNSS-IR under Loose Constraints Based on Local Mean Decomposition. Sensors. 2023; 23(14):6540. https://doi.org/10.3390/s23146540

Chicago/Turabian Style

Wei, Zhenkui, Chao Ren, Xingyong Liang, Yueji Liang, Anchao Yin, Jieyu Liang, and Weiting Yue. 2023. "Sea-Level Estimation from GNSS-IR under Loose Constraints Based on Local Mean Decomposition" Sensors 23, no. 14: 6540. https://doi.org/10.3390/s23146540

APA Style

Wei, Z., Ren, C., Liang, X., Liang, Y., Yin, A., Liang, J., & Yue, W. (2023). Sea-Level Estimation from GNSS-IR under Loose Constraints Based on Local Mean Decomposition. Sensors, 23(14), 6540. https://doi.org/10.3390/s23146540

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