Search

Recently, deep learning (DL) models, especially those based on long short-term memory (LSTM), have demonstrated their superior ability in resolving sequential data problems. This study investigated the performance of six models that belong to the supervised learning category to evaluate the performance of DL models in terms of streamflow forecasting. They include a feed-forward neural network (FFNN), a convolutional neural network (CNN), and four LSTM-based models. Two standard models with just one hidden layer—LSTM and gated recurrent unit (GRU)—are used against two more complex models—the stacked LSTM (StackedLSTM) model and the Bidirectional LSTM (BiLSTM) model. The Red River basin—the largest river basin in the north of Vietnam—was adopted as a case study because of its geographic relevance since Hanoi city—the capital of Vietnam—is located downstream of the Red River. Besides, the input data of these models are the observed data at seven hydrological stations on the three main river branches of the Red River system. This study indicates that the four LSTM-based models exhibited considerably better performance and maintained stability than the FFNN and CNN models. However, the complexity of the StackedLSTM and BiLSTM models is not accompanied by performance improvement because the results of the comparison illustrate that their respective performance is not higher than the two standard models—LSTM and GRU. The findings of this study present that LSTM-based models can reach impressive forecasts even in the presence of upstream dams and reservoirs. For the streamflow-forecasting problem, the LSTM and GRU models with a simple architecture (one hidden layer) are sufficient to produce highly reliable forecasts while minimizing the computation time.

Climate change currently affects the resilience and aquatic ecosystem. Climate change alters rainfall patterns which have a great impact on river flow. Annual flooding is an important hydrological characteristic of the Mekong River Basin (MRB) and it drives the high productivity of the ecosystem and biodiversity in the Tonle Sap floodplain and the Mekong Delta. This study aims to assess the impacts of climate change on river flow in the MRB and flood inundation in the Lower Mekong Basin (LMB). The changing impacts were assessed by a two-dimensional rainfall-runoff and inundation model (RRI model). The present climate (1979–2003) and future projected climate (2075–2099) datasets from MRI-AGCM3.2H and MRI-AGCM3.2S models were applied with a linear scaling bias correction method before input into the RRI model. The results of climate change suggested that flood magnitude in the LMB will be severer than the present climate by the end of the twenty-first century. The increment of precipitation between 6.6 and 14.2% could lead to increase extreme flow (Q5) 13–30%, peak inundation area 19–43%, and peak inundation volume 24–55% in the LMB for ranging of Representative Concentration Pathways (RCP) and sea surface temperature (SST) scenarios while there is no significant change on peak flood timing.

Precipitation is a crucial component of the water cycle and plays a key role in hydrological processes. Recently, satellite-based precipitation products (SPPs) have provided grid-based precipitation with spatiotemporal variability. However, SPPs contain a lot of uncertainty in estimated precipitation, and the spatial resolution of these products is still relatively coarse. To overcome these limitations, this study aims to generate new grid-based daily precipitation based on a combination of rainfall observation data with multiple SPPs for the period of 2003–2017 across South Korea. A Random Forest (RF) machine-learning algorithm model was applied for producing a new merged precipitation product. In addition, several statistical linear merging methods have been adopted to compare with the results achieved from the RF model. To investigate the efficiency of RF, rainfall data from 64 observed Automated Synoptic Observation System (ASOS) installations were collected to analyze the accuracy of products through several continuous as well as categorical indicators. The new precipitation values produced by the merging procedure generally not only report higher accuracy than a single satellite rainfall product but also indicate that RF is more effective than the statistical merging method. Thus, the achievements from this study point out that the RF model might be applied for merging multiple satellite precipitation products, especially in sparse region areas.

