ORIGINAL RESEARCH article

Front. Earth Sci., 23 April 2021

Sec. Environmental Informatics and Remote Sensing

Volume 9 - 2021 | https://doi.org/10.3389/feart.2021.657206

Using Generalized Regression Neural Network to Retrieve Bare Surface Soil Moisture From Radarsat-2 Backscatter Observations, Regardless of Roughness Effect

  • 1. Geomathematics Key Laboratory of Sichuan Province, Chengdu Technological University, Chengdu, China

  • 2. College of Water Conservancy and Civil Engineering, Inner Mongolia Agricultural University, Huhhot, China

  • 3. Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing, China

Abstract

The combined influence of surface soil moisture and roughness on radar backscatters has been limiting SAR’s application in soil moisture retrieval. In the past research, multi-temporal analysis and artificial neural network (ANN) inversion of physically based forward models were regarded as promising methods to decouple that combined influence. However, the former does not consider soil roughness change over a relatively longer period and the latter makes it hard to thoroughly eliminate the effect of soil roughness. This study proposes to use generalized regression neural network (GRNN) to derive bare surface soil moisture (BSSM) from radar backscatter observations regardless of the effect of soil roughness (GRNN inversion of backscatter observations). This method not only can derive BSSM from radar backscatters, provided soil roughness is unknown in any long period, but also can train models based on small-size sample data so as to reduce the manual error of training data created by simulation of physically based models. The comparison of validations between BSSM-backscatter models and BSSM-roughness-backscatter models both analyzed by GRNN shows that the incorporation of soil roughness cannot raise the prediction accuracy of models and, instead, even reduce it, indicating that the combined influence is thoroughly decoupled when being analyzed by GRNN. Moreover, BSSM-backscatter models by GRNN are recommended due to their good prediction, even compared to those related models in past publications.

Introduction

The verified combination of influence of bare surface soil moisture (BSSM) and roughness on backscattering coefficients of synthetic aperture radar (SAR) especially at the C-band, under given radar parameters and homogeneous soil textures, has been limiting its application as an operational source of soil moisture in hydrology, though it is regarded as the most suitable for monitoring surface soil moisture (SSM), due to its high sensitivity to water contents, its high spatial resolution in the order of tens of meters for the distributed soil moisture, and its ability to neglect the influence of the atmosphere (; ; ; ; ; ; ). How to minimize the auxiliary parameters so as to ably construct a direct relationship between backscattering coefficients and SSM is the key to break that limitation of application.

In past researches, multi-temporal analysis is regarded as the promising one to decouple SSM effect on radar backscatter from the effects of other surface soil parameters (e.g., soil roughness) (; ). That methodology is under the assumption that soil roughness is constant within a sufficiently short time interval or soil roughness changes over a longer time scale compared to BSSM change, and thus SAR backscatter change is just related to BSSM change during this period, under the given radar configurations and homogeneous soil textures. However, this method ignores cultivation practices during this period that may produce soil roughness change (; ), and yet have to choose the SAR sensors with a short repeat cycle so as to guarantee constant roughness within this period.

Artificial neural network (ANN) inversion of physically based forward electromagnetic models has also been realized ably to retrieve BSSM, provided that surface roughness is an unknown parameter in the training (e.g., ; ; , ; ). These physically based models account for the interactions between the microwave radiation and soil [e.g., the integral equation model (IEM) by and the advanced integral equation model (AIEM) by ] and thus ably simulate backscattering coefficients in terms of soil attributes (e.g., the dielectric constant and the surface roughness). However, those researches utilized the feedforward ANNs, which require a huge amount of datasets in the training to guarantee the precision and robustness of models, and thus had to create the simulated data by a physically based model to train ANNs, because the small size of sample data hardly satisfied that requirement of training. The proven obvious discrepancy between simulated and measured backscatters () is brought into the course of training in which ANN leads to some uncertainty of the retrieval model. For example, showed the resultant standard error of estimation (SEE) of 2.75% and the determination coefficient (R2) of 0.85 using ANN inversion of IEMs, compared to a SEE of 2.16% and an R2 of 0.91 by incorporating roughness to ANN inversion of IEMs. This contrast indicates that in spite of the good prediction accuracy for SSM retrieval by ANN inversion of physical-based forward models, this method does not thoroughly decouple the effect of BSSM on backscatter from that of soil roughness but just reduces the effect of soil roughness on BSSM retrieval from backscatter.

