E ﬀ ects of Aquifer Bed Slope and Sea Level on Saltwater Intrusion in Coastal Aquifers

: The quality of groundwater resources in coastal aquifers is a ﬀ ected by saltwater intrusion. Over-abstraction of groundwater and seawater level rise due to climate change accelerate the intrusion process. This paper investigates the e ﬀ ects of aquifer bed slope and seaside slope on saltwater intrusion. The possible impacts of increasing seawater head due to sea level rise and decreasing groundwater level due to over-pumping and reduction in recharge are also investigated. A numerical model (SEAWAT) is applied to well-known Henry problem to assess the movement of the dispersion zone under di ﬀ erent settings of bed and seaside slopes. The results showed that increasing seaside slope increased the intrusion of saltwater by 53.2% and 117% for slopes of 1:1 and 2:1, respectively. Increasing the bed slope toward the land decreased the intrusion length by 2% and 4.8%, respectively. On the other hand, increasing the bed slope toward the seaside increased the intrusion length by 3.6% and 6.4% for bed slopes of 20:1 and 10:1, respectively. The impacts of reducing the groundwater level at the land side and increasing the seawater level at the shoreline by 5% and 10% considering di ﬀ erent slopes are studied. The intrusion length increased under both conditions. Unlike Henry problem, the current investigation considers inclined beds and sea boundaries and, hence, provides a better representation of the ﬁeld conditions.


Introduction
Saltwater intrusion (SWI) is a serious and common problem in many coastal aquifers around the globe. The degree and severity of the problem depend on several natural and man-made parameters [1,2]. To investigate SWI problems, density dependent groundwater flow and solute transport models are employed to simulate water flow and concentration pattern in the transition zone between freshwater and saltwater bodies. Under certain circumstances, the saltwater body moves inland, and the problem becomes more pronounced. In severe cases, drinking and irrigation wells might be abandoned due to high water salinity. Over-pumping of groundwater from coastal aquifers disturbs the long-term hydraulic balance between freshwater and saltwater bodies and accelerates inland migration of saltwater [3,4].
The simulation of SWI in coastal aquifers started long-time back by [5,6]. Based on some experimental work, they developed an equation to analyze the problem. They assumed that the freshwater and saltwater are two immiscible fluids and hence a sharp interface will be developed between the two water bodies. This assumption, however, oversimplifies the problem and may not be valid in most cases. In reality, a transition (dispersion) zone is developed between the two water Hydrology 2020, 7, 5; doi: 10 bodies. In the transition zone, the flow of water is under the hydraulic gradient, while the transport of salt ions is governed by the hydrodynamic dispersion processes. For undisturbed coastal aquifers, saltwater intrusion problems can be studied analytically using the sharp interface assumption based on Ghyben-Herzberg relation and the single potential theory [7]. This model depends on the hydrostatic balance between fresh and saline water. In other words, the weight of a saltwater column is balanced with a weight of a freshwater column to locate the interface. The steady state Henry analytical solution [8] has exclusively served as a benchmark for testing SWI models for over five decades [8]. A comparison between the results of different codes applied to Henry's problem was conducted and a sensitivity analysis to identify the most dominant parameters that govern the severity of the problem was performed [9]. An example of the analytical solutions of SWI problems is provided in [10]. A numerical model considering the dispersion zone between freshwater and saltwater was developed [11]. The model was used to predict the effect of saltwater upconning on the salinity of pumped groundwater. Other researchers [12,13] investigated the seawater intrusion in the Nile Delta aquifer in the aerial view under different recharge conditions. The possible effects of climate change and possible seawater level rise on SWI in the Nile Delta aquifer of Egypt and the Madras aquifer of India were investigated using a numerical model [14]. Due to the geometric conditions and its wider exposure to the saline water body of the Mediterranean Sea, the Nile Delta aquifer is more vulnerable to the SWI problem under the conditions of climate change.
A numerical model using SEAWAT code to simulate SWI at Chaouia aquifer in Morocco was developed [15]. A laboratory study was conducted to assess the transport patterns of saltwater wedge in costal aquifers under different conditions [16]. The effect of the changes of groundwater and seawater levels on seawater intrusion in an unconfined coastal aquifer using SEAWAT was investigated [17]. The results showed that increasing seawater level will increase the inland migration of SWI [18] also used SEAWAT to simulate the SWI under different hydraulic parameters in coastal aquifers. The results indicated that the aquifer hydraulic properties have a significant effect on the dynamics and severity of SWI problems and hence might be altered to mitigate and retard the intrusion process.
Analytical methods and a numerical model (SUTRA) were used to study the impact of seawater level rise (SLR) and land-surface inundation (LSI) on SWI considering different recharge rates and variations of land-surface and aquifer bed slopes [19]. The results indicated that theses parameters have significant impact on SWI. A numerical model based on mesh free method to study SWI using Henry's problem was employed for the case of a confined aquifer [20]. They considered the effects of the aquifer hydraulic parameters, boundary conditions, pumping and recharge rates and aquifer depth [21] extended the semi-analytical solution of the dispersive Henry problem based on the Fourier-Galerkin method to simulate heterogeneous and anisotropic coastal aquifers. The results indicated that the modified analytical solution can be effectively used for verification of existing and future numerical models. The method can also be used to investigate SWI under the combined influence of stratification and anisotropy.
According to the available literature, many researchers used Henry problem as a benchmark for calibration of SWI models. The Henry problem assumed horizontal bed and vertical land and sea sides. Such conditions might not be applicable for the majority of the costal aquifers. Moreover, the groundwater and seawater levels were assumed constant. Unlike Henry problem, the current study aims to assess the effect of changing aquifer bed and seaside slops on SWI. The effects of groundwater and seawater levels on SWI are also considered. A better understanding of the SWI problem is, therefore, provided.

