Seung-bum Kim , Shaowu Huang*, Leung Tsang*, Joel Johnson**, Eni Njoku

1. Soil moisture retrieval over bare surfaces using time-series radar observations and a lookup table representation of forward scattering Seung-bum Kim, Shaowu Huang*, Leung Tsang*, Joel Johnson**, Eni Njoku Jet Propulsion Lab., California Inst. Technology * Univ. Washington ** The Ohio State Univ. IGARSS 2011, Vancouver, Canada

2. A ‘bare and sparsely vegetation surface covers 13% of the world’s land surface. Over vegetated surfaces, a soil moisture signal comes from the bare surface. Accurate soil moisture retrieval over the bare surface forms the basis of global soil moisture retrieval We focus on the soil surface where the roughness has the isotropic random distribution. We study the retrieval of soil moisture at L-band at 40deg incidence angle to apply to the Soil Moisture Active Passive (SMAP) mission NASA’s SMAP mission Global, high-resolution mapping of soil moisture (top 5cm) and its freeze/thaw state Three year mission, due for launch in 2015 1000 km-wide swath, enabling 2−3 day revisit Retrievals with radiometer (36km resolution), SAR (3km), and SAR/radiometer combined (9km) Multi-pol (HH, VV, HV) To achieve the global coverage, the number of single looks was compromised (60, worst case), leading to somewhat large speckle. Total radar measurement error ranges from 0.5dB (13%) to 0.7dB (17%) Objectives

3. Four parameters dominating the radar scattering from bare surface: surface roughness (s), soil moisture (Mv), correlation length of roughness (l), correlation function (F) Issues to resolve for accurate soil moisture retrieval With HH, VV, HV, we cannot determine all 4 unknowns Knowledge of correlation length and function is very inaccurate (50% error). Ambiguity: dry & rough soil σ0~ wet & smooth soil simga0  can lead to multiple solutions 13% to 17% radar measurement error Past literature Empirical retrieval model without considering l or F (Oh et al. 1992; Dubois et al. 1995) The exponential function fits data well (Shi et al. 1997; Mattia et al. 1997) Estimate s and l after assuming temporally static Using ancillary dry condition (Rahman et al. 2008) Using one time measurement (Joseph 2008) Statistical estimate (Verhoest et al. 2007) Using ancillary weather model (Mattia et al. 2009) Eliminate s during retrieval (Shi et al. 1997) Our goal: develop non-empirical and simple method that does not need ancillary information Issues with bare surface soil moisture retrieval

4. A forward model (Numerical Maxwell Model in 3-Dimension, NMM3D, Huang et al. 2010; 2011) is inverted using its lookup table representation. The NMM3D computes numerical solutions of Maxwell’s equations without approximate parameterizations or tuning. The NMM3D predictions compare well with in situ datasets representing wide ranges of roughness, soil moisture, and correlation length. Soil moisture retrieval method rmse=1.49dB (VV) 1.64dB (HH) Huang, Tsang, Njoku, Chen, TGRS, 2010 Results with correlation length/rms height (l/s) =4,7,10 is also available.

5. Retrieval with a lookup table • The retrieval based on a lookup table performs better. In comparison, the inversion of a sophisticated model (IEM) is not as accurate. Lookup table RMS Error = 0.0016 cm3/cm3 Courtesy: J.J. van Zyl Parameterized IEM based on 3 Mv and 3 ks simulations RMS Error = 0.0067 cm3/cm3

6. Two-polarization (HH and VV) are used for input. HV channel is set aside for vegetation information for future. Mv and roughness are retrieved. 2N independent input (HH and VV, N is the number of time-series) N+1 unknowns assuming roughness does not change in time. The retrieval of the soil moisture is accomplished by the least square minimization. Time-invariant roughness is estimated first. Then time-varying soil moisture is retrieved. Soil moisture retrieval method

7. The radar measurement noise (0.7dB, 17%) is modeled by a Gaussian random process. Errors in roughness estimate is smaller than 10%. Real part of the dielectric constant (εr) is retrieved first and Mv error is smaller than 0.06 cm3/cm3. A snapshot search of the same lookup table has no constraint on roughness, is subject to the ambiguity, and the rmse is very large. Monte-Carlo simulation Time-series search of lookup table Snap-shot search of lookup table

8. Truck-mounted radar measurements in Ypsilanti, Michigan were obtained over a two-month campaign (Oh et al. 2002). The LUT time-series and snap-shot performs comparably most likely because the radar measurement error (~0.4dB) is smaller than 0.7dB. Dubois method has outliers. Validation with in situ data

9. The retrievals from all 4 sites are combined in one scatter plot. Many retrievals of the Dubois retrieval become outliers. Overall the LUT time-series retrieval shows the best correlation with the in situ mv and the best rmse in soil moisture estimation. Validation with in situ data

10. The validation with the Michigan data is performed: (a) uses l (correlation length) / s (rms height) of 10 (= the same as forward model). (b) uses the truth. The choice of l/s does not affect the Mv retrieval. Even if the value of the cost function changes, the location of the minimum (=Mv retrieval) does not. Change of l/s merely adds a bias to the cost function Effects of correlation length +.. Retrieved Mv

11. If a single scattering process is dominant, the IEM shows that the roughness effect and the dielectric effect can be decoupled (Fung et al. 1992; Shi et al. 1997): σ0 (Mv, s, l/s) = f(Mv) + g(s, l/s) Then σ0 (Mv, s, l/s=a) - σ0 (Mv, s, l/s= b ) = g(s, l/s=a) - g(s, l/s=b) = independent of Mv. The bias offset is uniform wrt dielectric constant and polarization. Effects of correlation length l/s=4 l/s=15

12. Findings Time-invariant roughness and soil moisture are estimated using the time-series method (2N inputs solve N+1 unknowns). Does not require ancillary information. Simple search of a lookup table. Tested with in situ data: the error is 0.044 cm3/cm3 using 6-11 time-series inputs. The time-series method performs better than the other methods. Retrieval is mostly independent of the knowledge of correlation length  one degree of freedom is reduced. Discussion Increase in time-series window will further reduce the radar measurement noise  improve roughness estimate  improve Mv retrieval. Summary

13. backup

14. Soil moisture: retrieval performance