Runoff Simulation Research Articles

Proposal for Low-Cost Optical Sensor for Measuring Flow Velocities in Aquatic Environments

The ocean, with its intricate processes, plays a pivotal role in shaping marine life, habitats, and the Earth’s climate. This study addresses issues such as beach erosion, the survival of propagules from species like Posidonia oceanica, and nutrient distribution. To tackle these challenges, we propose an innovative sensor that quantifies hydrodynamic velocity by measuring the output voltage derived from detecting changes in light absorption and scattering using LEDs and LDRs. Our results not only demonstrate the effectiveness of the sensor but also the accuracy of the processing algorithm. Notably, the blue LED exhibited the lowest mean relative error of 7.59% in freshwater, while the yellow LED was most precise in chlorophyll-containing water, with a mean relative error of 6.80%. In a runoff simulation, we observed similar velocities with the blue, green, and white LEDs, 6.89 cm/s, 6.99 cm/s, and 7.05 cm/s, respectively, for nearly identical time intervals. It is important to highlight that our proposed sensor is not only effective but also highly cost-efficient, representing less than 0.43% of the cost of a Nortek Vector 6 MHz and 0.18% of the Teledyne Workhorse II 300 kHz Marine. This makes it a key tool for managing marine ecosystems sustainably.

Runoff simulation of the Kaidu River Basin based on the GR4J-6 and GR4J-6-LSTM models

Study regionThe Kaidu River Basin originates from the southern slope of the Tienshan Mountains in the Xinjiang Uygur Autonomous Region, China. Study focusAccurate runoff simulation and prediction significantly affect flood control, drought resilience, and water resource allocation decisions. This study establishes the GR4J-6 model (modèle du Génie Rural à 4 paramètres Journalier-6, including a snowmelt module) and integrates it with the LSTM (Long Short-Term Memory) model to construct the hybrid GR4J-6-LSTM model and enhance the simulation accuracy of snowmelt runoff. A case study is conducted in the Kaidu River Basin to demonstrate the applicability of these models in cold and arid regions. The accuracy of the GR4J-6, LSTM, and GR4J-6-LSTM models is evaluated using Nash-Sutcliffe efficiency (NSE), Kling-Gupta Efficiency (KGE), and Root Mean Squared Error (RMSE) metrics. In addition, the contributions of each feature variable in the models are analyzed using the SHapley Additive exPlanations (SHAP) method to enhance the reliability of the results. New hydrological insights for the regionThe GR4J-6 model demonstrated good applicability in the Kaidu River Basin, with NSE, KGE, and RMSE values of 0.69, 0.79, and 39.39 m3/s during the validation period, respectively. The hybrid model GR4J-6-LSTM exhibited the highest comprehensive accuracy among all the models, with NSE, KGE, and RMSE values of 0.84, 0.87, and 28.79 m3/s, respectively. In the LSTM model, temperature and precipitation were found to significantly influence the simulated runoff, indicating that higher temperature and precipitation lead to increased runoff. In the GR4J-6-LSTM model, Tmin (minimum temperature) and the hydrological feature variable Qsim exhibited a strong positive correlation with simulated runoff, as Tmin and Qsim increased, they promoted stronger flow production. This study provides a framework for runoff simulation in snowmelt river basins, offering a reference for projecting extreme hydrological events under climate change.

IRainSnowHydro v1.0: A distributed integrated rainfall-runoff and snowmelt-runoff simulation model for alpine watersheds

