The evaluation of soil bulk density (BD) is essential for studying soil processes like water retention, air circulation, and carbon sequestration. Conventional direct BD measurement methods are often time-consuming, costly, and laborious. In contrast, the use of pedotransfer functions (PTFs) based on readily available soil properties offers a cost-effective alternative approach that can provide indirect quick estimates of soil BD. Although several region-specific PTFs are available, none is tailored for use in the hilly North-Western Himalayas (NWH) area. Thus, this work developed a novel PTF to estimate the soil BD of 300 soil samples collected from the NWH area and analyzed by stepwise regression (SR), partial least-squares regression (PLSR), multiple linear regression (MLR), linear regression (LR), and advanced machine-learning algorithms, including artificial neural networks (ANN) and decision tree (DT). The site-specific SR-based equation performed better than the PLSR model by explaining up to 72.2% variability of BD and achieving a coefficient of determination (R2) of 0.73, a mean absolute error (MAE) of 0.0379 g cm⁻³, and a root mean square error (RMSE) of 0.0487 g cm⁻³ for the BD prediction. Soil organic matter and clay contents were the most significant contributors, showing the highest importance in the projection. Furthermore, the performance of the SR-based PTF was validated by comparison with 63 PTFs currently in use, which showed to perform poorly and vary widely in BD estimation (R2: 0.006–0.71, RMSE: 0.0649–0.405 gcm− 3, MAE: 0.0512–0.395 g cm− 3). The promising SR equation and PTF developed in this work showed high accuracy and excellent suitability for stakeholders and professionals in the study area and the scientific community as a reliable BD estimator for hilly soils. This study introduces proposed the first-ever site-specific PTF for BD estimation in hilly soils, contributing to scientific challenge by enhancing soil database and improving BD predictions in complex terrains worldwide. In terms of global significance, the outcomes of this study are essential for soil quality and vulnerability studies, resilience indicators estimates carbon sequestration initiatives, and digital soil mapping, and can also be used to construct a BD database for hilly soils across the globe.