Gaussian Process Regression and Cooperation Search Algorithm for Forecasting Nonstationary Runoff Time Series

: In the hydrology ﬁeld, hydrological forecasting is regarded as one of the most challenging engineering tasks, as runoff has signiﬁcant spatial–temporal variability under the inﬂuences of multiple physical factors from both climate events and human activities. As a well-known artiﬁcial intelligence tool, Gaussian process regression (GPR) possesses satisfying generalization performance but often suffers from local convergence and sensitivity to initial conditions in practice. To enhance its performance, this paper investigates the effectiveness of a hybrid GPR and cooperation search algorithm (CSA) model for forecasting nonstationary hydrological data series. The CSA approach avoids the premature convergence defect in GPR by effectively determining suitable parameter combinations in the problem space. Several traditional machine learning models are established to evaluate the validity of the proposed GPR-CSA method in three real-world hydrological stations of China. In the modeling process, statistical characteristics and expert knowledge are used to select input variables from the observed runoff data at previous periods. Different experimental results show that the developed GPR-CSA model can accurately predict nonlinear runoff and outperforms the developed traditional models in terms of various statistical indicators. Hence, a CSA-trained GPR model can provide satisfying training efﬁciency and robust simulation performance for runoff forecasting.


Introduction
Precise and timely runoff prediction is crucial for reducing flood damage, managing water resources and optimizing reservoir scheduling [1][2][3]. Moreover, multiple huge reservoirs in large rivers are put into operation in succession and, thus, the actual practical requirements for accurate runoff forecasting increase significantly [4][5][6]. Researchers and engineers are devoted to establishing a more comprehensive forecasting model that facilitates the scientific management of limited water resources under a changing environment [7][8][9]. However, the natural runoff process exhibits strong nonlinearity and nonstationarity owing to the combined influences of multiple factors, such as meteorological events, natural geography and watershed features [10][11][12][13]. Thus, accurate runoff forecasting remains an important but challenging research topic for hydrology experts and scientists.
In the last few years, various runoff forecasting methods have been successfully established, which can be broadly categorized into two different groups [14][15][16]: process-based approaches and data-driven approaches. Process-based approaches often involve intricate models, thorough knowledge of the physical mechanisms underlying runoff processes, sufficient hydrometeorological data, and scientific expertise judgment. The strict requirements pose application limitations that may lead to poor prediction performance and uncertainty. To address this problem, data-driven approaches have become the primary approach to produce reliable forecasting information for reservoir operation and hydropower energy management. Mathematical statistics techniques have been widely employed by many researchers to forecast nonstationary hydrological data series [16][17][18]. Although the linear relationship between predictors and predicted values can be well identified, these models cannot provide satisfying prediction results as the highly nonlinear characteristics inherent in runoff series are not given full consideration. With the advancement of intelligent computing, numerous artificial intelligence techniques have been widely used in runoff prediction [19,20]. Compared with the conventional regression approaches, the machine learning methods have demonstrated significant improvements in prediction accuracy.
As a classical machine learning approach, Gaussian process regression (GPR) relies on Bayesian theory and statistical learning theory [21][22][23]. By replacing the basis function used in Bayesian linear identification, GPR can address the complicated regression problems with small sample sizes and high dimensionality [24][25][26]. Compared to the traditional forecasting models, GPR has the advantages of easy implementation, adaptively acquiring hyperparameters, and producing probabilistic outputs. GPR has gained significant attention in regression problems, such as runoff, wind power, and solar power forecasting. Generally, the conjugate gradient method is often used to obtain GPR hyperparameters. However, it suffers from high dependence on the initial values and difficulty in determining the number of iterations. In other words, the GPR model may suffer from local convergence and strong parameter dependence, which obviously limit its practicality and interpretability in runoff prediction. To mitigate the research gaps, it is essential to find more practical methods to improve GPR's performance for forecasting nonlinear runoff series [27][28][29].
Recently, a novel meta-heuristic cooperation search algorithm (CSA) was developed to resolve intricate engineering optimization problems [30]. In the foundational concept of the CSA method, each solution can be viewed as a staff member in a teamwork setting and multiple solutions form the swarm for evolutionary computing. The swarm converges gradually towards the promising search regions around the optimal solution with three operators, including the team communication operator for improving global search ability, the reflective learning operator for facilitating local search ability, and the internal competition operator for ensuring the survival of elite solutions. The CSA approach has been used to resolve numerical optimization problems. Despite its potential, there is little research on applying the CSA to promote the GPR performance in runoff forecasting. Thus, the paper proposes a hybrid GPR-CSA model that leverages the CSA algorithm to enhance the generalization ability of GPR for runoff forecasting. The experiments show that compared with conventional models, the GPR-CSA approach offers better prediction accuracy for different applications. In summary, this paper contributes multiple effective models for forecasting the real-world runoff data series; moreover, a novel GPR-CSA method with better forecasting accuracy can be applied for nonlinear regression tasks, including runoff prediction, in various scenarios. The application of the GPR-CSA to three hydrological stations in China demonstrates that the GPR-CSA is able to fully identify the high-dimensional relationship between predictors and predicted values, providing an effective artificial intelligence model for addressing hydrologic forecasting problems.
The layout structure of this article is given as follows: Section 2 outlines the specifics of our methodology for runoff forecasting; Section 3 gives four evaluation criteria; Section 4 examines the applicability of the proposed method; and Section 5 gives the conclusion.

