Watershed prioritization is considered an important tool for soil and watershed management. This study focuses on the watershed prioritization of the Kailali district in terms of soil erosion, considering morphometric parameters and land use/landcover (LULC) datasets using GIS. ALOS DEM of 30mresolution was used to delineate sub-watersheds and calculate linear, areal, and relief morphometric parameters. Similarly, Esri LULC 2021 (Sentinel-2 imagery at 10mresolution) was used to calculate LULC parameters. An integrated approach of Principal Component Analysis (PCA) and Weighted Sum Analysis (WSA) was used for prioritization. PCA was used to reduce selected parameters, calculate the correlation matrix, and define the significant parameters. WSA was used to define weightage value, and Compound Value (CV) was calculated for the ranking of sub-watersheds. 22 sub-watersheds with at least 3rd order stream and 15 parameters were selected for prioritization. PCA integrated with WSA was found to be effective for prioritization. The findings showed that about 61.58% of the watershed area is in the high-priority category, suggesting those areas are at a higher risk of erosion. Therefore, different land rehabilitation programs and bioengineering techniques should be focused on the sub-watershed of high-priority categories followed by medium and low-priority categories to control further soil erosion. The adopted methodology of prioritization can also be performed for multi-hazard mapping.