Estimation of Spatial Groundwater Recharge and Surface Runoff in the Gaza Coastal Aquifer Using GIS-Based WetSpass Model

Estimation of Spatial Groundwater Recharge and Surface Runoff in the Gaza Coastal Aquifer Using GIS-Based WetSpass Model

Aish, A. M.,

Earth Sciences Department, Faculty of Science, Al Azhar University-Gaza, Palestine

E-mail : aaish@usa.com

Abstract

Assessment of the spatial distribution of water recharge and surface runoff is essential for water resource management and planning. This research evaluates spatial distributions of groundwater recharge, surface water flow, and evapo-transpiration for the Gaza coastal aquifer, Palestine in 2020 using the WetSpass model integrated with GIS. The ArcGIS tool is implemented to create grid maps of the WetSpass model data, such as temperature, rainfall, wind velocity, potential evapo-transpiration, soil map, land use, topography, gradient, and water depth. The model generates computerized maps of yearly groundwater recharge, surface water flow, and evapo-transpiration. The WetSpass model's calibration procedures were implemented by manually altering or changing the model parameters already included in the WetSpass model within a predetermined range of values. The simulated results revealed that the average yearly groundwater recharge was 42 mm/year, surface runoff was 38 mm/year and actual evapotranspiration was 268 mm/year. The research findings would help local stakeholders and decision-makers manage groundwater resources in a sustainable and effective manner.

Keywords: Water balance, GIS, Spatial Analysis, WetSpass

1. Introduction

The Gaza Strip (Palestine) is located on the southern end of the Mediterranean Sea, between latitudes 31° 16 and 31° 45 north and longitudes 34° 2 and 34° 25 east. Its coastal zone is 42 km long and ranges in width between 6 and 12 km, encompassing an area of 365 km ²[1]. Figure 1 illustrates a location map of the Gaza Strip. Elongated ridges and depressions characterize the topography. It extends NE–SW direction [2]. Ridges are composed of sandstone locally named as Kurkar [3]. The Kurkar deposits include silt, clay, and conglomerates [4]. The coastal aquifers are separated into four sub-aquifers close to the coastline. As depicted in Figure 2, the top sub-aquifer "A" is unconfined, while sub-aquifers "B1, B2, and C" are semi confined. Gaza coastal aquifer’s average thickness is about 120 meters at the coastline. The thickness along the northern Gaza border is about 5 to 10 meters and in the south is around 60 meters [5]. Rainfall naturally recharges the Gaza coastal aquifer, and irrigation return flow provides additional recharge. The water level is around 60 meters below the ground in the east, whereas in the west, it is just a few meters deep along the coast [6].

The soil is mainly composed of sands, clay, and loess. The sand is founded in the form of dunes covering the coastline from the north to the southern border of the Gaza Sector. Due to the dunes' shape, the thickness ranges from two to about 50 meters. Clay is found in the north-eastern of Gaza. Wadis are surrounded by loess soil of thickness ranging from 25 to 30 meters [7]. Temperature is affected by geographical conditions, altitude, marine exposure, etc. Over the year, the monthly mean temperature typically varies in the Gaza Strip as the highest temperature is recorded in August (summer) while the lowest temperature is recorded in January (winter) [9]. For the year 2020, the temperature average varied from 15.7 degrees Celsius in January and 27.9 degrees Celsius in August. The wind velocity peaks at noon and gradually decreases throughout the night. Almost all the winds originate in the southwest during the winter, and the strong northwest wind blows in summer [10]. The average monthly minimum wind velocity is 3.2 m/s in October and the average monthly maximum wind velocity is 5.6 m/s in February.

Figure 1: Map of the Gaza Strip's location and topography (Data source:PWA)

Figure 2: Hydrogeological cross section of the Gaza coastal aquifer [8]

The average monthly minimum humidity is 65.5 % for November and the average monthly maximum humidity is 85.6 % for July. The average monthly minimum sunshine hours is 4.4 for December and the average monthly maximum sunshine hours is 11.1 for July. The ArcGIS World Geocoding Service website was used to collect climatological information such as humidity, temperature, wind velocity, solar radiation, and sunshine hours, as illustrated in Table 1. Rainfall is considered the primary water source in the Gaza Sector; the most rainfall is observed in the winter, especially between October and March, with a significant annual variation. The rainfall increases from the south to the north [11]. The Palestinian Water Authority (PWA) and Ministry of Agriculture (MoA) gather daily data from 12 measuring rain gauge stations. Figure 3 shows the areal rainfall intensity distribution for the year 2020. The values vary from 579 millimeters annually in the northern area to 225 millimeters annually in the southern area, with an annual average is 348 millimeters.

2. Materials and Methods

2.1 Hydrological Modeling using WetSpass Model

Estimating groundwater recharge and surface runoff is essential for managing both surface and groundwater resources. Calculating groundwater recharge is difficult and complex because it depends on a variety of elements, including terrain, soil texture, land use and cover, groundwater depth, meteorological conditions, and other hydrologic variables. The WetSpass model was created as a physically-based approach to calculate the actual average evapotranspiration, groundwater recharge, and surface runoff [12] and [13]. A WetSpass Model is used in the current study to calculate the geographical distribution of water balance components. The model considers each raster cell's spatial distribution of topography, soil texture, slope, and climatic parameters. The WetSpass approach describes an area or basin as a regular grid cell arrangement. The water balance considerations of each locality are achieved by the vegetated, open-water, bare soil, while impervious fractions per grid pixel are calculated by the following equations [14].

