Increasing cadmium (Cd) accumulation in agricultural soils is undesirable due to its hazardous influences on human health. Thus, having more information on spatial variability of Cd and factors effective to increase its content on the cultivated soils is very important. Phosphate fertilizers are main contamination source of cadmium (Cd) in cultivated soils. Also, crop rotation is a critical management practice which can alter soil Cd content. This study was conducted to evaluate the effects of long-term consumption of the phosphate fertilizers, crop rotations, and soil characteristics on spatial variability of two soil Cd species (i.e., total and diethylene triamine pentaacetic acid (DTPA) extractable) in agricultural soils. The study was conducted in wheat farms of Khuzestan Province, Iran. Long-term (27-year period (1980 to 2006)) data including the rate and the type of phosphate fertilizers application, the respective area, and the rotation type of different regions were used. Afterwards, soil Cd content (total or DTPA extractable) and its spatial variability in study area (400,000 ha) were determined by sampling from soils of 255 fields. The results showed that the consumption rate of di-ammonium phosphate fertilizer have been varied enormously in the period study. The application rate of phosphorus fertilizers was very high in some subregions with have extensive agricultural activities (more than 95 kg/ha). The average and maximum contents of total Cd in the study region were obtained as 1.47 and 2.19 mg/kg and DTPA-extractable Cd as 0.084 and 0.35 mg/kg, respectively. The spatial variability of Cd indicated that total and DTPA-extractable Cd contents were over 0.8 and 0.1 mg/kg in 95 and 25 % of samples, respectively. The spherical model enjoys the best fitting and lowest error rate to appraise the Cd content. Comparing the phosphate fertilizer consumption rate with spatial variability of the soil cadmium (both total and DTPA extractable) revealed the high correlation between the consumption rate of P fertilizers and soil Cd content. Rotation type was likely the main effective factor on variations of the soil DTPA-extractable Cd contents in some parts (eastern part of study region) and could explain some Cd variation. Total Cd concentrations had significant correlation with the total neutralizing value (p < 0.01), available P (p < 0.01), cation exchange capacity (p < 0.05), and organic carbon (p < 0.05) variables. The DTPA-extractable Cd had significant correlation with OC (p < 0.01), pH, and clay content (p < 0.05). Therefore, consumption rate of the phosphate fertilizers and crop rotation are important factors on solubility and hence spatial variability of Cd content in agricultural soils.