Next Article in Journal
Evaluation of GPM and TRMM and Their Capabilities for Capturing Solid and Light Precipitations in the Headwater Basin of the Heihe River
Previous Article in Journal
Study on Performance and Operation Mechanism of a Separation Equipment for a PWR Steam Generator
Previous Article in Special Issue
General Method Based on Regressive Relationships to Parameterize the Three-Parameter Depth–Duration–Frequency Curve
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Intercomparison and Assessment of Stand-Alone and Wavelet-Coupled Machine Learning Models for Simulating Rainfall-Runoff Process in Four Basins of Pothohar Region, Pakistan

1
Department of Agricultural Engineering, Faculty of Agricultural Sciences and Technology, Bahauddin Zakariya University, Multan 60000, Pakistan
2
School of Engineering, University of Basilicata, 85100 Potenza, Italy
3
NUST Institute of Civil Engineering, National University of Sciences and Technology, Islamabad 44000, Pakistan
4
Department of Agricultural Engineering, MNSU Agriculture, Multan 60000, Pakistan
5
Department of Civil Engineering, The University of Lahore, Lahore 54590, Pakistan
*
Authors to whom correspondence should be addressed.
Atmosphere 2023, 14(3), 452; https://doi.org/10.3390/atmos14030452
Submission received: 14 December 2022 / Revised: 11 February 2023 / Accepted: 21 February 2023 / Published: 24 February 2023
(This article belongs to the Special Issue Climate Change Impacts on Urban Stormwater Management)

Abstract

:
The science of hydrological modeling has continuously evolved under the influence of rapid advancements in software and hardware technologies. Starting from simple rational formulae for estimating peak discharge and developing into sophisticated univariate predictive models, accurate conversion of rainfall into runoff and the assessment of inherent uncertainty has been a prime focus for researchers. Therefore, alternative data-driven methods have gained widespread attention in hydrology. Moreover, scientists often couple conventional machine learning models with data pre-processing techniques, i.e., wavelet transformation (WT), to enhance modelling accuracy. In this context, this research work attempts to explore the latent linkage between rainfall and runoff in Pothohar region of Pakistan by developing a novel linkage of five streamline techniques of machine learning, including single decision tree (SDT), decision tree forest (DTF), tree boost (TB), multilayer perceptron (MLP), and gene expression modeling (GEP), with a more sophisticated variant of WT, i.e., maximal overlap discrete wavelet transformation (MODWT), for boundary correction of the transformed components of timeseries data. This study also implements these machine learning models in a stand-alone mode for a more comprehensive comparative analysis of performances. Furthermore, the study uses a combined-basin approach that divides Pothohar region into two basins to compensate for the complex topographic division of the study area. The results indicate that MODWT-based DTF outperformed other stand-alone and hybrid models in terms of modeling accuracy. In the first scenario, considering the Bunha-Kahan River basin, MODWT-DTF yielded the highest NSE (0.86) and the lowest RMSE (220.45 mm) and R2 (0.92 at lag order 3 (Lo3)) when transformed with daubechies4 (db4) at level three. While in the Soan-Haro River basin, MODWT-DTF produced the highest accuracy modeling at lag order 4 (Lo4) (NSE = 0.88, RMSE = 21.72 m3/s, and R2 = 0.91). The highly accurate performance of 3- and 4-days lagged models reflects the temporal consistency in hydrological response of the study area. The comparison of simple and hybrid model performance indicates up to a 55% increase in modeling accuracy due to data pre-processing with wavelet transformation.

1. Introduction

Rainfall-Runoff (RR) modeling is one of the classical applications of hydrology. It has the purpose of simulating a river flow hydrograph in a given cross river section induced by an observed or a hypothetical rainfall forcing [1]. Depending on their complexity, rainfall-runoff models can also simulate the dynamics of water quality, ecosystems, and other dynamical systems related to water, therefore embedding laws of chemistry, ecology, social sciences, and other fields into the model. These models can be used for flood simulation, forecasting, and prevention [2,3,4,5,6], management of water resources [7,8], and water supply simulation and forecasting [9,10]. Moreover, the applications also include land use and land cover (LULC) impact assessment [11] and management practices and water quality evaluation [12].
The hydrology and water resources domains have witnessed growing interest in using data-driven models to improve the simulation of RR processes whether by directly replacing or being used in conjunction with classical hydrological models (HMs) [13].
The physics-based modeling methods, for example the Soil and Water Assessment Tool (SWAT), solve hydrological problems through computing the conservation equations of momentum, mass, and energy and mathematically quantifying real-world physics [14]. Meanwhile, data-driven modeling approaches, for example machine learning, depend primarily on numerical data to solve hydrological problems. The last decade has seen rapid developments and a trend shift toward data-driven hydrological modeling, primarily in order to overcome the uncertainty and high parametric requirement of physically-based models [15].
Among the many different methods of soft computing, Artificial Neural Networks (ANN) have been shown to be the method of preference when it comes to modeling the complicated rainfall-runoff phenomena. Neural networks were used by Tokar and Johnson [16] to estimate daily runoff as a product of daily precipitation, temperature, and snowmelt. According to the findings of the research, the ANN model had a higher degree of accurate prediction when compared with both regression and conceptual models. Wilby et al. [17] developed a conceptual and neural network rainfall-runoff model using the data they obtained from precipitation, evaporation, and discharge. To assess the degree to which neural networks can reliably ingest the hydrological processes, three separate experiments were carried out, each with a gradually diminishing quantity of information. According to the findings of the research, a neural network consisting of seven inputs and three hidden nodes could comprehend the behavior of the conceptual model.
Most recently, there has been a significant shift toward the use of hybrid soft computing strategies for the purpose of resolving difficulties faced in real life. The symbiotic relationship that exists between the various soft computing methodologies is seen to be the driving force for this expansion. The combination of methods helps to overcome the shortcomings of each of the separate approaches and ultimately results in the creation of reliable computational strategies [18]. For example, Okkan et al. [19] developed a novel nested hybrid rainfall-runoff modelling framework based on the confluence of machine learning (ANN and support vector machine (SVM)) and a conceptual rainfall-runoff model (dynamic water balance model (dynwbm)). The findings reflect that the nested hybrid model bested the standalone machine learning, and the conceptual models, and the coupled models. Likewise, Poonia and Tiwari [20] used two variants of ANN, radial basis function (RBF) and feed-forward back propagation (FFBP), to link rainfall events with runoff generation in the Hoshangabad basin of the Narmada River. The outcomes supported RBF as a highly accurate rainfall-runoff modeling technique (i.e., R2 = 0.9964) in the basin. Similarly, Gomes and Blanco [21] coupled conventional ANN models with MODWT to estimate daily rainfall in the Tocantins-Araguaia Hydrographic Region, Brazil. The results indicate an enhancement in the performance of hybrid MODWT-ANN models (i.e., NSE = 0.81–0.95) when compared to stand-alone ANN models. A similar study by Nourani et al. [22] formulated a wavelet-based M5 tree model (Wavelet-M5 model) and employed the model in a rainfall-runoff simulation of the Sardrud catchment in Iran. The authors also compared the novel hybrid model with a simple M5 tree model and a hybrid wavelet-ANN (WANN) model. The study revealed that wavelet coupling enhanced modeling performance of the M5 model by approximately 31%. In another study, Nourani et al. [23] systematically reviewed the application of wavelet-based hybrid Artificial Intelligence (AI) modelling in hydro-climatology. The authors primarily explored the wavelet-hybrid applications of AI models, including ANN, decision tree, and random forest models, in different contexts such as sediment modelling, flow forecasting, precipitation modelling, and rainfall-runoff modelling. Ouma et al. [24] compared the performance of a hybrid machine learning model (WNN) with a deep learning model long short-term memory (LSTM) in rainfall and runoff time-series trend analysis in the Nzoia basin. The results reflected the higher prediction accuracy of LSTM (i.e., R2 = 0.861) when compared to WNN models (0.782). Saifullah et al. [25] explored the rainfall-runoff linkage of the Jhelum River in Pakistan by using gene expression programming (GEP), a support vector machine coupled with radial biases function (RBF-SVM), and an M5-Model tree. The authors concluded that GEP performed better in light of observations made from flow duration curves. Similarly, Kavoosi and Khozeymehnnhad [26] also employed GEP to simulate rainfall-runoff process in the Halil River in Iran. The authors compared the performance with neural network and fuzzy logic models. The outcomes indicated that GEP was outperformed by a particle swarm optimized-adaptive neural fuzzy inference system (ANFIZ-PSO). Asadi et al. [27] used a genetic algorithm (GA) for the purpose of evolving the weights of the neural network that was utilized for modeling the rainfall-runoff process. To improve the accuracy of the model’s predictions, the data were pre-processed using several techniques, including data transformation, input variable selection, and data clustering. According to the findings of the research, using this technique results in a more expedient learning process, a high degree of accuracy, and good adaptability to the nonlinear functional connections that exist between rainfall and runoff. Similarly, Okkan and Serbes [28] employed three machine learning techniques–feed forward neural networks (FFNN), multiple linear regression (MLP), and least square support vector machine (LSSVM)–with discrete wavelet transformation (DWT) to model reservoir inflow of the Demirkopru Dam in Turkey. The authors observed that the DWT-FFNN hybrid model outperformed the other models used in the study. In another study, Okkan [29] proposed wavelet-based WNN and wavelet-multi linear regression (WREG) hybrid models to predict the monthly reservoir inflow of the Kemer Dam in Turkey. The outcomes indicated that integration of DWT with FFNN and REG models enhanced the prediction accuracy of hybrid models when compared to the standalone models.
Researchers have also used more advanced and sophisticated structures of artificial intelligence in hydrology and climatology problems. For instance, a deep learning (DL) model, the convolutional neural network (CNN), has been implemented to increase precipitation forecasting accuracy [30], rainfall-runoff modelling [31], and non-periodic flow prediction [32]. Similarly, basic and improved structures of long short-term memory (LSTM) models have been employed in solving hydro-climatological problems including precipitation forecasting [33], rainfall-runoff modelling [34], and wave height prediction [35]. The above cited works emphasize the importance of rainfall-runoff modeling and the application of machine learning techniques in developing high accuracy models for understanding this latent hydrological relationship. In this context, in a previous study [18], the authors applied five machine learning techniques and wavelet pre-processing to understand rainfall-runoff dynamics of the Soan River basin. However, the RR process of the other five basins, which play an equally important role in the hydrology of the Pothohar region, has not been analyzed. Indeed, the literature review highlights a significant research gap in terms of an unclear quantification of rainfall-runoff phenomenon in this region. Moreover, a comprehensive study on performance comparison of relevant machine learning techniques used for RR modeling on diverse and/or interrelated river basins is also required to demonstrate the applications of machine learning in hydrology. Lastly, the influence of using advanced wavelet transformation-based pre-processed input data on diverse machine learning techniques performance is still unclear.
This research work applies five of the major and most used machine learning techniques as stand-alone models and as hybrid models coupled with maximal overlap discrete wavelet transformation (MODWT) to evaluate the upshot of input data transformation on the modeling accuracy. The machine learning techniques can consist of logistic and linear regression models, naïve bayes, support vector machine, k-Nearest neighbor, and gradient boosting models but we have selected single decision tree (SDT), decision tree forest (DTF), tree boost (TB), multilayer perceptron (MLP), and gene expression modeling (GEP) owing to their contemporary modeling variants and inordinate and incessant applications in environmental and hydrological contexts. The above-mentioned machine learning methods have been discreetly explored in hydrological settings in previous studies [36,37,38,39,40]. However, this study will provide thorough research work on the application of machine learning in rainfall-runoff modeling along with employment of wavelet pre-processing for enhancing modeling accuracy. Moreover, the present study aims to contribute to the existing research gaps by studying the linkage between rainfall and runoff in the Pothohar region of Pakistan. The multi-basin study covering ML techniques and wavelet transformation in the Pothohar region is the original effort by the authors.
Hence, the main objectives of this study are (i) the clear determination of linkages between rainfall events and runoff generation in the Pothohar region of Pakistan, (ii) the stand-alone application and assessment of the selected ML approaches (DTF, SDT, MLP, TB, and GEP), (iii) coupling each of the above-mentioned technique with MODWT to quantify the influence of input pre-processing on model results, and (iv) the comparison of the performance of each machine learning technique in the basins of the Pothohar region.

2. Materials and Methods

2.1. The Study Area and Datasets

