Integration of GIS and F-Hydra Model for Aquifer Vulnerability Monitoring in the Afram Plains, Ghana

Aquifer vulnerability is a critical issue across the entire globe due to the depletion of groundwater and escalation of pollution levels, which poses a detrimental effect on the natural environment. To ascertain this contamination risk, an extensive study has been conducted to assess the aquifer vulnerability by using the F-hydra model. This paper presents the vulnerability technique for the theory and application of flow accumulation, land-use and hydraulic conductivity. The model was applied to a shallow aquifer in eastern Ghana’s Volta River Basin, with the results being compared to the standardised DRASTIC model. The model follows the aquifer vulnerability assessment concept of the source pathway receptor where flow accumulation represents the ponding areas with downward percolation of contaminant to the water table, land-use represents the human activities at the land surface, and hydraulic conductivity represents the driving force leading to the movement of contaminant. The results reveal that the moderate vulnerability region covers 51.55% (2598.12 km2) of the entire area. The high and low vulnerability regions cover a significant percentage of the area 1.13% (56.52 km2) and 47.32% (2384.93 km2), respectively. The final vulnerability index was largely influenced by the removal of the hydraulic conductivity and land-use parameters. The model was validated with nitrate concentration in drilled wells in the study area. The maps produced in this study could be utilised as a guide to vulnerability by policymakers, groundwater manager and planners aimed at preserving the integrity of this vulnerable resource. *All authors contributed equally to this work. #Corresponding author. How to cite this paper: Agyare, A., Kabo-Bah, A. T., Bayel, B. T., & Jalloh, A. S. (2021). Integration of GIS and F-Hydra Model for Aquifer Vulnerability Monitoring in the Afram Plains, Ghana. Journal of Geoscience and Environment Protection, 9, 222-248. https://doi.org/10.4236/gep.2021.93014 Received: January 25, 2021 Accepted: March 27, 2021 Published: March 30, 2021 Copyright © 2021 by author(s) and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access


