Next Article in Journal
Effects of Groin Type and Bed Properties on Flow in Groin Fields: Comparison of Fixed- and Mobile-Bed Experiments
Next Article in Special Issue
The University of West Florida Campus Ecosystem Study: Spatial and Temporal Variation in Water Quality at Thompson Bayou
Previous Article in Journal
Impact of Hydrological Conditions on the Isotopic Composition of the Sava River in the Area of the Zagreb Aquifer
Previous Article in Special Issue
Annual Evaluation of 17 Oestrogenic Endocrine Disruptors and Hazard Indexes in the Douro River Estuary—The Atlantic Discharge of the Highest-Flow River of Southwestern Europe
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Evaluation of À Trous-Based Record Extension Techniques for Water Quality Record Extension

1
Irrigation & Hydraulics Department, Faculty of Engineering, Ain-Shams University, Cairo 11566, Egypt
2
Civil Engineering Department, Faculty of Engineering (Mataria), Helwan University, Cairo 11795, Egypt
3
Construction and Building Engineering Department, College of Engineering and Technology, Cairo Branch, Arab Academy for Science, Technology and Maritime Transport, Cairo 11799, Egypt
4
National Water Research Center, Cairo 12622, Egypt
*
Author to whom correspondence should be addressed.
Water 2022, 14(14), 2264; https://doi.org/10.3390/w14142264
Submission received: 31 May 2022 / Revised: 15 July 2022 / Accepted: 17 July 2022 / Published: 20 July 2022
(This article belongs to the Special Issue Water Quality Modeling and Monitoring)

Abstract

:
Hydrological data in general and water quality (WQ) data in particular frequently suffer from missing records and/or short-gauged monitoring/sampling sites. Many statistical regression techniques are employed to substitute missing values or to extend records at short-gauged sites, such as the Kendall-Theil robust line (KTRL), its modified version (KTRL2), ordinary least squares regression (OLS), four MOVE techniques, and the robust line of organic correlation (RLOC). In this study, in aspiring to achieve better accuracy and precision, the À Trous-Haar wavelet transform (WT) was adopted as a data denoising preprocessing step prior to applying record extension techniques. An empirical study was performed using real WQ data, from the National WQ monitoring network in the Nile Delta in Egypt, to evaluate the performance of these eight record-extension techniques with and without the WT data preprocessing step. Evaluations included the accuracy and precision of the techniques when used for the restoration of WQ missing values and for the extension of the WQ short-gauged variable. The results indicated that for the restoration of missing values, the KTRL and WT-KTRL outperformed other techniques. However, for the extension of short-gauged variables, WT-KTRL2, WT-MOVE3, and WT-MOVE4 techniques showed more accurate and precise results compared with both other techniques and their counterparts without the WT.

1. Introduction

The availability of representative and accurate hydrological and water quality (WQ) data is a crucial part of long-term water resource management [1,2,3]. Water resources data in general and WQ data in particular suffer from missing records and/or short-gauged monitoring/sampling sites [4,5]. Many statistical regression approaches have been applied for the restoration of missing records and/or the extension of short-gauged water resources records. One of the widely used regression techniques for both the restoration of missing hydrological and WQ records and the extension of records at short-gauged sites is the simple linear regression technique (ordinary least-squares regression-OLS) [6]. Two of the OLS main assumptions are: the explanatory and response variables are linearly related (highly correlated); and the data used are representative. However, the OLS suffers from two major flaws:
  • It is sensitive to the existence of outliers (extreme/unusual records), as outliers’ existence seriously affects the estimation of the OLS intercept and slope parameters [7,8];
  • It generates extended records with an underestimated variance [4,9,10,11,12]. Producing extended records with underestimated variance results from a bias in the estimation of extreme values, which as a result produces a bias in the estimation of exceedance and non-exceedance probabilities [4,5]. For WQ management generally and particularly for WQ assessment, high percentiles and extreme values are critical for evaluating whether WQ is within accepted limits or standards [13].
Water quality data have unique features, like nonnegative values, positive skewness, nonnormal distribution, presence of censored values (e.g., below a detection limit), presence of outlier/extreme values, seasonal patterns, and autocorrelation. Two of the more common features are the positive skewness and existence of outliers; owing to these two features, WQ data often have a form approaching the lognormal distribution [5,14,15]. Given that these WQ data mainly characterize outliers’ existence, deviation from normal distribution, and the presence of censored values, a robust or nonparametric regression technique may be more appropriate.
Several robust regression techniques (e.g., least median of squares; least absolute deviations; Winsorized regression; and trimmed least square estimation) and nonparametric regression techniques (e.g., monotonic regression and Kendall–Theil robust line (KTRL) have been designed as analogs to OLS with the advantage of being insensitive to outliers’ existence [7,16]. Nevitt and Tam [7] studied the performance of different robust and nonparametric approaches compared to the performance of OLS in cases of outliers’ existence and deviation from normality. The results indicated that the KTRL is the most convincing analog to OLS regression, with the advantage of being insensitive to outliers’ existence and/or deviation from normality. The KTRL is widely employed not only for both record extension and the restoration of missing values/records but also for the assessment of monotonic trend in water resources data (e.g., [8,17,18,19,20]). Although the KTRL overcomes one of the main OLS flaws by being insensitive to outliers’ existence, it is similar to OLS in producing extended records with an underestimated variance.
MOVE1 & MOVE2 proposed by [9] and MOVE3 & MOVE4 proposed by [10] are examples of Maintenance of Variance Extension (MOVE) systems that give extended records with unbiased variance. Several studies have used MOVE approaches to extend stream-flow records (e.g., [21,22,23]), for missing precipitation values restoration (e.g., [24,25]), and for WQ record extension (e.g., [5]). Although the MOVE techniques overcome one of the major flaws of OLS, they also are sensitive to outliers’ existence.
Recently, three regression techniques to overcome not just one but the two major flaws of the OLS were proposed: Khalil et al. [4] provided a modified version of the KTRL (KTRL2), Khalil and Adamowski [4] proposed the Robust Line of Organic Correlation, and Khalil et al. [26] proposed a modified version of the MOVE-1 that involves the L-moments in the estimate of the model slope (LM-R). Several studies have been carried out to assess the performance of these three techniques (KTRL2, RLOC, and LM-R) using Monte Carlo and empirical experiments [4,5,27,28]. These studies showed that in the case where outliers exist, these three newly developed techniques (KTRL2, RLOC, and LM-R) outperform the four MOVE techniques in producing extended records with unbiased variance. However, for the restoration of missing records, the KTRL outperforms other techniques in the existence of outliers.
More recently, Nalley et al. [2] employed the À Trous-Haar wavelet transform (WT) as a data preprocessing step before applying record-extension techniques for streamflow record extension. In Nalley et al. [2], the WT was applied to the predictor (x) and response (y) variables to create an approximation and detailed components for each x and y; record extension techniques were then applied separately for the extension of each component, and the results to estimate the streamflow extended records were finally summed up. Nalley et al. [2] compared the performance of the OLS, KTRL, MOVE techniques, KTRL2, and RLOC, with and without the WT data preprocessing step, for the extension of streamflow records using streamflow data at 67 paired sites from Canada’s Reference Hydrometric Basin Network. The main results showed consistent improvements in the WT-KTRL2, WT-RLOC, WT-MOVE1, and WT-MOVE2 techniques’ precision and accuracy when compared to their traditional counterparts (without WT), especially for the extended records statistical parameters. However, as a restoration for missing streamflow values, the WT-based techniques showed inconsistent improvements.
In this study, the use of the “À Trous-Haar” WT as a data preprocessing step for WQ record extension and missing WQ values restoration was evaluated. In contrast to the approach of [2], in this study, the WT was applied only to the predictor, and only the predictor approximation component was used in the record extension technique to directly estimate the missing WQ values. In addition, unlike streamflow, WQ data have much shorter temporal coverage and suffer from irregular sampling. To put it another way, streamflow indicates a long time series, whereas WQ indicates a short time series. The major objective of this research is to explore the À Trous-Haar WT-record extension techniques for WQ record extension and missing value restoration. An empirical experiment that utilizes real WQ data obtained from the Egyptian National WQ Monitoring Network was applied. In this experiment, eight record-extension techniques, the OLS, KTRL, KTRL2, RLOC, and four MOVE techniques were examined with and without the À Trous-Haar WT data preprocessing step.

2. Materials and Methods

It is assumed that two highly correlated WQ variables, a WQ variable x that has n 1 + n 2 measured records, and a WQ variable y that has only n 1 concurrent measured records, are illustrated as the following forms:
  x 1 ; x 2 ; x 3 ; ; x n 1 ;   x n 1 + 1 ;   x n 1 + 2 ;   ; ... ; x n 1 + n 2
  y 1 ; y 2 ; y 3 ; ; y n 1 ;                                                                                                                                            
Records of the short-gauged WQ variable y can be estimated/extended for the period n 1 + 1 through n 1 + n 2 using statistical regression techniques. Another common case, where both WQ variables x and y have similar sizes (e.g., n 1 + n 2 ) of concurrent measured records, but the WQ variable y has some missing values (m), is illustrated as the following forms:
  x 1 ; x 2 ; x 3 ; ; ; x n 1 ;   x n 1 + 1 ;   x n 1 + 2 ;   ; ; ; ; ... ; x n 1 + n 2   y 1 ; y 2 ; y 3 ; . ; m ; ; m ;   x n 1 + 1 ;   x n 1 + 2 ;   m ; ; m ; ; ... ; x n 1 + n 2
Missing values (m) in the WQ variable y can also be substituted using regression techniques. Simple regression techniques that involve only one predictor (independent variable x) have the following form [29]:
y ^ i = c + s   x i + ε
where c is the regression intercept, s is the slope, x i is the observed x values, and y ^ i is the estimated y values for the period through i = n 1 + 1 ,…   n 1 + n 2 (record extension), or, for missing values (m) in WQ variable y, and ε is a random error (ε mean value is zero) [29]. The main difference between these simple regression techniques is the way in which c and s are estimated. The following subsections briefly illustrate the eight record-extension techniques examined in this study: the OLS, MOVE techniques, KTRL, KTRL2, and RLOC.

2.1. Ordinary Least Squares Regression (OLS)

The OLS regression technique portrays the covariation between the dependent (response) WQ variable (y) and independent (predictor) WQ variables (x). The OLS c O L S and s O L S estimates obtained by solving normal equations to minimize the squared error in the estimated y values are illustrated as follows [30]:
s O L S = r s y c s x c
  c O L S = y c ¯ s O L S   x c ¯
where y c ¯ , x c ¯ are the mean values of y c and x c , whereas s y c and s x c are the standard deviation values of y c and x c , which reflect the period of concurrent records, and r is their correlation coefficient for the series of concurrent records.
OLS has five assumptions: the dependent (y) and independent (x) variables are linearly related (highly correlated), the data are representative, and the model residuals are homoscedastic, independent, and normally distributed [29,31]. The OLS is unbiased with a small mean square error M S E = σ y 2 1 ρ x y for the estimation of a missing record [12], where σ y 2 is the y population variance and ρ x y is the population correlation coefficient between x and y. However, as a record-extension technique, it creates extended records with an underestimated variance [9,12,29,32]. In addition, as the OLS c and s are based on statistical parameters (mean and standard deviation) that are clearly altered by the existence of outliers. The presence of censored values, outliers and deviation from normality are three of the main WQ data characteristics that make OLS not the ideal technique for WQ record extension or restoration of missing values.

2.2. Maintenance of Variance Extension Techniques (MOVE)

MOVE approaches have the main advantage of producing extended records with unbiased variance, which overcome one of the OLS’s major flaws [9,10]. In MOVE1, c and s were created so that the whole estimated time series of y ^ i for i = 1 ; 2 ; . n 1 ; n 1 + 1 ;   n 1 + n 2 have a similar mean and standard deviation to those of the concurrent period, y c ¯ and s yc , respectively. MOVE1 c M 1 and s M 1 are defined as follows [9]:
  s M 1 = s i g n   r s y c s x c
  c M 1 = y c ¯ s M 1   x c ¯
where s i g n   r stands for the algebraic sign (+ or −) of the r . For MOVE2, proposed by [9], the c M 2 and s M 2 estimates were determined so that if MOVE2 is employed to estimate the entire sequence of y ^ i for i = 1 ; 2 ; n 1 ; n 1 + 1 ;   n 1 + n 2 , the unbiased population estimates of the mean μ ^ y and variance ( σ ^ y 2 ) introduced by [33] could be reproduced [9] , and the MOVE2 c M 2 and s M 2 estimates are defined as follows [9]:
s M 2 = σ ^ y s x c
  c M 2 = μ ^ y s M 2   x ¯
where the μ ^ y and σ ^ y 2 estimates are defined as follows [33]:
μ ^ y = y c ¯ + s i g n r n 2 n 1 + n 2 β ^ x E ¯ x c ¯
σ ^ y 2 = β ^ 2 s x E 2 + n 1 + n 2 3 n 1 3 n 1 + n 2 1 n 1 1 n 1 2 s y c 2 β ^ s x c 2
where x E ¯ is the mean value, s x E 2 is the estimate of variance based on the full series ( x i ) for i = n 1 + 1 ;   n 1 + n 2 , and β ^ is defined as follows:
β ^ = i = 1 n 1 x i x ¯ c y i y ¯ c i = 1 n 1 x i x ¯ c
Vogel and Stedinger [10] proved that the MOVE1 and MOVE2 strategies are unable to achieve their goal mainly because record extension techniques are only applied to extend records y ^ i only, for i = n 1 + 1 ;   n 1 + n 2 , and not for the entire series that includes the period of concurrent records. Vogel and Stedinger [10] proposed c M 3 and s M 3 estimates for the MOVE3 technique that can be employed to generate extended records so that the resulting sequence of n 1 + n 2 records { y 1 ,   y 2 , ……., y n 1 , y ^ n 1 + 1 , ……, y ^ n 1 + n 2 } has the unbiased μ ^ y and σ ^ y 2 proposed by [33], as follows [10]:
s M 3 = n 1 + n 2 1 σ ^ y 2 n 1 1 s y c 2 n 1 y ¯ c μ ^ y 2 n 2 a μ ^ 2 ( n 2 1 ) s xE 2
c M 3 = n 1 + n 2 μ ^ y + n 1 y ¯ c n 2 s M 3   x ¯ c  
Vogel and Stedinger [10] also showed that the μ ^ y and σ ^ y 2 proposed by [33] and used in MOVE2 and MOVE3 are not necessarily the best asymptotic estimators for small samples of interest, which are common to hydrology in general and WQ in particular. Vogel and Stedinger [10] provided improved population estimators μ ^ y * and σ ^ y * 2 for small samples, as follows:
μ ^ y * = y ¯ c + θ 1 n 2 n 1 + n 2 β ^ x ¯ E x ¯ c
where   θ 1 = ( n 1 1 ) ρ 2 n 1 4 ρ 2 + 1
σ ^ y * 2 = 1 θ 2 s y c 2 + θ 2 β ^ 2 s x 2 + 1 n 1 + n 2 3 n 1 3 n 1 + n 2 1 n 1 1 n 2 2 s y c 2 β ^ s xc 2
where θ 2 = ( n 1 4 ) ρ 2 n 1 8.5 ρ 2 + 4.5 where ρ is the population correlation coefficient. Thus, MOVE4 was proposed by [10] to create extended records so that the ensuing sequence of n 1 + n 2 records { y 1 ,   y 2 , ……., y n 1 , y ^ n 1 + 1 , ……, y ^ n 1 + n 2 } has a mean μ ^ y * and variance σ ^ y * 2 . The MOVE4 c M 4 and s M 4 estimates are obtained by replacing the [33] estimators μ ^ y and σ ^ y 2 in Equations (11) and (12) by μ ^ y * and σ ^ y * 2 , respectively.
MOVE slopes and intercepts are based on the x and y parameters, which are visibly affected by outliers, despite the fact that the main benefit of MOVE approaches is the preservation of variance in extended records [4,8,27,28,29]. As a result, the MOVE approaches, like OLS, are sensitive to the existence of outliers.

2.3. Kendall–Theil Robust Line (KTRL & KTRL2)

The KTRL slope ( s K T ) estimate is based on the Kendall rank correlation coefficient [4,27,28,29]. s K T is based on a pairwise comparison between each pair of records ( x i y i ) and all other pairs ( x j y j ) . For each pair of comparisons, a slope of Δy/Δx is calculated, a concurrent record of size n pairs results in n(n − 1)/2 calculated slopes, and the median of all the pairwise slopes is s K T [34]:
s K T = m e d i a n y j y i x j x i i < j   i = 1 ,   2 , n 1 1   a n d   j = 2 ,   3 , , n 1
The KTRL intercept ( c K T ) is defined as follows [35]
c K T = m e d i a n ( y c ) s K T * m e d i a n ( x c )
In the KTRL2 proposed by [4], the slope estimate ( s K T 2 ) was developed so that the extended records would have a cumulative distribution function (CDF) similar to the CDF of the measured records. The KTRL2 s K T 2 and c K T 2 estimates are defined as follows [4]:
s K T 2 = m e d i a n y j y i x j x i i < j   i = 5 t h ,   10 t h , 90 t h   a n d   j = 10 t h ,   15 t h , , 95 t h
c K T 2 = m e d i a n ( y c ) s K T 2 * m e d i a n ( x c ) .  
where y j and x j   . are the percentiles of y c and x c estimated during the period of concurrent records. Percentiles are calculated for the 5th, 10th..., and 95th percentiles, according to [4]. As a result, 171 [n (n1)/2 = 19 (19 − 1)/2] pairwise comparisons will come from a set of 19 ( x j , y j ) percentile pairs. A slope y x is produced for each of these comparisons, and the median of the 171 possible pairwise slopes is used as the slope estimate. As a result, rather than minimizing the error in the estimation of the y records, the goal was to minimize the error in the estimation of the y percentiles.

2.4. Robust Line of Organic Correlation (RLOC)

Khalil and Adamowski [5] introduced the RLOC as a modified version of MOVE1 with the benefit of being insensitive to the existence of outliers. The RLOC slope ( s R L ) and intercept ( c R L ) estimates are based on the median and interquartile range (IQR) ratios, as follows [5] :
s R L = y 75 y 25 x 75 x 25  
c R L = m e d i a n ( y c ) s R L * m e d i a n ( x c )
where y 75 ,   y 25 ,   x 75 ,   x 25 are the 75th and 25th percentiles of y and x measured concurrent records. When the data are normally distributed, [5] found that MOVE1 and RLOC have approximately the same performance, but the RLOC estimator becomes more accurate with a minor deviation from normality.

2.5. Wavelet Transform

The wavelet transform (WT) was initially developed in the mathematics community but has been proven to be a useful tool for analyzing nonstationary time series in hydrology and hydrogeology (e.g., [36,37,38,39,40,41]). The WT is a time and frequency domain multiresolution analysis that is an important derivative of the Fourier transforms (FT) [42]. The WT has advantages over the FT, as it can simultaneously gather information on the time, location, and frequency of a signal, whereas the FT only offers frequency information [42].
The À Trous-Haar WT was chosen as the most suitable discrete WT (DWT) for forecasting tests in this study [43]. Du et al. [44] have shown that although the DWT is quite popular, it is usually misemployed in hydrological forecasting as a data-preprocessing step. As with DWT, different algorithms (such as maximal overlap DWT-multiresolution analysis (MODWT-MRA) and DWT-MRA) use future values to compute the approximation and detail components [44]. However, the À Trous algorithm does not rely on future data in the decomposition process [45].
The WT (Ws) (details component) is calculated using the equations below [46] when C0 is the original time series and   C S is the approximation component at scale s:
W S k = C S 1 k C S k
C S k = l = + h l C s 1 k + 2 s l
where k is the wavelet transform calculation location (inside the time series), h(l) is the low pass filter, and l = (1/16, 1/4, 3/8, 1/4, 1/16). The À Trous-Haar wavelet transform proposed by [47] should be replaced with the standard À Trous WT to account for the requirement that future data values cannot be included in the wavelet transform calculation. The À Trous-Haar wavelet algorithm is identical to the À Trous wavelet method, with the exception that the low-pass filter h(l) is substituted with a simpler filter, with l = (0.5, 0.5). As a result, Equation (19) is [48]:
C S k = 1 2 ( C S 1 k C S 1 k 2 s )
As a result, the information after k is not used in calculating the wavelet coefficient at any time point k. In this study, the À Trous-Haar algorithm was applied to the predictor x as a data preprocess, and then the record-extension techniques were applied using the predictor approximation component to estimate the desired y records.

3. Empirical Experiment

The Edko drainage system is one of the Nile Delta’s major drainage systems. Beginning from Shubra-Kheitto free flow Lake Edko before reaching the Mediterranean, the Edko catchment area is approximately 96,000 hectares with a length of 48.8 km [49]. Since August 1997, monthly samples have been gathered at 11 water quality monitoring locations across the Edko drainage system (Figure 1).
The empirical experiment used ten years of monthly WQ measured records for Sodium (Na) and Total Dissolved Solids (TDS) measured at 10 out of the 11 monitoring locations. One location was excluded due to an incomplete Na dataset. It is assumed that 10 years of monthly records are representative of the Na and TDS measured at these 10 monitoring locations. It should be emphasized that the selection of the Na and TDS for this case study was due to their high correlation.
Table 1 shows descriptive statistics for the Na and TDS measured at the 10 monitoring locations. Figure 2 displays the box plots for Na and TDS, which reveal outliers and positive skewness. Table 2 displays the findings of the Kolmogorov–Smirnov goodness-of-fit test, which assesses the normality of the data under the null hypothesis that the sample was drawn from a normal distribution. The results show that none of the two WQ variables at any of the ten monitoring stations allow the test null hypothesis to be accepted (Table 2).
The TDS was employed as an explanatory variable (predictor or independent variable) to extend the Na records (response or dependent variable) using the eight record extension approaches. The effectiveness of the eight record-extension approaches, as well as the use of the À Trous-Haar WT data preprocessing step, were evaluated using a split-sample cross-validation method. In the split sample cross-validation, every two years of monthly records were eliminated n 2 = 24 from the ten years of available data and the remaining eight years were used to define the period of concurrent records. Four different sizes were considered for the period of concurrent records n 1 equal to 60, 70, 80, and 90 records.
The eight record-extension approaches were used for the estimation of Na using TDS as a predictor at each of the Edko drain 10 monitoring locations considered in this study. As a result, 200 different extended Na record realizations were constructed (10 locations × 5 different sample combinations x four different n 1 sizes = 200).
For each of the 200 distinct realizations studied, the correlation coefficient of y c and x c was consistently positive and ranged between 0.77 and 0.84.
The bias (BIAS) as accuracy metric and the root mean squared error (RMSE) as a precision metric were used to assess the performances of the eight record-extension approaches under the À Trous-Haar data preprocess step:
B I A S = 1 n t i = 1 n t S i ^ S i
R M S E = 1 n t i = 1 n t S i ^ S i 2
where S i and S i ^ are the measured and estimated statistics of the dependent variable for i = 1 ; n t , respectively, and n t is the number of trials in the empirical study.

4. Results and Discussion

Figure 3, Figure 4 and Figure 5 show the BIAS and RMSE values for the estimation of the extended records (Figure 3), extended records’ mean value (Figure 4), and extended records’ standard deviation (Figure 5).
For the restoration of missing records (Figure 3), the OLS and MOVE techniques showed comparable results, indicating that the WT step had no clear influence. However, for KTRL, KTRL2 and RLOC, the WT step deteriorates the accuracy and precision. These results are due to the high correlation between the two WQ variables (Na & TDS); high linear dependence is one of the assumptions of regression techniques, which may mask the influence of the data preprocessing step for the restoration of a missing record. High precision was generally observed for the KTRL and MOVE techniques, whereas high accuracy was generally observed for OLS, MOVE3, and MOVE4. The KTRL provided the most precise results for the period of concurrent records equal to 90 and 80, whereas MOVE2 was the most precise for concurrent records of sizes equal to 70 and 60. As the period of concurrent records increases, the data become more representative, and nonparametric techniques demonstrated their performance. This dynamic regarding the period of concurrent records is clear not only for the KTRL performance but also for all eight techniques with and without the WT data preprocessing step, as shown in Figure 3.
These results are in agreement with [2] results for the restoration of stream flow missing values, where the use of the À Trous-Haar WT as a preprocessing step did not show consistent improvement. Based on these results and results obtained by [2] and given the limited WQ data (10 years in this study) compared to streamflow data (40 years used by [2]), we may confirm that the use of the À Trous Haar WT as a data preprocessing step did not show clear improvement for the restoration of missing WQ values.
Figure 4 shows the BIAS and RMSE results for the estimation of the mean value of the extended records. In general, the OLS and MOVE techniques showed comparable results with their WT-based counterparts, which indicates that the WT step had no clear influence. However, for KTRL, KTRL2 and RLOC, the WT step deteriorates the accuracy and precision. Given that the regression techniques are proposing a line that passes by the mean (X) and mean (Y) for the OLS, and median (X) and median (Y) for the KTRL, the mean value of the estimated records should be accurate and precise. Given the linear dependence between the TDS and Na, the OLS and MOVE approaches provided more precise and accurate results with almost no influence of the WT data preprocessing step. The period of available records affects more the KTRL, KTRL2, and RLOC approaches; as the period of available records decreases, the accuracy and precision deteriorate. High precision was generally observed for the KTRL, MOVE3, and MOVE4 approaches, whereas high accuracy was generally observed for OLS, MOVE3, and MOVE4.
Figure 5 shows the BIAS and RMSE results for the estimation of the Na expanded records’ standard deviation values. Relatively high accuracy was obtained for the WT-MOVE3, WT-MOVE4, WT-KTRL2 and WT-RLOC techniques, whereas the WT-KTRL2 was the most precise for relatively large-size concurrent records (e.g., n1 = 90), and the WT-MOVE1 provided more precise results for small sizes (e.g., 60 records). Higher accuracy of the extended records’ standard deviation provided by the MOVE3 and MOVE4 approaches than the MOVE1 and MOVE2 is noted. This is attributed to the main difference between those techniques, where MOVE3 and MOVE4 were developed for small-size samples, unlike MOVE1 and MOVE2, which are based on population parameters. These results are in agreement with the results obtained by [10,27], where the four MOVE techniques were compared using streamflow and WQ data, respectively. This indicates the usefulness of the MOVE3 and MOVE4 for such small period of WQ records. The better accuracy provided by the KTRL2, RLOC, WT-KTRL2, and WT-RLOC approaches is due to the advantage of being insensitive to the existence of outliers and the ability to maintain variance of the extended records. These results are in agreement with results obtained by [27,28], which showed the better accuracy and precision provided by the KTRL2 and RLOC compared to MOVE techniques in the presence of outliers. For the extended records’ standard deviation, the WT data preprocessing step smooths the raw data into an approximation component that minimizes the influence of extreme values, which affects the extended records’ variance more than their mean value.
Figure 6 and Figure 7 show the BIAS and RMSE results for the Na extended records’ full length of percentiles from the 5th to the 95th percentile, respectively. Figure 8 and Figure 9 provide more focus on the BIAS and RMSE for the low and high percentiles, respectively. In general, Figure 6, Figure 7, Figure 8 and Figure 9 show that the accuracy and precision of the extended record percentiles increase as the period of concurrent records increases. This is attributed to the first assumption in regression techniques that the dependent and independent data are representative; as the period of concurrent records increases, the data become more representative, and their sample parameters become more accurate and precise.
For the four periods of concurrent records considered, the WT-KTRL provided extended records that are more accurate and precise in presenting the extended period low percentiles. The WT-KTRL showed either similar or slightly better performance than the KTRL. However, for high percentiles (e.g., 85th, 90th and 95th), both the KTRL and WT-KTRL techniques showed the least accuracy and precision (e.g., n1 = 60 and 70) or the second least accuracy and precision for n1 = 80 or 90.
These results are not in agreement with previous results provided by [4], which showed that the KTRL provided extended records that overestimated low percentiles and underestimated high percentiles. These different results are attributed to the existence of outliers in the independent variables (x) more/not in the dependent variable (y), as shown in Figure 2. The heavy existence of outliers in the independent variable relative to the dependent variable affects the variance and standard deviation, which becomes larger than it should be. A clear explanation of the OLS or MOVE1 slope estimates is that an overestimated (fault larger) standard deviation for the independent variable leads to an underestimation of the regression slope. Underestimation of the regression slope provides underestimation of the extended records, which leads to underestimating low as well as high values. For OLS, WT-OLS, KTRL, and WT-KTRL, underestimation of the regression slope leads to reduced overestimation of low values (low percentiles) but simultaneously increases underestimation of high values (high percentiles). Similarly, underestimating the MOVE slope results in underestimating high and low values, resulting in underestimating both high and low percentiles.
In evaluating the influence of the WT step, WT-KTRL and WT-OLS showed slightly more accurate and precise results than KTRL and OLS, respectively. For the MOVE techniques, for both the high and low percentiles, the WT-step did not show clear improvement but reduced the accuracy and precision. For KTRL2 and RLOC, accuracy and precision improvements were detected with the WT data preprocessing step.
For the four different periods of concurrent records, the highest accuracy in estimating high percentiles was obtained by WT-MOVE3 and WT-MOVE4, followed by the MOVE1-WT technique. The highest precision was obtained by the WT-MOVE1, WT-MOVE3, and WT-MOVE4 techniques, while WT-KTRL2 became the second highest precision technique when the period of concurrent records was 90.
Thus, when the period of concurrent records was large enough to be considered representative, KTRL2 and WT-KTRL2 showed accurate and precise estimates in representing their statistical parameters. However, for restoration of missing records, the KTRL showed more accurate and precise results for small-sample WQ data.
In brief, results showed that using the À Trous-Haar WT as a data preprocessing step did not improve the accuracy or precision of any of the eight regression approaches for the estimation of missing WQ records. However, for the extension of short-gauged WQ records, using the À Trous-Haar WT step improved the accuracy and precision of the extended records’ standard deviation and extreme percentiles. Thus, the use of the À Trous-Haar WT as a data preprocessing step is recommended if the objective is to extend WQ records at short-gauged sites. In addition, and in agreement with previous studies [4,28], results showed that the KTRL is preferable for the restoration of missing WQ values in case of the presence of outliers. For the extension of WQ records, the WT-MOVE3 or WT-MOVE4 is preferable, while in the presence of outliers, the WT-KTRL2 is desirable.

5. Conclusions and Recommendations

For the restoration of WQ missing data and the extension of WQ records at short-gauged locations, eight record-extension strategies were investigated in this study. The evaluation took into account a data preprocessing step that used the À Trous-Haar WT to denoise the predictor by just using the À Trous-Haar approximation component in the regression technique. Real WQ records from the Edko drainage system in Egypt’s Nile Delta were used in an empirical investigation. The results showed that adding the WT data preprocessing step did not improve the restoration of missing values significantly. However, for the extended records’ standard deviation and percentiles, it showed improvements in accuracy and precision.
It can be concluded that the selection of the appropriate record-extension technique is based on two main aspects: the existence of outliers; and the objective of the record substitution, either to estimate missing values or to extend short-gauged records. The existence of outliers should be checked carefully before selecting an applicable record-extension technique, especially if the objective is to provide extended records that preserve the statistical parameters. For the restoration of missing WQ values, either the OLS or the KTRL should be used, with said KTRL being preferred in the existence of outliers. For the extension of WQ records at short-gauged sites, any of the MOVE, KTRL2, or RLOC procedures can be used, with WT-KTRL2 and WT-RLOC preferable in the existence of outliers and WT-MOVE3 and WT-MOVE4 preferable in cases with small sample sizes.
The use of the À Trous-Haar WT as a data preprocessor for record-extension approaches requires more research, including case studies from various regions and an assessment of the impact of outliers, their position, intensity, and magnitudes. Furthermore, analyzing the impact of the period of concurrent records and the level of correlation between the dependent and independent variables on the use of the À Trous-Haar WT as a data preprocessing step would provide a clear assessment. In addition, as a natural extension of this work, other data decomposition techniques such as empirical mode decomposition (EMD) and ensemble EMD (EEMD) (e.g., [50,51,52,53]) would be used instead of WT and would allow comparing different data preprocessing methods.

Author Contributions

Conceptualization, S.A., B.K. and A.E.S.; methodology, S.A., B.K. and A.E.S.; software, S.A. and A.E.; investigation, S.A.; data curation, A.E.S.; writing—original draft preparation, S.A.; writing—review and editing, B.K., A.E.S., A.E. and M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The water quality data used in this study were provided by the National Water Research Center (NWRC) of Egypt. The authors are grateful to the Editor, the Associate Editor, and the three anonymous reviewers, whose comments and suggestions greatly improved the quality of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Halbe, J.; Pahl-Wostl, C.; Sendzimir, J.; Adamowski, J. Towards adaptive and integrated management paradigms to meet the challenges of water governance. Water Sci. Technol. 2013, 67, 2651–2660. [Google Scholar] [CrossRef]
  2. Nalley, D.; Adamowski, J.; Khalil, B.; Biswas, A. A comparison of conventional and wavelet transform based methods for streamflow record extension. J. Hydrol. 2019, 582, 124503. [Google Scholar] [CrossRef]
  3. Horne, J. Water Information as a Tool to Enhance Sustainable Water Management—The Australian Experience. Water 2015, 7, 2161–2183. [Google Scholar] [CrossRef] [Green Version]
  4. Khalil, B.; Ouarda, T.B.M.J.; St-Hilaire, A. Comparison of Record-Extension Techniques for Water Quality Variables. Water Resour. Manag. 2012, 26, 4259–4280. [Google Scholar] [CrossRef]
  5. Khalil, B.; Adamowski, J. Record extension for short-gauged water quality parameters using a newly proposed robust version of the Line of Organic Correlation technique. Hydrol. Earth Syst. Sci. 2012, 16, 2253–2266. [Google Scholar] [CrossRef] [Green Version]
  6. Hadzima-Nyarko, M.; Rabi, A.; Šperac, M. Implementation of Artificial Neural Networks in Modeling the Water-Air Temperature Relationship of the River Drava. Water Resour. Manag. 2014, 28, 1379–1394. [Google Scholar] [CrossRef]
  7. Nevitt, J.; Tam, H.P. A Comparison of Robust and Nonparametric Estimators under the Simple Linear Regression Model. In Proceedings of the Paper presented at the Annual Meeting of the American Educational Research Association, Chicago, IL, USA, 24–28 March 1997. [Google Scholar]
  8. Granato, G.E. United States, Office of the Natural and Human Environment, Geological Survey (U.S.). In Kendall-Theil Robust Line (KTRLine-Version 1.0): A Visual Basic Program for Calculating and Graphing Robust Nonparametric Estimates of Linear-Regression Coefficients between Two Continuous Variables; Techniques and methods of the, U.S. Geological Survey; Bibliogov: Washington, DC, USA, 2006; Volume 4, p. 31. [Google Scholar]
  9. Hirsch, R.M. A comparison of four streamflow record extension techniques. Water Resour. Res. 1982, 18, 1081–1088. [Google Scholar] [CrossRef]
  10. Vogel, R.M.; Stedinger, J.R. Minimum variance streamflow record augmentation procedures. Water Resour. Res. 1985, 21, 715–723. [Google Scholar] [CrossRef] [Green Version]
  11. Robinson, R.; Wood, M.S.; Smoot, J.L.; Moore, S.E. Parametric modeling of water quality and sampling strategy in a high-altitude appalachian stream. J. Hydrol. 2004, 287, 62–73. [Google Scholar] [CrossRef]
  12. Koutsoyiannis, D.; Langousis, A. Treatise on Water Science Chapter 27: Precipitation. Treatise Water Sci. 2011, 2, 27–78. [Google Scholar] [CrossRef]
  13. Canadian Council of Ministers of the Environment. 2015. Available online: https://www.ccme.ca/en/resources/water/water_quality.html (accessed on 1 June 2015).
  14. Lettenmaier, D.P. Multivariate Nonparametric Tests for Trend in Water Quality. J. Am. Water Resour. Assoc. 1988, 24, 505–512. [Google Scholar] [CrossRef]
  15. Berryman, D.; Bobee, B.; Cluis, D.; Haemmerli, J. Nonparametric metric tests for trend detection in water quality time series. J. Am. Water Resour. Assoc. 1988, 24, 545–556. [Google Scholar] [CrossRef]
  16. Dietz, E.J. A comparison of robust estimators in simple linear regression. Commun. Stat. Simul. Comput. 1987, 16, 1209–1227. [Google Scholar] [CrossRef]
  17. Albek, E. Estimation of point and diffuse contaminant loads to streams by non-parametric regression analysis of monitoring data. Water Air Soil Pollut. 2003, 147, 229–243. [Google Scholar] [CrossRef]
  18. Olsson, O.; Gassmann, M.; Wegerich, K.; Bauer, M. Identification of the effective water availability from streamflows in the Zerafshan river basin, Central Asia. J. Hydrol. 2010, 390, 190–197. [Google Scholar] [CrossRef]
  19. Dery, S.J.; Mlynowski, T.J.; Hernandez-Henriquez, M.A.; Straneo, F. Interannual variability and interdecadal trends in Hudson Bay streamflow. J. Mar. Syst. 2011, 88, 341–351. [Google Scholar] [CrossRef]
  20. Vorogushyn, S.; Merz, B. Flood trends along the Rhine: The role of river training. Hydrol. Earth Syst. Sci. 2013, 17, 3871–3884. [Google Scholar] [CrossRef] [Green Version]
  21. Calarullo, S.J.; Sullivan, S.L.; McHugh, A.R. Implementation of MOVE.1, Censored MOVE.1, and Piecewise MOVE.1 Low-Flow Regressions with Applications at Partial-Record Streamgaging Stations in New Jersey: U.S. Geological Survey Open-File Report 2018-1089. 20p. Available online: https://pubs.usgs.gov/of/2018/1089/ofr20181089.pdf (accessed on 6 November 2018).
  22. Jia, Y.; Culver, T.B. Bootstrapped artificial neural networks for synthetic flow generation with a small data sample. J. Hydrol. 2006, 331, 580–590. [Google Scholar] [CrossRef]
  23. Ryu, J.H.; Svoboda, M.D.; Lenters, J.D.; Tadesse, T.; Knutson, C.L. Potential extents for ENSO-driven hydrologic drought forecasts in the United States. Clim. Chang. 2009, 101, 575–597. [Google Scholar] [CrossRef]
  24. Raziei, T.; Saghafian, B.; Paulo, A.A.; Pereira, L.S.; Bordi, I. Spatial Patterns and Temporal Variability of Drought in Western Iran. Water Resour. Manag. 2008, 23, 439–455. [Google Scholar] [CrossRef] [Green Version]
  25. Raziei, T.; Bordi, I.; Pereira, L.S. An application of GPCC and NCEP/NCAR datasets for draught variability analysis in Iran. Water Resour. Manag. 2011, 25, 1075–1086. [Google Scholar] [CrossRef]
  26. Khalil, B.; Awadallah, A.G.; Adamowski, J.; Elsayed, A. A Novel Record-Extension Technique for Water Quality Variables Based on L-Moments. Water Air Soil Pollut. 2016, 227, 179. [Google Scholar] [CrossRef]
  27. Khalil, B.; Adamowski, J. Comparison of OLS, ANN, KTRL, KTRL2, RLOC, and MOVE as record-extension techniques for water quality variables. Water Air Soil Pollut. 2014, 225, 1966. [Google Scholar] [CrossRef]
  28. Khalil, B.; Adamowski, J. Evaluation of the performance of eight record-extension techniques under different levels of association, presence of outliers and different sizes of concurrent records: A Monte Carlo study. Water Resour. Manag. 2014, 28, 5139–5155. [Google Scholar] [CrossRef]
  29. Helsel, D.; Hirsch, R.M. Statistical Methods in Water Resources (USGS Numbered Series No. 04-A3), Statistical Methods in Water Resources, Techniques of Water-Resources Investigations; U.S. Geological Survey: Reston, VA, USA, 2002. [Google Scholar] [CrossRef]
  30. Draper, N.R.; Smith, H. Applied Regression Analysis; John Wiley: New York, NY, USA, 1966; 736p. [Google Scholar]
  31. Serinaldi, F.; Grimaldi, S.; Abdolhosseini, M.; Corona, P.; Cimini, D. Testing copula regression against benchmark models for point and interval estimation of tree wood volume in beech stands. Forstwiss. Centralblatt 2012, 131, 1313–1326. [Google Scholar] [CrossRef] [Green Version]
  32. Moog, D.B.; Whiting, P.J.; Thomas, R.B. Streamflow record extension using power transformations and application to sediment transport. Water Resour. Res. 1999, 35, 243–254. [Google Scholar] [CrossRef]
  33. Matalas, N.; Jacobs, B. A Correlation Procedure for Augmenting Hydrologic Data; (USGS Numbered Series No. 434-E), a Correlation Procedure for Augmenting Hydrologic Data, Professional Paper; U.S. Government Printing Office: Washington, DC, USA, 1964. [Google Scholar] [CrossRef]
  34. Theil, H. A rank-invariant method of linear and polynomial regression analysis I and II. Indag. Math. 1950, 12, 173. [Google Scholar]
  35. Groeneveld, R.A.; Conover, W.J. Practical Nonparametric Statistics, 2nd ed.; John Wiley and Sons: New York, NY, USA, 1980; 493p. [Google Scholar]
  36. Partal, T.; Kişi, Ö. Wavelet and neuro-fuzzy conjunction model for precipitation forecasting. J. Hydrol. 2007, 342, 199–212. [Google Scholar] [CrossRef]
  37. Partal, T.; Cigizoglu, K. Estimation and forecasting of daily suspended sediment data using wavelet-neural networks. J. Hydrol. 2008, 358, 317–331. [Google Scholar] [CrossRef]
  38. Adamowski, J.; Chan, H.F. A wavelet neural network conjunction model for groundwater level forecasting. J. Hydrol. 2011, 407, 28–40. [Google Scholar] [CrossRef]
  39. Pandey, B.K.; Tiwari, H.; Khare, D. Trend analysis using discrete wavelet transform (DWT) for long-term precipitation (1851–2006) over India. Hydrol. Sci. J. 2017, 62, 2187–2208. [Google Scholar] [CrossRef]
  40. Graf, R.; Zhu, S.; Sivakumar, B. Forecasting river water temperature time series using a wavelet–neural network hybrid modelling approach. J. Hydrol. 2019, 578, 124115. [Google Scholar] [CrossRef]
  41. Liu, Q.-J.; Shi, Z.-H.; Fang, N.-F.; Zhu, H.-D.; Ai, L. Modeling the daily suspended sediment concentration in a hyperconcentrated river on the Loess Plateau, China, using the Wavelet–ANN approach. Geomorphology 2013, 186, 181–190. [Google Scholar] [CrossRef]
  42. Tiwari, M.K.; Chatterjee, C. A new wavelet–bootstrap–ANN hybrid model for daily discharge forecasting. J. Hydroinform. 2010, 13, 500–519. [Google Scholar] [CrossRef] [Green Version]
  43. Mallat, S. A Wavelet Tour of Signal Processing, 2nd ed.; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  44. Du, K.; Zhao, Y.; Lei, J. The incorrect usage of singular spectral analysis and discrete wavelet transform in hybrid models to predict hydrological time series. J. Hydrol. 2017, 552, 44–51. [Google Scholar] [CrossRef]
  45. Renaud, O.; Starck, J.-L.; Murtagh, F. Wavelet-Based Combined Signal Filtering and Prediction. IEEE Trans. Syst. Man Cibern. Part B 2005, 35, 1241–1251. [Google Scholar] [CrossRef]
  46. Shensa, M. The discrete wavelet transform: Wedding the a trous and Mallat algorithms. IEEE Trans. Signal Process. 1992, 40, 2464–2482. [Google Scholar] [CrossRef] [Green Version]
  47. Zheng, G.; Starck, J.-L.; Campbell, J.; Murtagh, F. Multiscale transforms for filtering financial data streams. J. Comput. Intell. Financ. 1999, 7, 18–35. [Google Scholar]
  48. Murtagh, F.; Starck, J.L.; Renaud, O. On neuro-wavelet modeling. Decis. Support Syst. 2004, 37, 475–484. [Google Scholar] [CrossRef]
  49. El-Saadi, A. Economics and Uncertainty Considerations in Water Quality Monitoring Networks Design. Ph.D. Dissertation, Faculty of Engineering, Ain-Shams University, Cairo, Egypt, 2006. [Google Scholar]
  50. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.-C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  51. Nelsen, B.; Williams, D.A.; Williams, G.P.; Berrett, C. An Empirical Mode-Spatial Model for Environmental Data Imputation. Hydrology 2018, 5, 63. [Google Scholar] [CrossRef] [Green Version]
  52. Eze, E.; Halse, S.; Ajmal, T. Developing a Novel Water Quality Prediction Model for a South African Aquaculture Farm. Water 2021, 13, 1782. [Google Scholar] [CrossRef]
  53. Chu, T.-Y.; Huang, W.-C. Application of Empirical Mode Decomposition Method to Synthesize Flow Data: A Case Study of Hushan Reservoir in Taiwan. Water 2020, 12, 927. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Edko drainage system WQ monitoring locations.
Figure 1. Edko drainage system WQ monitoring locations.
Water 14 02264 g001
Figure 2. Box plots for the Na and TDS records.
Figure 2. Box plots for the Na and TDS records.
Water 14 02264 g002
Figure 3. The BIAS and RMSE results for the extended Na records for four different periods of concurrent records = 60; 70; 80; and 90.
Figure 3. The BIAS and RMSE results for the extended Na records for four different periods of concurrent records = 60; 70; 80; and 90.
Water 14 02264 g003
Figure 4. The BIAS and RMSE results for the extended Na records mean value for four different periods of concurrent records = 60; 70; 80; and 90.
Figure 4. The BIAS and RMSE results for the extended Na records mean value for four different periods of concurrent records = 60; 70; 80; and 90.
Water 14 02264 g004
Figure 5. The BIAS and RMSE results for the extended Na records standard deviation for four different concurrent sizes = 60; 70; 80; and 90.
Figure 5. The BIAS and RMSE results for the extended Na records standard deviation for four different concurrent sizes = 60; 70; 80; and 90.
Water 14 02264 g005
Figure 6. The BIAS results for the extended Na records’ percentiles.
Figure 6. The BIAS results for the extended Na records’ percentiles.
Water 14 02264 g006aWater 14 02264 g006b
Figure 7. The RMSE results for the extended Na records’ percentiles.
Figure 7. The RMSE results for the extended Na records’ percentiles.
Water 14 02264 g007aWater 14 02264 g007b
Figure 8. The BIAS (a) and RMSE (b) results for the extended Na records’ low percentiles for n1 = 60 and 90.
Figure 8. The BIAS (a) and RMSE (b) results for the extended Na records’ low percentiles for n1 = 60 and 90.
Water 14 02264 g008aWater 14 02264 g008b
Figure 9. The BIAS (a) and RMSE (b) results for the extended Na records’ high percentiles for n1 = 60 and 90.
Figure 9. The BIAS (a) and RMSE (b) results for the extended Na records’ high percentiles for n1 = 60 and 90.
Water 14 02264 g009
Table 1. Na and TDS descriptive statistics.
Table 1. Na and TDS descriptive statistics.
Monitoring
Locations
MinimumMaximumMeanSt. DeviationSkewness
NaTDSNaTDSNaTDSNaTDSNaTDS
(mg/l)
WE011.13203.0016.791411.004.14667.991.96179.362.840.61
WE021.20216.0058.714155.6121.452138.9510.40772.740.520.16
WE031.16213.0068.024734.709.171118.716.57450.356.384.46
WE050.31232.0033.052978.056.31950.243.09266.015.613.43
WE061.09197.0016.961602.077.01946.122.51207.511.310.43
WE071.13210.0037.713737.009.471141.436.14514.642.662.97
WE083.00562.0062.416390.0023.732313.4714.241153.320.660.64
WE101.00259.0011.611426.005.70835.951.88168.320.860.41
WE113.52642.0044.533556.006.78913.384.14306.546.866.20
WE211.27232.0043.883624.0017.361774.257.91577.390.980.82
Table 2. Kolmogorov–Smirnov goodness of fit test for Na and TDS.
Table 2. Kolmogorov–Smirnov goodness of fit test for Na and TDS.
Monitoring Locationsp Value
NaTDS
WE013.35 × 10−911.80 × 10−102
WE026.03 × 10−991.80 × 10−102
WE035.76 × 10−961.80 × 10−102
WE054.71 × 10−941.80 × 10−102
WE063.19 × 10−991.80 × 10−102
WE076.58 × 10−991.80 × 10−102
WE083.40 × 10−1021.80 × 10−102
WE108.20 × 10−971.80 × 10−102
WE216.03 × 10−991.80 × 10−102
WE112.00 × 10−1021.80 × 10−102
WE216.03 × 10−991.80 × 10−102
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Anwar, S.; Khalil, B.; Seddik, M.; Eltahan, A.; Saadi, A.E. An Evaluation of À Trous-Based Record Extension Techniques for Water Quality Record Extension. Water 2022, 14, 2264. https://doi.org/10.3390/w14142264

AMA Style

Anwar S, Khalil B, Seddik M, Eltahan A, Saadi AE. An Evaluation of À Trous-Based Record Extension Techniques for Water Quality Record Extension. Water. 2022; 14(14):2264. https://doi.org/10.3390/w14142264

Chicago/Turabian Style

Anwar, Samah, Bahaa Khalil, Mohamed Seddik, Abdelhamid Eltahan, and Aiman El Saadi. 2022. "An Evaluation of À Trous-Based Record Extension Techniques for Water Quality Record Extension" Water 14, no. 14: 2264. https://doi.org/10.3390/w14142264

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