Gaussian Process Regression (GPR)
Suppose that y and t are the dependent and independent variable belonging to the real number set R. The regression problem that further incorporates noisy information can be expressed as: Here, ε ∈ N 0, σ 2 and f (·) are the measurement error and an unidentified function. Based on the Gaussian process technique, f (·) is a stochastic function with both a mean function u(·) and a covariance function k(·, ·). The relationship between different instances can be formulated as below: where θ stands for the hyper-parameters requiring estimation.
For an observation set D = {(t 1 , y 1 ), (t 2 , y 2 ), · · · , (t n , y n )}, the equation can be expressed as the following equation: where {ε i } i=1,2,··· ,n are the random disturbances in the normal distribution with mean 0 and variance σ 2 . Thus, the joint distribution of y 1 , y 2 , · · · , y n obeys a multivariate normal distribution, which can be expressed as follows: where u i = u(t i ) is entries of average matrix µ. Ψ denotes an n × n matrix whose (i, j)th element is specified by where δ ij denotes the Kronecker delta. It is assumed that t * and y * are the testing point and possible response value, and for the training set D, the conditional distribution is characterized by a normal distribution where the average and variance are computed as follows: where ψ(t * ) = [k(t * , t 1 ; θ), · · · , k(t * , t n ; θ)] T is the covariance between f (t * ) and f = [ f (t 1 ), · · · , f (t n )] T . Ψ is the covariance matrix of (y 1 , y 2 , . . . , y n ) T . In GPR, the parameters are the hyperparameters θ in the covariance function, the noise variance σ 2 , and other coefficients (indicated generally by β) in the mean function u(·). These parameters are obtained by finding the maximum succeeding marginal log-likelihood function [31][32][33]:

Cooperation Search Algorithm (CSA)
The cooperation search algorithm (CSA) is a novel and effective meta-heuristic tool for tackling complex optimization problems [34][35][36]. In the CSA optimizer, the target problem is regarded as a growing company. Each solution is perceived as an employee, while a group of solutions forms a team. The supervisor committee consists of the personal best-known staff members, while the executive committee is composed of an elite staff member set that contains M global best-known solutions. The chairman-in-office is picked from the executive committee in a random manner. The CSA methodology starts with the random initialization of solutions, and then all solutions are dynamically updated to gradually discover high-quality solutions through three operators. Specifically, the team communication operator determines the solution's probability of being influenced by leader solutions, whereas the reflective learning operator determines whether to learn from its own best-known position or that of its supervisor. The internal competition operator selects solutions with better fitness values to compete for leadership positions, promoting distribution diversity and global search ability of the swarm. By the above procedures, the CSA method has been widely applied to resolve the complicated optimization problems in different engineering fields. Figure 1 shows the schematic diagram of the CSA method. Then, the technical details that are crucial for solving the multivariable optimization problems are as below.

