This study describes application of chemometric multi-way modeling approach to analyze the dataset pertaining to soils of industrial area with a view to assess the soil/sub-soil contamination, accumulation pathways and mobility of contaminants in the soil profiles. The three-way (sampling depths, chemical variables, sampling sites) dataset on heavy metals in soil samples collected from three different sites in an industrial area, up to a depth of 60 m each was analyzed using three-way Tucker3 model validated for stability and goodness of fit. A two component Tucker3 model, explaining 66.6% of data variance, allowed interpretation of the data information in all the three modes. The interpretation of core elements revealing interactions among the components of different modes (depth, variables, sites) allowed inferring more realistic information about the contamination pattern of soils both along the horizontal and vertical coordinates, contamination pathways, and mobility of contaminants through soil profiles, as compared to the traditional data analysis techniques. It concluded that soils at site-1 and site-2 are relatively more contaminated with heavy metals of both the natural as well as anthropogenic origins, as compared to the soil of site-3. Moreover, the accumulation pathways of metals for upper shallow layers and deeper layers of soils in the area were differentiated. The information generated would be helpful in developing strategies for remediation of the contaminated soils for reducing the subsequent risk of ground-water contamination in the study region.