Spatial variability of van Genuchten model parameters and its affecting factors in 0～20cm and 20～40cm soil layers were studied with multifractal and joint multifractal methods. The results showed that in 0～20cm and 20～40cm soil layers parameter had obvious multifractal characteristics, and its spatial varia b ility was strong caused by the distribution of lower and higher values respectiv ely. The multifractal characteristics of parameter n and θs were not obvious, and the spatial variability was weak. Correlation characteristics betwe en van Genuchten model parameters α and soil physical properties in 0～20 cm and 20～40cm soil layers were analyzed. At the observation scale, the corre lations between parameter n and sand and silt content were significant, the parameter n and sand and organic matter content had remarkable relationship, and θs was significantly related to sand content in 0～20cm soil layer. The relationship between parameters α and organic matter content was remarkabl e, and the parameter n and θs were significantly related to sand, sil t and organic matter content in 20～40cm soil layer. At the multiple scales, c orrelations between parameters α and sand and silt content were most obviou s in 0～20cm soil layer, and parameters α were remarkably related to bulk density and organic matter content in 20～40cm soil layer. 