Introduction
All groundwater is vulnerable to pollution from surface activities. However, the degree of vulnerability is dependent on several factors. Of prime importance is the nature of the aquifer system, i.e., its ability to act as a potential water resource (Stournaras et al., 2011;Lv et al., 2021;Mutiara et al., 2021;Xiang et al., 2021). In every community, water plays an important role. In addition to the fact that it is essential for life, groundwater facilitates a wide variety of ecological systems and contributes to the growth of the economy and recreational activities (Anornu et al., 2012;Egbi et al., 2019;Chakraei et al., 2021;Meran et al., 2021;Shen, 2021;Wang et al., 2021). Groundwater has been regarded as an essential and reliable resource for a water supply because of its generally low vulnerability to contamination in contrast with surface water and its capacity for storage (US EPA, 1985;Akter & Ahmed, 2021;Lv et al., 2021). One classical example is the Nubian Sandstone water and its area coverage of about 2 million km 2 . This aquifer is a vital water resource for groundwater level and therefore sustainable development plans are needed to protect these vulnerable resources (Nageswara Rao & Narendra, 2006).
Groundwater has a major contribution to agricultural, industrial, drinking, and other municipal uses, particularly in the region. Moreover, if it is also to be a dependable resource, capable to meet current and impending demand, then it has to be protected and managed effectively, not merely the quantity but also the quality. Thus, in defining sustainable water resources management plans at the regional or continental scale, assessments of groundwater resources and their associated impacts are strongly needed (Hasiniaina et al., 2010).
A thorough understanding of the hydrogeological characteristics of aquifer systems concerning their environmental vulnerability is necessary to effectively exploit and protect groundwater against deterioration from its pristine state. Aquifers are not only characterised by hydraulic conductivity, but also by transmissivity (product of aquifer thickness and hydraulic conductivity) and diffusivity (ratio of storage coefficient and transmissivity). Others include soil composition or characteristics, existing climatic conditions, topography, mode and source of recharge, pH, the resident time of water within the geological strata, the permeability of the soil cover and drainage area (Vrba & Zaporozec, 1994;Olorunfemi et al., 1999;Ouedraogo et al., 2016;Shirazi et al., 2012;Shirazi et al., 2013).
The culmination of unregulated urbanization and industrialisation (as evident in most developing countries) compromise the pristine quality of several urban groundwater resources. While assessing the degree of vulnerability of aquifers and their protection from pollution, it is essential to comprehend the intrinsic property and hydrogeological characteristics. The aforementioned properties rely on the vulnerability of the aquifer system to natural or human impacts (Vrba & Zaporozec, 1994).
Groundwater protection requires extensive information on aquifer vulnerability, notably mapping intrinsic properties of aquifers to pollution. For certain scenarios, a comprehensive vulnerability evaluation of the actual quantitative and qualitative status of a particular aquifer system is not feasible. This is attributable to inadequate knowledge for aquifer monitoring or the complexity of the aquifer systems (Worrall & Kolpin, 2003). Alternatively, aquifer vulnerability indices are established and mapped to indicate or predict the actual severity of human-induced deterioration in groundwater quantity and quality (Balaji et al., 2021;Li et al., 2021;Rani et al., 2021).
In data-limited regions, the issue of inadequate monitoring data is further exacerbated. Data limited regions are major regions with little or no documented scientific information for research applications (Oke & Vermeulen, 2020). Major regions in the African continent lack extensive geohydrological data owing to a reduced government budget on data collection and management of information. Most models of aquifer vulnerability assessments are regarded to be inapplicable to many regions of the African continent at times. This is not due to the scientific basis of the model, but primarily because of data unavailability.
The Afram Plains in the south-eastern region of Ghana's Volta River Basin, one of the Lower Voltaian sedimentary basins of West Africa (see Figure 1) is Journal of Geoscience and Environment Protection affected by the challenges of limited comprehensive hydrogeological data. Groundwater occurrence in the basin is found in a confined and unconfined state, depending on the sedimentary rock depositions that host the aquifer. The aquifers in the Afram Plains have also been tagged risky by Davies (2002), which is attributed to the extensive application of fertilizers and pesticides from the agricultural farmlands occupying about 70% of the entire landmass. Besides, other sources may include contaminants from septic or pit latrines, fuel stations, household and industrial waste disposal mainly from the developed areas (Agyare et al., 2017). Assessing the aquifer vulnerability in the basin, specifically, the unconfined aquifers would require evaluating the contributing factors responsible for groundwater protection. A breakdown by Agyare et al. (2017), indicates that more than half of the populace use water from tube wells (boreholes) while the remaining access water from large diameter hand-dug wells.
As mentioned earlier, various models for assessing groundwater vulnerability have been proposed and can also be useful in evaluating the aquifer systems in Afram Plains. However, due to the heterogeneity and complex nature of the aquifers and their protective cover coupled with insufficient hydrogeological data, there is the need of proposing other simplified models with fewer data needs. Thus, one way to develop an easily applicable model is to reduce the number of parameters (Gogu et al., 2003). The simplified models will be targeting the intrinsic properties of the aquifer protective cover and depicted with vulnerability maps. Natural protection investigations above the aquifer are required to promote laws and policies regarding land-use practises preventing contamination of groundwater. The objective of this study is to evaluate the vulnerability of aquifer systems in the Afram Plains of Ghana and determine the relationship be-tween the aquifer vulnerability indices and on-ground nitrate concentrations.