Covering a total area of nearly 22,500 km2 and housing a population of around 17,460,000, the Pothohar plateau is one of the hydrologically significant regions of Pakistan. Located at the upper boundary of Punjab province, this region spans from 32.5–34° N to 72–74° E. The region is surrounded by two medium-altitude mountain/hill ranges to the north (Margalla hills) and south (salt range) and is bounded by the Indus River in the west and the Jhelum River in the east. Soan and Haro are the two major rivers (or tributaries) in the region that meet the Indus River [41]. Moreover, there are a total of six different tributary rivers (Soan, Haro, and Reshi meeting the Indus River downstream of Tarbela dam, while Bunha, Kahan and Kanshi meet the Jhelum River near Mangla dam). As such, we have subdivided the Pothohar basin in to six sub-basins that are discussed below.
The study area could be characterized as a “Subtropical Triple Season Moderate Climate Zone” with a single rainy period from July to September that has a cooling effect on temperatures. Even though rain falls throughout the year, during this rainy period it rains particularly heavily and lasts from July through September. During the rainy season, roughly 600 mm of rain accumulates, contributing almost 60% of the yearly precipitation of around 1000 mm [42]. Furthermore, crop and vegetative production is highly unpredictable due to the late monsoon and inconsistent winter precipitation, both of which are prevalent occurrences. The geography of the hilly terrain with sharp elevation gradients aids the formation of multiple tributaries flowing at high speeds which degrade the arable soil. Due to steep slopes, precipitation water does not have time to infiltrate through the soil and form any ground groundwater reserve; this causes the rainfall to damage the soil, leading to land degradation. As a result, agricultural activity and production in this region is completely reliant on precipitation (rainfall), that can be scarce and non-homogenous at the time of irrigation. Drought is a common occurrence, as evidenced in recent years [43].
Figure 1 presents the location map of study area and its fragmentation into sub-basins. The locations of contributing rain gauge stations throughout the study area are also shown in Figure 1. The Pothohar region consists of multiple sub-basins: each sub-basin has a water level gauging site where surface runoff from all over the sub-basin enters the river. For the upper three sub-basins (Soan, Haro, and Reshi), runoff naturally flows towards the Indus River, entering the river at three different locations downstream of the Tarbela reservoir. Similarly, the three lower sub-basins (Bunha, Kahan, and Kanshi) meet the Jhelum River at three different locations downstream of the Mangla reservoir. Since runoff from the Pothohar region enters two different rivers at a minimum of six different gauged locations, it is not feasible to analyze rainfall and runoff data using single basin analysis. Moreover, the temporal longevity of daily data added to the spatial enormity of the study area, consisting of more than 20 rainfall stations, would have formed an enormous and hard to handle dataset; this would have also magnified the processing and data handling errors. Also, the processing time would have increased by orders of magnitude. In this light, two homogeneous areas were identified in the study: one consists of the two main basins (Soan and Haro) meeting the Indus River and the other consists of the two main basins (Bunha and Kahan) meeting the Jhelum River. The Reshi and Kanshi River basins were excluded due to inconsistent and/or unavailable rainfall data. Moreover, the two gauging sites located in the lower basins of each part, as shown in Figure 1, are considered in the runoff data for that part. This is because gauge reading at downstream stations represents the flow accumulation at the water level gauging stations of the upstream basin.

Data Acquisition

The daily time series data of rainfall and runoff spanning 17 years, observed from 1999 to 2016, were acquired from three government agencies: Water and Power Development Authority (WAPDA), Pakistan Meteorological Department (PMD), and Soil and Water Conservation Research Institute (SAWCRI). Data from a total of 15 rain gauge stations were accumulated to determine daily rainfall depth in the Soan-Haro River basin, while data from 7 rain gauge stations were used to formulate model inputs in the case of the Bunha-Kahan River basin. The preference of rain gauge data over satellite data is due to homogeneous spatial distribution, high temporal resolution data, and the intrinsic reliability of in situ rain gauge stations. Similarly, the daily runoff data (m3/s) observed at Makhad (an outlet of the Soan River basin) were acquired from WAPDA. Sixteen years of daily gauge reading data (converted from meters to millimeters for unit and error consistency of the analysis) gathered downstream of the Bunha-Jhelum junction near Darapur were used due to unavailability of the discharge data and rating curve at this gauging station. Recent observations (2017–2021) were excluded from the analysis due to various inconsistencies in data at multiple stations and unavailability at some stations. Including incomplete data would have affected the reliability of outcomes. However, no significant change in data trends were observed from 2017 to 2021. Statistical characteristics of the datasets are summarized in Table 1 and Table 2 presents the summary of rain gauge and water level gauging station attributes.
Due to the unavailability of rating curve information at the cross sections where gauging stations are located, it was decided to use gauge reading (water level elevation available in meters) in place of runoff (m3/s) for the Bunha and Kahan rivers. To ensure the consistency of analysis, units of gauge reading were converted from meters to millimeters to increase the range of the data and model’s response to them.
Since gauges of the Bunha and Kahan River basins are present in the path of flow of the Jhelum River, the rainfall-runoff (gauge reading) relationship was modeled aggregating the upstream and downstream sub-basins in a unique river basin. Likewise, rainfall data from stations located in both the Bunha and Kahan River basins was used to cover the whole combined basin area. Water level gauge reading data from the Bunha River basin gauging station (near Bhambar) were used in the analysis because it is downstream of the Kahan River basin gauging station (near Naugran). As such, it also represents the contribution of the flow from the Kahan River.
In the same way, all the rain gauge stations included in the Soan and Haro River basin analysis were included. Runoff data from the Soan River basin gauging station (Makhad) were used as it is located downstream of Gurriala (Haro River basin gauging station) on the Indus River in the Tarbela dam command area and represents the flow from upper basin.

2.2. Machine Learning Techniques

2.2.1. Single Decision Tree (SDT)

A technique for representing data in a tree-like organization is Single Decision Tree. The SDT is also classified as a logical data model since it presents a multi-directional splitting mechanism. Nodes, or the rectangular boxes that make up the SDT’s data flow structure, are used. Every node in the specified dataset reflects the number of rows that are present. Runoff is chosen as the goal variable and rainfall is chosen as the predictor variable. Because there is only one predictor variable employed in the current investigation, rainfall is also used as the root node. Figure 2 shows the general structure of SDT; details of the technique can be found in Khan et al. [18].

2.2.2. Tree Boost (TB)

In Breiman’s [44] tree boost methodology, a prediction algorithm is improved by applying it to several repetitions of the boosting procedure. The weight of each function is then calculated by adding the results of all the functions. This stage improves forecast accuracy and eliminates any errors that may have occurred earlier in the process. This method’s better trait is the division of independent and dependent components. A general structure of TB is shown in Figure 3; further details can be found in Khan et al. [18].

2.2.3. Decision Tree Forest (DTF)

A specific category-related machine learning technique is the DTF. It consists of multiple trees used to combine projected assessments to offer a broad estimation of the supplied data. The Decision Tree Forest approach allows for the development of several trees concurrently and independently of one another. Breiman’s [45] original notion for the random forest (RF) approach serves as the foundation for the DTF technique. DTF general structure is presented in Figure 4. Khan et al. [18] have discussed DTF in detail.
Figure 2. General structure of Single Decision Tree (modified from Raza et al. [46]).
Figure 2. General structure of Single Decision Tree (modified from Raza et al. [46]).
Atmosphere 14 00452 g002
Figure 3. General structure of Tree Boost (modified from Raza et al. [46]).
Figure 3. General structure of Tree Boost (modified from Raza et al. [46]).
Atmosphere 14 00452 g003

2.2.4. Multilayer Perceptron (MLP)

MLP is a type of fully connected feed forward artificial neural network (ANN) developed by Webros [47]. It is composed of layers of neurons. An input layer, one or more hidden layers, and an output layer are the minimum number of levels that make up an MLP. One neuron makes up the output layer that displays the output of the MLP ANN; in this example, the output is the estimated runoff in m3/s and gauge reading of the Bunha-Kahan basin. The number of inputs in the dataset is the same as the number of neurons in the input layers [48]. Figure 5 shows the general structure of MLP. More details can be found in the previous study [18].
Figure 4. General structure of Decision Tree Forest (modified from Raza et al. [46]).
Figure 4. General structure of Decision Tree Forest (modified from Raza et al. [46]).
Atmosphere 14 00452 g004

2.2.5. Gene Expression Programming (GEP)

GEP is a method that uses adaptive computing, influences spontaneous development, and produces computer programs and quantitative models. It uses a dataset to provide an outcome in the form of a tree-like structure, as proposed in [49]. To address the shortcomings of genetic programming (GP) and genetic algorithm (GA), the GEP methodology, as seen in Figure 6, integrates key ideas from both methods. Further details are given in Khan et al. [18].
Figure 5. General structure of Multilayer Perceptron (moified from Hammad et al. [40]).
Figure 5. General structure of Multilayer Perceptron (moified from Hammad et al. [40]).
Atmosphere 14 00452 g005
Figure 6. General structure of Gene Expression Modeling.
Figure 6. General structure of Gene Expression Modeling.
Atmosphere 14 00452 g006

2.3. Pre-Processing Techniques

Maximal Overlap Discrete Wavelet Transformation (MODWT)

As in previous research work [18], Maximal Overlap Discrete Wavelet Transformation (MODWT) was employed to decompose the input time series into scaled components. This advanced variant has a clear advantage over simple DWT as it deals with down samples and reduces the number of errors brought on by time series over-translation by rectifying the boundary of the converted data series [50,51,52]. Raw input (xn) is processed via two complimentary filters before being divided into two components: Approximations (a) and Details (d). The low pass filter ignores the time series’ high frequency components and only permits high scale low frequency elements; these are thought to be more significant since they capture the identity and structure of the signal [53]. The white noise and subtlety inherent in the rainfall data are accepted by high pass filters in the same way as they are accepted by low scale high frequency components (also known as detail components) at each level of decomposition. Figure 7 graphically presents the detail and approximation components of rainfall timeseries observed in the Soan-Haro River basin. It can be seen that the detail components at all three levels include the white noise feature of the timeseries data, whereas the approximation component at the highest level (i.e., at level 3) reflects the main structure of the signal.
By including the detail components of each level of decomposition and the highest-level approximation components, the original signal may be recreated [36,38]. The tendency of wavelet pre-processing to decompose time series data into high frequency and low frequency series helps machine learning models to train on the obscure and salient features of the data, ultimately enhancing modeling accuracy. In the current study, we decomposed rainfall time series data into approximation and details signals. Revealing the hidden multiscale information of the data helps machine learning models to perform better in quantifying the relationship between rainfall and runoff.

2.4. Model Selection

We have used the five most used and acclaimed machine learning methods (DTF, TB, SDT, GEP, and MLP) to link rainfall and runoff in the Pothohar region of Pakistan. However, rainfall at a given time (Rt) has also been included in the input variables to account for ongoing events. Data on daily precipitation and runoff were collected from 1999 to 2016 (a 17-year period). The training set and testing set were each given a separate subset of the data: 70% for models’ training and 30% for validation and testing. Two different model types were created for each of the abovementioned methods since the study also attempted to evaluate the impact of input pre-processing on the models. Rainfall data for time “t” were first used as an input and runoff time series served as the training target. Using original time series, the first models for each approach were trained and tested. The inputs were then lagged up to lag order (Lo) 10, which resulted in the generation of numerous models for each approach. Optimum lag order for a basin is a function of basin characteristics and the hydrology of the area. The order of lag that produces the highest accuracy of rainfall-runoff modeling reveals important hydrological information about the basins, such as lag time or basin lag.
Model input without any lag order is Lo0; this includes rainfall time series at time t as the only variable. Similarly, input at Lo1 has variables Rt, and Rt−1 that represent the rainfall events lagged by 1-day. Likewise, Lo2 having variables Rt, Rt−1, and Rt−2 represents the rainfall events lagged up to 2-days. In a similar manner, Lo10 contains rainfall events lagged from day 1 to day 10. To explore the temporal association of rainfall, time lag (measured in days) was introduced into the models. The optimum lag order for each model was selected based on a trial-and-error method and was further consolidated using a partial correlation method (PCM) [54]. Ribeiro et al. [55] suggested particle swarm optimization (PSO) as an alternative method of input optimization. The formulated input combinations make certain that the effects of earlier events are carried over to the subsequent lag, while the discrete combinations are limited to rainfall events occurring on that single day.
The model inputs were pre-processed by MODWT after rudimentary models for each approach were trained and evaluated. The initial MODWT-based model was modified using all the daubechies wavelet filters in order to choose the best wavelet function and amount of decomposition, with support in the literature from previous studies [38,56]. Other filters could not match the performance of the wavelet function db4. It was further discovered that MODWT transformation with db4 at level 3 was the best choice after comparing the performance at all decomposition levels. Inputs generated for each approach were then transformed in this way. Other wavelet families include Symlets, Haar, Spline, Sine, and Coiflet; however, Daubechies (db) is recognized to perform better than other wavelet families in similar hydrologic issues. As a result, the wavelet function (db4) is chosen while the db family is chosen based on performance.
The optimum parameters selected for the MLP, GEP, and decision tree-based machine learning techniques, including SDT, DTF, and TB, were obtained by simultaneously training multiple settings using the par-for loop function in MATLAB parallel computing environment. A potential range of each parameter was given to the computing tool and parallel models was trained on unique parametric combinations using the grid search function. The parameter combinations with the highest accuracy were selected again for training and testing the models in this study. The selection of optimum lag order range was based primarily on the trial-and-error method and further validated by the partial correlation method. The overview of all the parameters used throughout the training and testing of all the approaches employed in this study is shown in Table 3. The confluence of trial-and-error and input optimization strategies recommended by prior research was used to choose the best settings.

