Abstract: Based on grid sampling and laboratory analysis, spatial variability of surface soil nutrients was analyzed with GS⁺ and other statistics methods on the landslide area of Fenghuang Mountain, Leigu Town, Beichuan County. The results showed that except for high variability of available phosphorus, other soil nutrients exhibited moderate variability. The ratios of nugget to sill of the soil available phosphorus and soil organic carbon were 27.9% and 28.8%, respectively, showing moderate spatial correlation, while the ratios of nugget to sill of the total nitrogen (20.0%), total phosphorus (24.3%), total potassium (11.1%), available nitrogen (11.2%), and available potassium (22.7%) suggested strong spatial correlation. The total phosphorus had the maximum range (1232.7 m), followed by available nitrogen (541.27 m), total nitrogen (468.35 m), total potassium (136.0 m), available potassium (128.7 m), available phosphorus (116.6 m), and soil orga-nic  carbon (93.5 m). Soil nutrients had no significant variation with the increase of altitude, but gradually increased from the landslide area, the transition area, to the littleimpacted area. The total and available phosphorus contents of the landslide area decreased by 10.3% and 79.7% compared to that of the littleimpacted area, respectively. The soil nutrient contents in the transition area accounted for 31.1%-87.2% of that of the little-impacted area, with the least and largest reductions in total and available phosphorus contents, respectively. The disturbance by earthquake landslide was the dominant reason for the spatial variability of surface soil nutrients.