Research paperCurve number modifications and parameterization sensitivity analysis for reducing model uncertainty in simulated and projected streamflows in a Himalayan catchment
Introduction
Due to increments in the air temperature, the vapor carrying capacity is significantly affected across the world (Collins et al., 2013). The exponential rate of change of air temperature has been identified as one of the main causes of climate change (Taylor et al., 2012). The variations in the ratio of precipitation to snowfall and rainfall are increasing non-uniformly around the world (Chamaille-Jammes et al., 2007). A non-uniformity in the precipitation events and amount observed (Khadka et al., 2014). Extreme events will be more frequent in the Himalayan region in coming years (Shivam et al., 2016, Singh and Goyal, 2016, Singh et al., 2015). The Himalayan glacier ranges are the largest and most dynamic in the world (Pachauri et al., 2014, Collins et al., 2013). Several studies reported that the glaciers have been retreating since the end of the Little Ice Age (LIA) and they have been accelerated since 1970s (e.g. Collins et al., 2013). The earlier and accelerated melting rate of permanent snow covers and glaciers are contributing more water to the rivers (Collins et al., 2013). The earlier melting of snow-glaciers due to climate change have largely affected the seasonality and quantity of streamflow over Himalayan watersheds (Collins et al., 2013, Kulkarni et al., 2010). Many studies reveal that the Himalayan glaciers have been losing their mass and are shrinking (Pachauri et al., 2014, Khadka et al., 2014, Krishna, 2005) and that the flow of Himalayan rivers is increasing (Singh et al., 2015). Thus the downstream portions of the rivers are flooding (Radic et al., 2014, Kulkarni et al., 2010).
The eastern Himalayas is recognized as one of the most fragile ecosystems in the world. The poor and marginalized people living in this region face a pressing challenge of adapting to changing climate (Chettri et al., 2009). The shifts in rainfall and snowfall pattern over Himalayas due to climate change will affect the habitat of many species (Chakraborty et al., 2016). In particular, Chettri et al. (2009) and Huss (2011) conducted studies on the Himalaya and found that glaciers in this region have shrunk and the permanent snow cover is significantly decreased. The relative contribution of discharge from glaciers has decreased and the yearly peak discharge has shifted towards spring. Collins et al. (2013) worked on the climatic variations in Himalayan region and deduced that discharge of Ganges river, especially in Nepal, significantly fluctuated. However, the level of flow was usually preserved from the 1960s to 2000s. Bhutiyani et al. (2008) found that the discharge of the Satluj river significantly increased from the year 1991–2004 due to enhanced temperature warming. Bhambri et al. (2011) concluded that the snowmelt induced runoff will be higher in the 21st century over the Himalayan region and that could affect the aquatic ecosystems in Himalayas.
Distributed and deterministic eco-hydrological models have been effectively utilized for the management of water resources. The physical and real time meteorological parameters based hydrological models were found more reliable for the computation of water balance than other modeling approaches (Singh et al., 2013a., Singh et al., 2013b, Abbaspour, 2011, Arnold et al., 2012). The hydrological response unit (HRU) scale based semi-distributed hydrological model Soil and Water Assessment Tool (SWAT) and stochastic optimization tool SWAT Calibration and Uncertainty Program (SWATCUP) have already proven their utilities in various water balance studies around the world (Singh et al., 2013b, Abbaspour, 2011, Schuol et al., 2008, Gosain et al., 2006). To incorporate the present climate projections (e.g. Global Circulation Models: GCMs), we emphasized that the present modeling framework should be defined and explained (Li et al., 2015; Vermote, 2015; Slater et al., 2009; Molotch and Margulis, 2008). We highlighted the need to parameter optimization algorithms and uncertainty analysis to generate more accurate projection scenarios of main hydrological components such as streamflow, water yield and precipitation.
The Teesta river in the eastern Himalayas is characterized by extremely high elevation ranges (e.g. 1400 m to 7400 m, approximately) (Singh and Goyal, 2016). The topography is dominated by snow cover, glaciers and steep slopes (Singh and Goyal, 2016). In previous SWAT based studies, several authors reported unusual high peaks in modeled streamflows, especially when modelling hilly terrain. The high peak in the compute discharge may be occurring due to steep slopes. It is found that this problem is very common in the curve number (CN) based runoff computation (Rahman et al., 2013, Schilling et al., 2008). The unreliable high flow peaks may enhance the uncertainty level, which further makes calculating of water balance unreliable (Singh et al., 2015; Mishra et al., 2014). This will tend to increase an uncertainty level in the projection of water balance components such as water yield and snowfall/snowmelt (Schulz et al., 2008). This could affect the downstream movement of streamflow and may affect the aquatic life of the living species.
The SWAT model has proven its suitability for the calculation of in the projection of water balance components in large river basins across the world (Singh et al., 2015, Schuol et al., 2008, Gosain et al., 2006, Arnold et al., 2012). Gosain et al. (2006) used SWAT for simulating twelve large river basins in India to quantify the climate change impact on hydrology and concluded that the future discharge volume has been enhanced for those rivers originating from the Himalayan glaciers due to the higher melting rate of snow and glaciers. Schuol et al. (2008) used SWAT to analyze the hydrology of the entire Africa continent at sub-catchment spatial scales in the monthly time intervals. Schuol et al. (2008) found a significant variability in the discharge amount during historical and future time domains. In this study, the rainfall-runoff relationship based Soil Conservation Services (SCS) curve number (CN) method utilized for the computation of accurate streamflow discharge (Mishra et al., 2014; Neitsch et al., 2011).
As model uncertainties reduce the confidence level of projection scenarios, the prime objective of this study was to implement the modified CNs in SWAT at a HRU scale by the computation of fractional slopes. The fractional slopes are computed at each hydrological response unit (HRU) scale and then the CNs are adjusted for each HRU to avoid the slope steepness, especially over hilly catchment (Jain et al., 2010). The deterministic models do not compute randomness (Singh et al., 2013a., Singh et al., 2013b). This study incorporates stochastic advance optimization algorithms such as Sequential Uncertainty Parameter Fitting Approach (SUFI2) and Parametric Solution (ParaSol) (Abbaspour, 2011, Abbaspour et al., 2007). The results from both methods are compared to find out the suitability of the optimization methods for hilly catchments. The parameter based sensitivity analysis was done to find out the most influencing parameter over snow glacier induced Teesta river catchment (Singh et al., 2013a., Singh et al., 2013b, Abbaspour, 2011). The optimization algorithms are found helpful to reduce various kinds of model uncertainties (e.g. parameter uncertainty, input data related uncertainties and model structure uncertainty) in calibrated and projected modeling outcomes (Singh et al., 2015, Singh et al., 2013a., Singh et al., 2013b).
Another objective of this study was to evaluate the effects of climate change in the projected water yield and streamflows over the snow-glacier induced Teesta river Himalayan catchment. The projection scenarios were generated utilizing the Coupled Model Intercomparison Project Phase Five (CMIP5) Global Circulation Models (GCMs) with their Representative Concentration Pathway (RCP) experiments (Taylor et al., 2012). The projected streamflow and water yield are characterized at each sub-catchment scale to highlight the spatial variations due to climate change. The annual surface runoff and water yield were analyzed at each sub-catchment scale in the future time domain to assess increase and decrease in their amount. The percentage of change analysis was conducted to determine the magnitude of change of water yield at each sub-catchment scale, which may be useful to highlight the stress over watershed ecosystem and their components in a Himalayan region even at smaller scales. The Teesta river catchment has very diverse ecology and thus the any change in streamflow and water yield with respect to glaciered and non-glaciered portions may affect the biodiversity of the region.
Section snippets
Study area description
The Teesta river catchment (up to Chungthang gauge) part of Indus river basin has been selected for the analysis, which is a part of eastern Sikkim Himalaya India. The flow in the Teesta river originates from many small and large tributaries. Lachung river is one of the major tributaries of Teesta river joins the Teesta river at Chungthang gauge and also contributed a significant amount of flow in the downstream of Teesta river. The Teesta river originates from the Chhombo Chhu and Khangchun
Climate data
The daily observed gridded (0.5° × 0.5°) precipitation and temperature datasets for the year 1980–2005 (around 26 years) have been used for the analysis. These datasets were collected from the Indian Meteorological Department and Indian Institute of Tropical Meteorology. These gridded datasets are prepared using more than 1800 observed precipitation and temperature gauges distributed across the India. These gridded precipitation and temperature datasets significantly utilized in various research
Soil and water assessment tool (SWAT)
SWAT is defined as semi-distributed, deterministic, and continuous time scale based hydrological model, which is fully capable of simulating and forecasting various water balance components at Hydrological Response Unit (HRU) and sub-catchment scale. Around 2400 peer reviewed publications are available related to SWAT model applications (Neitsch et al., 2011). The modeling outcomes can be simulated on the daily, monthly and yearly time steps (Neitsch et al., 2011, Arnold et al., 2012). In SWAT
Results
Initially, the streamflow was computed as per the single slope (or constant slope) method and it figured highly uncertain at both catchment outlets shown in regression plots (Fig. 3a). Based on a single slope, the computed R2 is found less corresponding to observed and simulated stream flows and the fitted trend line look highly uncertain (Fig. 3a). Thus, multiple HRU slopes were computed and incorporated into SWAT modeling. The multiple HRU slope method improved the correlation between
Discussion
The CN method extensively used in various rainfall runoff and agricultural watershed management studies but CN method found uncertain in the streamflow simulation for hilly catchments. Many authors reported inconsistencies in the streamflow simulation over hilly catchments using SWAT inbuilt CN method (Singh et al., 2015; Zhang and Pan, 2014; Rahman et al., 2013; Huang et al., 2006). Previous Jha et al. (2014) and Huang et al. (2006) concluded that the CN cannot account the effect of slope. We
Conclusion
A hydrological model SWAT has been applied for the projection and evaluation of streamflows and water yield in a snow/glacier induced Teesta river Himalayan catchment. An advanced stochastic optimization tool SWATCUP inbuilt SUFI2 and ParaSol algorithms have been utilized for the model calibration, validation and parameter uncertainty analysis. The SWAT model performed well in a very complex hilly catchment, as all the objective functions (e.g. p-factor, r-factor, R2 and NSE) reached
Acknowledgments
This present research work has been carried out under DST research project No. SB/DGH-66/2013 and financial support is gratefully acknowledged. We also thank to GFDL, New Jersey USA, Earth Explorer NASA, CWC New Delhi India, IMD Pune India and IITM Pune India for providing the required datasets for the completion of this research work.
References (49)
- et al.
Modelling hydrology and water quality in the pre-alpine/alpine Thur watershed using SWAT
J. Hydrol.
(2007) - et al.
Predicting distribution of major forest tree species to potential impacts of climate change in the central Himalayan region
Ecol. Eng.
(2016) - et al.
Climatic variation and runoff from partially-glacierized Himalayan tributary basins of the Ganges
Sci. Total Environ.
(2013) - et al.
Development of a snowfall-snowmelt routine for mountainous terrain for the soil water assessment tool (SWAT)
J. Hydrol.
(2002) - et al.
Climate change impact on glacier and snow melt and runoff in Tamakoshi basin in the Hindu Kush Himalayan (HKH) region
J. Hydrol.
(2014) - et al.
Individual and combined impacts of future climate and land use changes on the water balance
Ecol. Eng.
(2017) - et al.
An adaptive inverse-distance weighting spatial interpolation technique
Comput. Geosci.
(2008) - et al.
Estimating the distribution of snow water equivalent using remotely sensed snow cover data and a spatially distributed snowmelt model: a multi-resolution, multi-sensor comparison
Adv. Water Resour.
(2008) - et al.
River flow forecasting through conceptual models part I—a discussion of principles
J. Hydrol.
(1970) - et al.
Estimation of freshwater availability in the West African sub-continent using the SWAT hydrologic model
J. Hydrol.
(2008)
Analysis and trends of precipitation lapse rate and extreme indices over north Sikkim eastern Himalayas under CMIP5ESM-2M RCPs experiments
Atmos. Res.
Use of the data depth function to differentiate between case of interpolation and extrapolation in hydrological model prediction
J. Hydrol.
Comment on “Estimating the distribution of snow water equivalent using remotely sensed snow cover data and a spatially distributed snowmelt model: a multi-resolution, multi-sensor comparison” by Noah P. Molotch and Steven A. Margulis (Adv. Water Resour. 31 (2008) 1503–1514)
Adv. Water Resour.
Predicting plausible impacts of sets of climate and land use change scenarios on water resources
Appl. Geogr.
Uncertainty of runoff projections under changing climate in Wami River sub-basin
J. Hydrol.: Reg. Stud.
An urban storm-inundation simulation method based on GIS
J. Hydrol.
SWAT-CUP4: SWAT Calibration and Uncertainty Programs—A User Manual
SWAT
Trans. ASABE
Mapping of debris-covered glaciers in the Garhwal Himalayas using ASTER DEMs and thermal data
Int. J. Remote Sens.
Changing streamflow patterns in the rivers of northwestern Himalaya: implications of global warming in the 20th century
Curr. Sci.
Climate-driven fluctuations in surface-water availability and the buffering role of artificial pumping in an African savanna: potential implication for herbivore dynamics
Austral Ecol.
Climate Change Impacts and Vulnerability in the Eastern Himalayas
Effective and efficient global optimization for conceptual rainfall-runoff models
Water Resour. Res.
Climate change impact assessment on hydrology of Indian river basins
Curr. Sci.
Cited by (13)
Glacier change and glacier runoff variation in the Himalayan Baspa river basin
2021, Journal of HydrologyCitation Excerpt :The R2 values are computed as 0.62 and 0.86 for TS1 (i.e. 2003–2010) and TS2 (i.e. 2011–2018), respectively. The computed R2 values on a daily scale are found satisfactory as compared to the previous studies (Singh and Goyal, 2017; Jain et al., 2010) The modeling accuracy of snow-glacier areas will majorly depend on how snow and glacier parameters are accurately input to the model (Abbaspour et al., 2015). The water balance components such as simulated flow (Sim flow), Rain flow (rainfall induced flow or runoff), Snow melt, Glacier melt and Base flow have been computed at the outlet i.e. Sangla gauge (Fig. 7), which shows an average annual distribution of each component during the years 2003 to 2018 in two-time series sets (i.e. TS1: 2003–2010 and TS2: 2011–2018 (total 16 years)).
Climate change impact on water balance and hydrological extremes in different physiographic regions of the West Seti River Basin, Nepal
2021, Ecohydrology and HydrobiologyCitation Excerpt :The low flows also indicate less water in the basin during the dry season, which may result in drought and that would affect meeting the water demand as well. These variations in the frequency of high flow and scarcity of water will hinder the volume and quality of proper drinking water and agricultural activity as well as ecosystem planning (Singh and Goyal, 2017). Therefore, proper watershed and water management interventions like afforestation, on-farm conservation, infiltration ponds, reservoirs operation policies, water harvesting measures, and water demand management should be promoted.
Activation soil moisture accounting (ASMA) for runoff estimation using soil conservation service curve number (SCS-CN) method
2020, Journal of HydrologyCitation Excerpt :As a result, the method is widely preferred by the hydrologists, engineers and watershed managers as an independent simple watershed model (Ponce and Hawkins, 1996; Garen and Moore, 2005; Singh et al., 2010; Hawkins, 2014; Bartlett et al., 2016), as well as the runoff estimating component in many complex process based model for water availability, soil erosion, flood control, and water quality such as Storm Water Management Model (SWMM) (Metcalf and Eddy, 1971); Hydrologic Engineering Center-1 (HEC-1) (HEC, 1981); Chemicals, Runoff, and Erosion from Agricultural Management Systems (CREAMS) (Smith and Williams, 1980), Agricultural Non-Point Source (AGNPS) (Young et al., 1989a,b) and Soil and Water Assessment Tool (SWAT) (Neitsch et al., 2002) to name a few of them. Singh and Goyal (2017) incorporated the modified CN relationships (Mishra et al., 2008) in to SWAT model to reduce the model uncertainty in simulated and projected streamflows in a Himalayan catchment. More recently, Zhang et al. (2019) incorporated rainfall intensity depended CN (ICN) relationships in the SWAT model (SWAT-ICN) to improve the rainfall-runoff simulation process.
Assessment of impact of climate change on streamflow and soil moisture in Pare watershed of Arunachal Pradesh, India
2024, Sustainable Water Resources ManagementProjected Discharge of Dudhnai River: A Tributary of the Brahmaputra River
2023, Journal of the Indian Society of Remote Sensing