Snowmelt runoff is an essential runoff component in alpine watersheds. On the Tibetan Plateau, the complex hydrometeorological and underlying surface conditions make a single runoff generation mode (either snowmelt-runoff or rainfall-runoff) cannot accurately simulate the runoff process. In this study, we developed a new method that combines the curve number, topographic index, and fractional snow cover to identify the sub-basin seasonal dominant runoff generation mode within the Jinsha River Basin. By constructing a surface ‘snow reservoir’ to depict snow melting impact on runoff generation, and quantitatively classifying the precipitation composition, an innovative integrated hydrological model named the distributed integrated Rainfall-runoff and Snowmelt-runoff simulation Hydrological model (iRainSnowHydro) is developed. With model, a method for identifying the seasonal varying dominant runoff generation mode is proposed. The results show that most sub-basins experience both snowmelt and rainfall driven runoff generation in spring, with snowmelt occurring earlier in regions of lower latitude and elevation. Besides, iRainSnowHydro performs well in daily runoff simulations at Zhimenda and Shigu stations with Nash coefficients of 0.81 and 0.85 in the calibration period, and 0.72 and 0.81 in the validation period. The correlation coefficient ranges from 0.92 to 0.96. Additionally, calculation through iRainSnowHydro indicates a noteworthy percentage of spring snowmelt water. Notably, the Zhimenda watershed, characterized by higher latitudes and elevations, displays an escalating trend from 56.6 % to 78.9 % of total precipitation for spring snowmelt water between 2014 and 2020, while the Shigu watershed maintains stable within 27 % ± 6 %. The methodologies outlined bear significance for simulating and predicting runoff in alpine watersheds and offers valuable insights into how snow cover responds to climate change on the Tibetan Plateau.

Comparative Evaluation of Evapotranspiration and Optimization Schemes for Green Roof Runoff Simulations Using HYDRUS-1D

The use of green roofs, a low-impact development practice, can be an effective means of reducing direct runoff in urban centers. Green roof modeling can enable efficient design by preliminarily grasping the behavior of the green roof system according to specific configurations. In this study, we aimed to find appropriate evapotranspiration and parameter optimization schemes for HYDRUS-1D, a commonly used modeling tool for green roofs. Comparative studies of this sort in the context of green roof runoff modeling have not been conducted previously and are important in guiding users to overcome the difficulties of choosing the right numerical schemes for an accurate prediction of runoff from a green roof. As a study site, the Portland Building Ecoroof in Portland, Oregon, USA, was chosen, as green roof configurations and observed data for climate and runoff were available. From the simulation results of the runoff volume, the Blaney–Criddle method, which was considered an alternative, was found to be appropriate for calculating evapotranspiration from a green roof (R2 = 0.82) relative to the Hargreaves method built in HYDRUS-1D (R2 = 0.46). In addition, this study showed that the optimization method using the harmony search algorithm, which was proposed as an alternative optimizer, was better (R2 = 0.95) than that of the HYDRUS-1D’s own optimization module (R2 = 0.82) in calibrating HYDRUS-1D for green roof runoff. The findings are thought to be useful in guiding modelers who are considering using HYDRUS-1D for green roof runoff simulations.

Applying Low-Impact Development Techniques for Improved Water Management in Urban Areas

Worldwide, the increase in impervious surfaces due to urbanization has led to significant water cycle issues such as groundwater depletion, urban heat islands, and flooding. To address these challenges, Low-Impact Development (LID) techniques are increasingly being applied in stormwater management. This study focuses on Ulsan, designated as a water cycle model city in South Korea, with a particular emphasis on the highly urbanized Okgyo drainage watershed. Using the Stormwater Management Model (SWMM) version 5.1, long-term runoff simulations were conducted to evaluate the effects of LID implementation on water cycle change rates and recovery rates. The model incorporates detailed hydrological and hydraulic parameters, including inflow, runoff, infiltration, and evapotranspiration for six subcatchments within the watershed. The SWMM was calibrated and validated using 30 years of historical rainfall data (1987–2016) from the Ulsan weather station. Calibration and validation processes used the NRCS-CN (Curve Number) method to ensure accuracy in simulating runoff patterns and water balance. The study specifically evaluated the effectiveness of two LID techniques: bioretention and permeable pavements. Three scenarios were modeled: bioretention applied to 5% of the area, permeable pavements applied to 5% of the area, and a combined application of both techniques. The results showed that the combined scenario provided the best outcome, with a 7.80% reduction in surface runoff and a 14.56% improvement in water cycle health. The LID application scenario confirmed the potential to achieve the water cycle management target of handling 25.5 mm of rainfall. These findings demonstrate that the introduction of LID techniques in public spaces can significantly enhance water management. This research provides insights into effective water cycle management methods tailored to specific urban land uses, laying a foundation for future urban planning and sustainable development.

Spring Runoff Simulation of Snow-Dominant Catchment in Steppe Regions: A Comparison Study of Lumped Conceptual Models

