Next Article in Journal
On-Demand Satellite Payload Execution Strategy for Natural Disasters Monitoring Using LoRa: Observation Requirements and Optimum Medium Access Layer Mechanisms
Next Article in Special Issue
Analysis and Mitigation of Crosstalk Effect on Coastal GNSS-R Code-Level Altimetry Using L5 Signals from QZSS GEO
Previous Article in Journal
Assessing a Prototype Database for Comprehensive Global Aquatic Land Cover Mapping
Previous Article in Special Issue
An Improved Method of Soil Moisture Retrieval Using Multi-Frequency SNR Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Kalman Filter Soil Moisture Inversion Model Using GPS SNR Data—A Dual-Band Data Fusion Approach

1
Institute of Space Science, Shandong University, Weihai 264209, China
2
College of Information Science and Engineering, Shandong Agricultural University, Tai’an 271018, China
3
School of Electronic and Information Engineering, Beihang University, Beijing 100191, China
4
School of Geological Engineering and Surveying and Mapping, Chang’an University, Xi’an 710054, China
5
China Association of Remote Sensing Application, Beijing 100094, China
6
Business School, Jinhua Polytechnic, Jinhua 321000, China
7
Laboratoire Géosciences Environnement Toulouse, Université Paul Sabatier, 31400 Toulouse, France
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(19), 4013; https://0-doi-org.brum.beds.ac.uk/10.3390/rs13194013
Submission received: 23 August 2021 / Revised: 29 September 2021 / Accepted: 5 October 2021 / Published: 7 October 2021
(This article belongs to the Special Issue Recent Advances in GNSS Reflectometry)

Abstract

:
This article aims to attempt to increase the number of satellites that can be used for monitoring soil moisture to obtain more precise results using GNSS-IR (Global Navigation Satellite System-Interferometric Reflectometry) technology to estimate soil moisture. We introduce a soil moisture inversion model by using GPS SNR (Signal-to-Noise Ratio) data and propose a novel Robust Kalman Filter soil moisture inversion model based on that. We validate our models on a data set collected at Lamasquère, France. This paper also compares the precision of the Robust Kalman Filter model with the conventional linear regression method and robust regression model in three different scenarios: (1) single-band univariate regression, by using only one observable feature such as frequency, amplitude, or phase; (2) dual-band data fusion univariate regression; and (3) dual-band data fusion multivariate regression. First, the proposed models achieve higher accuracy than the conventional method for single-band univariate regression, especially by using the phase as the input feature. Second, dual-band univariate data fusion achieves higher accuracy than single-band and the result of the Robust Kalman Filter model correlates better to the in situ measurement. Third, multivariate variable fusion improves the accuracy for both models, but the Robust Kalman Filter model achieves better improvement. Overall, the Robust Kalman Filter model shows better results in all the scenarios.

1. Introduction

Soil moisture is an essential component in the continental water–carbon cycle, and a key parameter for quantifying the energy and water exchange between land surface and atmosphere [1,2,3,4]. Precise soil moisture monitoring is an important foundation for achieving high yields of agriculture production. L-band microwaves have significant advantages in soil moisture remote sensing, for being little affected by clouds and atmosphere, and for its ability to penetrate vegetation. GNSS constellations can provide massive L-band signal sources for free [5,6,7]; a novel microwave remote sensing approach, GNSS-R (Global Navigation Satellite System-Reflectometry) technology was proposed by M. Martin-Neira in 1993 [8], which typically uses an up-looking antenna for receiving direct signals from GNSS satellites, and two down-looking antennas for receiving LHCP (left-hand circular polarization) and RHCP (right-hand circular polarization) reflected signals from the Earth’s surface. At first this method was mainly used for sea-states monitoring [9], and thereafter for land parameter observations such as soil moisture [10,11,12,13]. In 2008, K. Larson developed a new method called GNSS-IR [14], which uses a single RHCP antenna to receive both the direct and reflected signal simultaneously, with a Geodetic GNSS receiver to record the interference of the two signals in the SNR (Signal-to-Noise Ratio) file. Later, Zavorotny and Larson explained the GNSS-IR theoretically by developing a physical model of the interference mechanism together [15]. Compared to the conventional GNSS-R technology, GNSS-IR uses off-the-shelf commercial receivers and antennas, so that it has a great merit of low-cost, which makes it easier to build large scale in situ monitoring systems for soil moisture measurement.
Nievinski studied the forward model of GPS multipath signal for near-surface reflectometry thoroughly [16], and later built an open-source GPS multipath simulator in Matlab/Octave for demonstrating the impact of environmental observables such as soil moisture and snow density [17]. However, as the land surface is complex and varied, simulators may not fully simulate the real situation. Chew proposed a method for determining whether SNR data are significantly corrupted by vegetation and for correcting these effects [18]. Recently, leading GNSS-R groups are only just starting to grapple with the various complexities of land GNSS-R reflections and developing complementary models [19,20,21]. In addition, there are scholars that have combined remote sensing technology and GNSS-IR technology in estimating vegetation water content [22]. Han proposed a method to reconstruct direct and multipath signal from SNR data and then calculate the dielectric constant of soil [23]. Ting Yang retrieved the dielectric constant from BDS (BeiDou Navigation Satellite System) SNR data by using an analytical model and verified the applicability though experimentation [24].
R. E. Kalman introduced his famous discrete data filtering technique in 1960 [25]. The Kalman Filter algorithm has been extensively applied in GNSS [22], INS [26], data assimilation [27] and many other research fields for its ability to provide an efficient way of computing the least squares problem using a recursive method. As of recent years, the Kalman Filter has been applied in GNSS-R sea level monitoring [28] and sea wind retrieval [29,30,31], but little research has been conducted regarding soil moisture inversion.
The reflected signal of different frequency carriers has different information than the reflected surface, therefore the data fusion algorithm was introduced to GNSS-R applications. Nazi Wang proposed a Sea Level Estimation method base on GNSS dual-band carrier phase linear combinations and achieved altimetric accuracy <0.2 m. Wang et al. used a peak weighting method to fuse GPS L1 and L2 SNR for wind speed retrieval [32]. Most of the existing GNSS-IR soil moisture retrieval methods are focused on building an empirical model and only use a single-band of GNSS signal. However, a multi-band data fusion algorithm for GNSS-IR soil moisture measurement has rarely been studied directly. Therefore, in this study we have established a Robust Kalman Filter algorithm to retrieve soil moisture and to improve the robustness and accuracy of the retrieval.
This article aims to examine and evaluate the potential of data fusion algorithms for GNSS-R applications. First, an author-proposed robust regression model will be introduced, as well as an optimized novel Robust Kalman Filter model. For validation, a dataset collected at Lamasquère (South France) will be analyzed, and the third section presents the results obtained with both robust regression model and Robust Kalman Filter model. These model retrievals are compared to (1) a classical model obtain using an empirical regression model and to (2) measurements of in situ Theta Probe ML3 sensor. The last section concludes this study by highlighting pros and cons of the different methods.

2. Methods

