The efficient (site-specific) management of soil nutrients is possible by understanding the spatial variability in distribution of phyto-available nutrients (here after called available nutrients) and identifying the soil management zones (MZs) of agricultural landscapes. There is need for delineating soil MZs of agricultural landscapes of the world for efficient management of soil nutrients in order to obtain sustainability in crop yield. The present study was, therefore, undertaken to understand the spatial distribution pattern of available micronutrients (zinc (Zn), boron (B), iron (Fe), manganese (Mn) and copper (Cu)), available sulphur (S), and soil properties (soil acidity (pH), electrical conductivity (EC) and organic carbon (SOC) content) in soils of intensively cultivated Indo-Gangetic Plain (IGP) of India and to delineate soil MZs for efficient management of soil nutrients. Totally, 55101 soil samples from 0–15 cm depth were obtained from 167 districts of IGP during 2014 to 2017 and were analysed for different soil parameters. Soil pH, EC and SOC content varied from 4.44 to 9.80, 0.02 to 2.13 dS m-1 and 0.10 to 1.99%, respectively. The concentration of available Zn, B, Fe, Mn, Cu and S varied from 0.01 to 3.27, 0.01 to 3.51, 0.19 to 55.7, 0.05 to 49.0, 0.01 to 5.29 and 1.01 to 108 mg kg-1, respectively. Geostatistical analysis resulted in varied distribution pattern of studied soil parameters with moderate to strong spatial dependence. The extent (% area) of nutrient deficiencies in IGP followed the order: S > Zn > B > Mn > Cu > Fe. Principal component analysis and fuzzy c-means clustering produced six distinctly different soil MZs of IGP for implementation of zone-specific soil nutrient management strategies for attaining sustainability in crop yield. The developed MZ maps could also be utilized for prioritization and rationalization of nutrients supply in IGP of India.