Materials and Methods
In this study, the numerical model SEAWAT [22] is used to investigate the effects of different parameters on the seawater intrusion problem in coastal aquifers. The investigated parameters include the seaside and bed slopes, the seawater level, and the freshwater flux from the land side boundary. The results of the model were compared with a number of published cases and was then used to investigate the effects of the different parameters.

Numerical Model
The SEAWAT model is based on coupling MODFLOW and MT3DMS to solve fluid flow and solute transport equations. In this study, the model was validated using Henry problem. Figure 1 presents a vertical cross section of an idealized confined aquifer and the applied boundary conditions. The most popular analytical solution for saltwater intrusion phenomena was developed by [8]. A number of researchers used Henry problem for the verification of numerical models (e.g., [23][24][25][26][27][28]). Typically, the study domain is 200 cm in length, 100 cm in height and 10 cm in width. In this study, the domain was discretized using 4 rows, 80 columns and 44 layers with a cell dimension of 2.5 cm in each direction.

Numerical Model
The SEAWAT model is based on coupling MODFLOW and MT3DMS to solve fluid flow and solute transport equations. In this study, the model was validated using Henry problem. Figure 1 presents a vertical cross section of an idealized confined aquifer and the applied boundary conditions. The most popular analytical solution for saltwater intrusion phenomena was developed by [8]. A number of researchers used Henry problem for the verification of numerical models (e.g., [23][24][25][26][27][28]). Typically, the study domain is 200 cm in length, 100 cm in height and 10 cm in width. In this study, the domain was discretized using 4 rows, 80 columns and 44 layers with a cell dimension of 2.5 cm in each direction. The boundary conditions for head and concentration were set based on the following. At the land side, a hydrostatic freshwater head distribution was considered and a constant flux (Qin) of 0.5702 m 3 /d/m was implemented. The concentration (Cin) was sent constant along the boundary and equals to zero mg/L. At the seaside, a hydrostatic seawater head was considered to represent the density variation with depth at this boundary, with a hydraulic head of 1.0m above the upper boundary of the aquifer. A constant concentration of 35,000 mg/L (seawater concertation) was implemented. The upper and lower horizontal boundaries were set as no flow boundary for both the fluid and salt ions. As such, the stream lines should be parallel and the equipotential lines should be perpendicular to the upper and lower boundaries. The implement hydraulic parameters are presented in Table 1.