ET _m = a _b E _b + a _v ET_v + a _o E _o + a _i E _i

Equation 1

S _m = a _b S _b + a _v S_v + a _o S _o + a _i S _i

Equation 2

R _m = a _b R _b + a _v R_v + a _o R _o + a _i R _i

Equation 3

Where ET _m is the total evapotranspiration (mm), S _m the surface runoff (mm), R _m is groundwater recharge (mm), each having (b) bare soil, (v) vegetated, (o) open water and (i) impermeable region component. The terms a _b, a _v, a _o and a _i are the fraction region of bare soil, vegetated, open water and impervious region, respectively [15]. A detailed description of the WetSpass Model can be found in the manual available on the website of the Vrije University of Brussels, Belgium (VUB).

2.2 WetSpass Model Data Input

The model requires two different kinds of input data: GIS grid maps and parameter tables. GIS grid maps comprise climatological data (rainfall, potential evapotranspiration, average temperature, and wind speed), groundwater depth, slope, topography, soil type, and land use. The input data were set up as a raster map using the ESRI ASCII grid format with a cell size of 50 m × 50 m. The input parameters for the WetSpass Model are shown in Table 2.

3. Results and Discussion

The simulation results are organized by digital images of the spatial distribution of yearly mean values of groundwater recharge, surface runoff, and actual evapotranspiration for the year 2020. These maps are raster-shaped, and each pixel displays the size of the corresponding water balance component, represented as layer thickness in millimeters.

3.1 Groundwater Recharge

Groundwater recharge is important for evaluating groundwater resources, although it is difficult to estimate [16]. The WetSpass model determines the annual average areal recharge distributed as an area variable based on the terrain, slope, soil, land use, climate, etc. this is mainly to consider how the groundwater system is impacted by the geographical variability of the land surface [17]. The simulated annual mean of groundwater recharge of varied from 0 millimeters per year (minimum value) to 192 millimeters per year (maximum value), with a mean of 42 millimeters per year and a standard deviation of 40 millimeters per year as illustrated in Figure 4. Precipitation and runoff have an impact on the variability of groundwater recharge. It indicates that there is significant groundwater recharge in the northern governorate. It appears that there is little groundwater recharge in the Rafah governorate. Higher groundwater recharge values were simulated for areas with lower topography and porous soil, like the coastal regions in the north. Infiltration into the groundwater is further influenced by plant cover, slope, soil type, the existence or lack of clay lenses, and water table depth.

3.2 Surface Water Runoff

The WetSpass model uses the runoff coefficient technique to calculate surface runoff. The coefficient depends on the vegetation cover type, slope gradient, and soil composition [18]. Surface runoff in the Gaza coastal aquifer varies geographically with topography and other catchment characteristics. The simulated annual surface runoff varied from 0 millimeters per year (minimum value) to 240 millimeters per year (maximum value) with a mean of 38 millimeters per year and a standard deviation of 48 millimeters per year as illustrated in Figure 5. Surface runoff is relatively low in most of the Gaza Strip locations. All values are below 92 millimeters per year, except for urban areas like the city of Gaza and the camps of refugees, which have the highest runoff values.

3.3 Actual Evapotranspiration

The model estimates the overall real evapotranspiration as the result of the transpiration of the vegetative cover, the evaporation of water, and the evaporation from bare soil between plants [19].

Figure 4: Simulated groundwater recharge with WetSpass model

Figure 5: Simulated surface runoff with WetSpass model

Figure 6: Simulated actual evapotranspiration with WetSpass model

Table 3: Comparison of the Gaza coastal aquifer's water balance component

The simulated annual evapotranspiration of Gaza’s coastal aquifer ranges from 130 millimeters per year (minimum value) to 391 millimeters per year (maximum value) as the minimum with an average value of 268 millimeters per year and a standard deviation of 46 millimeters per year. The simulated annual transpiration extended from 0 millimeters per year (minimum value) to 172 millimeters per year (maximum value) with an average value of 64 millimeters per year. The annual interception varied from 0 millimeters per year (minimum value) to 22 millimeters per year (maximum value) with an average value of 8 millimeters per year. The soil evaporation yearly average varied from 0 millimeters per year (minimum value) to 249 millimeters per year (maximum value) with an average value of 196 millimeters per year as illustrated in Figure 7. The average of evapotranspiration is about 77 % of the annual average of rainfall. A variation in evapotranspiration that correlates with rainfall and vegetation cover can be noticed in the results. The city of Gaza and camps of refugees are primarily built-up areas with limited vegetation cover.

3.4 Comparison of Water Balance for the Gaza Coastal Aquifer with Earlier Studies

Water balance mainly describes the net impact of water inflow and outflow. The most crucial input component is precipitation. Surface runoff, evapotranspiration, and groundwater discharge are water balances' three most significant outflow elements.

Figure 7: Simulated transpiration, interception and soil evaporation with WetSpass model