A Rapid UV / Vis Spectrophotometric Method for the Water Quality Monitoring at On-Farm Root Vegetable Pack Houses

: Our research aim was to apply UV / Vis spectrophotometric techniques for the rapid monitoring of the quality of water sourced from on-farm root vegetable washing processes. To achieve this goal, the quality assessment of the washing water and wastewater at di ﬀ erent stages of the technological processes was performed using physicochemical, biological, and UV / Vis absorbance measurements as well as statistical methods, such as principal component analysis (PCA) and partial least squares (PLS) regression. Limit values of UV / Vis absorbance at speciﬁc wavelengths were predicted in order to adapt them for routine testing and water quality monitoring at the farm packhouses. Results of the lab analyses showed, that the main problems of the water quality were caused by suspended solids (470–3400 mg L − 1 ), organic substances (BOD 5 215–2718 mg L − 1 ; COD 540–3229 mg L − 1 ), nitrogen (3–52 mg L − 1 ), phosphorus (1–6 mg L − 1 ), and pathogenic microorganisms (TVC > 300 cfu mL − 1 , E. coli 5.5 × 10 3 –1.0 × 10 4 cfu mL − 1 , intestinal enterococci 2.8 × 10 2 –1.5 × 10 4 cfu mL − 1 , coliform bacteria 1.6 × 10 3 –2.0 × 10 4 cfu mL − 1 ). Suspended solids exceeded the limit values by 10–50 times, organic matter by 10–25 times, dissolved organic carbon by 3–5 times, nitrogen by 3–7 times, total phosphorus by 3–12 times, and microorganisms by 3–10 times. UV / Vis limit values calculated were as follows: A210 nm—3.997–4.009 cm − 1 , A240 nm—5.193–5.235 cm − 1 , A254 nm—4.042–4.047 cm − 1 , A320 nm—7.387–7.406 cm − 1 , and A660 nm—3.937–3.946 cm − 1 . UV / Vis measurements at A320 nm are proposed for the routine water quality monitoring. organic


Introduction
Washing water and wastewater quality control in on-farm packhouses is of great importance for the quality of agricultural production as well as the environment. Producers require rapid and efficient methods for water quality monitoring. Conventional standard analytical methods are expensive, time-consuming and require special laboratory conditions and equipment. Therefore they are inconvenient and impossible to use for the water quality monitoring at the on-farm packhouses. Ultraviolet/visual spectrometry (UV/Vis) can have significant advantages over other in situ methods for this application due to its wide availability and low application costs. Previously, various UV/Vis methods were used, for example, to determine the particle size of polystyrene in solution, the size distribution of quantum particles of a ZnO semiconductor, the size of gold-coated palladium clusters in solution, and to study the aggregation of b-carotene particles [1]. UV/Vis technology is rapid and Y = ε 1 c 1 l + ε 2 c 2 l + · · · + ε p c p l + E, (2) where Y is the response of a particular sample at wavelength λ, ε p is the molar extinction coefficient for the p sample, c p is the corresponding concentration and E is the spectral error. If ε is measured at several wavelengths, then the equation obtained is: where j is the index for samples, i is the index for wavelengths, p is the number of components to be determined and E y is spectral errors. In matrix notation, Equation (3) can be written as follows: where Y is the j·i matrix of calibration spectra, C is the j·l matrix of component concentration, K is the l·i matrix of absorptivity-path length products, and E Y is the j·i matrix of spectral errors.
The limitation of this model is that it can only be applied to up to three variables. Multivariate calibration methods enable analyzing those analytes that interact with each other and predict their concentration from spectral data [12][13][14]. Chemometrics enables exploring, analyzing, and modeling relationships among different water quality parameters by various statistical methods [13][14][15][16][17][18][19]. Generally, direct chemometrics models can only be used if all constituents' spectra are known and Beer-Lambert law is valid. However, wastewaters are complex mixtures of various soluble and insoluble organic and inorganic compounds interacting and transforming in time. Therefore, hidden isosbestic points and simple regression models can not completely describe the situation, especially if monitored substances chemically interfere with each other influencing the final spectra [19]. Moreover, specimens from 'real-life' sources may contain substances other than those of interest, contributing to the absorbance. Therefore, the use of direct chemometrics methods for wastewater characterization is limited [5,20].
Indirect chemometrics techniques, such as principal component analysis (PCA), multiple and multivariate regression, partial least square regression (PLS), and others, are based on estimating parameters from sample mixtures. Their are not dependent on the validity of Beer-Lambert law and systems of the even highly correlated spectra can be used for multi-component calibration [14,18,21,22]. Applying these techniques, it is possible to determine the UV/Vis absorption limits for wastewater quality according to the limit values of chemical analytes.
Therefore, in our study, chemometric techniques such as principal component analysis (pca) and partial least squares regression and projection to latent structures (PLS) were applied to analyze relationships between water quality chemical parameters (suspended solids (SS), electric conductivity (EC), dissolved organic carbon (DOC), chemical oxygen demand COD Cr , biochemical oxygen demand (BOD 5 )), and UV/Vis absorbance data. The aim was to establish a simple and reliable method for monitoring water quality which, sourced from the root vegetables washing and wastewater treatment processes within the on-farm packhouses.
Statistical models can increase the importance and scope of the obtained spectral data [23]. Partial least squares regression (PLS) is a multidimensional statistical method used to predict standard values of aggregated water quality parameters such as total organic carbon (TOC), chemical oxygen demand (COD) and total suspended solids (TSS) [24,25]. The use of regression algorithms must comply with certain assumptions, i.e., the data should be independent and have a normal distribution. Predictive abilities of PLS calibration models were evaluated in relation to the determination factor (R 2 ) [26].
Different statistics were used to assess the PLS model and performance indicators were proposed to classify the developed models. It is essential to consider whether UV/Vis spectra can develop robust TSS and COD estimators, with statistical regression carried out on the off-line spectra in the collected samples. This affects the reliability of the in situ TSS and COD forecast provided by UV/Vis spectrophotometer. Therefore, it is important to test several statistical parameters to evaluate the PLS model [26].