Model Verification
In this study, SEAWAT is employed to simulate Henry problem considering the above boundary conditions and hydraulic parameters. As shown in Figure 2a, the intrusion of 0.5 Isochlor (17,500 mg/L) reached a distance of 62.5 cm measured along the bottom boundary from the seaside. The results of this simulation represent the basic run. For model verification, the results are compared with other available simulations of Henry problem. Figure 2b shows a comparison between the results of current simulation by SEAWAT, other models, and the semi-analytical solution developed by [29]. It is clear that a good agreement between SEAWAT and other models has been achieved. As such, the developed model is considered calibrated and could be applied to study different settings of SWI problems considering different seaside and bed slopes toward the sea and the land.  The boundary conditions for head and concentration were set based on the following. At the land side, a hydrostatic freshwater head distribution was considered and a constant flux (Q in ) of 0.5702 m 3 /d/m was implemented. The concentration (C in ) was sent constant along the boundary and equals to zero mg/L. At the seaside, a hydrostatic seawater head was considered to represent the density variation with depth at this boundary, with a hydraulic head of 1.0 m above the upper boundary of the aquifer. A constant concentration of 35,000 mg/L (seawater concertation) was implemented. The upper and lower horizontal boundaries were set as no flow boundary for both the fluid and salt ions. As such, the stream lines should be parallel and the equipotential lines should be perpendicular to the upper and lower boundaries. The implement hydraulic parameters are presented in Table 1.

Model Verification
In this study, SEAWAT is employed to simulate Henry problem considering the above boundary conditions and hydraulic parameters. As shown in Figure 2a, the intrusion of 0.5 Isochlor (17,500 mg/L) reached a distance of 62.5 cm measured along the bottom boundary from the seaside. The results of this simulation represent the basic run. For model verification, the results are compared with other available simulations of Henry problem. Figure 2b shows a comparison between the results of current simulation by SEAWAT, other models, and the semi-analytical solution developed by [29]. It is clear that a good agreement between SEAWAT and other models has been achieved. As such, the developed model is considered calibrated and could be applied to study different settings of SWI problems considering different seaside and bed slopes toward the sea and the land.

Results
SEAWAT is applied on the Henry problem to study the effect of changing bed and seaside slopes on SWI considering a constant saltwater head of 100 cm at the seaside and a constant flux rate of (Qin) 0.5702 m 3 /d/m at the land side. The seawater level is expected to increase due to climate change, and groundwater levels in coastal aquifers are expected to decrease due to the expected increase of groundwater abstractions to substitute for the reduction of rainfall and natural recharge [13,14]. Three cases were considered to study the effect changing groundwater and seawater levels on saltwater intrusion including: (a) increasing seawater head, (b) decreasing groundwater level at the land side, and (c) a combination of both. Different settings of bed and seaside slopes were also examined, and all simulations were conducted under the steady state conditions to represent the ultimate case of inland migration of the seawater.

Effect of Changing Seaside Slope on SWI
A number of scenarios were considered to analyze the SWI problem under different slopes of seabed and aquifer bottom. In the basic run, the numerical model was employed to solve Henry problem using the same parameters and boundary conditions that were considered by other

Results
SEAWAT is applied on the Henry problem to study the effect of changing bed and seaside slopes on SWI considering a constant saltwater head of 100 cm at the seaside and a constant flux rate of (Q in ) 0.5702 m 3 /d/m at the land side. The seawater level is expected to increase due to climate change, and groundwater levels in coastal aquifers are expected to decrease due to the expected increase of groundwater abstractions to substitute for the reduction of rainfall and natural recharge [13,14]. Three cases were considered to study the effect changing groundwater and seawater levels on saltwater intrusion including: (a) increasing seawater head, (b) decreasing groundwater level at the land side, and (c) a combination of both. Different settings of bed and seaside slopes were also examined, and all simulations were conducted under the steady state conditions to represent the ultimate case of inland migration of the seawater.