Site Description
Afram Plains is located in the lower Volta River Basin in the Eastern Region of Ghana. It is sited between latitude 6.6080˚ and 7.1300˚ and longitude −0.0542˚ and −0.4926˚ (see Figure 1). The Afram Plains has a geographical area of 5446 km 2 , as the largest district in the region, with an undulating land surface ranging between 60 to 120 m above mean sea level. The area is bounded to the north by the Obosum River, to the east by the Volta Lake, and to the south-west and south by the Afram River (Agyare et al., 2017).
The average annual rainfall in the Afram Plains ranges between 1150 and 1650 mm. However, the highest rainfall ever recorded in the Afram Plains occurred in the year 1986 with an annual maximum depth of 1832 mm and a minimum depth of 546 mm, which also occurred in 1987. The wet season in the Afram Plains starts from April to October and the dry season from November to late March. Some water demands are met during the wettest rainy period from May to September. The climate within the year traverses between warm and hot with annual mean temperatures of about 36.6˚C to 36.8˚C which generally peaks up around February to late March before the start of the rainfall season. The relative humidity is generally between 70% to 80% throughout the year (Agyare et al., 2017).
The soil types found in the Afram Plains are plinthosols, fluvisols, regosols, lixisols, luvisols and leptosols. Lixisols are the most dominant soil type that is found in the area. The aquifer systems are characterized by semi-confined aquifers with low permeability and secondary permeability within fractured zones (Davies, 2002). The geologic formation in the north-eastern margin of the study area is mainly composed of massive coarse-grained conglomerates, sandstones and shales. The quartzitic sandstones cap the north-western and south-western margins. Also, the central part of the study area is underlain by cross-bedded yellowish-white, greyish-purple and red quartzitic medium-grained sandstones, fine-grained sandstones, siltstones and conglomerates. The south-eastern portion is underlain by cross-bedded dark greyish purple and reddish-brown medium-to-fine grained sandstones with interbedded thin gravel bands as illustrated in Figure 2.

Data Acquisition
A GIS spatial database was constructed for the hydraulic conductivity, land-use and aquifer recharge of the study area. All the data was processed using ArcGIS 10.6 TM . Besides, spatial data possess different spatial resolutions, hence, the data layers were resampled to fit the recommended resolution with a scale of 1:600,000. Thus, 0.5 km × 0.5 km resolution was selected for this study.
Data acquisition in Ghana is a great concern in Ghana as the required data is  not readily accessible by the public. Thus, the data were derived from both primary and secondary sources including field works, literature reviews and laboratory analyses for hydro-chemical parameters. The methodology deployed in the generation of the maps is summarised in Figure 3.

Sampling and Experimental Analysis
The hydrogeochemical characterisations of groundwater environments, in both time and space, are essential for understanding the variability in terms of the physical and chemical characteristics of groundwater. The chemical composition of groundwater is essential for assessing the groundwater quality for vulnerability purposes. Water rock contact, mode of recharge, and human impact are only a few of the factors that influence the chemical characteristics of groundwater (Chilton, 2006;Robins et al., 2007;van Beynen et al., 2012;Erostate et al., 2018;Jaunat et al., 2019).
Groundwater sampling was conducted from March 2020 to April 2020. This enabled groundwater samples to be collected in both dry and wet seasons, which is significant for groundwater elemental concentration monitoring. The collected samples include shallow unconfined borehole and hand-dug wells which are the prime targets for this study. The unfiltered water samples were collected in sterilised 100 mL plastic bottles per each site using a bailer and transported in a refrigerated box to the laboratory between 1˚C and 9˚C before analysis. The wells sampled for water quality analysis and each site visited in the study area are shown in Figure 4.

Hydro-Chemical Data Collection
Water samples were collected in 500 mL bottles, unpreserved, unfiltered, and stored below 9˚C for anion determination. Physico-chemical properties such as