Recently, machine learning (ML) is being applied to various fields, including hydrology and hydraulics. The numerical models based on ML algorithms have been widely used for forecasting water levels or flowrate in different timescales. Especially in estuary areas where the hydrodynamic regime becomes complicated, the water level forecast information in this area plays an essential role in the operation of tidal sluices. This study proposes an efficient approach using an ML model, long short-term memory (LSTM), to predict short-term water levels in tidal sluice gates from 6 to 48 hours ahead. The An Tho culvert located in the Bac Hung Hai irrigation system, the most extensive irrigation system in Vietnam, was selected as a case study station. The high accuracy of predictive results reveals LSTM models' effectiveness in different forecasting scenarios. In the first scenario using just water level data at the prediction station, the Kling-Gupta efficiency (KGE) coefficient ranges from nearly 0.89 to 0.96. Meanwhile, in the second scenario, the combination of observed data of three gauge stations exhibited better performance with KGE coefficients ranging from just under 0.93 to 0.98 for eight forecasted cases. The findings of this study highlight the performance of LSTM models in providing high-accuracy short-period water level forecasts for areas near estuaries. These obtained results can play a vital role in the management and operation of tidal sluices in the Bac Hung Hai irrigation system, as well as a reference for the operation of other irrigation systems around the world.

In Korea, approximately 70% of the country is mountainous, with steep slopes and heavy rainfall in summer from June to September. Korea is classified as a high-risk country for soil erosion, and the rate of soil erosion is rapidly increasing. In particular, the operation of Doam dam was suspended in 2001 because of water quality issues due to severe soil erosion from the upstream areas. In spite of serious dam sediment problems in this basin, in-depth studies on the origin of sedimentation using physic-based models have not been conducted. This study aims to analyze the spatial distribution of net erosion during typhoon events using a spatially distributed physics-based erosion model and to improve the model based on a field survey. The spatially uniform erodibility constants of the surface flow detachment equation in the original erosion model were replaced by land use erodibility constants based on benchmarking experimental values to reflect the effect of land use on net erosion. The results of the upgraded model considering spatial erodibility show a significant increase in soil erosion in crop fields and bare land, unlike the simulation results before model improvement. The total erosion and deposition for Typhoon Maemi in 2003 were 36,689.0 and 9893.3 m3, respectively, while the total erosion and deposition for Typhoon Rusa in 2002 were 142,476.6 and 44,806.8 m3, respectively, despite about twice as much rainfall and 1.2 times as high rainfall intensity. However, there is a limitation in quantifying the sources of erosion in the study watershed, since direct comparison of the simulated net erosion with observed spatial information from aerial images, etc., is impossible due to nonperiodic image photographing. Therefore, continuous monitoring of not only sediment yield but also periodic spatial detection on erosion and deposition is critical for reducing data uncertainty and improving simulation accuracy.

Spatiotemporal precipitation data is one of the essential components in modeling hydrological problems. Although the estimation of these data has achieved remarkable accuracy owning to the recent advances in remote-sensing technology, gaps remain between satellite-based precipitation and observed data due to the dependence of precipitation on the spatiotemporal distribution and the specific characteristics of the area. This paper presents an efficient approach based on a combination of the convolutional neural network and the autoencoder architecture, called the convolutional autoencoder (ConvAE) neural network, to correct the pixel-by-pixel bias for satellite-based products. The two daily gridded precipitation datasets with a spatial resolution of 0.25° employed are Asian Precipitation-Highly Resolved Observational Data Integration towards Evaluation (APHRODITE) as the observed data and Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks-Climate Data Record (PERSIANN-CDR) as the satellite-based data. Furthermore, the Mekong River basin was selected as a case study, because it is one of the largest river basins, spanning six countries, most of which are developing countries. In addition to the ConvAE model, another bias correction method based on the standard deviation method was also introduced. The performance of the bias correction methods was evaluated in terms of the probability distribution, temporal correlation, and spatial correlation of precipitation. Compared with the standard deviation method, the ConvAE model demonstrated superior and stable performance in most comparisons conducted. Additionally, the ConvAE model also exhibited impressive performance in capturing extreme rainfall events, distribution trends, and described spatial relationships between adjacent grid cells well. The findings of this study highlight the potential of the ConvAE model to resolve the precipitation bias correction problem. Thus, the ConvAE model could be applied to other satellite-based products, higher-resolution precipitation data, or other issues related to gridded data.