Spring Runoff Simulation of Snow-Dominant Catchment in Steppe Regions: A Comparison Study of Lumped Conceptual Models

Assessing the Impact of Rainfall Inputs on Short-Term Flood Simulation with Cell2Flood: A Case Study of the Waryong Reservoir Basin

This study explored the impacts of various rainfall input types on short-term runoff simulations using the Cell2Flood model in the Waryong Reservoir Basin, South Korea. Six types of rainfall data were assessed: on-site gauge measurements, spatially interpolated data from 39 Automated Synoptic Observing System (ASOS) and 117 Automatic Weather System (AWS) stations using inverse distance weighting (IDW), and Hybrid Surface Rainfall (HSR) data from the Korea Meteorological Administration. The choice of rainfall input significantly affected model accuracy across the three rainfall events. The point-gauged ASOS (P-ASOS) data demonstrated the highest reliability in capturing the observed rainfall patterns, with Pearson’s r values of up to 0.84, whereas the radar-derived HSR data had the lowest correlations (Pearson’s r below 0.2), highlighting substantial discrepancies. For runoff simulation, the P-ASOS and ASOS-AWS combined interpolated dataset (R-AWS) achieved relatively accurate predictions, with P-ASOS and R-AWS exhibiting Normalized Peak Error (NPE) values of approximately 0.03 and Peak Time Error (PTE) within 20 min. In contrast, the HSR data produced large errors, with NPE up to 4.66 and PTE deviations exceeding 200 min, indicating poor temporal accuracy. Although input-specific calibration improved performance, significant errors persisted because of the inherent uncertainty of rainfall data. These findings underscore the importance of selecting and calibrating appropriate rainfall inputs to enhance the reliability of short-term flood modeling, particularly in ungauged and data-sparse basins.

A Novel Approach for Identifying and Assessing MOR‐Based CMIP6 Model for Hydrological Analysis in an Ungauged Watershed

ABSTRACTThe identification of the onset and retreat dates of the monsoon season is a crucial and intricate phenomenon, given its annual spatiotemporal variability. The monsoon season contributes significantly to rainfall, replenishing water sources and hydrological systems but causes hydrological extremes, especially for the high‐altitude watersheds in Southeast Asia. Global Circulation Model (GCM)‐Coupled Model Intercomparison Project Phase 6 (CMIP6)‐based rainfall and temperature data are helpful for adequately representing present and future climate scenarios. However, the usability of uncorrected GCM‐CMIP6 datasets needs to be assessed regionally. This study focuses on identifying the best‐suited GCM‐CMIP6 based on the monsoon onset (MO) and retreat (MR) dates along with other climatological temporal parameters. A numerical definition for MO and MR has been formulated to find the best‐suited GCM‐CMIP6 (i.e., CMCC‐ESM2). In this context, runoff simulation is carried out using the meteorological inputs of the monsoon onset‐retreat (MOR)‐based best‐suited GCM to evaluate its usability. A multi‐model simulation approach has been carried out for runoff estimation based on observed datasets to find a better‐suited hydrological model. The proposed overall methodology is tested in a hydrological extreme‐prone ungauged watershed (i.e., Ranikhola). CMCC‐ESM2 and SSP2‐4.5 has been identified as best‐suited SSP based on statistical evolution (R2 [0.693], NSE [0.662] and RSR [0.581]) for future daily runoff prediction. Future hydrological analysis shows that the average monsoon peak runoff magnitude will increase from the calibrated period (2015–2020) by 18.01% in the coming years (2021–2049).

Impacts of LULC changes on runoff from rivers through a coupled SWAT and BiLSTM model: A case study in Zhanghe River Basin, China