Furthermore, if BSSM sample data and their corresponding SAR’s backscatter data are used in ANNs’ training so as not to produce the manual error from simulated BSSM data by physically based models, it will probably explain how ANN techniques work on separating the combined effect of BSSM and soil roughness on backscatter. GRNN, involving a single-pass learning and without backpropagation so as to reduce times of iterations, has a high-accuracy estimation due to its Gaussian kernel and can handle noises in the inputs. More important is that GRNN can train data in a multidimensional measurement space even with sparse datasets (). Recently, GRNN is increasingly exploited in the field of estimating soil moisture, due to its proven outperformance over backpropagation neural network (BPNN). For example, trained GRNN on sparse ground-based measurements by passive microwave, and trained GRNN using 84 ground-station measurements and their corresponding multi-source remote sensing data (optical and microwave).

The study uses a local case in China to construct backscatter-BSSM and backscatter-roughness-BSSM relationships by GRNN based on 147 datasets of samples. The 10-fold cross-validations of backscatter-BSSM and backscatter-roughness-BSSM models by GRNN are compared to show the null or even negative effect of soil roughness on BSSM retrieval from backscatters here. Moreover, as a contrast, this study uses parametric regression method to analyze backscatter-BSSM and backscatter-roughness-BSSM relationships, whose 10-fold cross-validations are also compared to confirm the combined effect of BSSM and soil roughness on backscatter. The contrast shows that soil roughness works totally differently when using nonparametric analysis of GRNN and parametric analysis of traditional regression.

Materials and Methods

The study area is located in Hetao Irrigation Plain, Inner Mongolia, China (Figure 1), covering a region of 457.7 km2, with a dry and cold climate and an annual average temperature of 6.3–7.7°C. The average annual precipitation and evaporation are approximately 139.4 and 2,070 mm, respectively.

FIGURE 1

Radarsat-2 Observations

A Radarsat-2 scene was acquired using fine quad-polarization mode and a mean incidence angle of 31.32°, with SAR sensor of C-band (∼5.33 GHz), multiple imaging modes, and different polarization channels. The acquisition date of this scene is on April 11, 2016, when the surface of the study region was bare, immediately prior to spring irrigation. The spatial resolution of the image was 8 m, and the coverage was 25 km × 25 km. Standard processing to extract full-polarization backscatter coefficients was applied to the Radarsat-2 image with ENVI 5.3.1 software, including multilook, single image filtering, geocoding, and radiometric calibration. Terrain correction was not included in pre-processing due to the flat terrain in the study area, and the resultant principal parameters of backscattering coefficients were georeferenced using multiple ground control points. Resultantly, two co-polarization backscatters (transmitting and receiving directions are the same, e.g., hh or vv) and two cross-polarization backscatters (transmitting and receiving directions are not the same, e.g., hv or vh) are described in Table 1 and Supplementary Table 1.

TABLE 1

ParameterRangeMean
Backscattering coefficients (dB)hh-polarization −15.20 to (−4.84)−9.56
vv-polarization −15.32 to (−5.92)−9.80
hv-polarization −28.72 to (−13.50)−21.85
vh-polarization −30.74 to (−13.82)−21.94
Radar configurationsFrequency (GHz)5.335.33
Mean incidence angle (°)31.3231.32

Radarsat-2 data description.

Ground Measurements

Ground measurements of BSSM in the study area were conducted simultaneously to Radarsat-2 imagery acquisition, approximately on April 11–14, 2016. A total of 147 sampling plots were collected (Figure 1) on the bare tillable surface. There was neither precipitation nor agricultural practices during the whole campaign. A portable global positioning system device was used to record sampling locations.

