As the field observations of K_s and K are time-consuming, and K_s and K are sensitive to soil texture (Williams et al. 1983; Kaczmarek et al. 2016), the pedo-transfer function method (Sun et al. 2016) was applied to get the values of K_s and K. The soil property estimation module integrated in the SPAW hydrological model (Saxton & Rawls 2006) was applied to calculate K_s. K was calculated with the erosion productivity impact calculator (Williams et al. 1983). Since there are two layers for each soil type recorded in HWSD, the K_s and K values for each soil type are depth weighted. For the water bodies, the values of K_s and K were set as 1 to exclude them from soil. The statistics of the spatial distributions of K_s, K, and K_s·K in YLX, WJB, and HQ are shown in Table 1.