Next Article in Journal
Transport and Fate of Nitrate in the Streambed of a Low-Gradient Stream
Previous Article in Journal
Assessment of Alternative Agricultural Land Use Options for Extending the Availability of the Ogallala Aquifer in the Northern High Plains of Texas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Short-Term Water Demand Forecasting Model Combining Variational Mode Decomposition and Extreme Learning Machine

Department of Constructional and Environmental Engineering, Kyungpook National University, Sangju 37224, Korea
*
Author to whom correspondence should be addressed.
Hydrology 2018, 5(4), 54; https://doi.org/10.3390/hydrology5040054
Submission received: 14 August 2018 / Revised: 19 September 2018 / Accepted: 25 September 2018 / Published: 27 September 2018

Abstract

:
Accurate water demand forecasting is essential to operate urban water supply facilities efficiently and ensure water demands for urban residents. This study proposes an extreme learning machine (ELM) coupled with variational mode decomposition (VMD) for short-term water demand forecasting in six cities (Anseong-si, Hwaseong-si, Pyeongtaek-si, Osan-si, Suwon-si, and Yongin-si), South Korea. The performance of VMD-ELM model is investigated based on performance indices and graphical analysis and compared with that of artificial neural network (ANN), ELM, and VMD-ANN models. VMD is employed for multi-scale time series decomposition and ANN and ELM models are used for sub-time series forecasting. As a result, ELM model outperforms ANN model. VMD-ANN and VMD-ELM models outperform ANN and ELM models, and the VMD-ELM model produces the best performance among all the models. The results obtained from this study reveal that the coupling of VMD and ELM can be an effective forecasting tool for short-term water demands with strong nonlinearity and non-stationarity and contribute to operating urban water supply facilities efficiently.

1. Introduction

Many urban areas around the world are confronted with stresses associated with water supply due to natural and social factors including economic growth, overpopulation, and climate change [1,2,3,4,5]. To solve the problems, accurate water demand forecasting and the expansion and efficient operation of water supply and distribution facilities are essential. Especially, reliable and accurate water demand forecasting is essential to develop reliable water supply expansion strategies at the lowest cost. However, water demand forecasting is still a challenging task due to the availability of data, various influencing factors, various forecasting periods, and the nonlinearity and non-stationarity of data [1,6].
Water demand forecasting can be classified as short-term, medium-term, and long-term forecasting in terms of forecast horizon although there is no general rule for the classification. According to Tiwari and Adamowski [7], hourly forecasting (up to 48-h lead time), daily forecasting (up to 14-day lead time), and weekly forecasting (up to 26-week lead time) are classified as short-term forecasting. Medium-term forecasting includes monthly forecasting with up to 24-month lead time, whereas annual and decadal forecastings are considered as long-term forecasting. Decision problems in tactical planning level including revenue forecast and investment planning are treated based on medium-term forecasting, whereas long-term forecasting is used for decision problems such as capacity expansion in strategic planning level. This study is focused on short-term forecasting based on daily water demand time series. Short-term forecasting can be utilized for dealing with decision problems related to the operational management and optimization of water supply system in operational planning level [6].
Water demand forecasting methods can be classified into unit water demand analysis, qualitative methods, univariate time series methods, exponential smoothing models, moving average models, time series regression models, stochastic process models, scenario-based approaches, decision support system, machine learning models, and composite models [6]. Especially, machine learning models and composite models have been recently applied in order to improve the accuracy of short-term water demand forecasting [7,8,9,10]. For example, Adamowski et al. [8] proposed an artificial neural network (ANN) coupled with discrete wavelet transform (DWT) for urban water demand forecasting. They found that the DWT-ANN model forecasted water demands more accurately than other models including multiple linear regression (MLR), multiple nonlinear regression (MNLR), autoregressive integrated moving average (ARIMA), and ANN models. Tiwari and Adamowski [7] proposed bootstrap-based ANN (BNN) and wavelet-bootstrap-neural network (WBNN) models for short-term urban water demand forecasting. They revealed that the WBNN model reduced the forecast uncertainty and produced more accurate and reliable water demands and confidence intervals. Multi-scale relevance vector regression (MSRVR) was proposed by Bai et al. [9] for urban water demand forecasting. They used stationary wavelet transform (SWT) in order to decompose water supply time series into multi-scale components and then trained relevance vector regression (RVR) using wavelet coefficient for each scale. They found that the MSRVR model was able to be an effective tool for improving daily urban water demand forecasting. Brentan et al. [10] proposed a hybrid regression combining support vector regression (SVR) and adaptive Fourier series (AFS) in order to forecast short-term water demands. They revealed that the combination of forecasting (SVR) and post-processing (AFS) models were effective to near real-time water demand forecasting.
Furthermore, many other studies on water demand forecasting have been conducted recently. The combination of seasonal autoregressive integrated moving average (SARIMA) models and data assimilation for short-term water demand forecasting was presented by Arandia et al. [11]. They found that for offline mode, the models with weekly seasonality outperformed ones with daily periodicity, and for online mode, the forecasting accuracy was able to be improved significantly by data assimilation. Gagliardi et al. [12] proposed a Markov chain-based short-term water demand forecasting model. They showed that homogeneous Markov chain model was more efficient than non-homogeneous Markov chain model and produced better forecasting accuracy compared with ANN and naïve models. A moving window-based model for short-term water demand forecasting was presented by Pacchin et al. [13]. They demonstrated that the proposed model was able to provide good accuracy over the entire forecasting time. To assess the predictive uncertainty of hourly water demand forecasting, Alvisi and Franchini [14] applied model conditional processor (MCP) which can combine water demand forecasts obtained from two or more forecasting models. They revealed that MCP was able to provide better forecasting accuracy compared with single forecasting models (ANN and periodic pattern-producing models). An overview of short-term water demand forecasting methods including autoregressive (AR), moving average (MA), autoregressive moving average (ARMA), autoregressive moving average exogenous (ARMAX), feed-forward back-propagation neural network (FFBP-NN), and hybrid models (combination of ARMA and FFBP-NN) was presented by Anele et al. [15]. From comparative analysis utilizing a common dataset, they concluded that ARMA, ARMAX, and hybrid models may be the best candidates for assessing the predictive distribution of future water demands. Anele et al. [16] applied an MCP approach for water demand forecasting and compared it with ARMA, FFBP-NN, and hybrid models (coupling of ARMA and FFBP-NN). They revealed that MCP marginally outperformed the compared models and may be an effective tool for real-time short-term water demand forecasting.
This study proposes an extreme learning machine (ELM) coupled with variational mode decomposition (VMD) in order to enhance the accuracy of short-term water demand forecasting. VMD, which is a relatively novel time series decomposition technique in hydrological fields, is employed in this study for decomposing an original water demand time series into sub-time series. VMD has been successfully applied for the analysis of hydrological variables and exhibited better performance compared with conventional methods including DWT and ensemble empirical mode decomposition (EEMD) [17,18]. DWT has drawbacks that mother wavelets, decomposition level, and edge effect should be considered for time series decomposition, and empirical mode decomposition (EMD) has disadvantages including mode mixing, stopping criteria, and end-point effect. EEMD was developed for mitigating the problems of EMD (especially the mode mixing) but the problems were not completely resolved. On the other hand, VMD decomposes an original signal into multiple modes and then updates them based on Wiener filtering. By utilizing Wiener filtering, VMD can be more robust to sampling and noise and yield narrow-banded modes (see Dragomiretskiy and Zosso [19] and Theorem 1 of Polyak and Pearlman [20]). This helps to alleviate the effect of mode mixing and extract the time-frequency features accurately. In addition, VMD can be implemented faster and more efficiently compared with EMD and EEMD since VMD is a non-recursive algorithm based on augmented Lagrangian method and alternate direction method of multipliers (ADMM) (see Algorithm 1 of Dragomiretskiy and Zosso [19]).
Meanwhile, ELM is selected as a forecasting model for the decomposed water demands due to its ease of modeling and excellent performance [18,21,22]. Conventional gradient descent-based ANN models require many iterative learning steps for obtaining the optimal learning performance. This increases the computational time and the possibility of being trapped in the local minima. Recurrent neural network (RNN) is known as a powerful model for time series modeling. However, in practice, RNN has the disadvantage that it is difficult to be trained properly due to the vanishing gradient and exploding gradient problems [23]. Deep learning models require many training datasets and very long computational time since they have many weights and biases that should be trained due to their structural complexity. Further, optimization algorithms need to be used for determining the optimal learning parameters since deep learning models have many learning parameters that should be selected in advance. Moreover, deep learning models require many computational resources and long computational time, even using graphics processing unit (GPU)-based deep learning models. In practice, one may get just as good performance with non-deep learning models using far less computational time and resources although this is dependent on the given problems and datasets [24]. On the other hand, for the ELM model, the input weights and biases are assigned randomly, and the smallest norm least-squares solution for the output weights is calculated analytically using Moore-Penrose generalized inverse (see Definitions 2.1–2.2 and Theorem 2.1 of Huang et al. [25]). This can increase the learning speed extremely and provide better generalization compared with conventional ANN models. Thus, in this study, the performance of VMD-ELM model for short-term water demand forecasting in six cities (Anseong-si, Hwaseong-si, Pyeongtaek-si, Osan-si, Suwon-si, and Yongin-si), South Korea is investigated and compared with that of ANN, ELM, and VMD-ANN models, based on performance indices and graphical analysis.

2. Materials and Methods

2.1. Data Used

Figure 1 shows the locations of water purification plant and water supply areas. The Suji water treatment plant, which is located in the north of study area, supplies water to six cities, Anseong-si, Hwaseong-si, Pyeongtaek-si, Osan-si, Suwon-si, and Yongin-si, which are located in the southern region of Gyeonggi-do, South Korea. The plant has the capacity of 916,000 m3/day and has been operated by the K-water (http://www.kwater.or.kr, accessed on 12 August 2018). Table 1 shows the area, population, and water demand of study area. The total area is 2460.51 km2 with a total population of 3,840,907. Especially, Suwon-si is a city with the largest population among local governments, South Korea and has the largest population density in the study area. Yongin-si is the fastest growing city in South Korea, which has the second largest population in the study area. Suwon-si and Yongin-si have the largest water demand among the cities in the study area. Domestic water demand is higher in Osan-si, Suwon-si, and Yongin-si, whereas agricultural water demand is higher in Anseong-si, Hwaseong-si, and Pyeongtaek-si.
In this study, daily water demand data between 2008 and 2017 provided by K-water were used for developing water demand forecasting models. For training the models efficiently [26], the water demand data were scaled into the range of [0, 1]. There is no clear criterion to divide the entire data into training and testing data. However, the length of the training data is typically set at 70–90% of the total data length [27,28]. Thus, in this study, the scaled data were then partitioned into training (2008–2014, data length = 2557 (70%)) and testing datasets (2015–2017, data length = 1096 (30%)).

2.2. Variational Mode Decomposition (VMD)

VMD, which is developed by Dragomiretskiy and Zosso [19], is a non-recursive and adaptive time-frequency analysis method. The VMD decomposes an original time series f into K sub-time series called intrinsic mode functions (IMFs). The IMFs can be obtained by updating Equations (1)–(3) until convergence [19].
u ^ k n + 1 ( ω ) = f ^ ( ω ) i < k u ^ i n + 1 ( ω ) i > k u ^ i n ( ω ) + λ ^ n ( ω ) / 2 1 + 2 α ( ω ω k n ) 2  
ω k n + 1 = 0 ω | u ^ k n + 1 ( ω ) | 2 d ω 0 | u ^ k n + 1 ( ω ) | 2 d ω  
λ ^ n + 1 ( ω ) = λ ^ n ( ω ) + τ [ f ^ ( ω ) k u ^ k n + 1 ( ω ) ]  
where u ^ k n = the kth mode in nth iteration, f ^ = the Fourier transform of f, ω k n = the kth center frequency in nth iteration, λ ^ n = the Lagrange multiplier in nth iteration, α = the quadratic penalty factor, τ = the time step of dual ascent. For the theoretical details of VMD, one can refer to Dragomiretskiy and Zosso [19].

2.3. Artificial Neural Network (ANN)

ANN is an artificial intelligence (AI) model for solving various problems including classification and regression [29]. Multilayer perceptron (MLP) is a kind of feedforward ANN organized into layers with multiple neurons [30], which has been widely applied in hydrological fields. Figure 2 represents the general structure of three-layer MLP. The neurons of the input layer receive the input vectors (lagged water demands), and output values (forecasted water demands) are yielded from the neuron of the output layer. The hidden layer, which is the middle layer of MLP, connects the neurons of the input and output layers. The MLP can be represented as Equation (4) [31]:
y ^ j = i = 1 L β i h ( w i x j + b i ) + β 0 , j = 1 , 2 , , N  
where x j and y ^ j = the input and output vectors, w i and β i = the weights for the hidden and output layers, b i and β 0 = the biases for the hidden and output layers, L = the number of the neurons of hidden layer, N = the length of time series data, and h = the transfer function. The determination of the parameters (weights and biases) of MLP, which is called model learning, is conducted through learning algorithms such as backpropagation (BP). The detailed theoretical background of ANN can be found in Alpaydin [32].

2.4. Extreme Learning Machine (ELM)

ELM is a least square-based single-hidden layer feed-forward neural networks (SLFNs) [33]. In terms of training speed and generalization capability, ELM is known to be superior to conventional ANN [34]. The SLFN with the weights and biases generated randomly can be expressed as Equation (5) and written compactly as Equation (6).
i = 1 L β i h ( w i x j + b i ) = t j , i = 1 , 2 , , L , j = 1 , 2 , , N  
H β = T  
where H = h ( w i x j + b i ) = the matrix comprised of the outputs from hidden neurons, h = the transfer function, w i = [ w i 1 , w i 2 , , w i n ] T = the ith weight vector between input and hidden neurons, x j = [ x j 1 , x j 2 , , x j n ] T = the jth input vectors, b i = the ith bias for hidden neurons, L and N = the number of hidden neurons and the data length, respectively, β = [ β 1 T , β 2 T , , β L T ] L × m T = the matrix comprised of weights between hidden and output neurons, β i = [ β i 1 , β i 2 , , β i m ] T = the ith weight vectors between hidden and output neurons, and T = [ t 1 T , t 2 T , , t N T ] N × m T = the matrix comprised of target vectors t j m .
Unlike conventional ANN, the β is analytically determined by the method of least squares as Equation (7).
β ^ = H T  
where H = P H T is the Moore-Penrose generalized inverse of H, and P = ( H T H ) 1 is the inverse of the covariance matrix of H. Figure 3 represents the general structure of ELM model used in this study. For detailed theoretical background on the ELM, one can refer to Huang et al. [33].

2.5. VMD-Based Water Demand Forecasting

VMD-based machine learning models (VMD-ANN and VMD-ELM) for water demand forecasting are to couple VMD and single machine learning models (ANN and ELM), respectively. VMD was adopted as a time series decomposition method for decomposing water demand time series into IMFs, and ANN and ELM models were used for model learning and forecasting for each IMF. Figure 4 shows the flowchart of VMD-based water demand forecasting consisting of the following four steps:
Step 1.
Decomposition: Water demand time series is decomposed into IMFs using VMD.
Step 2.
Model learning: ANN and ELM models are learned for each IMF.
Step 3.
IMF forecasts: ANN and ELM models produce forecasted values for each IMF.
Step 4.
Final forecasts: Summing the forecasted IMFs produces the final water demand forecasts.

2.6. Performance Evaluation Indices

In this study, multiple performance indices were used for evaluating the model performances. Mathematical expressions and ranges for the performance indices are presented in Table 2. Since there is no universal performance index, it is desirable to choose the indices that are appropriate for a particular application. Furthermore, it is common to use multiple performance indices since there are pros and cons in each index [35,36]. The performance indices used in this study are categorized as absolute, relative, and dimensionless errors. The absolute errors include the mean absolute error (MAE), the root mean squared error (RMSE), and the fourth root mean quadrupled error (R4MS4E). RMSE and R4MS4E are sensitive to forecasting errors for peak and high data values, and thus can be good performance indices for high data values. In contrast, MAE is evaluated based on all deviations without being weighted to lower and higher data values. MAE, RMSE, and R4MS4E have the advantage that they can represent the size of a typical error effectively since they are evaluated in the same unit as the original data. For a perfect model, MAE, RMSE, and R4MS4E would be zero [36]. The mean absolute relative error (MARE) and median absolute percentage error (MdAPE) belong to relative errors. MARE is a good index for low data values since it is more affected to forecasting errors for low data values. MdAPE is similar to MARE but it has the advantage that it is less sensitive to skewed error distributions and outliers. For a perfect model, MARE and MdAPE would be zero [36]. The dimensionless errors include the modified versions of the coefficient of efficiency (CE) and the index of agreement (IOA). The CE and IOA have the drawback that they can overestimate the model performance for peak data values. The modified CE (MCE) and modified IOA (MIOA) for j = 1 can reduce the overestimation for peak data values and represent better overall evaluation, whereas larger values of j can be used for the performance evaluation of high data values [35]. For a perfect model, MCE and MIOA would be one [35,36].

3. Results and Discussion

3.1. Development of Single and VMD-Based Forecasting Models

To develop VMD-based water demand models, water demand time series should first be decomposed using the VMD. The decomposition results depend on the quadratic penalty factor (α) and the number of IMFs (K). In this study, the values of α and K were determined based on the following trial-and-error method:
Step 1.
For K = {2, 3, ......, 20} and α = {5, 10, 20, 50, 70, 100, 150, 200, 500, 1000, 2000}, the corresponding IMFs are generated from water demand time series.
Step 2.
For each set of K and α values, the IMFs are summed to reconstruct the water demand time series.
Step 3.
Correlation coefficients (r) between original and reconstructed water demand time series are estimated.
Step 4.
The sets of K and α values corresponding to r 1 are selected.
Step 5.
The optimal values of K and α are selected based on the performances of VMD-based forecasting models.
Figure 5 represents the correlation coefficients between the water demand time series and the reconstructed time series for different K and α values. Based on Figure 5, the sets of K and α values corresponding to r 1 were selected and then K = 4 and α = 5 producing the best performance of VMD-based forecasting models were selected as the optimal values. Using the K and α values, the original water demand time series was decomposed into four sub-time series (IMF 1-IMF 4) as seen in Figure 6.
For developing forecasting models using ANN and ELM, the potential influencing input variables should be selected in advance. Various lag times for daily water demand time series were considered in this study. The optimal lag time was determined by autocorrelation function (ACF), partial autocorrelation function (PACF), and average mutual information (AMI) [37]. The structure of input and target variables for the ANN and ELM models is as follows:
Q t + i = f ( Q t 5 , Q t 4 , Q t 3 , Q t 2 , Q t 1 , Q t )  
where Q t j = the input water demand time series lagged by j days (j = 0, 1, ......, 5), Q t + i = the target water demand time series for lead time i day (i = 1, 2, ......, 7), and f = the ANN and ELM models. In the same manner, the structure of input and target variables for VMD-ANN and VMD-ELM models is determined as follows:
I M F 1 t + i = f ( I M F 1 t 9 , I M F 1 t 8 , , I M F 1 t 1 , I M F 1 t )  
I M F 2 t + i = f ( I M F 2 t 1 , I M F 2 t )  
I M F 3 t + i = f ( I M F 3 t 1 , I M F 3 t )  
I M F 2 t + i = f ( I M F 2 t 1 , I M F 2 t )  
Figure 7 shows correlation coefficients between IMFs obtained from the entire water demand time series and ones from the partial water demand time series with different lengths. From Figure 7, it was seen that time series decomposition by VMD was dependent on the length of water demand time series. When the length of the partial time series was more than 500, the decomposition result was almost the same as the result for the entire time series. Thus, for real-time forecasting, VMD should be executed whenever new observations are acquired, and at least 500 previous time series data should be used for time series decomposition.
In ANN and ELM modeling, selecting the optimal number of hidden neurons is a critical step since it affects the model performance. A trial-and-error approach [37] was utilized for determining the optimal number in this study. In most hydrological fields, the logistic sigmoid function has been widely applied to ANN modeling in order to compute the output of each neuron [26]. The function is also applied to ANN and ELM modeling in this study. Furthermore, the ANN model was learned by BP algorithm and the default values of learning parameters (learning and momentum rates) presented by Zell et al. [38] were used. When considering both processing speed and accuracy, it is more efficient to use the appropriate small values of learning parameters than to determine them using a trial-and-error approach [39]. More information on this can be found in Rumelhart et al. [40], Pao [41], and Seo et al. [18].

3.2. Performance Evaluation

In this study, absolute error indices (MAE, RMSE, and R4MS4E), relative error indices (MARE and MdAPE), and dimensionless error indices (MCE and MIOA) were employed in order to evaluate and compare model performances for testing data. MCE and MIOA with j = 1 (MCE1 and MIOA1) were used to evaluate overall performances, and MCE and MIOA with j = 2 and 3 (MCE2, MCE3, MIOA2, and MIOA3) were employed to evaluate performances for high water demands. The results of performance evaluation for testing data are summarized in Table 3.
For single forecasting models (ANN and ELM), the ELM model yielded lower MAE and higher MCE1 and MIOA1 than ANN model. ANN model produced negative MCE1. These indicated that ELM model represented better overall performance compared with ANN model, and ANN model performed worse than a “no knowledge” model. ELM model yielded lower values of MARE and MdAPE compared with ANN model. ELM model produced better performance than ANN model for low water demands. Furthermore, the ELM model yielded lower values of RMSE and R4MS4E and higher values of MCE2, MIOA2, MCE3, and MIOA3 than ANN model. ANN model produced negative values of MCE2 and MCE3. These indicated that ELM model represented better performance for higher water demands compared with ANN model, and ANN model produced worse performance than a “no knowledge” model. Consequently, it was found from the results of single forecasting models that ELM model performed better over the entire range of testing data compared with ANN model.
For VMD-based forecasting models (VMD-ANN and VMD-ELM), VMD-ELM model yielded better performance indices than VMD-ANN model. For example, the MAE values of VMD-ELM model were lower than those of VMD-ANN model, and the MCE1 and MIOA1 values of VMD-ELM model were higher than those of VMD-ANN model. These represented that the overall performance of VMD-ELM model was better compared with VMD-ANN model. Furthermore, the values of MARE and MdAPE for VMD-ELM model were smaller than those for VMD-ANN model, and VMD-ELM produced lower values of RMSE and R4MS4E and higher values of MCE2, MIOA2, MCE3, and MIOA3 compared with VMD-ANN models. These results represented that VMD-ELM model performed better than VMD-ANN model for low and high water demands. Consequently, from the results of VMD-based forecasting models, it was found that VMD-ELM model outperformed VMD-ANN model over the entire range of testing data.
These results are due to the theoretical characteristics of ANN and ELM models. For the ANN model, gradient descent-based learning algorithms such as BP algorithm adjust all the parameters iteratively and require many iterative steps for obtaining satisfactory performance. Thus, the learning process is very slow, and it is difficult for the ANN model to reach the minimum training error due to infinite training iteration and local minima. In addition, the ANN model may be over-trained by gradient descent-based learning and produce worse generalization performance. Gradient descent-based learning algorithms only ty to reach the smallest training error, and do not consider the magnitude of the weights. These characteristics make it difficult for the ANN model to achieve the best generalization performance [25]. On the other hand, the ELM model is considered as a linear system after the input weights and the biases of the hidden layer are assigned randomly. Namely, the input weights and the biases of the hidden layer do not need to be tuned and any iterative learning process is also not required. The output weights of the ELM model can be determined analytically by finding the smallest norm least-squares solution for the linear system using the Moore-Penrose generalized inverse. Thus, the ELM model can produce the minimum training error, the smallest norm of weights, and the best generalization performance [25].
From the comparison of single and VMD-based forecasting models, VMD-ANN and VMD-ELM models represented better performance compared with ANN and ELM models. For instance, VMD-ANN and VMD-ELM models yielded lower MAE and higher MCE1 and MIOA1 values than ANN and ELM models, which represented better overall performance for VMD-ANN and VMD-ELM models. Furthermore, VMD-ANN and VMD-ELM models yielded lower MARE and MdAPE values than ANN and ELM models. For VMD-ANN and VMD-ELM models, the values of RMSE and R4MS4E were lower and the values of MCE2, MIOA2, MCE3, and MIOA3 were higher compared with ANN and ELM models. These results represented that VMD-ANN and VMD-ELM models produced better performance than ANN and ELM models for low and high water demands. Therefore, from the comparison of single and VMD-based forecasting models, it was found that VMD-ANN and VMD-ELM models outperformed ANN and ELM models over the entire range of testing data, VMD-ELM model produced the best performance, and the performances of ANN and ELM models were able to be improved by VMD.
In general, water demand time series is a combination of different time-frequency features and has strong nonlinearity and non-stationarity. In this study, VMD decomposes an original water demand time series into multiple modes (low- and high-frequency modes), and then the modes are updated by Wiener filtering [42]. Since the modes represent simpler patterns than the original time series, the ANN and ELM models can be learned more effectively using the modes as training datasets than using the original time series. Especially, since low-frequency mode (IMF 1), which represents the trend component of water demand time series, can be modeled very accurately by the single forecasting models, the performance of water demand forecasting can be improved significantly by VMD. Thus, for water demand time series with strong nonlinearity and non-stationarity, VMD-based forecasting models can produce better performance than single forecasting models.
Figure 8 shows scatter diagrams for forecasted and observed water demands. It was observed from Figure 8 that the scatter points of VMD-ANN and VMD-ELM models were closer to 1:1 slope line than those of ANN and ELM models. Based on scatter diagrams, VMD-ELM model produced the most accurate result among all other models. Although VMD-ANN and VMD-ELM models represented similar dispersion around 1:1 slope line, VMD-ANN model yielded underestimated values for high water demands. ELM model represented better accuracy than ANN model. Especially, the ANN model produced greatly underestimated result for high water demands. Although VMD was able to improve the underestimation of ANN model for high water demands, the underestimation was not resolved completely even in VMD-ANN model. Consequently, VMD-ELM model represented better accuracy compared with all other models, in terms of under- and overestimation and dispersion around 1:1 slope line. It was also found from these results that the accuracy of ANN and ELM models was able to be improved significantly by VMD.
Figure 9 shows Taylor diagram for single and VMD-based forecasting models. The diagram provides a graphical summary for the agreement of observed and forecasted patterns, based on correlation coefficient (CC), standard deviation (SD), and centered root-mean-square difference (CRMSD) [43]. In Figure 9, the gray contour lines and black arc represent the values of CRMSD and SD for observed pattern, respectively. The closer a model is to the black hollow circle marker, the closer the forecasted pattern is to the observed pattern. From Figure 9, it was observed that VMD-ANN and VMD-ELM models yielded higher CC and lower CRMSD than ANN and ELM models, and ELM and VMD-ELM models produced SD values closer to the observed pattern compared with ANN and VMD-ANN models. Especially, it was shown that VMD-ELM models were located the closest to the black hollow circle marker among all other models. Consequently, based on Taylor diagram, it can be said that VMD-ANN and VMD-ELM models outperformed ANN and ELM models, and VMD was able to improve the accuracy of ANN and ELM models.
Figure 10 and Figure 11 show residual boxplots and time series plots for single and VMD-based forecasting models, respectively. It was observed from Figure 10 that the ranges (maximum—minimum) and interquartile ranges (third quartile—first quartile) of residuals for VMD-ANN and VMD-ELM models were smaller than those for the ANN and ELM models, and the median values of residuals for the ANN and ELM models were biased in the negative direction from zero compared with VMD-ANN and VMD-ELM models, respectively. Furthermore, it was observed from Figure 11 that VMD-ANN and VMD-ELM models forecasted the observed pattern more accurately than ANN and ELM models. Especially, compared with VMD-ANN and VMD-ELM models, the ANN and ELM models showed greater over- and underestimation for low and high water demands, respectively. These results indicated that VMD was able to reduce the residuals of ANN and ELM models significantly, improve under- and overestimation represented in ANN and ELM models, and thus enhance the accuracy of ANN and ELM models.
From the results obtained from this study, VMD-based forecasting models produced better performances than single forecasting models in daily water demand forecasting. Similar results have also been reported from previous studies dealing with DWT-based hybrid modeling for water demand [7,8,9] and VMD-based hybrid modeling for river stage [17] and rainfall-runoff [18]. VMD decomposes water demand time series, which is a mixture of complex and irregular components, into multiple sub-time series (IMFs). Machine learning-based forecasting models are then constructed for each sub-time series. Since the sub-time series represent simpler temporal patterns than the original time series, the machine learning models can be trained more effectively, and thus water demand time series with strong nonlinearity and non-stationarity can be forecasted more accurately. Furthermore, by applying Wiener filtering, VMD can be more robust to noise affecting the decomposition accuracy of water demand time series negatively and yield narrow-based modes. Due to these characteristics, the VMD not only alleviates the mode mixing but also provides accurate time-frequency features. In addition, since VMD is a non-recursive signal decomposition algorithm, error propagation which inevitably occurs in iterative algorithms does not occur. These characteristics enable the VMD to produce more accurate time-frequency features compared with other signal decomposition methods such as EMD [42,44].
Meanwhile, VMD-based water demand forecasting models can be a more reliable forecasting tool by combining with resampling techniques including simple bootstrap, block bootstrap, Gaussian process regression bootstrap, Bayesian bootstrap, and maximum entropy bootstrap. Hybrid modeling similar to this concept was carried out successfully by Tiwari and Adamowski [7]. They developed a hybrid model (WBNN) coupling simple bootstrap, DWT, and ANN in order to forecast short-term urban water demand and concluded that WBNN model was able to decrease the forecasting uncertainty and thus provide more accurate and reliable water demand forecasts. Therefore, the combined model of resampling techniques, VMD, and machine learning models can be suggested as a future study for improving short-term water demand forecasting.

4. Conclusions

This study investigates the performances of VMD-based water demand forecasting models. VMD and two machine learning models, the ANN and ELM, were employed for the decomposition of water demand time series and the forecasting of sub-time series, respectively. The performances of single forecasting models (ANN and ELM) and VMD-based forecasting models (VMD-ANN and VMD-ELM) are evaluated based on performance indices and graphical analysis. For single forecasting models, the ELM model outperforms the ANN model. VMD-ANN and VMD-ELM models perform better compared with ANN and ELM models, and the VMD-ELM model achieves the best performance among all the models. From the results of this study, it is revealed that VMD can reduce forecasting error and improve model performance in machine learning-based water demand forecasting. Therefore, VMD-based machine learning models can be effective forecasting tools for short-term water demand with strong nonlinearity and non-stationarity and contribute to operate urban water supply facilities efficiently. This study is limited to short-term water demand forecasting for a specific area. Further studies need to be performed with regard to water demand forecasting for other areas which have demographic, economic, and weather conditions similar to the area where forecasting models are developed.

Author Contributions

Conceptualization, Y.S.; Methodology, Y.S.; Formal Analysis, Y.S.; Investigation, Y.S.; Data Curation, Y.S.; Writing-Original Draft Preparation, Y.S.; Writing-Review & Editing, Y.S., S.K. and Y.C.; Supervision, Y.C.; and Project Administration, Y.C.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. House-Peters, L.A.; Chang, H. Urban water demand modeling: Review of concepts, methods, and organizing principles. Water Resour. Res. 2011, 47, W05401. [Google Scholar] [CrossRef]
  2. Schuetze, T.; Santiago-Fandiño, V. Quantitative assessment of water use efficiency in urban and domestic buildings. Water 2013, 5, 1172–1193. [Google Scholar] [CrossRef]
  3. Hao, L.; Sun, G.; Liu, Y.; Qian, H. Integrated modeling of water supply and demand under management options and climate change scenarios in Chifeng city, China. J. Am. Water Resour. Assoc. 2015, 51, 655–671. [Google Scholar] [CrossRef]
  4. Arsiso, B.K.; Tsidu, G.M.; Stoffberg, G.H.; Tadesse, T. Climate change and population growth impacts on surface water supply and demand of Addis Ababa, Ethiopia. Clim. Risk Manag. 2017, 18, 21–33. [Google Scholar] [CrossRef]
  5. Lee, J.S.; Kim, J.W. Assessing strategies for urban climate change adaptation: The case of six metropolitan cities in South Korea. Sustainability 2018, 10, 2065. [Google Scholar] [CrossRef]
  6. Donkor, E.A.; Mazzuchi, T.A.; Soyer, R.; Roberson, J.A. Urban water demand forecasting: Review of methods and models. J. Water Res. Plan. Manag. 2014, 140, 146–159. [Google Scholar] [CrossRef]
  7. Tiwari, M.K.; Adamowski, J. Urban water demand forecasting and uncertainty assessment using ensemble wavelet-bootstrap-neural network models. Water Resour. Res. 2013, 49, 6486–6507. [Google Scholar] [CrossRef] [Green Version]
  8. Adamowski, J.; Chan, H.F.; Prasher, S.O.; Ozga-Zielinski, B.; Sliusarieva, A. Comparison of multiple linear and nonlinear regression, autoregressive integrated moving average, artificial neural network, and wavelet artificial neural network methods for urban water demand forecasting in Montreal, Canada. Water Resour. Res. 2012, 48, W01528. [Google Scholar] [CrossRef]
  9. Bai, Y.; Wang, P.; Li, C.; Xie, J.; Wang, Y. A multi-scale relevance vector regression approach for daily urban water demand forecasting. J. Hydrol. 2014, 517, 236–245. [Google Scholar] [CrossRef]
  10. Brentan, B.M.; Luvizotto, E., Jr.; Herrera, M.; Izquierdo, J.; Pérez-García, R. Hybrid regression model for near real-time urban water demand forecasting. J. Comput. Appl. Math. 2017, 309, 532–541. [Google Scholar] [CrossRef]
  11. Arandia, E.; Ba, A.; Eck, B.; McKenna, S. Tailoring seasonal time series models to forecast shot-term water demand. J. Water Res. Plan. Manag. 2016, 142, 1–10. [Google Scholar] [CrossRef]
  12. Gagliardi, F.; Alvisi, S.; Kapelan, Z.; Franchini, M. A probabilistic short-term water demand forecasting model based on the Markov chain. Water 2017, 9, 507. [Google Scholar] [CrossRef]
  13. Pacchin, E.; Alvisi, S.; Franchini, M. A short-term water demand forecasting model using a moving window on previously observed data. Water 2017, 9, 172. [Google Scholar] [CrossRef]
  14. Alvisi, S.; Franchini, M. Assessment of predictive uncertainty within the framework of water demand forecasting using the Model Conditional Processor (MCP). Urban Water J. 2017, 14, 1–10. [Google Scholar] [CrossRef]
  15. Anele, A.O.; Hamam, Y.; Abu-Mahfouz, A.M.; Todini, E. Overview, comparative assessment and recommendations of forecasting models for short-term water demand prediction. Water 2017, 9, 887. [Google Scholar] [CrossRef]
  16. Anele, A.O.; Todini, E.; Hamam, Y.; Abu-Mahfouz, A.M. Predictive uncertainty estimation in water demand forecasting using the model conditional processor. Water 2018, 10, 475. [Google Scholar] [CrossRef]
  17. Seo, Y.; Kim, S.; Singh, V.P. Comparison of different heuristic and decomposition techniques for river stage modeling. Environ. Monit. Assess. 2018, 190, 392. [Google Scholar] [CrossRef] [PubMed]
  18. Seo, Y.; Kim, S.; Singh, V.P. Machine learning models coupled with variational mode decomposition: A new approach for modeling daily rainfall-runoff. Atmosphere 2018, 9, 251. [Google Scholar] [CrossRef]
  19. Dragomiretskiy, K.; Zosso, D. Variational mode decomposition. IEEE Trans. Signal Process. 2014, 62, 531–544. [Google Scholar] [CrossRef]
  20. Polyak, N.; Pearlman, W.A. Stationarity of the Gabor basis and derivation of Janssen’s formula. In Proceedings of the IEEE-SP International Symposium on Time-Frequency and Time-Scale Analysis, Victoria, BC, Canada, 4–6 October 1992. [Google Scholar] [CrossRef]
  21. Deo, R.C.; Şahin, M. An extreme learning machine model for the simulation of monthly mean streamflow water level in eastern Queensland. Environ. Monit. Assess. 2016, 188, 90. [Google Scholar] [CrossRef] [PubMed]
  22. Yaseen, Z.M.; Jaafar, O.; Deo, R.C.; Kisi, O.; Adamowski, J.; Quilty, J.; El-Shafie, A. Stream-flow forecasting using extreme learning machines: A case study in a semi-arid region in Iraq. J. Hydrol. 2016, 542, 603–614. [Google Scholar] [CrossRef]
  23. Bengio, Y.; Simard, P.; Frasconi, P. Learning ling-term dependencies with gradient descent is difficult. IEEE Trans. Neural Netw. 1994, 5, 157–166. [Google Scholar] [CrossRef] [PubMed]
  24. Brink, H.; Richards, J.W.; Fetherolf, M. Real-World Machine Learning; Manning Publications Co.: Shelter Island, NY, USA, 2017; ISBN 978-1617291920. [Google Scholar]
  25. Huang, G.B.; Zhu, Q.Y.; Siew, C.K. Extreme learning machine: A new learning scheme of feedforward neural networks. In Proceedings of the 2004 IEEE International Joint Conference on Neural Networks, Budapest, Hungary, 25–29 July 2004. [Google Scholar] [CrossRef]
  26. Dawson, C.W.; Wilby, R.L. Hydrological modelling using artificial neural networks. Prog. Phys. Geogr. 2001, 25, 80–108. [Google Scholar] [CrossRef]
  27. Partal, T.; Cigizoglu, H.K. Estimation and forecasting of daily suspended sediment data using wavelet-neural networks. J. Hydrol. 2008, 358, 317–331. [Google Scholar] [CrossRef]
  28. Mohanta, A.; Patra, K.C.; Sahoo, B.B. Anticipate Manning’s coefficient in meandering compound channels. Hydrology 2018, 5, 47. [Google Scholar] [CrossRef]
  29. ASCE Task Committee. Artificial neural networks in hydrology. I: Preliminary concepts. J. Hydrol. Eng. 2000, 5, 115–123. [Google Scholar] [CrossRef]
  30. Zahmatkesh, Z.; Goharian, E. Comparing machine learning and decision making approaches to forecast long lead monthly rainfall: The city of Vancouver, Canada. Hydrology 2018, 5, 10. [Google Scholar] [CrossRef]
  31. Lima, A.R.; Cannon, A.J.; Hsieh, W.W. Forecasting daily streamflow using online sequential extreme learning machines. J. Hydrol. 2016, 537, 431–443. [Google Scholar] [CrossRef]
  32. Alpaydin, E. Introduction to Machine Learning, 2nd ed.; The MIT Press: Cambridge, MA, USA, 2010; ISBN 978-0262028189. [Google Scholar]
  33. Huang, G.B.; Zhu, Q.Y.; Siew, C.K. Extreme learning machine: Theory and application. Neurocomputing 2006, 70, 489–501. [Google Scholar] [CrossRef]
  34. Yu, L.; Dai, W.; Tang, L. A novel decomposition ensemble model with extended extreme learning machine for crude oil price forecasting. Eng. Appl. Artif. Intell. 2016, 47, 110–121. [Google Scholar] [CrossRef]
  35. Krause, P.; Boyle, D.P.; Bäse, F. Comparison of different efficiency criteria for hydrological model assessment. Adv. Geosci. 2005, 5, 89–97. [Google Scholar] [CrossRef]
  36. Dawson, C.W.; Abrahart, R.J.; See, L.M. HydroTest: A web-based toolbox of evaluation metrics for the standardized assessment of hydrological forecasts. Environ. Model. Softw. 2007, 22, 1034–1052. [Google Scholar] [CrossRef]
  37. Seo, Y.; Kim, S.; Kisi, O.; Singh, V.P. Daily water level forecasting using wavelet decomposition and artificial neural intelligence techniques. J. Hydrol. 2015, 520, 224–243. [Google Scholar] [CrossRef]
  38. Zell, A.; Mamier, G.; Mache, M.V.N.; Hübner, R.; Dörin, S.; Hermann, K.U. SNNS Stuttgart Neural Network Simulator v. 4.2, User Manual. University of Stuttgart/University of Tübingen. 1998. Available online: http://www.ra.cs.uni-tuebingen.de/downloads/SNNS/SNNSv4.2.Manual.pdf (accessed on 30 June 2018).
  39. Dai, H.; MacBeth, C. Effects of learning parameters on learning procedure and performance of a BPNN. Neural Netw. 1997, 10, 1505–1521. [Google Scholar] [CrossRef]
  40. Rumelhart, D.; Hinton, G.E.; Williams, R.J. Learning representations by back-propagating errors. Nature 1986, 323, 533–536. [Google Scholar] [CrossRef]
  41. Pao, Y.H. Adaptive Pattern Recognition and Neural Networks; Addison-Wesley Publishing Company, Inc.: New York, NY, USA, 1988; ISBN 978-0201125849. [Google Scholar]
  42. Shi, P.; Yang, W. Precise feature extraction from wind turbine condition monitoring signals by using optimised variational mode decomposition. IET Renew. Power Gener. 2017, 11, 245–252. [Google Scholar] [CrossRef]
  43. Taylor, K.E. Summarizing multiple aspects of model performance in a single diagram. J. Geophys. Res. Atmos. 2001, 106, 7183–7192. [Google Scholar] [CrossRef] [Green Version]
  44. Sun, G.; Chen, T.; Wei, Z.; Sun, Y.; Zang, H.; Chen, S. A carbon price forecasting model based on variational mode decomposition and spiking neural networks. Energies 2016, 9, 54. [Google Scholar] [CrossRef]
Figure 1. Water supply areas and location of water purification plant.
Figure 1. Water supply areas and location of water purification plant.
Hydrology 05 00054 g001
Figure 2. General structure of multilayer perceptron (MLP)-type artificial neural network (ANN).
Figure 2. General structure of multilayer perceptron (MLP)-type artificial neural network (ANN).
Hydrology 05 00054 g002
Figure 3. General structure of extreme learning machine (ELM).
Figure 3. General structure of extreme learning machine (ELM).
Hydrology 05 00054 g003
Figure 4. Flowchart for variational mode decomposition (VMD)-based water demand forecasting.
Figure 4. Flowchart for variational mode decomposition (VMD)-based water demand forecasting.
Hydrology 05 00054 g004
Figure 5. Correlation relationship between original water demand time series and reconstructed time series for different K and α values.
Figure 5. Correlation relationship between original water demand time series and reconstructed time series for different K and α values.
Hydrology 05 00054 g005
Figure 6. Water demand time series and IMFs.
Figure 6. Water demand time series and IMFs.
Hydrology 05 00054 g006aHydrology 05 00054 g006b
Figure 7. Correlation coefficients between IMFs obtained from the entire water demand time series and ones from the partial water demand time series with different lengths.
Figure 7. Correlation coefficients between IMFs obtained from the entire water demand time series and ones from the partial water demand time series with different lengths.
Hydrology 05 00054 g007
Figure 8. Scatter diagrams for single and VMD-based forecasting models.
Figure 8. Scatter diagrams for single and VMD-based forecasting models.
Hydrology 05 00054 g008
Figure 9. Taylor diagram for single and VMD-based forecasting models.
Figure 9. Taylor diagram for single and VMD-based forecasting models.
Hydrology 05 00054 g009
Figure 10. Residual boxplots for single and VMD-based forecasting models
Figure 10. Residual boxplots for single and VMD-based forecasting models
Hydrology 05 00054 g010
Figure 11. Time series plots for single and VMD-based forecasting models.
Figure 11. Time series plots for single and VMD-based forecasting models.
Hydrology 05 00054 g011aHydrology 05 00054 g011b
Table 1. Area, population, and water demand of study area.
Table 1. Area, population, and water demand of study area.
CitiesArea (km2)Population (People)Population Density (People/km2)Water Demand (103 m3/year)
TotalDomesticIndustrialAgricultural
Anseong-si553.39182,294329.4245,00041,35317,873185,774
Hwaseong-si693.92729,9391051.9381,37596,07530,260255,040
Pyeongtaek-si458.08489,0811067.7384,01380,82546,221256,967
Osan-si42.73218,6355116.737,40324,16064746769
Suwon-si121.051,203,2859940.4582,384520,38514,33947,660
Yongin-si591.341,017,6731721.0705,685418,1712733284,781
Source: Water Resources Management Information System (http://www.wamis.go.kr (accessed on 12 August 2018)).
Table 2. Mathematical expressions and ranges for performance indices.
Table 2. Mathematical expressions and ranges for performance indices.
Performance IndicesEquationsRanges
Absolute errorsMAE M A E = 1 N i = 1 N | Q i Q ^ i | [ 0 , ]
RMSE R M S E = i = 1 N ( Q i Q ^ i ) 2 N [ 0 , ]
R4MS4E R 4 M S 4 E = 1 N i = 1 N ( Q i Q ^ i ) 4 4 [ 0 , ]
Relative errorsMARE M A R E = 1 N i = 1 N | Q i Q ^ i | Q i [ 0 , ]
MdAPE M d A P E = M e d i a n ( | Q i Q ^ i Q i | × 100 ) [ 0 , ]
Dimensionless errorsMCE M C E j = 1 i = 1 N | Q i Q ^ i | j i = 1 N | Q i Q ¯ | j with j [ , 1 ]
MIOA M I O A j = 1 i = 1 N | Q i Q ^ i | j i = 1 N ( | Q ^ i Q ¯ | + | Q i Q ¯ | ) j with j [ 0 , 1 ]
Q ^ i : forecasted water demand, Q i : observed water demand; Q ¯ : mean of observed water demand; and N: data length.
Table 3. Performance evaluation for testing data.
Table 3. Performance evaluation for testing data.
ModelsLead Times (Days)MAE (m3/day)RMSE (m3/day)R4MS4E (m3/day)MAREMdAPEMCE1MIOA1MCE2MIOA2MCE3MIOA3
ANN118,69424,36634,7000.0282.330−0.1410.576−0.2840.797−0.4220.901
221,32227,65839,6360.0322.592−0.3620.516−0.8310.728−1.5330.835
322,71429,46842,8120.0352.851−0.4600.491−1.0860.692−2.0880.792
421,92429,20343,9840.0332.616−0.4310.498−1.1060.688−2.2960.775
522,04429,44745,1900.0342.606−0.4320.496−1.1260.683−2.4210.761
623,44430,66945,4190.0362.824−0.5140.478−1.2820.666−2.6100.752
725,25232,16645,9930.0393.245−0.6390.452−1.5510.644−3.0650.739
ELM114,62019,92830,4460.0221.6220.4210.6930.6480.8940.7830.961
216,70323,18437,6960.0251.8490.3380.6430.5230.8460.6250.920
317,89024,93240,9000.0272.0030.2910.6180.4490.8220.5270.898
419,80926,84142,1270.0302.2850.2160.5580.3610.7670.4520.857
521,00628,23543,9370.0312.5100.1690.5330.2940.7400.3720.832
619,58527,32243,6230.0292.1610.2260.5940.3390.7930.3980.875
720,37527,51643,0730.0312.4160.1950.5610.3300.7700.4120.858
VMD-ANN16417871913,4550.0090.7170.7460.8640.9330.9810.9820.997
2885211,89917,3490.0131.0040.6490.8080.8740.9610.9570.993
312,16616,15423,8470.0181.4320.5180.7280.7690.9230.8910.978
414,29118,21025,5810.0211.8240.4340.6820.7060.9020.8570.972
513,70118,04626,7350.0201.6550.4580.6860.7110.8990.8490.967
616,22521,09430,6590.0241.9870.3590.6310.6060.8600.7650.946
717,05823,01335,3250.0261.9450.3260.6090.5310.8270.6640.916
VMD-ELM13637492775190.0050.4030.8560.9270.9780.9940.9970.999
26327864312,9990.0100.6960.7490.8680.9340.9810.9830.997
310,15613,82421,8050.0151.1610.5980.7840.8300.9490.9250.987
410,06813,69121,9450.0151.1590.6010.7840.8340.9490.9260.987
511,88816,01025,4080.0181.3930.5300.7450.7730.9300.8840.979
612,83617,51327,6040.0191.4600.4930.7280.7280.9170.8470.972
714,40520,07731,5530.0221.6340.4310.6960.6430.8900.7660.956

Share and Cite

MDPI and ACS Style

Seo, Y.; Kwon, S.; Choi, Y. Short-Term Water Demand Forecasting Model Combining Variational Mode Decomposition and Extreme Learning Machine. Hydrology 2018, 5, 54. https://doi.org/10.3390/hydrology5040054

AMA Style

Seo Y, Kwon S, Choi Y. Short-Term Water Demand Forecasting Model Combining Variational Mode Decomposition and Extreme Learning Machine. Hydrology. 2018; 5(4):54. https://doi.org/10.3390/hydrology5040054

Chicago/Turabian Style

Seo, Youngmin, Soonmyeong Kwon, and Yunyoung Choi. 2018. "Short-Term Water Demand Forecasting Model Combining Variational Mode Decomposition and Extreme Learning Machine" Hydrology 5, no. 4: 54. https://doi.org/10.3390/hydrology5040054

APA Style

Seo, Y., Kwon, S., & Choi, Y. (2018). Short-Term Water Demand Forecasting Model Combining Variational Mode Decomposition and Extreme Learning Machine. Hydrology, 5(4), 54. https://doi.org/10.3390/hydrology5040054

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