GNSS reflectometry works like Bi-static Radar systems, which consider the GNSS satellites and receiver as the radar transmitter and receiver, respectively. The main idea of GNSS-IR technology is to use a single RHCP antenna to receive both the direct and reflected signal simultaneously [11]. In the scenario of in situ observation with low antenna height (the reflected over-path must be less than L1 C/A code wavelength ~293 m), the difference of Doppler shift between direct and reflected signal can be ignored because our application scenario is ground based measurement. If it is being used with space-borne GNSS-R receivers, the Doppler shift cannot be ignored. At the same time, if the transmit path’s difference is less than one chip, so that the two signals are coherent, then the two signals are interfering each other at the center of antenna. The scenario of interference generation is demonstrated in Figure 1.
The reflected signal contains more RHCP component in the low elevation angle scenario, so that the interference is more significant in that case. The interference phenomenon is recorded in the SNR data stored by the geodetic GNSS receiver. Following [16,17], SNR can be expressed as a combination of direct and reflected signal as shown in Equation (1):
S N R 2 = A d 2 ( θ ) + A r 2 ( θ ) + 2 A d ( θ ) A r ( θ ) cos ( φ )
where A d , A r is the amplitude of direct and reflected signal, respectively, φ is the phase difference between reflected and direct signal, and θ is the elevation angle of GNSS satellite. When the elevation angle θ changes, the phase difference φ also changes, by which the oscillation in the SNR amplitude is created.
GNSS-IR technology normally utilizes geodetic receivers and antennas. In that case, the antenna’s gain pattern is designed to suppress the multipath signal coming from the bottom side for better positioning accuracy, i.e., the amplitude of SNR is mainly contributed by the direct signal (see antenna gain pattern, Figure 1).
Supposing δ is the difference between transmit paths of the reflected signal and direct signal, then φ can be expressed as Equation (2):
φ = 2 π λ δ = 4 π h λ sin θ
where h is the effective height of the antenna, which refers to the distance from antenna phase center to the reflecting plane, and λ is the wavelength of GNSS carrier wave.
Due to the microwave that penetrates the soil a few centimeters or decimeters (depending on the ground composition and the soil moisture) when reflected by soil surface, the reflecting plane is beneath the soil surface, so the effective antenna height varies in a range of few centimeters/decimeters with respect to the soil moisture, a key parameter in determining the penetration depth of the microwave. As the moisture usually does not change much in a few hours without precipitation, the effective antenna height can be treated as a constant to simplify the analysis.
The modified oscillation frequency can be derived by Equation (3):
f x = d φ d ( sin θ ) = 4 π h λ
For soil moisture measurement the direct signal is not of interest and typically can be removed by a second-order polynomial fitting and then derive the reflected component only, which refers to the multipath component that can be expressed by Equation (4):
S N R m ( θ , ε ) = A m cos ( 4 π h λ sin θ + ϕ ( θ , ε ) )
where ε is the dielectric constant of soil, A m and ϕ are the amplitude and phase of multipath oscillation, respectively, and h is the effective antenna height. Typical data processing method is deriving the frequency f x and h by using Lomb–Scargle Spectral Analysis (LSSA), and then deriving A m and ϕ by using least mean square estimation, which corresponds to a conventional empirical model used to retrieve soil moisture.
Existing studies show that normally A m and ϕ have higher correlation to soil moisture than f x and h [13,14], so in this study we use A m and ϕ to build soil moisture inversion models. Here, we define jth day amplitude and phase observable on Lk band as a j k , p j k (k = (1, 2)), and then the observables’ time series vectors are A k , P k which are constituted by daily observables a j k , p j k , respectively. In the same way the soil moisture time series vector can be defined as Y R +   N which is constituted by daily soil moisture y j . In general, a GNSS-IR soil moisture retrieve model is a map function from A k , P k to Y .
In this section, the conventional linear regression model and an author-proposed robust regression model will be introduced. Furthermore, a novel Kalman Filter model will be proposed for soil moisture inversion.

2.1. Conventional Linear Regression Soil Moisture Inversion Model

The general idea of conventional linear regression is to build a map from independent variable to dependent variable matrix, as shown as Equation (5):
Y = β X + v
where β is different coefficients being calculated for each satellite at each frequency band and v is the residual vector. Y is the soil moisture time series for the whole measuring period. X is observation or its combination.
In this study we use three different strategies to make the map for GNSS-IR linear regression.
In the first strategy, we used a single-band univariate linear regression. In this case, X is the single vector A k or P k . In the second strategy we used dual-bands data fusion univariate linear regression; in this case we can define two matrices A = [ A 1 , A 2 ] , P = [ P 1 , P 2 ] and X can be either A or P . In the last strategy we used dual-bands data fusion multivariate linear regression. In this case we build a joint observable matrix [ A P ] .
In the classical process of linear regression, the least mean square method is generally used to determine β and v , but the disturbance of outlier noise cannot be effectively eliminated. Therefore, a robust regression model and a Robust Kalman Filter model are proposed to solve this issue.

2.2. Robust Regression Soil Moisture Inversion Model

