Design Aspects, Energy Consumption Evaluation, and O ﬀ set for Drinking Water Treatment Operation

: Drinking water treatment, wastewater treatment, and water distribution are energy-intensive processes. The goal of this study was to design the unit processes of an existing drinking water treatment plant (DWTP), evaluate the associated energy consumption, and then o ﬀ set it using solar photovoltaics (PVs) to reduce carbon emissions. The selected DWTP, situated in the southwestern United States, utilizes coagulation, ﬂocculation, sedimentation, ﬁltration, and chlorination to treat 3.94 m 3 of local river water per second. Based on the energy consumption determined for each unit process (validated using the plant’s data) and the plant’s available landholding, the DWTP was sized for solar PV (as a modeling study) using the system advisor model. Total operational energy consumption was estimated to be 56.3 MWh day − 1 for the DWTP including water distribution pumps, whereas energy consumption for the DWTP excluding water distribution pumps was 2661 kWh day − 1 . The results showed that the largest consumers of energy—after the water distribution pumps (158.1 Wh m − 3 )—were the processes of coagulation (1.95 Wh m − 3 ) and ﬂocculation (1.93 Wh m − 3 ). A 500 kW PV system was found to be su ﬃ cient to o ﬀ set the energy consumption of the water treatment only operations, for a net present value of $0.24 million. The net reduction in carbon emissions due to the PV-based design was found to be 450 and 240 metric tons CO 2 -eq year − 1 with and without battery storage, respectively. This methodology can be applied to other existing DWTPs for design and assessment of energy consumption and use of renewables.


Introduction
Energy conservation and sustainability is one of the pressing issues being faced by the water industry today as drinking water treatment and wastewater treatment are highly energy-intensive processes [1]. These processes account for about 2% of the United States' total energy use, resulting in the emissions of about 45 million tons of greenhouse gases (GHG) [2]. One of the ways to reduce carbon emissions is to incorporate renewables into existing water infrastructure. The objective of this study was evaluation of the energy consumption of a large-scale drinking water treatment plant (DWTP), associated carbon emissions, and then offsetting that energy consumption by utilizing photovoltaics (PVs).
Energy use related to drinking water treatment may include energy for water conveyance, energy for unit operations of the DWTP, for facility related equipment (lighting, heating and ventilation), and energy for water distribution [3,4]. Energy use for water treatment depends on various factors such as raw water quality, water source, age of water delivery system, conveyance distance, water storage of solar energy [53,54]. If the landholdings of the treatment plant are not sufficient, then additional land area would need to be acquired.
Design and performance analysis of solar PV can be achieved using available modeling tools. The system advisor model (SAM) developed by the National Renewable Energy Laboratory in 2005 [55] is a performance and cost-modeling tool and facilitates in decision making for renewables of solar, wind, and geothermal and biomass combustion. Various studies have used SAM for analyzing solar technologies, particularly photovoltaics [56][57][58][59]. Ref. [59] utilized SAM to generate hourly performance simulations for the validation of their solar PV design. Solar PV was used to generate power for the aerators, utilized for the treatment of acid mine drainage at Hwangji Youchang facility, located in Korea. The 30.1 kW PV system, using a factor of safety of 5, resulted in electricity production of 3016 kWh month −1 , to meet the 342.39 kWh month −1 electricity demand [59]. Ref. [56] utilized SAM to generate life cycle costs per kWh for utility-scale development of solar PV as well as for using solar PV to meet the electricity demands for residential and commercial sectors in the United States. Ref. [58] employed SAM for analyzing life cycle energy production and associated costs of a solar PV system and solar water heater for a residential dwelling unit in Houston, Texas.
Solar energy has great potential as a clean source of energy for water treatment processes [7,[60][61][62], wastewater treatment processes [63][64][65][66], and desalination [67][68][69][70][71]. Studies evaluating the potential of existing DWTPs to incorporate solar energy technologies are limited. Ref. [61] made performance comparisons between a conventional DWTP and a wind-solar-powered nano-filtration pilot plant located in Spain. Ref. [62] explored the potential of using wind and solar energy for a large-scale water treatment plant in the Netherlands. The goal of this study was to: (a) Design and determine the energy consumption and energy intensity of each unit operation (validated using plant's data) of an existing DWTP, (b) Conduct a modeling study and size the DWTP for solar PV, to offset the energy consumption of the plant, based on available landholdings and economic assessment, (c) Determine the net reduction in carbon emissions due to solar PV installation compared to non-PV based design.
The approach used in this study can be utilized as a guide to design DWTPs and analyze their energy consumption. This is useful for engineers and researchers as generally consulting firms do not publicly disclose such information. For a meaningful design, data was collected from an existing DWTP. Further, sustainability is incorporated into the design of the plant by using solar PV, with the long-term goal of energy independence and reduced carbon emissions for an existing DWTP.

Study Area
For this study, the selected treatment plant is located in southwestern United States, serving a city with a population of 0.24 million. The primary water supply for the city is a river originating from a lake, enduring a 161 km course through different canyons and ultimately having an outfall in another lake. Only 6% of the water flowing through the river is used by the community, and in addition, almost half of that 6% is returned to the river after the generated wastewater is treated at the wastewater treatment plants. Hence, only 3% of the water flowing through the river is used in a non-drought year, while 8% is used in a drought year. The DWTP has the capacity to treat 3.94 m 3 of source water per second. For confidentiality reasons, the location of the plant cannot be disclosed, as required by the DWTP's management, since data was collected from the DWTP's managers for validation of the results generated in this study. However, for the purpose of analysis of the PV system, it was assumed that the plant was located in Nevada (direct beam insolation level 6.3 kWh m −2 day −1 ), having insolation levels favorable for solar deployment.

Data Sources
Data sources for this study are as follows: Water 2020, 12, 1772 4 of 27

Process Flow Diagram
A process flow diagram was acquired from the treatment plant managers. The process flow diagram discussed in this section for the DWTP is based on the information provided. Figure 1 shows the process flow diagram for the plant, which utilizes a conventional treatment train to treat 3.94 m 3 of river water per second. Raw water from the river at first undergoes pre-treatment and passes through bar screens, raw water basins, and then automatic screens for removal of larger debris, settlement of heavier grit, and removal of smaller debris, respectively. Then, the water is coagulated and flocculated. Next, the water is filtered through dual-media filters. Lastly, the water is disinfected by application of sodium hypochlorite and then stored for distribution throughout the city.
Water 2020, 12, x FOR PEER REVIEW 4 of 30

Data Sources
Data sources for this study are as follows:

Process Flow Diagram
A process flow diagram was acquired from the treatment plant managers. The process flow diagram discussed in this section for the DWTP is based on the information provided. Figure 1 shows the process flow diagram for the plant, which utilizes a conventional treatment train to treat 3.94 m 3 of river water per second. Raw water from the river at first undergoes pre-treatment and passes through bar screens, raw water basins, and then automatic screens for removal of larger debris, settlement of heavier grit, and removal of smaller debris, respectively. Then, the water is coagulated and flocculated. Next, the water is filtered through dual-media filters. Lastly, the water is disinfected by application of sodium hypochlorite and then stored for distribution throughout the city.

Water Quality Report
A raw water quality report was acquired from the DWTP's managers, and any missing data was reproduced from [43]. The design for each unit process is dependent upon the raw water characteristics as well as the corresponding maximum contaminant levels (MCL) and secondary maximum contaminant levels (SMCL) determined by the U.S. Environmental Protection Agency (USEPA), as shown in Table 1 [72].

Treatment Plant Design Criteria
Sizing was achieved by using industry accepted design criteria [43,[73][74][75] and discussed in detail in Section 4.0.

Water Quality Report
A raw water quality report was acquired from the DWTP's managers, and any missing data was reproduced from [43]. The design for each unit process is dependent upon the raw water characteristics as well as the corresponding maximum contaminant levels (MCL) and secondary maximum contaminant levels (SMCL) determined by the U.S. Environmental Protection Agency (USEPA), as shown in Table 1 [72].

Treatment Plant Design Criteria
Sizing was achieved by using industry accepted design criteria [43,[73][74][75] and discussed in detail in Section 4.

Carbon Emissions
Data sources and parameters for carbon dioxide emissions and southwestern state's energy-source mix for electricity generation is shown in Table 2. Carbon emissions data in units of gCO 2 -eq kWh −1 , for different energy sources for electricity generation, was obtained from [76]. Nevada's electricity source mix for various energy sources was obtained from U.S. Energy Information Administration (USEIA) [77].  [76] and Nevada's electricity source mix [77].

PV System
Parameters used for sizing solar PV, related to the PV module, inverter, and battery storage are shown in Table 3. The PV system was sized as 2-axis tracking and consists of a mono-crystalline silicone module assumed to be ground mounted. A 2-axis tracking system allows higher energy output compared to single-axis trackers or fixed systems. The module type selected was Helio USA 7T2 305 [78]. Inverters are used to convert the direct current output of the solar modules into alternating current, so that the generated electricity can be utilized by the DWTP or the excess generation can be utilized by the electricity grid. The inverter type selected was Fronius Symo 10.0 240V. The battery type chosen was lead acid flooded [79,80] (Table 3).  Various financial parameters incorporated within the economic model of SAM, are listed in Table 4. These parameters were based upon the review of published literature in the years 2016 and 2017 [81][82][83][84][85][86][87][88].  [49]. These incentives were applied in this study for economic assessment. Federal investment tax credit (ITC) represents a percentage of investments that the owners are allowed to deduct from their taxes, dollar-to-dollar. Solar projects of a commercial, industrial, or agricultural nature, implemented before or during the year 2019, can claim 30% ITC, while 26%, 22%, and 10% ITC can be claimed by projects implemented during 2020, 2021, and 2022 and onwards, respectively. If federal taxes owed are less than 30% of installed costs, solar PV owners owe no taxes for that year, and any unused tax credit is carried over to the next year. A state property tax exemption incentive, applicable for the design life of solar PV, was also incorporated into the economic analysis. This incentive cannot be claimed if another state tax abatement or exemption is claimed by the same building. Therefore, only these two incentives were incorporated in the analysis.

Methodology
The first step of the current research focused on sizing the DWTP by designing the various unit processes involved. Sizing was achieved by using industry accepted design criteria [43,73,74]. Design of the plant and estimation of energy consumption were determined for the maximum day demand flow anticipated for the design life of the plant (3.94 m 3 s −1 ). Thus, the analysis is representative of the worst-case scenario for the DWTP. Other parameters chosen were also reflective of the extreme conditions.
Each unit process was designed and its energy driving units were identified; then, energy consumption (kWh day −1 and kWh m −3 ) for each unit process (coagulation, flocculation, sedimentation, filtration, chlorination, residual management) was individually determined and validated using the plant's data (Section 4.1). Computations were made using Microsoft Excel Spreadsheet. Next, the system advisor model (SAM) was used to conduct a modeling study to determine the techno-economic performance of using solar PV (by sizing solar PV and battery storage based on economic metrics) to offset the energy consumption of the DWTP (Section 4.2). Land requirements of the PV deployment were compared against available landholdings of the plant to check for the sufficiency of the available area for PV deployment (Section 4.3). Net reduction in carbon emissions was determined by comparing a PV-based design (modeling study) to a non-PV based design (actual plant operation) (Section 4.4).

Pre-Sedimentation
Pre-sedimentation initiates suspended solids removal and equalizing of variable turbid loadings entering a DWTP. Typically, plants treating river water are equipped with raw water basins, as river water is high in turbidity [92]. Two tanks should be provided at a minimum to incorporate redundancy in the system. Surface loading rates vary between 200 and 400 m 3 m −2 day −1 , while detention time varies between 6 and 15 min [43].

Coagulation
Application of a coagulant assists in destabilizing a colloidal suspension, causing collides to agglomerate and form flocs. A coagulant's dose and type is dependent upon the results of a jar test. The selected DWTP utilized ferric sulphate Fe 2 (SO 4 ) 3 as a coagulant. After the coagulant is added, rapid or flash mixing ensures uniform dispersion of the coagulant throughout the basin. In this case, a static mixer was used because mixing was achieved in the conveyance pipe between the coagulant intake pipe and flocculation basin. At least two mixers should be provided. A Reynolds number >2000 implies turbulence and good coagulation performance [75]. Velocity gradient (G) remains a significant parameter for the design of mixers. The equation to compute the G value of a static mixer [75] is as follows: where G = root-mean-square velocity gradient (s −1 ), P = power imparted to water (N-m s −1 ), V = mixer volume = product of pipe area and length of mixer (m 3 ), and µ = dynamic viscosity (N-s m −2 ). Energy consumption (E) in the units of kWh day −1 of the mixers was found using the following equation [75].
where t m = motor run time (h day −1 ), and f = motor efficiency factor. Typically, motor efficiency is in the range 0.7-0.9 [93]. For this study, a value of 0.8 was used. Since the plant operates 24 h per day, a motor run time of 24 h was used. Energy consumption of the metering pump was determined using the brake horsepower (BHP) equation [94], which is as follows: where γ = specific weight (kN m −3 ), Q = flow rate (m 3 s −1 ), H = total dynamic head (m), and e = wire-to-water efficiency. Wire-to-water efficiency can be estimated as the product of pump, motor, and drive efficiency. A jet diffuser pump was used to flash-mix coagulant aid polymer. Power (P) estimations were made using following equation [75].
where v jet is the jet velocity emerging from the orifice (m s −1 ), and g = acceleration due to gravity (9.81 m s −2 ). As in Equation (3), Q = flow rate, and γ = specific weight of water.

Flocculation
After coagulation, the flocculation process assists in the formation of flocs, which are later separated through the processes of sedimentation and filtration. In contrast to the coagulation process, which involves flash mixing to disperse the coagulant within the feed water, the flocculation process involves slow mixing to enhance the contact between the coagulant and the feed water and for the subsequent formation of flocs. Paddle wheels are used for the flocculation process (rotational speed 0.3-3 rev min −1 ). The coagulated water flows into the flocculation basins, where flocculation is achieved in stages with varying velocity gradient (5-40 s −1 ) [43]. The following equation was used to estimate the energy consumption of paddle wheel flocculators.
where n f = number of flocculation basins, Q f = flow rate in each flocculation basin (m 3 min −1 ), and t s = detention time per stage (minute). The variables G, µ, f, and t m are the same in as in Equation (2).

Sedimentation
After the water is coagulated and flocculated, it moves into the sedimentation basins so that heavier flocs settle out of suspension. Compared to conventional gravity sedimentation, parallel plate settlers provide enhanced solids removal capabilities in a smaller area, even during periods of extremely high solid loadings to the plant. Parallel plate sellers were designed based on the criteria provided by [43,75,95]. The energy consumer for the parallel plate settlers were sludge transfer pumps, and their brake horsepower can be estimated using Equation (3).

Filtration
A dual-media (anthracite and sand) filter is employed to trap the flocculated water, using a rapid filtration process, to provide a finished water quality of ≤0.3 NTU and 4 log removal of cryptosporidium. Parameters used for the filter design [74] include filtration rate (14.5 m h −1 ), effective size of filter media (anthracite 1 mm, sand 0.5 mm), filter media depth (anthracite 1.4 m, sand 0.25 m), and net available head (1.8-3.0 m).
Typically, to arrive at the plant the water has to be pumped from the water source to the plant. For the selected DWTP, raw water flows by gravity from the river, passes through bar screens and flows by gravity through the rest of the plant. The raw water usually flows by gravity between one unit to the next, in a typical water treatment facility. In the selected plant, after undergoing the processes of coagulation, flocculation, and sedimentation, raw water flows by gravity into the filtration unit. The filtered water is chlorinated and then moved into the finished water reservoir. The energy-consuming part of the filtration process is the backwash system (backwash pumps and air scour). Air blowers operate by developing a pressure differential to move air between the entrance and exit points [75]. Air scour of the filter media was achieved for 4 min before the application of water backwash. The power estimation for the air blower was made using the following equation [96].
where w = air mass flow (kg s −1 ), R = universal gas constant (kJ k −1 mole • K), T = air temperature at inlet ( • K), P 1 = absolute pressure at entrance (Pascal), P 2 = absolute pressure at exit (Pascal), e = efficiency = 0.8, n = (k−1) (k −1 )= 0.283 for air, and 8.41 is constant for air (kg k −1 mole). The filter is backwashed when the head loss exceeds the available head or increased turbidity is displayed in the filter effluent. The backwash system was designed [43] using the data obtained from the DWTP (backwash duration 9 min, backwash frequency 24 h). The energy consumption of the backwash pumps can be determined using Equation (3).

Chlorination
Disinfection can be accomplished through chlorination, ultra-violet radiation treatment, or ozonation. For the DWTP, filtered water was disinfected through chlorination using sodium hypochloride (NaOCl) applied in liquid form. The chlorination design was achieved using the design criteria provided by [93,[97][98][99]. The surface water treatment rule was applied for 4 log inactivation of viruses and 3 log inactivation of giardia. A minimum one-hour detention at average design flow or 30 min detention at peak hourly flow was used, whichever was greater. The energy consumers for the chlorination process were the metering pumps; the energy consumption can be determined using Equation (3). Finished water was distributed throughout the city using highly energy-intensive zone-1, zone-2, and zone-3 pumps.

Residual Management
Filtration backwash wastewater and sedimentation sludge generated during filtration and sedimentation is treated using decant basin and water recovery basins, designed using the same criteria as pre-sedimentation basins. At these facilities, plant waste solids are concentrated and dewatered to recover water and to facilitate ultimate solids disposal. Decanted water is recycled through the plant, while the decant basin settled solids are pumped by solids transfer pumps to water recovery basins for secondary solids settling. Dewatering of drinking water sludge will vary depending upon coagulant type used, solids loading, and pH [100].
Energy consumption was evaluated for two main scenarios for the selected treatment plant. Scenario 1 (S-1) represents analysis for the DWTP including water distribution pumps, whereas Scenario 2 (S-2) represents the analysis for the DWTP without including the water distribution pumps.

System Advisor Model
Energy consumption for each unit process of the DWTP was determined (Section 4.1) and used as an input for SAM, a modeling tool utilized to evaluate the techno-economic performance of using solar PV to offset the energy consumption of the DWTP.

Technical Analysis
SAM evaluates energy performance and cost for grid-connected systems based on location and system design parameters. Electric load of the DWTP, electric rate, typical meteorological year (TMY3) data of the location, as well as the parameters related to PV panels were the primary inputs for SAM. An analysis period of 25 years was used.
The energy consumption estimate determined using the steps outlined in Section 4.1 was used as an electric load input for sizing solar PV. Electricity generation by the PV system meets the facility's electrical load that would otherwise be met by electricity from the grid and hence reduces the electric bill. Time of use (TOU) electricity rates for large general service were applied that mainly consist of energy charges, demand charges, fixed monthly charges, and facilities charges. Summer months comprised of July, August, and September, whereas the remaining months were treated as winter months. Electricity rates are different for summer and winter months. Rate charges are also based upon voltage levels categorized as secondary (≤600 V), primary (levels between 600-25,000 V), and high voltage transmission (>25,000 V). Rates for secondary voltage levels were applied for the DWTP.
The PV system was sized using SAM by using the inputs of desired array size and a DC-to-AC ratio. Using the parameters of the selected module and inverter, the number of inverters and modules were determined. Nameplate capacity was determined as the product of the module's maximum rated power and total number of modules. Some of the equations used to size solar PV are as follows: M parrallel = C I pmax M string (8) M n = M string × M parallel (10) where M string = modules per string, M parallel = module strings in parallel, V mppt-max = maximum MPPT voltage, V mppt-min = minimum MPPT voltage, V max = maximum power voltage, C = array nameplate capacity (kW), I pmax = module maximum power, I n = number of inverters, M mp = Module maximum power, I map = Inverter maximum AC power, and M n = total number of modules. Battery storage was also provided, primarily to meet energy demands at nighttime since the treatment plant operates for 24 h per day. SAM allows the analysis of three types of batteries: lead-acid, lithium-ion, and vanadium redox flow. For this study, a lead-acid battery was used. The battery was assumed to be connected on the AC side. Battery bank size can be specified as either the number of battery cells in a series and battery strings in parallel, or it is estimated using SAM using the inputs of desired bank capacity and bank voltage. (14) where B ns = number of battery cells in series, B dbv = desired bank voltage, B ncv = battery nominal cell voltage, B nbv = nominal bank voltage, B np = number of battery strings in parallel, B cc = battery cell capacity, B nbc = nominal bank capacity, and B db c = desired bank capacity. Nominal cell voltage was the manufacturer reported voltage (volt) of a single cell in a battery, which was used to size the battery storage. C-rate governed the charge and discharge rate of the battery. Max C-rate of charge and discharge of 0.12 per hour was used. SAM utilized a lifetime model that considered battery charge cycles as the main reason for capacity degradation. Capacity degradation was simulated as a function of depth of discharge and number of charge cycles, where depth of discharge described the state of the battery's remaining charge. The capacity losses were incorporated into the maximum battery capacity after each charge cycle had elapsed. SAM does not incorporate the thermal behavior of the battery into the lifetime model; hence, it was assumed that the battery was stored in an air-conditioned room, at constant temperature.

Financial Analysis
Economic performance was analyzed using SAM by employing the parameters of net present value (NPV) and levelized cost of electricity (LCOE). NPV for a system examines the difference between the present value of revenues and costs. A positive NPV means that the present expenses incurred are predicted to provide a return greater than the initial investment. If the NPV is negative, then the cash inflows are predicted to be less than the cash outflows and the investment is not financially beneficial. LCOE was also estimated, which is the present value of the life cycle cost of the PV system. Real LCOE is adjusted for inflation, whereas nominal LCOE is not. Net savings in $ year −1 are estimated by taking the difference between the costs for purchased electricity, with and without the solar system. For this study, NPV was used as the primary metric for financial assessment of the PV system. This is because the analysis includes the input parameters of taxes and incentives, which may complicate the estimation of more simplistic metrics such as payback period. In addition, payback period is not recommended when investment is to be considered explicitly. NPV gives value to the size of investment, while the metrics of LCOE and payback period do not [101].
Various financial parameters incorporated within the economic model of SAM are listed in Table 4, while others are discussed in the following sections. An analysis period of 25 years was used based upon the warranty period of the PV panels.
Net capital cost of the solar system was calculated using SAM, by subtracting cash incentives from total installed costs (sum of direct and indirect cost) ( Table 4). Annual costs incurred for the operation, maintenance, and repairs of the solar system, after it has been installed, constitute the O&M cost ($18 kW −1 year −1 ), which increased each year based on the inflation rate for this analysis. Any savings through the application of state or federal incentives was also taxable. Savings achieved due to the installation of a solar system were not taxable; however, O&M costs were taxable. Debt percentage, which is a fraction of the net capital costs, was specified as 100%, and the debt payments were assumed constant. Principal amount was estimated using SAM as the product of total installed cost and debt percentage. A loan term of 25 years and interest rate of 3.5% were used as inputs to SAM [87]. Salvage values were taken as 20% of installed costs [91].
If land required for installation of a PV system needs to be purchased, the land cost can be incorporated in SAM's financial model. This study assumed that the existing empty land acreage of the treatment plant would be utilized for the installation of the PV system. Hence, land purchase costs were taken as zero.
In SAM, the inflation rate (2.5%) applied to the system costs, electricity rates, and cash incentives, among others, was the yearly mean consumer price index data between the years 1991 and 2012. The real discount rate (8%) is the nominal discount rate minus the inflation rate and is used for the estimation of NPV and LCOE using SAM. In SAM, yearly insurance cost and property tax cost is tax deductible for commercial projects and is a component of the annual operating costs. The insurance rate was taken as 0.25% of the installed costs [88]. For the first year, the insurance amount was a product of the insurance rate and total installed cost, which is then increased by an inflation rate of 2.5% for later years. The property tax amount is the product of the property tax rate and assessed value of the property. The property tax rate was taken as zero to incorporate the property tax exemption incentive.
Depreciation reduces the federal and state taxable income. Different depreciation methods are available. Sometimes the method used to estimate depreciation is dictated by the Federal Government or the state. MACRS (modified accelerated cost recovery system), a depreciation method to help recover capital costs and reduce tax liability, was utilized for this analysis since qualifying systems are those that claim the 30% ITC. Commercial PV is eligible to depreciate over a five-year period [49] as determined by the Internal Revenue Service at depreciation rates of 20%, 32%, 19.2%, 11.52%, 11.52%, and 5.76% for year 1, 2, 3, 4, 5, and 6, respectively. The renewable energy system's taxable basis can be reduced by one half of the investment tax credit of 30%, hence 85% of the taxable basis can be reduced. A bonus depreciation of 50%, 40%, and 30% can be further applied for PV systems placed in service before the year 2018, during 2018, and during 2019, respectively.
Sales tax amount was estimated as the product of sales tax % and direct cost. Since SAM incorporates the sales tax amount as a part of the total installed costs of the solar system, the sales tax amount influences the estimates of depreciation, debt amount, and debt interest payments, as well as the debt interest payment deductible estimate from federal and state income tax. The sales tax is partially exempted based on Nevada state incentives, but the selected DWTP was not eligible because Nevada's property tax exemption was incorporated in the analysis. A sales tax of 8.1% was used.
Scenarios were simulated to reveal the effects of battery storage, governmental incentives, and geographical location on the techno-economic performance of the PV system. A sensitivity analysis was also performed.

Land Area Requirements
Land area requirements for the PV system were estimated using SAM by dividing the total module area (product of a single PV module area and number of modules) with the ground coverage ratio (GCR). GCR is the ratio of the side view diagonal length of one row to bottom length spacing between two rows [102]. The GCR value ranges between 0.01 and 0.99; a value close to zero implies the PV arrays are spaced further apart compared to a GCR value close to 1. For the analysis, a GCR value of 0.3 was selected based on a self-shading analysis [103] for tracking type PV modules. Other than the reason of minimizing self-shading, spacing between the modules is required for access and maintenance work. Land area availability for the development of solar PV was determined by estimating the available landholdings of the DWTP using ArcGIS 9.0 software.

Carbon Emissions
Carbon emissions were estimated by making use of the data provided in Table 2. Nevada's electricity source distribution was used to estimate carbon emissions in the case of the non-PV-based design. Carbon emissions related to the PV-based design were subtracted from those of the non-PV-based design to determine net reduction in the units of carbon emissions per kWh. Carbon emissions generated due to PVs are those generated during its transportation, construction, and disposal phase, but during operation of a solar installation, there are negligible emissions [104].

DWTP Design and Energy Consumption
The results for the analysis are discussed as follows.

Pre-Sedimentation
Typically, to arrive at the plant the water has to be pumped from the water source to the plant. For the selected DWTP, raw river water flows by gravity from the river canal, passing through bar screens, and flows by gravity through the rest of the plant. The bar screens are cleaned manually. Next, the water flows by gravity through the pre-sedimentation basins and the screening facility. Two pre-sedimentation basins (44.5 m × 10 m × 4 m) are provided to avoid the settling of suspended solids in the DWTP's piping. For this study, a surface loading rate of 380 m 3 m −2 day −1 was used. The raw water basins are cleaned manually using wheel loaders. Next, the water flows into the screening facility, which is designed to remove floatable solids from the raw water. The fine-screen opening size is 9.53 mm, and each screen area is 1.4 m 2 . This screen accomplishes 90% of solid removal larger than 4.7 mm. The screens are cleaned using backwashing jets (551.6 kPa) for 15 min day −1 ( Table 5).

Coagulation
Screened raw water flows into the coagulation facility where in-line rapid mixers and a jet diffuser are used to flash-mix coagulant and coagulant aid polymer, respectively. Two in-line mixers are provided. Each 12.5 kW mixer is housed in a 1.5 m diameter pipe and operated for 24 h a day. The jet diffuser pump utilizes a diffuser with a diameter of 0.01 m. The energy consumption was found to be 202.8 kWh day −1 (Table 5).

Flocculation
Downstream of rapid mixing, treated water flows through six flocculation basins where the slow mixing of the coagulated water is achieved through paddle wheels ( Table 5). The detention time for the basins is 25 min. Three stages, with velocity gradients of 70 s −1 , 50 s −1 and 30 s −1 are provided for each flocculation basin of volume 985.8 m 3 .

Sedimentation
Effluent from six flocculation basins flows through submerged openings into six parallel plate sedimentation basins. The peak hydraulic capacity of each basin is 15 MGD, with a volume of 1105.6 m 3 . The flow pattern within the sedimentation basins is across and upward through packs of parallel plates (1.7 m × 0.75 m). Spacing provided between the plates is 60 mm, inclined at an angle of 60 degrees. The Reynolds number was determined to be <1000 to ensure laminar flow, while the Froude number was >10 −5 . Water collected in the troughs at the top of the plate settlers discharges the effluent to the sedimentation basin's effluent channel and flows by gravity to the filters.

Filtration
Following the removal of most of the suspended solids in the sedimentation basins, water is filtered through twelve dual-media filters, including one filter for redundancy. Water enters the filter influent distribution channel from the sedimentation basin, enters each filter and passes through the filter media and underdrain system to the filter flume channel and is finally conveyed to the finished water reservoir. The area of each filter basin is 89 m 2 . Clean bed head loss was determined to be 0.77m. The net available operating head for the DWTP is 3.0 m. The filter is backwashed when the head loss exceeds the available head or limiting head.
Filters are backwashed using both air and water. Air scour is provided to agitate the filter media for about 4 min, before backwashing the filters with water. An airflow rate of 1.4 m 3 sec −1 is used. The filters are backwashed every 24 h, while a filter-to-waste duration of 15 min is used. Backwash water supply is pumped from the plant pump station to the filters by backwash water supply pumps. Recovery is 97%.

Chlorination
Sodium hypochlorite is used to disinfect filtered water at the two finished water reservoirs, each having a volume of 9460 m 3 . The surface water treatment rule is applied for 4 log inactivation of viruses and 1.5 log inactivation of giardia for a CT value of 165 mg L −1 -min at pH 8. Disinfection of cryptosporidium is achieved during filtration, hence, this was not included in the chlorination design. A contact time of 88 min was used in this study. Residual chlorine concentration was 1.9 mg L −1 . A length to width ratio of 47:1 was used for the baffles (hence greater than 40:1 ratio). A soda ash system was designed to prepare and convey soda ash solution to the rapid mixer (0.6 kW), designed for a velocity gradient of 1000 s −1 and detention time of 20 sec. Soda ash feed pumps convey the slurry from the mixing tank to the points of application at the finished water pump station, where it is used to increase pH and alkalinity. Finally, the water is pumped to different areas designated as zone-1, zone-2, and zone-3.

Residual Management
Filtration backwash wastewater and sedimentation sludge generated during filtration and sedimentation is conveyed to a decant basin and water recovery basins, respectively, via a 1.2 m pipeline. Backwash wastewater flows by gravity from filters to two decant basins (37 m × 8 m × 4 m) for primary solids settling. While one is online, the other is used as a backup during maintenance. After a short detention period of 13 min, decanted water is pumped to a raw water junction box, where it is recycled through the plant. Decant basin settled solids are pumped to a water recovery basin (53 m × 8.5 m × 4 m) for secondary solids settling. An additional water recovery basin is provided as a backup basin. The basins were designed for a surface loading rate of 375 m 3 m −2 day, and a detention time of 15 min. Wheel loaders are used for sludge removal. Where sludge drying beds are feasible due to favorable weather, typically filter presses and centrifuge etc. are not used, as in this case.

Energy Consumption of the DWTP
The selected DWTP is unique concerning raw water intake. Raw river water flows by gravity from the river canal into the DWTP; thus, no energy is utilized for water intake. For scenario S-1, the total operational energy consumption was 56.3 MWh day −1 (165.3 Wh m −3 ) for the DWTP including water distribution pumps, whereas energy consumption for scenario S-2, i.e., for the DWTP excluding water distribution pumps was 2661 kWh day −1 (7.8 Wh m −3 ). The energy consumption estimated by the current study (Tables 5 and 6, Figures 2 and 3) was validated by comparing the estimated motor sizes to the plant motor sizes. As shown by Tables 5 and 6, the values are in good agreement with each other. Water 2020, 12, x FOR PEER REVIEW 16 of 30 Figure 2. Results for energy consumption by percentage for water intake, water treatment units, and finished water distribution pumps for scenario S-1 (energy consumption of water treatment facility including distribution pumps) and S-2 (energy consumption of water treatment facility excluding distribution pumps).  Table 5). Pumping operations within the DWTP and the mixing operations (during coagulation, flocculation, and soda ash slurry mixing) utilized about 54% and 40% of the total operational energy consumption of the DWTP, respectively (S-2). Backwashing jet pumps for screens were the smallest consumers of energy within the DWTP, consuming about 0.004 Wh m −3 (Figure 3). Figure 2. Results for energy consumption by percentage for water intake, water treatment units, and finished water distribution pumps for scenario S-1 (energy consumption of water treatment facility including distribution pumps) and S-2 (energy consumption of water treatment facility excluding distribution pumps).  Further, unit process-based quantitative assessment of energy intensities as computed in this study may help to compare and improve the performance of treatment plants. Plants of similar capacities may not be utilizing similar paths of treatment. Hence, the energy consumption may be recomputed based on the unit processes being utilized by the plants, resulting in a fairer energy performance comparison ( Figure 3). The unit process-based quantitative assessment of energy intensities may help compare the performance of different plants more accurately and may identify opportunities for improvement.
Sensitivity analysis was performed for the energy consumption estimates related to the pumping operation for the treatment plant. Wire-to-electric efficiencies of the pumps were increased by 5% and 10%, which resulted in a 2.6% (2.59 MWh day −1 ) and 4.9% (2.53 MWh day −1 ) decrease in the total operational energy consumption (2.66 MWh day −1 ) of the water treatment units, respectively, whereas decreasing wire-to-water efficiencies by 5% (2.74 MWh day −1 ) and 10% (2.82 MWh day −1 ) resulted in increased total operational energy consumption by 2.8% and 6%.
Various studies have been conducted to evaluate the energy consumption of treatment plants. Ref. [62] determined energy consumption for a large-scale water treatment plant in the Netherlands and reported a total energy consumption of 172 Wh m −3 for the treatment plant, including intake pumps and water distribution pumps. Ref. [105] reported that in the U.S. the average electricity use for water supply systems, including water pumping, treatment, and distribution, amounted to about 369.8 Wh m −3 , whereas that of New York State was about 186.2 Wh m −3 .
The methodology developed in this work is applicable to estimating energy consumption for plants treating either surface or groundwater sources. Because the energy computations are based on unit operations, one can simply lay out the units in any desirable configuration. For example, for a groundwater source, typically only filtration and disinfection units would be needed. In this study, the computation does not involve advanced treatment units such as ozonation, membrane filtration, or UV disinfection. This is because these units were not part of the example DWTP, which uses conventional treatment and does not employ/require these advanced techniques for the treatment of their raw water. However, energy considerations of advanced treatment units have been recently published in other studies [106].
Usually, water treatment plants do not have energy consumption instrumentation installed in individual unit processes; rather, they have totalizers for large sessions of the plant. The reason is the high cost associated with instrumentation. Indeed, for this reason most water or wastewater treatment plants do not have this information and cannot benchmark their energy consumption when contemplating switching to renewable sources of energy. In sustainability studies, and in consideration of energy savings, many plants have struggled with lack of data. Even in southwestern U.S., where the plants are relatively new and modern, such data do not exist. The methodology proposed in this study, to overcome the lack of energy consumption data in treatment plants, considers every individual unit of operation present in a treatment plant, and thus equipment-level energy consumption was computed. In addition, the method developed uses motors (considering their respective efficiencies) as the major energy-consuming unit. The design equations proposed, based on typical engineering design parameters and equipment manufacturer data, can be used by any plant to estimate their energy consumption, independent of the existence of electricity consumption meters installed in their unit operations. The need to establish acceptable approaches to compute energy consumption in treatment plants was also recently evidenced by [107].
The plant uses a conventional treatment train (utilizing the processes of coagulation, flocculation, sedimentation, filtration, and disinfection) for water treatment. Other typical treatment trains employed for drinking water treatment are direct filtration (utilizing the processes of coagulation, flocculation, filtration, and disinfection) and inline filtration (utilizing the processes of coagulation, filtration, and disinfection). Thus, the current work covers all the basic treatment processes used by a typical DWTP for water treatment. The design criteria followed by this study are based on industry-accepted design criteria and can be applied with some modifications for any treatment plant. This work is unique because it provides information about the detailed design and the energy consumption evaluation, along with references for replication of this work. Further, the provided information is valuable and meaningful and can help design engineers and researchers alike for application of this work to their treatment plant. The energy consumption evaluation and the potential for solar energy application on plants of various sizes and complexity using the same methodology have been determined by other studies [6,7]. Ref. [6] determined the energy consumption of a plant, which treats 0.425 m 3 of groundwater per second by utilizing the processes of coagulation, filtration, and disinfection (inline filtration). Ref. [7] determined the energy consumption of a treatment plant that treats 1 Mm 3 of raw river water per day by using direct filtration process. Thus, the methodology developed in this work can be applied for design and energy consumption evaluation of other treatment plants.

System Advisor Model
Results achieved using SAM for the PV design and economic analysis are summarized as follows. The techno-economic evaluation was only conducted for scenario S-2, i.e., to offset the energy consumption of the DWTP excluding distribution pumps. A 500 kW PV system was sized, with large storage capacity incurring a net capital cost of $14.5M (Table 7). Monthly performance of the PV system is shown in Figure 4. There is a yearly electricity bill of $6444 (Table 7); since the system is grid-connected, these are the fixed yearly costs for grid connection. Table 7. Results for technical and financial analysis of solar PV system using the system advisor model for scenario S-2 (analysis for energy consumption of water treatment plant excluding distribution pumps) for locations of Nevada and New York.

Effect of Battery Storage
The DWTP operates 24 h a day, hence, to meet the nighttime energy requirements, as well to ensure a balanced supply of energy, battery storage was provided. The battery type chosen was lead acid flooded. Results showed that the battery storage comprised over 80% of the capital costs of the PV system (Table 7). Even though the NPV was positive, the value is small compared to the net capital costs incurred. A higher NPV is better than a lower NPV. A scenario was simulated by

Effect of Battery Storage
The DWTP operates 24 h a day, hence, to meet the nighttime energy requirements, as well to ensure a balanced supply of energy, battery storage was provided. The battery type chosen was lead acid flooded. Results showed that the battery storage comprised over 80% of the capital costs of the PV system (Table 7). Even though the NPV was positive, the value is small compared to the net capital costs incurred. A higher NPV is better than a lower NPV. A scenario was simulated by decreasing the battery prices from $157.7 kWh −1 to $15.77 kWh −1 . This resulted in net capital costs of $3M and a positive NPV for the PV system. The results show that for solar energy to become an attractive prospect for investors, and a reliable source of electricity throughout the 24 h, prices for battery storage must decrease substantially. The DWTP was also analyzed without incorporating battery storage (assuming all other parameters remain the same), for comparison purposes and to reflect the effect of battery storage on cost and performance of the PV system. In this case, the PV system would rely on grid connection to meet the nighttime energy requirements for the plant. It was determined that a PV system without storage has the potential to offset 2.2% (500 kW system, net capital costs $1.7 M, NPV $0.16 M) of the electric load, respectively.
It can be seen that the cost incurred when battery storage was not incorporated was much lower, compared to when storage was provided. The cost associated with battery storage, comprised over 80% of the PV system costs. Successful cost recovery of the PV system with battery storage during the lifetime of the PV would require application of novel approaches for the production of low-cost battery storage systems [108], which would greatly promote the role of solar PV systems as a viable and competitive source of electricity generation for the 24-h duration of the day.
Various studies have explored the use of renewables for water treatment. Ref. [61] compared the performance of a conventional DWTP and a wind-solar-powered nano-filtration pilot plant located in Spain. The study concluded that the application of a nano-filtration process resulted in higher quality water, and usage of renewables not only assisted in making an energy-intensive process sustainable but also helped in making the water treatment process more efficient. Ref. [62] explored the potential of using wind and solar energy for a water treatment plant located in the Netherlands. The study found that about 70-96% of energy independence can be achieved by using 8 MW wind turbines and 5.6 MW solar panels, but due to high costs for large battery storage, 100% energy independence could not be achieved. Solar energy was used in conjunction with wind energy because the site possessed large but erratic wind potential while solar energy using photovoltaics provided a more balanced supply. Ref. [109] determined solar PV to be viable for meeting the electricity needs of two villages located in Bangladesh.

Effect of Governmental Incentives
Financial analysis heavily relies on governmental incentives. Two incentives were applied, federal investment tax incentive and the state property tax exemption incentive for this study for the Nevada site, resulting in a nominal LOCE of 2.65 cents kWh −1 . If either of these incentives were removed from the analysis, the development of solar PV would no longer be profitable as it would result in a negative NPV and increased LCOE. Ref. [82] reported LCOE values for commercial PV systems in the U.S. in the range of 9-12 cents kWh −1 when not incorporating federal ITC and in the range of 6-8 cents kWh −1 when incorporating it. For the current study, when the 30% ITC was not incorporated, the NPV became negative while the nominal LCOE value was estimated as $0.39 kWh −1 .
Based on the renewable portfolio standard (RPS) for the state, 25% of electricity generation by the year 2025 must be achieved using renewables [49]. Portfolio energy credits earned due to the deployment of solar PVs at treatment plants can assist in meeting RPS goals. Incorporating renewables into the existing water infrastructure assists in reducing the dependency of water facilities on the traditional energy sources, thus generating reduced emissions. This helps in making water facilities energy independent, which may lead to cost savings in the long run. Policy makers for water and energy need to participate in collaborative planning and implement policies that provide incentives to existing and future DWTPs for using renewables to offset their energy requirements.

Effect of Location Change
A hypothetical scenario was examined by changing the location of the treatment plant from southwest Nevada (6.3 kWh m −2 day −1 ) to the U.S. east coast. The location selected was New York (3.5 kWh m −2 day −1 ). The financial parameters used are shown in Table 4. The electric rate input, downloaded from the SAM database for electric rates, was based upon the rates of the New York Power Authority for large general service. The results showed that the PV system (500 kW PV, 75 MWh battery) in the Nevada location (Table 7) when located in New York was able to offset 85% of the total load (2661 MWh day −1 ) analyzed. Forty-seven percent of the load offset was achieved using battery storage. It was further determined that a 850 kW PV system using a battery storage capacity of 82 MWh was able to offset 100% of the 2661 MWh day −1 electric load, for the New York location, however, the NPV was negative (Table 7). Fifty-eight percent of the load offset was achieved using battery storage.

Sensitivity Analysis
Sensitivity analysis for the financial parameters that were changed when analyzing for the Nevada and New York locations showed that property tax rates had the largest effect. For the New York location, the NPV changed from negative to positive, when the property tax rate was changed from 2% to 1.2%, respectively. The property tax rate may vary depending upon the location of the property. A property tax exemption or partial property rate incentives have the potential of greatly promoting the development of solar facilities in the state.
Sensitivity analysis was also performed for various financial parameters for the solar system sited in Nevada and showed that a discount rate of 7% or less resulted in a negative NPV. For the current study, a debt percentage of 100% was used. A debt percentage of 96% or less rendered a negative NPV and higher LCOE. Reducing the loan term decreases the NPV and increases the LCOE. Using a loan term of 23 years or less, or loan interest rates of 4% or higher, resulted in a negative NPV. For the current study, an inflation rate of 2.5% was used. Using an inflation rate of 1.9% or less resulted in a negative NPV. It can be seen that for this particular study, the viability of the project is sensitive to the changes in discount rate, inflation rate, debt percentage, loan interest rates, and loan term. For the same PV system, if battery storage is not provided, a debt % of 85% or less, loan term of 18 years or less, discount rate of 4% or less, interest rate of 6% or higher, and inflation rate of zero or less rendered a negative NPV and higher LCOE values.
Ref. [46] determined the economic feasibility of using solar PVs to fulfill the electricity needs of a rural community consisting of 300 homes, located in Gusau, Nigeria. The study found that solar PVs was highly sensitive to changes in electricity rates and inflation rates but slightly sensitive to changes in loan term and loan interest rate. Changes in year-round insolation levels due to change in location may also significantly alter the energy output of PV systems and thus the results of financial analysis as well [46,110]. Ref. [110] conducted a techno-economic assessment of using solar and wind energy with storage systems for locations in the Kingdom of Saudi Arabia (KSA), Canada, and Australia. The study determined the cost of electricity to be lowest for the Yanbu area in KSA (0.609 $ kWh −1 ). Further, the study determined the renewable potential of the KSA to be better than the selected locations in Canada and Australia due to higher year-round levels of solar irradiance.

Land Area Requirements
The land coverage of the PV system for the DWTP (S-2) was 0.01 km 2 ( Table 7). The empty real estate acreage of the DWTP was about 0.4 km 2 , estimated using ArcGIS. Hence, it can be seen that ample land area was available for development of solar PV, and additional costs need not be incurred for the purpose of purchasing land for solar development. Moreover, ample empty land area is available near the DWTP that can be acquired for the installation of PV.
Various studies have analyzed the land area demands of solar installations. Ref. [111] determined that about 15% of the total land area of Iran provided excellent conditions for solar installations. Ref. [112] determined optimal sites in Afghanistan for the development of wind and solar energy systems, resulting in potential generations of 0.34 million GWh per year and 0.15 million GWh per year, respectively.

Carbon Emissions
Net reduction in carbon emissions was found to be 450 metric tons CO 2 -eq year −1 (equivalent to one year emissions for: 95 passenger cars driven, electricity consumption of 67 homes, or carbon sequestration due to 530 acres of U.S. forests). Without battery storage, the net reduction in carbon emissions was found to be 240 metric tons CO 2 -eq year −1 (equivalent to one year emissions for: 50 passenger cars driven, electricity consumption of 35 homes, or carbon sequestration due to 280 acres of U.S. forests). The equivalencies were calculated using the USEPA GHG Equivalencies Calculator [113].
Ref. [114] calculated the reduction in carbon emissions due to residential installations of solar PV in Australia and determined it to be 3.7 million ton CO 2 -eq in the year 2013 and 8 million ton CO 2 -eq in the year 2020. Ref. [115] estimated a GHG reduction of 1.7 ton year −1 for a zero energy solar home in São Paulo, Brazil. Ref. [116] developed a simulation model incorporating estimated future costs for renewables of wind and solar and determined an 80% reduction in carbon emissions by the year 2030 relative to the levels of the year 1990, without any rise in LCOE.
In general, this study evaluated the energy consumption and integration of renewable energy sources into existing water treatment plants for energy sufficiency purposes. Further, the gains achieved in public health and environmental impacts due to reduction in carbon emissions makes the development of solar PVs a good investment [117].

Conclusions
The objective of this study was to utilize solar PVs to meet the energy demands of an existing DWTP treating 3.94 m 3 of river water per second by (a) designing and determining the energy consumption of each unit operation of the treatment plant, (b) sizing the DWTP for solar PVs based on available land holdings and economic analysis, and (c) determining the net reduction in carbon emissions due to the development of solar PVs. The analysis was successfully conducted for scenario S-1 (100% of electric load) and scenario S-2 (5% of electric load, i.e., excluding water distribution pumps), with battery storage and without storage.
For scenario S-1, total operational energy consumption was 56.3 MWh day −1 (165.3 Wh m −3 ) for the DWTP including water distribution pumps, whereas energy consumption for S-2, i.e., for the DWTP excluding water distribution pumps, was 2661 kWh day −1 (7.8 Wh m −3 ). For S-1, water distribution pumps utilized 95% of the total operational energy consumption, whereas the water treatment only operations (S-2) utilized the remaining 5% of the total operational energy consumption of the DWTP. Based on the results of the quantitative analysis of energy consumption, a techno-economic assessment was conducted for solar PVs, to offset the energy consumption of the plant for scenario S-2 only. The development of a 500 kW PV system with battery storage of 75 MWh was found to have a positive net present value of $0.24M, only after the federal and state incentives were incorporated in the analysis. Further, the study showed that reduction in battery prices can tremendously help solar PVs to become a viable source of electricity generation for the entire 24-h duration of the day. Economic feasibility was dependent upon the selection of financial parameters, based on the review of existing literature published between the years 2016 and 2017. Changes in the values of these parameters can greatly affect the results of cost analysis. The economics of solar PVs was found to be sensitive to changes in debt %, inflation rate, loan term, loan interest rates, and discount rates. The real estate acreage of the DWTP was found to be sufficient for the development of solar PVs, which was found to utilize a land area of 0.01 km 2 . The net reduction in carbon emissions because of a solar-based design, with and without the provision of battery storage, was found to be 450 and 240 metric tons CO 2 -eq year −1 , respectively. It was also determined that changing the plant's location from the southwest (higher insolation levels) to the east coast (lower insolation levels) would require a larger PV system size and battery storage capacity and thus increased costs. These results were found to be most sensitive to property tax rate, among the parameters changed between the two locations (sales tax, state income tax, property tax, electric rates). Property tax exemption or partial property tax rate incentives were shown to have promising implications for solar development.
In this study, determination of energy consumption was limited to drinking water treatment operation. Future work may include evaluation of energy consumption related to drinking water distribution within the community.
The design approach utilized for the plant, and equations used, highlighted important design factors that may be adjusted to reduce the energy needs of the plants. Energy intensity values generated rigorously for various processes for the plant may allow a more accurate comparison of a plant's performance against that of other plants. This method can be used for other plants by making suitable alterations. Studies like these can provide the right motivation for increasing the efforts towards a low-carbon water sector or carbon-neutral water infrastructure. Further, the detailed aspects of the design of actual plants are generally not shared publically, and thus this work provides valuable insight for engineering professionals and decision makers. This methodology can be applied to other DWTPs with some alterations for attainment of sustainability goals. The deployment of solar PVs will help achieve the effects of increased public and environmental health and climate benefits.