Effect of Changing Seaside Slope on SWI
A number of scenarios were considered to analyze the SWI problem under different slopes of seabed and aquifer bottom. In the basic run, the numerical model was employed to solve Henry problem using the same parameters and boundary conditions that were considered by other researchers. In this study, the intrusion length is represented by the 0.5 equi-concentration (isochlor) line and is measured along the bottom boundary from the seaside. Based on the results of the basic run, the intrusion length was 62.5 cm measured along the bottom boundary, as shown in Figure 2a. In the first and second scenarios, two seaside slopes (1h:1v and 2h:1v) were considered as shown in Figure 3, and the same flux of the Henry problem was used. All other hydraulic and transport parameters remained constant as given in Table 1. The intrusion length of 0.5 equi-concentration line for the two cases reached 95.75 and 136.0 cm, respectively, measured along the bed as shown in Figure 4a,b. The relation between intrusion length and seaside slope was developed to determine the intrusion length (X T ) for different seaside slopes as shown in Figure 5. The results indicate that the intrusion length increases with the increase of the horizontal to vertical ratio of the seaside slope. researchers. In this study, the intrusion length is represented by the 0.5 equi-concentration (isochlor) line and is measured along the bottom boundary from the seaside. Based on the results of the basic run, the intrusion length was 62.5 cm measured along the bottom boundary, as shown in Figure 2a.
In the first and second scenarios, two seaside slopes (1h:1v and 2h:1v) were considered as shown in Figure 3, and the same flux of the Henry problem was used. All other hydraulic and transport parameters remained constant as given in Table 1. The intrusion length of 0.5 equi-concentration line for the two cases reached 95.75 and 136.0 cm, respectively, measured along the bed as shown in Figure  4a,b. The relation between intrusion length and seaside slope was developed to determine the intrusion length (XT) for different seaside slopes as shown in Figure 5. The results indicate that the intrusion length increases with the increase of the horizontal to vertical ratio of the seaside slope.

Effect of Changing Aquifer Bed Slope on SWI
In scenarios three and four, the effect of changing the aquifer bed slope on SWI is investigated. Two bed slopes were considered. In scenario three, the aquifer bed slope was taken as 20h:1v researchers. In this study, the intrusion length is represented by the 0.5 equi-concentration (isochlor) line and is measured along the bottom boundary from the seaside. Based on the results of the basic run, the intrusion length was 62.5 cm measured along the bottom boundary, as shown in Figure 2a.
In the first and second scenarios, two seaside slopes (1h:1v and 2h:1v) were considered as shown in Figure 3, and the same flux of the Henry problem was used. All other hydraulic and transport parameters remained constant as given in Table 1. The intrusion length of 0.5 equi-concentration line for the two cases reached 95.75 and 136.0 cm, respectively, measured along the bed as shown in Figure  4a,b. The relation between intrusion length and seaside slope was developed to determine the intrusion length (XT) for different seaside slopes as shown in Figure 5. The results indicate that the intrusion length increases with the increase of the horizontal to vertical ratio of the seaside slope.

Effect of Changing Aquifer Bed Slope on SWI
In scenarios three and four, the effect of changing the aquifer bed slope on SWI is investigated. Two bed slopes were considered. In scenario three, the aquifer bed slope was taken as 20h:1v representing a mild increase in the aquifer depth toward the seaside. In scenario four, the bed slope was taken as 10h:1v ( Figure 6). All other parameters and boundary conditions remained unchanged, scenario three, Isochlor line 0.5 reached a distance of 61.25 cm measured from the seaside boundary, Figure 7a. Under scenario four, the same Isochlor line reached a distance of 59.5 cm, measured from the seaside boundary, Figure 7b. The relation between intrusion length and bed slope toward the land is presented in Figure 8. A linear equation is developed to determine the intrusion length (XT) for different bed slops (Figure 8). The results reveled that for the same hydraulic and boundary conditions, the intrusion length would decrease if the inland depth of the aquifer decreases with the distance from the shore boundary.

