Soil Moisture Estimation for the Chinese Loess Plateau Using MODIS-derived ATI and TVDI

Timely and effective estimation and monitoring of soil moisture (SM) provides not only an understanding of regional SM status for agricultural management or potential drought but also a basis for characterizing water and energy exchange. The apparent thermal inertia (ATI) and Temperature Vegetation Dryness Index (TVDI) are two widely used indices to reflect SM from remote sensing data. While the ATI-based model is routinely used to estimate the SM of bare soil and sparsely vegetated areas, the TVDI-based model is more suitable for areas with dense vegetation coverage. In this study, we present an iteration procedure that allows us to identify optimal Normalized Difference Vegetation Index (NDVI) thresholds for subregions and estimate their relative soil moisture (RSM) using three models (the ATI-based model, the TVDI-based model, and the ATI/TVDI joint model) from 1 January to 31 December 2017, in the Chinese Loess Plateau. The initial NDVI (NDVI0) was first introduced to obtain TVDI value and two other thresholds of NDVIATI and NDVITVDI were designed for dividing the whole area into three subregions (the ATI subregion, the TVDI subregion, and the ATI/TVDI subregion). The NDVI values corresponding to maximum R-values (correlation coefficient) between estimated RSM and in situ RSM measurements were chosen as optimal NDVI thresholds after performing as high as 48,620 iterations with 10 rounds of 10-fold cross-calibration and validation for each period. An RSM map of the whole study area was produced by merging the RSM of each of the three subregions. The spatiotemporal and comparative analysis further indicated that the ATI/TVDI joint model has higher applicability (accounting for 36/38 periods) and accuracy than the ATI-based and TVDI-based models. The highest average R-value between the estimated RSM and in situ RSM measurements was 0.73 +/- 0.011 (RMSE-root mean square error, 3.43 +/- 0.071% and MAE-mean absolute error, 0.05 +/- 0.025) on the 137th day of 2017 (DOY-day of the year, 137). Although there is potential for improved mapping of RSM for the entire Chinese Loess Plateau, the iteration procedure of identifying optimal thresholds determination offers a promising method for achieving finer-resolution and robust RSM estimation in large heterogeneous areas.