Soil samples were collected at three evenly distributed points at a depth of 0–10 cm in each plot under the temperature of 10.4°C and later analyzed for gravimetric moisture content and bulk density measurements using the oven drying method. Surface roughness parameters [root-mean-square (RMS) height and correlation length] were measured on location in each sampling plot using the meshboard method and were alternatively taken in south–north and east–west directions. The meshboard used is 1 m long, and the sampling interval is 1 cm. As a result, the measured volumetric soil moisture contents and roughness parameters are described in Table 2 and Supplementary Table 1.

TABLE 2

ParameterRangeMeanUnit
RMS height (S)0.26–1.680.78(cm)
Correlation length (L)1.87–70.8223.38(cm)
BSSM10–2618.00(vol.%)
Percentage of clay8.278.27(vol.%)
Percentage of sand40.4440.44(vol.%)
Percentage of silt51.2951.29(vol.%)
Temperature10.410.4(°C)

Description of the collected of soil parameters from 147 sampling plots.

Methodology

GRNN Algorithm

Compared to standard feedforward neural network (e.g., ; , ; ; ), the structure of GRNN is relatively simple and static with four layers: input, pattern, summation, and output layers (). Once the input goes through each unit in the pattern layer, the relationship between the input and the response would be stored in the unit. Therefore, the regression performed by GRNN is actually to output the most probable scalar , given specified input vector X:

where D_i is the distance between the input vector X and the training sample input vector Xi stored in the pattern unit (or cluster centre), and σ, called “spread” or “smoothing parameter,” is the only unknown parameter in the network and affects the fitness in GRNN architecture that needs optimization (). The spread represents the standard deviation of the input vector X from the training sample input vector Xi in the Gaussian kernel of , and GRNN can produce the optimal results as long as the spread is no more than the standard deviation (). Therefore, the optimized σ is theoretically not unique, and it was often suggested to use the cross-validation to estimate the optimized spread (; ). In this study, the GRNN package integrated in the MATLAB software is used.

There are 15 GRNN architectures due to 15 different kinds of input, respectively, for BSSM-backscatter and BSSM-roughness-backscatter relationships (Figure 2).

FIGURE 2

Model Evaluation

In this study, a 10-fold cross-validation technique is chosen to test the models’ predictive capabilities due to a limited number of soil samples (; ; ). The mean absolute error (MAE) and the root-mean-square-error (RMSE) () are used as the metrics to assess the performances of GRNN algorithms for soil moisture retrieval in terms of prediction accuracy. These metrics mentioned above are defined as follows:

where is the predicted bare surface soil moisture content of sample i, BSSMi is the in situ measurement of sample i, is the mean value of soil moisture contents of all in situ measurements, and n is the number of all samples.

The coefficient of determination was often used in different versions of mathematical definition in past publications (), which makes it difficult to compare. Moreover, the value of R2 is also related to the number of observations: for example, R2 could be as high as 0.97 in fitting a linear relationship between an independent variable X and a normally distributed dependent variable Y based on three observations, even though X is unrelated to Y, but with 100 observations, the R2 of 0.07 is enough to establish statistical significance at the 1% level (). Therefore, the statistical significance that is reflected by the value of R2 in different models with different observations is also difficult to compare. Instead, the significance level of F test () is used here as a complementary reference to the MAE and RMSE metrics, and the statistical significance package integrated in the MATLAB software is used.

Results

Fifteen BSSM-backscatter models are constructed by GRNN here, based on 15 combinations of full-polarized backscattering coefficients: { }, { }, {}, {}, {, }, {, }, {, }, {, }, {, }, {, }, {, , }, {, , }, {, , }, {, , }, and {, , , }. The 10-fold cross-validations are used to determine the appropriate spread parameters of these models. The empirical knowledge in the past publications indicates that the spread σ definitely was more than zero and it was probably to produce a good result by a value between 0.01 and 1 (), so here the range of the spread σ is explored from 0.01 to 1. This study utilizes MAE and RMSE to estimate the models, and utilize the significance level (P value) to confine the MAE and RMSE. The MAE and RMSE, which have P values of above 5%, will be eliminated, and thus those ranges not circled in a green dashed line in Figure 3 are eliminated first.