Effect of Changing Aquifer Bed Slope on SWI
In scenarios three and four, the effect of changing the aquifer bed slope on SWI is investigated. Two bed slopes were considered. In scenario three, the aquifer bed slope was taken as 20h:1v representing a mild increase in the aquifer depth toward the seaside. In scenario four, the bed slope was taken as 10h:1v ( Figure 6). All other parameters and boundary conditions remained unchanged, and the simulations were allowed to continue until reaching the steady state conditions. Under scenario three, Isochlor line 0.5 reached a distance of 61.25 cm measured from the seaside boundary, Figure 7a. Under scenario four, the same Isochlor line reached a distance of 59.5 cm, measured from the seaside boundary, Figure 7b. The relation between intrusion length and bed slope toward the land is presented in Figure 8. A linear equation is developed to determine the intrusion length (X T ) for different bed slops (Figure 8). The results reveled that for the same hydraulic and boundary conditions, the intrusion length would decrease if the inland depth of the aquifer decreases with the distance from the shore boundary.
Hydrology 2019, 6, x 6 of 14 and the simulations were allowed to continue until reaching the steady state conditions. Under scenario three, Isochlor line 0.5 reached a distance of 61.25 cm measured from the seaside boundary, Figure 7a. Under scenario four, the same Isochlor line reached a distance of 59.5 cm, measured from the seaside boundary, Figure 7b. The relation between intrusion length and bed slope toward the land is presented in Figure 8. A linear equation is developed to determine the intrusion length (XT) for different bed slops (Figure 8). The results reveled that for the same hydraulic and boundary conditions, the intrusion length would decrease if the inland depth of the aquifer decreases with the distance from the shore boundary.   Figure 7a. Under scenario four, the same Isochlor line reached a distance of 59.5 cm, measured from the seaside boundary, Figure 7b. The relation between intrusion length and bed slope toward the land is presented in Figure 8. A linear equation is developed to determine the intrusion length (XT) for different bed slops (Figure 8). The results reveled that for the same hydraulic and boundary conditions, the intrusion length would decrease if the inland depth of the aquifer decreases with the distance from the shore boundary.    In scenarios five and six, the bed slope was set as 20h:1v and 10h:1v, respectively ( Figure 9). The seawater represented by 0.5 Isochlor reached 64.75 and 66.5 cm, for scenarios five and six, respectively, as given in Figure 10a,b. A linear relation between the intrusion length and bed slope is established as illustrated in Figure 11. The results indicated that the intrusion length increases with the increase of the landward increase of the aquifer depth.    In scenarios five and six, the bed slope was set as 20h:1v and 10h:1v, respectively ( Figure 9). The seawater represented by 0.5 Isochlor reached 64.75 and 66.5 cm, for scenarios five and six, respectively, as given in Figure 10a,b. A linear relation between the intrusion length and bed slope is established as illustrated in Figure 11. The results indicated that the intrusion length increases with the increase of the landward increase of the aquifer depth.  In scenarios five and six, the bed slope was set as 20h:1v and 10h:1v, respectively ( Figure 9). The seawater represented by 0.5 Isochlor reached 64.75 and 66.5 cm, for scenarios five and six, respectively, as given in Figure 10a,b. A linear relation between the intrusion length and bed slope is established as illustrated in Figure 11. The results indicated that the intrusion length increases with the increase of the landward increase of the aquifer depth.   In scenarios five and six, the bed slope was set as 20h:1v and 10h:1v, respectively ( Figure 9). The seawater represented by 0.5 Isochlor reached 64.75 and 66.5 cm, for scenarios five and six, respectively, as given in Figure 10a,b. A linear relation between the intrusion length and bed slope is established as illustrated in Figure 11. The results indicated that the intrusion length increases with the increase of the landward increase of the aquifer depth.

Effect of Changing Sea Level on SWI
One of the most probable and important application of SWI problems is the effect of seawater level rise. Prior to 1990, available records from tide gauges showed a 1 mm sea level rise per year over a period of 2 centuries. After 1990, both satellites and tide gauges have revealed a rise of 3.2 mm per year [30]. Due to global warming trends and the associated melting of ice sheets and expansion of water in oceans, it is estimated that the seawater level will continue to rise even if the emission of greenhouse gases is controlled [31].
The possible seawater rise was simulated by increasing water level at the seaside by 5% and 10%, respectively. In scenario seven, the seawater level was set at 105 cm, while in scenario eight, it was set at 110 cm above the aquifer bottom. As before, a constant recharge 0.5702 m 3 /d of freshwater was considered at the land side. Under the steady state conditions, the 0.5 Isochlor reached a distance of 70.5 and 76.0 cm, respectively, for the two scenarios as shown in Figure 12a,b measured along the horizontal bed from the shore line. The relations between seawater level and intrusion length (XT) for different bed and seaside slopes are shown in Figure 13. A summary of the intrusion length for different cases is provided in Table 2.

Effect of Changing Sea Level on SWI
One of the most probable and important application of SWI problems is the effect of seawater level rise. Prior to 1990, available records from tide gauges showed a 1 mm sea level rise per year over a period of 2 centuries. After 1990, both satellites and tide gauges have revealed a rise of 3.2 mm per year [30]. Due to global warming trends and the associated melting of ice sheets and expansion of water in oceans, it is estimated that the seawater level will continue to rise even if the emission of greenhouse gases is controlled [31].
The possible seawater rise was simulated by increasing water level at the seaside by 5% and 10%, respectively. In scenario seven, the seawater level was set at 105 cm, while in scenario eight, it was set at 110 cm above the aquifer bottom. As before, a constant recharge 0.5702 m 3 /d of freshwater was considered at the land side. Under the steady state conditions, the 0.5 Isochlor reached a distance of 70.5 and 76.0 cm, respectively, for the two scenarios as shown in Figure 12a,b measured along the horizontal bed from the shore line. The relations between seawater level and intrusion length (XT) for different bed and seaside slopes are shown in Figure 13. A summary of the intrusion length for different cases is provided in Table 2.

Effect of Changing Sea Level on SWI
One of the most probable and important application of SWI problems is the effect of seawater level rise. Prior to 1990, available records from tide gauges showed a 1 mm sea level rise per year over a period of 2 centuries. After 1990, both satellites and tide gauges have revealed a rise of 3.2 mm per year [30]. Due to global warming trends and the associated melting of ice sheets and expansion of water in oceans, it is estimated that the seawater level will continue to rise even if the emission of greenhouse gases is controlled [31].
The possible seawater rise was simulated by increasing water level at the seaside by 5% and 10%, respectively. In scenario seven, the seawater level was set at 105 cm, while in scenario eight, it was set at 110 cm above the aquifer bottom. As before, a constant recharge 0.5702 m 3 /d of freshwater was considered at the land side. Under the steady state conditions, the 0.5 Isochlor reached a distance of 70.5 and 76.0 cm, respectively, for the two scenarios as shown in Figure 12a,b measured along the horizontal bed from the shore line. The relations between seawater level and intrusion length (X T ) for different bed and seaside slopes are shown in Figure 13. A summary of the intrusion length for different cases is provided in Table 2.        Increasing seawater level by 5% and 10% combined with different seaside and bed slopes affected the intrusion length as shown in Table 2 and Figure 13. For a seaside slope of 1:1, the intrusion of 0.5 Isochlor reached to 106.0 and 115.25 cm, respectively, while it reached 149 and 158.25 cm for a seaside slope of 2:1, respectively. For a bed slope of 1:20 toward the land, the intrusion reached to 66.6 and 72.25 cm, respectively, while it reached 65.5 and 71.0 cm, respectively, for a bed slope of 1:10. For bed slope of 1:20 toward the sea, the intrusion reached to 71.50 and 77.0 cm, respectively, while it reached 72.25 and 78.75 cm, respectively for a bed slope of 1:10 toward the sea. The bed slope toward the land gave less intrusion than the slope toward the sea for the two seawater levels. However, the seaside slope resulted in a bigger intrusion length for the two cases.

Effect of Reducing the Landside Recharge
In this case, the effect of excessive pumping from the aquifer and hence reduction of the specific recharge at the landside were considered. The case was simulated by decreasing the flux of freshwater by 5% and 10% from the original rate (Q in ) of 0.5702 m 3 /d to be 0.5417 m 3 /d and 0.5132 m 3 /d with a constant saltwater head, at the seaside of 1m. Figure 14a and b indicate that the intrusion of 0.5 Isochlor reached 66.70 and 68.15 cm, measured from the shore line along the bottom boundary, respectively.
The relationships between the recharge rate and intrusion length (X T ) for different slopes are shown in Figure 15. Decreasing the freshwater landside recharge rate by 5% and 10% combined with different seaside and bed slopes has affected the intrusion length with different values as shown in Figure 15. Increasing seawater level by 5% and 10% combined with different seaside and bed slopes affected the intrusion length as shown in Table 2 and Figure 13. For a seaside slope of 1:1, the intrusion of 0.5 Isochlor reached to 106.0 and 115.25 cm, respectively, while it reached 149 and 158.25 cm for a seaside slope of 2:1, respectively. For a bed slope of 1:20 toward the land, the intrusion reached to 66.6 and 72.25 cm, respectively, while it reached 65.5 and 71.0 cm, respectively, for a bed slope of 1:10. For bed slope of 1:20 toward the sea, the intrusion reached to 71.50 and 77.0 cm, respectively, while it reached 72.25 and 78.75 cm, respectively for a bed slope of 1:10 toward the sea. The bed slope toward the land gave less intrusion than the slope toward the sea for the two seawater levels. However, the seaside slope resulted in a bigger intrusion length for the two cases.

Effect of Reducing the Landside Recharge
In this case, the effect of excessive pumping from the aquifer and hence reduction of the specific recharge at the landside were considered. The case was simulated by decreasing the flux of freshwater by 5% and 10% from the original rate (Qin) of 0.5702 m 3 /d to be 0.5417 m 3 /d and 0.5132 m 3 /d with a constant saltwater head, at the seaside of 1m. Figure 14a and b indicate that the intrusion of 0.5 Isochlor reached 66.70 and 68.15 cm, measured from the shore line along the bottom boundary, respectively.
The relationships between the recharge rate and intrusion length (XT) for different slopes are shown in Figure 15. Decreasing the freshwater landside recharge rate by 5% and 10% combined with different seaside and bed slopes has affected the intrusion length with different values as shown in Figure 15.

Combined Effect of Increasing Sea Level and Decreasing Landside Recharge
This case represents the combination of the two previous scenarios; increasing seawater head by 5% and 10% and decreasing the flux of freshwater at the land side by 5% and 10%. The resulted intrusion of sweater is shown in Figure 16a,b. The intrusion of 0.5 Isochlor reached 72.9 and 81.5 cm, respectively. More inland intrusion has been observed in this scenario. Figure 17 shows the relation between the combined scenario and intrusion length (XT) for different slopes. Increasing seawater head and decreasing recharge rate by 5% and 10% combined with different seaside and bed slopes has increased the intrusion length with different levels as shown in Figure 17. For a seaside slope of 1:1, the intrusion length advanced to 108.5 and 120.5 cm, and for a seaside slope of 2:1, the corresponding intrusion lengths were, 151.0 and 162.25 cm, respectively.  Figure 17 indicates that the presence of seaside slope would increase the total intrusion length, while the presence of a bed slope toward the land would reduce the intrusion length for the studied cases.

Combined Effect of Increasing Sea Level and Decreasing Landside Recharge
This case represents the combination of the two previous scenarios; increasing seawater head by 5% and 10% and decreasing the flux of freshwater at the land side by 5% and 10%. The resulted intrusion of sweater is shown in Figure 16a,b. The intrusion of 0.5 Isochlor reached 72.9 and 81.5 cm, respectively. More inland intrusion has been observed in this scenario. Figure 17 shows the relation between the combined scenario and intrusion length (X T ) for different slopes. Increasing seawater head and decreasing recharge rate by 5% and 10% combined with different seaside and bed slopes has increased the intrusion length with different levels as shown in Figure 17.  Figure 17 indicates that the presence of seaside slope would increase the total intrusion length, while the presence of a bed slope toward the land would reduce the intrusion length for the studied cases.

Combined Effect of Increasing Sea Level and Decreasing Landside Recharge
This case represents the combination of the two previous scenarios; increasing seawater head by 5% and 10% and decreasing the flux of freshwater at the land side by 5% and 10%. The resulted intrusion of sweater is shown in Figure 16a,b. The intrusion of 0.5 Isochlor reached 72.9 and 81.5 cm, respectively. More inland intrusion has been observed in this scenario. Figure 17 shows the relation between the combined scenario and intrusion length (XT) for different slopes. Increasing seawater head and decreasing recharge rate by 5% and 10% combined with different seaside and bed slopes has increased the intrusion length with different levels as shown in Figure 17. For a seaside slope of 1:1, the intrusion length advanced to 108.5 and 120.5 cm, and for a seaside slope of 2:1, the corresponding intrusion lengths were, 151.0 and 162.25 cm, respectively. For a bed slope of 1:20 toward the land, the intrusion lengths reached to 69.0 and 77.0 cm, and for a bed slope of 1:10, the corresponding intrusion lengths were 68.0 and 75.25 cm, respectively. For a bed slope 1:20 of toward the sea, the intrusion lengths reached 74.0 and 83.0 cm, and for a bed slope of 1:10, the corresponding intrusion distances were 75.0 and 85.5 cm, respectively. Figure 17 indicates that the presence of seaside slope would increase the total intrusion length, while the presence of a bed slope toward the land would reduce the intrusion length for the studied cases.

Conclusions
Groundwater quality in coastal aquifers is vulnerable to deterioration due to the seawater intrusion problem. The seawater level is expected to rise due to the anticipated climate change, and the freshwater flux from the landside is expected to decline due to limited natural recharge and overabstraction of groundwater. This will help accelerate the seawater intrusion problem. In this paper, the effects of aquifer geometric slopes, including seaside slope and bed slopes toward the land and the sea, on seawater intrusion were presented and discussed. Moreover, the effects of changing groundwater and seawater levels on seawater intrusion were investigated. The SEAWAT code was applied to Henry problem considering different seaside slopes of 1:1 and 2:1 and bed slopes toward the sea and the land of 1:20 and 1:10. Different conditions of seawater levels and land side recharge were considered with different boundary and geometric slopes. The results revealed that increasing the bed slope toward the sea would increase the intrusion length, while increasing the bed slope toward the land would decrease the intrusion of seawater. Increasing the seaside slope would increase the total intrusion length by 50% for a slope of 1:1 and 100% for a slope of 2:1. Increasing seawater level and decreasing the land side freshwater flux would also increase the intrusion length. The concurrence of these two scenarios gave additional inland intrusion. The study demonstrated the sensitivity of the seawater intrusion problem to a number of parameters including slopes of the bed and sea boundary, freshwater flux at the land side and seawater level.

Conclusions
Groundwater quality in coastal aquifers is vulnerable to deterioration due to the seawater intrusion problem. The seawater level is expected to rise due to the anticipated climate change, and the freshwater flux from the landside is expected to decline due to limited natural recharge and overabstraction of groundwater. This will help accelerate the seawater intrusion problem. In this paper, the effects of aquifer geometric slopes, including seaside slope and bed slopes toward the land and the sea, on seawater intrusion were presented and discussed. Moreover, the effects of changing groundwater and seawater levels on seawater intrusion were investigated. The SEAWAT code was applied to Henry problem considering different seaside slopes of 1:1 and 2:1 and bed slopes toward the sea and the land of 1:20 and 1:10. Different conditions of seawater levels and land side recharge were considered with different boundary and geometric slopes. The results revealed that increasing the bed slope toward the sea would increase the intrusion length, while increasing the bed slope toward the land would decrease the intrusion of seawater. Increasing the seaside slope would increase the total intrusion length by 50% for a slope of 1:1 and 100% for a slope of 2:1. Increasing seawater level and decreasing the land side freshwater flux would also increase the intrusion length. The concurrence of these two scenarios gave additional inland intrusion. The study demonstrated the sensitivity of the seawater intrusion problem to a number of parameters including slopes of the bed and sea boundary, freshwater flux at the land side and seawater level.

Conclusions
Groundwater quality in coastal aquifers is vulnerable to deterioration due to the seawater intrusion problem. The seawater level is expected to rise due to the anticipated climate change, and the freshwater flux from the landside is expected to decline due to limited natural recharge and over-abstraction of groundwater. This will help accelerate the seawater intrusion problem. In this paper, the effects of aquifer geometric slopes, including seaside slope and bed slopes toward the land and the sea, on seawater intrusion were presented and discussed. Moreover, the effects of changing groundwater and seawater levels on seawater intrusion were investigated. The SEAWAT code was applied to Henry problem considering different seaside slopes of 1:1 and 2:1 and bed slopes toward the sea and the land of 1:20 and 1:10. Different conditions of seawater levels and land side recharge were considered with different boundary and geometric slopes. The results revealed that increasing the bed slope toward the sea would increase the intrusion length, while increasing the bed slope toward the land would decrease the intrusion of seawater. Increasing the seaside slope would increase the total intrusion length by 50% for a slope of 1:1 and 100% for a slope of 2:1. Increasing seawater level and decreasing the land side freshwater flux would also increase the intrusion length. The concurrence of these two scenarios gave additional inland intrusion. The study demonstrated the sensitivity of the seawater intrusion problem to a number of parameters including slopes of the bed and sea boundary, freshwater flux at the land side and seawater level. Funding: The support of the National Water Center, UAE University towards the open access publication of this paper is acknowledged.