2.5. Model Performance Indices

Numerous analytical techniques may be used to evaluate the performance of the models developed for the rainfall-runoff modeling in this study. The selection of performance evaluation indicators is significantly influenced by the range of computational models that are employed in the research, as well as their outcomes. We included empirical measurements like Coefficient of Determination (R2), Root Mean Squared Error (RMSE), and Nash Sutcliffe Efficiency (NSE) [57]. The following equations mathematically represent the selected performance indices.
R 2 = [ n   [ i = 1 n   ( Q o b s * Q m o d ) ( i = 1 n Q o b s )   ( i = 1 n Q m o d )   ]   ] 2 [   n i = 1 n Q o b s 2 ( Q m o d ) 2   ] [   n i = 1 n Q m o d 2 ( Q m o d ) 2   ]
R M S E = 1 n i = 1 n ( Q o b s Q m o d ) 2
N S E = 1 i = 1 n ( Q o b s Q m o d ) 2 i = 1 n ( Q o b s Q o b s ¯ )   2
where Q o b s and Q m o d denote the observed and modelled runoff using the developed models, respectively, and Q o b s ¯ is the average of observed runoff while n is number of observations.
The Value of R2 suggests the degree of similarity across the actual and predicted values. Higher numbers (near to 1) suggest more resemblance, while lower scores (closer to −1) reflect greater dissimilarity. The square root of the average value of squared errors, or RMSE, is a measurement of the mean gap between data points and the regression curve. The NSE is a metric that assesses a model’s capability to predict observed data. NSE has a spectrum of scores from −∞ to 1 [58]. A greater NSE score (approaching 1) indicates that the model can correctly predict, and vice versa. Generally, NSE scores are high and RMSE values are small in a well-performing model.

3. Results and Discussion

3.1. Bunha and Kahan River Basins

In the subsequent sections, a comprehensive description of the performance of the application of the selected ML techniques used as stand-alone models (Figure 8 for models training and Figure 9 for model testing) and in coupling with MODWT (Figure 10 refers to models train and Figure 11 to models test) for the combined Bunha and Kahan basin have been provided alongside an intercomparison of their performance assessment (shown in Figure 12, Figure 13 and Figure 14).

3.1.1. Single Decision Tree (SDT)

As observable in Figure 8, Figure 9, Figure 10 and Figure 11, MODWT-based SDT models surpassed the simple SDT model at all lag orders (from Lo0 to Lo10). However, it can also be observed that the difference between testing RMSE of MODWT-SDT models and stand-alone SDT models at lower lag orders was lower when compared to that at higher lag orders. Likewise, the gap between testing the NSE of stand-alone SDT models and MODWT-SDT models continued to decrease with an increase in lag order. At Lo0, the RMSE of the MODWT-SDT model is 0.58 times the RMSE value obtained from SDT models. At Lo10, the difference between the RMSE obtained from both models is insignificant (0.89 times). This finding indicates the probability of stand-alone SDT models outperforming MODWT-SDT at higher lag orders. Despite this, the maximum NSE in testing was observed in the MODWT-SDT-Lo4 model, which equaled 0.32. Moreover, Lo4 also generated outputs with the lowest RMSE value, equal to 344.48 mm. The mean NSE for SDT models was 0.21 in training and 0.19 in testing; for MODWT-SDT, it was 0.32 and 0.27, respectively. The average training and testing RMSE of SDT models was 498.72 mm and 539.86 mm, respectively. For MODWT-SDT models, the average training and testing RMSE was equal to 363.67 mm and 395.02 mm, respectively, indicating the effectiveness of wavelet pre-processing of raw data. The above statistics indicate that neither SDT nor MODWT-SDT models could produce satisfactory outcomes.

3.1.2. Decision Tree Forest (DTF)

DTF model performance is graphically expressed in Figure 8 and Figure 9. Figure 10 and Figure 11 demonstrate the performance of MODWT-DTF models, where it can be observed that MODWT-DTF bested DTF. Although simple DTF models yielded high modeling accuracy, averaging training NSE at 0.83 and testing NSE at 0.75, the performance indicators of MODWT-DTF models enhanced insignificantly due to less room for improvement in DTF models performance, having average training RMSE as low as 235.4 mm and average testing RMSE equal to 266.17 mm. The average NSE of MODWT-DTF in training was observed to be 0.88 and 0.86 in testing, while maximum efficiency (NSE = 0.86) was observed at Lo3. As observed, modeling efficiency remained consistent in a narrow range of 0.84 to 0.86 at lag orders above Lo3. As such, it is concluded that MODWT-DTF yielded the best performance at Lo3, and training NSE was also considerably higher (NSE = 0.87). However, it was also observed that higher lag orders yielded lower RMSE values (testing RMSE at Lo8 = 203 mm) when compared to lower lag orders (testing RMSE at Lo3 = 211 mm). However, this difference is negligible when compared to the increase in the number of input variables at Lo8. Maximum testing NSE, equal to 0.86, for stand-alone DTF models were also obtained for Lo6.
In conclusion, input pre-processing through MODWT improved DTF modelling accuracy and improved testing NSE by a factor of 1.14 on average, while in training, average NSE increased from 0.83 to 0.88 by a factor of 1.06 only. Graphical representations of performance indices illustrate that the trends of NSE values and RMSE had greater consistency during training and testing of MODWT-DTF models when compared to that of DTF models. The line graphs of coefficient of determination (R2) depict that during both training and testing of simple DTF, values of R2 were considerably low at Lo0. Starting from 0.65 in training and 0.18 in testing at Lo0, the values gradually spiked to a maximum of 0.87 in training and 0.86 in testing at Lo6 and showed no significant difference after lag order 4. Contrastingly, for MODWT-DTF models, testing values of R2 started at 0.85 at Lo0 and jumped to 0.96 at Lo3, consolidating the outcomes of NSE and RMSE at Lo3. The R2 values showed no spike or recession after Lo3. Therefore, MODWT-DTF-Lo4 is the best performing model among DTF and MODWT-DTF models. The comparative analysis of the best performing models from each technique also supported the observation that MODWT-DTF-Lo3 yielded the highest modeling accuracy for the Bunha-Kahan River basin among the 110 other models formulated and run on the study area sub-basin.

3.1.3. Tree Boost (TB)

After SDT and DTF, TB and MODWT-TB models were employed and performance indices were calculated from the modelled and target values. Figure 8 and Figure 9 show the training and testing performance of TB models in terms of NSE, RMSE (mm), and coefficient of determination (R2). The performance indices of MODWT-TB models are graphically presented in Figure 10 and Figure 11.
By comparing the performance of TB and MODWT-TB models, a considerable difference between the performance of pre-processed and un-processed models can be observed, as hybrid MODWT-TB models yield slightly greater accuracy in terms of NSE and R2 in comparison to stand-alone TB models. The highest efficiency of MODWT-TB in training was 0.54 at Lo1 and a maximum testing NSE equal to 0.46 was observed at Lo1 and Lo3. The average training and testing NSE for TB models was 0.28 and 0.25, respectively. Similarly, the average training and testing NSE for MODWT-TB models was 0.44 and 0.40, respectively; substantially higher than simple TB models. However, an insignificant decrease in RMSE was observed due to MODWT integration, as training average RMSE decreased from 479.892 mm in TB to 451.55 mm in MODWT-TB and testing mean RMSE decreased from 502.76 mm to 485.61 mm due to MODWT coupling with TB models. Minimum RMSE for training and testing of MDOWT-TB, equal to 381.33 mm and 417 mm, respectively, was observed at Lo1, while a negligible increase or decrease in RMSE was observed at higher lags. Therefore, Lo1 is the best performing model among TB and MODWT-TB models as a result of the lowest RMSE value in training and testing, the highest training and testing NSE values, and highest training and testing R2 values (training R2 = 0.56 and testing R2 = 0.47). The trends of R2 in the training and testing of TB and MODWT-TB models follow an inconsistent trend. R2 in training of TB is high at lower lags and decreases to a minimum of 0.20 at Lo8, while R2 values in testing also follow the similar trend throughout the lag order window. Similarly, in MODWT-TB models, training R2 values showed an inconsequential spike at Lo1, from 0.48 to 0.56, while testing values jumped from 0.41 to 0.47. However, both training and testing NSE MODWT-TB reach maximum R2 values at Lo1. TB and MODWT-TB models with high NSE and R2 values yielded low RMSE and models with low NSE and R2 yielded high RMSE values. As with increases in training accuracy (in terms of NSE), the testing NSE value of MODWT-TB models increased substantially. Hence, wavelet transformation enhanced the performance of stand-alone TB models, However, in comparison to other techniques, these models did not perform well.

3.1.4. Multi-Layer Perceptron (MLP)

Figure 8, Figure 9, Figure 10 and Figure 11 represent the stand-alone MLP and hybrid MODWT-MLP model performance in accordance with the performance indices. The NSE and RMSE values confirm that MODWT-MLP models marginally outperformed standalone MLP models, as the average training and testing NSE for stand-alone MLP models was 0.17 and 0.14, respectively, compared to that of 0.22 and 0.15 for MODWT-MLP, respectively.
MODWT coupling with MLP models could not enhance performance efficiency in training and testing up to desired limits. Moreover, at higher lag orders, the testing NSE of simple MLP models is either equal to or higher than that of MODWT-MLP models. However, at lag orders, such as Lo6, where MODWT-MLP testing NSE is slightly better than that of MLP, the RMSE of the former is significantly lower than that of the latter. For example, the testing NSE of MLP and MODWT-MLP is 0.07 and 0.10 at Lo6, respectively, But the testing RMSE at same lag order is 843 mm and 691 mm, respectively. Moreover, simple MLP models produced the highest testing NSE at L04, equal to 0.27, while MODWT-MLP yielded the best testing efficiency at Lo1, equal to 0.29. The highest training NSE value (0.33) was also observed at the same lag order for the MODWT-MLP model. Likewise, the lowest testing RMSE value (426 mm) was observed at Lo1. The average RMSE for MLP in training and testing was 625.22 mm and 659.44 mm compared to 501.34 mm and 581 mm for MODWT-MLP models, respectively. Likewise, the highest testing R2 value for MODWT-MLP was observed at Lo1 (0.31). Therefore, in this case, Lo1 showed the best performance. Additionally, wavelet pre-processing improved the performance of MLP models, although it was still insufficient.

3.1.5. Gene Expression Programming (GEP)

The performance of stand-alone GEP and hybrid MODWT-GEP models is compared graphically in Figure 8, Figure 9, Figure 10 and Figure 11, which amply demonstrate that MODWT-GEP hybrid models outperformed stand-alone GEP models. In testing, the average NSE for the wavelet-based hybrid GEP models was 0.44, as opposed to 0.29 for stand-alone GEP models. Additionally, MODWT-GEP yielded a testing RMSE that was on average lower (444.95 mm) than for GEP models (610.16 mm). The maximum NSE value measured during testing was 0.62 at Lo2 for MODWT-MLP, while the lowest RMSE, equal to 334.19 mm, was also measured at Lo2. The average RMSE for GEP in training and testing was 559.68 mm and 610.16 mm as compared to 423.60 mm and 444.95 mm for MODWT-MLP models, respectively. The average training and testing R2 value for GEP models was 0.36 and 0.29, respectively. Similarly, average training and testing R2 value for MODWT-GEP models was 0.49 and 0.44, respectively, considerably higher than simple GEP models. Gene expression programming, in contrast to TB, MLP, and SDT, generated results that were acceptable, while wavelet processing showed potential for improving accuracy.

3.1.6. Comparative Analysis