Description of the Root-Vegetables Washing and Wastewater Treatment Systems
Four typical farms were chosen that represent different root vegetables, such as carrots and beetroots, and use different washing and wastewater treatment systems. All farms transferee the washing water used to the settling ponds (A-C) or settling tank (D). Settled water continuously is used for the first wash. Groundwater from a well (depth 80-150 m) is used for the final wash of vegetables (Figures 1-3). As a certain amount of excess water consistently generates, it is stored in the ponds. Excess water and sediments are spread out on the fields periodically according to legislation.
Appl. Sci. 2020, 10, 9072 3 of 15 relationships between water quality chemical parameters (suspended solids (SS), electric conductivity (EC), dissolved organic carbon (DOC), chemical oxygen demand CODCr, biochemical oxygen demand (BOD5)), and UV/Vis absorbance data. The aim was to establish a simple and reliable method for monitoring water quality which, sourced from the root vegetables washing and wastewater treatment processes within the on-farm packhouses. Statistical models can increase the importance and scope of the obtained spectral data [23]. Partial least squares regression (PLS) is a multidimensional statistical method used to predict standard values of aggregated water quality parameters such as total organic carbon (TOC), chemical oxygen demand (COD) and total suspended solids (TSS) [24,25]. The use of regression algorithms must comply with certain assumptions, i.e., the data should be independent and have a normal distribution. Predictive abilities of PLS calibration models were evaluated in relation to the determination factor (R 2 ) [26].
Different statistics were used to assess the PLS model and performance indicators were proposed to classify the developed models. It is essential to consider whether UV/Vis spectra can develop robust TSS and COD estimators, with statistical regression carried out on the off-line spectra in the collected samples. This affects the reliability of the in situ TSS and COD forecast provided by UV/Vis spectrophotometer. Therefore, it is important to test several statistical parameters to evaluate the PLS model [26].

Description of the Root-Vegetables Washing and Wastewater Treatment Systems
Four typical farms were chosen that represent different root vegetables, such as carrots and beetroots, and use different washing and wastewater treatment systems. All farms transferee the washing water used to the settling ponds (A-C) or settling tank (D). Settled water continuously is used for the first wash. Groundwater from a well (depth 80-150 m) is used for the final wash of vegetables (Figures 1-3). As a certain amount of excess water consistently generates, it is stored in the ponds. Excess water and sediments are spread out on the fields periodically according to legislation.

C
D washing system ( Figure 3)-carrots (10 t h −1 ): 1st wash-recirculating water from settling pond (70 m 3 t −1 ); polishing and 2nd wash-groundwater from well (7 m 3 t −1 ); excess flow from settling tank is transferred to the constructed system of three biological ponds. Firstly, excess water is discharged to the biological pond. Then it is supplied to surface flow wetland/macrophytes wetland. After that water is discharged into the 3 d pond and, finally, to the river basin ( Figure 3).

D
Settling tank Macrophytes wetland 3 d Pond Quality of the washing water and wastewater were monitored monthly over a 12-month period at the marked points (Figures 1-3). Results of the parameters monitored are presented in the same order as samples taken from marked points.

Sampling and Laboratory Analyses
The results of parameters monitored are presented in the same order as samples taken from marked sampling points (Figures 1-3). Temperature (T,°C) was monitored in-situ. In all cases, water samples were collected by the grab method and kept for one day at 4 °C prior to laboratory analyses. The collected samples were analyzed for the water quality parameters, such as color, suspended solids (SS), total dissolved solids (TDS), pH, electrical conductivity (EC), biological oxygen demand in 5 days (BOD5), chemical oxygen demand (CODCr and CODMn), dissolved organic carbon (DOC), Kjeldahl nitrogen (NKj), total phosphorus Ptot, and inorganic nitrogen N(NH4-N+ NO2-N + NO3-N). Deionized water samples were analyzed as standards.
All analytical methods are in accordance with standard methods [27] of the water and wastewater analysis, such as ISO 11923:1997, ASTM D5907-13, ISO 10523:2008, ISO 7888:1985, ISO   Figure 2. Flowchart of carrots washing and processing system C.
D washing system ( Figure 3)-carrots (10 t h −1 ): 1st wash-recirculating water from settling pond (70 m 3 t −1 ); polishing and 2nd wash-groundwater from well (7 m 3 t −1 ); excess flow from settling tank is transferred to the constructed system of three biological ponds. Firstly, excess water is discharged to the biological pond. Then it is supplied to surface flow wetland/macrophytes wetland. After that water is discharged into the 3 d pond and, finally, to the river basin ( Figure 3).

D
Settling tank Macrophytes wetland 3 d Pond Quality of the washing water and wastewater were monitored monthly over a 12-month period at the marked points (Figures 1-3). Results of the parameters monitored are presented in the same order as samples taken from marked points.

Sampling and Laboratory Analyses
The results of parameters monitored are presented in the same order as samples taken from marked sampling points (Figures 1-3). Temperature (T,°C) was monitored in-situ. In all cases, water samples were collected by the grab method and kept for one day at 4 °C prior to laboratory analyses. The collected samples were analyzed for the water quality parameters, such as color, suspended solids (SS), total dissolved solids (TDS), pH, electrical conductivity (EC), biological oxygen demand in 5 days (BOD5), chemical oxygen demand (CODCr and CODMn), dissolved organic carbon (DOC), Kjeldahl nitrogen (NKj), total phosphorus Ptot, and inorganic nitrogen N(NH4-N+ NO2-N + NO3-N). Deionized water samples were analyzed as standards.
D washing system ( Figure 3)-carrots (10 t h −1 ): 1st wash-recirculating water from settling pond (70 m 3 t −1 ); polishing and 2nd wash-groundwater from well (7 m 3 t −1 ); excess flow from settling tank is transferred to the constructed system of three biological ponds. Firstly, excess water is discharged to the biological pond. Then it is supplied to surface flow wetland/macrophytes wetland. After that water is discharged into the 3 d pond and, finally, to the river basin ( Figure 3).
Quality of the washing water and wastewater were monitored monthly over a 12-month period at the marked points (Figures 1-3). Results of the parameters monitored are presented in the same order as samples taken from marked points.

Sampling and Laboratory Analyses
The results of parameters monitored are presented in the same order as samples taken from marked sampling points (Figures 1-3). Temperature (T, • C) was monitored in-situ. In all cases, water samples were collected by the grab method and kept for one day at 4 • C prior to laboratory analyses. The collected samples were analyzed for the water quality parameters, such as color, suspended solids (SS), total dissolved solids (TDS), pH, electrical conductivity (EC), biological oxygen demand in 5 days (BOD 5 ), chemical oxygen demand (COD Cr and COD Mn ), dissolved organic carbon (DOC), Kjeldahl nitrogen (N Kj ), total phosphorus P tot , and inorganic nitrogen N(NH 4 -N + NO 2 -N + NO 3 -N). Deionized water samples were analyzed as standards.
All For the estimation of the uncertainties in the lab analysis, a triplicate analysis was performed, and the mean results, as well as standard deviation, were calculated. Measurements of the UV/Vis absorbance were performed with UV/Vis spectrophotometer Helios γ (Thermo Scientific, USA) in the range of 210 to 660 nm in a quartz cell of 10 mm path length after appropriate dilution with deionised water in order to overcome signal saturation in the UV region (absorbance values above 1 a.u.). In addition, blank samples of deionized water were measured to adjust the 100% transmittance signal. The data set was composed of 672 spectra of washing water and wastewater samples collected over a period of 12 months and of deionized water (control) at specific wavelengths of the UV/Vis spectra. Absorbance at A210 nm is used for roughly evaluation of nitrates, A225 nm-detergents, A240 nm-organic matter OM + suspended solids SS, A254 nm-dissolved organic carbon DOC, A320 nm-suspended solids SS, A390 nm-color, A660 nm-turbidity [2,6,28]. The dissolved organic carbon (DOC) quality was evaluated by specific ultraviolet absorbance (SUVA). SUVA 254 , which indicates the proportion of aromatic compounds in the DOC, was calculated according to equation [29]: Normally SUVA index of surface waters ranges from 1 to 6. As SUVA index increases, a larger concentration of aromatic compounds are in the solution when SUVA < 3, mainly hydrophilic low molecular weight compounds, are present. When SUVA > 4, organic matter is hydrophobic and high in aromatic compounds. Also, the bigger SUVA index can be resulted by iron, colloidal, or other impurities [7,30,31].
In order to use A254 nm as alternative parameters to estimate BOD 5 and COD, suspended solids have to be removed by filtering with a 0.45 µm pore size membrane filter. However, it is inconvenient, time-consuming and expensive. Chevakidagarn [2] studies of wastewater UV/Vis measurements showed a strong linear correlation between absorbance at A550 and the difference of filtered and non-filtered UV absorbance at A260 nm. Some filtered and unfiltered samples were measured at A240 nm, A320 nm, A254 nm and A660 nm and strong linear relationships were obtained (R 2 0.7-0.9). Therefore, measurements of unfiltered samples were made without the interference of suspended solids using the equation: Each sample's UV/Vis spectra were composed of three recorded spectra that were averaged after manual removal of erroneous spectra. Initial data were processed and the data arrangement of the data set was performed using the software Excel 2016 (Microsoft, Redmond, WA, USA).

Chemometrics Methods
The experimental data set was properly arranged in data matrices and auto-scaled (column mean-centered and scaled) to eliminate offsets and changes in measurement scales. Exploratory analysis of the autoscaled data set was conducted by principal component analysis PCA [22,[32][33][34][35][36][37]. The main objectives of the principal component analysis (PCA) are the transformation of the original data's coordinate system into a more relevant one and reduction of the original system dimension using a limited number of the principal components (PC's) which reflect the inherent structure of the data [38]. A principal component model approximates a given data matrix as follows [28]: where n represents the number of objects (samples), p-the number of variables (wavelengths), and d-the number of principal components (PC's), T-the scores matrix, L-the loadings matrix, LT-the transposed loadings matrix, and E-the residuals matrix.
The number of PC's was determined by the eigenvalue-one criterion analysis of the covariance matrix, according to which only the PC's with eigenvalues greater than one are considered relevant. The rotation of the components extracted was converted into six iterations by Varimax method with Kaiser normalization in order to extract optimal amount of PC's. Tests of PC's derived were performed by correlation and regression means. Two types of plots were obtained from the application of PCA. Loadings plots which describe and map the relationships between water quality parameters and the extracted principal components, and scores plots, which describe and map the mean samples UV/Vis spectra in the new axes defined by the principal components. One-vector scores plot to interpret the spectral PCA results was performed [28].
Multivariate calibration was performed using principal component regression analysis (PCR). The multiple linear regression (MLR) models assume that in addition to the p independent x-variables, a response variable y is measured, which can be explained as a linear combination of the x-variables. The prediction of the factor of interest y j can be described as [14]: where b i is the computed coefficient, x i represents the absorbance at each wavelength and e i,j , is the error. One drawback of PCA technique is that it captures only the characteristics of the X-vector or predictive variables. No importance is given of each predictive variable relation to the dependent or the target variable. Therefore, it is an unsupervised dimension reduction technique. Partial least square (PLS) allows achieving this balance and providing an alternative approach to PCA technique. Therefore, advanced modeling was performed by partial least square regression (PLS) statistical techniques [13,16,22,32,34,37,39,40]. PLS takes the X-Block (physico-chemical properties) and correlates them with Y-Block (spectral data) information to calculate latent variables (LV) for PLS calibration and prediction models. The calibration model is built by means of PLS, outlier tests and validation procedures to reach both high correlation quality and robustness. The result of the calibration procedure is a recovery function. The recovery function is obtained by plotting the reference measurements (actual targets) vs. the predicted values (estimated targets).
Chemometric models are used to correlate the concentrations of required determinants to spectra. Direct chemometric models are not used in wastewater. This is due to the lack of linear absorbance dependence of a single determinant required by Lambert-Beer law and is related to the parameters used for wastewater analysis. Replacement parameters are commonly used (e.g., COD and TSS), and from the specific properties of wastewater-a strong correlation between different parameters (e.g., COD is strongly correlated with filtered COD and TSS). Therefore, there is a need for intermediate chemometric models. They are based on estimating the calibration parameters from the calibration mixtures. They were using intermediate chemometric models. The validity of the Lambert-Beer law is no longer a prerequisite. Background modeling in the main component becomes feasible. Even systems with highly correlated spectra can also be used for multi-component calibrations [41].
The PLS model used takes into account in particular concentration-spectral dependencies and results in the most robust calibrations today. In the case of multifactorial calibration, the entire spectrum assessment is used. The calibration model is built using PLS, various validation procedures and outlier tests. The aim is to achieve high-quality correlations and quality of the analysis. The spectra and reference measurements must be standardized to perform PLS [41]. The procedure includes multi-factor calibration with PLS, multiple cross-validation and detection and elimination of outliers. The whole procedure is repeated until a good calibration result is obtained. The result of the calibration procedure is a recovery function. The recovery function is obtained by plotting actual cells using reference measurements against detailed cells by predicted values [21,42].
The wide range of wavelengths available allows great flexibility in selecting the best correlating wavelengths for the calibration function. This avoids cross-sensitivity. This is an advantage compared to systems that provide absorption measurements of only one or two wavelengths. The UV/Vis spectrometer provides global calibration for given parameters, but recalibration with local reference samples will improve performance and precision for quantitative measurements. For most cases, the reference data is a critical part of the whole calibration procedure. Therefore, to get good calibration results, it is necessary to guarantee the reference measurements' quality. They should relate to the reference analysis method, sample identity, sampling randomization, measuring range and sampling errors [42].
Quality assessment of the obtained results was performed by comparison of predicted values versus measured values, both for calibration and validation data sets. Model validation, cross-validation procedures, outlier tests and prediction errors RMSEC and RMSEP were performed to reach both high correlation quality and robustness [13,19,21,32].

Data of the Water Quality Monitoring
Descriptive statistics of the physicochemical parameters of the washing water quality are summarized in Table 1. Determination of the dissolved organic carbon (DOC) with TOC analyzer is the modern and rapid method. However, limit values of DOC are still under consideration by some researchers [2,5,7,10,30,45]. DOC can also be described by means of COD or BOD 5 . Results of this research showed the average fit of the linear regression models (Figure 4).
Appl. Sci. 2020, 10, 9072 7 of 15 reference data is a critical part of the whole calibration procedure. Therefore, to get good calibration results, it is necessary to guarantee the reference measurements' quality. They should relate to the reference analysis method, sample identity, sampling randomization, measuring range and sampling errors [42]. Quality assessment of the obtained results was performed by comparison of predicted values versus measured values, both for calibration and validation data sets. Model validation, crossvalidation procedures, outlier tests and prediction errors RMSEC and RMSEP were performed to reach both high correlation quality and robustness [13,19,21,32].

Data of the Water Quality Monitoring
Descriptive statistics of the physicochemical parameters of the washing water quality are summarized in Table 1. Determination of the dissolved organic carbon (DOC) with TOC analyzer is the modern and rapid method. However, limit values of DOC are still under consideration by some researchers [2,5,7,10,30,45]. DOC can also be described by means of COD or BOD5. Results of this research showed the average fit of the linear regression models (Figure 4). − It can be seen from the model, that the DOC limit values of the root vegetable washing water are in the range between 100 to 122 mg L −1 .
Contamination of municipal wastewater depends of its origin and load, particle sizes, rainfalls, storm events, and other conditions [2,3,5,9]    It can be seen from the model, that the DOC limit values of the root vegetable washing water are in the range between 100 to 122 mg L −1 .
Limits of wastewater discharge are based on a condition that no discharge should compromise the current state of a receiving water body [50]. They vary between countries. As shown in Table 1, the average SS concentration exceeded limit values for up to 40-90 times, BOD 5 -38-86 times, COD Cr -26 times, COD Mn -88 times, N Kj -1.8, P tot -2.5-19.6 and color-2 times.
The main reasons for the root vegetables washing water contamination and its' restriction of possible reuse in the packhouses are suspended solids (SS), biochemical oxygen demand (BOD 5 ), chemical oxygen demand (COD), Kjeldahl nitrogen N Kjt and total phosphorus P tot . Washing water of the carrots is more contaminated than that of beetroots, especially after polishing and processing. These findings show that the root vegetables' washing water is greatly contaminated and must be treated before discharge or further usage.
SUVA index was in the range of wastewater and natural waters. The third pond of the wastewater treatment system D increased sharply as DOC content was reduced greatly due to sedimentation and microbial decomposition. It is dominated by high molecular weight humic substances as it is typical for natural standing waters [7,30,31,51].
UV/Vis spectra were in accordance with typical absorbance spectra of wastewaters [2,6,28]. Two peaks at A240 and A320 indicate increased amounts of organic compounds and suspended solids. In addition, the effect of wastewater treatment of the pond in system D are obvious. It is obvious that absorbance decreases and the spectra of the water in the sampling spot 5 at the D site are close to the spectra of the deionized water ( Figure 5).

PCA and PCR Results of the Data Set
In order to prove the hypothesis, factor analysis was performed using the PCA method. Kaiser-Meyer-Olkin measure of sampling adequacy criterion (KMO) was 0.649, a significance level of Bartlett's test of sphericity Sig. < 0.0001. All measures of sampling adequacy (MSA) were >0.520. Extraction communalities were high and varied from 0.415 to 0.957. This indicates that the extracted components are a good representation of the variables (four extracted PCs explained 80.2% of the total data variance [52]). Principal component PC1 (39.5%) named as "Spectral data" was explained by UV/VIS spectral data excluding absorbance at the wavelength A 210 nm. Principal component PC2 (22.3%) named as "Dissolved solids" is explained by the N Kj , P tot , TDS, and color. Principal component PC3 (11.0%) named as "Organic substances" is explained by COD Cr , BOD 5 , and DOC. PC4 (7.5%), named as "Electrochemical properties" is explained by the pH, EC, and A210 nm, which is the measure of nitrates and is related to the electrochemical properties.
UV/Vis spectra were in accordance with typical absorbance spectra of wastewaters [2,6,28]. Two peaks at A240 and A320 indicate increased amounts of organic compounds and suspended solids. In addition, the effect of wastewater treatment of the pond in system D are obvious. It is obvious that absorbance decreases and the spectra of the water in the sampling spot 5 at the D site are close to the spectra of the deionized water ( Figure 5). Figure 5. UV/VIS fingerprints of the washing water and wastewater monitored at the sampling spots 1, 2, 3, 4 and 5 (Figures 1-3), C-deionized water.

PCA and PCR Results of the Data Set
In order to prove the hypothesis, factor analysis was performed using the PCA method. Kaiser-Meyer-Olkin measure of sampling adequacy criterion (KMO) was 0.649, a significance level of Bartlett's test of sphericity Sig. <0.0001. All measures of sampling adequacy (MSA) were >0.520. Extraction communalities were high and varied from 0.415 to 0.957. This indicates that the extracted components are a good representation of the variables (four extracted PCs explained 80.2% of the total data variance [52]). Principal component PC1 (39.5%) named as "Spectral data" was explained UV/Vis scores plot of spectral data revealed the highest correlation at 270-320 nm showing the highest DOC and SS contribution to PC1 while the lowest correlation was at 240 nm and decreased evenly from 320 nm to 660 nm. These findings confirm some speculations that DOC correlates better with UV/Vis spectra in comparison to organic matter (OM) while SS interferes with them due to scattering [5,21,30]. The right-side variables in Figure 6b are positively correlated one with another UV/Vis scores plot of spectral data revealed the highest correlation at 270-320 nm showing the highest DOC and SS contribution to PC1 while the lowest correlation was at 240 nm and decreased evenly from 320 nm to 660 nm. These findings confirm some speculations that DOC correlates better with UV/Vis spectra in comparison to organic matter (OM) while SS interferes with them due to scattering [5,21,30]. The right-side variables in Figure 6b are positively correlated one with another and with chemical analytes of the wastewater. Two clusters of the chemical properties and spectral data were established, while pH seems to be an outlier. Significant Pearson correlations (p < 0.05) were observed between PC1 (spectral data), PC3 (organic matter) (r = 0.86) and PC4 (electrochemical properties) (r = −0.79) while the correlation of PC2 (dissolved solids) was wake and statistically insignificant. The Chi-Square test (Asymp. Sig. 0.102 > 0.05) confirmed that our model consisting of four PCs represents a sufficient description of the chemical data. One-Sample Kolmogorov-Smirnov Test (Asymp. Sig < 0.05) confirmed the hypothesis about compatibility between normal empirical and theoretical distributions of PCs. Monte Carlo significance p > 0.05 let to reject H0 hypothesis and conclude that the PCs' distributions are compatible with normal distribution. Data of the Cook's test < 4 showed that there are not outliers in the model. Mahalanobis test identified outlier sample only of deionized water, which was used as background for spectral analysis.

PLS Regression and Prediction Results of the Data Set
Partial least square regression model (PLS2) of the mean-centered spectral data with three latent variables (LV) explained 85% of the total variance while LV1 explained 32% after calibration, cross-validation, and outliers' removal. The best fitting to the model showed data at absorbance of A320 ( Figure 7). PLS2 regression model confirmed PCA findings of two uncorrelated clusters ( Figure 5b) and showed the great influence of suspended solids (SS), dissolved solids (TDS and EC), and organic matter (BOD 5 , COD Cr , and DOC) while the influence of pH, N Kj , P tot and color was negligible. In Table 2, the limit values of the contamination of the root vegetables washing water and wastewater and absorbance spectra of the UV/Vis predicted by PLS regression model are presented.
As A320 data have the greatest impact on the model (Figure 5a), it was proposed to use this parameter to monitor water quality in the plants for washing the root vegetables.

PLS Regression and Prediction Results of the Data Set
Partial least square regression model (PLS2) of the mean-centered spectral data with three latent variables (LV) explained 85% of the total variance while LV1 explained 32% after calibration, crossvalidation, and outliers' removal. The best fitting to the model showed data at absorbance of A320 ( Figure 7). PLS2 regression model confirmed PCA findings of two uncorrelated clusters ( Figure 5b) and showed the great influence of suspended solids (SS), dissolved solids (TDS and EC), and organic matter (BOD5, CODCr, and DOC) while the influence of pH, NKj, Ptot and color was negligible. In Table

Conclusions
Chemical analysis of the investigated water samples showed that washing water of root vegetables exceeded the limit values by 10-50 times for suspended solids, 10-25 times for organic matter, 3-5 times for dissolved organic carbon, 3-7 times for total nitrogen, 3-12 times for total phosphorus, and by 3-10 times for microorganisms. Therefore, washing water quality monitoring and proper wastewater treatment is of great importance.
UV/Vis absorbance spectra indicated increased amounts of organic compounds and suspended solids and effects of the water treatment in the ponds system. PCA and PLS regression methods have been applied to analyze relationships between chemical properties and UV/Vis spectral data. The PCA and PLS modeling results showed that UV/Vis at 320 nm could be a sufficient measure for the water pollution with suspended, dissolved solids and organic matter monitoring. According to these materials' chemical limit values, limit value of the absorption at A320 nm is 7.39-7.41 cm −1 . This limit value will be validated at the feature research.