Changes in river runoff have a significant impact on the sustainable use of water resources in a watershed, and these changes are closely linked to variations in land use/land cover (LULC). This research explores an innovative approach in the Zhang River Basin (ZRB), China, by coupling a concept-based hydrological model, the Soil and Water Assessment Tool (SWAT), with a deep-learning model, the Bidirectional Long Short-Term Memory Network (Bi-LSTM), to improve the accuracy of river runoff simulations. By analyzing LULC changes in 2002, 2012, and 2022, this study developed three SWAT models and three coupled SWAT-BiLSTM models to quantitatively assess the impacts of these changes on river runoff through eight LULC scenarios. The findings revealed significant LULC changes from 2002 to 2022, with cropland and grassland areas decreasing while forest and urban land areas increased. The total area of grassland, forest, and cropland made up over 93 % of the basin, indicating active land type conversions. Calibration and validation results demonstrated that the SWAT-BiLSTM model outperformed the conventional SWAT model, yielding higher accuracy in runoff simulations. Specifically, the SWAT-BiLSTM model achieved R2 values of 0.89 and 0.90 during calibration and validation, compared to the SWAT model's R2 values of 0.76 and 0.79. Scenario analyses indicated that expansions in farmland, grassland, and urban areas were correlated with increased river runoff, while an expansion in forested areas led to reduced runoff. Notably, urban land changes had the most pronounced impact on runoff, emphasizing the need for careful runoff management and flood risk mitigation in urban planning. By combining SWAT and Bi-LSTM models, this study provides an innovative assessment of the impact of LULC changes on water resources in the ZRB. The results offer valuable insights for water resource management, LULC optimization, and flood risk management, highlighting the potential application of deep learning techniques in hydrological simulation. This research serves as a scientific basis for policy-making and sustainable land use planning in the ZRB and similar regions.

Study on Runoff Simulation with Multi-source Precipitation Information Fusion Based on Multi-model Ensemble

Study on Runoff Simulation with Multi-source Precipitation Information Fusion Based on Multi-model Ensemble

Snow depth retrieval method for PolSAR data using multi-parameters snow backscattering model

Snow depth (SD) is a crucial property of snow, its spatial and temporal variation is important for global change, snowmelt runoff simulation, disaster prediction, and freshwater storage estimation. Polarimetric Synthetic Aperture Radar (PolSAR) can precisely describe the backscattering of the target and emerge as an effective tool for SD retrieval. The backscattering component of dry snow is mainly composed of volume scattering from the snowpack and surface scattering from the snow-ground interface. However, the existing method for retrieving SD using PolSAR data has the problems of over-reliance on in-situ data and ignoring surface scattering from the snow-ground interface. We proposed a novel SD retrieval method for PolSAR data by fully considering the primary backscattering components of snow and through multi-parameter estimation to solve the snow backscattering model. Firstly, a snow backscattering model was formed by combining the small permittivity volume scattering model and the Michigan semi-empirical surface scattering model to simulate the different scattering components of snow, and the corresponding backscattering coefficients were extracted using the Yamaguchi decomposition. Then, the snow permittivity was calculated through generalized volume parameters and the extinction coefficient was further estimated through modeling. Finally, the snow backscattering model was solved by these parameters to retrieve SD. The proposed method was validated by Ku-band UAV SAR data acquired in Altay, Xinjiang, and the accuracy was evaluated by in-situ data. The correlation coefficient, root mean square error, and mean absolute error are 0.80, 4.49 cm, and 3.95 cm, respectively. Meanwhile, the uncertainties generated by different SD, model parameters estimation, solution method, and underlying surface are analyzed to enhance the generality of the proposed method.

Modelacion del impacto del deshielo en lo caudales en la cuenca del Lago Mansfiled Hollow en Connecticut, USA

Storm runoff predictions are essential for minimizing flood hazards and increasing resilience to extreme weather events. In this study, an analysis was conducted to simulate snowmelt runoff in the Mansfield Hollow Lake Watershed, which is a tributary of the Thames River watershed in Connecticut, New England. The United States Army Corp of Engineers (USACE) model HEC-HMS was applied to simulate snowmelt runoff during the winter-spring of 2010 and 2019. The Mansfield Hollow Lake Watershed is composed of three main tributaries, namely the Fenton, Mount Hope, and Natchaug Rivers. These runoff simulations and the watershed response to snowmelt are crucial for evaluating the potential impacts of watershed management decisions, particularly during high-flow periods. The HEC-HMS model was calibrated during the 2010 event and validated for the 2019 events. The study found that for the snow storms during 2010 and 2019 events, HEC-HMS model provided highly accurate predictions of snowmelt runoff with R-squared and, Nash - Sutcliffe correlation values exceeding 0.76. These findings highlight the efficacy of HEC-HMS model for simulating snowmelt runoff and demonstrate the utility of such model in predicting and managing flood risks. The results of this study provide valuable insights into the potential impacts of snowmelt runoff and will inform future watershed management decisions in the Mansfield Hollow Lake Watershed and similar regions.

