Reuse of Agriculture Drainage Water in a Mixed Land-Use Watershed

Water demand in Korea has triggered the need for fresh water to be used for agriculture. Agricultural drainage water (ADW) is a way of coping with the growing demand for fresh water for agriculture. In this study, a water quality model (WQM), and an algorithm were used in order to determine the water quality and optimize the water reuse quantity in the Osan stream drain, South Korea. The water quality associated with the drain was stimulated using the QUAL2Kw model and the uncertainty analysis and sensitivity analysis with the use of Monte Carlos Simulation was performed to determine the performance of the WQM. Jaya algorithm technology was used as an optimization tool to find optimal ADW reuse quantities at particular withdrawal points. For calibration and validation, the model was applied twice for both summer and winter seasons. The results show that the reuse quantities represent 77.2% and 49.8% of the available ADW in the study area for summer and winter, respectively, representing 49.1% and 54.5% of seasonal canal delivery. The utilization of the simulation-optimization model is usually well suited for decision support leading to near-optimum reuse assortment of ADW for irrigation.


Introduction
In the last century, water use increased worldwide at more than twice the population rate.In arid regions, water use increasingly reaches the limit where water services can be delivered sustainably.About one-fifth of the total population, live in shortage regions.Another 1.6 billion people, or about a quarter of the total population, face financial water shortages [1].About 20 hectares of agricultural land is irrigated by wastewater, and about 10% of the world's population consumes food products from farmland [2,3].
South Korea faces many water-related challenges, including an increasing population and related industrial and agricultural activities that will raise the demand for water to levels reaching the limits of the accessible supply.In 2007, an aggregate volume of 15.9 billion m 3 water was utilized with regard to agriculture in Korea, which involved around 48% of the aggregate yearly water use in Korea [4,5].In line with the Ministry of Environment [6], the total wastewater treatment limit in Korea is usually 6.6 billion m 3 for each year, which often comes to 20% associated with the aggregate yearly water use.An ongoing nationwide overview on the future water request and source of Korea announced an anticipated deficiency of more than 4.4 billion m 3 of water by 2020 [4,7].
The Korean Government has been implementing wastewater reclamation projects for agricultural land, mainly for paddy fields, since 2005 and the reused amount in 2011 reached 45 million m 3 .Most wastewater reuse in agriculture is indirect, which means that the wastewater is disposed of into streams and the diluted wastewater is used for irrigation [8].Therefore, farms near Wastewater Treatment Plants (WWTPs), which are significantly affected by the effluent quality of a WWTP, particularly in the dry season [9], are likely to be expanded.In both direct and indirect wastewater reuse, disposal levels of WWTPs affect irrigation water quality.In the case of reuse of indirect wastewater, the hydrological conditions under which the treated wastewater is diluted also affect irrigation water in addition to the level of disposal.Therefore, before any analysis of the effects of indirect wastewater reuse, it is important to address how irrigation water is affected by WWTP effluent on irrigation water.Optimizing the water quality and quantity of reused drainage water is, therefore, the primary issue of any strategic plan for the effective management of water in South Korea.
Through the years, SWQMs have increasingly already been recognized to support decisionmaking and water quality management.Recently several SWQMs were developed to evaluate and predict the quality of water in streams.QUAL2Kw is a model with different parameters for flexible and applicable water quantity and quality simulation in different settings [10][11][12][13].The complex relationship between waste loads and the resulting water quality is best described in the QUAL2Kw model in the receiving water bodies [14,15].
Regardless of the wealth of studies on wastewater irrigation, few studies have been conducted on the influence of WWTP effluents on irrigation water in situations where indirect wastewater reuse accounts for the majority of wastewater reuse.Other studies have been done with QUAL2Kw and Genetic algorithm, but none have assessed QUAL2Kw and Jaya algorithm which is a new algorithm.
Water quality for irrigation will have declined, due to mixed land use along the Osan stream.This study aims to simulate ADW reuse by integrating a WQM (QULA2Kw) in the mixed land use watershed.The research used the QUAL2Kw model to study the water quality of the Osan drain and its environs.Uncertainty and sensitivity analyses have been carried out to ensure the reliability of the ADW quality and quantity simulation using QUAL2Kw.The Jaya Algorithm technique was then applied as an optimization tool to determine the optimal reuse quantities of ADW.

Study Area
The QUAL2Kw model was applied to the Osan River watershed situated in central Korea with latitude and longitude which range from N 37 • 05 56" to N 37 • 14 05" and from E 127 • 01 29" to E 127 • 09 38", respectively (Figure 1).The Osan River watershed is suitable for the study watershed since it includes the Osan WWTP, which has a treatment capacity of 140,000 m 3 d −1 , downstream from the main water body.A pumping station that irrigates diluted effluent discharge to paddy fields is situated approximately 1 km downstream from the WWTP.In this study, the main use irrigation water source is drainage water reuse, hence Tap-dong WWTP water is mixed with fresh Osan stream irrigation water.The canals are not labeled and the flow rates are not measured.The mixed water is used by the farmers in the Osan watershed for farming.Land use data for 2009 of the watershed is available at the Ministry of Environment (MOE).It is characterized, as follows-32.3%forest, 26.1% built-up, 18.3% paddy, 8.1% bare soil, 6.2% upland, and 3.5% water.The average effluent water quality from the WWTPs at the time of 2011 was presented as 10.1 mg L −1 for total nitrogen (T-N), 0.8 mg L −1 for total phosphorus (T-P), and 155 MPN 100 mL −1 for total coliform.Additionally, it satisfied the effluent standard for WWTPs in Korea.On the other hand, it was reported that some agricultural land remained irrigated with untreated wastewater [14]   In this study, the average monthly records of water quality for January 2015 and July 2016 were collected for 6 locations from the mainstream.The average monthly records of dissolved oxygen (DO), temperature, pH, biological oxygen demand (BOD), total dissolved solutes (TDS), total nitrogen (TN), chemical oxygen demand (COD), total phosphorous (TP) and total suspended solids (TSS) and so on were measured according to [16].Water quality variables for the 6 monitoring locations, as shown in Table 1 below.In South Korea, the Ministry of Environment [17] proclaimed standards for water quality for treated wastewater are based their specific purpose of use.In this study, the average monthly records of water quality for January 2015 and July 2016 were collected for 6 locations from the mainstream.The average monthly records of dissolved oxygen (DO), temperature, pH, biological oxygen demand (BOD), total dissolved solutes (TDS), total nitrogen (TN), chemical oxygen demand (COD), total phosphorous (TP) and total suspended solids (TSS) and so on were measured according to [16].Water quality variables for the 6 monitoring locations, as shown in Table 1 below.In South Korea, the Ministry of Environment [17] proclaimed standards for water quality for treated wastewater are based their specific purpose of use.

Methods
In this study, the simulation optimization model was developed to achieve near-optimal mixing ratios between ADW and fresh (canal) water at the reused sites.The Osan stream drain was simulated by QUAL2Kw.The SWQM model yields the water quality parameters which include biological oxygen demand (BOD), dissolved oxygen (DO), and total suspended solids (TSS) along the drainage network.Uncertainty and sensitivity analyses, based on MCS were used to evaluate the accuracy of the model and Jaya Algorithm to optimize the total reuse of the proposed reuse sites.In the simulation optimization framework, the SWQM and the optimization methods are combined.

Water Quality Modeling
The SWQM is a one-dimensional and steady stream flow/water quality model.The SWQM can simulate up to 4 user-specified water quality components.Simulated here were two non-conservative constituents (BOD and DO concentrations) and two conservative constituents (TSS and TN concentrations).QUAL2Kw is a stationary and one-dimensional stream modeling designed to represent a modernized version of the U.S. EPA standard river water quality model: QUAL2E [18] to QUAL2Kw.The framework also includes a number of new features that can be applied to shallow and upland streams, as well as more current science.QUAL2Kw simulates the transport and fate of conventional pollutants.The model represents the river as a one-dimensional channel with a steady, non-uniform flow and simulates both point and non-point pollutants.
The QUAL2Kw model simulates the transport and various constituents, such as temperature, wind speed, dew point, temperature, dissolved oxygen, phytoplankton, carbonaceous biochemical oxygen demand, phosphorus, and nitrogen etc.It also simulates several other constituents not typically included in the generally available software, and it is free software.The model simulates pathogenic bacteria, pH, alkalinity, inorganic suspended solids, and ground algae in particular.The inclusion of bottom algae is essential for the simulation of shallow streams so that these algae have a new feature of the nitrogen and phosphorus nutrient variable stoichiometry.
The QUAL2Kw model has other functions that distinguish it from all other models.In other words, oxygen and nutrient flows are calculated on the basis of particulate organic matter, sediment reactions and soluble forms concentrations in overlying waters.The hyporheic zone is modeled, therefore it is the area under the stream bed where water flows through spaces between rocks and cobbles.It is also a necessary feature to simulate shallow streams.
QUAL2Kw has been installed in Microsoft Excel.It is programmed in Visual Basic for Applications (VBA).Input data, model execution, and output viewing was a graphical user interface that was performed in excel.In a model run, a compiled FORTRAN 95 program run by the Excel VBA program performs numerical integration.Constituent concentration (c i ) of mass balance in the water column of a reach is written as where C i is the concentration component of water quality, Q i is flow, V i is volume, E i is dispersion coefficient and W i is the outer load component of reach i. S i is the sinks and sources of the component, which are due to a wide number of transformation mechanisms and reactions within reach i. Q out , i is the flow abstraction from reach i.
The source/sink term (S i in Equation ( 1)) requires a large number of parameters for each state variable (e.g., maximum growth rate of bottom algae) to be specific.The user may select which parameters are held at constant values and which are to be optimized by Jaya Algorithm.The QUAL2Kw Documentation and User's Manual describes a comprehensive discussion of the model theory [19] and provides a detailed description of the model [20].
The Osan Stream drain was divided into 5 reaches with different lengths and subdivided into a total of 47 elements of 0.350 km length (Figure 2).The main criteria for this division were the similarity of hydraulic characteristics, including slope and width, in the reach of a canal.Hydraulic characteristics were determined by measuring geometries and water velocities at each sampling location.Manning's equation was used to determine water velocity and the depth of the stream.For all reaches, a Manning roughness coefficient of 0.04 was used [21].The input parameters for water quality included in the model are total coliform, TSS, TDS, DO, BOD 5 , COD NH 4 -N, NO 3 -N and NO 2 -N.Phytoplankton, detritus, and pathogens were not measured and were left empty.Algae and bottom oxygen demand were assumed to be 40% and 100%, respectively.The thickness of the sediment/hyporheic zone was assumed to be 10 cm and the calculation time step was set at 11.25 min to prevent model instability.Five days was the number of days of simulation period while Euler's method was done by integration.characteristics were determined by measuring geometries and water velocities at each sampling location.Manning's equation was used to determine water velocity and the depth of the stream.For all reaches, a Manning roughness coefficient of 0.04 was used [21].The input parameters for water quality included in the model are total coliform, TSS, TDS, DO, BOD5, COD NH4-N, NO3-N and NO2-N.Phytoplankton, detritus, and pathogens were not measured and were left empty.Algae and bottom oxygen demand were assumed to be 40% and 100%, respectively.The thickness of the sediment/hyporheic zone was assumed to be 10 cm and the calculation time step was set at 11.25 min to prevent model instability.Five days was the number of days of simulation period while Euler's method was done by integration.The QAUL2Kw model was designed for a population size of 100 with 100 generations of evolution [22].The calibrated parameters values will be the default values for now or same as presented in [23].For calibration, parameters for the Osan drain QUAL2Kw model using data sets from January 2015 to July 2015.The calibrated model was used for validation with a completely different data set from January 2016 and July 2016 (Water Information system) to aid to the smooth running of the model and to predict the quality and quantity of water under several conditions.

Uncertainty Analysis
An uncertainty analysis was executed to assess the uncertainties related to input variable values (BOD, DO, TSS and TN) through QUAL2Kw in order to obtain the statistical measurements of the distributions of the resulting outputs.Random outputs of the QUAL2Kw model within the Osan drain resulting in the magnification of errors in the outputs of the model.It is important to understand the uncertainty in the ability of the QUAL2Kw model hence to predict the drain's behavior and in managing reuse activities to validate the simulated drain.MCS was used to perform Uncertainty analysis in QUAL2Kw model.Monte Carlos Simulation (MCS) is a technique that handles non-normal distributions, complex algorithms, and correlations between input factors for the model in question.This procedure randomly samples input variables from probability distributions without or with correlation, hence our model gives different outputs for each selected variable [13,24,25].The various outputs have been analyzed statistically to assess the uncertainty of predicting the model [13,26,27].The MCS will be perform using YASAIw add-in, is an open source and freely available framework for MCS in Excel [19].
YASAIw is an amendment to the original YASAI add-in developed at Rutgers University [28].The basic YASAI functionality and theory is detailed [28,29].YASAI provides Excel functions that return random numbers with specified distributions.They will generally return different, randomly chosen values every time you recalculate a spreadsheet [29].YASAI is implemented as a single file in the Visual Basic for Applications (VBA) language.Any user on any computer equipped with a reasonably recent version of Excel should be able to use it, provided they are allowed to run Excel macros [30].YASAI provides only six distributions with no graphics and the graphics can easily be added with great difficulty.In this study, an uncertainty analysis was performed by using the QUAL2Kw model with YASAIw tool to evaluate the input variables in the Osan drain to obtain uncertainty outputs that confirms the QUAL2Kw water quality model.For each output, the model can produce histograms, frequencies, error bars, and probability density functions.

Sensitivity Analysis
Sensitivity analysis examines the importance of imprecision in model inputs in modeling or decision-making.In this study, outputs from Monte Carlos Simulation (MCS) produces sensitivity analysis which estimates the correlation among the predicted variables and the chosen input variable hence the contribution to the variance of the predicted variable from the chosen variables in the MCS.Four water quality outputs (BOD, DO, TN and TSS) at one major location (Osan Stream 2-1, 5.15 km) was chosen for and sensitivity analyses, because of the presence of a Wastewater treatment plant (WWTP) for both summer and winter respectively.A total of 1000 MCSs were performed to generate 1000 outputs for the four water quality variables.The Monte Carlos Simulation is implemented by simulating the QUAL2Kw model to identify parameters that have the greatest impact on water quality outputs.A total of 66 input variables were chosen, which include BOD, DO, TN and TSS from headwater and QUAL2Kw model calibrated rate constants.In this study, sensitivity analysis produces the contribution to the variance and Spearman's correlation coefficients from each output.Spearman's rank correlation is a nonparametric measure of the correlation that uses the rank of observations in its calculation, rather than the original numeric values.It measures the monotonic relationship between two variables.Table 2 presents the mean water quality for summer and winter 2016 at proposed reused stations.Table 3 presents the kinetics parameters adapted from the literature and used in the model calibration process [18,22].

Jaya Algorithm
Jaya algorithm is a new optimization algorithm which has no algorithm-specific controlling parameter for solving constrained and unconstrained optimization problems.In contrast to other population-based heuristic algorithms, only the two ordinary control parameters, population size (m) and generation number (Gn), are included [31].The basic Jaya algorithm has only one phase of the above concept and is, therefore, a simple technique of optimization.Using binary coding, parameters can represent their values with a high degree of accuracy [31].
To further explain the Jaya Algorithm, let Q j,k,i is the objective function.At any iteration i, assume that there are 'm' number of design variables (i.e., j = 1, 2, . . ., m), 'n' number of candidate solutions (i.e., population size, k = 1, 2, . . ., n).Let Q j,best,i and Q j,worst,i be the best and worst candidate solutions respectively.If Q j,k,i stands for the value of the variable j th for the kth candidate during the i th iteration.The formula is modified and written as the following Equation (2).
where, Q j,best,i stands for the j th variable for the best candidate and Q j,worst,i stands for the j th variable for the worst candidate.Q j,k,i is the updated value of Q j,k,i and r 1,j,i and r 2,j,i are the two random numbers in the range [0, 1].This "r 1,j,i (Q j,best,i − |Q j,k,i |)" shows the tendency to approach to the best solution and the term "−r 2,j,i (Q j,worst,i −|Q j,k,i |)" shows the tendency to avoid the worst solution.Q j,k,i is accepted if it's better than the previous candidate or it's maintained [31].
The Jaya algorithm does need any algorithm-specific parameters, which requires extensive calibration before conducting the actual computational experiments, if not done correctly, leads unavoidable and unwanted convergence.The controlling of the algorithm-specific parameters is not as easy as it looks.Also, controlling the parameters in each iteration is most of the time difficult and time-consuming.This entire process does not belong to the Jaya algorithm.The main reason is the victorious nature of the Jaya algorithm, which makes it more powerful than any other algorithm.

Jaya Algorithm Application
The Jaya algorithm is used as an optimization tool to maximize the ADW reuse quantities at the proposed reuse sites unless the standards for the use of mixed water in irrigation have been violated.In the simulation-optimization framework, the Jaya algorithm, and the QUAL2Kw model will be combined to achieve the optimum use of ADW in the Osan drain.
The aim is to maximize the removal of ADW from the drain to the mixing sites on the canal (Figure 3).
Subject to where Q t is the total ADW reuse drain discharges.TSS c = TDS concentration (mg/L) for the irrigation water downstream from the mixing location in a canal c; TSS c;max = maximum permissible concentration of TSS (mg/L) for irrigation water; BOD c = BOD concentration (mg/L) for the irrigation water downstream from the mixing location; BOD c;max = maximum permissible BOD concentration (mg/L) for irrigation water; DOc = concentration of DO (mg/L) for the irrigation water downstream from the mixing location; and DO c;min = minimum permissible concentration of DO (mg/L) for irrigation water.The Osan stream drain is composed of many individual outlets from the stream made by the farmers.To optimize this drain we assumed five proposed reuse locations due to the availability of water quality parameters and due to the unavailability of data on flow and water quality parameters, the flow rate of the stream for the past 8 years was used to extrapolate the flow rate at various locations.The water quality parameters were extrapolated from the water quality data at the monitoring points.

Results and Discussion
This optimization model was applied to the 2016 summer and winter Osan for the maximum and minimum reuse of the ADW.Water quality and quantity data were therefore used in the winter (January 2016) and summer (July 2016) to model the drainage water quality and quantity of the Osan drain.QUAL2Kw model calibrated parameters were not changed.Our choice of the month January to represent the winter period was due to the low rainfall index and the great influence of the contribution of the diffuse load for the water quality.Similarly, our choice of the month of July to represent the summer period was due to the higher rainfall indices.In two subsequent Jaya Algorithm experiments, QUAL2Kw models were used for summer and winter.Given that agricultural water consumption accounts for more than 47% of total water consumption and the trend towards reuse in Korea is generally increasing [25,32], it is important to categorize the culture of reuse of wastewater and to introduce representative cases at the local level for a safe and sustainable environment.

Results and Discussion
This optimization model was applied to the 2016 summer and winter Osan for the maximum and minimum reuse of the ADW.Water quality and quantity data were therefore used in the winter (January 2016) and summer (July 2016) to model the drainage water quality and quantity of the Osan drain.QUAL2Kw model calibrated parameters were not changed.Our choice of the month January to represent the winter period was due to the low rainfall index and the great influence of the contribution of the diffuse load for the water quality.Similarly, our choice of the month of July to represent the summer period was due to the higher rainfall indices.In two subsequent Jaya Algorithm experiments, QUAL2Kw models were used for summer and winter.Given that agricultural water consumption accounts for more than 47% of total water consumption and the trend towards reuse in Korea is generally increasing [25,32], it is important to categorize the culture of reuse of wastewater and to introduce representative cases at the local level for a safe and sustainable environment.

Water Quality Simulation
The simulation model was executed twice by using completely different data sets from summer and winter 2016 to represent seasons with maximum and minimum ADW reuse requirements.In order to simulate this drainage system, it is divided into five reaches on the basis of its hydraulic characteristics.The model predictions along the Osan stream drain were plotted and summarized in Figure 4 below.
wastewater from the Tap-dong WWTP add high organic and inorganic materials resulted low DO.The water quality standards for agricultural water are classified into four classes of water quality [33].
Furthermore, the results showed that the TSS and T-N values in many drain segments meet the standards (i.e., below 1000 and 15 mg/L respectively).The TSS concentration for summer increases sharply after 1 km and decreases downstream hence the presence of pathogens available in the drain while in winter it's in a steady state from headwater to downstream.T-N concentrations increase along the drain, so the impact of fertilizer application in the agricultural drainage basin areas can be attributed.

Uncertainty Analysis Outcome
The results suggested that input data were efficiently simulated in the Osan stream drain.The frequency histograms and error bars from the uncertainty analysis show that the data are clumped around the mean for all the simulated variables at the Osan Stream 2-1, 5.15 km for winter and summer respectively (Figure 5).At Osan Stream 2-1, the mean for the four variables examined suggests that both summer and winter the BOD, DO, TN and TSS predictions are more reliable, as they have a small standard deviation bar.Corresponding values are 0.026, 8.1, 9.3 and 18.35 for summer, and 3.15, 12.37, 4.2 and 7.89 for winter, respectively.The results showed that MCS simulated The direct flow of ADW into these drains is simulated by the distributed flow along the Osan stream.For winter and summer simulations, this value was found to be 1.34 m 3 /s and 1.85 m 3 /s.TSS, TN, BOD and DO along the drainage canal ranged from 0.92, 1.22, 0.75 and 5.26 mg/L for winter and 3.30, 9.46, 2.08 and 8.0 mg/L for summer, respectively, in the calibration (January 2015 to December 2015) to 4.20, 7.88, 1.12 and 12.37 mg/L for winter and 17.88, 9.34, 0.03 and 8.28 mg/L for summer, respectively, in the validation (January 2016 to December 2016).
DO is always greater than 2 mg/L, so it meets the 4th grade criterion for streams, and BOD meets the criteria of 8 mg/L.In the case of T-N and T-P, there is no standard for the lower level.Therefore, it is compared with the water quality recommendation standard [5].
Along the drain, BOD ranged from 0.10 to 0.4 mg/L for summer and from 0.84 to 2.20 mg/L for winter, as shown in Figure 4. Similarly, the DO ranged from 7.73 to 8.28 mg/L in summer and from 12.50 to 12.69 mg/L winter, for reusing drainage water in irrigation in South Korea [15] hence wastewater from the Tap-dong WWTP add high organic and inorganic materials resulted low DO.The water quality standards for agricultural water are classified into four classes of water quality [33].
Furthermore, the results showed that the TSS and T-N values in many drain segments meet the standards (i.e., below 1000 and 15 mg/L respectively).The TSS concentration for summer increases sharply after 1 km and decreases downstream hence the presence of pathogens available in the drain while in winter it's in a steady state from headwater to downstream.T-N concentrations increase along the drain, so the impact of fertilizer application in the agricultural drainage basin areas can be attributed.

Uncertainty Analysis Outcome
The results suggested that input data were efficiently simulated in the Osan stream drain.The frequency histograms and error bars from the uncertainty analysis show that the data are clumped around the mean for all the simulated variables at the Osan Stream 2-1, 5.15 km for winter and summer respectively (Figure 5).At Osan Stream 2-1, the mean for the four variables examined suggests that both summer and winter the BOD, DO, TN and TSS predictions are more reliable, as they have a small standard deviation bar.Corresponding values are 0.026, 8.1, 9.3 and 18.35 for summer, and 3.15, 12.37, 4.2 and 7.89 for winter, respectively.The results showed that MCS simulated mean for BOD, DO, TN and TSS increased downstream hence the presence of the wastewater pumping station is located downstream.The results indicated that the likelihood of variables in water quality, namely BOD, DO, TN, and TSS, deviating from observed values was low.The QUAL2Kw model is used in the study and can be used as an intervention analysis for DO, BOD, TN, and TSS.The model can, therefore, be used effectively to manage water quality.

Uncertainty Analysis Outcome
The results suggested that input data were efficiently simulated in the Osan stream drain.The frequency histograms and error bars from the uncertainty analysis show that the data are clumped around the mean for all the simulated variables at the Osan Stream 2-1, 5.15 km for winter and summer respectively (Figure 5).At Osan Stream 2-1, the mean for the four variables examined suggests that both summer and winter the BOD, DO, TN and TSS predictions are more reliable, as they have a small standard deviation bar.Corresponding values are 0.026, 8.1, 9.3 and 18.35 for summer, and 3.15, 12.37, 4.2 and 7.89 for winter, respectively.The results showed that MCS simulated mean for BOD, DO, TN and TSS increased downstream hence the presence of the wastewater pumping station is located downstream.The results indicated that the likelihood of variables in water quality, namely BOD, DO, TN, and TSS, deviating from observed values was low.The QUAL2Kw model is used in the study and can be used as an intervention analysis for DO, BOD, TN, and TSS.The model can, therefore, be used effectively to manage water quality.

Sensitivity Analysis Results
The results of the sensitivity analysis showed that the contribution to variance parameter makes up about 95% of the BOD, DO, TN and TSS variance at Osan Stream 2-1 both in summer and in winter.Table 4 showed that Spearman's Rho values of −0.93, −0.85, −0.92 and 0.97 for winter and −91.1, 0.80, −0.71 and −0.75 for summer were the highest input for BOD, DO, TN and TSS outputs at the Osan Stream 2-1 reuse station.Their contribution to variance were 81.6%, 71.5%, 80.7%, and 88.1% for winter and 76.6%, 64.6%, 55.2% and 48.6% for summer.This is because of the effect of the drain, which releases its pollution loads into the stream.

Sensitivity Analysis Results
The results of the sensitivity analysis showed that the contribution to variance parameter makes up about 95% of the BOD, DO, TN and TSS variance at Osan Stream 2-1 both in summer and in winter.Table 4 showed that Spearman's Rho values of −0.93, −0.85, −0.92 and 0.97 for winter and −91.1, 0.80, −0.71 and −0.75 for summer were the highest input for BOD, DO, TN and TSS outputs at the Osan Stream 2-1 reuse station.Their contribution to variance were 81.6%, 71.5%, 80.7%, and 88.1% for winter and 76.6%, 64.6%, 55.2% and 48.6% for summer.This is because of the effect of the drain, which releases its pollution loads into the stream.

Optimization of ADW Reuse
The simulation optimization model using Jaya Algorithm was run independently for summer and winter, representing seasons with maximum and minimum ADW reuse requirements.These two experiments in the Jaya algorithm identified the water quality, the mixing ratio and optimal reuse withdrawals between fresh water and drainage at each mixing point for the experiments.Total reuse quantities were 45.73 and 6.80 m 3 /s for summer and winter, respectively (Table 5).The water quantity in the irrigation canals upstream from the mixing points for summer and winter was 47.35 and 5.75 m 3 /s so that the average mixing ratio between drainage and fresh water was 2.0 in summer and 2.2 in winter.In addition, the criteria for concentrations of BOD, DO, TSS and TN in canal water after mixing at all mixing points met the criteria for irrigation.
The water requirements for irrigation downstream of the reused points were 93.08 m 3 /s for summer and 12.55 m 3 /s for winter (Figure 6).The reused amounts of ADW, therefore, represent 49.1% and 54.2% of the summer and winter requirements for fresh water for irrigation.In addition, the total demand for irrigation water for summer and winter was 59.2 and 13.65 m 3 /s, with recycled water volumes accounting for 77.2% and 49.8% of the available ADW for summer and winter in the study area.Therefore, the amount of ADW released is higher in summer than in winter.The optimization results showed that the Jaya algorithm was used successfully in the management of ADW reuse.Therefore, the amount of ADW released is higher in summer than in winter.The optimization results showed that the Jaya algorithm was used successfully in the management of ADW reuse.

Optimization of ADW Reuse
The simulation optimization model using Jaya Algorithm was run independently for summer and winter, representing seasons with maximum and minimum ADW reuse requirements.These two experiments in the Jaya algorithm identified the water quality, the mixing ratio and optimal reuse withdrawals between fresh water and drainage at each mixing point for the experiments.Total reuse quantities were 45.73 and 6.80 m 3 /s for summer and winter, respectively (Table 5).The water quantity in the irrigation canals upstream from the mixing points for summer and winter was 47.35 and 5.75 m 3 /s so that the average mixing ratio between drainage and fresh water was 2.0 in summer and 2.2 in winter.In addition, the criteria for concentrations of BOD, DO, TSS and TN in canal water after mixing at all mixing points met the criteria for irrigation.
The water requirements for irrigation downstream of the reused points were 93.08 m 3 /s for summer and 12.55 m 3 /s for winter (Figure 6).The reused amounts of ADW, therefore, represent 49.1% and 54.2% of the summer and winter requirements for fresh water for irrigation.In addition, the total demand for irrigation water for summer and winter was 59.2 and 13.65 m 3 /s, with recycled water volumes accounting for 77.2% and 49.8% of the available ADW for summer and winter in the study area.Therefore, the amount of ADW released is higher in summer than in winter.The optimization results showed that the Jaya algorithm was used successfully in the management of ADW reuse.Therefore, the amount of ADW released is higher in summer than in winter.The optimization results showed that the Jaya algorithm was used successfully in the management of ADW reuse.This simulation-optimization model can be applied to other drainage basins with similar characterizes of ADW in South Korea hence it may aid to estimate the permissible waste load to be released into a drain.Moreover, to attain suitable reuse water for irrigation and to manage the total pollution load in the drain.
Figure 7 shows how the objective function decreases at each iteration with Jaya algorithm.This provides evidence in favor of optimality of the Jaya Algorithm as the graph converges to the minimum with the given accuracy.This simulation-optimization model can be applied to other drainage basins with similar characterizes of ADW in South Korea hence it may aid to estimate the permissible waste load to be released into a drain.Moreover, to attain suitable reuse water for irrigation and to manage the total pollution load in the drain.
Figure 7 shows how the objective function decreases at each iteration with Jaya algorithm.This provides evidence in favor of optimality of the Jaya Algorithm as the graph converges to the minimum with the given accuracy.

Conclusions
A simulation optimization methodology was developed and applied to Osan Stream, South Korea, for the reuse of ADW in irrigation.This study simulates the water quality of the Osan stream drain using surface water quality model, QUAL2Kw and an optimization algorithm, Jaya algorithm.Jaya algorithm was used to determine the optimal reuse of ADW at five proposed reuse sites.
Results show that ADW from the Osan stream drainage basin can be used directly for irrigation at the five reuse locations with only a few violations of Korean standards and FAO guidelines.Thus, mixing wastewater with freshwater increases its suitability for irrigation and makes it suitable for irrigation.The QUAL2Kw model was calibrated and applied using data obtained during the period 2014-2016.Water quality results show accepted levels of BOD, DO, TN and TSS.Uncertainty and sensitivity analyses were applied to QUAL2Kw to verify the input variable.Uncertainty analysis predictions are reliable for both summer and winter while the sensitivity analysis showed that the contribution to variance parameter makes up about 95% of the sensitivity analysis showed the variance contributions of the parameters accounted for up to 95% of the variance.The optimization results showed that the reuse quantities can reach 49.1% and 54.5% ADW available in the summer and winter drain.These quantities represent 80% and 54.6% of the demand for fresh irrigation water in the drainage basin for summer and winter.This approach can be extended to irrigated systems to other drain irrigation systems with similar characteristics.Finally, the simulation-optimization model may improve the decision making on agricultural drainage water for irrigation for mixed watersheds.
and limited research has been carried out into various domestic examples of wastewater irrigation.The study watershed has a temperate monsoon climate.The average temperature is approximately 12 • C, and the yearly average precipitation is approximately 1312 mm.Two-thirds of annual rainfall in the summer season lasts from June to August.The study watershed furthermore has an area of 98.3 km 2 , drainage density of 0.52 km −2 , main river length of 16.49 km, and drainage slope of 10.21%.At the Tap-Dong streamflow gauge station, which is managed by the Ministry of Land, Infrastructure, and Transport (MOLIT), daily streamflow from the watershed was measured in real-time.The National Institute of Environmental Research (NIER) monitored the water quality, including TDS, DO, BOD, NH 4 -N, NO 3 -N, NO 2 -N, and COD at the watershed outlet (Osan stream 1, 1.2, 1.3, 2, 2.1 and 3 stations) on a monthly basis.

Figure 1 .
Figure 1.Location of the study area.

Figure 1 .
Figure 1.Location of the study area.

Figure 3 .
Figure 3. Flowchart of the optimization model.

Figure 3 .
Figure 3. Flowchart of the optimization model.

Figure 4 .
Figure 4. Water quality simulation 2016 results for the Osan stream drain.

Figure 4 .
Figure 4. Water quality simulation 2016 results for the Osan stream drain.Water quality properties (TSS, TN, BOD and DO concentrations) and flow rates measured on 25 January 2016 and 24 July 2016 are used to create models for the summer and winter seasons.

Figure 4 .
Figure 4. Water quality simulation 2016 results for the Osan stream drain.

Figure 5 .
Figure 5. Error bars and frequency distribution from uncertainty analysis for BOD, DO, TN and TSS for summer and winter at Osan stream 2-1.

Figure 5 .
Figure 5. Error bars and frequency distribution from uncertainty analysis for BOD, DO, TN and TSS for summer and winter at Osan stream 2-1.

Figure 6 .
Figure 6.Water demand, availability and reuse in the Osan stream drain.

Figure 6 .
Figure 6.Water demand, availability and reuse in the Osan stream drain.

Figure 7 .
Figure 7. Objective function vs iteration for Jaya Algorithm.

Table 1 .
Water quality average at the monitoring points of Osan stream drain from January 2015 to July 2016.

Table 2 .
Mean water quality for summer and winter 2016 at proposed reused stations.

Table 3 .
System parameters in the Osan Stream.

Table 4 .
Monte Carlo Simulation summary for BOD, DO, TSS, and TN at proposed reuse station.

Table 4 .
Monte Carlo Simulation summary for BOD, DO, TSS, and TN at proposed reuse station.

Table 5 .
Jaya Algorithm results at proposed reuse stations during summer and winter 2016.