 Astronomy
 Biomedical Optics & Medical Imaging
 Defense & Security
 Electronic Imaging & Signal Processing
 Illumination & Displays
 Lasers & Sources
 Micro/Nano Lithography
 Nanotechnology
 Optical Design & Engineering
 Optoelectronics & Communications
 Remote Sensing
 Sensing & Measurement
 Solar & Alternative Energy
 Sign up for Newsroom EAlerts
 Information for:
Advertisers




Remote Sensing
Enhancing remotely sensed soil moisture estimations using in situ ancillary information
Matias Barber, Francisco Grings and Haydee Karszenbaum
Developing algorithms used to generate soil moisture maps of agricultural, bare soils from synthetic aperture radar systems benefits from including ancillary information within a Bayesian retrieval scheme.
7 February 2013, SPIE Newsroom. DOI: 10.1117/2.1201301.004710
Soil moisture is a key variable related to crop development and yield. Too much or too little moisture can each have devastating effects: persistently dry soil causes plants to wilt and diminishes their ability to transpire and grow. By contrast, excess moisture results in poor seed germination, inadequate nutrient uptake, and soil compaction. Consequently, the early assessment of soil moisture reserves, and monitoring moisture changes, is particularly important prior to—and throughout—the seeding season.
In response to this demand for information, a systematic effort has been made to develop maps of soil moisture over agricultural areas. Orbiting microwave synthetic aperture radar (SAR) systems offer the opportunity to monitor moisture content at different scales and under almost any weather condition. SAR systems leverage the known sensitivity that the backscattered signal (which originated in the sensor and was scattered by the target) exhibits to soil parameters, including soil moisture and soil roughness.^{1} However, the relationship between the backscattered signal and soil parameters is not straightforward, and consequently no operational SARderived soil moisture maps are yet available.
The major obstacles to developing such maps include the difficulty in modeling the scattering processes that relate backscattering to soil properties (moisture and roughness),^{2} the speckle noise,^{3} and the difficulty in measuring soil roughness in the field. The former two mainly relate to the SAR imaging system, whereas the latter relates to soil heterogeneity.^{4, 5} Moreover, several combinations of surface parameters can usually produce the same SAR observation. As a consequence, the retrieval of soil parameters remains challenging, and soilmoisture maps derived from remotely sensed SAR data remain fairly inaccurate.
Figure 1. Diagram showing the block elements comprising the Bayesian retrieval algorithm for soil parameters. SAR: Synthetic aperture radar. M: Model.
Many approaches have been explored to develop retrieval methodologies to infer soil condition, for example, changedetection procedures^{6–8} and radar backscatter modeling (theoretical and semiempirical).^{2, 9–11} Changedetection methods exploit the availability of temporal series of SAR acquisitions whereby variations in surface backscatter are expected to reflect changes in soil moisture, because other parameters affecting radar backscatter can be considered fairly constant. Radar backscatter modeling deals with the dependency of soil parameters to backscattered signal from a deterministic, physical point of view. By contrast, only Bayesian approaches^{12–14} are able to incorporate ancillary information about the soil in a systematic, natural way. We have been working on such a kind of Bayesian approach, whereby better performance is reached against a ‘blind’ (i.e., without ancillary data) estimation.
Roughly speaking, given a set of measured SAR backscattering coefficients, one tries to assign the best set of soil parameters (i.e., soil moisture and roughness) that account for those measurements. A Bayesian inference scheme for soil parameters is shown in Figure 1. Typically, the scheme comprises both the forward and error model (or ‘likelihood’ in Bayes' jargon), ancillary information (‘prior’) and output (or ‘posterior’). The likelihood measures the degree of compatibility between a certain SAR measurement and certain soil parameters constrained to some given forward model. The higher the values of the likelihood, the more likely that the SAR measurement comes from that specific combination of soil parameters. The prior involves all the information available about the soil parameters and can originate from historical or cadastral records, estimation from other sensors, in situ data, and/or contextual information. The pointbypoint multiplication of the likelihood by the prior leads to the posterior. For this reason, the posterior is a ‘modulated’ version of the likelihood function owing to the prior envelope. Finally, estimates of soil parameters are computed from the posterior in a number of ways (e.g., the mean or the maximum a posteriori).
Figure 2. Optical image from the study area. Plots are labeled as north (N) or south (S) followed by a site number (1–10). Image courtesy of the Argentine Microwave Observation Satellite (SAOCOM) Mission, Argentine Space Agency (CONAE).
Figure 3. Absolute difference between estimated soil moisture and in situ measurements using (a) our Bayesian approach and (b) the minimization approach. Dotted lines indicate a ‘desirable’ 0.05error bound (a measure of the desired maximum error). RMSE: Root mean squared error.
We recently modeled speckle noise and terrain heterogeneity within a Bayesian retrieval methodology.^{15} As a natural advantage of the Bayesian approach, prior information about soil condition can be easily included, enhancing the performance of the retrieval. Also, the results indicate that the Bayesian model enlarges the validity region of a standard minimizationbased procedure.
With this model in hand, we carried out an experimental study to assess the performance of the Bayesian approach against a traditional minimization procedure, estimates from which can be compared to in situ measurements. The study area consisted of a set of 20, 120×50m agricultural plots in Centro Espacial Teófilo Tabanera (CETT), near Falda del Carmen, Córdoba, Argentina (see Figure 2). The test site belongs to the Argentine Space Agency (CONAE). Tillage treatments are presented on the plots, leading to different roughness conditions and, therefore, to ‘expected’ roughness (prior information). Over the study site, 16 plots were notillage whereas only four were ploughed. Prior information for roughness was collected on a previous field campaign using a bidimensional laser profiler.^{16} The hydrological soil condition was dry, with a cumulative rainfall of 3.4mm in the 60day period preceding the fieldwork. Soil moisture was retrieved from a full polarimetric Lband (wavelength of 23cm) image from CONAE's airborne SAR system SARAT over the 20 plots.
The forward model employed in this study is a semiempirical one developed by Oh.^{11} After decorrelating the pixels to diminish the speckle noise, the mean values of the backscattered signal are computed on each plot and used as inputs to the retrieval algorithm. Thus, 20 soil moisture estimates (with their associated errors) can be computed and then compared with in situ data. Figure 3 shows the absolute difference between the estimated soil moisture and in situ measurements. Both retrieval approaches are shown: Bayesian—see Figure 3(a)—and minimization: see Figure 3(b). Correspondence between our Bayesian model and in situ measurements was achieved with a root mean squared error (RMSE) of 0.03cm^{3}/cm^{3} and a maximum absolute error (MAE) of 0.06cm^{3}/cm^{3}. For the minimization procedure, the RMSE was calculated at ∼0.04cm^{3}/cm^{3} with a MAE of 0.07cm^{3}/cm^{3}. Between the models, the Bayesian one performs better with less spread of the errors.
In summary, rendering quantitative soil moisture maps is of critical relevance for the management of the agricultural sector. Such maps are feasible using SAR images acquired by satellites or airborne systems. A Bayesian retrieval scheme can deal with the speckle noise introduced by the imaging system and the heterogeneity of the soil parameters. A remarkable feature of the Bayesian approach is that it includes information ancillary to the retrieval algorithm to enhance the output estimates. Our next step will be to assess the errors in the output estimates caused by the heterogeneity of the soil parameters.
The authors thank CONAE for providing the radar image and the staff at the Centro Espacial Teófilo Tabanera who supported the field campaign.
Matias Barber, Francisco Grings, Haydee Karszenbaum
Institute of Astronomy and Space Physics (IAFE), Argentina
Buenos Aires, Argentina
Matias Barber received his degree in physics from the University of Buenos Aires, Argentina, in 2009. He is currently pursuing a PhD, working on surface scattering models at the IAFE. He has been conducting field work since 2009.
References:
1. F. T. Ulaby, R. K. Moore, A. K. Fung, Microwave Remote Sensing: Active and passive. Radar Remote Sensing and Surface Scattering and Emission Theory, AddisonWesley, 1982.
2. A. K. Fung, Microwave Scattering and Emission Models and Their Applications, Artech House, 1994.
3. C. Oliver, S. Quegan, Understanding Synthetic Aperture Radar Images, SciTech Publishing, 2004.
4. A. W. Western, R. B. Grayson, G. Blöschl, G. R. Willgoose, T. A. McMahon, Observed spatial organization of soil moisture and its relation to terrain indices, Water Resour. Res. 353, p. 797810, 1999. doi:10.1029/1998WR900065
5. N. E. C. Verhoest, H. Lievens, W. Wagner, J. AlvarezMozos, M. S. Moran, F. Mattia, On the soil roughness parameterization problem in soil moisture retrieval of bare surfaces from synthetic aperture radar, Sensors 8(7), p. 42134248, 2008. doi:10.3390/s8074213
6. A. Balenzano, F. Mattia, G. Satalino, M. W. J. Davidson, Dense temporal series of C and Lband SAR data for soil moisture retrieval over agricultural crops, IEEE J. Sel. Topics Appl. Earth Observ. Remote Sens. 4(2), p. 439450, 2011. doi:10.1109/JSTARS.2010.2052916
7. F. Mattia, G. Satalino, V. R. N. Pauwels, A. Loew, Soil moisture retrieval through a merging of multitemporal Lband SAR data and hydrologic modelling, Hydrol. Earth Syst. Sci. 13, p. 343356, 2009. doi:10.5194/hess133432009
8. Y. Kim, J. J. van Zyl, A timeseries approach to estimate soil moisture using rolarimetric radar data, IEEE Trans. Geosci. Remote Sens. 47(8), p. 25192527, 2009. doi:10.1109/TGRS.2009.2014944
9. Y. Oh, K. Sarabandi, F. T. Ulaby, An empirical model and an inversion technique for radar scattering from bare soil surfaces, IEEE Trans. Geosci. Remote Sens. 30(2), p. 370381, 1992. doi:10.1109/36.134086
10. P. C. Dubois, J. van Zyl, T. Engman, Measuring soil moisture with imaging radars, IEEE Trans. Geosci. Remote Sens. 33(4), p. 915926, 1995. doi:10.1109/36.406677
11. Y. Oh, Quantitative retrieval of soil moisture content and surface roughness from multipolarized radar observations of bare soil surfaces, IEEE Trans. Geosci. Remote Sens. 42(3), p. 596601, 2004. doi:10.1109/TGRS.2003.821065
12. Z. S. Haddad, P. D. Dubois, J. J. van Zyl, Bayesian estimation of soil parameters from radar backscatter data, IEEE Trans. Geosci. Remote Sens. 34(1), p. 7682, 1996. doi:10.1109/36.481895
13. C. Notarnicola, F. Posa, Bayesian algorithm for the estimation of the dielectric constant from active and passive remotely sensed data, IEEE Trans. Geosci. Remote Sens. 1(3), p. 179183, 2004. doi:10.1109/LGRS.2004.827461
14. F. Mattia, G. Satalino, L. Dente, G. Pasquariello, Using a priori information to improve soil moisture retrieval from ENVISAT ASAR AP data in semiarid regions, IEEE Trans. Geosci. Remote Sens. 44(4), p. 900912, 2006. doi:10.1109/TGRS.2005.863483
15. M. Barber, F. Grings, P. Perna, M. Piscitelli, M. Maas, C. Bruscantini, J. JacoboBerlles, H. Karszenbaum, Speckle noise and soil heterogeneities as error sources in a Bayesian soil moisture retrieval scheme for SAR data, IEEE J. Sel. Topics Appl. Earth Observ. Remote Sens. 5(3), p. 942951, 2012. doi:10.1109/JSTARS.2012.2191266
16. M. Barber, F. Grings, P. Perna, C. Bruscantini, H. Karszenbaum, Field Campaign CETT 2011: analysis of roughness measurements, tech. rep., 2011.


