Subodh Dhakal
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Subodh Dhakal has a total of 8 co-authors. Here are the top co-authors by number of shared publications:
Bharat Prasad Bhandari
Central Department of Environmental Science, Tribhuvan University, Kirtipur, Nepal
1 shared publication- Bharat Prasad Bhandari
Central Department of Environmental Science, Tribhuvan University, Kirtipur, Nepal
1 shared publication - Ching-Ying Tsou
Faculty of Agriculture and Life Science, Hirosaki University, Hirosaki 036-8561, Japan
1 shared publication - Gopi Krishna Sedhain
Nepal Environment Protection Center, Kathmandu 44600, Nepal
1 shared publication - Muna Neupane
Central Department of Environmental Science, Tribhuvan University, Kathmandu 44600, Nepal
1 shared publication - Niraj Bal Tamang
River and Catchment Dynamics, WSP NZ, Christchurch, New Zealand
1 shared publication
Latest Publications
- Open Access
- 24 March 2026
- Published by Nepal Journals Online (JOL) in Nepal Journal of Environmental Science
- Vol. 14 (1) , 51-61
Geological contribution on landslide distribution in Siwaliks of parts of central and eastern Nepal
The Siwalik is the youngest among all the geological units within the Nepal Himalaya, and it is inherently weak and fragile. Conservation of the Siwalik region in Nepal has become a national challenge, and the government is currently committed to bring breakthrough in controlling landslides. Therefore, a study was conducted in the Siwaliks of central and eastern Nepal to understand the spatial distribution of landslides in different geological formations. A geological formation distribution map is prepared by collecting lithological and structural geological data from the field, and the spatial distribution of the landslide is delineated from satellite images as well as field visits. It is found that the Siwaliks in the study area are divisible into Lower, Middle, and Upper Siwaliks. In terms of landslide distribution, the Middle Siwaliks is more hazardous in terms of landslide occurrence than the other two. In an average, 60% of the total landslides in the study area fall within the Middle Siwaliks, 22% fall within Upper Siwaliks, and 13% landslides have occurred in Lower Siwaliks. It is also found that highly weathered rock beds in Lower and Upper Siwaliks favor shallow failures, whereas the easily erodible mudstone layers interbedded with comparatively less erodible sandstone layers in the Middle Siwaliks make it unstable and more prone to rock slides and rock falls. It indicates the control of geology in landslide occurrences and their types. This study will provide insights into the possible geological controls for landslides throughout the Siwaliks and support for preparing a geohazard management strategy in the study area.
- 24 March 2026
- Published by Nepal Journals Online (JOL) in Nepal Journal of Environmental Science
- Vol. 14 (1) , 51-61
Open AccessGeological contribution on landslide distribution in Siwaliks of parts of central and eastern Nepal
The Siwalik is the youngest among all the geological units within the Nepal Himalaya, and it is inherently weak and fragile. Conservation of the Siwalik region in Nepal has become a national challenge, and the government is currently committed to bring breakthrough in controlling landslides. Therefore, a study was conducted in the Siwaliks of central and eastern Nepal to understand the spatial distribution of landslides in different geological formations. A geological formation distribution map is prepared by collecting lithological and structural geological data from the field, and the spatial distribution of the landslide is delineated from satellite images as well as field visits. It is found that the Siwaliks in the study area are divisible into Lower, Middle, and Upper Siwaliks. In terms of landslide distribution, the Middle Siwaliks is more hazardous in terms of landslide occurrence than the other two. In an average, 60% of the total landslides in the study area fall within the Middle Siwaliks, 22% fall within Upper Siwaliks, and 13% landslides have occurred in Lower Siwaliks. It is also found that highly weathered rock beds in Lower and Upper Siwaliks favor shallow failures, whereas the easily erodible mudstone layers interbedded with comparatively less erodible sandstone layers in the Middle Siwaliks make it unstable and more prone to rock slides and rock falls. It indicates the control of geology in landslide occurrences and their types. This study will provide insights into the possible geological controls for landslides throughout the Siwaliks and support for preparing a geohazard management strategy in the study area.
- Top Cited PapersOpen Access
- 31 January 2026
- Published by Elsevier in Geosystems and Geoenvironment
- Vol. 5 (2)
Landslide susceptibility assessment and slope stability analysis of cut slopes along Dumre-Paudhi Bazar road section, western Nepal
- 31 January 2026
- Published by Elsevier in Geosystems and Geoenvironment
- Vol. 5 (2)
Top Cited PapersOpen AccessLandslide susceptibility assessment and slope stability analysis of cut slopes along Dumre-Paudhi Bazar road section, western Nepal
- Open Access
- 25 March 2025
- Published by Nepal Journals Online (JOL) in Journal of Institute of Science and Technology
- Vol. 30 (1) , 89-99
Landslide Susceptibility Assessment in the Marin Khola Watershed of the Sub Himalaya, Central Nepal
Nepal is facing the threat of landslides each year causing huge loss of lives and properties. Landslide prediction and susceptibility assessments help in identifying the potential zones of landslide occurrences and provide opportunities to treat them prior to their occurrence. Among different methods of landslide susceptibility mapping, the InfoVal method is one of the simple and useful methods. In this study, this method is used to study the landslide susceptibility in the Marin Khola Watershed within the Siwaliks of central Nepal as this area comprises of the weak geological formations that contribute to high potentialities of landslides, yet there are no studies for predicting landslides. A total of 217 landslides were taken for the study and they were divided into two groups: working landslides and validating landslides. 75% of these total landslides were selected as working landslides and the remaining 25% were selected for validating landslides. Spatial relationships of the landslide distribution with different causative factors including topographic factors, hydrologic factors, geological factors and landuse factors were employed and analyzed. The results depict that very high, high, moderate, low and very low susceptibility classes cover 1.15%, 49.93%, 30.17%, 11.48%, and 11.28% area, respectively. The Middle Siwaliks are most susceptible to landslides compared to the Upper Siwaliks and Lower Siwaliks. The accuracy values are found to be affected by the difference in the landslide characteristics and types occurring in the study area. The model accuracy remains at 66% and predictive accuracy at 75%.
- 25 March 2025
- Published by Nepal Journals Online (JOL) in Journal of Institute of Science and Technology
- Vol. 30 (1) , 89-99
Open AccessLandslide Susceptibility Assessment in the Marin Khola Watershed of the Sub Himalaya, Central Nepal
Nepal is facing the threat of landslides each year causing huge loss of lives and properties. Landslide prediction and susceptibility assessments help in identifying the potential zones of landslide occurrences and provide opportunities to treat them prior to their occurrence. Among different methods of landslide susceptibility mapping, the InfoVal method is one of the simple and useful methods. In this study, this method is used to study the landslide susceptibility in the Marin Khola Watershed within the Siwaliks of central Nepal as this area comprises of the weak geological formations that contribute to high potentialities of landslides, yet there are no studies for predicting landslides. A total of 217 landslides were taken for the study and they were divided into two groups: working landslides and validating landslides. 75% of these total landslides were selected as working landslides and the remaining 25% were selected for validating landslides. Spatial relationships of the landslide distribution with different causative factors including topographic factors, hydrologic factors, geological factors and landuse factors were employed and analyzed. The results depict that very high, high, moderate, low and very low susceptibility classes cover 1.15%, 49.93%, 30.17%, 11.48%, and 11.28% area, respectively. The Middle Siwaliks are most susceptible to landslides compared to the Upper Siwaliks and Lower Siwaliks. The accuracy values are found to be affected by the difference in the landslide characteristics and types occurring in the study area. The model accuracy remains at 66% and predictive accuracy at 75%.
- Top Cited PapersOpen Access
- 2 March 2024
- Published by MDPI AG in Sustainability
- Vol. 16 (5) , 2092
Assessing the Prediction Accuracy of Frequency Ratio, Weight of Evidence, Shannon Entropy, and Information Value Methods for Landslide Susceptibility in the Siwalik Hills of Nepal
The main objective of this study is to assess the prediction and success rate based on bivariate frequency ratio (FR), weight of evidence (WoE), Shannon entropy (SE), and information value (IV) models for landslide susceptibility in the sedimentary terrain of Nepal Himalaya, as the area is facing threat for sustainable development as well as sustainable resource management. This study also seeks to evaluate the causative factors for landslide susceptibility. Initially, a landslide inventory map was created, consisting of 1158 polygons. These polygons were randomly divided into two sets, with a ratio of 70% for training and 30% for testing data. The multicollinearity approach was evaluated to assess the relevance of selected conditioning variables and their inclusion in the model construction process. The area under the curve (AUC) and other arithmetic evaluation methods were employed to validate and compare the outcomes of the models. In comparison, the predictive accuracy of the FR model surpasses that of the IV and SE models. The success rates, ranked in descending order, are as follows: WoE (79.9%), FR (75.3%), IV (74.4%), and SE (73.2%). Similarly, the success rates of four distinct models, namely WoE, FR, IV, and SE, are 85%, 78.75%, 78.57%, and 77.2%, correspondingly. All models have an accuracy and prediction rate exceeding 70%, making them suitable for assessing landslide susceptibility in the Siwalik Hills of Nepal. Nevertheless, the weight of evidence model provides more precise outcomes than other models. This study is expected to provide important information for road and settlement sustainability in the study area.
- 2 March 2024
- Published by MDPI AG in Sustainability
- Vol. 16 (5) , 2092
Top Cited PapersOpen AccessAssessing the Prediction Accuracy of Frequency Ratio, Weight of Evidence, Shannon Entropy, and Information Value Methods for Landslide Susceptibility in the Siwalik Hills of Nepal
The main objective of this study is to assess the prediction and success rate based on bivariate frequency ratio (FR), weight of evidence (WoE), Shannon entropy (SE), and information value (IV) models for landslide susceptibility in the sedimentary terrain of Nepal Himalaya, as the area is facing threat for sustainable development as well as sustainable resource management. This study also seeks to evaluate the causative factors for landslide susceptibility. Initially, a landslide inventory map was created, consisting of 1158 polygons. These polygons were randomly divided into two sets, with a ratio of 70% for training and 30% for testing data. The multicollinearity approach was evaluated to assess the relevance of selected conditioning variables and their inclusion in the model construction process. The area under the curve (AUC) and other arithmetic evaluation methods were employed to validate and compare the outcomes of the models. In comparison, the predictive accuracy of the FR model surpasses that of the IV and SE models. The success rates, ranked in descending order, are as follows: WoE (79.9%), FR (75.3%), IV (74.4%), and SE (73.2%). Similarly, the success rates of four distinct models, namely WoE, FR, IV, and SE, are 85%, 78.75%, 78.57%, and 77.2%, correspondingly. All models have an accuracy and prediction rate exceeding 70%, making them suitable for assessing landslide susceptibility in the Siwalik Hills of Nepal. Nevertheless, the weight of evidence model provides more precise outcomes than other models. This study is expected to provide important information for road and settlement sustainability in the study area.
- 8 April 2020
- Published by Springer Nature in SN Applied Sciences
- Vol. 2 (5) , 1-17
Spatio-temporal dynamics of landslides in the sedimentary terrain: a case of Siwalik zone of Babai watershed, Nepal
- 8 April 2020
- Published by Springer Nature in SN Applied Sciences
- Vol. 2 (5) , 1-17
Spatio-temporal dynamics of landslides in the sedimentary terrain: a case of Siwalik zone of Babai watershed, Nepal
- Open Access
Climatic Variability and Land Use Change in Kamala Watershed, Sindhuli District, Nepal
This study focuses on the land use change and climatic variability assessment around Kamala watershed, Sindhuli district, Nepal. The study area covers two municipalities and eight Village Development Committees (VDCs). In this paper, land use change and the climatic variability are examined. The study was focused on analyzing the changes in land use area within the period of 1995 to 2014 and how the climatic data have evolved in different meteorological stations around the watershed. The topographic maps, Google Earth images and ArcGIS 10.1 for four successive years, 1995, 2005, 2010, and 2014 were used to prepare the land use map. The trend analysis of temperature and precipitation data was conducted using Mann Kendall trend analysis and Sen’s slope method using R (3.1.2 version) software. It was found that from 1995 to 2014, the forest area, river terrace, pond, and landslide area decreased while the cropland, settlement, and orchard area increased. The temperature and precipitation trend analysis shows variability in annual, maximum, and seasonal rainfall at different stations. The maximum and minimum temperature increased in all the respective stations, but the changes are statistically insignificant. The Sen’s slope for annual rainfall at ten different stations varied between −38.9 to 4.8 mm per year. Land use change and climatic variability have been analyzed; however, further study is required to establish any relation between climatic variability and land use change.
Open AccessClimatic Variability and Land Use Change in Kamala Watershed, Sindhuli District, Nepal
This study focuses on the land use change and climatic variability assessment around Kamala watershed, Sindhuli district, Nepal. The study area covers two municipalities and eight Village Development Committees (VDCs). In this paper, land use change and the climatic variability are examined. The study was focused on analyzing the changes in land use area within the period of 1995 to 2014 and how the climatic data have evolved in different meteorological stations around the watershed. The topographic maps, Google Earth images and ArcGIS 10.1 for four successive years, 1995, 2005, 2010, and 2014 were used to prepare the land use map. The trend analysis of temperature and precipitation data was conducted using Mann Kendall trend analysis and Sen’s slope method using R (3.1.2 version) software. It was found that from 1995 to 2014, the forest area, river terrace, pond, and landslide area decreased while the cropland, settlement, and orchard area increased. The temperature and precipitation trend analysis shows variability in annual, maximum, and seasonal rainfall at different stations. The maximum and minimum temperature increased in all the respective stations, but the changes are statistically insignificant. The Sen’s slope for annual rainfall at ten different stations varied between −38.9 to 4.8 mm per year. Land use change and climatic variability have been analyzed; however, further study is required to establish any relation between climatic variability and land use change.
- Open Access
- 28 December 2016
- Published by MDPI AG in Condensed Matter
- Vol. 2 (1) , 2
Experimental Study of Particle Interactions in Moderate to Dense Granular Shear Flows of Disks
To find the status of multiple collisions in transitional granular flow, moderate to densely packed monosized plastic disks were sheared in experimental 2D shear flow apparatus by applying shear strain rate in the range of 14.8 s−1 to 34.6 s−1. Application of high speed video camera and subsequent image processing techniques precisely measure the spatial positions of the particles involved in making the flow. Collision detection and contact duration finding algorithms were formulated to detect inter particle collision in each time step. The proportion of binary and multiple collisions was quantified in each time step. It is found that the contribution of multi particle collision is 13% for the lowest normalized solid fraction ( = 0.50) and 68% for the highest normalized solid fraction ( = 0.85) under consideration. The inter particle collision time is also found greater than binary collision time in all the flows under consideration. The group size is determined for each time step as an additional length scale associated with multi-particle interactions. Number of groups was made of more than two particles in all the flows under consideration, and the amount of such groups increased on increasing the solid fraction irrespective of the shear rate. In the dense cases (solid fraction > 0.60), rapid increase in the group size and the formation of occasional force chains should be attributed for the change in the stress generation pattern.- 28 December 2016
- Published by MDPI AG in Condensed Matter
- Vol. 2 (1) , 2
Open AccessExperimental Study of Particle Interactions in Moderate to Dense Granular Shear Flows of Disks
To find the status of multiple collisions in transitional granular flow, moderate to densely packed monosized plastic disks were sheared in experimental 2D shear flow apparatus by applying shear strain rate in the range of 14.8 s−1 to 34.6 s−1. Application of high speed video camera and subsequent image processing techniques precisely measure the spatial positions of the particles involved in making the flow. Collision detection and contact duration finding algorithms were formulated to detect inter particle collision in each time step. The proportion of binary and multiple collisions was quantified in each time step. It is found that the contribution of multi particle collision is 13% for the lowest normalized solid fraction ( = 0.50) and 68% for the highest normalized solid fraction ( = 0.85) under consideration. The inter particle collision time is also found greater than binary collision time in all the flows under consideration. The group size is determined for each time step as an additional length scale associated with multi-particle interactions. Number of groups was made of more than two particles in all the flows under consideration, and the amount of such groups increased on increasing the solid fraction irrespective of the shear rate. In the dense cases (solid fraction > 0.60), rapid increase in the group size and the formation of occasional force chains should be attributed for the change in the stress generation pattern. - Open Access
Climate Change Impact and Adaptation Practices in Agriculture: A Case Study of Rautahat District, Nepal
This study identifies the impact of climate change and adaptation practices on agriculture in the Rautahat district of central Nepal by analyzing the atmospheric temperature, rainfall pattern, soil moisture, and direct field survey. The impact and status of crop production systems are emphasized. Primary data on crop production system were collected through household surveys and adaptation practices in crop production were collected through focus group discussions, key informants’ interviews, and direct observations. Time series data on key climatic variables and productivity of major crops were collected from the government sources. Mann–Kendall trend analysis and Sen’s Slope methods were used for the analysis and quantification of temperature and rainfall trends. Spearman’s rank correlation analysis was performed to find the relation of seasonal rainfall with the crop yields. The study showed that the annual average rainfall was decreasing at the rate of 10.21 mm/year and the annual mean temperature had increased at a rate of 0.02 °C/year over the last 30 years, but their variations were found to be statistically insignificant. Seasonal rainfall also increased, except for the post-monsoon rain. Estimation of Soil Moisture Index through remote sensing technique indicates that it has been reduced considerably over the past 15 years at the beginning of the monsoon. It was observed that farmers have been using different adaptation measures like the use of high-yielding varieties of crops, enhanced irrigation system, switching to hybrid seed, and increased access to pesticides. As a result, the yields of major crops including rice, maize, wheat, sugarcane, potatoes, and pulses all showed increasing trends during 1999 to 2014. However, the total costs of production of all crops have increased many fold as a consequence of the cost associated with the arrangements for such adaptation measures and shifts towards market-based commodities.
Open AccessClimate Change Impact and Adaptation Practices in Agriculture: A Case Study of Rautahat District, Nepal
This study identifies the impact of climate change and adaptation practices on agriculture in the Rautahat district of central Nepal by analyzing the atmospheric temperature, rainfall pattern, soil moisture, and direct field survey. The impact and status of crop production systems are emphasized. Primary data on crop production system were collected through household surveys and adaptation practices in crop production were collected through focus group discussions, key informants’ interviews, and direct observations. Time series data on key climatic variables and productivity of major crops were collected from the government sources. Mann–Kendall trend analysis and Sen’s Slope methods were used for the analysis and quantification of temperature and rainfall trends. Spearman’s rank correlation analysis was performed to find the relation of seasonal rainfall with the crop yields. The study showed that the annual average rainfall was decreasing at the rate of 10.21 mm/year and the annual mean temperature had increased at a rate of 0.02 °C/year over the last 30 years, but their variations were found to be statistically insignificant. Seasonal rainfall also increased, except for the post-monsoon rain. Estimation of Soil Moisture Index through remote sensing technique indicates that it has been reduced considerably over the past 15 years at the beginning of the monsoon. It was observed that farmers have been using different adaptation measures like the use of high-yielding varieties of crops, enhanced irrigation system, switching to hybrid seed, and increased access to pesticides. As a result, the yields of major crops including rice, maize, wheat, sugarcane, potatoes, and pulses all showed increasing trends during 1999 to 2014. However, the total costs of production of all crops have increased many fold as a consequence of the cost associated with the arrangements for such adaptation measures and shifts towards market-based commodities.
- 10 November 2014
- Published by Springer Nature
- p. 53-72
Evolution of Geomorphologic Hazards in Hindu Kush Himalaya
- 10 November 2014
- Published by Springer Nature
- p. 53-72
Evolution of Geomorphologic Hazards in Hindu Kush Himalaya