Tannery affected surface soils from 72 sampling sites from industrial area of Sialkot district, Pakistan were collected and analyzed for nine physico-chemical parameters, nine heavy metals and four macro-nutrients. Most of the soils were poor in organic matter (0.11-2.98%), basic in nature with pH (7.1-10.6) and electrical conductivity (1.2-17.9 mS/cm). Mean concentration of total dissolved solids (TDS), Cl1-, alkalinity, NO31-N, salinity and PO43- were 3093 mg/L, 6587, 3929, 301.3, 10.3, and 1.7 mg/kg. The results showed that concentration of macro-nutrients were in the order: Na>Mg>K>Ca whereas heavy metals followed the order: Cr>Fe>Ni>Mn>Cu>Zn>Co>Pb>Cd. Factor Analysis based on Principal Component Analysis (FA/PCA), Cluster Analysis (CA) and correlation analysis identified contribution of metals from tannery effluents, agrochemicals, automobiles exhaust and natural weathering processes. Tannery affected soils were enriched with Cd followed by Cr, Pb, Ni, Cu, Co, Zn and Mn. Geo-accumulation index (Igeo) classified the soil samples in unpolluted to moderately polluted categories. Metal Pollution Index (MPI) provided better estimation of heavy metal pollution as compared to Pollution Load Index (PLI). Ecological risk index (RI) showed high potential ecological risk associated with Cd and Cr with mean concentrations above respective average shale/background values. The results are useful for heavy metals source identification, enrichment, risk assessment, management of tannery affected soils and can contribute to monitoring programs at regional levels.