This study aimed to assess the suitability of the parameters of a physically based, distributed, grid-based rainfall-runoff model. We analyzed parameter sensitivity with a dataset of eight rainfall events that occurred in two catchments of South Korea, using the Sobol’ method. Parameters identified as sensitive responded adequately to the scale of the rainfall events and the objective functions employed. Parameter sensitivity varied depending on rainfall scale, even in the same catchment. Interestingly, for a rainfall event causing considerable runoff, parameters related to initial soil saturation and soil water movement played a significant role in low flow calculation and high flow calculation, respectively. The larger and steeper catchment exhibited a greater difference in parameter sensitivity between rainfall events. Finally, we found that setting an incorrect parameter range that is physically impossible can have a large impact on runoff simulation, leading to substantial uncertainty in the simulation results. The proposed analysis method and the results from our study can help researchers using a distributed rainfall-runoff model produce more reliable analysis results.

The hydrological model assessment and development (hydromad) modeling package is an R-based package that can be applied to simulate hydrological models and optimize parameters. As the hydromad package is incompatible with hydrological models outside the package, the parameters of such models cannot be directly optimized. Hence, we proposed a method of optimizing the hydrological-model parameters by combining the executable (EXE) file of the hydrological model with the shuffled complex evolution (SCE) algorithm provided by the hydromad package. A physically based, spatially distributed, grid-based rainfall–runoff model (GRM) was employed. By calibrating the parameters of the GRM, the performance of the model was found to be reasonable. Thus, the hydromad can be used to optimize the hydrological-model parameters outside the package if the EXE file of the hydrological model is available. The time required to optimize the parameters depends on the type of event, even for the same catchment area.

The choice of the computational time step (dt) value and the method for setting dt can have a bearing on the accuracy and performance of a simulation, and this effect has not been comprehensively researched across different simulation conditions. In this study, the effects of the fixed time step (FTS) method and the automatic time step (ATS) method on the simulated runoff of a distributed rainfall–runoff model were compared. The results revealed that the ATS method had less peak flow variability than the FTS method for the virtual catchment. In the FTS method, the difference in time step had more impact on the runoff simulation results than the other factors such as differences in the amount of rainfall, the density of the stream network, or the spatial resolution of the input data. Different optimal parameter values according to the computational time step were found when FTS and ATS were used in a real catchment, and the changes in the optimal parameter values were smaller in ATS than in FTS. The results of our analyses can help to yield reliable runoff simulation results.

This paper presents an algorithm that can efficiently simulate a grid-based one-dimensional distributed rainfall-runoff model by performing parallel computations using flow accumulation values for individual grid cells, which are calculated through an eight flow direction method. This parallel computation algorithm uses information about flow accumulation to automatically find parallel computation target grid cells within the overall area and perform parallel computation on the grid by unit. The Microsoft.NET Parallel class was used to apply and evaluate the parallel computation algorithm independently on two machines. The results showed that the time reduction effect of parallel computation differed for each target domain, because flow accumulation values varied depending on the domain. Parallel computation reduced computation time by around 40% to 78% in virtual domains and around 63% in the real domain compared to sequential computation. The results of this study can be utilized to reduce the computation time of distributed models.