Abstract: To provide a theoretical basis for simulating soil hydraulic processes, classical statistics and geo-statistics were used to characterize the spatial distribution of saturated soil hydraulic conductivity (K_s) at a slope scale in the northern Loess Plateau. Furthermore, the first order autoregressive state-space model was applied to simulate the distribution of K_s. The magnitude of K_s variability was moderate according to the coefficient of variation. The degree of spatial dependence was moderate with range of 42 m. Bulk density, sand content, silt content, and clay content were auto-correlated. Cross-correlation existed between K_s and other factors at different lag distance. Bulk density and soil particles were the key factors to affect the spatial distribution of K_s on the sloping land. A combination of bulk density and soil particles was the best to explain the spatial variation of K_s(R²>0.9). The first order autoregressive state-space model is suitable for simulating the K_s distribution at a slope scale.