After a rigorous performance assessment, a comprehensive comparison analysis of the best and most consistent stand-alone and/or hybrid machine learning techniques was developed. For the thorough comparison analysis, one model with the best performance indicators was chosen from each approach employed in this study. The statistical indicators and wavelet parameters of the models are summarized in Table 4. It should be emphasized that basic models combined with wavelet pre-processing surpassed counter-models with raw inputs in each approach. Additionally, the performance statistics demonstrate that these models usually succeeded at lower lag orders, specifically Lo1, Lo2, Lo3, and Lo4. The wavelet-coupled decision tree forest (MODWT-DTF) model was shown to perform better than all other models in this experiment, according to subsequent evaluation. The Bunha-Kahan River basin’s highest testing NSE value (0.86) was obtained at lag order 3 (Lo3) by the hybrid MODWT-DTF model. This statistic reflects a minimal modeling error of 220.45 mm and an accuracy of 86% in terms of RMSE. As observable in Figure 12, the comparison between the performance accuracy of stand-alone machine learning techniques in Bunha-Kahan River basin clearly indicates that DTF-Lo6 outperformed the best models of other stand-alone techniques. GEP-Lo6 was second highest with an NSE value of 0.40, and SDT-Lo5 had the lowest NSE value (above 0.20) among all best performing models. Similarly, the intercomparison of best models from each stand-alone and hybrid ML technique is graphically shown in Figure 13. The figure shows that highest accuracy enhancement was observed in the case of GEP, where a 55% increase in NSE was observed due to MODWT transformation. However, the smallest NSE increase was observed in the DTF model, which yielded above 0.8 NSE in both stand-alone and hybrid modes. Interestingly, it was also observed that MODWT transformation decreased the optimum lag order when compared to stand-alone models. The evaluation of models with the best accuracy for each approach is shown visually in Figure 14. The detailed description of performance indices of stand-alone and hybrid machine learning models is presented in Table A1 in Appendix A.
As was previously mentioned, the MODWT-DTF model generated the highest modelling accuracy (around 86%) at Lo3. The regression plot of the MODWT-DTF-Lo3 model, displayed in Figure 15, demonstrates a high R2 value of 0.92, that was previously listed in Table 4. Figure 16 compares the output of the model with the actual gauge reading recorded at the gauge station. MODWT-DTF-Lo3, represented by the orange-colored solid line, closely replicated the measured runoff, represented by the black-colored solid line. The model was successful at capturing the patterns of actual runoff at both frequent and extreme levels (outliers).

3.2. Soan and Haro River Basins

The following sections provide a thorough explanation of the modelling outcomes obtained from the application of stand-alone (Figure 17 and Figure 18 refer to models’ training and testing, respectively) and MODWT-hybrid (Figure 19 for models’ training and Figure 20 for models’ testing) machine learning models for the combined Soan and Haro basin have been provided. The intercomparison of their performance assessment is also presented in the subsequent sections (Figure 21, Figure 22 and Figure 23).

3.2.1. Single Decision Tree (SDT)

It can be seen in Figure 17, Figure 18, Figure 19 and Figure 20 that MODWT-based SDT models outperformed simple SDT model at all lag orders (from Lo0 to Lo10). However, it can also be observed that the difference between the testing RMSE of MODWT-SDT models and stand-alone SDT models at lower lag orders was lower when compared to that at higher lag orders. Likewise, the gap between the testing NSE of stand-alone SDT models and MODWT-SDT models continued to decrease with increase in lag order. At Lo0, the RMSE of MODWT-SDT model is 0.76 times that of the RMSE value obtained from SDT models. Meanwhile, at Lo10, the difference between the RMSE obtained from both models is insignificant (0.91 times). In this case, the MODWT-SDT-Lo3 model had the highest testing NSE, equivalent to 0.31. Additionally, Lo3 produced outputs with the lowest RMSE value of 63.02 m3/s. In training, the mean NSE for SDT models was 0.23, while in testing it was 0.18. However, MODWT-SDT models produced 0.27 NSE in testing and 0.32 in training. SDT models had an average training and testing RMSE of 72.64 mm and 76.76 mm, respectively. The average training and testing RMSE for MODWT-SDT models was 61.26 mm and 65.47 mm, respectively, demonstrating the usefulness of wavelet pre-processing of original data. However, neither the SDT nor the MODWT-SDT models could produce adequate results.

3.2.2. Decision Tree Forest (DTF)

Figure 17 and Figure 18 illustrate the performance of DTF models graphically, while Figure 19 and Figure 20 depict the efficiency of MODWT-DTF models. These figures show that MODWT-DTF models performed better than DTF models. The performance metrics of MODWT-DTF models noticeably improved even though simple DTF models produced acceptable modelling accuracy, with an average training NSE of 0.63 and an average testing NSE of 0.56; these models also had an average training RMSE as low as 42.81 m3/s and an average testing RMSE equal to 53.21 m3/s. This was because there was still an opportunity for enhancing the effectiveness of DTF models. An average NSE of MODWT-DTF in training was observed to be 0.87 and 0.84 in testing, while maximum efficiency (NSE = 0.89) was observed at Lo5. However, Lo4 also yielded a high NSE value, equal to 0.88, marginally less than the NSE at Lo5. Moreover, Lo4 produced the lowest testing RMSE value, equal to 21.71 m3/s, when compared to Lo5 and other lag orders. This justifies the selection of MODWT-DTF-Lo4 as the best performing model among DTF and MODWT-DTF models in the Soan-Haro River basin. As observed, modeling efficiency remained consistent in a narrow range of 0.88 to 0.89 at lag orders above Lo4; as such, it is concluded that the MODWT-DTF model yielded the best performance at Lo4 with an enhancement in the NSE performance in testing and training by a factor of 1.49 and 1.37 on average, respectively.
Graphical representations of performance indices show that the trends of NSE values and RMSE had greater consistency during the training and testing of MODWT-DTF models when compared to simple DTF models. The line graphs of coefficient of determination (R2) depict that during both training and testing of the simple DTF model, values of R2 were considerably low at Lo0. Starting at 0.40 in training and 0.33 in testing at Lo0, the values gradually spiked to a maximum of 0.75 in training and 0.68 in testing at Lo5. They showed no significant difference after lag order 5. Contrastingly, for MODWT-DTF models, training, and testing values of R2 started at 0.67 and 0.57 at Lo0 and jumped to 0.95 and 0.91 at Lo4, respectively, consolidating the outcomes of NSE and RMSE at Lo4. The R2 values showed no significant spike or recession after Lo4. Therefore, the MODWT-DTF-Lo4 model is the best performing model among DTF and MODWT-DTF models. The comparative analysis of best performing models from each technique also supported the observation that MODWT-DTF-Lo4 yielded the highest modeling accuracy in the Soan-Haro River basin among all 110 other models formulated and run on the study area sub-basin.

3.2.3. Tree Boost (TB)

Figure 17 and Figure 18 show the training and testing performance of TB models in terms of NSE, RMSE (m3/s), and coefficient of determination (R2), respectively. The performance indices of MODWT-TB models are similarly presented in Figure 19 and Figure 20. By comparing the performance of TB and MODWT-TB models, a significant difference between the performance of pre-processed and un-processed models can be observed, as MODWT-TB models yielded greater accuracy in terms of NSE and R2 when compared to simple TB models. The highest performance efficiency of MODWT-TB in training was 0.66 at Lo2 and the highest testing NSE, equal to 0.54, was also observed at Lo2. The average training and testing NSE for TB models was 0.28 and 0.18, respectively. Similarly, the average training and testing NSE for MODWT-TB models was 0.52 and 0.40, respectively, significantly greater than simple TB models. Moreover, a similar trend of a decrease in RMSE was observed due to MODWT integration, as training average RMSE decreased from 54.51 m3/s in TB to 35.67 m3/s in MODWT-TB and testing mean RMSE decreased from 51.25 m3/s to 29.46 m3/s due to MODWT coupling with TB models. However, minimum RMSE for training and testing of MDOWT-TB, equal to 29.49 m3/s and 35.67 m3/s, respectively, was observed at the highest lag order (Lo10), while at Lo2, a slightly higher testing RMSE, equal to 42.34 m3/s, was observed. Therefore, Lo2 is the best performing model among TB and MODWT-TB models due to a low RMSE value in training and testing, highest training and testing NSE values, and second-highest training and testing R2 values (training R2 = 0.67 and testing R2 = 0.60). The trends of R2 in training and testing of TB and MODWT-TB models follow an inconsistent trend. R2 in the training of the TB model is higher at lower lags and decreases to a minimum of 0.23 at Lo8, while R2 values in testing also follow a similar trend throughout the lag order window. Similarly, in MODWT-TB models, training R2 values show a significant spike at Lo1, from 0.54 to 0.73, while testing values jump from 0.48 to 0.64. Also, both training and testing values for the NSE MODWT-TB model reach the highest R2 value at Lo1. TB and MODWT-TB models with high NSE and R2 values yielded low RMSE and models with low NSE and R2 yielded high RMSE values. As with the increase in training accuracy (in terms of NSE), the testing NSE value of MODWT-TB models increased substantially. As such, wavelet transformation improved the performance of simple TB models. In comparison to other techniques, these models performed comparatively well.

3.2.4. Multi-Layer Perceptron (MLP)

Figure 17, Figure 18, Figure 19 and Figure 20 represent the MLP and MODWT-MLP models’ performance in accordance with the performance evaluation parameters, respectively. The graphs show that MODWT-MLP models performed marginally better than stand-alone MLP models. This is supported by NSE and RMSE values. The average training and testing NSE for standalone MLP models was 0.17 and 0.14, compared to 0.22 and 0.15 for MODWT-MLP, while stand-alone MLP models had NSE of 0.17 and 0.14, respectively.
MODWT coupling with MLP models could not enhance performance efficiency in training and testing up to desired limits. Moreover, at higher lag orders, the testing NSE of simple MLP models was either equal to or higher than that of MODWT-MLP models. However, at lag orders such as Lo6, where MODWT-MLP testing NSE was slightly better than that of MLP, the RMSE of the former was significantly less than that of the latter. For example, the testing NSE of MLP and MODWT-MLP is 0.07 and 0.10 at Lo6, respectively. But the testing RMSE at the same lag order is 843 mm and 691 mm, respectively. Moreover, simple MLP models produced the highest testing NSE at L04, equal to 0.27, whereas MODWT-MLP yielded the best testing efficiency at Lo1, equal to 0.29. The highest training NSE value (0.33) was observed at the same lag order for the MODWT-MLP model. Likewise, the lowest testing RMSE value (426 mm) was observed at Lo1. The average RMSE for MLP in training and testing was 625.22 mm and 659.44 mm, as compared to 501.34 mm and 581 mm for MODWT-MLP models, respectively. Likewise, the highest testing R2 value for MODWT-MLP was observed at Lo1 (0.31). We may draw the conclusion that in this case, Lo1 showed the best performance. Additionally, wavelet pre-processing improved MLP model performance, although it was still unsatisfactory.

3.2.5. Gene Expression Programming (GEP)

The modelling accuracy of GEP and MODWT-GEP models is shown graphically in Figure 17, Figure 18, Figure 19 and Figure 20. This clearly shows that MODWT-GEP models surpassed GEP models. Wavelet pre-processed GEP models provided runoff simulations with greater accuracy than that of basic GEP models. The average NSE in testing for the former was 0.22 as opposed to 0.15 for the latter. Additionally, MODWT-GEP generated a testing RMSE that was lower on average (63.55 m3/s) when compared to GEP models (76.18 m3/s). While lag order 1 produced runoff values with the highest RMSE values, the lowest RMSE (59.58 m3/s) was recorded at lag order 4 for the same model. The maximum NSE value in testing was obtained at Lo1 for MODWT-MLP, equal to 0.30. The average RMSE for GEP in training and testing was 62.34 m3/s and 76.18 m3/s, as compared to 55.66 m3/s and 63.55 m3/s mm for MODWT-MLP models, respectively. The average training and testing R2 value for GEP models was 0.45 and 0.33, respectively. Similarly, average training and testing R2 values for MODWT-GEP models was 0.61 and 0.50, respectively. These were considerably higher than simple GEP models. Gene expression programming produced acceptable outcomes, in contrast to TB, MLP, and SDT, while wavelet processing showed the potential to improve modelling accuracy.

3.2.6. Comparative Analysis