Proposed Runoff Forecasting Method
The GPR method has satisfying performance in learning efficiency and satisfactory generalization ability compared to the traditional forecasting models. However, the standard GPR method may yield undesirable hydrologic forecasting outcomes in actual application scenarios because unsatisfying parameter combinations will lead to the local convergence problem. To efficiently overcome this problem, this paper establishes a hybrid evolutionary artificial intelligence model where the first-rank parameter combinations of the GPR model is determined by the CSA method, and then the optimized GPR model to predict runoff data series in the coming periods. By linking the advantages of both GPR and CSA, this study can offer a more robust runoff forecasting model with higher compactness than the traditional GPR model. As showed in Figure 2, the specific process of the proposed model is given as follows:  (1) Team building phase. The initial positions of all staff members in the team are determined by Equation (8). Based on the fitness value, M elite solutions are used to form the exterior leader set.
where I is the number of solutions. The jth value of the ith solution at the kth iteration is represented by x k i,j . A uniformly distributed random number in [L, U] is denoted with φ(L, U). x j and x j are the lower and upper bound of the jth variable, respectively.
(2) Team communication operator. Each solution has the opportunity to acquire fresh insights from leader staff members. As showed in Equation (9), the team communication operator uses three components: the expertise A from the chairman, the cumulative knowledge B from the leader staff members in board of directors, and the combined knowledge C from leader staff members in the board of supervisors. The chairman is selected randomly from M global best-known solutions, whereas all directors and supervisors are assigned equal roles. The detailed equations are given as below: where the jth element of the ith group agent at iteration k + 1 is denoted with u k+1 i,j . The jth element of the ith agent's best-known position at iteration k is represented by pBest k i,j . The jth element of the indth global best-known agent is represented by gBest k ind,j . ind is the integer randomly chosen from {1,2, . . . ,M}. The expertise affected by the chairman is denoted with A k i,j . The mean expertise affected by M global best-known staff members and I personal best-known staff members are denoted as B k i,j and C k i,j . α and β are the learning parameters.
(3) Reflective learning operator. In addition to studying from elite agents, each agent can also acquire new information by considering their own experiences and observations, which can be represented as follows: where the jth element of the ith reflective agent at iteration k + 1 is denoted as v k+1 i,j .
(4) Internal competition operator. By guaranteeing retention of the high-quality agents, the competitiveness of the swarm can be gradually enhanced by the following equation: where F(x) is the specific fitness score associated with staff x. Based on the fitness value, the board of directors and the board of supervisors are updated by the following equation: where pBest k+1 i is the ith personal best-known solution at iteration k + 1. gBest k+1 is the board of directors containing M global best-known solutions.

Proposed Runoff Forecasting Method
The GPR method has satisfying performance in learning efficiency and satisfactory generalization ability compared to the traditional forecasting models. However, the standard GPR method may yield undesirable hydrologic forecasting outcomes in actual application scenarios because unsatisfying parameter combinations will lead to the local convergence problem. To efficiently overcome this problem, this paper establishes a hybrid evolutionary artificial intelligence model where the first-rank parameter combinations of the GPR model is determined by the CSA method, and then the optimized GPR model to predict runoff Water 2023, 15, 2111 6 of 17 data series in the coming periods. By linking the advantages of both GPR and CSA, this study can offer a more robust runoff forecasting model with higher compactness than the traditional GPR model. As showed in Figure 2, the specific process of the proposed model is given as follows:

Proposed Runoff Forecasting Method
The GPR method has satisfying performance in learning efficiency and satisfactory generalization ability compared to the traditional forecasting models. However, the standard GPR method may yield undesirable hydrologic forecasting outcomes in actual application scenarios because unsatisfying parameter combinations will lead to the local convergence problem. To efficiently overcome this problem, this paper establishes a hybrid evolutionary artificial intelligence model where the first-rank parameter combinations of the GPR model is determined by the CSA method, and then the optimized GPR model to predict runoff data series in the coming periods. By linking the advantages of both GPR and CSA, this study can offer a more robust runoff forecasting model with higher compactness than the traditional GPR model. As showed in Figure 2, the specific process of the proposed model is given as follows:  Step 1: Preparatory work. The computation parameters of the proposed GPR-CSA method are set before calculation, such as the maximum number of iterations K, solutions I, leaders M in the CSA method, and the kernel function in the GPR model.
Step 2: Parameter optimization. Based on the training data, the detailed procedures of the CSA method to determine the GPR parameters are given as follows: Step 2.1: Define the counter k = 1. Then, use Equation (8) in the team building phase to create the initial population in the feasible zone.
Step 2.2: Evaluate the fitness value of all staff members to update the optimal position of each staff member and the globally best-known position of the swarm.
Step 2.3: Use the team communication operators defined in Equations (9)-(12) to enhance global exploration as while the reflective learning operators in Equations (13)- (16) to improve local exploitation. Then, the internal competition operator in Equation (17) is used to select better solutions for the iteration k + 1.
Step 2.4: Increment the counter k by 1. If k is smaller than the maximum iteration, go to Step 2.2; otherwise, the globally best-known position of the CSA represents the ideal GPR model parameters.
Step 3: Operational prediction. Utilize the optimized GPR model to forecast the potential predicted values of the new predictors in the testing dataset.