The F-Hydra Model
The subjective vulnerability model that was developed to address this gap is the F-hydra model. This model is focused on the diffuse recharge through a specified soil thickness with the assumption of precipitation as the driving mechanism in the vertical infiltration and subsequent percolation that replenishes the groundwater. The F-hydra model employs the objective as well as the subjective based vulnerability models. The F-hydra model was developed by Kabo-bah et al. (2014), which considers 3 core parameters as critical in the determination of the extent of vulnerability of the aquifer to contamination. The main condition in the use of the F-hydra model is to understand the basic parameters, namely flow accumulation (fa), hydraulic conductivity (hyd), and landcover (lan).
The F-hydra model stands on the basis that the amount of precipitation infiltration and land-use activities at the surface must be known, and the dissolved solute contaminants in the infiltrating water have the same mobility rate as water. The contributing parameters are assigned a rating of 1 to 10 based on a characteristic range of values. The assigned weights vary from 30 to 40, whereas the highest weight denotes the most critical parameter. The total summation of the product of the ratings and weights of the respective parameters gives the vulnerability index. Table 1 and Table 2 presents the data required for the model and the assigned weights for each parameter, respectively. Equation (1) shows an addition formula of the F-hydra vulnerability index expressed as a percentage.  The concept of the model was originally defined by the European Groundwater Vulnerability Approach, origin-pathway-target (Daly et al., 2002;European Commission COST Action 620, 2003). In general, the fa parameter is used to define the drainage pattern of terrain and provides information about the potential availability of groundwater recharge from rainfall. The hyd parameter is a function of the recharge rate, aquifer thickness, and the length of effective groundwater drainage, while the lan parameter, on the other hand, is a function of the soil type, rock formation and slope (Kabo-bah et al., 2014). map from the computation of the flow direction will automatically be used to derive the flow accumulation map of the terrain. The output map contains cumulative hydrologic flow values that signify the number of input pixels that contribute any water to the outlets (Step 4). The flow accumulation rating is presented in Figure 6 and Table 3.

Hydraulic Conductivity Parameter
This parameter measures permeability. It depends on the intrinsic permeability    of the material and the degree of saturation. This is a critical factor that controls the contaminant migration and dispersion from the injection point within the unsaturated zone through to the saturated zone. Therefore, the higher conductivity, the more susceptible the aquifer. The calculated hydraulic conductivity for the study area was evaluated for consistency with other literature estimates. In this study, the hydraulic conductivity was inferred from the global permeability map. The data were given in log permeability (logk). Since hydraulic conductivity is a critical parameter in assessing aquifer vulnerability, the permeability (k) was converted into conductivity values using Equation (2).
However, it was revealed that literature estimates of 0.01 -0.05 m/day by Lutz et al. (2007). The estimates from Equation (2) vary between 0.01 -0.03 m/day.
For is well known, that estimating accurately the hydraulic conductivity is quite difficult, yet the overall consistency of this estimate was found to be satisfactory.
The most prevalent value covering about 34% of the area (R = 6), followed by 35% of the area (R = 4) and the remaining 31% (R = 2) as presented in Table 4 and

Landcover Parameter
The landcover parameter reveals the change or modification to the natural environment as a result of the activities on the Earth surface. These activities may cause a detrimental effect on both the quality and quantity of the aquifer. The following are the definition of the various landcover classes identified as shown in Figure 8.  Water bodies: This class represents the stream networks (Lake Volta and its tributaries), the reservoirs and any identifiable water bodies on the surface.  Evergreen Forest: It characterizes a group of trees holding their leaves permanently throughout the year. The decrease of such forest in the study area is also significant probably due to deforestation.  Shrublands: It characterizes vegetation dominated by shrubs including grasses, herbs, and geophytes.  Grasslands: Represents all grasses, shrubs and all types of small vegetation.  Croplands: These types of lands are mainly used for crop cultivation. The agricultural activities in the area are mainly rain-fed agriculture. Besides, some irrigated areas are identified and are all classified under croplands. Thus, one of the major non-point source pollutions of the aquifer is caused by nitrate in the agriculture areas. Hence, high contamination risk should be anticipated due to the intense agricultural activities and input such as fertilizer and pesticides. In this study, 15 series of land-use were identified and re-grouped into 5 classes namely: water bodies (R = 3), evergreen broadleaf forest (R = 1), shrublands (R = 4), grasslands (R = 4) and croplands (R = 10) as illustrated in Figure  8. Subsequently, the weight and ratings assigned to the land-use parameter are presented in Table 5.