In this part, a comprehensive comparison analysis of the best and most consistent models for each approach is developed after a rigorous performance assessment of each technique with and without wavelet transformation. For the thorough comparison analysis, one model with the best performance indicators is chosen from each approach employed in the Soan-Haro River basin. The statistical indicators and wavelet parameters of the models are summarized in Table 5. It should be noted that basic models combined with wavelet pre-processing surpassed counter-models with raw inputs in each approach. Furthermore, the findings reflect that these models excelled at lower lag orders. The MODWT-DTF model was shown to perform better than all other models in this study, according to extensive examination. At lag order 4 (Lo4), the MODWT-DTF model generated the highest testing NSE value (0.88) for rainfall-runoff modelling of the Soan-Haro River basin. The model had a minimum modeling error of 21.72 m3/s and a maximum modeling accuracy of 88%. As can be seen in Figure 21, the comparison between the performance accuracy of stand-alone machine learning techniques in the Soan-Haro River basin clearly reflects that DTF-Lo5 yielded the highest NSE value (training = 0.71, testing = 0.65) among best models of other stand-alone techniques. TB-Lo2 was second highest. With a testing NSE value of 0.45, and MLP-Lo2 had the lowest NSE value (0.20) among the best performing models. The intercomparison of the best models from each stand-alone and hybrid ML technique is graphically shown in Figure 22. The figure shows that highest NSE increment was observed in case of GEP, despite being second lowest overall; a 43% increase in NSE was observed due to MODWT transformation. However, despite the high percentage increase in NSE value due to MODWT transformation, the MODWT-GEP-Lo1 model did not yield satisfactory or acceptable modelling accuracy. Furthermore, the lowest NSE increase was observed in the TB model, which yielded a testing NSE value of 0.45 in stand-alone mode and 0.54 in hybrid application. As in the first scenario, MODWT transformation decreased the optimum lag order when compared to stand-alone models in the case of the Soan-Haro River basin.
The comparison of the models with the best accuracy for each approach is represented visually in Figure 23. Table A2 in Appendix B provides the tabulation of detailed statistical metrics of model performance in the Soan-Haro River basin.
As discussed above, the MODWT-DTF model produced the highest modeling accuracy (88%) at Lo4. Figure 24 represents the regression plot of the MODWT-DTF-Lo4 model in the Soan-Haro River basin while Figure 25 compares the model’s rainfall-runoff plot to actual runoff. It can be seen that the MODWT-DTF-Lo4 model (orange line) effectively mirrored the measured runoff (black line). The model accurately depicted both recurrent frequencies and patterns of observed runoff at extreme levels (outliers).

4. Discussion

The performance of five widely used machine learning approaches is comprehensively compared in this work alongside the impact of MODWT pre-processing on rainfall-runoff modeling in the Pothohar region of Pakistan. This research work is a follow up to our previous work [18], in which the above-mentioned techniques were applied to model rainfall-runoff relationship in the Soan River basin only. Here, the authors have studied the links between rainfall and runoff in the entire Pothohar region.
In particular, the study revealed that the MODWT-DTF model performed best with an NSE equal to 0.88 and R2 equal to 0.91 at Lo4 in the Soan-Haro River basin. In contrast to our study, Hussain et al. [58] modeled the rainfall-runoff relation of this region with a maximum R2 value of only 0.66. Ghani et al. [59] used HEC-HMS and HEC-GeoHMS software to map the rainfall-runoff connection of the Pothohar region and the Soan River basin with 92% accuracy. The results of our study offer quantitative evidence for the conclusions regarding machine learning model efficiency increases when combined with wavelet pre-processing.
Additionally, the Soan-Haro River basin’s optimal lag order was found to be four, meaning that the models worked effectively when rainfall data was delayed up to four days with input variables Rt−1, Rt−2, Rt−3, and Rt−4. The time it takes for runoff generated due to rainfall in the farthest part of the basin to reach the basin outlet is related to this optimum lag number. Examination of the rainfall and runoff dataset revealed that the highest runoff (1407.69 m3/s) occurred four days after the highest rainfall event (68.70 mm). This was in contrast to the testing period, when the highest runoff was discovered to occur after three days and lasted up to four days after the maximum rainfall event. Previous investigations [36,37,60] have also noted this relationship between the theoretical idea of lag order, sometimes known as time lag or lag number, and the physical concept of time of concentration (Tc).
The single decision tree models are largely prone to lose information in tasks requiring multiple categorizations. This makes SDT less accurate in learning and predicting the relationship between continuous variables [61]. Moreover, the decision tree-based models, including SDT and TB, are highly unstable and overly sensitive to minor changes in time series data. As such, modelling performance is significantly impacted by the range of time series data. In contrast, DTF combines the outcomes from several decision trees and the outcome is based on an ensembled result of variable categorization in all the trees. Moreover, the decision tree forest also creates more diverse sets of decision trees that capture different features in the data [62]. This makes DTF less likely to be sensitive to a minor variation and overfitting problem and therefore more accurate in prediction tasks when compared to SDT and TB. Similar trends have been observed in the findings of our study, as DTF and MODWT models clearly outperformed stand-alone and MODWT-based SDT and TB models in both scenarios. Similarly, gene expression programming (GEP) and multi-layer perceptron (MLP) models can sometimes lack generalization and interpretability of data features in prediction tasks [63,64]. Moreover, both GEP and MLP are prone to overfitting in cases of diverse hyperparameter settings. Lower NSE values of stand-alone and hybrid MLP and GEP models are consistent with these caveats of the techniques.
This work contributed to the body of knowledge by demonstrating the usefulness of DTF and MODWT-based DTF in the hydrological modeling of the Pothohar region, while other earlier studies [28,38,40,51,65] support the conclusions of our study regarding improving modeling accuracy using wavelet transformation.
Nevertheless, it may be noted that, although the Pothohar region is one of the main basins in the Indus River basin, it does not entirely represent the hydrology of the Indus River plane, since there exist other portions of the Indus basin.
A combination of evapotranspiration, rainfall, and infiltration data can also be used in any future research work to better comprehend the association between these hydrological processes. Moreover, the hydrological response of catchment basins is expected to become more complicated under the influence of changing climate trends. Therefore, multivariate quantification of runoff would be helpful in solving prospective hydrological problems. Furthermore, this experiment only utilizes rainfall and pre-processed rainfall data to retrieve obscured details and additional information from the time series. Long short-term memory (LSTM), a deep learning approach, may be useful in examining the novel aspects of the hydrological parameters of the basin since deep learning applications in hydrology have progressed significantly.
These recommendations for further study effort are highly likely to be implemented by the authors in subsequent studies.

5. Conclusions

This study employed five machine learning techniques (SDT, TB, DTF, MLP, and GEP) as stand-alone models and hybrid models coupled with maximal overlap discrete wavelet transformation (MODWT) to garner the relationship between rainfall and runoff in the Pothohar region of Pakistan. The study area consisted of six sub-basins, from which four major sub-basins were selected as the focused study area. The analysis was run in two different homogeneous areas including the combined Bunha-Kahan River basin and the Soan-Haro River basin. The rainfall data of each combined river basin were aggregated against the water level data (converted from m to mm) from each gauging station, located in the downstream river basin.
In each case, the analysis was initially run using standalone ML technique (SDT, TB, DTF, MLP, and GEP) settings and then coupled with MODWT. The findings reflected that the MODWT-based DTF model at lag order 3 (Lo3) performed best in the Bunha-Kahan River basin and at Lo4 in the Soan-Haro River basin. The modeling accuracy ranged from 0.85 to 0.90, in terms of NSE, and from 0.90 to 0.92, in terms of R2. Additionally, it was proven that wavelet pre-processing improved the modeling accuracy of stand-alone machine learning models.
In the Bunha-Kahan River basin, the DTF mode, in combination with MODWT, performed best at Lo3 with an NSE value of 0.86. RMSE and R2-values were 220.45 and 0.92, respectively. The MODWT-GEP model showed the second-best performance, with an NSE value of 0.62 at Lo2, followed by an NSE value of 0.46 for the MODWT-TB model at Lo1. In the Soan-Haro River basin, the MODWT-based DTF model again performed best, with an 88% accuracy at Lo4. The MODWT-TB model recreated runoff data with 54% accuracy at lag order 2 (Lo2). The Decision Tree Forest (DTF) model performed well in each of the sub-basins when compared to standalone counter-models. When coupled with MODWT, the performance of all ML models was enhanced. However, other than the MODWT-DTF model, the modeling accuracy of SDT, TB, MLP, and GEP was generally unsatisfactory, except MODWT-GEP-Lo2 in the Bunha-Kahan River basin (NSE = 0.62).
This performance improvement is a result of the wavelet transformation’s propensity to recover latent details from the univariate time series during input signal pre-processing. Therefore, MODWT may be applied with confidence to issues requiring many variables to provide better results while using less parametric information and sparse modeling. The type, geography, and characteristics of the basin might affect the wavelet technique’s findings. Thus, it is important to choose the right wavelet family and type as well as decomposition level in order to obtain the best results.

Author Contributions

Conceptualization, M.S.; methodology, M.S., M.T.K. and R.A.; software, M.T.K., H.S. and M.H.; validation, M.S. and M.A.I.; formal analysis, M.T.K., M.H., M.U.A. and H.S.; investigation, M.S. and M.T.K.; resources, M.T.K., M.S., R.A., S.A. and H.S.; data curation, H.S., M.H. and M.K.U.; writing—original draft preparation, M.T.K., M.H. and H.S.; writing—review and editing M.T.K., R.A. and S.A.; visualization, R.A., S.A., M.A.I., S.H. and M.U.A.; supervision, M.S., R.A. and M.A.I.; project administration, R.A.; funding acquisition, R.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the project POR FESR On DEmand Services for Smart Agriculture (ODESSA) held by Raffaele Albano.

Data Availability Statement

The authors are grateful to Pakistan Meteorological Department (PMD) for the provision of data for the study. The data may be obtained directly from PMD.

Acknowledgments

The authors are grateful to the reviewers and editor for their insightful comments, and Pakistan Meteorological Department (PMD) for the provision of data for the study.

Conflicts of Interest

The authors report no conflict of interest.

Appendix A. Model Performance Indices in the Bunha-Kahan River Basin

Table A1. Performance indices of stand-alone and hybrid machine learning models in the Bunha-Kahan River basin.
Table A1. Performance indices of stand-alone and hybrid machine learning models in the Bunha-Kahan River basin.
Lag Order (Lo) SDT MODWT-SDT
Training Testing Training Testing
NSERMSE *MAE **R2NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2
Lo00.16635.00300.020.170.13756.00357.190.140.29393.00185.680.280.24446.00210.720.25
Lo10.17597.65283.750.190.16614.25291.630.120.30379.90180.370.300.29369.10175.240.28
Lo20.21483.81230.830.220.19517.26246.780.180.39292.23139.420.350.28382.29182.390.27
Lo30.19534.74256.380.200.18546.00261.780.190.35327.06156.810.330.27394.66189.220.27
Lo40.20508.00244.770.230.13706.00340.170.200.33344.48165.980.340.32339.72163.680.31
Lo50.22472.40228.740.230.22452.45219.080.230.33345.36167.230.320.28382.29185.110.29
Lo60.24416.31202.590.240.22446.73217.400.220.29393.00191.250.270.26414.97201.940.26
Lo70.20508.00248.450.190.20481.90235.690.130.30379.90185.800.290.28382.29186.970.29
Lo80.22461.82227.010.230.20481.90236.880.220.31372.58183.140.330.30356.80175.390.31
Lo90.25406.40200.780.280.21468.00231.210.230.29393.00194.160.300.23465.39229.930.24
Lo100.22461.82217.140.240.21468.00220.050.210.30379.90178.620.290.26411.69193.570.25
Lag Order (Lo) DTF MODWT-DTF
Training Testing Training Testing
NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2
Lo00.65362.00171.030.670.18489.50231.270.180.87210.2099.310.980.85223.00105.360.92
Lo10.75306.00144.580.820.59344.20162.620.730.87206.4097.520.980.85217.00102.530.93
Lo20.81262.00124.390.900.69300.90142.860.840.86206.3097.950.980.85211.00100.180.93
Lo30.85227.00108.300.950.82233.00111.160.920.87205.1097.850.980.86220.00104.960.92
Lo40.86212.00101.640.970.85231.00110.750.920.87205.6098.570.990.86218.00104.520.92
Lo50.86207.5099.980.980.85230.00110.820.930.88203.5098.050.980.86213.00102.630.93
Lo60.86205.0099.260.980.86225.00108.950.930.88204.9099.220.980.86215.00104.110.93
Lo70.87202.9098.740.980.86223.00108.520.940.89204.8099.660.980.86206.00100.250.93
Lo80.87202.0098.800.980.86219.00107.110.930.89201.0098.310.980.85203.0099.280.94
Lo90.87201.6099.100.980.86216.00106.180.940.89204.20100.380.980.86205.00100.770.93
Lo100.87200.7099.160.980.86216.00106.710.950.89204.00100.790.980.86204.00100.790.93
Lag Order (Lo) TB MODWT-TB
Training Testing Training Testing
NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2
Lo00.23539.00253.430.240.19594.00279.290.200.44468.00220.050.480.38505.00237.440.41
Lo10.35384.20181.520.370.28423.07199.890.300.54381.33180.170.560.46417.17197.100.47
Lo20.36374.36177.740.350.27438.00207.950.280.51403.76191.690.500.45426.44202.460.47
Lo30.35384.20183.300.360.29409.17195.210.310.48392.00187.020.500.46417.17199.030.46
Lo40.26506.81242.990.290.23510.70244.850.240.44417.00199.930.450.41468.05224.410.41
Lo50.33405.67195.460.340.28424.52204.540.290.49411.00198.030.510.44436.14210.140.42
Lo60.25500.25242.230.360.24515.88249.800.250.41453.00219.350.430.39492.05238.260.39
Lo70.35372.00181.030.360.33374.20182.100.310.37516.00251.110.380.35548.29266.820.35
Lo80.20649.85317.830.240.18647.00316.440.200.39528.00258.240.400.34564.41276.040.36
Lo90.22567.43278.920.230.21583.50286.820.210.40492.00241.840.410.37518.65254.940.36
Lo100.21594.30293.610.200.20610.33301.530.200.37505.00249.490.370.35548.29270.880.35
Lag Order (Lo) MLP MODWT-MLP
Training Testing Training Testing
NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2
Lo00.18598.00293.610.200.16647.00304.210.170.25448.00448.190.260.18503.18236.590.20
Lo10.23469.57281.170.220.22472.73223.350.200.33369.00369.620.350.29426.45201.480.31
Lo20.21514.29221.860.220.20520.00246.880.210.28385.00385.460.300.23474.22225.150.20
Lo30.26415.38244.170.280.20520.00248.090.200.29404.00404.770.270.20496.76237.000.22
Lo40.29372.41198.180.320.27385.19184.680.230.39452.00452.710.400.24515.63247.220.22
Lo50.14671.43178.550.180.09655.56315.860.070.20528.00528.140.240.11632.94304.960.15
Lo60.11781.82323.510.100.07842.86408.130.050.23576.00576.740.260.10691.16334.670.12
Lo70.13730.77378.570.140.10868.00420.300.090.11548.00548.10.120.09676.42327.540.11
Lo80.09772.63353.850.110.05811.36394.840.110.14566.005660.120.10600.35292.150.10
Lo90.14671.43375.990.150.11645.45315.680.120.09645.00645.50.100.04723.59353.900.07
Lo100.10880.00328.390.090.07885.71428.880.070.13590.00590.350.110.11655.12317.220.07
Lag Order (Lo) GEP MODWT-GEP
Training Testing Training Testing
NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2
Lo00.37498.00242.350.400.33525.00255.490.340.49406.00197.920.510.44430.00209.370.44
Lo10.26666.38322.670.280.21767.86371.810.250.51394.12190.840.510.46415.65201.270.47
Lo20.35520.74253.410.340.32538.28261.950.260.64334.38162.720.660.62334.19162.630.58
Lo30.43442.47216.410.430.40450.63220.400.370.45433.33211.930.510.41454.15222.120.50
Lo40.57358.35176.150.580.44418.75205.840.400.60350.00172.040.620.55364.00178.930.49
Lo50.39477.59235.950.400.36489.58241.880.400.45433.33214.090.440.40463.00228.740.39
Lo60.44434.68204.380.430.40450.63211.880.370.51394.12185.310.490.50390.40183.560.47
Lo70.28625.93295.730.300.20701.25331.320.220.48412.50194.890.480.46415.65196.380.44
Lo80.30590.87286.110.300.17725.00351.060.190.45433.33209.830.440.44430.00208.210.44
Lo90.19875.05425.830.210.15935.00455.010.150.31583.87284.130.360.24705.00343.080.28
Lo100.26666.38325.920.250.23709.78347.140.190.39484.62237.020.400.37492.43240.840.35
* RMSE is in 000 mm. ** MAE is in 000 mm.