FIGURE 3

Furthermore, the independent test dataset is used to assess the final models’ performance and the metrics are shown in Figure 4. Based on the results of both 10-fold cross-validations and independent-test metrics, those ranges of the spread parameters in translucent blue are recommended to use for GRNN models. Table 3 lists the recommended backscatter-BSSM models, their respective appropriate spread ranges, and their corresponding MAE and RMSE values. Generally, the nine recommended backscatter-BSSM models have excellent MAE and RMSE values ranging, respectively, from 1.43 to 1.84 and 1.83 to 2.4, which are even better than the related models in past publications.

FIGURE 4

TABLE 3

Backscatter combinationRecommended the spreadMAERMSE
0.1–11.43–1.601.95–2.06
0.1–165–1.991.98–2.4
0.8∼11.51–1.551.83–1.87
(., .)0.4–11.69–2.112.01–2.47
(, )0.7–11.56–1.691.97–2.11
(, )0.9–11.56–1.602.1
(, )0.8–11.58–1.611.95–1.98
(, , )0.9–11.72–1.782.10∼2.15
(, , )0.7–11.63–1.842.08–2.28

The recommended backscatter models, spread parameters, and their performance metrics.

Comparisons and Discussion

Comparison Between Backscatter-BSSM and Backscatter-Roughness-BSSM Models by GRNN

Comparisons are made between backscatter-BSSM and backscatter-roughness-BSSM models constructed by GRNN in order to be clear whether soil roughness is thoroughly decoupled from the BSSM effect on SAR backscatters (Figure 5). Here, the MAE and RMSE are compared in terms of trend lines fitted by power function, where the solid lines denote trend lines of 15 backscatter-BSSM models and the dashed lines denote those in Figure 5, which are, respectively, in blue and magenta. The translucent ranges in Figure 5 are the recommended spread parameters. Generally, the solid lines are lower than their corresponding dashed lines for the 15 pairs of comparison, which implies an overall better prediction of 15 backscatter-BSSM models than 15 backscatter-roughness-BSSM models by GRNN. In the recommended ranges of spread parameters, backscatter-BSSM models have better prediction than backscatter-roughness-BSSM models.

FIGURE 5

Comparison Between Backscatter-BSSM and Backscatter-Roughness-BSSM Models by Parametric Regression

The comparison here is in order to confirm the combined influence of BSSM and surface roughness on SAR backscatters by parametric regression in this study. The backscatter-BSSM and backscatter-roughness-BSSM models by parametrical regression are defined based on the research (), which are written as follows:

In Figure 6, the red lines and markers symbolize the MAE and RMSE metrics of the 15 BSSM-backscatter regression analyses, whereas the blue lines and markers symbolize the MAE and RMSE metrics of the 15 BSSM-roughness-backscatter regression analyses. In general, the MAE and RMSE of the regression analyses with roughness for the 15 combinations of backscatters are lower than the corresponding values of the regression analyses without roughness, which means better prediction by considering the soil roughness effect when using the parametric analysis methods.

FIGURE 6

Discussion

The first comparison (Figure 5) shows that the incorporation of auxiliary roughness into BSSM retrieval from SAR backscatters by GRNN does not raise the prediction accuracy of models and even, most of the time, reduces the prediction accuracy, which indicates that the soil roughness is an invalid effect in BSSM retrieval from SAR backscatters by GRNN. The second comparison (Figure 6) shows that the incorporation of roughness into BSSM retrieval from SAR backscatters by parametric regression raises the prediction accuracy of models, which indicates the positive effect of soil roughness in BSSM retrieval from SAR backscatters by parametric regression. GRNN analysis signifies high-dimensional analysis, whereas parametric regression denotes low-dimensional analysis. The contrast between the two comparisons may be a hint that the relationships among BSSM, soil roughness, and backscattering coefficients are distinct in the analytic spaces with different dimensions.

Conclusion

This study derives BSSM directly from SAR backscatter observations using small-size sample data to train GRNN. This method is more promising to remove the effect of soil roughness on soil moisture retrieval from SAR backscatter than both multi-temporal analysis and feedforward ANN inversion of physically based forward models. Firstly, this method can avoid the inconvenience of temporal characteristics of roughness changes from multi-temporal analysis. Secondly, it minimizes the error of data to train ANNs using samples data, compared to feedforward ANN inversion that uses simulated data from physically based forward models. Furthermore, the models’ validations show that the incorporation of soil roughness into BSSM retrieval from backscatter by GRNN does not raise prediction accuracy and even reduces it (Figure 5), indicating thoroughly decoupling the effect of BSSM on backscatter from the effect of soil roughness. By contrast, as mentioned in Section 1, ANN inversion of IEMs (e.g., ) actually does not thoroughly decouple the combined effect of BSSM and soil roughness on SAR backscatter. However, the comparison of models’ validations between BSSM-backscatter and BSSM-roughness-backscatter models by parametric regression still shows the combined influence of soil roughness and BSSM on SAR backscatter. Conclusively, although ANNs are rather advantageous to recognize the superbly complex nonlinear relationships compared to traditional parametric analysis methods, the selection of ANNs is still crucial and it is better to choose the ANNs that can work well even with small-size data (e.g., GRNN).

In the 15 backscatter-BSSM models by GRNN, #x03C3;hh-BSSM, σvv-BSSM, σhv-BSSM, (, )-BSSM, (, )-BSSM, (, )-BSSM, (, )-BSSM, (, , )-BSSM, and (, , )-BSSM models are recommended (Figure 4 and Table 3), because they have small MAE (1.43 to 2.11 vol.%) and RMSE values (1.83 to 2.47 vol.%) under the spread range of statistical significance (5%). Especially recommended is the σhh-BSSM model whose MAE ranges from 1.43 to 1.60 vol.% and whose RMSE ranges from 1.95 to 2.06 vol.% over the 0.1–1 spread of statistical significance (5%).

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Author contributions

LZ and QL: conceptualization. LZ: formal analysis and methodology, software, validation, and writing—original draft. QL, LJ, and JF: funding acquisition. QL: investigation, project administration, and resources. QL and LJ: supervision. JF: visualization. LZ and LJ: writing—review and editing. All authors contributed to the article and approved the submitted version.

Funding

This work was funded by the National Science Foundation of China (No. 52069020, No. 1212011085468, No. 42002294, and No. 41972308), the Research Initiation Fund for Teacher Development from Chengdu University of Technology (No. 10912-2019KYQD07430), the State Key Research Project in 13th Five-Year (No. 51569018), the Science and Technology Major Project of Xinjiang Uygur Autonomous Region, and the China Postdoctoral Science Foundation (No. 2017M623307XB).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2021.657206/full#supplementary-material

References

Summary

Keywords

soil moisture, surface roughness, generalized regression neural network, SAR backscatter, high-dimensional analyses

Citation

Zeng L, Liu Q, Jing L, Lan L and Feng J (2021) Using Generalized Regression Neural Network to Retrieve Bare Surface Soil Moisture From Radarsat-2 Backscatter Observations, Regardless of Roughness Effect. Front. Earth Sci. 9:657206. doi: 10.3389/feart.2021.657206

Received

22 January 2021

Accepted

04 March 2021

Published

23 April 2021

Volume

9 - 2021

Edited by

Jie Dou, Nagaoka University of Technology, Japan

Reviewed by

Qiyu Chen, China University of Geosciences Wuhan, China; Xiaoning Song, University of Chinese Academy of Sciences, China

Updates

Copyright

*Correspondence: Quanming Liu,

This article was submitted to Environmental Informatics and Remote Sensing, a section of the journal Frontiers in Earth Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics