Economics of alternative water resources with an emphasis on aquifer storage and recovery

This paper demonstrates a comprehensive methodology for assessing the comparison of unit water production cost (UWPC) between alternative water resources including desalination, freshwater reservoirs, single-purpose dams, underground dams and two indirect water in take technologies – riverbank filtration and aquifer storage and recovery (ASR). This study considers the Monte Carlo simulation as the only viable solution to tackle this critical question, which can be used to evaluate the economics of diverse water supply schemes incorporating those alternatives and prepare long-term water supply planning. Built upon actual and conceptual cost data for each alternative, total project cost and operation and management cost estimation models for each alternative were developed and used for generating mean UWPC information using the Monte Carlo simulation approach. The mean UWPC differences between alternative water supply schemes were found to be statistically significant and the simulation results revealed that ASR is the lowest-cost option to provide drinkable water for both cases when a conventional water treatment plant (WTP) and advanced WTP were used as a connected post-treatment process. This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/). doi: 10.2166/ws.2017.137 om https://iwaponline.com/ws/article-pdf/18/2/612/206678/ws018020612.pdf 2019 Jae-ho Choi (corresponding author) Young-Gyoo Shim DongA University, Busan, South Korea E-mail: jaehochoi@dau.ac.kr Miroslaw Skibniewski Department of Civil and Environmental Engineering, University of Maryland at College Park, College Park, USA


INTRODUCTION
Alternative water resources are usually perceived as new water supply methods, particularly in islands and mountainous areas, where the direct in take of water through conventional methods such as dams, surface water, and direct use of groundwater is difficult to adopt (Park ).
Despite no clear definition of alternative water resources, desalination (DES), underground dams (UGDs), and riverbank filtration (RBF) are classified as alternative water resource projects because they can secure water sources, increase the cost effectiveness of producing water, and meet other site-specific purposes (K-water a).Another alternative method that is now gaining greater attention as a new source of water and under investigation in several countries including Korea, China, Qatar, and a few more Middle Eastern countries is aquifer storage and recovery (ASR).The escalating costs and environmental challenges of conventional water supplies that use mostly surface water have encouraged water professionals to explore ASR.
ASR was defined by Pyne () as 'the storage of water in a suitable aquifer through a well during times when water is available, and recovery of the water from the same well during times when it is needed.'The vadose zone plays a filtering, absorbing, and biologically mediated reacting role in removing or reducing the concentrations of suspended solids, biodegradable organic matter, nutrients, metals, and pathogenic microorganisms.The beneficial uses of ASR for coping with various water-related issues, such as securing drinking water resources, preparing water resources for disaster, preventing jurisdictional water disputes, preventing seawater intrusion in coastal regions, and so forth, are gradually being recognized throughout the world.
More than 134 ASR well-fields (over 544 ASR wells) are operational in 22 states in the USA, and other countries, including England, Australia, The Netherlands, Namibia, South Africa, United Arab Emirates, India, Canada, and Israel, have either experimental or operational well-fields (Bloetscher et al. ; Pyne ).
The principal reason for the increased interest in ASR technology from many countries in the last decade can be attributed to its cost effectiveness.In almost all situations where the need exists for alternative or supplemental water supplies, ASR wells can provide this water at less than half the capital cost of alternative water sources.The operating costs are not well known but are believed to be similar to, and perhaps slightly higher than, typical production well-field operating costs (Pyne ).However, Maliva () questioned the lesser degree of its adoption in the market and attributed the reason to a lack of sound economic cases for decision makers such as utility managers, water management agency officials, and political leaders.
In 2013, a 5-year-long national research project to develop technology packages for the practical use of largescale ASR was funded by Korea's Ministry of Land, Infrastructure, and Transport (Choi et al. ).One of the most sensitive and critical aspects of this project, particularly its market adoption, is its cost.This study serves to provide the economics of a drinkable water supply scheme using ASR by comparing its unit water production cost (UWPC) to those of other schemes using different water resources (hereinafter 'alternatives') such as RBF, freshwater reservoir (FWR), single-purpose dam (SPD), UGD, and DES.FWR and SPD are mainly used for drinkable water production purposes in Korea.The UWPCs of a water supply scheme using an alternative water resource project reflect all the capital and operation and management (O&M) costs from constructing and operating alternatives, in take facility, water transmission pipeline, and post-water-treatment utility.Therefore, building a thorough understanding of the water supply project schemes, water treatment processes, construction methods, and work scopes of an alternative water resource project is desirable to form the same baseline for the UWPC determination.
Previous studies (K-water ; Carollo Engineers ; CDM ) that investigated the economics of alternative water supply systems have either used cost opinion information or focused on one specific alternative of interest, which was compared with a few other conventional water supply methods.Collectively, most of the related studies used cost opinions with a wide range of accuracy and focused on deterministic UWPC calculation, while no studies have examined the uncertainties and wide ranges of variables that affect the UWPC.Given the inherent uncertainty of input variables and estimates of UWPC, the only way to compare the UWPCs between alternatives is through the means of cost uncertainty analyses that are based on Monte Carlo simulation including the Latin hypercube sampling method (hereinafter 'MCS').In this regard, the major contribution of this study is a direct comparison of UWPCs between water supply schemes with different alternatives based on actual as-designed and as-built project cost data by using statistical simulation.

METHODS
The objective can be achieved by a four-step procedure based on an extensive compilation of total project costs (TPCs) and O&M costs in relation to the various water resources and water treatment projects in the context of Korea as follows.
Step (1): identify alternative water supply schemes consisting of alternative water resource projects such as DES, SPD, FWR, UGD, RBF, and ASR, and postwater-treatment utilities such as conventional water treatment plant (WTP) and advanced WTP (hereinafter 'AWTP'), and collect their TPC and O&M costs based on real project documents and conceptual feasibility reports.
Step (2): develop TPC and O&M estimation models for each alternative by using regression analysis with a single independent variable (water development capacity, m 3 /d).
TPC estimation and O&M estimation models for conventional and advanced WTPs were also developed.Step (3): set up variables that influence the UWPC estimation of different water supply schemes using each alternative as source water by using MCS, including the interest rate, service life, pipeline extension, pavement material cost, and post-treatment utility types.Step (4): compare UWPCs among these alternative water supply schemes by using different variable conditions, followed by a series of statistical methods to validate the mean difference between any pair of the schemes and the overall assessment of the comparison results.

Database compilation for UWPC estimation
All the alternatives produce source water that is transported to a standard WTP or an AWTP except DES, which mostly uses a simplified post-treatment process.The authors compiled all the data required for determining UWPC for each alternative water supply scheme by using real as-design and as-constructed project data and references (K-water , b; Lee ; Moon et al. ).Since ASR has not yet been implemented in Korea, cost information for a total of 11 ASR well-fields and four ASR O&M cost data reported by CDM () were examined and applied for this study.This study only used four cases in which the aquifer was recharged with either the treated surface water or partly treated surface water, and the recovered water from the aquifer was used for public water systems.
Figure 1 shows detailed information on the number of projects (4,86,17,3,18,19 in the order of items on the X-axis) and the relative distributions of the capacity ranges and TPCs.The numbers of alternatives were not evenly distributed, and the capacity ranges and TPCs varied between alternatives.The largest proportion belonged to DES, followed by UGD, SPD, FWR, ASR, and RBF.The medians of DES and FWR in terms of the capacity range were much lower, whereas the medians of all the alternatives in terms of the TPC were not much different, except for SPD.
These results imply that DES and FWR were planned and operated on relatively smaller scales, while SPDs demanded higher project costs.The base time-point for comparing the UWPCs between various water supply schemes using the selected alternatives was set to August 2006, which was

TPC estimation and O&M cost estimation models
TPC estimation models for each alternative were developed based on the compiled dataset from the previous section, as shown in Table 1.Each model produced a TPC for August 2006 and could be used only within the capacity range in the last column, from which the models were developed.
The goal of this study is to determine the unit cost of producing drinkable water that is treated by either conventional WTP or AWTP that sources water from an alternative, so the TPC models of the WTP and the AWTP are also required for the UWPC calculation.Cost data for a total of 25 rapid sand filter WTPs were compiled to develop a reliable TPC model (Equation (1)) by using regression analysis.Other types of WTPs that use slow sand filters and membranes were not considered because their numbers were very limited; in particular, they have been used on much smaller scales below 15,000 m 3 /d development rates.In addition, most FWR, UGD, and SPD projects in our data were highly associated with rapid sand filter WTPs.
Total annual O&M cost estimation models for each alternative were based on O&M cost data from real projects or O&M cost determination criteria (K-water ; KDI ), as shown in Table 1.Annual ASR O&M costs were based on two real project data that were provided by CDM (), which explains why R 2 equals 1.The annual O&M costs of each alternative water supply scheme using FWR, UGD, and SPD are assumed to be a combination of the O&M cost of each alternative resource infrastructure and the O&M cost of the following conventional WTP.The annual total O&M cost, for example, of a WTP that used FWR as its primary source of water was the combination of y FWP and y WTP .This study also assumed that the annual O&M costs of UGD and SPD are 0.5% of the construction cost according to the guideline published by the Ministry of Construction and Transportation (MOCT ), while FWRs were set to 0.3% because this alternative has fewer ancillary structures than UGD and SPD.The annual O&M cost estimation model (i.e.y WTP ) of a conventional WTP that sources water from alternatives was developed based on 26 real project O&M cost data by using regression analysis (Equation ( 2)).

Factors that affect the UWPC modelling
Once the TPC and annual O&M cost estimation models of each alternative and the TPC and O&M cost estimation models of WTPs were prepared, we could determine reasonable and commonly accepted values of factors that affect the UWPCs of various water supply schemes that incorporate each alternative.We subdivided the input factors into four groupseconomic, construction, operation and utility system parameter groupsas shown in Table 2. MCS was utilized to accommodate the probabilistic distributions of input factors and different water supply schemes to mimic the real data conditions.Economic factors consider the time value when determining the equivalent uniform annualized cost (EUAC), which is defined by the cost per year of owning and operating an asset over its entire lifespan.In particular, the EUAC is more effective to decision makers because organizations typically report their annual activities and need to compare alternatives with unequal lifespans with greater computational efficiency (Jones & Smith ). where x unit: 1,000 m 3 /day (applicable development range: 1,000-1,000,000 m 3 /d), y unit: 1,000 USD (August 2006). where x unit: 1,000 m 3 /day (applicable development range: 1,000-1,000,000 m 3 /d), This study applied the EUAC method to compare the cost effectiveness between different water supply schemes utilizing different water resource alternatives.The total EUAC of a water supply scheme is the sum of the EUAC of all the infrastructures comprising the water supply scheme including an alternative, transmission pipeline,  Finally, the utility system parameter is related to the quality of the source water and integration of alternatives into the following water utility.

Post-water-treatment utility selection
The With regard to the O&M cost of AWTP, the processes that are added onto conventional WTP cause a 5.8-6.5% increase in O&M costs, which are modelled as a Gaussian distribution in a MCS, as shown in where x unit: 1,000 m 3 /day (applicable development range: 28,000-85,000 m 3 /d), y unit: 1,000 USD (August 2006).

RESULTS AND DISCUSSION
MCS was conducted to simulate real data conditions and determine the means of UWPCs for each alternative based on probabilistic and deterministic input variables.MCS first needs to set up variables that influence the UWPC estimation of different water supply schemes using each alternative as source water, including the interest rate, service life, pipeline extension, pavement material cost, and post-treatment utility types.Each variable is assumed to follow a Gaussian distribution, and every random value is generated through the Latin hypercube sampling.Figure 2 is the plot of the means of 10,000 simulations for different water supply schemes using each alternative.This study also performed two simulations by setting the post-watertreatment utility of RBF and ASR to option 1 (conventional WTP) and option 2 (AWTP), respectively.The other alternatives are linked only to conventional WTPs, whereas DES was connected to the simplified post-treatment option.
Figure 2 shows the simulated UWPC means of each alternative from the first MCS with 10,000 iterations over the water development capacity.Conventional water treatment processes were used for RBF and ASR in this simulation, and the UWPCs for each alternative were valid and drawn only within the capacity range, from which the O&M cost estimation models were developed.
All the UWPCs and the differences in the means between alternatives became smaller as the water development capacity increased as seen in Figure 2. The overall trends of the UWPCs are divided into higher UWPC alternatives and lower UWPC alternatives with a dark dotted line.
Some points exist where the UWPCs of two different alternatives overlap (i.e., ⓐ-ⓒ), and the largest numbers of these alternatives including ASR can be compared (i.e., ⓓ).UWPC ranking in descending order from 40,000 to 45,000 m 3 /d was DES, SPD, RBF, UGD, and ASR although the mean differences between alternatives could not be statistically proven.At point ⓒ, similar to points ⓐ and ⓑ, the mean difference between two alternatives (i.e., DES and SPD) was also very close and thus needed to be validated with statistical follow-up procedures.
The recommended approach to find statistical mean differences between more than two groups (i.e., alternatives) is to use analysis of variance, or ANOVA.This approach tests all the means simultaneously while considering within-group variability to determine whether at least one reliable difference exists between any of the group means or a linear combination of group means in the set of means.However, this method requires follow-up procedures, which are called multiple comparisons, to precisely pinpoint sources of significant differences

CONCLUSIONS
The MCS results practically demonstrate that DES, SPD, and
and WTP, and the sum of the annual O&M cost of the infrastructures.The construction parameter is related to the material type and extension of water transmission pipelines and the pavement type.This study used the cost information from the research report (MOE ) in relation to materials and construction of the pipeline, which are classified by pavement type(i.e., soil, asphalt, concrete).The operation parameter determines the O&M cost of each alternative.
selection of a post-water-treatment utility for the two representative indirect alternatives, RBF and ASR, are subject to changes according to the amount, the location, and particularly the quality of the source in take water.Determining the UWPCs of water supply schemes that use these alternatives therefore must reflect added project cost components and O&M cost variations depending on the type of connected post-water-treatment utility.The post-water-treatment utility of the three RBF projects that were used to build our estimation models used advanced water treatment processes.Most of the advanced water treatment processes are added to conventional WTPs, thus increasing capital and O&M costs.This study analysed a total of 18 AWTP expansion projects, in which four cases added the powdered activated carbon process, eight cases added the post-ozone and biologically activated carbon (BAC) processes, and the rest added pre-and post-ozone and BAC processes to the conventional water treatment processes.This study selected the project expansion cost models that were built based on the first four cases by using regression analysis because the cost model was only applicable to the capacity range of interest in this study, as shown in Equation (3).

The
first section, which has a development capacity of 1,000 m 3 /d-5,000 m 3 /d, illustrates that the lowest-cost alternative (i.e., DES) at 1,000 m 3 /d becomes the highestcost option at 5,000 m 3 /d.The relative UWPC ranking orders from the 5,000 m 3 /d and 15,000 m 3 /d points until where the UWPCs of FWRs are available are DES, SPD, FWR, and UGD in descending order, although the UWPCs between SPD and FWR at point ⓐ are difficult to discern.Before and after point ⓑ, where RBF and UGD have similar UWPCs, their relative orders are reversed and the UWPC of RBF becomes higher than that of UGD because the RBF cost estimation model was developed based on the low-cost vertical well option, which is usually used with lower development capacities (e.g., 20,000 m 3 /d), and the high-cost collection well option, which is usually required for higher development capacities (e.g., 60,000 m 3 /d).The relative

Figure 2 |
Figure 2 | MCS for estimating UWPCs for each alternative (at 10,000 iterations for RBFs and ASRs with a conventional WTP connected).

(
Carlberg ).Point ⓓ was tested to determine whether any significant difference existed between the mean values for both cases: CASE 1 had a difference at 40,000 m 3 /d with a conventional water treatment utility, and CASE 2 had a difference at 45,000 m 3 /d with an advanced water treatment utility.Figure3shows the probability density results from the MCS for both cases.DES shows a narrower density distribution because no water transmission pipelines and connected water treatment utilities were present.One well-regarded multiple comparison procedure for controlling Type I error, or the incorrect rejection of a true null hypothesis, is Bonferroni's adjustment.The Bonferroni adjustment is often considered to be overly conservative because this method simply divides the Type I error rate (0.05) by the number of tests (15 for 40,000 m 3 /d, 10 for 45,000 m 3 /d).
utilities that were connected to RBF, the p-values at points ⓐ, ⓑ, and ⓒ were 0, 0.00046 and 0, respectively, which all indicate a significant difference in the mean values for those pairs of alternatives at a significance level of 0.01.

Figure 3 |
Figure 3 | Probability density for alternative water supply schemes: (a) water development capacity at 40,000 m 3 /d for a conventional water treatment utility with RBF and ASR (FWRs are assumed to be available); (b) water development capacity at 45,000 m 3 /d for an advanced water treatment utility with RBF and ASR (no FWR information available).

FWR
are considered higher UWPC alternatives, whereas UGD, RBF, and ASR are considered lower UWPC alternatives.DES costs less than FWR and UGD at a development capacity of 1,000 m 3 /d but is the most expensive alternative over all the development capacities mainly because of higher O&M costs.FWR is the best option for lower development capacities from 5,000 to 15,000 m 3 /d among the higher UWPC alternatives.RBF is a highly competitive option at relatively small water development capacities; for instance, RBF is the cheapest option at 20,000 m 3 /d, but its position is overtaken by UGD over 25,000 m 3 /d.However, these findings must be carefully adjusted and used in reality by considering the availability, quality, and recovery capacity of injected source water into an aquifer and by projecting the outer environment, such as the availability of a nearby WTP and transmission pipeline extension.As was indicated by two articles written by Bloetscher et al. () and Pyne (), ASR projects have a higher operational risk mainly due to clogging issues resulting in an earlier closure of operation wells than planned.If the risk of failure in satisfying the life expectancy of ASR is reflected in the UWPC simulation, the statistical result could be changed.It also needs to address the limitation of applying the findings produced in the previous section because the results are based on a limited number of projects, particularly related to ASR and RBF.For a complete uncertainty analysis of UWPC estimation for the alternatives, not only the variabilities of the input factors summarized in Table2, but also the types and confidence levels of the developed regression models need to be considered in a hybrid approach combining both sensitivity analysis and Monte Carlo analysis for a future research direction.Considering the level of water availability particularly for RBF in the simulation process may lead to a different result since RBF is highly dependent on river flow rates, which are susceptible to future climate change.UWPC information that considers various input variables and water supply schemes could be used for policy or decision makers to swiftly measure the viability of selected water resource options in a specific project environment.Future studies should include incorporating the operational risk and water security issues into the UWPC simulation of alternative water supply systems and monetizing the benefits of ASR to conduct benefit-cost analyses.

Table 2 |
Parameter assumptions for UWPC MCS a Average and standard deviation based on three pavement types: asphalt, concrete, soil.b Steel pipeline (900 mm), topsoil depth: 1.5 m. c Rapid sand filter WTP.

Table 2
. In a similar way, recovered water from ASR sites may require further treatment to meet drinking water standards, although the quality is improved through the removal of pathogens and nutrients(Steinel ).ASR projects that are designed to produce public water are generally known to use a simpli-

Table 3
RBF, UGD, and ASR.For conventional water treatment

Table 3 |
Pairwise T-test comparisons with the Bonferroni adjustment at 40,000 m 3 /d and at 45,000 m 3 /d