Enhanced runoff simulation by precise capture of snowmelt variation signals with satellite-based snow products in a high-elevation basin

Hydrological models stand as a pivotal instrument for runoff simulation, while encountering notable uncertainties in output due to intricate terrain conditions and limited ground-based observations, especially in high-elevation basins. Leveraging satellite-based images presents a promising avenue for deciphering the hydrological model’s state variables. In pursuit of refining runoff simulation, this study developed a two-step enhancement strategy, which involved (1) calibrating snow-related parameters in the hydrological model using remotely sensed snow cover area (SCA) and snow water equivalent (SWE) and (2) capturing snowmelt runoff through the hydrological model and image-based products. Coupled with the Variable Infiltration Capacity (VIC) model, we adopted this strategy as a case study in the Dadu River Basin, China. The results indicated (1) the daily Nash-Sutcliffe Efficiency (NSE) of runoff simulation reached 0.84 by the enhancement strategy, markedly surpassing the strategy reliant on soil parameters with a single calibrated reference (daily NSE of 0.66), and exhibited comparability to the strategy incorporating snow parameters but calibrated solely based on observed discharge (daily NSE of 0.83). (2) the enhancement strategy demonstrated hydrological consistency with snowmelt information derived from imagery. Specifically, multi-year average contributions of model’s and image-based snowmelt calculations were 31.8% and 33.4%, respectively. Additionally, simulated dates of snow accumulation and ablation appeared an average deviation of approximately one week compared to the imagery results. This study elucidates a potential methodological approach for offering valuable insights into hydrological processes within analogous high-elevation basins globally.

Evaluation of SWAT Model in Runoff Simulation Using Rainfall and Temperature Derived From Satellite Images

Evaluation of SWAT Model in Runoff Simulation Using Rainfall and Temperature Derived From Satellite Images

Distributed modelling of snow and ice melt in the Naltar Catchment, Upper Indus basin

