A pioneering study employed a holistic geostatistical approach to predict the spatial variability of a non sampled area in the Chenab River, Pakistan, using kriging interpolation for organochlorine pesticide (OCP)-polluted risk zones. The Present research intended to investigate the carcinogenic and non-carcinogenic human health risks, contamination levels, and spatial variation of OCPs in the Chenab River, Pakistan. The residual OCP content in sediment samples (n = 120) ranged from 0.056 to 32.14 ng/g. DDE and α-HCH were prevalent among all the samples analyzed, with mean concentrations of 15.84 ± 8.02 and 12.45 ± 6.72 ng/g, respectively. The order of magnitude of OCPs in sediment samples was DDTs > α-HCH > chlorothalonil > heptachlor > endosulfan > aldrin > dieldrin. The findings of the single (SPI) and Nemerow (Nel) pollution index of α-HCH, heptachlor, and aldrin depicted the Chenab River as a serious pollution risk zone. The outcomes of the Pearson correlation coefficient analysis represent the positive correlation among all OCPs, revealing the common origin. Distribution trends showed substantially higher (p 10-4 illustrated a substantial cancer health risk posed by α-HCH, heptachlor, aldrin, and dieldrin in the downstream zone. We recommend the urgent cessation of the ongoing discharge of OCPs into the Chenab River, which needs to be highlighted owing to the significant cancer risk to public health to ensure the good health and wellbeings.