Characterization of the Outfall Area of a Multi-Stage-Flash Desalination Plant in Bahrain

In Bahrain, like the other Gulf Cooperation Council (GCC) countries, desalination is inevitable to meet the escalating municipal water demands. However, desalination is associated with many environmental effects, which need to be minimized to their lowest possible limits. One of the major environmental concerns of desalination in the Arabian Gulf region is the local and regional effects of the outfall areas on the marine environment. In this study, the outfall area of a government-owned MSF desalination plant is characterized in terms of temperature and salinity. The spatial extent of the plume of the desalination plant’s effluent is mapped by a field survey conducted during the winter season around the plant’s outfall area at 25 cm and 1 m below the water surface and at low and high tide. The results of the characterization indicated that the temperature of the brine discharged to the outfall was 37 ̊C, higher than the ambient seawater temperature by 16.5 ̊C at high tide and 17.5 ̊C at low tide, and that the extent of the mixing zone area was found at about 260 m and 1 km from the outfall point at high tide and low tide, respectively. The results also showed that brine thermal discharge is not in compliance with the standard limits (<3 ̊C from ambient within 100 m of shoreline) both at high and low tides with differences reaching more than 10 ̊C. In terms of salinity, the brine discharged salinity was 56.2 parts per trillion (ppt) compared to an ambient seawater salinity of 43.2 ppt. The maximum salinity measured near the outfall point was 56 ppt at low tide and 51 ppt at high tide, both at 1 m below the surface water column. It is found that the current design structure consisting of two jetties to isolate the desalination plant outfall area from its surroundings is not environmentally sound, as the current surface/inter-tidal outfall location is susceptible to significant increases in salinity and temperature around the outfall area due to the limited flushing it experiences. Therefore, the current design of the outfall area needs to be reviewed to ensure meeting brine discharge regulations and mitigate its impact on the surrounding marine area. The spatial extent of the brine plume can be How to cite this paper: Al-Zubari, W.K., El-Sadek, A.A. and Khadim, M.J. (2018) Characterization of the Outfall Area of a Multi-Stage-Flash Desalination Plant in Bahrain. Journal of Water Resource and Protection, 10, 287-303. https://doi.org/10.4236/jwarp.2018.103017 Received: February 22, 2018 Accepted: March 16, 2018 Published: March 19, 2018 Copyright © 2018 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access W. K. Al-Zubari et al. DOI: 10.4236/jwarp.2018.103017 288 Journal of Water Resource and Protection minimized by building a discharge area further offshore at a sub-tidal location where turbulent flow exists to minimize the spatial extent and intensity of the brine plume. It is recommended that this characterization be extended to all desalination plants in Bahrain, and a regular monitoring program, which should also include selected biological communities and organisms of ecological relevance, be established around the desalination plants outfall areas.


Introduction
Like the rest of the GCC countries, Bahrain has extremely poor endowment of water resources.It has one of the lowest per capita renewable freshwater resources in the world that continue to decline rapidly due to escalating population growth (Figure 1).Overall per capita renewable freshwater availability in Bahrain has been rapidly declining from about 500 cubic meters per year (m 3 /yr) in 1970 to about 80 m 3 /yr in 2010, considerably below the acute water poverty line of 500 m 3 /yr, where water becomes a major constraint for development impacting the standard of living, health and the environment [1].However, despite the rapid increase in water demands and the limitation of its conventional freshwater resources, Bahrain has done well in providing drinking water for its rapidly expanding population by resorting to desalination.Per capita desalination capacity share has been increasing from about 11 m 3 in 1980 to about 170 in 2010.This share continues to increase with time.In 2012 per capita desalination capacity share in Bahrain was about 180 m 3 .Figure 2 illustrates the trend in desalination water production to meet municipal water demands in the country.At present, municipal water supply in Bahrain relies mainly on desalinated water, which is used either directly or blended with groundwater, which has significantly positive implication on the quality of the supplied water.The quality of the supplied municipal water meets the Bahrain/GCC high drinking water quality standards in terms of salinity (<500 mg/L) and beyond.
In Bahrain there are five major operating desalination plants, all of which located at the eastern coast of the country (Figure 3).The details of these desalination plants are indicated in Table 1.In 2014, the total desalination capacity was about 870 m 3 /d (315 Mm 3 /yr).Like the GCC countries, the current trend in desalination expansion in Bahrain is expected to continue in the future [2].
The primary desalination process used in Bahrain is the thermal process      However, while Bahrain has been able to meet the rising municipal water demands by the expansion of desalination plants production, this has been associated with enormous costs manifested by Al-Zubari (2014) [3]: 1) the required energy (oil and gas) for desalinated water production (including its opportunity cost and in-situ value); 2) financial and energy/electricity cost of every stage in the operation of the water cycle system (i.e., production, transmission, and distribution); and 3) environmental costs in terms of thermal brine and other waste products discharge by desalination plants and their impacts on the surrounding coastal and marine environment, as well as air pollution by burned fossil fuel and their impacts on human health and the environment.
Depending on the physical and ecological characteristics of the receiving waters, these substances can have a harmful impact on the local environment.Especially vulnerable are areas such as mangroves, salt marshes, coral reefs, or generally, low energy intertidal areas; enclosed seas, such as the Arabian Gulf, have limited water exchange capacities and are generally shallow and less energetic, thus more sensitive to effluent discharges [4].Discharge of waste products chemicals (i.e., biocides, chlorination, and de-scaling chemicals) can lead to chronic toxicity and small-scale alterations to community structure in marine environments, particularly for corals [5].Moreover, as desalination effluents has high salinity and denser than seawater, they sink to the seabed and slowly circulates causing harm to sea grasses and other ecosystems on which a large range of aquatic life depend [6] [7].Changes to salinity can play a significant role in the growth and size of aquatic life and the marine species disturbance [8].
In light of the current population growth and municipal water consumption trends, expansion in desalination is inevitable to meet the escalating municipal water requirements in Bahrain.Therefore, management efforts have to be made on both the supply and demand sides, to reduce its associated costs to the minimum levels possible.This can be achieved by improving water efficiency (use, supply, recycle, reuse) in the municipal water supply sector to reduce water requirements, and thus reducing desalination production and its associated financial, economic and environmental costs.
Yet, desalination externalities will continue to impact the marine environment and ecology, and will need to be mitigated to the minimum possible levels by technical and technological means.For new desalination plants, basic knowledge of the resulting concentration distributions allows for an impact assessment and design optimization.The concentration distribution depends on the sitting of the outfall, the amount of mixing and the transport capacities of the prevailing currents [9].On the other hand, for an existing desalination plant, a typical environmental impact assessment and mitigation procedures would involve outfall site characterization and modeling approach to improve the design of the discharge such that effluent impacts are minimized.Such approach would require the followings steps: the impact of a given desalination plant on the surrounding marine environment is characterized in the vicinity of the brine discharge area; then, the results of the characterization stage are used to calibrate and develop a hydrodynamic simulation model (e.g., CORMIX) for the desalination plant and its surrounding marine environment; once the simulation model is calibrated to satisfactorily represent the existing system of the desalination plant and its surrounding marine area, it is used to investigate the effectiveness of various proposed mitigation options.
Mitigation options can be made either within the desalination plant itself (e.g., to dilute brine with power plant cooling waters [10], dilute brine with seawater prior to discharge [11] [12], and development of effective anti-scalants with no biological effects [13] [14], or at the outfall area, such as submerged discharge via pipeline and nozzle or diffuser further offshore [15] [16] [17] [18].
In this research, the outfall area of an MSF desalination plant in Bahrain is characterized in terms of temperature and salinity, and the horizontal and vertical extent of the plume of the desalination plant's effluent is investigated.Such basic knowledge of the spatial distribution of the salinity and temperature of the effluent is necessary for an environmental impact assessment, hydrodynamic modeling of the outfall area and for evaluating the effectiveness of any proposed mitigation options.

Location and Description of the Desalination Plant
The study was conducted on the outfall area of Sitra Power and Water Station (SPWS).SPWS's Phase I was commissioned in 1975 and was designed to meet the growing demand for electricity and drinking water in Bahrain.Consisting of four sets of boilers, steam turbo generators and two Multi-stage flash (MSF) seawater desalination units, the facility has a design production capacity of 100 Mega Watts (MW) of electricity and 22,725 m 3 /day (5MGPD) of distilled water.In 1984, phase II of SPWS was commissioned and consisted of a 25 MW gas Turbine with a waste heat recovery boiler, an auxiliary boiler, and a 22,725 m 3 /day (5MGPD) capacity MSF unit.Phase II, or unit 5 is a self-contained plant with separate support infrastructure including: gas supply, seawater Intake, outfall culvert, and other Auxiliary System.In order to meet further growth in demand for electricity and water, Phase III was commissioned in 1984/85.Phase III consisted of three identical MSF units, each with an installed capacity of 22,725 m 3 /day (5MGPD) when operated at 90˚C Top Brine Temperature (TBT), but capable of producing 30,906 m 3 /day (6.8MGPD)When operated at 110˚C (TBT).
SPWS has two seawater intakes and four outfall culverts.A flow diagram of the intake and outfall is shown in Figure 4.The outfall area is being constrained by two jetties.These Jetties, shown in Figure 5, were constructed to minimize the spread of brine plumes and thus minimize areas of ecological impacts, and in the case of the north jetty to limit the intrusion of brines into seawater intake areas [19].

Field Survey and Data Acquisition
The locations of the sampling points are shown in Figure 5. Overall eighteen sampling points were used to characterize the outfall area of SPWS.The locations of the sampling points were determined using a GPS instrument (Garmin).
Temperature (in degree centigrade) and salinity (in total dissolved Solids (TDS)) Journal of Water Resource and Protection were measured at these spot sampling points in the field using Marine Water Quality Monitor (YSI) instrument.The design of the sampling points was made in the form of a grid to ensure full coverage of the spatial distribution of the outfall plume between the two barriers.Sampling points X1 to X16 represent the outfall area, while both points X17 (desalination plant feed-water side) and X18 (outfall side) represent the ambient conditions of the area.
The outfall area and its surrounding are relatively shallow and does not exceed 1.5 meters, with some areas appear at the surface during low tide (Figure 6).Hence, it was felt that it is important to investigate the changes in the salinity and temperature during low and high tides.Sampling and measurements were taken at 25 cm and 1 m below the surface of the water column during both high and low tide (1.65 m and 0.96 m, respectively).The sampling for low and high tide was made on the same day (21/02/2013, 08:30 am and 14:42 pm).Table 2 displays the measurement results of the survey during high and low tide at 25 cm and 1m below the surface of water column for the 18 sampling points.

Spatial Data Representation
Data interpolation and contouring was made using the inverse weighted distance  and the 25 cm below the surface water column planes for both low and high tides were also prepared by subtracting the produced maps using the same software.

Temperature
The ambient temperature of the study area was 20.5˚C at high tide and was 19.5˚C at low tide.The temperature of the brine discharged at the outfall was 37˚C; i.e., the temperature of the brine discharged to the outfall was higher than the ambient temperature by 16.5˚C at high tide and 17.5˚C at low tide.(Figure 7 and Figure 8) show the spatial distribution of temperature at 25 cm and 1 meter below the water surface, respectively, for both high and low tide.In the vertical direction, temperatures are higher at 1 m below surface water column than at 25 cm.This is illustrated in Figure 9, which shows the temperature difference between the 1 m and the 25 cm below the surface water column Journal of Water Resource and Protection planes.The maximum difference temperature was 5.5˚C at low tide and 5˚ at high tide.
The measured temperatures were compared with mixing zone water quality standards of the Kingdom of Bahrain, which states that there should be no direct heat addition within 100 m of shoreline and no thermal alteration, which would cause temperature to deviate from ambient temperature by more than 3˚C (Environmental Legislative Decree 21 for the year 1996).The survey results showed that brine thermal discharge is not in compliance with the standard limits at high tide and low tide (refer to samples 1, 2, 5, and 6 in Table 2).At low tide, the difference reached more than 10˚C.It is observed that the impact of the temperature of the brine discharge at low tide is generally more than that at high tide.This is attributed to the relatively shallow area of the discharge zone, where mixing becomes even more restricted at low tide.

Salinity
The average ambient salinity at high tide and low tide was 43.2 ppt.The brine  salinity discharged at the outfall was 56 ppt; i.e., the salinity of the brine discharged to the outfall was higher than the ambient salinity by 13 ppt.Figure 10 and Figure 11 show the spatial distribution of TDS at high and low tide for 25 cm and 1m sampling depths, respectively.The highest TDS value in the surrounding area of the brine discharge occurred very close to the mouth of the outfall, with the maximum salinity measured near the outfall point was 56 ppt at low tide and 51 ppt at high tide, both at 1m below the surface water column.The relatively high increase in the salinity at low tide than at high tide is attributed to the shallow areas of the outfall at during low tide which leads to less mixing and dilution of the brine.The two figures indicate that, as was the case of the temperature, at high tide, the salinity plume is concentrated within the discharge area, while at low tide the salinity plume is spread towards the open sea.
Vertically, the two figures indicate that the salinity of the seawater is higher at 1m below surface water column than at 25 cm.This can be clearly observed in Figure 12, which shows the salinity difference between 1 m and 25 cm below the  surface water column.The maximum difference in salinity at low tide was 7.5 ppt and it was 4.4 ppt at high tide.This difference in salinity at different level is attributed to the brine density is higher than seawater density leading it to sink toward the sea bed.
Bahrain has no regulation for the salinity at the outfall area (i.e., Environmental Legislative Decree 21 for the year 1996), and compliance analysis cannot be made.However, if compared with the standards of other Gulf countries, for example the Emirate of Dubai, where the standard is set to have the salinity not to exceed 5% of the ambient seawater salinity within 100 m from the outfall area, then this criteria has been exceeded at all times.

Conclusions and Recommendations
The outfall area of the MSF desalination plant of Sitra Power and Water Station was characterized in terms of the spatial distribution of the brine salinity and  The current design of the outfall area needs to be reviewed to ensure meeting brine discharge regulations and mitigate its impact on the surrounding marine area.In order to do so, optimized high efficiency mixing designs are needed for the brine discharge as part of a sustainable concentrate management plan.The current surface/inter-tidal location is susceptible to significant increases in salinity and temperature around the outfall area due to the limited flushing it experiences.The spatial extent of the brine plume can be minimized by building a discharge area further offshore at a sub-tidal location where turbulent flow exists to minimize the spatial extent and intensity of the brine plume.Low Tide Salinity difference between 1 m and 25 cm below Sea level

Figure 1 .
Figure 1.Trends of availability of annual per capita renewable freshwater and desalination capacity in Bahrain, 1970-2010.
renewable freshwater per capita annual desalination capacity W. K. Al-Zubari et al.DOI: 10.4236/jwarp.2018.103017289 Journal of Water Resource and Protection Desalination technology was introduced in Bahrain in 1975 and has developed very rapidly to counteract the shortage and quality deterioration in groundwater resources and to meet the qualitative requirements for drinking water standards.

(
64%), namely the Multi-Stage-Flash (MSF) distillation and the Multiple Effect Distillation (MED) technologies.The MSF technology is an established technology and is combined with co-generation of electricity which greatly improves the economics of desalination.It also exhibits significant economies of scale which are critical for large scale production.In addition, MSF plants have a useful life of about 25 years that can be nearly doubled through proper plant maintenance and refurbishment.The MED combined with thermal vapor compression, is more energy-efficient even for smaller desalination plants than MSF, and has been increasing in Bahrain as well as in the Arabian Gulf region in the past few years.On the other hand, Reverse Osmosis (RO) technology, both seawater and brackish water, have been adopted recently for relatively large plant (i.e., Al-Dur Seawater RO).

Figure 2 .
Figure 2. Trends in municipal water supply and its sources in Bahrain, 1980-2011 (data source: Electricity and Water Authority).

Figure 6 .
Figure 6.Discharge area of SPWS located at the intertidal zone.
The highest temperature value in the surrounding area of the brine discharge occurred, as expected, very close to the mouth of the outfall; the maximum temperature measured near the outfall point at high tide was 30.3˚C at 1 m below surface water column and was 36.3˚C at low tide.The relatively high increase in temperature at low tide than at high tide is attributed to the extreme reduction of the water column during low tide restricting mixing.The figures indicate that high temperatures are concentrated toward the discharge point at high tide, while at low tide the high temperatures spread toward the open sea.The extent of the mixing zone area was estimated at about 260 m from the outfall point at high tide, reaching about 1km at low tide.

Figure 7 .
Figure 7. Observed Temperature at SPWS outfall 25 cm below sea level, in degrees Centigrade: a) High tide; and b) Low tide.
at 25 cm below Sea level W. K. Al-Zubari et al.DOI: 10.4236/jwarp.2018.103017297 Journal of Water Resource and Protection

Figure 8 .
Figure 8. Observed Temperature at SPWS outfall 1 m below sea level, in degrees centigrade: a) High tide; and b) low tide.

35
Journal of Water Resource and Protection

Figure 9 .
Figure 9. Observed Temperature Difference at SPWS outfall between 1 m & 25 cm below sea level, in degrees centigrade a) High tide; and b) Low tide.

5
Journal of Water Resource and Protection

Figure 10 .
Figure 10.Observed TDS at SPWS outfall 25cm below sea level, in ppt: a) High tide; and b) Low tide.
Journal of Water Resource and Protection

Figure 12 .
Figure 12.Observed TDS Difference at SPWS outfall between 1m & 25cm below sea level, in ppt: a) High tide; and b) Low tide.

Table
. Daily production capacity of desalination plants in Bahrain.

Table 2 .
Temperature and Total Dissolved Solids Measurements at high tide and low tide at 25 cm and 1 meter below sea level.