In establishing adequate climate change policies regarding water resource development and management, the most essential step is performing a rainfall-runoff analysis. To this end, although several physical models have been developed and tested in many studies, they require a complex grid-based parameterization that uses climate, topography, land-use, and geology data to simulate spatiotemporal runoff. Furthermore, physical rainfall-runoff models also suffer from uncertainty originating from insufficient data quality and quantity, unreliable parameters, and imperfect model structures. As an alternative, this study proposes a rainfall-runoff analysis system for the Kratie station on the Mekong River mainstream using the long short-term memory (LSTM) model, a data-based black-box method. Future runoff variations were simulated by applying a climate change scenario. To assess the applicability of the LSTM model, its result was compared with a runoff analysis using the Soil and Water Assessment Tool (SWAT) model. The following steps (dataset periods in parentheses) were carried out within the SWAT approach: parameter correction (2000&ndash, 2005), verification (2006&ndash, 2007), and prediction (2008&ndash, 2100), while the LSTM model went through the process of training (1980&ndash, 2100). Globally available data were fed into the algorithms, with the exception of the observed discharge and temperature data, which could not be acquired. The bias-corrected Representative Concentration Pathways (RCPs) 4.5 and 8.5 climate change scenarios were used to predict future runoff. When the reproducibility at the Kratie station for the verification period of the two models (2006&ndash, 2007) was evaluated, the SWAT model showed a Nash&ndash, Sutcliffe efficiency (NSE) value of 0.84, while the LSTM model showed a higher accuracy, NSE = 0.99. The trend analysis result of the runoff prediction for the Kratie station over the 2008&ndash, 2100 period did not show a statistically significant trend for neither scenario nor model. However, both models found that the annual mean flow rate in the RCP 8.5 scenario showed greater variability than in the RCP 4.5 scenario. These findings confirm that the LSTM runoff prediction presents a higher reproducibility than that of the SWAT model in simulating runoff variation according to time-series changes. Therefore, the LSTM model, which derives relatively accurate results with a small amount of data, is an effective approach to large-scale hydrologic modeling when only runoff time-series are available.

The Mekong River plays an extremely important role in Southeast Asia. Flowing through six countries, including China, Myanmar, Thailand, Laos PDR, Cambodia, and Vietnam, it is a site of great biological and ecological diversity and the habitat of numerous species of fish. It also supports a very large population that lives along the river basin. Therefore, much attention has been focused on the giant Mekong River Basin, particularly, its soil erosion and sedimentation problems. In fact, many methods have been used to calculate and simulate these problems. However, in the case of the Mekong River Basin, the available data is limited because of the extreme size of the area (about 795,000 ㎢) and lack of equipment systems in the countries through which the Mekong River flows. In this study, we applied the Universal Soil Loss Equation (USLE) model in a GIS (Geographic Information System) framework to calculate the amount of soil erosion and sediment load during the selected period, from 1951 to 2007. The result points out dangerous areas, such as the Upper Mekong River Basin and 3S Basin (containing the Sekong, Sesan, and Srepok Rivers) that are suffering the serious consequences of soil erosion problems. Moreover, the present model is also useful for supporting river basin management in the implementation of sustainable management practices in the Mekong River Basin and other basins.

An accurate landslide-susceptibility assessment is fundamental for preventing landslides and minimizing damage. In this study, a new time-variant slope-stability (TiVaSS) model for landslide prediction is developed. A three-dimensional (3D) subsurface flow model is coupled with the infinite slope-stability model to consider the effect of horizontal water movement in the subsurface. A 3D Richards' equation is solved numerically for the subsurface flow. To overcome the massive computational requirements of the 3D subsurface flow module, partially implicit temporal discretization and the simplification of first-order spatial discretization are proposed and applied in TiVaSS. A graphical user interface and two-dimensional data visualization are supported in TiVaSS. The model is applied to a 2011źMt. Umyeon landslide in the Republic of Korea, and its overall performance is satisfactory. A new time-variant slope-stability model for predicting rainfall-induced landslides is developed.The three-dimensional subsurface unsaturated flow and an infinite slope-stability model are combined.A graphical user interface including input and output visualization is supported.

A rainfall induced slope failure is a common natural hazard in mountainous areas worldwide. Sudden and rapid failures which have a high possibility of occurrence in a steep slope are always the most dangerous due to their suddenness and high velocities. Based on a series of experiments this study aimed to determine a critical angle which could be considered as an approximate threshold for a sudden failure. The experiments were performed using 0.42 mm mean grain size sand in a 200 cm long, 60 cm wide and 50 cm deep rectangular flume. A numerical model was created by integrating a 2D seepage flow model and a 2D slope stability analysis model to predict the failure surface and the time of occurrence. The results showed that, the failure mode for the entire material will be sudden for slopes greater than 67°; in contrast the failure mode becomes retrogressive. There is no clear link between the degree of saturation and the mode of failure. The simulation results in considering matric suction showed good matching with the results obtained from experiment. A subsequent discarding of the matric suction effect in calculating safety factors will result in a deeper predicted failure surface and an incorrect predicted time of occurrence.

Many researchers have evaluated the influence of vegetation cover on slope stability. However, due to the extensive variety of site conditions and vegetation types, different studies have often provided inconsistent results, especially when evaluating in different regions. Therefore, additional studies need to be conducted to identify the positive impacts of vegetation cover for slope stabilization. This study used the Transient Rainfall Infiltration and Grid-based Regional Slope-stability Model (TRIGRS) to predict the occurrence of landslides in a watershed in Jinbu-Myeon, Pyeongchang-gun, Korea. The influence of vegetation cover was assessed by spatially and temporally comparing the predicted landslides corresponding to multiple trials of cohesion values (which include the role of root cohesion) and real observed landslide scars to back-calculate the contribution of vegetation cover to slope stabilization. The lower bound of cohesion was defined based on the fact that there are no unstable cells in the raster stability map at initial conditions, and the modified success rate was used to evaluate the model performance. In the next step, the most reliable value representing the contribution of vegetation cover in the study area was applied for landslide assessment. The analyzed results showed that the role of vegetation cover could be replaced by increasing the soil cohesion by 3.8 kPa. Without considering the influence of vegetation cover, a large area of the studied watershed is unconditionally unstable in the initial condition. However, when tree root cohesion is taken into account, the model produces more realistic results with about 76.7% of observed unstable cells and 78.6% of observed stable cells being well predicted.

This paper proposes an extended model evaluation method that considers not only the model performance but also the model structure and parameter uncertainties in hydrologic modeling. A simple reservoir model (SFM) and distributed kinematic wave models (KWMSS1 and KWMSS2 using topography from 250-m, 500-m, and 1-km digital elevation models) were developed and assessed by three evaluative criteria for model performance, model structural stability, and parameter identifiability. All the models provided acceptable performance in terms of a global response, but the simpler SFM and KWMSS1 could not accurately represent the local behaviors of hydrographs. Moreover, SFM and KWMSS1 were structurally unstable; their performance was sensitive to the applied objective functions. On the other hand, the most sophisticated model, KWMSS2, performed well, satisfying both global and local behaviors. KMSS2 also showed good structural stability, reproducing hydrographs regardless of the applied objective functions; however, superior parameter identifiability was not guaranteed. A number of parameter sets could result in indistinguishable hydrographs. This result indicates that while making hydrologic models complex increases its performance accuracy and reduces its structural uncertainty, the model is likely to suffer from parameter uncertainty.

This paper addresses the effects of extreme rainfall on the stability of cut slopes in Yen Bai city, Northern Viet Nam. In this area, natural slopes are excavated to create places for infrastructures and buildings. Cut slopes are usually made without proper site investigations; the design is mostly based on experience. In recent years, many slope failures have occurred along these cuts especially in rainy seasons, resulting in properties damaged and loss of lives. To explain the reason that slope failure often happens during rainy seasons, this research analyzed the influence of extreme rainfalls, initial ground conditions, and soil permeability on the changes of pore water pressure within the typical slope, thereafter determining the impact of these changes on the slope stability factor of safety. The extreme rainfalls were selected based on all of the rainfalls triggering landslide events that have occurred over the period from 1960 to 2009. The factor of safety (FS) was calculated using Bishop's simplified method. The results show that when the maximum infiltration capacity of the slope top soil is less than the rainfall intensity, slope failures may occur 14 hours after the rain starts. And when this happens, the rainfall duration is the deciding factor that affects the slope FS values. In short, cut slopes in Yen Bai may be stable in normal conditions after the excavation, but under the influence of tropical rain storms, their stability is always questionable.

This study aims to assess the slope stability variation of Jeonbuk drainage areas by RCM model outputs based on A1B climate change scenario and infinite slope stability model based on the previous research by Choi et al.(2013). For a large-scale slope stability analysis, we developed a GIS-based database regarding topographic, geologic and forestry parameters and also calculated daily maximum rainfall for the study period(1971~2100). Then, we assess slope stability variation of the 20 sub-catchments of Jeonbuk under the climate change scenario. The results show that the areal-average value of safety factor was estimated at 1.36(moderately stable) in spite of annual rainfall increase in the future. In addition, 7 sub-catchments became worse and 5 sub-catchments became better than the present period(1971~2000) in terms of safety factor in the future.

Catchments soil erosion, one of the most serious problems in the mountainous environment of the world, consists of a complex phenomenon involving the detachment of individual soil particles from the soil mass and their transport, storage and overland flow of rainfall, and infiltration. Sediment size distribution during erosion processes appear to depend on many factors such as rainfall characteristics, vegetation cover, hydraulic flow, soil properties and slope. This study involved laboratory flume experiments carried out under simulated rainfall in a 3.0 m long × 0.8 m wide × 0.7 m deep flume, set at 17° slope. Five experimental cases, consisting of twelve experiments using three different sediments with two different rainfall conditions, are reported. The experiments consisted of detailed observations of particle size distribution of the out-flow sediment. Sediment water mixture out-flow hydrograph and sediment mass out-flow rate over time, moisture profiles at different points within the soil domain, and seepage outflow were also reported. Moisture profiles, seepage outflow, and movement of overland flow were clearly found to be controlled by water retention function and hydraulic function of the soil. The difference of grain size distribution of original soil bed and the out-flow sediment was found to be insignificant in the cases of uniform sediment used experiments. However, in the cases of non-uniform sediment used experiments the outflow sediment was found to be coarser than the original soil domain. The results indicated that the sediment transport mechanism is the combination of particle segregation, suspension/saltation and rolling along the travel distance.

Flood forecasting is an essential requirement in integrated water resource management. This paper suggests a Long Short-Term Memory (LSTM) neural network model for flood forecasting, where the daily discharge and rainfall were used as input data. Moreover, characteristics of the data sets which may influence the model performance were also of interest. As a result, the Da River basin in Vietnam was chosen and two different combinations of input data sets from before 1985 (when the Hoa Binh dam was built) were used for one-day, two-day, and three-day flowrate forecasting ahead at Hoa Binh Station. The predictive ability of the model is quite impressive: The Nash&ndash, Sutcliffe efficiency (NSE) reached 99%, 95%, and 87% corresponding to three forecasting cases, respectively. The findings of this study suggest a viable option for flood forecasting on the Da River in Vietnam, where the river basin stretches between many countries and downstream flows (Vietnam) may fluctuate suddenly due to flood discharge from upstream hydroelectric reservoirs.

There has been an increase for the landslide areas and restoration expenses due, in large part, to the increased locally heavy rains caused by recent climate change as well as the reckless development. This study carried out a slope stability analysis by the application of distributed wetness index, using the GIS-based infinite slope stability model, which took the root cohesion effect into consideration, for part of Mt. Umyeon in Seoul, where landslide occurred in July 2011, in order to compensate the defects of existing analysis method, and subsequently compared its result with the case on the exploitation of lumped wetness index. In addition, this study estimated the distributed wetness index by methodology, applying three methods of specific catchment area calculation: single flow direction (SFD), multiple flow direction (MFD), and infinity flow direction (IFD), for catchment area, one of the variables of distributed wetness indices, and finally implemented a series of comparative analysis for slope stability by methodology. The simulation results showed that most unstable areas within the study site were dominantly located in cutting-area surroundings along with the residential area and the mountaintop and unstable areas of IFD and lumped wetness index method were similar while SFD and MFD provided smaller unstable areas than the two former methods.

This paper presents the effect of spatially-distributed rainfall on both rainfall-sediment-runoff and erosion or deposition in the experimental Cheoncheon catchment: upstream of Yongdam dam basin. The rainfall fields were generated by three rainfall interpolation techniques (Thiessen polygon: TP, Inverse Distance Weighting: IDW, Kriging) based only on ground gauges and two radar rainfall synthetic techniques (Gauge-Radar ratio: GR, Conditional Merging: CM). Each rainfall field was then assessed in terms of spatial feature and quantity and also used for rainfall-sediment-runoff and erosion-deposition simulation due to the spatial difference of rainfall fields. The results showed that all the interpolation methods based on ground gauges provided very similar hydrologic responses in spite of different spatial pattern of erosion and deposition while raw radar and GR rainfall fields led to underestimated and overestimated simulation results, respectively. The CM technique was acceptable to improve the accuracy of raw radar rainfall for hydrologic simulation even though it is more time consuming to generate spatially-distributed rainfall.

This study aims to demonstrate the relationship between various factors and soil erosion or deposition, simulated from distributed rainfall-sediment-runoff model applications. We selected area, overland flow length, local slope as catchment representative characteristics among many important geographic factors and also used the grid-based accumulated rainfall as a representative hydro-climatic factor to assess the effect of these two different types of factors on erosion and deposition. The study catchment was divided based on the Strahler`s stream order method for analysis of the relationship between area and erosion or deposition. Both erosion and deposition increased linearly as the catchment area became larger. Erosion occurred widely throughout the catchment, whereas deposition was observed at the grid-cells near the channel network with short overland flow lengths and mild slopes. In addition, the relationship results between grid-based accumulated rainfall and soil erosion or deposition showed that erosion increased gradually as rainfall amount increased, whereas deposition responded irregularly to variations in rainfall. Within the context of these results, it can be concluded that deposition is closely related with the geographic factors used in this study while erosion is significantly affected by rainfall.

The purpose of this research is to understand the change of runoff characteristics by estimated spatial rainfall. Therefore, this paper largely composed of two parts. First, we compared the simulated result according to estimation method, ID(Inverse Distance Method, ID2(Inverse Square Distance Method), and Kr(General Covariance Kriging Method), after letting miss rainfall data to the observed data. Second, we reviewed the runoff characteristics of the distributed runoff model according to the estimated spatial rainfall. On the basis of Yuseong water level station, we select the target basin as Gabchun watershed. We assumed 1 point or 2 point of the 6 rainfall gauge stations in watershed were missed. We applied the spatial rainfall distributed by Kr to Hy-GIS GRM, distributed runoff model. When 1 point rainfall data is missed, Kr is superior to others in point rainfall estimation and runoff estimation of Hy-GIS GRM. However, in case rainfall data of 2 points is missed, all of three methods did not give suitable result for them. In conclusion, Kr showed better applicability than other estimated methods if rainfall’s data less than 2 points is missed.Key Words :General covariance Kriging, HyGIS-GRM, Missing rainfall data

In hydrologic modeling, prediction uncertainty generally stems from various uncertainty sources associated with model structure, data, and parameters, etc. This study aims to assess the parameter uncertainty effect on hydrologic prediction results. For this objective, a distributed rainfall-sediment yield-runoff model, which consists of rainfall-runoff module for simulation of surface and subsurface flows and sediment yield module based on unit stream power theory, was applied to the mesoscale mountainous area (Cheoncheon catchment; 289.9 ). For parameter uncertainty evaluation, the model was calibrated by a multi-objective optimization algorithm (MOSCEM) with two different objective functions (RMSE and HMLE) and Pareto optimal solutions of each case were then estimated. In Case I, the rainfall-runoff module was calibrated to investigate the effect of parameter uncertainty on hydrograph reproduction whereas in Case II, sediment yield module was calibrated to show the propagation of parameter uncertainty into sedigraph estimation. Additionally, in Case III, all parameters of both modules were simultaneously calibrated in order to take account of prediction uncertainty in rainfall-sediment yield-runoff modeling. The results showed that hydrograph prediction uncertainty of Case I was observed over the low-flow periods while the sedigraph of high-flow periods was sensitive to uncertainty of the sediment yield module parameters in Case II. In Case III, prediction uncertainty ranges of both hydrograph and sedigraph were larger than the other cases. Furthermore, prediction uncertainty in terms of spatial distribution of erosion and deposition drastically varied with the applied model parameters for all cases.

Accelerated soil erosion due to extreme climate change, such as increased rainfall intensity, and human-induced environmental changes, is a widely recognized problem. Existing soil erosion models are generally based on the gross erosion concept to compute annual upland soil loss in tons per acre per year. However, such models are not suitable for event-based simulations of erosion and deposition in time and space. Recent advances in computer geographic information system (GIS) technologies have allowed hydrologists to develop physically based models, and the trend in erosion prediction is towards process-based models, instead of conceptually lumped models. This study aims to propose an effective and robust distributed rainfall-sediment yield-runoff model consisting of basic element modules: a rainfall-runoff module based on the kinematic wave method for subsurface and surface flow, and a runoff-sediment yield-runoff model based on the unit stream power method. The model was tested on the Cheoncheon catchment, upstream of the Yongdam dam using hydrological data for three extreme flood events due to typhoons. The model provided acceptable simulation results with respect to both discharge and sediment discharge even though the simulated sedigraphs were underestimated, compared to observations. The spatial distribution of erosion and deposition demonstrated that eroded sediment loads were deposited in the cells along the channel network, which have a short overland flow length and a gentle local slope while the erosion rate increased as rainfall became larger. Additionally, spatially heterogeneous rainfall intensity, dependant on Thiessen polygons, led to spatially-distinct erosion and deposition patterns.

A natural rainfall-runoff process is conceptualized by means of mathematical form reflecting a physical law and hydrologist’s perception and it is typically calibrated and verified based on streamflow data, the commonly used catchment response. The streamflow data is obviously required, but is not sufficient to identify a conceptual parameter of a hydrologic model since numerous parameter combinations can often result in quite similar numerical values in terms of objective function and even indistinguishable simulated hydrographs. One of the efficient techniques to resolve this impasse is a combination of the model with streamflow data augmented by other kinds of hydrological information relevant to the prediction. In this study, internal dynamic response of catchment, spatiotemporally separated components of hydrograph, to the various mimic parameter sets is presented. It is concluded that all of simulations based on plausible parameter combinations sampled by using deterministic and stochastic automatic optimization algorithms, are performed equally well in terms of both model performance measure and visual comparison between observed and simulated hydrograph while internal behaviors of catchment show totally different aspects. The spatiotemporal information presented here can be utilized as one of the complementary constraints capable of filtering out non-physical parameter set(s).

No matter how sophisticated and accurate hydrologic models may be, prediction uncertainty is unavoidable problem in rainfall-runoff modeling and it stems from various components. Therefore, it is essential to identify and reduce sources of uncertainty for more precise agreement between model prediction and observations in the real system. In the context of uncertainty issues, the problem of model structural uncertainty or stability is an issue of increasing interest in recent research. This paper examines a nature of model structural inadequacy using a single-objective global optimization method in hydrological modeling and proposes a framework to assess model structural stability through a comparison of two hydrologic models.