Performance Evaluation Criteria
In this section, four evaluation criteria are used to evaluate the performance of the developed models in hydrologic forecasting, including root mean squared error (RMSE), mean absolute error (MAE), correlation coefficient (R), and Nash-Sutcliffe efficiency (NSE). In practical applications, a reliable and robust model can be capable of producing lower where y i and y i are the ith point in the recorded and estimated dataset. The average of all the recorded and estimated points are denoted with y avg and y avg . The set of the data being evaluated is represented by n.

Engineering Background
The presented GPR-CSA model is used to predict the nonlinear and nonstationary runoff data series of three representative hydrological stations in China, i.e., SX station with daily runoff data, LYX station and TNH station with weekly runoff data. Next, this study conducts a comparative analysis utilizing the runoff data of the SX, LYX and TNH stations. Figure 3 offers a comprehensive overview and statistical information of the recorded runoff data. The recorded runoff data are partitioned into two distinct subsets, with the first 70% reserved for training and validation, while the last 30% are reserved for testing.

Engineering Background
The presented GPR-CSA model is used to predict the nonlinear and nonstationary runoff data series of three representative hydrological stations in China, i.e., SX station with daily runoff data, LYX station and TNH station with weekly runoff data. Next, this study conducts a comparative analysis utilizing the runoff data of the SX, LYX and TNH stations. Figure 3 offers a comprehensive overview and statistical information of the recorded runoff data. The recorded runoff data are partitioned into two distinct subsets, with the first 70% reserved for training and validation, while the last 30% are reserved for testing.

Model Development
Several runoff forecasting models were applied to check the effectiveness of the proposed model, including linear regression (LR), artificial neural network (ANN), recurrent neural network (RNN), Gaussian process regression (GPR), and long short-term memory network (LSTM). It is worth noting that appropriate inputs play an important place in improving the performance of machine learning models. Thus, expert knowledge and

Model Development
Several runoff forecasting models were applied to check the effectiveness of the proposed model, including linear regression (LR), artificial neural network (ANN), recurrent neural network (RNN), Gaussian process regression (GPR), and long short-term memory network (LSTM). It is worth noting that appropriate inputs play an important place in improving the performance of machine learning models. Thus, expert knowledge and partial autocorrelation functions were implemented to identify suitable inputs. For the SX, LYX and TNH stations, the input variables are all three antecedent runoffs at period t − 1, t − 2, and t − 3, which are used to forecast the runoff at period t + τ, where τ is the forecasting period. Moreover, the parameters for the models were determined as follows: for the ANN, RNN and LSTM models, the activation function was set as the sigmoid function, and the Adam optimizer was used for tuning parameters; for GPR, the combined RBF and rational quadratic kernel was used while the default parameter configurations in the scikit-learn toolbox in Python 3.10 were adopted for tunning parameters.

Case 1: One-Step-Ahead Prediction Outcomes
This study develops several runoff forecasting models and then compares their performances at different forecasting horizons by four statistical evaluation indicators, namely RMSE, MAE, R, and NSE. Table 1 lists the detailed statistical indicators of one-step-ahead predicting results by the GPR-CSA method and several control methods. The data illustrate that the data at the SX station present the largest forecasting biases at both the training and testing sets, and the TNH station and LYX station present the lower forecasting biases than the SX station. This phenomenon lies in the higher runoff values of the SX station, showing the significant impacts of the runoff dataset at hand. Table 1 also supports the following interesting conclusions: (a) the GPR-CSA model displays superior forecasting ability owing to its lowest values of RMSE and MAE, coupled with the highest R and NSE indicators; (b) the standard GPR method outperforms the LR model in terms of the fitting ability and overall performance, highlighting the importance of the employed model structure; (c) the comparison of the GPR and GPR-CSA models highlights the validity of the CSA algorithm in identifying feasible parameters, demonstrating the superior ability of CSA in optimizing multivariable combinations. To illustrate the prediction ability of the developed model, Figure 4 shows the onestep-ahead predicting results for the analyzed runoff data series. It shows that the runoff predicting curves of the GPR-CSA model approximate the original runoff curves, better than the other models at the three hydrological stations. Figure 5 depicts the radar plots of R and NSE at the one-step-ahead forecasting results. It can be seen that the values of R and NSE of the GPR-CSA method in the radar plots are the farthest to the center zone, proving the effectiveness of GPR-CSA in providing satisfactory predicting results. Thus, the GPR-CSA model can trace the complex features of runoff data, thereby leading to highly satisfactory forecasting outcomes. These findings serve to prove the feasibility of the model for hydrological forecasting.

Case 2: Two-Step-Ahead Prediction Outcomes
The refined model boasts adequate forecasting accuracy in the above one-step-ahead runoff prediction. In real-world scenarios, the forecasting model's performance at various horizons is also critical to promote water resource utilization. Consequently, the two-stepahead runoff predicting results are compared. As outlined earlier, Table 2 gives the statistical indicators of the predicting outcomes by various models. For both the training and   The refined model boasts adequate forecasting accuracy in the above one-step-ahead runoff prediction. In real-world scenarios, the forecasting model's performance at various horizons is also critical to promote water resource utilization. Consequently, the two-stepahead runoff predicting results are compared. As outlined earlier, Table 2 gives the statis-

Case 2: Two-Step-Ahead Prediction Outcomes
The refined model boasts adequate forecasting accuracy in the above one-step-ahead runoff prediction. In real-world scenarios, the forecasting model's performance at various horizons is also critical to promote water resource utilization. Consequently, the twostep-ahead runoff predicting results are compared. As outlined earlier, Table 2 gives the statistical indicators of the predicting outcomes by various models. For both the training and testing data, the statistical data fully highlight the superiority of the GPR-CSA model compared to other control methods. Thus, this section provides further evidence of the engineering feasibility of the hydrological forecasting approach.  Figure 6 illustrates the scatter plots of two-step-ahead predicting results for the testing dataset derived through several techniques. It shows that the proposed model exhibits superior prediction accuracy in comparison to other models, as it attains the largest correlation between the recorded and predicted runoff in all simulations. Figure 7 shows the bar graphs of the RMSE and MAE for the two-step-ahead predicting results at the testing phase. It shows that GPR-CSA has smaller RMSE and MAE values compared with other forecasting methods, demonstrating the superiority of CSA in forecasting nonstationary runoff series. Thus, incorporating artificial intelligence and metaheuristic optimization can effectively meet the practical needs of hydrological forecasting tasks.  Table 3 displays the statistical metrics of three-step-ahead predicting results obtained through multiple models. Figure 8 illustrates the statistical indicators of the three-stepahead estimated results for the testing dataset. The data show that the standalone models    Table 3 displays the statistical metrics of three-step-ahead predicting results obtained through multiple models. Figure 8 illustrates the statistical indicators of the three-stepahead estimated results for the testing dataset. The data show that the standalone models    Table 3 displays the statistical metrics of three-step-ahead predicting results obtained through multiple models. Figure 8 illustrates the statistical indicators of the three-stepahead estimated results for the testing dataset. The data show that the standalone models yield limited forecasting results, while the evolutionary algorithm substantially improves the achieved outcomes in different cases. Compared to other models, the developed GPR-CSA model attains the best forecasting performance for the testing datasets at the three stations. Thus, the introduced parameter optimization strategy can significantly enhance the forecasting effectiveness of a standalone model for runoff forecasting. yield limited forecasting results, while the evolutionary algorithm substantially improves the achieved outcomes in different cases. Compared to other models, the developed GPR-CSA model attains the best forecasting performance for the testing datasets at the three stations. Thus, the introduced parameter optimization strategy can significantly enhance the forecasting effectiveness of a standalone model for runoff forecasting.