Appendix B. Model Performance Indices in the Soan-Haro River Basin

Table A2. Performance indices of stand-alone and hybrid machine learning models in the Soan-Haro River basin.
Table A2. Performance indices of stand-alone and hybrid machine learning models in the Soan-Haro River basin.
Lag Order (Lo) SDT MODWT-SDT
Training Testing Training Testing
NSERMSE *MAE **R2NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2
Lo00.0976.9323.430.090.0881.5024.820.080.2356.0217.060.230.2061.6318.770.20
Lo10.2273.8822.570.220.1876.7723.450.180.2864.5419.720.280.2766.4720.310.27
Lo20.2373.4722.520.230.2075.7623.220.200.3061.0718.720.300.2764.9519.900.27
Lo30.2569.0621.230.250.1975.6923.270.190.3359.3418.240.330.3163.0319.380.31
Lo40.2970.2521.670.290.2573.6022.700.250.3463.5419.600.340.3067.0420.680.30
Lo50.2473.4522.720.240.1877.8824.100.180.3558.7018.160.350.2965.2920.200.29
Lo60.2473.2022.720.240.1877.3124.000.180.3459.6418.510.340.2764.2019.930.27
Lo70.1776.5323.830.170.1477.7424.210.140.3460.2418.760.340.2663.6119.810.26
Lo80.2074.1823.170.200.1676.1223.780.160.3360.9019.020.330.2664.8520.260.26
Lo90.3269.4521.760.320.2275.6123.690.220.3264.9120.340.320.2569.6521.830.25
Lo100.3268.6020.830.320.2376.3423.180.230.3464.9219.710.340.2669.4921.100.26
Lag Order (Lo) DTF MODWT-DTF
Training Testing Training Testing
NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2
Lo00.3868.8820.980.400.3276.9723.440.330.6635.3210.760.670.5738.5711.750.57
Lo10.4959.9618.260.500.4168.8720.970.400.8221.906.670.770.7824.387.420.70
Lo20.5347.2614.440.550.4361.7418.860.490.8532.659.970.890.7836.0211.000.83
Lo30.6440.7012.470.630.5653.0716.260.600.9019.886.090.950.8721.766.670.90
Lo40.6737.7311.090.680.5950.1115.410.620.9018.885.800.950.8821.726.610.91
Lo50.7136.0911.130.750.6544.6013.760.680.9018.355.660.950.8922.306.880.90
Lo60.7135.8211.080.740.6445.8414.180.680.9018.315.660.950.8922.576.980.92
Lo70.7136.3211.270.740.6545.9814.270.680.9118.065.610.960.8922.376.940.91
Lo80.7036.3811.330.740.6546.0214.330.680.9018.105.640.960.8922.416.980.91
Lo90.7035.4911.190.740.6445.7314.280.670.9018.085.650.950.8922.186.930.92
Lo100.7136.2411.360.740.6446.4014.540.680.9118.105.670.960.8922.347.000.91
Lag Order (Lo) TB MODWT-TB
Training Testing Training Testing
NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2
Lo00.0652.8416.040.400.0256.6017.180.340.1335.6810.830.540.1037.8411.490.48
Lo10.4749.2214.990.700.3851.5515.700.610.5240.3512.290.730.4942.2912.880.64
Lo20.4748.9714.960.660.4549.4715.110.620.6640.2612.300.670.5442.3412.940.60
Lo30.2153.5316.410.410.1055.7917.100.310.5031.369.610.640.4135.8010.970.46
Lo40.2748.3414.860.350.0955.6517.110.300.5031.409.650.600.3636.1211.100.40
Lo50.2750.9015.700.390.1155.4917.110.300.4132.369.980.550.3736.2911.190.44
Lo60.2751.2215.850.430.1654.8816.980.290.5330.519.440.570.4235.9111.110.44
Lo70.2848.6615.100.400.1254.9817.060.260.5431.089.650.550.4035.7411.090.43
Lo80.2249.4715.400.400.1155.0417.140.230.6029.429.160.590.3536.1111.240.36
Lo90.3451.0315.940.470.2654.3016.960.380.6830.429.500.620.4735.9111.220.49
Lo100.2751.2516.060.470.2154.5117.080.350.6029.499.240.580.4535.6711.180.48
Lag Order (Lo) MLP MODWT-MLP
Training Testing Training Testing
NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2
Lo00.1965.4519.870.190.1485.2325.880.150.2469.2521.020.250.2172.8822.130.21
Lo10.2858.5617.830.280.1981.3424.770.200.3166.0320.110.310.2578.9424.040.25
Lo20.2171.3921.810.220.2081.5524.910.200.2967.7720.700.290.2479.6124.320.25
Lo30.2464.2119.680.250.1584.5625.910.150.2662.1519.050.290.2167.8620.800.25
Lo40.2463.4919.520.250.1584.1325.860.160.2860.6018.630.310.2365.8820.250.25
Lo50.1671.6222.090.160.1285.6626.420.130.2457.2517.660.270.2068.3121.070.24
Lo60.1873.2922.680.180.1683.8425.940.170.2462.0219.190.260.1969.4621.490.22
Lo70.1770.8121.910.180.1484.9526.280.150.2857.4917.790.290.2066.7820.660.25
Lo80.1969.3921.540.190.1285.3026.470.150.2558.3118.100.260.2167.5220.960.22
Lo90.1974.1023.070.200.1484.1526.200.160.2459.3018.460.250.2168.3321.280.21
Lo100.1776.3423.620.180.1484.7426.220.150.2661.1918.930.250.2268.9121.320.22
Lag Order (Lo) GEP MODWT-GEP
Training Testing Training Testing
NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2NSERMSEMAER2
Lo00.2457.4517.830.540.1676.4623.730.330.2762.1119.280.550.2465.2920.260.48
Lo10.2754.0216.710.560.2073.6622.790.420.3260.2218.630.650.3070.1621.710.60
Lo20.2164.9820.170.450.2173.6122.850.440.3360.0818.650.660.2870.8021.970.58
Lo30.2259.6818.580.460.1676.2023.730.330.2856.0017.440.630.2261.4519.130.54
Lo40.2358.9218.400.480.1875.3423.530.380.2956.0217.500.640.2459.5818.610.54
Lo50.1864.4220.190.350.0978.8124.700.210.2751.3416.090.610.2161.8319.380.49
Lo60.2165.9820.030.410.1774.9622.760.370.2954.9516.680.610.2361.1618.570.52
Lo70.2063.3219.280.420.1376.9923.450.300.3051.8515.790.630.1961.4018.700.48
Lo80.2162.4519.320.420.1476.7723.750.330.2852.2016.150.590.1962.4619.320.41
Lo90.2166.4020.610.450.1277.3524.010.280.2852.8316.400.590.2161.7019.150.44
Lo100.2068.0621.190.420.1277.8224.230.270.2554.7017.030.530.1663.2819.700.38
* RMSE in m3/s. ** MAE in m3/s.

References

  1. Montanari, A.; Young, G.; Savenije, H.H.G.; Hughes, D.; Wagener, T.; Ren, L.L.; Belyaev, V. “Panta Rhei-Everything Flows”: Change in hydrology and society-The IAHS Scientific Decade 2013–2022. Hydrol. Sci. J. 2013, 58, 1256–1275. [Google Scholar] [CrossRef]
  2. Ramly, S.; Tahir, W. Application of HEC-GeoHMS and HEC-HMS as Rainfall–Runoff Model for Flood Simulation. In International Symposium on Flood Research and Management [Internet]; Tahir, W., Abu Bakar, P., Wahid, M., Mohd Nasir, S., Lee, W., Eds.; Springer: Singapore, 2016. [Google Scholar] [CrossRef]
  3. Masseroni, D.; Cislaghi, A.; Camici, S.; Massari, C.; Brocca, L. A reliable rainfall-runoff model for flood forecasting: Review and application to a semi-urbanized watershed at high flood risk in Italy. Hydrol. Res. 2017, 48, 726–740. [Google Scholar] [CrossRef]
  4. Młyński, D.; Petroselli, A.; Wałęga, A. Flood frequency analysis by an event-based rainfall-runoff model in selected catchments of Southern Poland. Soil Water Res. 2018, 13, 170–176. [Google Scholar] [CrossRef] [Green Version]
  5. Zeng, Q.; Chen, H.; Xu, C.Y.; Jie, M.X.; Hou, Y.K. Feasibility and uncertainty of using conceptual rainfallrunoff models in design flood estimation. Hydrol. Res. 2016, 47, 701–717. [Google Scholar] [CrossRef] [Green Version]
  6. Talchabhadel, R.; Shakya, N.M. Rainfall Runoff Modelling for Flood. NEAJ Newsl. 2015, 1, 23–27. [Google Scholar]
  7. Sharghi, E.; Nourani, V.; Najafi, H.; Molajou, A. Emotional ANN (EANN) and Wavelet-ANN (WANN) Approaches for Markovian and Seasonal Based Modeling of Rainfall-Runoff Process. Water Resour. Manag. 2018, 32, 3441–3456. [Google Scholar] [CrossRef]
  8. Chau, K.W. Use of meta-heuristic techniques in rainfall-runoffmodelling. Water 2017, 9, 186. [Google Scholar] [CrossRef] [Green Version]
  9. Sitterson, J.; Knightes, C.; Parmar, R.; Wolfe, K.; Avant, B. An Overview of Rainfall-Runoff Model Types An Overview of Rainfall-Runoff Model Types. In Proceedings of the 9th International Congress on Environmental Modelling and Software [Internet], Fort Collins, CO, USA, 24–28 June 2018; pp. 1–29. Available online: https://scholarsarchive.byu.edu/iemssconference/2018/Stream-C/41/?utm_source=scholarsarchive.byu.edu%2Fiemssconference%2F2018%2FStream-C%2F41&utm_medium=PDF&utm_campaign=PDFCoverPages (accessed on 17 November 2022).
  10. Merabtene, T.; Kawamura, A.; Jinno, K.; Olsson, J. Risk assessment for optimal drought management of an integrated water resources system using a genetic algorithm. Hydrol. Process. 2003, 16, 2189–2208. [Google Scholar] [CrossRef]
  11. Sajikumar, N.; Remya, R.S. Impact of land cover and land use change on runoff characteristics. J. Environ. Manag. 2015, 161, 460–468. [Google Scholar] [CrossRef]
  12. Liu, Y.; Bralts, V.F.; Engel, B.A. Evaluating the effectiveness of management practices on hydrology and water quality at watershed scale with a rainfall-runoff model. Sci. Total Environ. 2015, 511, 298–308. [Google Scholar] [CrossRef]
  13. Rezaali, M.; Quilty, J.; Karimi, A. Probabilistic urban water demand forecasting using wavelet-based machine learning models. J. Hydrol. 2021, 600, 126358. [Google Scholar] [CrossRef]
  14. Fraga, I.; Cea, L.; Puertas, J. Effect of rainfall uncertainty on the performance of physically based rainfall–runoff models. Hydrol. Process. 2019, 33, 160–173. [Google Scholar] [CrossRef] [Green Version]
  15. Young, C.C.; Liu, W.C.; Wu, M.C. A physically based and machine learning hybrid approach for accurate rainfall-runoff modeling during extreme typhoon events. Appl. Soft Comput. J. 2017, 53, 205–216. [Google Scholar] [CrossRef]
  16. Tokar, A.S.; Johnson, P.A. Rainfall-Runoff Modeling Using Artificial Neural Networks. J. Hydrol. Eng. 1999, 4, 223–239. [Google Scholar] [CrossRef]
  17. Wilby, R.L.; Abrahart, R.J.; Dawson, C.W. Detection of conceptual model rainfall-runoff processes inside an artificial neural network. Hydrol. Sci. J. 2003, 48, 163–181. [Google Scholar] [CrossRef] [Green Version]
  18. Khan, M.T.; Shoaib, M.; Hammad, M.; Salahudin, H.; Ahmad, F.; Ahmad, S. Application of Machine Learning Techniques in Rainfall—Runoff Modelling of the Soan River Basin, Pakistan. Water 2021, 13, 3528. [Google Scholar] [CrossRef]
  19. Okkan, U.; Ersoy, Z.B.; Ali Kumanlioglu, A.; Fistikoglu, O. Embedding machine learning techniques into a conceptual model to improve monthly runoff simulation: A nested hybrid rainfall-runoff modeling. J. Hydrol. 2021, 598, 126433. [Google Scholar] [CrossRef]
  20. Poonia, V.; Tiwari, H.L. Rainfall-runoff modeling for the Hoshangabad Basin of Narmada River using artificial neural network. Arab. J. Geosci. 2020, 13, 944. [Google Scholar] [CrossRef]
  21. Gomes, E.P.; Blanco, C.J.C. Daily rainfall estimates considering seasonality from a MODWT-ANN hybrid model. J. Hydrol. Hydromech. 2021, 69, 13–28. [Google Scholar] [CrossRef]
  22. Nourani, V.; Tajbakhsh, A.D.; Molajou, A. Data mining based on wavelet and decision tree for rainfall-runoff simulation. Hydrol. Res. 2019, 50, 75–84. [Google Scholar] [CrossRef] [Green Version]
  23. Nourani, V.; Hosseini Baghanam, A.; Adamowski, J.; Kisi, O. Applications of hybrid wavelet-Artificial Intelligence models in hydrology: A review. J Hydrol. 2014, 514, 358–377. [Google Scholar] [CrossRef]
  24. Ouma, Y.O.; Cheruyot, R.; Wachera, A.N. Rainfall and runoff time-series trend analysis using LSTM recurrent neural network and wavelet neural network with satellite-based meteorological data: Case study of Nzoia hydrologic basin. Complex Intell. Syst. 2021, 8, 213–236. [Google Scholar] [CrossRef]
  25. Kavoosi, M.; Khozeymehnehad, H. Review and compare performance of 4 modeling methods LS-SVM, NN, GEP and ANFIS-PSO in Simulation of Rainfall–Runoff (Study Area: Halil River—Jiroft Dam). J. Irrig. Water Eng. 2021, 11, 96–110. [Google Scholar]
  26. Saifullah, M.W.M.; Hashim, S.; Shoaib, M.; Naseem, A.; Khan, M. Modelling of Rainfall-runoff process by GEP, RBF-SVM and M5 model tree in Jhelum River Basin, Pakistan. In Proceedings of the Internation Conference on Hydrology and Water Resources (ICHWR-2021), Lahore, Pakistan, 25 March 2021; pp. 105–112. [Google Scholar]
  27. Asadi, S.; Shahrabi, J.; Abbaszadeh, P.; Tabanmehr, S. A new hybrid artificial neural networks for rainfall-runoff process modeling. Neurocomputing 2013, 121, 470–480. [Google Scholar] [CrossRef]
  28. Okkan, U.; Serbes, Z.A. The combined use of wavelet transform and black box models in reservoir inflow modeling. J. Hydrol. Hydromech. 2013, 61, 112–119. [Google Scholar] [CrossRef] [Green Version]
  29. Okkan, U. Wavelet neural network model for reservoir inflow prediction. Sci. Iran. 2012, 19, 1445–1455. [Google Scholar] [CrossRef] [Green Version]
  30. Notarangelo, N.M.; Hirano, K.; Albano, R.; Sole, A. Transfer learning with convolutional neural netwoks for rainfall detection in single images. Water 2021, 13, 588. [Google Scholar] [CrossRef]
  31. Van, S.P.; Le, H.M.; Thanh, D.V.; Dang, T.D.; Loc, H.H.; Anh, D.T. Deep learning convolutional neural network in rainfall-runoff modelling. J. Hydroinf. 2020, 22, 541–561. [Google Scholar] [CrossRef] [Green Version]
  32. Xu, Y.; Liu, Y.; Jiang, Z.; Yang, X.; Wang, X.; Zhang, Y.; Qin, Y. Improved Convolutional Neural Network and its Application in Non-Periodical Runoff Prediction. Water Resour. Manag. 2022, 36, 6149–6168. [Google Scholar] [CrossRef]
  33. Poornima, S.; Pushpalatha, M. Prediction of rainfall using intensified LSTM based recurrent Neural Network with Weighted Linear Units. Atmosphere 2019, 10, 668. [Google Scholar] [CrossRef] [Green Version]
  34. Li, W.; Kiaghadi, A.; Dawson, C. High temporal resolution rainfall–runoff modeling using long-short-term-memory (LSTM) networks. Neural Comput. Appl. 2021, 33, 1261–1278. [Google Scholar] [CrossRef]
  35. Luo, Q.-R.; Xu, H.; Bai, L.-H. Prediction of significant wave height in hurricane area of the Atlantic Ocean using the Bi-LSTM with attention model. Ocean Eng. 2022, 266, 112747. [Google Scholar] [CrossRef]
  36. Shoaib, M.; Shamseldin, A.Y.; Khan, S.; Sultan, M.; Ahmad, F.; Sultan, T.; Dahri, Z.H.; Ali, I. Input Selection of Wavelet-Coupled Neural Network Models for Rainfall-Runoff Modelling. Water Resour. Manag. 2019, 33, 955–973. [Google Scholar] [CrossRef]
  37. Shoaib, M.; Shamseldin, A.Y.; Melville, B.W. A comparison between wavelet based static and dynamic neural network approaches for runoff prediction. J. Hydrol. 2016, 535, 211–225. [Google Scholar] [CrossRef]
  38. Shoaib, M.; Shamseldin, A.Y.; Melville, B.W. Comparative study of different wavelet based neural network models for rainfall-runoff modeling. J. Hydrol. 2014, 515, 47–58. [Google Scholar] [CrossRef]
  39. Shoaib, M.; Shamseldin, A.Y.; Melville, B.W. Hybrid Wavelet Neural Network Approach. In Network Modelling, Studies in Computational Intelligence; Springer International Publishing: Cham, Switzerland, 2016; pp. 127–143. [Google Scholar]
  40. Hammad, M.; Shoaib, M.; Salahudin, H.; Baig, M.A.I.; Khan, M.M.; Ullah, M.K. Rainfall forecasting in upper Indus basin using various artificial intelligence techniques. Stoch. Environ. Res. Risk Assess. 2021, 35, 2213–2235. [Google Scholar] [CrossRef]
  41. Mahmood, G.G.; Rashid, H.; Anwar, S.; Nasir, A. Evaluation of climate change impacts on rainfall patterns in pothohar region of pakistan. Water Conserv. Manag. 2019, 3, 1–6. [Google Scholar] [CrossRef]
  42. Tallat, Q.; Siddiqui, M.; Hashmi, H.N.; Ghumman, A.R. Runoff Modeling for a Watershed in Pothowar Region of Pakistan. Mehran Univ. Res. J. Eng. Technol. 2012, 31, 379–394. [Google Scholar]
  43. Siddiqui, Q.T.M.; Hashmi, H.N.; Mughal, H.R. Development of Small Dams in Pothowar Plateu of Punjab (Pakistan). In Proceedings of the ICID 21st International Congress on Irrigation and Drainage, Tehran, Iran, 15–23 October 2011; pp. 121–129. [Google Scholar]
  44. Breiman, L. Arcing the Edge; University of California: Berkeley, CA, USA, 1998. [Google Scholar]
  45. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  46. Raza, A.; Shoaib, M.; Khan, A.; Baig, F.; Faiz, M.A.; Khan, M.M. Application of non-conventional soft computing approaches for estimation of reference evapotranspiration in various climatic regions. Theor. Appl. Climatol. 2020, 139, 1459–1477. [Google Scholar] [CrossRef]
  47. Webros, P.J. Beyond Regression: New Tools for Prediction and Analysis in the Behavioral Sciences; Harvard University: Cambridge, MA, USA, 1974. [Google Scholar]
  48. Goodfellow, I.; Bengio, Y.; Courville, A. Deep Learning; MIT Press: Cambridge, MA, USA, 2016. [Google Scholar]
  49. Ferreira, C. Gene Expression Programming: A New Adaptive Algorithm for Solving Problems. Complex Syst. 2001, 13, 87–129. [Google Scholar]
  50. Percival, D.B.; Walden, A.T. The Maximal Overlap Discrete Wavelet Transform. In Wavelet Methods for Time Series Analysis; Cambridge University Press: Cambridge, UK, 2000; pp. 159–205. [Google Scholar]
  51. Quilty, J.; Adamowski, J. Addressing the incorrect usage of wavelet-based hydrological and water resources forecasting models for real-world applications with best practices and a new forecasting framework. J. Hydrol. 2018, 563, 336–353. [Google Scholar] [CrossRef]
  52. Walden, A.T. Wavelet Analysis of Discrete Time Series. Eur. Congr. Math. 2001, 2, 627–641. [Google Scholar]
  53. Ramana, R.V. Monthly Rainfall Prediction Using Wavelet Neural Network Analysis. Water Resour. Manag. 2013, 27, 3697–3711. [Google Scholar] [CrossRef] [Green Version]
  54. May, R.J.; Dandy, G.C.; Maier, H.R.; Nixon, J.B. Application of partial mutual information variable selection to ANN forecasting of water quality in water distribution systems. Environ. Model. Softw. 2008, 23, 1289–1299. [Google Scholar] [CrossRef]
  55. Ribeiro, G.H.T.; De Neto, P.S.G.M.; Cavalcanti, G.D.C.; Tsang, I.R. Lag selection for time series forecasting using Particle Swarm Optimization. In Proceedings of the International Joint Conference on Neural Networks, San Jose, CA, USA, 31 July–5 August 2011; pp. 2437–2444. [Google Scholar]
  56. Maheswaran, R.; Khosa, R. Comparative study of different wavelets for hydrologic forecasting. Comput. Geosci. 2012, 46, 284–295. [Google Scholar] [CrossRef]
  57. Nash, J.E.; Sutcliffe, J.V. River Flow Forecasting Through Conceptual Models Part I—A Disclission of Principles*. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  58. Hussain, F.; Wu, R.; Nabi, G.; Hussain, B. A Study on Rainfall, Runoff and Sediment Yield Relations in Small Watersheds of Pothwar Region Pakistan. In Proceedings of the 2017 Conference on Sustainable Development for Environmental Resources, Taiwan, 24–26 November 2017. [Google Scholar]
  59. Ghani, M.W.; Arshad, M.; Shabbir, A.; Shakoor, A.; Mehmood, N.; Ahmad, I. Investigation of potential water harvesting sites at potohar using modeling approach. Pakistan J. Agric. Sci. 2013, 50, 723–729. [Google Scholar]
  60. Jhong, Y.D.; Chen, C.S.; Lin, H.P.; Chen, S.T. Physical hybrid neural network model to forecast typhoon floods. Water 2018, 10, 632. [Google Scholar] [CrossRef] [Green Version]
  61. Song, Y.Y.; Lu, Y. Decision tree methods: Applications for classification and prediction. Shanghai Arch. Psychiatry 2015, 27, 130–135. [Google Scholar]
  62. Lin, N.; Noe, D.; He, X. Tree-based methods and their applications. In Springer Handbook of Engineering Statistics, 49th ed.; Pham, H., Ed.; Springer: London, UK, 2006; pp. 551–570. [Google Scholar]
  63. Danial Mohammadzadeh, S.; Kazemi, S.F.; Mosavi, A.; Nasseralshariati, E.; Tah, J.H.M. Prediction of compression index of fine-grained soils using a gene expression programming model. Infrastructures 2019, 4, 26. [Google Scholar] [CrossRef] [Green Version]
  64. Tu, J.V. Advantages and disadvantages of using artificial neural networks versus logistic regression for predicting medical outcomes. J. Clin. Epidemiol. 1996, 49, 1225–1231. [Google Scholar] [CrossRef]
  65. Adamowski, J.; Sun, K. Development of a coupled wavelet transform and neural network method for flow forecasting of non-perennial rivers in semi-arid watersheds. J. Hydrol. 2010, 390, 85–91. [Google Scholar] [CrossRef]
Figure 1. Location map of river basins within the Pothohar Region.
Figure 1. Location map of river basins within the Pothohar Region.
Atmosphere 14 00452 g001
Figure 7. Graphical illustration of detail components (at level 1, 2, and 3) and approximation component at level 3 for rainfall timeseries of the Soan-Haro River basin.
Figure 7. Graphical illustration of detail components (at level 1, 2, and 3) and approximation component at level 3 for rainfall timeseries of the Soan-Haro River basin.
Atmosphere 14 00452 g007
Figure 8. Training NSE score of stand-alone machine learning models in Bunha-Kahan River basin.
Figure 8. Training NSE score of stand-alone machine learning models in Bunha-Kahan River basin.
Atmosphere 14 00452 g008
Figure 9. Testing NSE score of stand-alone machine learning models in Bunha-Kahan River basin.
Figure 9. Testing NSE score of stand-alone machine learning models in Bunha-Kahan River basin.
Atmosphere 14 00452 g009
Figure 10. Training NSE score of hybrid machine learning models in Bunha-Kahan River basin.
Figure 10. Training NSE score of hybrid machine learning models in Bunha-Kahan River basin.
Atmosphere 14 00452 g010
Figure 11. Testing NSE score of hybrid machine learning models in Bunha-Kahan River basin.
Figure 11. Testing NSE score of hybrid machine learning models in Bunha-Kahan River basin.
Atmosphere 14 00452 g011
Figure 12. Comparative analysis of best performing models of stand-alone ML techniques in Bunha-Kahan River basin.
Figure 12. Comparative analysis of best performing models of stand-alone ML techniques in Bunha-Kahan River basin.
Atmosphere 14 00452 g012
Figure 13. Comparative analysis of best performing models of stand-alone and MODWT-hybrid techniques in Bunha-Kahan River basin.
Figure 13. Comparative analysis of best performing models of stand-alone and MODWT-hybrid techniques in Bunha-Kahan River basin.
Atmosphere 14 00452 g013
Figure 14. Comparative analysis of best models in Bunha-Kahan River basin.
Figure 14. Comparative analysis of best models in Bunha-Kahan River basin.
Atmosphere 14 00452 g014
Figure 15. Regression plot (MODWT-DTF-Lo3) in Bunha-Kahan River basin.
Figure 15. Regression plot (MODWT-DTF-Lo3) in Bunha-Kahan River basin.
Atmosphere 14 00452 g015
Figure 16. Rainfall-Runoff modeling of Bunha-Kahan River basin (MODWT-DTF-Lo3).
Figure 16. Rainfall-Runoff modeling of Bunha-Kahan River basin (MODWT-DTF-Lo3).
Atmosphere 14 00452 g016
Figure 17. NSE score of stand-alone machine learning models in training in the Soan-Haro River basin.
Figure 17. NSE score of stand-alone machine learning models in training in the Soan-Haro River basin.
Atmosphere 14 00452 g017
Figure 18. NSE score of stand-alone machine learning models in testing in the Soan-Haro River basin.
Figure 18. NSE score of stand-alone machine learning models in testing in the Soan-Haro River basin.
Atmosphere 14 00452 g018
Figure 19. NSE score of hybrid machine learning models in training in the Soan-Haro River basin.
Figure 19. NSE score of hybrid machine learning models in training in the Soan-Haro River basin.
Atmosphere 14 00452 g019
Figure 20. NSE score of hybrid machine learning models in testing in the Soan-Haro River basin.
Figure 20. NSE score of hybrid machine learning models in testing in the Soan-Haro River basin.
Atmosphere 14 00452 g020
Figure 21. Comparative analysis of best performing models of stand-alone ML techniques in the Soan-Haro River basin.
Figure 21. Comparative analysis of best performing models of stand-alone ML techniques in the Soan-Haro River basin.
Atmosphere 14 00452 g021
Figure 22. Comparative analysis of best performing models of stand-alone and MODWT-hybrid techniques in the Soan-Haro River basin.
Figure 22. Comparative analysis of best performing models of stand-alone and MODWT-hybrid techniques in the Soan-Haro River basin.
Atmosphere 14 00452 g022
Figure 23. Comparative analysis of best models in the Soan-Haro River basin.
Figure 23. Comparative analysis of best models in the Soan-Haro River basin.
Atmosphere 14 00452 g023
Figure 24. Regression plot (MODWT-DTF-Lo4) in the Soan-Haro River basin.
Figure 24. Regression plot (MODWT-DTF-Lo4) in the Soan-Haro River basin.
Atmosphere 14 00452 g024
Figure 25. Rainfall-Runoff Modeling of the Soan-Haro River basin (MODWT-DTF-Lo4).
Figure 25. Rainfall-Runoff Modeling of the Soan-Haro River basin (MODWT-DTF-Lo4).
Atmosphere 14 00452 g025
Table 1. Statistical features of rainfall and runoff datasets of the two selected homogeneus areas.
Table 1. Statistical features of rainfall and runoff datasets of the two selected homogeneus areas.
Min.Max.MeanModeS.D.
Soan and Haro River basins
Training (1999 to 2011)
Rainfall (mm)068.702.8605.06
Runoff (m3/s)2.161407.6925.494.9074.00
Testing (2012 to 2016)
Rainfall (mm)072.843.3705.71
Runoff (m3/s)1.961559.7433.425.8972.81
Bunha and Kahan River basins
Training (1999 to 2011)
Rainfall (mm)065.432.3706.18
Gauge Reading (m)790.40870.10856.958575.41
Testing (2012 to 2016)
Rainfall (mm)0140.722.6107.10
Gauge Reading (m)857.00868.40857.218570.61
Table 2. Summary of attributes related to gauge stations showed in Figure 1.
Table 2. Summary of attributes related to gauge stations showed in Figure 1.
Serial No.Station NameBasin(s)Latitude (°N)Longitude (°E)Altitude (m)
Raingauge Stations
1MurreeSoan & Haro33.907073.39432025
2RawalpindiSoan33.565173.0169540
3Kotli SattianSoan33.808273.52551352
4ChakwalSoan & Bunha32.932872.8630522
5Fateh JungSoan & Haro33.563572.6375514
6TalagangSoan32.917272.4081457
7Gujjar KhanSoan33.261673.3058458
8PendighebSoan33.245272.266310
9TaxilaSoan & Haro33.746372.8397549
10Khanpur DamSoan33.801872.9305545
11NathiagalliSoan34.068073.39222410
12JabbriHaro33.904573.1733923
13Wahh CanttHaro33.784372.7388471
14HattarHaro33.852172.8501513
15AbbotabadHaro34.149573.19951220
16NaraKuhan & Bunha33.983373.21661627
17Mulhal MughlanKuhan & Bunha73.149733.0027523
18Khokhar BalaKuhan & Bunha72.814632.7716853
19Mohra SharifKuhan & Bunha73.68333.6671500
20WaghKuhan & Bunha32.734373.3966272
21SohawaKuhan & Bunha33.117173.4149440
22DinaKuhan & Bunha33.030673.6069281
Runoff Gauging Stations
1Makhad (gauging station)Soan33.028171.7393252
2Bhambar (gauging station)Bunha32.769473.3634263
Table 3. Parametric summary of adopted techniques.
Table 3. Parametric summary of adopted techniques.
ML TechniqueParameters
SDTMin. rows in a nodeMin. size node to splitMax. tree levelsCross-validation treesSmooth min. spikes-
51010103-
TBMax. trees in seriesMin. size node to splitTrimming factorMin. trees in seriesSmooth min. spikesRandom percent
400100.0110520
DTFTrees in forestMin. size node to splitMax. tree levelsRandom predictor control-
200250square root of total predictors-
MLPHidden layersMin. neuronsMax. neuronsCross-validationMax. iterationAlg.
1220V-fold10,000LV- Marquardt
GEPPopulation sizeMax. initial pop.Max. generationGene per chromosomeGene head lengthFitness threshold
5010,0002000481
Table 4. Summary of best performing model’s comparison in Bunha-Kahan River Basin.
Table 4. Summary of best performing model’s comparison in Bunha-Kahan River Basin.
ML TechniqueLag OrderWavelet FilterLevelTesting
NSERMSE (mm)R2
MODWT-SDTLo4db430.32339.720.31
MODWT-DTFLo3db430.86220.450.92
MODWT-TBLo1db430.46417.220.47
MODWT-MLPLo1db430.29425.740.31
MODWT-GEPLo2db430.62332.120.58
Table 5. Summary of best performing model’s comparison in the Soan-Haro River basin.
Table 5. Summary of best performing model’s comparison in the Soan-Haro River basin.
ML TechniqueLag Order (Lo)Wavelet FilterLevelTesting
NSERMSER2
MODWT-SDTLo3db430.3163.030.31
MODWT-DTFLo4db430.8821.720.91
MODWT-TBLo2db430.5442.340.60
MODWT-MLPLo1db430.2560.570.25
MODWT-GEPLo1db430.3070.160.60
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Khan, M.T.; Shoaib, M.; Albano, R.; Inam, M.A.; Salahudin, H.; Hammad, M.; Ahmad, S.; Ali, M.U.; Hashim, S.; Ullah, M.K. Intercomparison and Assessment of Stand-Alone and Wavelet-Coupled Machine Learning Models for Simulating Rainfall-Runoff Process in Four Basins of Pothohar Region, Pakistan. Atmosphere 2023, 14, 452. https://doi.org/10.3390/atmos14030452

AMA Style

Khan MT, Shoaib M, Albano R, Inam MA, Salahudin H, Hammad M, Ahmad S, Ali MU, Hashim S, Ullah MK. Intercomparison and Assessment of Stand-Alone and Wavelet-Coupled Machine Learning Models for Simulating Rainfall-Runoff Process in Four Basins of Pothohar Region, Pakistan. Atmosphere. 2023; 14(3):452. https://doi.org/10.3390/atmos14030452

Chicago/Turabian Style

Khan, Muhammad Tariq, Muhammad Shoaib, Raffaele Albano, Muhammad Azhar Inam, Hamza Salahudin, Muhammad Hammad, Shakil Ahmad, Muhammad Usman Ali, Sarfraz Hashim, and Muhammad Kaleem Ullah. 2023. "Intercomparison and Assessment of Stand-Alone and Wavelet-Coupled Machine Learning Models for Simulating Rainfall-Runoff Process in Four Basins of Pothohar Region, Pakistan" Atmosphere 14, no. 3: 452. https://doi.org/10.3390/atmos14030452

APA Style

Khan, M. T., Shoaib, M., Albano, R., Inam, M. A., Salahudin, H., Hammad, M., Ahmad, S., Ali, M. U., Hashim, S., & Ullah, M. K. (2023). Intercomparison and Assessment of Stand-Alone and Wavelet-Coupled Machine Learning Models for Simulating Rainfall-Runoff Process in Four Basins of Pothohar Region, Pakistan. Atmosphere, 14(3), 452. https://doi.org/10.3390/atmos14030452

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop