Quantiﬁcation of Submarine Groundwater Discharge in the Gaza Strip

: Gaza Strip has suffered from seawater intrusion during the past three decades due to low rainfall and high abstraction from the groundwater resource. On a yearly basis, more than 170 million m 3 of groundwater is abstracted, while the long-term average recharge from rainfall is 24.4 million m 3 /year. Submarine groundwater discharge (SGD) has never been studied in the Gaza Strip, due to lack of experience in this ﬁeld, next to the ignorance of this subject due to the seawater intrusion process taking place. Continuous radon measurements were carried out in six sites along the Gaza Strip to quantify the SGD rate. The ﬁnal result shows SGD to occur in all sampled sites. The range of SGD rates varies from 0.9 to 5.9 cm · day − 1 . High values of SGD are found in the south (Rafah and Khan Younis governorates). The high values are probably related to the shallow unconﬁned aquifer, while the lowest values of SGD are found in the middle of Gaza Strip, and they are probably related to the Sabkha formation. In the north of Gaza Strip, SGD values are in the range of 1.0 to 2.0 cm · day − 1 . Considering that SGD would occur with the measured rates in a strip of 100 m wide along the whole coast line, the results in a quantity of 38 million m 3 of groundwater being discharged yearly to the Mediterranean Sea along Gaza coast. Nutrient samples were taken along Gaza Strip coastline, and they were compared to the onshore wells, 600 m away from the Mediterranean Sea. The results show that SGD has higher NO 3 − + NO 2 − than nutrient-poor seawater, and that it is close to the onshore results from the wells. This conﬁrms that the source of SGD is groundwater, and not shallow seawater circulation. In a coastal strip of 100 m wide along the Gaza coast, a yearly discharge of over 400 tons of nitrate and 250 tons of ammonium occurs from groundwater to the Mediterranean Sea.


Introduction
Although Submarine groundwater discharge (SGD) may not play a significant role in the global water balance, there are reasons to believe that the geochemical cycles of some major and minor elements may be strongly influenced either by the direct discharge of fresh groundwater into the sea, or by chemical reactions that occur during the recirculation of seawater through a coastal aquifer system [1][2][3]. Groundwater contamination, being a wide-spread problem, SGD may bring pollution to coastal seawater [3,4]. SGD affects the nutrient balance of the sea near the coast, causing harmful algal blooms and changing the flora and fauna of coastal waters [5][6][7][8][9][10].
Several papers have emphasized the important role of recirculated seawater in the transport of solutes from aquifers to coastal water. Moore [2] proposed the term "subterranean estuary (STE)" for the aquifer zone, where recirculating seawater mixes with fresh groundwater, and where water-rock interaction affects the mobility of constituents, including nutrients, towards the sea [11]. Lebbe [12], as well as Werner and Lockington [13], have demonstrated that tidal fluctuations can induce saltwater recirculation in the upper aquifer, producing a saline water table in the coastal fringe.
Measuring SGD by traditional hydrogeological or water balance estimates may be off by several orders of magnitude, largely because of difficulties in constraining hydraulic conductivities. Yet, quantifying SGD is important, as concentrations of dissolved constituents in SGD are often greater than in surface waters, resulting in significant groundwater-derived solute contributions [10,14]. One potential means of evaluating groundwater pathways and fluxes into the coastal zone more accurately is through the use of natural tracers, where 222 Rn is an excellent tracer [3,[15][16][17][18]. The very large enrichment of 222 Rn concentration in groundwaters over surface waters (usually 1000-fold or greater), and its unreactive nature and short half-life (t 1/2 = 3.83 days) make 222 Rn an excellent tracer to identify areas of significant groundwater discharge [3,19]. Continuous radon monitoring (RAD 7-device with RAD-AQUA) could provide reasonably high-resolution data on the radon concentration of coastal seawater at one location over time [20].
In the past three decades, researchers were focusing on the seawater intrusion phenomenon in Gaza Strip, while no one has tackled the reverse flow direction, to establish whether it may occur in the hydrogeological setting. This study focuses on the submarine groundwater discharge from different places along the Gaza coastline using continuous radon measurements with RAD7 and RAD-AQUA devices, beside nutrient analysis along the shoreline and onshore in monitoring wells.

Hydrogeological Background
Gaza Strip is a coastal area along the eastern Mediterranean Sea and lies at latitude 31 • 25 59 N and longitude 34 • 22 34 E. Gaza Strip forms a transition zone between the semi-humid coastal zone in the north, the semi-arid zone in the east, and the Sinai desert in the south. Gaza Strip has a surface area of 365 km 2 where more than 1.8 million inhabitants are living. Groundwater is the only water source for the population of Gaza [21].
The Kurkar Group stratigraphy of Gaza Strip near the coast subdivides the aquifer into four sub aquifers (sub-aquifers A, B1, B2 and C), separated by marine clay layers, whereas in the east, there is only one aquifer ( Figure 1). Those sub-aquifers are semi-confined near the coast, except for the upper layer. The bottom of the aquifer is an impermeable layer known as the Saqiya Group (Miocene-Pliocene age), which is a very thick sequence of marls, marine shales, and clay stones [22]. The Saqiya Group pinches out at around 10 km to 15 km from the coast in the South of Gaza Strip, and the coastal aquifer rests directly on Eocene-age chalk. The saturated thickness of the Gaza coastal aquifer near the coast is around 180 m in the North West, while it is around 40 m at the North Eastern border, and in the South East, it is only few meters (5 m to 10 m; [23]). The Gaza Strip has more than 5000 wells distributed all over the area; 200 wells are used for the water supply system, and the rest are either agricultural or private wells [24].
As a result of aquifer exploitation and imbalance between recharge and abstraction in the past few decades, the groundwater level has dropped to more than 10 m below mean sea level in the southern part of Gaza Strip [21,24]. The yearly groundwater decline rate varies from north to south: in the north, it ranges between 10 to 30 cm, and in the south, it reaches 70 cm [25].

Measurement of Nutrients in Groundwater
Coastal aquifers facilitate the interaction of nutrient-rich groundwater with nutrient-poor seawater. This is especially important in oligotrophic settings, such as the eastern Mediterranean, where the (NO3 − + NO2 − ) content of seawater is typically <100 µM [11,27]. Nutrient sampling of shallow groundwater was done to find indications of SGD, in order to select the best locations for radon measurements. A total of 232 water samples have been collected, 228 of them are groundwater samples taken along the Gaza coastline (51 samples in the north, 68 in the middle and 109 in the south), with a distance ranging from 10 to 20 m from the high-tide line, after digging for 1 to 1.2 m in the beach sand. Four samples were collected from shallow groundwater wells within 600 m from the coastline, which penetrate a few meters in the groundwater (in the south of the Gaza Strip). All collected samples were analyzed for nitrate (NO3 − ), nitrite (NO2 − ), ammonia (NH4 + ), phosphate (PO4 3− ), and salinity. Figure 2A shows the coastal samples' location (232 groundwater samples) and the grouping based on the sample location. Each group of coastal samples were averaged based on location: north, middle, and south of the Gaza Strip.

Measurement of Nutrients in Groundwater
Coastal aquifers facilitate the interaction of nutrient-rich groundwater with nutrient-poor seawater. This is especially important in oligotrophic settings, such as the eastern Mediterranean, where the (NO 3 − + NO 2 − ) content of seawater is typically <100 µM [11,27]. Nutrient sampling of shallow groundwater was done to find indications of SGD, in order to select the best locations for radon measurements. A total of 232 water samples have been collected, 228 of them are groundwater samples taken along the Gaza coastline (51 samples in the north, 68 in the middle and 109 in the south), with a distance ranging from 10 to 20 m from the high-tide line, after digging for 1 to 1.2 m in the beach sand. Four samples were collected from shallow groundwater wells within 600 m from the coastline, which penetrate a few meters in the groundwater (in the south of the Gaza Strip). All collected samples were analyzed for nitrate (NO 3 − ), nitrite (NO 2 − ), ammonia (NH 4 + ), phosphate (PO 4 3− ), and salinity. Figure 2A shows the coastal samples' location (232 groundwater samples) and the grouping based on the sample location. Each group of coastal samples were averaged based on location: north, middle, and south of the Gaza Strip.

Measurement of Nutrients in Groundwater
Coastal aquifers facilitate the interaction of nutrient-rich groundwater with nutrient-poor seawater. This is especially important in oligotrophic settings, such as the eastern Mediterranean, where the (NO3 − + NO2 − ) content of seawater is typically <100 µM [11,27]. Nutrient sampling of shallow groundwater was done to find indications of SGD, in order to select the best locations for radon measurements. A total of 232 water samples have been collected, 228 of them are groundwater samples taken along the Gaza coastline (51 samples in the north, 68 in the middle and 109 in the south), with a distance ranging from 10 to 20 m from the high-tide line, after digging for 1 to 1.2 m in the beach sand. Four samples were collected from shallow groundwater wells within 600 m from the coastline, which penetrate a few meters in the groundwater (in the south of the Gaza Strip). All collected samples were analyzed for nitrate (NO3 − ), nitrite (NO2 − ), ammonia (NH4 + ), phosphate (PO4 3− ), and salinity. Figure 2A shows the coastal samples' location (232 groundwater samples) and the grouping based on the sample location. Each group of coastal samples were averaged based on location: north, middle, and south of the Gaza Strip.  Samples locations were verified, to make sure that sewage outflow to the sea was avoided. Six locations were selected based on high nutrient concentrations to demonstrate SGD using RAD7 with RAD-AQUA for the continuous measurement of radon. Three maps of NO 3 − + NO 2 − , NH 4 + , and salinity were included to show the selection of the SGD locations (see Figure 2B-D).

Measurement of Radon Concentrations in Groundwater
Groundwater is enriched in radon compared to the surface water. While most of the radon ( 222 Rn) pore water is produced from particle-surface-bound radium ( 226 Ra), and the accumulation of the radium is likely regulated by the presence of manganese (hydr)oxides [28].
Radon has greater utility than radium as a SGD tracer when the discharge is fresh or a mixture of fresh and recirculated seawater, since any groundwater in the aquifer (independent of its salinity) is enriched in radon, due to its contact with the sediments and rocks. Radon-enriched fluid is transported to coastal waters, and a radon mass balance can be used to calculate total submarine groundwater discharge [28]. The main sources of uncertainty in the radon mass balance model are associated with: (1) quantifying 222 Rn loss by evasion to the atmosphere [29], (2) quantifying 222 Rn loss via mixing with offshore waters [30] and (3) characterizing the groundwater end-member radon activity that supplied the measured excess Rn in coastal waters [31][32][33].
Since radon is an inert gas, we may expect that the sediment and/or rock uranium/radium content is the main source of radon; it is not produced by chemical reactions or pore water chemistry [28]. The sole source of radon ( 222 Rn with a half time of 3.8 days) in groundwater is the radioactive decay of radium ( 226 Ra with a half time of 1600 years). Radium, the parent of radon, can be retained on grain surface coatings, bound in the mineral lattice of the aquifer matrix, and dissolved in pore water [2]. The short half-life of 222 Rn makes it an excellent tracer to identify areas of significant groundwater discharge [34]. The automated radon system analyses 222 Rn from a constant stream of water produced by a submersible pump and passed through an air-water exchanger, which distributes radon from the running flow of water to a closed air loop. The air stream is fed to the radon-in-air monitor. The radon concentration in the water can easily be calculated from the known temperature dependence of the radon distribution at equilibrium between the air and water phases [20].
Measurements of radon concentrations in the water column have been accomplished by standard oceanographic sampling and analysis techniques (radon emanation) for measurement of 222 Rn taking the special care required for trace gas sampling [3,35,36]. Recently, "continuous" radon monitoring has been described, which can provide high-resolution data on the radon concentration of coastal seawater at one location over time [3,20].
The operation principle of the RAD AQUA device is the closed air loop in which the radon concentration reaches equilibrium during water supply flow. The radon concentration in the water compared to radon concentration in the air is governed by a temperature coefficient of 0.25 at room temperature [3], meaning that four times more radon is present in the air, compared to water. This has a favorable impact on measuring accuracy.
Six locations along the Gaza Strip coastline were measured for continuous radon, to assess the SGD ( Figure 3): two locations in the north, two in the middle, and two in the south of the Gaza Strip.

Field Work
To measure the submarine groundwater discharge using radon as a groundwater tracer, RAD7 and RAD-AQUA were used together to measure continuous radon in water in six locations in the Gaza Strip (Figure 3). A submersible pump was placed at the sea floor (at a depth ranging from 0.83 to 1.05 m), pumping water from the seafloor through a hose to the RAD-AQUA device for spraying, while the RAD7 measures the produced radon from the sprayed water. We initially intended to measure temporal changes of SGD, but due to the situation in the Gaza Strip, we were unable to do this. The time duration for sampling was different in each sampling location according to the field situation and obstacles found in the field. The duration ranged from 1.5 hr to 2 days. In the middle part of Gaza Strip (Deir El Balah governorate), two short measurements were carried out, ranging from 1.5 to 3 hr; due to some practical constraints, we were not able to spend more time for measuring at these sites. Table 1 shows the location and duration of each measurement.

SGD Model from Radon Tracing
Estimating groundwater discharge via radon is based on a mass balance approach. Inventories are measured, either as a snapshot or continuously over time, and these inventories are converted to input fluxes after making allowances for losses, due to decay, atmospheric evasion, and net coastal

Field Work
To measure the submarine groundwater discharge using radon as a groundwater tracer, RAD7 and RAD-AQUA were used together to measure continuous radon in water in six locations in the Gaza Strip (Figure 3). A submersible pump was placed at the sea floor (at a depth ranging from 0.83 to 1.05 m), pumping water from the seafloor through a hose to the RAD-AQUA device for spraying, while the RAD7 measures the produced radon from the sprayed water. We initially intended to measure temporal changes of SGD, but due to the situation in the Gaza Strip, we were unable to do this. The time duration for sampling was different in each sampling location according to the field situation and obstacles found in the field. The duration ranged from 1.5 h to 2 days. In the middle part of Gaza Strip (Deir El Balah governorate), two short measurements were carried out, ranging from 1.5 to 3 h; due to some practical constraints, we were not able to spend more time for measuring at these sites. Table 1 shows the location and duration of each measurement.

SGD Model from Radon Tracing
Estimating groundwater discharge via radon is based on a mass balance approach. Inventories are measured, either as a snapshot or continuously over time, and these inventories are converted to input fluxes after making allowances for losses, due to decay, atmospheric evasion, and net coastal "mixing" terms. Although changing radon concentrations in coastal waters could be in response to a number of processes, the advective transport of radon-rich groundwater (pore water) through sediment is often the dominant process [37]. Thus, if one can measure or estimate the radon concentration in these advecting fluids, the 222 Rn fluxes may be easily converted to water fluxes. 222 Rn is also produced by 226 Rn in the sampled water, and this quantity has to be subtracted to obtain "Excess 222 Rn", which is relevant for our purposes. The inventory refers to the total amount of excess 222 Rn ( 222 Rn in water minus 226 Rn in pore water) per unit area. Decay is not considered because the fluxes are evaluated on a very short (1-2 h) time scales that are relative to the half-life of 222 Rn [38].
The following data are required to perform the analysis: continuous measurement of total 222 Rn activities (dpm·m −3 ) in the coastal water column, water depth measurement, water and air temperatures, wind speed, atmospheric 222 Rn concentrations, and 226 Ra in the coastal water.
RAD7 device measures the total 222 Rn (dpm·m −2 ) in the water column every 30 min. 226 Ra could not be measured during our research; hence, it has been taken from Weinstein et al. [39]. They measured 226 Ra several times (Carmel coast in Israel, 60 km to the north of Gaza Strip) by gamma spectrometry, to show 226 Ra activity in the range of 200-240 dpm/kg in the sand. On the other hand, they found 222 Rn in pore water to be in a range of 340 to 390 dpm·L −1 for the Kurkar groundwater [39].
Due to un-accessibility to measure the 222 Rn activity offshore, we used the value of 1000 dpm·m −3 taken from Peterson et al. [40]. The wind speed at the time of measurement was taken from http://www.wunderground.com.
We used the various assumptions that are inherent in the application of a 222 Rn box model as stated in Burnett and Dulaiova [3], Dulaiova et al. [30] and Burnett et al. [38], to derive a rate of submarine groundwater discharge.

Nutrients
High contents of (NO 3 − + NO 2 − ) in the coastal samples may be representing SGD ( Figure 4A).
The average groundwater concentration is 440 µM in wells, while in the coastal samples, the concentration ranges from 204.6 to 320.1 µM. A trend of decreasing (NO 3 − + NO 2 − ) ( Table 2) is shown from the north to the south of Gaza Strip. Literature states that fresh groundwater concentration for (NO 3 − + NO 2 − ) is around 330 µM [11]. The analyzed concentrations in the coastal samples of our study, thus indeed indicating a substantial contribution of SGD. The average concentration in the sampled groundwater wells was 4.9 µM of PO 4 3− , while the coastal samples were in the range of 5.6 to 11 µM (Table 2 and Figure 4B). Phosphate in the groundwater wells is lower than the measured phosphate in the coastal samples.   The average concentration of NH 4 + in the onshore wells is 360 µM, while the coastal samples are ranging from 362.5 to 380.8 µM (Table 2 and Figure 4C). Hence, the fluxes to the sea have higher NH 4 + concentrations than seawater that is poor in NH 4 + (the ammonia concentrations in seawater vary from <0.1 to 15.3 µM [41]). The nutrient analytical results for the three groups (north, middle, south) were tested against the significance of their difference, and it was found that, for nitrate, the three groups have significantly different averages at the 5% significance level. This was not always the case for the other parameters.

Radon
Seven measurements for continuous radon monitoring using RAD7 and RAD-AQUA were carried out in six locations ( Figure 3 and Table 1).
The extracted data were processed following the Burnett and Dulaiova [3] model. The results processed by the model are summarized in Table 3, while Figure 5 shows the model results for the SGD flux per unit area based on each sample's time series. Figure 6 illustrates the relationship between the ebb/flood tide (water depth) and the SGD rate at the Furosia location, where positive advection fluxes represent groundwater discharge into the Mediterranean Sea and tend to occur during ebb tide.
In the north, three measurements have been carried out, one in Furosia and two in Waha sites. The Waha site should have one continuous measurement but due to pump failure after 9.5 h, we stopped the measurement and continued with a second measurement after replacing the portable submersible pump. The results show 0.92 ± 0.27 dpm·L −1 for radon in water at the Furosia site, which results in 1.4 cm·day −1 as an average SGD flux per unit area for two days of measurements. At the Waha site, radon was in the range of 0.37 ± 0.17 to 0.45 ± 0.18 dpm·L −1 and the produced SGD rate is 1.5 cm·day −1 as an average for both measurements. The average SGD rate in the north of Gaza Strip is 1.5 cm·day −1 .
In the south of the Gaza Strip, higher radon concentrations were found: 3.51 ± 0.67 and 8.2 ± 1.2 dpm·L −1 in Rafah and Khan Younis, respectively, which produced SGD fluxes per unit area of 5.1 cm·day −1 and 5.9 cm·day −1 in Rafah and Khan Younis, respectively (average of 5.5 cm·day −1 ).
The calculated SGD in the middle part of Gaza Strip showed identical values of 0.9 cm·day −1 at both measuring sites, which gave us some confidence that the produced fluxes, in spite of the short time during which they were measured, were still fine. Further monitoring is recommended to enhance the results.
A sensitivity analysis was carried out to check the effect of the value for 222 Rn offshore taken from Peterson et al. [40] (1000 dpm·m −3 ) to the overall SGD rate calculations, by dividing and multiplying the value for 222 Rn offshore by 2, and re-calculating the SGD rate for each measured location (Table 4).

Discussion
This is the first study that has been carried out to quantify submarine groundwater discharge in the Gaza Strip, whereas seawater intrusion has been studied in this area, due to extensive abstraction and low groundwater recharge. Gaza Strip groundwater abstraction amounts to 170 million m 3 [24,42], while the long term rainfall recharge is 24.4 million m 3 /year [21], leading to lowering of the groundwater level by more than 10 m in the south, and 5 m in the north.
Nutrient analysis is a good way to identify potential locations of SGD occurrence, taking into consideration any other nutrient source (i.e., sewage outflow to the sea and/or agricultural activities close to the sea). Seawater is mainly poor in nutrients; hence, any high level of nutrients along the shoreline would be from groundwater. In this study, we used the nutrient analysis to indicate the potential sites for SGD, and to identify potential radon measuring sites. We found that the average (NO3 − + NO2 − ) in all groups had increased values. Also, considering the high NH4 + concentration in coastal samples compared to the seawater to be within the range of the onshore groundwater wells, this allowed us to conclude that SGD was indeed taking place.
Phosphate is as low as 5 µM in the groundwater wells, probably due to phosphate removal by water-rock interaction, while the coastal samples showed higher phosphate concentration in all groups, except the samples in Rafah governorate.

Discussion
This is the first study that has been carried out to quantify submarine groundwater discharge in the Gaza Strip, whereas seawater intrusion has been studied in this area, due to extensive abstraction and low groundwater recharge. Gaza Strip groundwater abstraction amounts to 170 million m 3 [24,42], while the long term rainfall recharge is 24.4 million m 3 /year [21], leading to lowering of the groundwater level by more than 10 m in the south, and 5 m in the north.
Nutrient analysis is a good way to identify potential locations of SGD occurrence, taking into consideration any other nutrient source (i.e., sewage outflow to the sea and/or agricultural activities close to the sea). Seawater is mainly poor in nutrients; hence, any high level of nutrients along the shoreline would be from groundwater. In this study, we used the nutrient analysis to indicate the potential sites for SGD, and to identify potential radon measuring sites. We found that the average (NO3 − + NO2 − ) in all groups had increased values. Also, considering the high NH4 + concentration in coastal samples compared to the seawater to be within the range of the onshore groundwater wells, this allowed us to conclude that SGD was indeed taking place.
Phosphate is as low as 5 µM in the groundwater wells, probably due to phosphate removal by water-rock interaction, while the coastal samples showed higher phosphate concentration in all groups, except the samples in Rafah governorate.

Discussion
This is the first study that has been carried out to quantify submarine groundwater discharge in the Gaza Strip, whereas seawater intrusion has been studied in this area, due to extensive abstraction and low groundwater recharge. Gaza Strip groundwater abstraction amounts to 170 million m 3 [24,42], while the long term rainfall recharge is 24.4 million m 3 /year [21], leading to lowering of the groundwater level by more than 10 m in the south, and 5 m in the north.
Nutrient analysis is a good way to identify potential locations of SGD occurrence, taking into consideration any other nutrient source (i.e., sewage outflow to the sea and/or agricultural activities close to the sea). Seawater is mainly poor in nutrients; hence, any high level of nutrients along the shoreline would be from groundwater. In this study, we used the nutrient analysis to indicate the potential sites for SGD, and to identify potential radon measuring sites. We found that the average (NO 3 − + NO 2 − ) in all groups had increased values. Also, considering the high NH 4 + concentration in coastal samples compared to the seawater to be within the range of the onshore groundwater wells, this allowed us to conclude that SGD was indeed taking place. Phosphate is as low as 5 µM in the groundwater wells, probably due to phosphate removal by water-rock interaction, while the coastal samples showed higher phosphate concentration in all groups, except the samples in Rafah governorate.
Only one method of SGD quantification has been applied in the study area, while other tools could not be carried out to enhance the accuracy of the outcomes. Seepage meters collecting in situ pore water in a closed system, are expected to provide more reliable results on the local scale.
The results show a variation in SGD quantity from the north to the south of Gaza Strip, where the southern part has higher fluxes than the north, while the middle area of Gaza Strip has a low quantity of SGD.
In the northern part of the Gaza Strip, where three continuous radon measurements were performed, the SGD rates were in the range of 1.0 to 2.0 cm·day −1 (average of 1.5 cm·day −1 ). On the other hand, the SGD rates were higher in the south (range of 5.1 to 5.9 cm·day −1 ), while they were minimal in the middle part (0.9 cm·day −1 ). It was calculated ( Table 5) that, considering solely a strip of 1 m wide along the coastline, in the north, a groundwater quantity of 97 × 10 3 m 3 /year is discharged into the Mediterranean Sea, while this is only 31 × 10 3 m 3 /year in the middle part, and a large flux of 247 × 10 3 m 3 /year in the south. The results thus show a variation in SGD quantity from the north to the south of the Gaza Strip, where the southern part shows higher fluxes than the north, while the middle area shows low quantities of SGD. Calculating SGD quantity in a 1 m wide strip along the whole Gaza coastline results in total SGD of 0.38 million m 3 /year, where the SGD quantity in the south is more than two times higher compared to the north. SGD may be expected to occur over a much wider area, although the rates can be expected to decrease with increasing distance from the coast. For this reason, we refrain from taking the whole outcrop area of Kurkar Group at the sea bottom into account, as the measured SGD may be expected to represent the maximum rates, which are not representative for the whole submarine outcrop area. Considering that SGD would occur with this rate in a strip of 100 m wide (which still appears to be a conservative estimate), results in a quantity of 38 million m 3 of groundwater being discharged yearly to the Mediterranean Sea along Gaza coast.
The reason behind having higher SGD in the south compared to the north, while the recharge rate from rainfall is higher in the north than in the south [21], probably refers to the amount of abstraction quantities (in the north, 61.6% of total groundwater abstraction (for domestic water supply) is located, while it is only 16.4% in the middle area and 22% in the southern part of Gaza Strip; [24]). More abstraction in the north leads to less groundwater discharge to the sea. On the other hand, the low SGD rate in the middle part of Gaza Strip (0.9 cm·day −1 on average) is probably related to the presence of Sabkhas, which act as discharge areas for groundwater and intercept the groundwater before it can reach the sea.
In a coastal strip of only 1 m wide, a yearly discharge of over 4 tons of nitrate and 2.5 tons of ammonium to the Mediterranean Sea occurs ( Table 5). Assuming that SGD with the measured rates would occur in a strip of 100 m wide, would result in nutrient discharges that are one hundred times larger than calculated in Table 5, that is, over 400 tons of nitrate and 250 tons of ammonium that are yearly discharged by Gaza Strip groundwater into the sea. The obtained results from using different 222 Rn offshore values (Table 4) showed no significant difference to the overall calculated SGD rate, which gave us confidence to use 222 Rn offshore, as stated in Peterson et al. [40].

Conclusions
Using nutrient analytical results along the shoreline helps us to identify potential sites of submarine groundwater discharge, after eliminating sites where nutrients may be due to other potential sources (such as sewage outflow to the shoreline and/or agricultural practices that are close to the shoreline). In our study, the raised nutrient concentrations in shallow groundwater collected near the shoreline, helped us to confirm that SGD is taking place, and that it can be further investigated by the means of 222 Rn measurements.
Gaza Strip has suffered from seawater intrusion in the past few decades, even though groundwater discharge to the Mediterranean Sea is still occurring, with a range of up to 5.9 cm·day −1 (in the selected sites the minimum SGD rate was 0.9 cm·day −1 ). SGD rate was found to be maximal in the south of Gaza Strip (average of 5.5 cm·day −1 ), while being much lower in the north (average of 1.5 cm·day −1 ) in spite of the substantially higher groundwater recharge in the north. This must be related to groundwater exploitation, which is much higher in the north. In central Gaza, measured SGD rate is minimal (average of 0.9 cm·day −1 ), which must be due to the presence of Sabkhas, acting as groundwater discharge areas, and intercepting the groundwater before it reaches the sea. The sensitivity analysis carried out for 222 Rn offshore showed no major effect on the calculated SGD rate in our study.
The analysis showed a groundwater discharge to the sea of about 0.38 million m 3 /year in a 1 m wide strip along the full length of Gaza coast. Assuming that the same SGD rate occurs in a 100 m wide strip, this results in a yearly SGD quantity of 38 million m 3 /year along the Gaza coastline. This SGD transports over 400 tons of nitrate and 250 tons of ammonium to the Mediterranean Sea.
Regular measurements of SGD will enhance our understanding of water balance in the Gaza Strip, and describe the temporal and spatial variability of SGD along the Gaza coastal aquifer.