S1–S3 Figs, S1–S7 Tables.

<div> This study assessed the risks of soil pollution by heavy metals in Chilmari Upazila, northern Bangladesh, using the static environmental resilience (Pi) model of soil. Geostatistical modeling and self-organizing maps (SOM) identified pollution areas and spatial patterns, while a positive matrix factorization (PMF) model revealed pollution sources. The results showed that the average concentrations of Cr, Pb and As were well above background levels. Agricultural and industrial soils were mainly contaminated with Cr, Pb and As according to the Nemerow Pollution Index (NPI), Ecological Risk (ER) and Pi Index. Over 70% of the sites were contaminated with Pb and Cr, while co-contamination was particularly high. A one-way ANOVA showed significant correlations between Pb, Cu and Zn levels and human activities. The PMF analysis revealed that industrial effluents, agrochemicals and lithogenic sources were the main contributors to soil contamination with 16%, 41% and 43%, respectively. The SOM analysis revealed three distinct spatial patterns (Pb-Zn, Cr-Cu-Ni and Co-Mn-As), which are consistent with the PMF results. These results emphasize the need for stringent measures to reduce industrial emissions and remediate soil contamination in order to improve soil quality and food security. </div>