In thesis [33], the authors proposed a robust regression soil moisture inversion model (in Chinese). For convenience, we present now a brief description of this method in this study.
Massive studies and experimental campaigns show that the soil moisture has a near-linear relation to the observables of GNSS SNR multipath oscillation. Regression is a statistical method to build the relationship between the independent variables and the dependent variables. Robust regression is a proper algorithm to suppress the effect of the environmental noise and the thermal noise induced by GNSS receiver and antenna.
The general form of regression model is shown as Equation (5) in Section 2.1. The residual vector v of the equation includes the constant bias of the linear regression function and also random noise.
The concept of robust regression is based on M-estimate by utilizing the Iterative Reweighed Least Squares (IRLS) method for regression coefficient estimating [34].
A function Q is defined as Equation (6):
Q = j = 1 N g ( v j )
where g ( v j ) is the Huber robust error function defined as Equation (7):
g ( v j ) = { 1 2 ( v j ) 2 i f | v j | c c | v j | 1 2 c 2 i f | v j | > c
where c is harmonic coefficient, which generally is an empirical value. According to article [35], when c is set to 1.345 the regression can achieve 95% efficiency with high robustness, and this value is also the default value of Matlab function robustfit, so in this study we also set c is to 1.345.
Then, the regression problem turns to an M optimization problem defined as Equation (8) for each frequency:
arg min H j N Q , s . t . v = Y β X
when Q is minimalized, Q / β = 0 , then Equation (9) is denoted:
j = 1 n f ( y j x j β j ) x j = 0
where f ( v j ) = d g ( v j ) / d v j , then Equation (10) is denoted:
f ( v j ) = { c i f   v j c v i f   v j c c i f   v j c
For achieving higher robustness, a scaled estimation value s j for each satellite is introduced to standardize the residuals. Following [34], we obtain Equation (11):
s j = m e d | v j m e d ( v j ) | 0.6745
where 0.6745 is median absolute deviation proposed by Hample [36] to guarantee the unbiased estimation under normal distribution, and m e d ( a ) is the for denoting the median absolute deviation. Therefore, we can normalize the residuals as Equation (12):
μ j = v j / s j
The weight of the observation on jth day can be defined as Equation (13):
W j = f ( μ j ) μ j W = d i a g [ W j ] X T W ( Y β X ) = 0
Then, we can obtain the iteration formulation of robust regression defined as Equation (14):
β ~ = ( X T W X ) 1 X T W Y
The main idea of the robust regression model is to assert different weights for different points depending on its residual—the smaller the residual of one point, the greater weight it will have. Then, the weights are optimized by iterating a weighted algorithm for 15 times and the first three observables are not used.

2.3. Robust Kalman Filter Soil Moisture Inversion Model

In our second model we developed a processing chain using a Kalman Filter. The idea of a Kalman Filter is to use recursion of input and output values to calculate and update by least mean square error estimation of the state. In this section we establish the state equation and observation equation. Then, we need to ameliorate the observation equation using the Huber-M estimation method. The standard Kalman Filter model assumes the true state at time j is evolved from the states at time j 1 [26], as shown in Equation (15). At time j , daily soil moisture y j of the true state x j ( a j , p j ) is made as Equation (16):
x j = F j x j 1 + ω j
y j = H j x j + ε j
where F j is the state transition matrix which is applied to the previous state x j 1 ; ω j is the process noise with a value of zero which is assumed to be drawn from a zero-mean multivariate Gaussian distribution χ , with covariance Ω j : ω j ~ χ ( 0 , Ω j ) .
H j is the observation model which maps the true state space into the observing space and ε j is the observation noise which is assumed to be zero-mean Gaussian white noise with covariance R j : ε j ~ χ ( 0 , R j ) .
Here is given a recursive process. We estimate the state at time j through the state at time j − 1, then calculate the error correlation matrix P and Kalman gain K , update the state variable, and then output y . The recursive algorithm of Kalman Filter is demonstrated as follows.
First, set both the F 1 and H 1 to be an identity matrix as the entrance of the recursive algorithm. Then, calculate each state x at time j predicted by state x at time j 1 defined as Equation (17):
x j | j 1 = F j x j 1 | j 1
Calculate the error covariance matrix that is between predicted and true values defined as Equation (18):
P j | j 1 = F j P j 1 | j 1 F j T + Ω j
According to Equation (18), Kalman gain can be computed as Equation (19):
K j = P j | j 1 H j T ( R j + H j P j | j 1 H j T ) 1
Equations (17) and (19) are used to calculate the estimated value of x at time j , defined as Equation (20):
x j | j = x j | j 1 + K j ( y j H j x j | j 1 )
Then, compute the error covariance matrix that is between estimated and true values defined as Equation (21):
P j | j = ( I K j H j ) P j | j 1 ( I K j H j ) T + K j R j K j T
where I is the identity matrix.
To obtain the optimal Kalman gain, the error covariance matrix further simplifies to Equation (22):
P j | j = ( I K j H j ) P j | j 1
Normally, Kalman Filter is more effective for suppressing Gaussian distributed noise. Otherwise, the robustness will be impacted. Therefore, we try to use the Huber-M estimation described in Section 2.2 to reconstitute the observables for improving the robustness. The relationship between state truth value x j and its predicted value x j | j 1 at time j is shown as Equation (23):
x j = x j | j 1 + δ j
where δ is the error of prediction. A linear regression model is constructed by combining Equation (16) as Equation (24):
[ y j x j | j 1 ] = [ H j I ] x j + [ ε j δ j ]
Then, we define Equation (25):
D j = [ R j 0 0 P j | j 1 ] z j = D j   1 / 2 [ y j x j | j 1 ] M j = D j   1 / 2 [ H j I ] γ j = D j   1 / 2 [ ε j δ j ]
Combined with Equation (2), we can define Equation (26):
z j = M j x j + γ j
By combining Equation (13), iteration x j is solved as described in Equation (27):
M j T W ( z j M j x j ) = 0
where M j T is the transpose of M j , W is the matrix of weight, and z j is the observation for true value x j   .
The result of the i + 1 th iteration as Equation (28):
x j i + 1 = ( M j T W M j ) 1 M j T W z j
At the end of the iteration, the variance is obtained as Equation (29):
P j | j = ( M j T W M j ) 1
Introducing Equations (28) and (29) above results into the updating of the observations in Equations (20) and (22), and the Huber-M estimation Robust Kalman Filter algorithm is computed.

3. Results

3.1. Experiment Campaign

To validate the proposed models, we chose to use collected ground thrust data from Lamasquère area (Southwest of France, see Figure 2). The GNSS dataset was collected at Lamasquère, France (43°29′14.45″N, 1°13′44.11″) from 5th of February 2014 to 15th of March 2014. The details of the experiment campaign were thoroughly demonstrated in article [22], so in this study we only brief the key information. The site is a soya field equipped with a Leica GR25 multifrequency receiver and an AR10 antenna. At the same time, two Theta Probe ML3 soil moisture sensors are installed to collect ground truth measurements. The sensors measured soil moisture every 2 min at 2-cm and 5-cm depths with few meters away from the receiver. The antenna height is about 1.7 m. The soil type is silt clay, which consists of 18% sand, 41% clay and 41% silt. During dry conditions L-band waves do not penetrate more than a few centimeters (maximum observed 5.2 cm). The soil temperature varies between 5 ℃ and 16 ℃. During the experiment campaign, the field had no vegetation cover and was considered as bare soil.
In this study, the SNR data used for the validation came only from low-elevation satellites with angles ranging from 2° to 30° of elevation. First, the data was processed by the method mentioned in Section 2.1 to obtain the observables such as frequency, h, amplitude, and phase of SNR multipath oscillation. The second step was solving Equation (5) for the classical model and Equation (14) for the robust regression model. Then, the conventional regression model, robust regression model and the Robust Kalman Filter model were applied on the observables. For making a thorough comparison of those three models, including the conventional model, robust regression model and Robust Kalman Filter model with in situ measurement corresponding to the ground thrust data used for validation, we applied them using three scenarios:
(1) Single-band univariate regression, by using only one observable.
(2) Dual-band data fusion univariate regression.
(3) Dual-band data fusion multivariate regression
The data processing flow is demonstrated by Figure 3.

3.2. Comparison of Results of Single-Band Univariate Regression on GPS L1—Scenario One

In this section, we inverted the soil moisture by conventional linear regression model, robust regression model and robust Kalman model by using one of two observables, such as amplitude and phase. The comparisons of the results are shown as follows.

3.2.1. Using Amplitude Observables

First, we demonstrated the relationship between in situ soil moisture and amplitude observables by the making of PRN1 as an example in Figure 4.
We computed the correlation coefficients between the observables of amplitude and in situ soil moisture on L1 and L2 bands, as shown in Figure 4, where RA1 is the correlation coefficients of amplitude on L1, and RA2 is the correlation coefficients of amplitude on L2. It is clear to see that there is a linear correlation between amplitude and in situ soil moisture. For this satellite the correlation on the L2 band is better than L1. After training and predicting, we obtained the inversion results of the amplitude observable on GPS L2. We show the soil moistures retrievals in Figure 5 of our time series (47 days) for satellite PRN1 using the amplitude observable as input of the inverting model.
We compared the results of the three models by using the amplitude observable with in situ soil moisture as our reference dataset. For the initialization step of all methods the first solution is coarse in comparison with the reference dataset. During days 5 until 20 one can see a gradual increase of the soil moisture for in situ data (from ~28*10−2 m3/m3 to ~29*10−2 m3/m3). As demonstrated in Figure 5, the results of the Robust Kalman Filter model show much better agreement to the increase of in situ soil moisture after day 10 than the other two models. The second increase (days 20 to 30) of the in situ soil moisture is clearly identified by the models (lower amplitudes, Figure 5). After day 35 the models differ more from real measurements. However, we can see clearly that the Robust Kalman Filter model improved the estimation results and is better than robust/classical regressions.

3.2.2. Using Phase Observables

Figure 6 demonstrates the relationship between observed phase and in situ soil moisture.
We computed the correlation coefficients between the observables of amplitude and in situ soil moisture on L1 and L2 bands, as Figure 6 shows, where RP1 is the correlation coefficients of amplitude on L1, RP2 is the correlation coefficients of amplitude on L2. According to Figure 6, when the in situ soil moisture increased, the observed phase on L1 and L2 showed a trend of decrease and vice versa, and both of them have a negative correlation. The inversion results of soil moisture using GPS L1 phase observable are shown in Figure 7.
We compared the soil moisture results of all the three models using phase with in situ soil moisture. The in situ time series changes were well reconstructed by the three models, therefore the observed data in the first 10 days are well transcribed by the three models. We then have a sudden decrease until day 13 for the two linear regression models which show significant differences with in situ soil moisture. Until day 23 the Robust Kalman model continue to decrease which is more consistent with in situ data. For the three models, at day 30 they are all well reconstructed, especially the robust regression model, and the main decrease from day 30 to 45 is only visible with Robust Kalman model. The regression models are strongly divergent from the field data after day 35. The Robust Kalman model improves more significantly than the robust regression and classical ones.
We now compare phase results to those using amplitude as input. The results of the phase give better correlation and lower difference with respect to in situ soil moisture. As with the amplitude, the most robust model remains the one that uses the Kalman Filter, as it clearly follows the validation time series and presents small deviations from in situ time series even after day 37. The regression models, whether classic or robust, do not show this strong final decrease; we can see, however, that in using the phase the deviations from the validation data is smaller than for amplitude.
In general, in the single-band univariate scenario, compared with conventional method and the robust regression model, the Robust Kalman model can improve the precision of soil moisture inversion, and the Robust Kalman Filter model achieves much higher performance than regression models.

3.3. Comparison between Dual-Band Data Fusion Classical Regression Model, Robust Regression Model and Robust Kalman Model—Scenario Two

In scenario two, we invert the soil moisture using dual-band data fusion methods. For the four observables, amplitude and phase of both L1 and L2 bands, the comparisons of the results are shown as follows.

3.3.1. Using Amplitude

Figure 8 shows the results using the amplitude observable on GPS L1 and L2.
We compared the soil moisture results of all the three models with dual-band using amplitude with in situ soil moisture. The in-situ time series changes are well reconstructed by the three models, therefore the observed data in the first 5 days are well transcribed by the Robust Kalman Filter models, while the other two regression models have a strange difference. We then have a gradual increase from day 5 to day 20 with in situ soil moisture. The Robust Kalman Filter model has the lowest difference with in situ data but has an oscillation which is not in agreement with the decreasing trend of in situ soil moisture from day 40 to 47, while the other two models have more difference. The second increase (days 20 to 30) of the in situ soil moisture is clearly identified by all the models, and after day 35 the two models (robust regression and conventional regression) strongly differ from real measurements. However, we can clearly see that the Robust Kalman Filter model improved the estimation results better than the robust/classical regressions.

3.3.2. Using Phase

The inversion results of dual-band data fusion methods by using amplitude observables on GPS L1 and L2 are demonstrated in Figure 9.
We compared the results of the three models by using the phase observable with in situ soil moisture as our reference dataset. The first 5 days for three models are more precise than the single-band scenario. During days 5 until 20 one can see a gradual increase of the soil moisture for in situ data (from ~28*10−2 m3/m3 to ~29*10−2 m3/m3). For the two models (robust regression and conventional regression) in general this growth has a dramatic reduction on day 13. The second increase (days 20 to 30) as compared to scenario one of the in situ soil moisture is clearly identified by the models. The robust regression model has a lower difference in general regarding this growth, and after day 30 the models strongly decrease. However, we can see clearly that the Robust Kalman Filter model improved the estimation results better than the robust/classical regressions. The two linear regression models differ with in situ soil moisture from day 35, but the Robust Kalman Filter model improved the estimation results and agreed with the decreasing trend.
In general, for the dual-band data fusion scenario, compared with the conventional method and robust regression method, the Robust Kalman model can improve the precision of soil moisture inversion when compared with the data fusion conventional method and robust regression method. Meanwhile, the data fusion Robust Kalman improves more obviously with all the observables than the robust regression method if comparing them in a single-band scenario.

3.4. Multivariate Variable Dual-Band Data Fusion—Scenario Three

In this section we invert the soil moisture by conventional linear regression model, robust regression model and Robust Kalman model by using all of the observables on both GPS L1 and L2 bands simultaneously. The comparison results are shown in Figure 10.
We compared the results of the three models by using multivariate variable dual-band with in situ soil moisture as our reference dataset. In the first 5 days, the Robust Kalman Filter model shows minimal differences with in situ data, and it is better than scenario one and two. The other two models have much more differences with in situ data. During days 5 until 20 one can see a gradual increase of the soil moisture for in situ data (from ~28*10−2 m3/m3 to ~29*10−2 m3/m3). For the Robust Kalman Filter model in general this growth is much faster and reaches a maximum intensity after 10 days. The second increase (days 20 to 30) of the in situ soil moisture is clearly identified by all the three models, and after day 35 models of conventional method and robust regression method differ from real measurements, but they are better than the other two scenarios. However, we can see clearly that the Robust Kalman Filter model improved the estimation results better than the robust/classical regressions, and it is also better than the previous scenarios.

4. Analysis and Discussion

It has been demonstrated that GNSS-IR has the ability to retrieve land surface parameters, and especially the observable of soil moisture. After K. Larson first retrieved soil moisture from GPS SNR data [14], massive studies have developed models for GNSS-IR soil moisture inversion. Although most studies focus on low elevation angles, in article [37] which was mentioned above, the authors take the pseudo-dynamic of the surface into account and obtain a significantly improved and also utilized signal with high elevation angles. How to apply the Robust Kalman Filter model with signal of high elevation angles is a main interest in future work.
In the experimental campaign of Lamasquère, the land is bare soil, but how to eliminate the impact of vegetation is also important. In article [38], the author proposed a multivariate adaptive regression spline method, considering the impact of the vegetation moisture content, with a correlation coefficient of 0.916 and a root-mean-squared error of 0.021 m3/m3. In article [39], the authors performed a 15-month observation which covered an entire growing cycle by two antennas and developed an inversion model on GPS L2C and L5 SNR, achieving a precision of 0.035 m3/m3 for the whole meadow growing cycle, and of 0.018 m3/m3 after grass cutting.
For every model, the correlation coefficients are calculated for each satellite on each band or dual-band combination as Figure 11. Satellites in which correlation coefficients are over 0.5 are considered as effective cases.
Compared with the conventional model, the robust regression model increases on average 28.77% for L1, 18.33% for L2 and 5% for dual-band. The correlation coefficient of the Robust Kalman Filter model increases on average 32.33% for L1, 28.14% for L2 and 19.10% for dual-band. As the statistical data shows, the Robust Kalman Filter model achieves the highest precision.
The correlation coefficients of the three models show that the Robust Kalman Filter model gives a better correlation than regression models, except for satellites PRN 6, 7, 8, 9, 19 and 24, for which the correlation coefficient is very weak, reflecting a non-correlation (R < 0.5). For demonstration, we counted the number of effective satellites for all the methods, as shown in Table 1.
There are 7 satellites for the robust regression model and 4 for conventional regression model on the L1 band, but the Robust Kalman Filter model has 13 satellites for soil moisture inversion. On the other hand, data fusion has also had a positive impact on the increase in the number of effective satellites. If we compare the results of the conventional method for dual-band and single-band, the correlation coefficient of the L1 band with an average increment is 48.66% and a maximum grow of 90.38%. Meanwhile, the correlation coefficient of L2 with an average increment is 26.93% and a maximum grow of 98.21%. For the robust regression model, the correlation coefficient of the L1 band with an average increment is 34.45% and a maximum grow of 80.36%. Meanwhile the correlation coefficient of L2 with an average increment is 14.06% and a maximum grow of 91.07%. The correlation coefficient of the Robust Kalman Filter model on the L1 band with an average increment is 34.80% and a maximum grow of 95.45%. Meanwhile, the correlation coefficient of L2 with an average increment is 18.50% and a maximum grow of 93.33%.
Then, we analyzed the correlation between the phase observable and in situ soil moisture as shown in Figure 12.
The correlation of phase is not better than amplitude and there are few satellites effective for soil moisture estimation. However, we can see that there are some satellites which are not effective for conventional methods but effective for robust regression and Robust Kalman Filter method. It shows both the robust regression and Robust Kalman Filter models are capable of increasing the number of effective satellites, such as PRN5 and PRN26 for the L1 band, and PRN3 and PRN14 for the L2 band. The figure also shows that the dual-band data fusion method makes improvement of correlation between the phase observable and in situ soil moisture.
Compared with the conventional model, the robust regression model increased by an average of 24.59% for L1, 35.22% for L2 and 24.63% for dual-band. The correlation coefficient of the Robust Kalman Filter model increased by an average of 33.96% on L1, 43.92% on L2 and 35.29% for dual-band. The Robust Kalman Filter can greatly improve the inversion accuracy of the model regardless of single frequency or dual-band fusion.
If we look the precision of the three models, i.e., correlation coefficients, one can see that the Robust Kalman Filter model gives a better correlation than regression models, except for satellites PRN 6, 9, 10, 12, 18, 19, 21, 22, 23 and 29, where the correlation coefficient is very weak and reflects a non-correlation (R < 0.5). For demonstration, we counted the number of effective satellites for all the methods, as shown Table 2.
The table shows that there are more effective satellites on the L1 band than L2 band and the dual-band data fusion also has a positive impact for increasing the number of effective satellites.
If we compare the results of the conventional method for dual-band and single-band, the correlation coefficient of the L1 band with an average increment is 25.68% and a maximum grow of 95.45%. Meanwhile, the correlation coefficient of L2 with an average increment is 36.46% for L2 and a maximum grow of 93.24%. For the robust regression model, the correlation coefficient of the L1 band with an average increment is 17.38% and a maximum grow of 93.84%. Meanwhile, the correlation coefficient of L2 with an average increment is 28.60% and a maximum grow of 98.75%. The correlation coefficient of the Robust Kalman Filter model on the L1 band with an average increment is 19.74% and a maximum grow of 98.61%. Meanwhile, the correlation coefficient of L2 with an average increment is 31.76% and a maximum grow of 97.67%.
Compared with the correlation coefficient of the amplitude observable, dual-frequency fusion has a better effect on the L2 band and a greater improvement. However, the amplitude is reversed.
For every model, the root-mean-squared error (RMSE) is also calculated for each satellite on each band or dual band combination. The results are shown in Figure 13 and Figure 14. All the discussion for RMSE is based on valid satellites defined in line 498~499.
First, we discuss the results with amplitude. Compared to the conventional method, the RMSE of the robust regression model decreases by an average of 2.54% on L1 and 2.81% on L2, meanwhile the RMSE of the Robust Kalman Filter model decreases by an average of 10.19% on L1 and 12.39% on L2. Over half of the satellites’ RMSE are between 1% m3/m3 and 2% m3/m3 for the Robust Kalman Filter model.
After dual-band data fusion, the RMSE of the conventional model decreases 12.10% more than L1 and 3.46% more than L2. The RMSE of robust regression model decreases 14.30% more than L1 and 0.72% more than L2. The RMSE of the Robust Kalman Filter model decreases 18.23% more than L1 and 7.28% more than L2. As the statistical data show, the Robust Kalman Filter model demonstrates the highest precision.
Moreover, we discuss the results with phase. Compared to the conventional method, the RMSE of the robust regression model decreases by an average of 1.32% on L1 and 6.11% on L2, meanwhile the RMSE of the Robust Kalman Filter model decreases by an average of 11.49% on L1 and 8.87% on L2. Over half of satellites’ RMSE are between 1% m3/m3 and 2% m3/m3 for the Robust Kalman Filter model.
After dual-band data fusion, the RMSE of conventional model decreases 5.45% more than L1 and 46.5% more than L2. The RMSE of the robust regression model decreases 2.08% more than L1 and 45.01% more than L2. The RMSE of the Robust Kalman Filter model decreases 1.24% more than L1 and 47.18% more than L2. As the statistical data show, the Robust Kalman Filter model demonstrates the highest precision.
In the end, we discuss dual-band multivariate data fusion models. The correlation coefficients and RMSE are shown in Figure 15.
As demonstrated in Figure 15, the correlation between the inverted and the in situ soil moisture has been significantly enhanced for all the models, regardless of single-band or dual-band data fusion approach. The conventional model has 22 effective satellites un-der the multivariate fusion scenario, while the robust regression and Robust Kalman Filter methods have 25 and 28 satellites, respectively, by which means we have more effective satellites to estimate soil moisture than univariate models.
Compared with the results of the univariate models, the correlation coefficient of dual-band data fusion multivariate conventional method self-improves for all the GPS satellites with an average increase of 44.97% for phase and an average increase of 29.43% for amplitude; meanwhile, the RMSE decreases by an average of 37.00% for amplitude. The correlation coefficient of dual-band data fusion multivariate robust regression model self-improves with an average increase of 29.43%, meanwhile the RMSE decreases for all satellites with an average reduction of 29.42% for amplitude. The correlation coefficient of dual-band data fusion multivariate Robust Kalman Filter self-improves with an average increase of 20.7%, meanwhile the RMSE decreases for all satellites with an average reduction of 31.42%.

5. Conclusions

The inversion models of GNSS-IR soil moisture remote sensing have been investigated in this paper. First, we introduced the conventional regression model and the robust regression model and proposed a Robust Kalman Filter model. Then, the validation of the models was carried out in three different scenarios, including the “single-band, single-observable” scenario, the “dual-band, single-observable” scenario and the “dual-band, dual-observables” scenario. In the single-observable scenarios, the results retrieved using the phase observable had better performance than those retrieved using the amplitude observable, and dual-band data fusion performed better than single-band retrieval. For the dual-band scenarios, the results obtained using dual observables achieved better performance than those retrieved using a single observable. Furthermore, the results provided by the proposed model achieved superior performance to those obtained by the other two models in each scenario. This indicated that the proposed model was available for the data fusion of both dual-band and dual-observable cases, which was significant to take full use of the multi observables provided by multi-band GNSS signals to improve the retrieval performance in GNSS-IR soil moisture remote sensing.

Author Contributions

Conceptualization, L.J., F.G. and L.Y.; methodology, L.Y. and D.Y.; software, L.J., W.Y., H.R. and X.H.; data collection, J.D.; validation, L.J., B.S., W.Y. and Y.L.; writing—original draft preparation, L.Y. and L.J.; writing—review and editing, T.X., X.H., J.D. and N.W.; project administration, L.Y.; funding acquisition, L.Y. and H.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number 31971781, and Zhejiang Provincial Basic Public Welfare Research Project Foundation of China, grant number LGN19D040001.

Acknowledgments

The authors would like to thank Nicolas Roussel for data collection and Mutian Han for the help of data processing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Munoz-Martin, J.; Onrubia, R.; Pascual, D.; Park, H.; Pablos, M.; Camps, A.; Rüdiger, C.; Walker, J.; Monerris, A. Single-Pass Soil Moisture Retrieval Using GNSS-R at L1 and L5 Bands: Results from Airborne Experiment. Remote Sens. 2021, 13, 797. [Google Scholar] [CrossRef]
  2. Choudhary, R.; Athira, P. Effect of root zone soil moisture on the SWAT model simulation of surface and subsurface hydrological fluxes. Environ. Earth Sci. 2021, 80, 1–16. [Google Scholar] [CrossRef]
  3. Xu, Y.; Dong, K.; Jiang, M.; Liu, Y.; He, L.; Wang, J.; Zhao, N.; Gao, Y. Soil moisture and species richness interactively affect multiple ecosystem functions in a microcosm experiment of simulated shrub encroached grasslands. Sci. Total Environ. 2022, 803, 149950. [Google Scholar] [CrossRef] [PubMed]
  4. Saadatabadi, A.R.; Izadi, N.; Karakani, E.G.; Fattahi, E.; Shamsipour, A.A. Investigating relationship between soil moisture, hydro-climatic parameters, vegetation, and climate change impacts in a semi-arid basin in Iran. Arab. J. Geosci. 2021, 14, 1–18. [Google Scholar] [CrossRef]
  5. Yang, Y.; Mao, Y.; Sun, B. Basic performance and future developments of BeiDou global navigation satellite system. Satell. Navig. 2020, 1, 1–8. [Google Scholar] [CrossRef] [Green Version]
  6. Hein, G.W. Status, perspectives and trends of satellite navigation. Satell. Navig. 2020, 1, 22. [Google Scholar] [CrossRef]
  7. Lu, J.; Guo, X.; Su, C. Global capabilities of BeiDou Navigation Satellite System. Satell. Navig. 2020, 1, 27. [Google Scholar] [CrossRef]
  8. Senyurek, V.; Lei, F.; Boyd, D.; Kurum, M.; Gurbuz, A.C.; Moorhead, R. Machine Learning-Based CYGNSS Soil Moisture Estimates over ISMN sites in CONUS. Remote Sens. 2020, 12, 1168. [Google Scholar] [CrossRef] [Green Version]
  9. Soulat, F. Sea state monitoring using coastal GNSS-R. Geophys. Res. Lett. 2004, 31, 133–147. [Google Scholar] [CrossRef] [Green Version]
  10. Carreno-Luengo, H.; Luzi, G.; Crosetto, M. Impact of the Elevation Angle on CYGNSS GNSS-R Bistatic Reflectivity as a Function of Effective Surface Roughness Over Land Surfaces. Remote Sens. 2018, 10, 1749. [Google Scholar] [CrossRef] [Green Version]
  11. Martin-Neira, M. A Passive Reflectometry and Interferometry System (PARIS): Application to ocean altimetry. ESA J. 1993, 17, 331–355. [Google Scholar]
  12. Zavorotny, V.U.; Voronovich, A.G. Bistatic GPS signal reflections at various polarizations from rough land surface with moisture content. Proc. IEEE/IGARSS 2000, 7, 2852–2854. [Google Scholar] [CrossRef]
  13. Katzberg, S.J.; Torres, O.; Grant, M.S.; Masters, D. Utilizing calibrated GPS reflected signals to estimate soil reflectivity and dielectric constant: Results from SMEX02. Remote Sens. Environ. 2006, 100, 17–28. [Google Scholar] [CrossRef] [Green Version]
  14. 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. 2007, 12, 173–177. [Google Scholar] [CrossRef]
  15. Zavorotny, V.U.; Larson, K.M.; Braun, J.; Small, E.E.; Gutmann, E.D.; Bilich, A. A Physical Model for GPS Multipath Caused by Land Reflections: Toward Bare Soil Moisture Retrievals. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2010, 3, 100–110. [Google Scholar] [CrossRef]
  16. Nievinski, F.G.; Larson, K.M. Forward modeling of GPS multipath for near-surface reflectometry and positioning applications. GPS Solut. 2013, 18, 309–322. [Google Scholar] [CrossRef]
  17. Nievinski, F.G.; Larson, K.M. An open source GPS multipath simulator in Matlab/Octave. GPS Solut. 2014, 18, 473–481. [Google Scholar] [CrossRef]
  18. Chew, C.; Small, E.E.; Larson, K.M. An algorithm for soil moisture estimation using GPS-interferometric reflectometry for bare and vegetated soil. GPS Solut. 2015, 20, 525–537. [Google Scholar] [CrossRef]
  19. Kim, H.; Lakshmi, V.; Kwon, Y.; Kumar, S.V. First attempt of global-scale assimilation of subdaily scale soil moisture estimates from CYGNSS and SMAP into a land surface model. Environ. Res. Lett. 2021, 16, 074041. [Google Scholar] [CrossRef]
  20. Chew, C.; Small, E. Estimating inundation extent using CYGNSS data: A conceptual modeling study. Remote Sens. Environ. 2020, 246, 111869. [Google Scholar] [CrossRef]
  21. Al-Khaldi, M.M.; Johnson, J.T.; O’Brien, A.J.; Balenzano, A.; Mattia, F. Time-Series Retrieval of Soil Moisture Using CYGNSS. IEEE Trans. Geosci. Remote Sens. 2019, 57, 4322–4331. [Google Scholar] [CrossRef]
  22. Pan, Y.; Ren, C.; Liang, Y.; Zhang, Z.; Shi, Y. Inversion of surface vegetation water content based on GNSS-IR and MODIS data fusion. Satell. Navig. 2020, 1, 21. [Google Scholar] [CrossRef]
  23. Han, M.; Zhu, Y.; Yang, D.; Hong, X.; Song, S. A Semi-Empirical SNR Model for Soil Moisture Retrieval Using GNSS SNR Data. Remote Sens. 2018, 10, 280. [Google Scholar] [CrossRef] [Green Version]
  24. Yang, T.; Wan, W.; Chen, X.; Chu, T.; Hong, Y. Using BDS SNR Observations to Measure Near-Surface Soil Moisture Fluctuations: Results from Low Vegetated Surface. IEEE Geosci. Remote Sens. Lett. 2017, 14, 1308–1312. [Google Scholar] [CrossRef]
  25. Kalman, R.E. A New Approach to Linear Filtering and Prediction Problems. J. Basic Eng. 1960, 82, 35–45. [Google Scholar] [CrossRef] [Green Version]
  26. O’Driscoll, C.; Petovello, M.G.; Lachapelle, G. Choosing the coherent integration time for Kalman Filter-based carrier-phase tracking of GNSS signals. GPS Solut. 2011, 15, 345–356. [Google Scholar] [CrossRef]
  27. Pham, T.M. Kalman Filter mechanization for INS airstart. IEEE AES Syst. Mag. 1992, 7, 3–11. [Google Scholar] [CrossRef]
  28. Ambadan, J.T.; Tang, Y. Sigma-Point Kalman Filter Data Assimilation Methods for Strongly Nonlinear Systems. J. Atmos. Sci. 2009, 66, 261–285. [Google Scholar] [CrossRef]
  29. Strandberg, J.; Hobiger, T.; Haas, R. Real-time sea-level monitoring using Kalman Filtering of GNSS-R data. GPS Solut. 2019, 23, 61. [Google Scholar] [CrossRef] [Green Version]
  30. Li, Z.; Zhang, B.; Yang, D.; Zhang, G. Kalman Filter Model for Spaceborne GNSS-R Ocean Wind Retrieval. J. Navig. Position. 2020, 8, 37–38. [Google Scholar] [CrossRef]
  31. Huang, F.; Garrison, J.L.; Rodriguez-Alvarez, N.; O’Brien, A.J.; Schoenfeldt, K.M.; Ho, S.C.; Zhang, H. Sequential Processing of GNSS-R Delay-Doppler Maps to Estimate the Ocean Surface Wind Field. IEEE Trans. Geosci. Remote Sens. 2019, 57, 10202–10217. [Google Scholar] [CrossRef]
  32. Wang, J.; Xu, T.; Wang, N.; He, Y.; Gao, F. Research on Sea Surface Height Measurement Based on GNSS-IR Dual Frequency Data Fusion. China Satellite Navigation Conference (CSNC), Chengdu, China, 22–25 November 2020; pp. 153–165. [Google Scholar]
  33. Jing, L. Surface Soil Moisture Estimated from Dual-Frequency GNSS Interferometric Reflectometry; Shandong Agricultural University: Tai’an, China, 2019. [Google Scholar]
  34. Huber, P.J.; Ronchetti, E.M. Robust Statistics. Wiley Series in Probability and Statistics, 2nd ed.; John Wiley &Sons: Hoboken, NJ, USA, 2009. [Google Scholar] [CrossRef]
  35. Huber, P. Robust Statistics; John Wiley: New York, NY, USA, 1981. [Google Scholar]
  36. Hample, F.R. Robust Statistics: The Approach Based on Influence Function; John Wiley &Sons: New York, NY, USA, 1986. [Google Scholar]
  37. Roussel, N.; Frappart, F.; Ramillien, G.; Darrozes, J.; Baup, F.; Lestarquit, L.; Ha, M.C. Detection of Soil Moisture Variations Using GPS and GLONASS SNR Data for Elevation Angles Ranging From 2° to 70°. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 4781–4794. [Google Scholar] [CrossRef]
  38. 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. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 13, 2442. [Google Scholar] [CrossRef]
  39. Zhang, S.; Calvet, J.-C.; Darrozes, J.; Roussel, N.; Frappart, F.; Bouhours, G. Deriving surface soil moisture from reflected GNSS signal observations from a grassland site in southwestern France. Hydrol. Earth Syst. Sci. 2018, 22, 1931–1946. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Scenario of Interference Generating of GNSS-IR, where θ is the elevation angle of the considered GNSS satellite, blue line represents the antenna gain pattern, and h is the antenna effective height. DS is the direct signal and RS the reflected one.
Figure 1. Scenario of Interference Generating of GNSS-IR, where θ is the elevation angle of the considered GNSS satellite, blue line represents the antenna gain pattern, and h is the antenna effective height. DS is the direct signal and RS the reflected one.
Remotesensing 13 04013 g001
Figure 2. Lamasquère Experiment location: (a) map of France with a black box locating the region surrounding Lamasquère; (b) zoom on the Lamaquère plots where the measurements were taken (yellow dot); (c) zoom on the installed instrumentation with 2 ML3 Theta probes and Leica GR25 receiver and AR10 antenna.
Figure 2. Lamasquère Experiment location: (a) map of France with a black box locating the region surrounding Lamasquère; (b) zoom on the Lamaquère plots where the measurements were taken (yellow dot); (c) zoom on the installed instrumentation with 2 ML3 Theta probes and Leica GR25 receiver and AR10 antenna.
Remotesensing 13 04013 g002
Figure 3. Data processing workflow: the in situ measurements of soil moisture are used as reference and are compared to the robust regression model and the Robust Kalman Filter model.
Figure 3. Data processing workflow: the in situ measurements of soil moisture are used as reference and are compared to the robust regression model and the Robust Kalman Filter model.
Remotesensing 13 04013 g003
Figure 4. Normalization of in situ soil moisture and observables (amplitude of PRN1).
Figure 4. Normalization of in situ soil moisture and observables (amplitude of PRN1).
Remotesensing 13 04013 g004
Figure 5. Estimation results of methods by using amplitude on L2 (PRN1), we highlighted for validation data the main moisture peaks (P1 day 17; P2 day 29) and holes (H1 day 6; H2 day 20; and H3 day 47). We can note that the difference between H1 and P1 is less than 1*10−2 m3/m3 and between H2 and P2 it is ~1*10−2 m3/m3 which is low and will allow us to test the precision of the inversions and to show the capacity of this technique to differentiate between very weak soil moisture fluctuations.
Figure 5. Estimation results of methods by using amplitude on L2 (PRN1), we highlighted for validation data the main moisture peaks (P1 day 17; P2 day 29) and holes (H1 day 6; H2 day 20; and H3 day 47). We can note that the difference between H1 and P1 is less than 1*10−2 m3/m3 and between H2 and P2 it is ~1*10−2 m3/m3 which is low and will allow us to test the precision of the inversions and to show the capacity of this technique to differentiate between very weak soil moisture fluctuations.
Remotesensing 13 04013 g005
Figure 6. Normalization of in situ soil moisture and observables (phase on L1 for PRN1).
Figure 6. Normalization of in situ soil moisture and observables (phase on L1 for PRN1).
Remotesensing 13 04013 g006
Figure 7. Estimation results of methods by using phase on L1 (PRN1), we highlighted for validation data the main moisture peaks (P1 day 17; P2 day 29) and holes (H1 day 6; H2 day 20; and H3 day 47). We can note that the difference between H1 and P1 is less than 1*10−2 m3/m3 and between H2 and P2 it is ~1*10−2 m3/m3 which is low and will allow us to test the precision of the inversions and to show the capacity of this technique to differentiate between very weak soil moisture fluctuations.
Figure 7. Estimation results of methods by using phase on L1 (PRN1), we highlighted for validation data the main moisture peaks (P1 day 17; P2 day 29) and holes (H1 day 6; H2 day 20; and H3 day 47). We can note that the difference between H1 and P1 is less than 1*10−2 m3/m3 and between H2 and P2 it is ~1*10−2 m3/m3 which is low and will allow us to test the precision of the inversions and to show the capacity of this technique to differentiate between very weak soil moisture fluctuations.
Remotesensing 13 04013 g007
Figure 8. Estimation results of dual-band data fusion methods by using amplitude (PRN1), we highlighted for validation data the main moisture peaks (P1 day 17; P2 day 29) and holes (H1 day 6; H2 day 20; and H3 day 47). We can note that the difference between H1 and P1 is less than 1*10−2 m3/m3 and between H2 and P2 it is ~1*10−2 m3/m3 which is low and will allow us to test the precision of the inversions and to show the capacity of this technique to differentiate between very weak soil moisture fluctuations.
Figure 8. Estimation results of dual-band data fusion methods by using amplitude (PRN1), we highlighted for validation data the main moisture peaks (P1 day 17; P2 day 29) and holes (H1 day 6; H2 day 20; and H3 day 47). We can note that the difference between H1 and P1 is less than 1*10−2 m3/m3 and between H2 and P2 it is ~1*10−2 m3/m3 which is low and will allow us to test the precision of the inversions and to show the capacity of this technique to differentiate between very weak soil moisture fluctuations.
Remotesensing 13 04013 g008
Figure 9. Estimation results of dual-band data fusion methods by using phase, we highlighted for validation data the main moisture peaks (P1 day 17; P2 day 29) and holes (H1 day 6; H2 day 20; and H3 day 47). We can note that the difference between H1 and P1 is less than 1*10−2 m3/m3 and between H2 and P2 it is ~1*10−2 m3/m3 which is low and will allow us to test the precision of the inversions and to show the capacity of this technique to differentiate between very weak soil moisture fluctuations.
Figure 9. Estimation results of dual-band data fusion methods by using phase, we highlighted for validation data the main moisture peaks (P1 day 17; P2 day 29) and holes (H1 day 6; H2 day 20; and H3 day 47). We can note that the difference between H1 and P1 is less than 1*10−2 m3/m3 and between H2 and P2 it is ~1*10−2 m3/m3 which is low and will allow us to test the precision of the inversions and to show the capacity of this technique to differentiate between very weak soil moisture fluctuations.
Remotesensing 13 04013 g009
Figure 10. Estimation results of multivariate variable dual-band data fusion methods (PRN1), we highlighted for validation data the main moisture peaks (P1 day 17; P2 day 29) and holes (H1 day 6; H2 day 20; and H3 day 47). We can note that the difference between H1 and P1 is less than 1*10−2 m3/m3 and between H2 and P2 it is ~1*10−2 m3/m3 which is low and will allow us to test the precision of the inversions and to show the capacity of this technique to differentiate between very weak soil moisture fluctuations.
Figure 10. Estimation results of multivariate variable dual-band data fusion methods (PRN1), we highlighted for validation data the main moisture peaks (P1 day 17; P2 day 29) and holes (H1 day 6; H2 day 20; and H3 day 47). We can note that the difference between H1 and P1 is less than 1*10−2 m3/m3 and between H2 and P2 it is ~1*10−2 m3/m3 which is low and will allow us to test the precision of the inversions and to show the capacity of this technique to differentiate between very weak soil moisture fluctuations.
Remotesensing 13 04013 g010
Figure 11. Correlation coefficients of three methods on each frequency with amplitude.
Figure 11. Correlation coefficients of three methods on each frequency with amplitude.
Remotesensing 13 04013 g011
Figure 12. Correlation coefficients of three methods on each frequency with phase.
Figure 12. Correlation coefficients of three methods on each frequency with phase.
Remotesensing 13 04013 g012
Figure 13. RMSE of three methods on each frequency with amplitude.
Figure 13. RMSE of three methods on each frequency with amplitude.
Remotesensing 13 04013 g013
Figure 14. RMSE of three methods on each frequency with phase.
Figure 14. RMSE of three methods on each frequency with phase.
Remotesensing 13 04013 g014
Figure 15. Correlation coefficients and RMSE of dual-band with multivariate.
Figure 15. Correlation coefficients and RMSE of dual-band with multivariate.
Remotesensing 13 04013 g015
Table 1. Number of effective satellites used amplitude on each band.
Table 1. Number of effective satellites used amplitude on each band.
MethodL1L2Dual-Band
Conventional model 41419
Robust regression model 71619
Robust Kalman model131824
Table 2. Number of effective satellites used phase on each band.
Table 2. Number of effective satellites used phase on each band.
MethodL1L2Dual-Band
Conventional model 9412
Robust regression model 12614
Robust Kalman model141218
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jing, L.; Yang, L.; Yang, W.; Xu, T.; Gao, F.; Lu, Y.; Sun, B.; Yang, D.; Hong, X.; Wang, N.; et al. Robust Kalman Filter Soil Moisture Inversion Model Using GPS SNR Data—A Dual-Band Data Fusion Approach. Remote Sens. 2021, 13, 4013. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13194013

AMA Style

Jing L, Yang L, Yang W, Xu T, Gao F, Lu Y, Sun B, Yang D, Hong X, Wang N, et al. Robust Kalman Filter Soil Moisture Inversion Model Using GPS SNR Data—A Dual-Band Data Fusion Approach. Remote Sensing. 2021; 13(19):4013. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13194013

Chicago/Turabian Style

Jing, Lili, Lei Yang, Wentao Yang, Tianhe Xu, Fan Gao, Yilin Lu, Bo Sun, Dongkai Yang, Xuebao Hong, Nazi Wang, and et al. 2021. "Robust Kalman Filter Soil Moisture Inversion Model Using GPS SNR Data—A Dual-Band Data Fusion Approach" Remote Sensing 13, no. 19: 4013. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13194013

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