Raed Bashitialshaaer*, Kenneth M. Persson
Department of Water Resources Engineering, Lund University, PO Box 118, SE-221 00 Lund, Sweden
Abstract
In this paper water and salt mass balances for the Dead Sea were modeled. Precipitation, evaporation, river discharges, groundwater flows, input/output from potash companies and salt production, and brine discharge were included in the models. The mixing time in the Dead Sea was modeled using a single-layer (well-mixed) a two-layer (stratified) system. Using the single-layer approach the water level was predicted to change from 411 m below mean sea level (bmsl) (in 1997) to 391 m and 479 m bmsl (in 2097) based on water mass balances including and excluding brine discharge, respectively, and to reach 402 m and 444 m for the two cases based on a salt mass balance. In the two-layer approach, the water level after 100 years was predicted to change from 411 m bmsl (1997) to 397 m and 488 m for a water mass balance including and excluding brine discharge, respectively, and to reach 387 m and 425 m for the two cases using a salt mass balance. The water mixing time using the single-layer description increased from 58 to 116 years when excluding brine discharge. Using the two-layer approach the exchange or mixing time increased in both layers when adding brine discharge to the system, from 1.2 to 1.7 years and 11 to 15.3 years in the upper and lower layers, respectively. Good agreement was found between the models and historical data.
Keywords: Water-Salt balance; Red Sea-Dead Sea Canal (RSDSC); Single-Layer and Two-Layer system; mixing time.
Introduction
The Dead Sea is a salt lake in southwestern Asia. It is bounded on the west by Israel and the West Bank and on the east by Jordan, forming part of the Israeli, Jordanian, and Palestinian Authority border (Fig. 1). The Dead Sea is fed mainly by the Jordan River, which enters from the north. Several smaller
streams also discharge into the Dead Sea, primarily from the east. The Dead Sea has no outlet. Fresh water is transported away by evaporation, which is rapid in the hot desert climate. As a result of large-scale projects by both Israel and Jordan to divert water from the Jordan River for irrigation and other water needs, the level of the Dead Sea has been falling for at least the past 50 years (Gaster, 1997–2006). The level of the Dead Sea was 408.5 m below mean sea level (bmsl) in 1995, being the lowest water surface on Earth (Al-Weshah, 2000). At the end of 1997, the water level was 411 m bmsl and the surface area 640 km2, (Gavrieli, 1997 & 2000). At the beginning of the last century, the level was about 390 m bmsl and had a surface area of 950 km2. In 1966, the Dead Sea covered an area of 940 km2 (76% of
the lake was in the northern basin), and its total length was 76 km, with an average width of 14 km. The total volume of the water in the Dead Sea was estimated to be 142 km3 (only 0.5% being in the southern basin). The surface area continues to decrease due to the high rate of evaporation and decreasing water inflow. The Dead Sea is located in the northern part of the Great Rift Valley. On the east side of the lake the Moab plateau rises about 1,340 m above sea level, and on the west side, the Judea plateau rises to approximately half that height (Asmar, & Ergenzinger, 2002a). *Raed Bashitialshaaer, Tel.: +46 46 222 4367 Fax: ++46 222 4435; E-mail: Raed.Alshaaer@tvrl.lth.se © 2010 International Association for Sharing Knowledge and Sustainability DOI:10.5383/swes.02.02.001 Originally, the Dead Sea consisted of two basins (a northern and a southern basin) divided by the Lisan Peninsula, which extends from the eastern shore out into the lake. The maximum depth in the northern basin was 749 m bmsl, whereas the average depth in the southern basin was only 10 m (Asmar & Ergenzinger, 2002a). The level of the Dead Sea has fallen continuously since the early 1930s, at increasing rates (0.5, 0.7, 0.90, and 0.95 m/y reported in (Steinhorn, & Gat, 1983; Al-Weshah, 2000; JRV, 1996a; Abu-Jaber, 1998), respectively). A summary of
estimated rates of evaporation as a function of salinity in the Dead Sea is presented in Table 1, based on a compilation for the period between the late 1950s and 1980 (Levy, 1984b). The salinity (Total dissolved solids, TDS) increased during this period from 225 g/l to 279 g/l. During the past ten years, the level of the Dead Sea has fallen by more than 25 meters and reached a level of about 416 m bmsl in 2003. Thus, the system is not presently sustainable and various projects aimed at supplying seawater to the Dead Sea have been suggested. The Red Sea-Dead Sea Canal (RSDSC) is being considered as one alternative.
The proposed layout, which may include a desalination plant, is shown in Fig. 1. The inflow of seawater (or reject brine from the desalination plant) into the Dead Sea will have a major impact on its limnology, geochemistry, and biology. The inflow of seawater from the proposed project (RSDSC) (density ≈ 1030
kg/m3) to the highly saline Dead Sea (density ≈ 1240 kg/m3) can be approximated by the inflow of fresh water to the seawater. The impact of such an inflow can be estimated by comparison with the inflow of fresh water during the rainy winter of 1991/2 when the level of the Dead Sea rose by 2 meters and the surface brine was diluted by up to 30% (Beyth et al., 1993).
Table 1. Estimated evaporation rates as a function of surfac salinity in the Dead Sea (Levy, 1984b)
TDS g/l | Level mbmsl | Period | Datasource | Evapt. rate(m/y) | |
225 | 393 | 1942/46 | Neumann,1958 | 1.70–1.75 | |
240 | 395 | 1959/60 | Neev,1967 | 1.47–1.65 | |
256-279 | 401 | 1979/80 | Anatiet al.,1987 | 1.30–1.54 |
When the runoff goes into the Dead Sea, the salt content can fall from its usual 35% to 30% or lower (Salt Works, 2001-2006). Under normal conditions, the salinity of the Dead Sea is about nine times greater than the average salinity of the oceans. Depending on the season, the salinity of the uppermost 35 m of the Dead Sea ranges between 30% and 40%, and the temperature between 19°C and 37°C. Below a transition zone, the bottom layer of the Dead Sea consistently has a temperature of about 22°C, with complete saturation of halite (or sodium chloride, NaCl). Because the water is saturated, salt precipitates out of
solution onto the sea floor. The Dead Sea is presently also saturated or oversaturated with respect to aragonite (CaCO3) and anhydrite (CaSO4) (Gavrieli et al., 1989). Mixing between the calcium-rich Dead Sea brine and sulfate-rich seawater will result in gypsum precipitation (CaSO4·2H2O) (Gavrieli et al., 2005). The Dead Sea salts are commercially important and their production increases the evaporation of water. Wisniak (2002) presented a description of the properties of the Dead Sea together with a chemical analysis of the processes utilized to exploit it commercially. Several numerical models have been developed to determine the water balance (Al-Weshah, 2000; Abu-Jaber, 1998). Horizontal variations in temperature and salinity may occur in the Dead Sea (Krumgalz et al., 1995 & 2000), but these are usually neglected in modeling. Due to the particular chemical composition of the water, the accepted definition of “salinity” is not useful for the Dead Sea and it has been replaced by an equivalent salinity (Anati, 1999), referred to as ‘‘quasi-salinity’’. Carnallite (KMgCl3.H2O) is next to crystallize out, however, this is only
expected to happen when the brine reaches a specific gravity of 1300 kg/m3 (Steinhorn & Gat, 1983). Table 2 presents a comparison of elemental analyses of water from the Dead Sea, the Jordan River, and the Mediterranean Sea based on more than 40 years of data. Opinions vary slightly as to which salt content the water of the Dead Sea has. Gavrieli et al., (2005) reported a density and salinity of the Dead Sea of about 1237 kg/m3 and 342.4 g/l, respectively, whereas Vengosh and Rosenthal (1994) reported the Dead Sea salinity to be about 332.06 g/l. Mixing of the water in the Dead Sea due to low waves and wind is slow compared to that in other water bodies e.g. Seas and Oceans. The Dead Sea may thus be considered a stratified water body, based on 44 available data sets on potential temperature, quasi-salinity, and potential density. These data indicate that the depth of the upper layer of the Dead Sea is about 10% of the maximum depth; defined as the average between June 1998 and December 2007 at the Ein-Gedi 320 station (see Figure 2) (ISRAMAR, 2009).
The chemistry of the Dead Sea
When the runoff goes into the Dead Sea, the salt content can fall from its usual 35% to 30% or lower (Salt Works, 2001-2006). Under normal conditions, the salinity of the Dead Sea is about nine times greater than the average salinity of the oceans. Depending on the season, the salinity of the uppermost 35 m of the Dead Sea ranges between 30% and 40%, and the temperature between 19°C and 37°C. Below a transition zone, the bottom layer of the Dead Sea consistently has a temperature of about 22°C, with complete saturation of halite (or sodium chloride, NaCl). Because the water is saturated, salt precipitates out of
solution onto the sea floor. The Dead Sea is presently also saturated or oversaturated with respect to aragonite (CaCO3) and anhydrite (CaSO4) (Gavrieli et al., 1989). Mixing between the calcium-rich Dead Sea brine and sulfate-rich seawater will result in gypsum precipitation (CaSO4·2H2O) (Gavrieli et al., 2005). The Dead Sea salts are commercially important and their production increases the evaporation of water. Wisniak (2002) presented a description of the properties of the Dead Sea together with a chemical analysis of the processes utilized to exploit it commercially. Several numerical models have been developed to determine the water balance (Al-Weshah, 2000; Abu-Jaber, 1998). Horizontal variations in temperature and salinity may occur in the Dead Sea (Krumgalz et al., 1995 & 2000), but these are usually neglected in modeling. Due to the particular chemical composition of the water, the accepted definition of “salinity” is not useful for the Dead Sea and it has been replaced by an equivalent salinity (Anati, 1999), referred to as ‘‘quasi-salinity’’. Carnallite (KMgCl3.H2O) is next to crystallize out, however, this is only
expected to happen when the brine reaches a specific gravity of 1300 kg/m3 (Steinhorn & Gat, 1983). Table 2 presents a comparison of elemental analyses of water from the Dead Sea, the Jordan River, and the Mediterranean Sea based on more than 40 years of data. Opinions vary slightly as to which salt content the water of the Dead Sea has. Gavrieli et al., (2005) reported a density and salinity of the Dead Sea of about 1237 kg/m3 and 342.4 g/l, respectively, whereas Vengosh and Rosenthal (1994) reported the Dead Sea salinity to be about 332.06 g/l. Mixing of the water in the Dead Sea due to low waves and wind is slow compared to that in other water bodies e.g. Seas and Oceans. The Dead Sea may thus be considered a stratified water body, based on 44 available data sets on potential temperature, quasi-salinity, and potential density. These data indicate that the depth of the upper layer of the Dead Sea is about 10% of the maximum depth; defined as the average between June 1998 and December 2007 at the Ein-Gedi 320 station (see Figure 2) (ISRAMAR, 2009).
Results of the previous studies
The decreasing water level in the Dead Sea has been of concern for many years, and the idea of diverting seawater from the Mediterranean Sea or the Red Sea has been discussed several times. Two scenarios for the RSDSC project were studied by Al Weshah (2000) using a water balance model for the coming 50 years. The first scenario assumes a diversion capacity of 70 m3/s and shows that the Dead Sea will recover its design level of 395 m bmsl in about 40 years. The second scenario assumes a diversion capacity of 60 m3/s and predicts a level of 400.5 m bmsl after about 40 years. The annual net inflow to the Dead Sea
from the River Jordan falls from 175 to 60 millions cubic meter (MCM) due to the fact that Jordan is permitted to utilize 20 MCM of the river water directly and to store another 20 MCM in Lake Tiberias during the winter for consumption in the summer months (JRV, 1996a). Gavrieli and Bein (2006) assumed a diversion capacity of 60 m3/s and a starting elevation in 1995 of about 408.5m bmsl. After 40 years of discharge from a reverse osmosis plant, the new elevation was 400.5 m bmsl (Gavrieli & Bein, 2006). In
another scenario, it was assumed that there was no RSDSC project, and the model showed that the level of the Dead Sea would fall to 444.4 m bmsl. Asmar and Ergenzinger (2002b) modeled a period of 100 years between 1989 and 2088 assuming that the current conditions would prevail. However, the industrial intake increased in the year 2000 by a projected amount of about 25%, leading to a change in water level from 406 to 460m bmsl during the study period in question.
Objectives of the current study
The extreme decrease in the Dead Sea water level, which is mainly caused by the reduced inflow from the Jordan River, high evaporation, low precipitation, and the intake of water to potash companies around, has resulted in substantial damage to the development of the area around the Dead Sea. The RSDSC project has been considered to be one of the most important elements in a development strategy to provide desalinated drinking water to three countries in the region. The brine discharge from this project is an important alternative to fill the Dead Sea that may significantly raises the level. The aim of the present study was to investigate methods for understanding the variations of water level and volume of the Dead Sea. However, it is also of interest to predict the impact of the water discharge depending on release point (discharge location for each term) and density values in each term of input and output. Another task is to evaluate the mixing time (i.e. theoretical time for Dead Sea surface water to mix with the entire water body) for the Dead Sea depending on density and salinity of the input and output. It roughly expresses the amount of time that is needed for waters introduced into the Dead Sea to
be totally mixed. Two models were developed, namely for a single-layer (well-mixed) system and for a two-layer (stratified) system, to estimate the mixing time when taking changes in the Dead Sea level and surface area into consideration.
Overview and background
The principal objective of the RSDSC project was initially to provide desalinated drinking water for the inhabitants of the areas surrounding (Palestinian Authority, Israel, and Jordan) (JRV, 1996b). In contrast, the Mediterranean-Dead Sea Canal, first studied in 1973, focused on the generation of electricity and
this project was intended to produce about 800 MW during peak hours [26] (Mediterranean-Dead Sea co., 1984). At the end of the 1980s and the beginning of the 1990s, the major development goals for the region were re-evaluated. A project was initiated by the Israeli Ministry of Energy and Infrastructure involving hydrostatically driven desalination with reverse osmosis (RO) membranes to produce fresh water from Red Sea water (Tahal, 1994; Beyth, 2002). The main objectives of the revised project are to provide a sustainable source of fresh water for the region (Palestinian Authority, Israel, and Jordan) and to halt the fall in the level of the Dead Sea. The total available hydrostatic head near the Dead Sea end of the proposed canal ranges from 331 ton 545 m, depending on the alignment and the pumping head provided at the Red Sea end (Fig. 3) (Al-Weshah, 2000). The capacity of such a conveyance system has been studied for different diversion flow rates ranging from 40 to 70m3/s (JRV, 1996a).
Fig. 3: Schematic layout of the Red Sea-Dead Sea Canal (Al Weshah, 2000).
The desalination plant is proposed to produce fresh water using the RO method. The estimated annual fresh water production of this plant depends on the capacity of the conveyance system. In one design it could be as high as 850 MCM annually, where the design and economy are limited by the diversion capacity. It is expected that two-thirds of this production will go to Jordan and the rest will go to Israel and the Palestinian Authority (JRV, 1996a). The total dissolved solids (TDS) from the desalination plant are expected to be between 200 and 300 mg/l. It is planned to restore the level of the Dead Sea to about 400 m bmsl after approximately 50 years. Table 3 summarize some properties of the Dead Sea and its output/input, where the river inflow represents the total water received from rivers, springs, and infiltration. Dead Sea level was 390 bmsl, surface area 950 km2 and volume of 155 km3 according to the year 1997 and precipitation was found between 70-90 mm and evaporation between 1300-1600mm/y (Gavrieli, & Bein, 2006; Gavrieli, 2000 and 1997).
Potash and salt works
The salt balance in the Dead Sea is significantly affected by potash and salt exploitation. In the early part of the 20th century, the Dead Sea began to attract interest from chemists who concluded that the lake was a natural deposit for potash and bromine. The Palestinian Authority Potash Company was chartered in 1929 and two more plants were established after 1934; the first plant on the north shore of the Dead Sea at Kalia and the second in the Sodom area south of the Lashon region. Israel produced 1.77 million tons potash, 206,000 tons elemental bromine, 44,900 tons of caustic soda, 25,000 tons of magnesium metal, and sodium chloride from the Dead Sea brine in 2001. On the Jordanian side of the Dead Sea, Arab Potash Company (A-C), formed in 1956, produces 2.0 million tons of potash annually, as well as sodium chloride and bromine. The power plant on the Israeli side allows the production of magnesium metal (by a subsidiary, Dead Sea Magnesium Ltd.).
Environmental impact
One of the major aims of the RSDSC project is to compensate for the negative water balance, which strongly affects the Dead Sea water level (see Fig. 4). Figure 4 also presented the main parameters of the Dead Sea to the left at the beginning of the century, level is 390 bmsl, surface area is 950 km2 and volume of 155 km3 and to the right the level is 411 bmsl, surface area is 640 km2 and volume of 131 km3 is presented at 1997 (Gavrieli, 1997 and 2000). The level recommended in the plane for the Katif alignment was 390 m bmsl (Mediterranean-Dead Sea Co., 1984), and only 12 years later for the RSDSC project it was 400 m bmsl (JVA, 1996c). However, the inflow of seawater or concentrated brines could build up an upper layer of lower density, which may create a hypolimnion with unfavorable environmental conditions and high precipitation of gypsum, at the same time causing changes in the biological environment of the epilimnion (Gavrieli, 1997). The water level has been decreasing over the past 10 years at a rate of about 0.9 m/y, representing an annual water loss of about 600 million cubic meters (MCM). This sharp decrease in water level is due to the withdrawal of fresh water by neighboring countries (over 1000 MCM/y of freshwater), which in the past supplied the Dead Sea. Approximately 200-250 MCM/y of this, corresponding to a fall in water level of about 35 cm/y, is attributed to the activities of the Israeli and Jordanian potash and salt works. These industries together pump 400-450 MCM/y from the Dead Sea into evaporation ponds located in the southern basin, where halite and carnallite (KMgCl3·H2O) are precipitated. At the end of the process, less than 200 MCM/y of concentrated brines (density = 1.35 kg/l; TDS = 500 g/l) are returned to the Dead Sea (Gavrieli, & Bein, 2006).
Methodology
The single-layer system
The system will be described firstly as a single-layer (also called a well-mixed system). The mathematical model used in this study is based on the LOICZ Biogeochemical Modeling Guidelines and has been validated by comparing its performance with other modeling studies of the Dead Sea (Gordon et al., 1996). The model was first employed to describe the dynamic behavior of the Dead Sea using data available in 1997 as initial conditions and simulating the evolution over a 100-year period. Historical data from the period 1976 to 2006 were then used for comparison with simulations obtained with the model. The Dead Sea is not in a steady-state but was assumed to be close to a steady state during the first year. Water and salt balances may have internal inputs and outputs, but this is only of concern in the two-layer approach. A general mass balance for the Dead Sea system yields:
dM/dt = Σ Inputs – Σ Outputs + Σ (sinks)
where dM/dt represents the change in mass of any particular material in the system with respect to time, and Σsinks is the mass of chemical precipitated za per unit time. The amount of salt produced has been found in previous study to be approximately zz0.1 m annually, and this was included as a sink in the calculations (Lensky et al., 2005). The various inputs and outputs are illustrated in Figure 5, where the total mass inflow of water from rivers, springs, and infiltration is denoted by QRi, and direct precipitation enters the system by QP (probably the major source of fresh water in the system). The water used by the Potash Company and other industrial works is denoted QInd-o and the brine discharged by these companies QInd-i (also called industrial outtake and industrial brine input for the system). The brine discharge from the desalination plant (QRO), which will come through the RSDSC project is an important factor in this study. Evaporation represents (QE) a significant output of fresh water, and is considered to be the most important factor causing the decrease in the Dead Sea level. In principle, the net flow (QN), also called the ‘‘residual flow’’ (QN = QRi + QP + QInd-i + QRO – QInd-o – QE- sinks), can be either to or from the system. It is treated algebraically as an input or an output. The equation describing the water budget of the system can now be rewritten as:
(dM/dt)N = Σ(QRO +QP +QRi + QInd-i) – Σ(QInd-o + QE) – Σ sinks
where “(dM/dt)N” is the same as (QN), i.e. the change in the net mass of water with time. This equation describes the water mass budget for the system of interest. It is also possible to write an equation describing the salt mass budget. To obtain the salt budget, an average salinity is assigned to each of the water inputs and outputs. For some of the terms, it is sufficiently accurate to assume that the salinity is zero (e.g., precipitation and evaporation (Lensky et al., 2005)). The equation the salt budget is:
QN SN = Σ(QRO SRO+QP SP+QRi SRi QInd-i SInd-i) -Σ(QInd-o SInd-o+QE SE)-Σ sinks
This equation can be simplified by leaving out the terms in which the salinity is likely to be close to zero:
QN SN = Σ(QRO SRO + QRi SRi + QInd-i SInd-i) Σ(QInd-o SInd-o) – Σ sinks
where, SN = (Ssyst – Save ) / 2, the average salinity of the whole system (SN) is the salinity near the boundary between the output and the input to the system, and it is often adequate to assign this salinity as the average of the adjacent water and the system water. The residual flow (or exchange flow) determines the water mixing time (τ), which is obtained by dividing the system volume (Vsyst) by the absolute value of the net surface flow, which is estimated from the net flow rate through the receiving water as a result of the water mass balance. Thus, the water mixing time is given by:
τ = VSyst / (|QN|)
Mixing time in the Dead Sea is a calculated quantity expressing the mean time that waters of inputs and outputs (having different composition and salinity) spend to let the system return to its characteristics. At its simplest this is equal to the result of dividing the Dead Sea volume by the flow in or out (net flow,
QN) of the system. It roughly expresses the amount of time taken for waters introduced into the Dead Sea to be totally mixed (if possible). This water may be a result of yearly accumulated waters that are added to the Dead Sea, subsequently changing the characteristics of the Dead Sea. It is important that the mixing time of the system is less than one year in order to avoid less dense (e.g. brine water) input from floating in the upper layer, as less dense fluid implies a higher evaporation rate. One
example can be mentioned here. When fish is carried into the Dead Sea by the River Jordan, it dies as soon as the mixing with the saltier waters takes place. This mixing has a specific time. The freshwater doesn’t mix immediately with the saltwater and can float for a considerable time on the surface of the Dead Sea. Until the mixing takes place, the freshwater fish and other biota can continue to live in the floating water. Mixing is interesting.
The two-layer system
The Dead Sea shows relatively strong vertical stratification that can be assumed to resemble a two-layer system (also called a stratified system) (Asmar, & Ergenzinger, 2002a, 2002b & 2003). Such a system is likely to develop in the case of a relatively large freshwater input. Figure 6 illustrates a modified water budget for a two-layer system, which includes entrainment flow (QDeep’), the vertical exchange volume (QZ), and a modified exchange time (τ). The calculations are quite similar to those in the single-layer system, although three further assumptions are needed for the two-layer system: The outflow volume associated with freshwater inputs (QN) occurs in the surface layer, displacing water in the surface layer of the system (SSys-1). River and industrial flows (QRi, QInd-i) enter the deep layer, flow upwards into the surface layer, and out again from the surface layer. The combined outflow from the surface layer (Qsurf + QRi) has positive sign, whereas both QN and Qsurf have negative signs. The salt balance is maintained through a vertical exchange flow (Qz) between the surface and layer and the deep layer (the epilimnion and the hypolimnion). The calculation of QDeep’ is possible only if there is a salinity difference between the inflowing deep water and the surface water of the system, QSys-1, where QSys-1 is about 10% of the total volume of the system. The top layer of the Dead Sea extends approximately 30 m of the total depth of 300 m. Note that this geometry creates a vertical loop circulation, with QDeep’ flowing upwards within the system (it is also called entrainment flow), carrying water of deep system salinity (SSys-2) to the surface, which then spread over the surface. To balance the salt, there is an additional term, the vertical mixing term (QZ), which exchanges surface and deep water within the system. This vertical mixing term only contributes if there is a vertical salinity difference between the surface and deep water. The salt balance of the surface layer is given by:
QROSRO + QPSP + QESE +QRiSRi +QInd-oSInd-o+QDeep’ SDeep’ –Σ sinks +QZ (SSys-2-SSys-1) = 0
To find the vertical exchange volume (QZ) the salinity terms which are equal to or close to zero must be eliminated. The system inflow is balanced by a vertical flow to the surface layer (Qz); and does not enter into the calculations because its magnitude is the same in both vertical directions,
QZ = (-QInd-oSInd-o+QROSRO+QRiSRi –Σ sinks +QDeep’ SDeep’)/(SSys-2-SSys-1)
where QDeep’ is the flow that enters the deep layer, flows upwards into the surface layer, and out again from the surface layer (circulation):
QDeep’ = QN S1 / S Deep’
where, S Deep’ = SSys-2. The water mixing times for the two layers are:
τ1 = V1 / (QInd-0 + |QZ|)
τ2 = V2 / (|QDeep’| + |QZ|)
Results and Discussion
Data for the Dead Sea were excerpted and collected from varying research to assure the accuracy of the results. These data were collected and calculated considering the changes in meteorological and hydrographical parameters in the last 30 years (Gavrieli and Bein, 2006; Gavrieli, 1997 & 2000; Asmar & Ergenzinger, 2002b). This detail of input and output was used for the model validation in the 30-year historical data.
Single-layer versus two-layer
Two different mathematical models were developed to describe the dynamic behavior of the Dead Sea based on the above-discussed models were mass balance equations for the single-layer and two-layer systems. The models were employed using data from 1997 and forecasts were made for a 100-year period. The two modelswere formulated following the LOICZ Biogeochemical Modeling Guidelines and comparing the results with previous studies of the Dead Sea. Based on the current conditions, the models can be used to predict the future behavior of the Dead Sea for different scenarios, including the proposed RSDSC project. The first model employed encompassed a single layer for which the water and salt mass balances were derived. Salinity variations and water discharged from the desalination plant were taken into account with and without the proposed project. Significant variations were predicted between the salinities of the output from and input into the system, especially when RSDSC project is considered. Tables 4 and 5 presents a single box model for a well-mixed system used to calculate the water and salt mass balances, giving the residual flow (QN) and the mixing time (τ) for: a) salinity variations including RO discharge, and b) salinity variations excluding RO discharge.Considering the significant differences in the salinities and densities of the input and output and the Dead Sea itself with respect to depth, a two-layer system was judged to provide a better description of the conditions than the single-layer system. The upper layer constitutes in average of about 10% of the total depth and the rest of the lake constitutes a rather homogeneous lower layer, as shown in Figure 2. The mixing time, (τ), between the two layers, determined by the flow that enters the deep layer (QDeep’) and the vertical exchange volume (QZ), is an important parameter in characterizing the conditions in the Dead Sea. Table 4 shows a comparison between the calculation of the results obtained with the single-layer model and the two-layer model based on previously collected data for 1997 (Gavrieli, 2000). The water and salt mass balances included salinity variations with and without RO discharge. In these calculations, the brine disposal from the desalination plant through the RSDSC was included together with other inputs and outputs (evaporation, rainfall, industrial intake, river inflow, and industrial brine). Values of volume, surface area, elevation, and cumulative levels of the Dead Sea for a 100-year period predicted by the single-layer and the two-layer are presented in Figure 7.
Comparisons between single and two-layer
In a single-layer, over a period of 100 years, the elevation was predicted to change from 411 m bmsl (in 1997) to 391.4 m and 479 m bmsl in the cases with and without RO discharge based on a water mass balance and to reach 401.9 m and 443.6 m bmsl in the two cases based on a salt mass balance. The
discharge from the RO plant affected the level and volume of the Dead Sea significantly. The net residual flow volumes predicted were 292.2 (with RO discharge) and 807.8MCM/y (without RO discharge), and the mixing time increased by 48.2% when RO discharge was not included in the system. The net residual flow volumes obtained with the salt mass balance were 132 (with RO discharge) and 427 MCM/y (without RO
discharge) and the mixing time increased from 58 to 116 years, which is equivalent to about 50% when RO discharge was not included in the system. In the two-layer, after the same period the calculated elevation had changed from 411 m bmsl (in 1997) to 397 m and 488 m bmsl in the two cases (with and without RO discharge, respectively) based on the water mass balance, and to 387 m and 425 m bmsl based on the salt mass balance, respectively. In this model, the exchange times for the upper mixed layer are 1.2 and 1.7 years, including RO discharge and excluding RO discharge, respectively, and 11.0 and 15.3 years for the lower layer with and without RO discharge, respectively. It is found that the exchange time decreased by 29.8% in the upper layer RO discharge was included. When the RO discharge was not included, the exchange time increased by 28% in the lower layer, as can be seen in Table 5. The two-layer model predicted that the amount of (QZ) needed were about 10.56 and 7.28 (x 109m3) with and without RO discharge, respectively, which represents an increase of 31% that will probably increase the exchange time.
Historical data comparison
Comparison was also made between historical data and model results using an initial level of 420.6 m bmsl in 2006. Simulation results covering 30 years, based on water and salt mass balances, using both the single-layer and two-layer models, are presented in Figure 8. The historical data display significant variations during some years. For example, in 1991 and 1992 there was a large amount of rainfall in the Dead Sea catchment area resulting in a great deal of runoff into the Dead Sea. The water and salt mass balances using the single-layer model yielded water levels in 1976 of about 398.4 and 398.5 m bmsl, respectively, compared with the 399 m bmsl from historical data. Employing the two-layer model and the water and salt mass balance, a water level of about 401.9 and 400.9 m bmsl, respectively, were predicted for 1976, compared with the observed level 399 m bmsl. Better agreement was thus obtained when comparing historical data with the single-layer model, with both the water and salt mass balance calculations. Small differences in the results were found with the two-layer model depending on the locations of the output and input and the fact that water with zero or close to zero salinity was excluded, thus omitting the impact of the water with low salinity on the sea level but including all input/output in the water mass balance. Differences may also be caused by uncertainties in the production of the potash company the amount of salts extracted from the Dead Sea and poor control of the input and output in the catchment. The annual amount of salt produced was estimated to be approximately 0.1m based on previous studies. It was very difficult to find the exact amount of the huge rainfall at the years 1991-92. Beyth et al. (1993), they reported for the same period that the amount of 1.5 billion cubic meters of water reached the Dead Sea raising it level by about 2 meters. Therefore, the model results have significant differences in that period because in some the level was added to the Dead Sea was unclear and reported in varying values between 1 to 2 meters. The affected thickness from the additional water formed an upper layer with relatively low salinity that was not more than 20 m deep (see Fig. 4 of Beyth et al., 1993).
Conclusions
Two models, namely a single-layer model and a two-layer model, were used to predict the behavior of the Dead Sea over periods of various lengths. It was shown that the results depend strongly on differences in the salinity of the input and output water, and the possible disposal of brine from a planned RO desalination plant, which is a part of the Red Sea-Dead Sea Canal project. This study included simulations of the conditions in the Dead Sea with the models over a period of 100 years, assuming that the current conditions continue or that brine discharge from the RO desalination plant is added to the system. Mixing times for part of the system with the whole system were significantly different: the two-layer model predicted much lower values than the single-layer model when the RSDSC project was included. Based on the results of the two models, the two-layer description seems to give better results than the single-layer system because the two-layer describes the condition of the Dead Sea. Compared to previous studies, the single-layer and two-layer models proved to be robust alternatives to the traditional water and salt balance techniques, enabling successful calculations of the water exchange through a relatively simple description of a complex and dynamic system such as the Dead Sea. The single-layer model predicted 1.4% and 2% higher water levels than the two-layer model using the water mass balance, with and without RO discharge, respectively, while the two-layer model yielded 3.7% and 4% higher values than the single-layer system using the salt mass balance, with and without RO discharge, respectively. Uncertainties in the amount of salt extracted and the production of other chemicals by industrial works that counted as an outtake from the Dead Sea could be one reason for the discrepancies.