Instance(s) of F-Hydra Model Application and Limitations
The F-hydra model has not been applied widely in aquifer vulnerability studies. Kabo-bah et al. (2014), utilized this newly developed F-hydra index to evaluate the pollution potential of the Densu Basin. The results indicated that more than 65% of the entire landmass is highly susceptible to pollution once they occur while less than 2% of the total areas were classified as low.
Notwithstanding, the F-hydra as compared to other models has certain practical constraints. Several drawbacks of the F-hydra model include:   2) Complexity in updating maps (Gogu & Dassargues, 2000); and; 3) Inaccurate in predicting the extent of pollution or vulnerability to some degree in the subsurface region (Kabo-bah et al., 2014).

Sensitivity Analysis
Sensitivity analysis presents a thorough detail on the influence of weights or ratings assigned to the individual parameters considered for the model (Krishna et al., 2015). There are 2 types of sensitivity analysis namely; single-parameter analysis and map removal sensitivity analysis. The single-parameter removal sensitivity analysis test points out the influence of each parameter on the final vulnerability index. Thus, the effective weight of each parameter can be estimated using Equation (3). The second category represents the sensitivity of the final vulnerability map to the removal of one or more parameter map layers and determined using Equation (4). Before the analysed results, authors can ascertain that their assigned weight was accurate or there is a need for modification. Estimation of the variation index demonstrates the influence of each map layer on the overall vulnerability index when removed.
where S i is the sensitivity index associated with the removal of one parameter map; V i and V xi are the unperturbed and perturbed vulnerability indices respectively; i signifies the cell number; N and n signify the number of parameters used to compute V i and V xi ; W i signifies the effective weight; P r,i and P w,i are the rating and the weight of each parameter map layer respectively; V i refers to the overall vulnerability index; V ar is the variation index of the removal parameter map layer. Thus, the unperturbed vulnerability index is considered to be the overall vulnerability index obtained using all the 3 parameter map layers and the perturbed vulnerability index refers to the index calculated using fewer parameter map layers.

The F-Hydra Vulnerability Index
The underlying principle of this study is that all aquifer systems are vulnerable to some extent. As a result, this present work aimed to evaluate the aquifer vulnerability to pollution in the Afram Plains District. To achieve this goal, the F-hydra model was applied for this evaluation and compared to the widely known DRASTIC model to ascertain the dependability of the results. The final F-hydra vulnerability map was produced by combining the 3 vital parameters. This was however after classification through numerical rating and creating a raster file from the numerical values. The F-hydra index was computed using Equation (1). The aquifer vulnerability index was sub-divided into 3 classes. These classes were: Low (140 -420), Moderate (421 -580), High (581 -840). The vulnerability zones of the study area were obtained with the indication of indexes as illustrated in Figure 9.
The vulnerability obtained is low to moderate vulnerability, and this is according to the indication of indices. The percentages of the area covered under the various classes are 47.32%, 51.55% and 1.13% for low, moderate and high, respectively are presented in Table 6. The vulnerability distribution of the study  area shows low vulnerability in areas with low recharge and low hydraulic conductivity. This is due to geologic formation hosting the aquifer in these areas are composed of tertiary clay and clayey-sand deposits which are impermeable (Al-Adamat et al., 2003;Neukum & Azzam, 2009;Al Hallaq & Elaish, 2012). The final F-hydra output map was compared with the DRASTIC output map of the study area developed by Agyare et al. (2017). The rationale behind this comparison is because DRASTIC model has been recognized as a standardised model that offers the user the ease to explain the vulnerability index from their expert knowledge based on the hydrogeological setting of the terrain (Gogu et al., 2003).
According to Figure 10, the DRASTIC index range between 71 and 197. The aquifer vulnerability derived was classified from low to high based on the indexes. The percentages of the area covered under the various classes are 20.71%, 49.78% and 29.50% for low, moderate and high, respectively, as presented in Table 7. The distribution of the aquifer vulnerability by the DRASTIC model exhibits high vulnerability in regions with low topography, depth-to-water, and relatively high hydraulic conductivity, while the low vulnerability was mapped for regions having a high vadose thickness, high clayey content in soils and likely the lowest precipitation region of the basin (Agyare et al., 2017).
Comparatively, both models yielded different results in terms of percentages of the mapped zones (see Table 7). These differences are quite significant with exception of the moderate vulnerability class which can be argued off as minor. However, the F-hydra model mapped the central and eastern margins as moderately vulnerable whilst the DRASTIC model assessed the region as highly vulnerable. The difference between the moderate and high vulnerability in the DRASTIC map is greatly influenced by the depth-to-water table.
Obviously, the depth-to-water table parameter in the DRASTIC model determines whether a zone will be classed as moderate or high vulnerability. Nonetheless, this parameter is not incorporated in the F-hydra model for aquifer vulnerability assessment (Gogu et al., 2003;Gumma & Pavelic, 2013). This could be the reason for the differences in the percentages of mapped zones, especially between the moderate and high vulnerability regions.

Sensitivity Parameter Analysis
A summary of the parameter relationship analysis for the F-hydra parameters is presented in Table 8. The result of the analysis reveals that the land-use   parameter has the highest rate value (mean = 5.15). Thus, a high contamination risk of the aquifer is anticipated to emanate from this parameter. Moreover, the flow accumulation parameter (mean = 5.12) and hydraulic conductivity (mean = 4.24) have a moderate risk of contamination. A high influence on the variation of the final vulnerability index is anticipated by the coefficient of variability of the hydraulic conductivity parameter (CV = 65.23%), while moderate influence is anticipated by the variability of land-use (CV = 51.89%) and flow accumulation (CV = 46.21%) parameters.

Map Removal Sensitivity Analysis
The summary of the map removal sensitivity calculated by removing one or more map layer(s) is presented in Table 9. The result indicates that the hydraulic conductivity and land-use parameters influence the final vulnerability index greatly according to their mean variation (V i = 4.08% and 2.44%). This is as a result of the high weight assigned to these parameters (W T = 40 and 30). Besides, it was observed that the influence on the final vulnerability index was moderate when the flow accumulation parameter was removed (V i = 1.83%).

Single Parameter Sensitivity Analysis
The results obtained demonstrated some variations from the theoretical weights as presented in Table 10. The result reveals that the effective weight of parameter flow accumulation (W eff = 30.00%) is greater than the theoretical weight (W T = 34.90%). Also, the flow accumulation parameter was considered the most effective parameter for this analysis. This result reaffirms the reason why this parameter is considered a critical parameter in this model. However, other parameters such as land-use and conductivity have their mean effective weight value (W eff = 29.51% and 15.09%) to be less than the theoretical weight (W T = 30.00% and 40.00%).

Vulnerability Map Validation
The relationship between various indexes of aquifer vulnerability and nitrate concentrations in wells are discussed in this section. The justification behind the selection of this chemical parameter was due to the improper use of fertilizers and manure from the agricultural sectors. This model, however, presumes that the contaminant has to migrate with water. Nitrates are generally very stable,  A. Agyare et al. dissolve quickly in water and can migrate several kilometres from their source fulfil this assumption. Consequently, production drilled wells from 30 separate locations were sampled and analysed for nitrate concentrations in a laboratory. The handheld Global Positioning System (GPS) Garmin etrex 20 was deployed to record the geographical locations of the sampled drilled wells.

The Spatial Nitrate Concentration
The laboratory assessment on the sampled wells revealed that the nitrate concentrations vary from 0.21 mg/L to 34.89 mg/L as presented in Figure 11. Table  11 presents a statistical summary of the sampling distribution. The acceptable threshold for nitrate concentration in drinking water is 50 mg/L according to the World Health Organization (WHO) guideline value. Also, many researchers have specifically limelight that nitrate concentrations greater than 10 mg/L in production wells suggest pollution from human sources and can lead to serious and fatal health consequences, such as methemoglobinemia (Sener et al., 2009;Adimalla et al., 2021;Emeka et al., 2021;Kühling et al., 2021;Tsagris & Tzouvelekas, 2021). Graphically, Figure 12 presents the nitrate concentrations of the sampled wells in the study area. A layer of the sampled wells was generated and overlayed on the intrinsic vulnerability map using ArcGIS. This was conducted to ascertain the number of wells with low, moderate and high nitrate concentrations falls within their corresponding vulnerability classes for validation purposes.

Regression Analysis
The model was evaluated using one of the most widely used statistical tests in geohydrology. The coefficient of determination (R 2 ) determines the correlation between the observed and simulated values. For this analysis, the F-hydra index was utilised as a replacement of the vulnerability map and plotted against the concentrations of on-ground nitrates. Awotwi et al. (2015), explicitly revealed that the coefficient of determination should be greater than 0.6 for a model to be accepted or approved. According to Figure 13, it can be established that the concentration of nitrate was closely related to the vulnerability indexes. The correlation of coefficient was found to be R 2 = 0.6792, which is slightly above the acceptable threshold (R 2 = 0.6) for any model validation. Nonetheless, the DRASTIC model exhibited a positive correlation with the coefficient of determination (R 2 ) of 0.9006 which is far above the acceptable threshold and can predict with 90% accuracy (Agyare et al., 2017). Besides, the F-hydra model can only predict with 67% accuracy. Thus, the F-hydra proves an innovative way of evaluating aquifer vulnerability, however, the results also show that further validation is strongly recommended using more measured data.

Conclusion
Aquifer vulnerability assessment is essential for decision-makers to prepare, prioritize, and plan monitoring programs. Groundwater vulnerability can be assessed using several models. The DRASTIC model is the most widely known and applicable in most geological settings. The major setback with this model is that  it requires a lot of data, making it unsuitable for areas where is scarce. In this study, an attempt has been made to evaluate the vulnerability of aquifer systems to contamination in Afram Plains. This objective was achieved using the newly established F-hydra model. The study area is classified into 3 regions on the premise of a pollution index and both models show a close correlation to the mapped regions of low and moderate, except the high. Although, the index is still under development and may not be able to accurately predict the vulnerability of aquifers in specific environments. This model proved functional and tangible to the capabilities of the DRASTIC model. The following are the overall observation and conclusions arrived from this study.
A moderate vulnerability index implies that the aquifer is susceptible to several contaminants if only there is continuous leaching or discharging. The results reveal that the moderate vulnerability region covers 51.55% (2598.12 km 2 ) of the entire area. The high and low vulnerability regions cover a significant percentage of the area 1.13% (56.52 km 2 ) and 47.32% (2384.93 km 2 ), respectively.
The sensitivity analysis results reveal that the removal of the hydraulic conductivity and land-use influences a large variation in the final vulnerability index. This result confirms how these parameters are critical to this model.
The model was validated with nitrate concentrations in production drilled wells in the study area. This study has proven that nitrate concentration increases with increasing vulnerability index. A cross plot between the mean F-hydra vulnerability indexes and the on-ground nitrate concentration exhibited a relatively high correlation of R 2 = 0.6792 which is higher than the acceptable threshold for aquifer vulnerability modelling.

Recommendation
Groundwater is the primary water supply that provides more than half of the demand for water from the populace. An investigation into aquifer vulnerability ought to be an indispensable component of the country's pursuit of setting up a system for water policy development. Although the static water table is relatively shallow, the findings recommend that topmost priority should be given to new developments which are yet to be sited in the moderate vulnerable regions. The vulnerability map developed in this study could be used by policymakers, groundwater managers and planners as a general guide to vulnerability evaluation when implementing measures to protect this fragile resource.
The F-hydra model established an innovative way of evaluating the vulnerability of the aquifer to contamination in areas with limited data accessibility. However, this model requires a colossal advancement and validation to ascertain its applicability and dependability of the results, most especially a more standardised rating of the parameters.