Energy balance distributed modelling in High Mountain Asia (HMA) is important to examine glaciological and hydrological processes and assess changes in streamflow in the current and future climate. In this study, the Physically based Distributed Snow Land and Ice Model (PDSLIM) using detailed observed meteorological data at hourly scale is employed to simulate the hydrological response of the Naltar catchment, 242.62 km2 in size, (in the Karakoram region in Pakistan to simulate its glaciers’ mass balance as well as daily runoff. The results exhibited overall satisfactory performance in terms of coefficient of determination (R2 = 0.96) and Nash-Sutcliffe Efficiency (NSE=0.95) modelled against satellite-based snow cover areas, for internal model verification, in eight years. The results of runoff simulations compared for external model verification, with observed daily discharge resulted in NSE 0.90 and 0.89 for calibration and validation period respectively. Flow composition analysis revealed that the streamflow regime of Naltar catchment is composed to 40 % by glacier runoff, 42 % by sub-surface runoff and 18 % by surface runoff. The eight year mean value of net mass balance exhibited a slightly negative mass balance (−0.810 ± 0.31 m w.e. a-1) less pronounced than that observed globally in several continental glaciers distinct from the Greenland and Antarctic ice sheets in the current climate. It seems that the ‘Karakoram anomaly’, i.e. the balanced to slightly positive glacier budgets observed in the region in the recent decades, a unique dynamics worldwide, has a moderate impact in the central Karakoram. Overall, the distributed energy-balance model PDSLIM, so far tested in the Alps, results to be a suitable tool to estimate energy and mass balance in the glacierized catchments of Karakoram and Himalaya and to better understand snow and ice melt runoff dynamics and floods in highly complex and glacierized mountain basins in the current and, in our research perspective, in the future climate.

Combining Multiple Machine Learning Methods Based on CARS Algorithm to Implement Runoff Simulation

Runoff forecasting is crucial for water resource management and flood safety and remains a central research topic in hydrology. Recent advancements in machine learning provide novel approaches for predicting runoff. This study employs the Competitive Adaptive Reweighted Sampling (CARS) algorithm to integrate various machine learning models into a data-driven rainfall–runoff simulation model. We compare the forecasting performance of different machine learning models to improve rainfall–runoff prediction accuracy. This study uses data from the Maduwang hydrological station in the Bahe river basin, which contain 12 measured flood events from 2000 to 2010. Historical runoff and areal mean rainfall serve as model inputs, while flood data at different lead times are used as model outputs. Among the 12 flood events, 9 are used as the training set, 2 as the validation set, and 1 as the testing set. The results indicate that the CARS-based machine learning model effectively forecasts floods in the Bahe River basin. Under the prediction period of 1 to 6 h, the model achieves high forecasting accuracy, with the average NSE ranging from 0.7509 to 0.9671 and the average R2 ranging from 0.8397 to 0.9413, though the accuracy declines to some extent as the lead time increases. The model accurately predicts peak flow and performs well in forecasting high flow and recession flows, though peak flows are somewhat underestimated for longer lead times. Compared to other machine learning models, the SVR model has the highest average RMSE of 0.942 for a 1–6 h prediction period. It exhibits the smallest deviation among low-, medium-, and high-flow curves, with the lowest NRMSE values across training, validation, and test sets, demonstrating better simulation performance and generalization capability. Therefore, the machine learning model based on CARS feature selection can serve as an effective method for flood forecasting. The related findings provide a new forecasting method and scientific decision-making basis for basin flood safety.

Runoff simulation modeling method integrating spatial element dynamics and neural network for remote sensing precipitation data

Effective management of water resources amidst global climate change necessitates advancements in runoff simulation techniques. This study introduces a novel methodology that integrates a spatial element dynamic model with a neural network designed for satellite precipitation product (SPP) and remote sensing data to address the complex dynamics of precipitation-runoff processes. Implemented in the Songhua River basin of China, the methodology features a runoff estimation module that operates across both high-resolution and SPP grid scales. This module synergizes spatial attributes with fuzzy classification for refined confluence optimization, and yielding simulations through long short-term memory network. The results indicate that the proposed multi-scale production-confluence scheme achieved Nash-Sutcliffe efficiency (NSE) scores of 0.87 and 0.84 during calibration (2011–2015) and validation (2017 and 2018) periods, respectively, representing an 8.15% improvement in peak percentage error (PPE10%) over conventional hydrologic models. Additionally, the scheme demonstrated superior simulation performance at tributary sites set compared to main channel sites set, with differences of 0.08 and 0.07 in NSE and Kling-Gupta efficiency (KGE), respectively. The introduction of the surface Ⅱ operator for adjusting precipitation intensity effectively improved runoff estimation accuracy, reducing the average error by 25.48% across both site sets. This study provides a scientific foundation for applying SPP in hydrological modeling and offers decision support for managing water resources.

An increase in the spatial extent of European floods over the last 70 years

Abstract. Floods regularly cause substantial damage worldwide. Changing flood characteristics, e.g., due to climate change, pose challenges to flood risk management. The spatial extent of floods is an important indicator of potential impacts, as consequences of widespread floods are particularly difficult to mitigate. The highly uneven station distribution in space and time, however, limits the ability to quantify flood characteristics and, in particular, changes in flood extents over large regions. Here, we use observation-driven routed runoff simulations over the last 70 years in Europe from a state-of-the-art hydrological model (the mesoscale Hydrologic Model – mHM) to identify large spatiotemporally connected flood events. Our identified spatiotemporal flood events compare well against an independent flood impact database. We find that flood extents increase by 11.3 % on average across Europe. This increase occurs over most of Europe, except for parts of eastern and southwestern Europe. Over northern Europe, the increase in flood extent is mainly driven by the overall increase in flood magnitude caused by increasing precipitation and snowmelt. In contrast, the increasing trend in flood extent over central Europe can be attributed to an increase in the spatial extent of heavy precipitation. Overall, our study illustrates the opportunities to combine long-term consistent regional runoff simulations with a spatiotemporal flood detection algorithm to identify large-scale trends in key flood characteristics and their drivers. The detected change in flood extent should be considered in risk assessments as it may challenge flood control and water resource management.

Runoff Prediction of Tunxi Basin under Projected Climate Changes Based on Lumped Hydrological Models with Various Model Parameter Optimization Strategies

Runoff is greatly influenced by changes in climate conditions. Predicting runoff and analyzing its variations under future climates are crucial for ensuring water security, managing water resources effectively, and promoting sustainable development within the catchment area. As the key step in runoff modeling, the calibration of hydrological model parameters plays an important role in models’ performance. Identifying an efficient and reliable optimization algorithm and objective function continues to be a significant challenge in applying hydrological models. This study selected new algorithms, including the strategic random search (SRS) and sparrow search algorithm (SSA) used in hydrology, gold rush optimizer (GRO) and snow ablation optimizer (SAO) not used in hydrology, and classical algorithms, i.e., shuffling complex evolution (SCE-UA) and particle swarm optimization (PSO), to calibrate the two-parameter monthly water balance model (TWBM), abcd, and HYMOD model under the four objective functions of the Kling–Gupta efficiency (KGE) variant based on knowable moments (KMoments) and considering the high and low flows (HiLo) for monthly runoff simulation and future runoff prediction in Tunxi basin, China. Furthermore, the identified algorithm and objective function scenario with the best performance were applied for runoff prediction under climate change projections. The results show that the abcd model has the best performance, followed by the HYMOD and TWBM models, and the rank of model stability is abcd > TWBM > HYMOD with the change of algorithms, objective functions, and contributing calibration years in the history period. The KMoments based on KGE can play a positive role in the model calibration, while the effect of adding the HiLo is unstable. The SRS algorithm exhibits a faster, more stable, and more efficient search than the others in hydrological model calibration. The runoff obtained from the optimal model showed a decrease in the future monthly runoff compared to the reference period under all SSP scenarios. In addition, the distribution of monthly runoff changed, with the monthly maximum runoff changing from June to May. Decreases in the monthly simulated runoff mainly occurred from February to July (10.9–56.1%). These findings may be helpful for the determination of model parameter calibration strategies, thus improving the accuracy and efficiency of hydrological modeling for runoff prediction.

Runoff Simulation in Data-Scarce Alpine Regions: Comparative Analysis Based on LSTM and Physically Based Models

Runoff simulation is essential for effective water resource management and plays a pivotal role in hydrological forecasting. Improving the quality of runoff simulation and forecasting continues to be a highly relevant research area. The complexity of the terrain and the scarcity of long-term runoff observation data have significantly limited the application of Physically Based Models (PBMs) in the Qinghai–Tibet Plateau (QTP). Recently, the Long Short-Term Memory (LSTM) network has been found to be effective in learning the dynamic hydrological characteristics of watersheds and outperforming some traditional PBMs in runoff simulation. However, the extent to which the LSTM works in data-scarce alpine regions remains unclear. This study aims to evaluate the applicability of LSTM in alpine basins in QTP, as well as the simulation performance of transfer-based LSTM (T-LSTM) in data-scarce alpine regions. The Lhasa River Basin (LRB) and Nyang River Basin (NRB) were the study areas, and the performance of the LSTM model was compared to that of PBMs by relying solely on the meteorological inputs. The results show that the average values of Nash–Sutcliffe efficiency (NSE), Kling–Gupta efficiency (KGE), and Relative Bias (RBias) for B-LSTM were 0.80, 0.85, and 4.21%, respectively, while the corresponding values for G-LSTM were 0.81, 0.84, and 3.19%. In comparison to a PBM- the Block-Wise use of TOPMEDEL (BTOP), LSTM has an average enhancement of 0.23, 0.36, and −18.36%, respectively. In both basins, LSTM significantly outperforms the BTOP model. Furthermore, the transfer learning-based LSTM model (T-LSTM) at the multi-watershed scale demonstrates that, when the input data are somewhat representative, even if the amount of data are limited, T-LSTM can obtain more accurate results than hydrological models specifically calibrated for individual watersheds. This result indicates that LSTM can effectively improve the runoff simulation performance in alpine regions and can be applied to runoff simulation in data-scarce regions.