Identifying the primary source of heavy metals (HMs) pollution and the key pollutants is crucial for safeguarding eco-health and managing risks in industrial vicinity. For this purpose, this investigation was carried out to investigate the pollution area identification with soil static environmental capacity (QI), receptor model-oriented critical sources, and Monte Carlo simulation (MCS) based probabilistic environmental and human health hazards associated with HMs in agricultural soils of Narayanganj, Bangladesh. The average concentration of Cr, Ni, Cu, Cd, Pb, Co, Zn, and Mn were 98.67, 63.41, 37.39, 1.28, 23.93, 14.48, 125.08, and 467.45 mg/kg, respectively. The geoaccumulation index identified Cd as the dominant metal, indicating heavy to extreme contamination in soils. The QI revealed that over 99% of the areas were polluted for Ni and Cd with less uncertain regions whereas Cr showed a significant portion of areas with uncertain pollution status. The positive matrix factorization (PMF) model identified three major sources: agricultural (29%), vehicular emissions (25%), and industrial (46%). The probabilistic assessment of health hazards indicated that both carcinogenic and non-carcinogenic risks for adult male, adult female, and children were deemed unacceptable. Moreover, children faced a higher health hazard compared to adults. For adult male, adult female, and children, industrial operations contributed 48.4%, 42.7%, and 71.2% of the carcinogenic risks, respectively and these risks were associated with Ni and Cr as the main pollutants of concern. The study emphasizes valuable scientific insights for environmental managers to tackle soil pollution from HMs by effectively managing anthropogenic sources. It could aid in devising strategies for environmental remediation engineering and refining industry standards.