In order to explore the spatial distribution and variation characteristics of soil water in coal mining subsidence in Mu Us sandy land, and provide theoretical basis for the restoration of the mining area, experiments based on a linear sampling and classic statistical and geostatistical methods were conducted. Spatial distribution characteristics and variation law of soil water in the typical 0 to 100 cm dune area in the subsidence and the non-collapse (control) were studied. The results showed that in the typical sand dune location of non-collapse (control), the probability distribution curves of soil water changes in all layers along vertical and horizontal directions were all normal distribution, and it was consistent with the temporal and spatial variation characteristics of soil water in conventional dunes in Mu Us sandy land. In contrast, two years after the coal mine collapsed, the variations of soil water in different layers along vertical and horizontal directions were different, and soil water loss was more serious than that of control dune by nearly 10% to 30%, and the standard deviation varied from 0.54 to 1.05, increasing by 52.08% compared with the non-subsidence area (control). The probability of positive and negative deviation greater than 1 was over 50%, and the coefficient of variation varied from 0.14 to 0.28, which was 80% higher than that of non-collapse area (control). After collapsing, the average level of soil water, standard deviation, variance and variation coefficient had great changes, and influence of coal mining subsidence on soil water was the most in the middle layer (30—70 cm), and was not obvious in the surface (0—20 cm) and lower (80—100 cm). In coal mining subsidence, the dispersion degree of soil water in different layers along the vertical and horizontal direction was greatly improved, which enhanced spatial variation of soil water.