Abstract: Abstract: The soil As pollution in agricultural land has received increasing attention globally due to its significant potential harm to the environment and human health. Knowing the source as well as the characteristics of pollutants is the premise for their control and treatment. But now, sampling for understanding pollutants in soil is difficult and costly resulting in limited number of soil samples from the As contaminated land. For the problem that it is hard to accurately and efficiently map the spatial distribution of soil As concentration based on limited soil samples. In this study, we collected a total of 104 top soil samples (0~20cm) from agricultural land, in a southeast of Guangdong province. A combined approach of spatial analysis and multiple regression modeling was developed to recognize the statistical characteristic of As concentration in sampling soil, and to reveal their formation causes of spatial variation and associated spatial patterns in vegetable field, paddy and orchard. Results showed that in general, As concentration of soil samples was largely not exceeding the national standard in the agricultural land of the studied area. However, As concentration of few soil samples was still exceeding the National Environment Quality Standard for soil (GB15618-1995), with the maximum value was 137.80 mg/kg. Among three different types of agricultural lands, the soil As concentration in vegetable field was the highest, followed by the paddy and orchard, with mean values of 11.04, 9.89 and 2.54 mg/kg, respectively. The sources of soil As contamination were variable in different agricultural land. In vegetable field, soil As concentration was closely interrelated with shortest distances between samples and chimneys and rivers, and the slope of sample sites. Significance value of the soil As concentration simulation model in vegetable field was 0.010 (P<0.05). In paddy, the dominant factor of soil As concentration was the shortest distances between samples and rivers. Significance value of the soil As concentration simulation model in paddy was 0.044 (P<0.05). However, the multiple linear regression models of soil As concentration simulation in orchard were not significant (P>0.05) in this study. The results suggested that the soil As contamination in vegetable field mainly was affected by the river sewage irrigation and surface runoff of polluted gas from chimneys. The area with high soil As concentration was distributed in the northern part, and lower in other regions. In paddy, the soil As contamination was significantly only affected by the sewage irrigation, with higher As contamination in the north and south, and lower in other regions. But no obvious source of soil As pollution was found in orchard soil. This study would provide a scientific method for accurately mapping the spatial distribution of heavy metal pollution in agricultural land and consequently guide the policy making of precise prevention and control based on sparse sampling data.