Analysis of the Kernel Function
The experiments are executed to show the influences of different kernel functions on the predicting results at the three stations. Table 4 gives the statistical indicators of onestep-ahead predicting results using different kernel functions, where kernel1 represents the radial basis function kernel, kernel2 represents the rational quadratic kernel, and kernel3 represents the compound kernel with the radial basis function kernel and rational quadratic kernel. Figure 9 draws the correlation coefficients of the predicting results with different kernels at the SX station. The following phenomenon can be observed: (1) compared with the standard GPR model, GPR-CSA achieves better prediction results regardless of the employed kernel functions; (2) for the same station, the predicting results change with the kernel functions, and thus, it is necessary to carefully select the kernel function according to the actual runoff situation by experiments; (3) with the extension of the forecasting period, the prediction performances of three kernel functions decrease gradually at the SX station. In all cases, the GPR-CSA model is always superior to the GPR model, demonstrating the superiority of the CSA method in finding suitable computation parameters of the GPR model. Thus, the GPR-CSA method is an effective tool to provide accurate hydrological forecasting information.

Analysis of the Model Robustness
To analyze the robustness of the GPR-CSA model, the GPR-CSA model was run 50 times for the one-step-ahead prediction at the three stations. Figure 11 shows the RMSE of the one-step-ahead predicting results from the GPR-CSA and ANN model at the testing phase in different runs. Compared with the conventional ANN method, the predicting results of the GPR-CSA model at the testing dataset are stable with little fluctuations, which demonstrates the reliable forecasting performance of the GPR-CSA method. Thus, the GPR-CSA model has an outstanding parameter optimization ability and a robust performance for runoff prediction.

Analysis of the Model Robustness
To analyze the robustness of the GPR-CSA model, the GPR-CSA model was run 50 times for the one-step-ahead prediction at the three stations. Figure 11 shows the RMSE of the one-step-ahead predicting results from the GPR-CSA and ANN model at the testing phase in different runs. Compared with the conventional ANN method, the predicting results of the GPR-CSA model at the testing dataset are stable with little fluctuations, which demonstrates the reliable forecasting performance of the GPR-CSA method. Thus, the GPR-CSA model has an outstanding parameter optimization ability and a robust performance for runoff prediction.

Analysis of the Model Robustness
To analyze the robustness of the GPR-CSA model, the GPR-CSA model was run 50 times for the one-step-ahead prediction at the three stations. Figure 11 shows the RMSE of the one-step-ahead predicting results from the GPR-CSA and ANN model at the testing phase in different runs. Compared with the conventional ANN method, the predicting results of the GPR-CSA model at the testing dataset are stable with little fluctuations, which demonstrates the reliable forecasting performance of the GPR-CSA method. Thus, the GPR-CSA model has an outstanding parameter optimization ability and a robust performance for runoff prediction.
times for the one-step-ahead prediction at the three stations. Figure 11 shows the RMSE of the one-step-ahead predicting results from the GPR-CSA and ANN model at the testing phase in different runs. Compared with the conventional ANN method, the predicting results of the GPR-CSA model at the testing dataset are stable with little fluctuations, which demonstrates the reliable forecasting performance of the GPR-CSA method. Thus, the GPR-CSA model has an outstanding parameter optimization ability and a robust performance for runoff prediction.

Conclusions
Accurate hydrological prediction is critical for the effective management of water energy resources. For addressing practical demands, this article proposes a hybrid artificial intelligence model for predicting runoff under uncertainty. For the first time, the cooperation search algorithm (CSA) is used to find the suitable parameter combinations of the classical Gaussian process regression (GPR) model. Through three well-designed operators, the CSA tool effectively overcomes the local convergence defects associated with the traditional gradient-based methods. To validate its efficacy, the GPR-CSA model is used to predict the nonlinear runoff data of three hydrological stations. The simulations indicate that the CSA method achieves a balance between global search and local search by optimizing the computational parameters of the traditional GPR model. Moreover, the results of the GPR-CSA method are better than several control models in both training and testing datasets. Thus, a novel hybrid artificial intelligence model is proposed for forecasting nonstationary and nonlinear streamflow. Combining the advantages of the GPR model and the CSA algorithm, the proposed model shows superior forecasting ability and robust prediction results in addressing complex hydrologic forecasting tasks. The research presented in this paper has resulted in innovative outcomes for the application of artificial intelligence methods in the field of hydrological forecasting. The findings have facilitated early warnings of flood disasters in river basins and the efficient utilization of water resources, making them a significant contribution to this field.

Data Availability Statement:
The data supporting the findings of this paper are available from the corresponding author upon reasonable request.