Application of an Index-Based Approach in Geospatial Techniques for the Mapping of Flood Hazard Areas: A Case of Cape Coast Metropolis in Ghana

: Floods remain one of the disasters that destroy properties, livelihoods, and in extreme situations, take lives. As a way of prevention, geospatial applications have been employed in many cities to map ﬂood zones and predict ﬂoods. For a country such as Ghana, ﬂoods have been ranked as the second fatal disaster after epidemics leading to several kinds of research to resolve them. To date, the Cape Coast Metropolis (CCM) has received little attention in terms of research, though ﬂood cases in the area continue to escalate. This study, therefore, examines the use of geospatial techniques as tools in addressing ﬂood problems in the CCM of Ghana. From a Digital Elevation Model, hydrologic variables were generated using the ArcGIS software (Esri, Redlands, CA, USA). The soil drainage classiﬁcation for the study was generated from a downloaded African Soil Grid Drainage map, while other important factors that inﬂuenced ﬂooding in the CCM were obtained from Landsat 8 imagery. Over 21% of the CCM was classiﬁed as high ﬂood hazard zones with areas around the river Kakum estuary being ﬂood hotspots. It is, therefore, recommended that the CCM Assembly fund dredging of streams / rivers and promote a ﬀ orestation along river banks to reduce the risk of ﬂooding within the metropolis.


Introduction
Globally, floods remain among the most common and dangerous disasters that do not only displace and kill people but also destroy livelihoods, personal properties, and public infrastructure [1][2][3]. Floods often affect victims' mental health and disrupt the social functioning of communities [4]. However, research has shown that the devastating impact of floods can be reduced through the application of proactive technology [5].
Modern technological applications, notably remote sensing (RS) and geographic information system (GIS), together known as geospatial techniques, are emerging as effective methods for preventing flood hazards [5]. In contemporary societies, for instance, geospatial techniques can be used for mapping areas of flood hazard and for predicting flood occurrences [2,[6][7][8]. Accurate prediction of areas of flood hazard remains an important priority in flood risk reduction initiatives [9]. Even before a

Study Area
The Cape Coast Metropolis (CCM), a coastal city in the Central Region of Ghana, is located within longitudes 1 • 11 -1 • 41 west and latitudes 5 • 7 -5 • 20 north (see Figure 1). The CCM with a land size of nearly 122 km 2 shares boundaries with the Komenda-Edina-Eguafo-Abrem Municipality in the west, Abura-Asebu-Kwamankese District in the east, Twifo-Hemang-Lower-Denkyira District in the north, and the Gulf of Guinea in the south [34,35]. The Metropolis has a coastal stretch of about 13 km long [34].
Water 2020, 12, x FOR PEER REVIEW 3 of 16 a valuable reference document for flood control authorities and experts in the fight against flooding in the Cape Coast Metropolis and Ghana as a whole.

Study Area
The Cape Coast Metropolis (CCM), a coastal city in the Central Region of Ghana, is located within longitudes 1°11′-1°41′ west and latitudes 5°7′-5°20′ north (see Figure 1). The CCM with a land size of nearly 122 km 2 shares boundaries with the Komenda-Edina-Eguafo-Abrem Municipality in the west, Abura-Asebu-Kwamankese District in the east, Twifo-Hemang-Lower-Denkyira District in the north, and the Gulf of Guinea in the south [34,35]. The Metropolis has a coastal stretch of about 13 km long [34]. The CCM records temperatures between 24 and 32 °C with relative humidity between 60% and 80% [34]. The area experiences a double maxima rainfall receiving annual rainfall between 90 and 110 mm along its coastal stretch to 160 mm up north with the harmattan setting in between November and February [34]. The vegetative cover in the area is made up of secondary forest, comprising mainly thickets and shrubs. The CCM has a hilly topography with water bodies flowing in between them. The Kakum River and streams, such as Siwere, are located within the area [34].
The CCM with an estimated population size of 169,894 and a growth rate of 3.1% is one of the oldest districts in Ghana [36]. It was upgraded to a municipality in 1987 and subsequently received metropolitan status in 2007 due to its population size [35]. Cape Coast, the capital of CCM, historically served as the administrative capital of the Gold Coast (currently Ghana) and became a major manufacturing and commercial hub for migrants before 1877 [37]. This vibrant economic activity The CCM records temperatures between 24 and 32 • C with relative humidity between 60% and 80% [34]. The area experiences a double maxima rainfall receiving annual rainfall between 90 and 110 mm along its coastal stretch to 160 mm up north with the harmattan setting in between November and February [34]. The vegetative cover in the area is made up of secondary forest, comprising mainly thickets and shrubs. The CCM has a hilly topography with water bodies flowing in between them. The Kakum River and streams, such as Siwere, are located within the area [34].
The CCM with an estimated population size of 169,894 and a growth rate of 3.1% is one of the oldest districts in Ghana [36]. It was upgraded to a municipality in 1987 and subsequently received metropolitan status in 2007 due to its population size [35]. Cape Coast, the capital of CCM, historically served as the administrative capital of the Gold Coast (currently Ghana) and became a major manufacturing and commercial hub for migrants before 1877 [37]. This vibrant economic activity slowed down as a result of the relocation of the capital to Accra. However, its tourism and educational Water 2020, 12, 3483 4 of 16 potentials have made the metropolis a center of attraction for both foreigners and indigenes. The area has a lot of small-scale entrepreneurs and a host of government agencies [28], while farmers and fishermen constitute about 60% of the active workforce [34].

Data Sources and Processing
We downloaded a global Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) digital elevation model (DEM) at 1 arc-second spatial resolution (https://search.earthdata. nasa.gov/) for the study. A single downloaded image tile was enough to cover the region of interest. Further processing of the DEM was carried out in ArcGIS by filling all gaps in the DEM.
The Africa Soil Grid Drainage Classes map was obtained freely from (https://files.isric.org/public/ afsis250m/af_DRAINFAO_T__M_250m.tif) the International Soil Reference and Information Center (ISRIC) website. It was resampled from a 250 m resolution to a 30 m resolution to match other spatial data used in the analysis. Later, the projection was changed from the World Geodetic System (WGS) 1984 Lambert Azimuthal Equal Area to Universal Transverse Mercator (UTM) 30 N.
We acquired a Landsat 8 satellite image via the Earth Explorer website (https://earthexplorer.usgs. gov) as additional material for the study. Applying selection parameters, only one image with less than 5% land and scene cloud cover was acquired. The authors applied an atmospheric correction to improve image quality.

Setting up Data for Use in ArcGIS
We created a file geodatabase for easy saving of files under a single coordinate system. The UTM 30N Projected Coordinate System was used since Ghana falls within this zone. Therefore, all files for the study including shapefiles and images were imported into this geodatabase and were analyzed using the ArcGIS 10.4.1 software.

Digital Elevation Model (DEM)
We used a standard digital dataset that is available at the official data source. According to the works of Senevirathne and Willgoose, and Prodanovic et al. [38,39], DEMs contain topographic depressions which can take the form of a single cell (pits) or contiguous areas (depressions), and these can interrupt continuous flow paths, so we created a depressionless DEM before computing our parameters. This was achieved using the Fill tool in ArcGIS. The following variables were obtained using tailor-made Spatial Analyst tools that have in-built algorithms for processing DEMs. The flow direction raster of CCM was first created, which then became an input material for extracting the flow accumulation image. The surface drainage network was extracted using the flow accumulation raster. The stream density map was produced from generated drainage networks, while the distance from streams was also generated using the Euclidean distance tool. The elevation and slope of the study area were also generated from the filled DEM (see Figure 2).

ISRIC World Soil Information
From the downloaded Africa Soil Grid Drainage Classes' map, a subset for the CCM was generated using the study area's shapefile. This map had classified soils in the study area based on their drainage qualities on a scale from "1" being "very poor" to "5" being "well-drained".

Landsat 8 Imagery
From the Landsat 8 image, the Normalized Difference Vegetation Index (NDVI) map for the study area was produced using the NDVI tool in ArcGIS. The Landsat 8 imagery after undergoing atmospheric correction was used for land use/land cover (LULC) classification. The different land cover types were derived using the interactive supervised classification method after creating training samples. Four LULC classes were discovered, namely water, dense vegetation, sparse vegetation, and towns (i.e., built-up).
To establish whether the classification results could be used for this study, the classified map was validated against a reference image, which is a Google Earth natural color-based image. In performing the accuracy assessment, a total of 801 testing samples were used. They represent pixels of the classified image that were randomly selected and uniformly distributed across the entire study area. A pixel to pixel comparison of the testing samples and the reference image was conducted. This

Landsat 8 Imagery
From the Landsat 8 image, the Normalized Difference Vegetation Index (NDVI) map for the study area was produced using the NDVI tool in ArcGIS. The Landsat 8 imagery after undergoing atmospheric correction was used for land use/land cover (LULC) classification. The different land cover types were derived using the interactive supervised classification method after creating training samples. Four LULC classes were discovered, namely water, dense vegetation, sparse vegetation, and towns (i.e., built-up).
To establish whether the classification results could be used for this study, the classified map was validated against a reference image, which is a Google Earth natural color-based image. In performing the accuracy assessment, a total of 801 testing samples were used. They represent pixels of the classified image that were randomly selected and uniformly distributed across the entire study area. A pixel to pixel comparison of the testing samples and the reference image was conducted. This helped us to generate an error matrix that formed the input for calculating the overall accuracy and kappa coefficient (see outcome in the Results section).

Factors Included in the Flood Hazard Index (FHI)
We produced a flood hazard index of the study setting using the extracted factors (elevation, slope, flow accumulation, stream density, and distance from streams; Normalized Difference Vegetation Index (NDVI), soil drainage capacity, and LULC). These have been used in similar flood mapping studies and proven to be successful [40]. Analyzing these factors would not have been possible without first standardizing them in the same unit [41]. Therefore, using a scale of 5-20, all aforementioned flood causal factors were reclassified with "5" indicating the least possibility for a factor to influence flooding and "20" assigned for the highest likelihood. For instance, with elevation, areas with high values were given the value of 5 and vice versa because high elevation values are associated with low flooding likelihood and vice versa ( Table 1). We employed a Modified Flood Hazard Index due to the scarcity of data on CCM and the ability to modify the tool to the context of the study area [40]. For instance, stream density was considered a causal factor for floods in this study, although it was absent in the original model. The flood hazard index was derived by overlaying all the reclassified causal factors using the Raster Calculator. To promote objectivity, the authors made no judgment on the relative importance of different factors. Weights were not personally assigned to determine which factor influenced flood hazards more than the other. The following equation was used to produce the FHI.
where: Reclassification of the Index was categorized into four flood hazard zones-high, medium, low, and no-using the natural breaks method [13,40]. We verified the authenticity of the reclassified index using historical flood records of the area to arrive at the flood hazard map.

Digital Elevation Model
Elevation was higher in the northeastern section of the Cape Coast Metropolis (CCM), while areas closer to wetlands recorded the least elevation values towards the southwestern portion (see Figure 3a). The coastal stretch of the city also had low elevation values.

Slope
The slope map showed high percentage values in the northernmost and southeastern portions of the metropolis, while the southwestern portion had the lowest slope value of zero (Figure 3b). Slope and elevation maps of the study area show a similar output since the slope is a derivative of elevation.

Flow Accumulation
Flow accumulation shows the cell in an output raster where all flowing water will gather downslope. A higher value is a measure of intense flow hence higher flood hazard and vice versa [40]. The categorization revealed flow accumulation values from 0 to 1500, being the lowest, to the highest which is over 12,000 (Figure 3c). The highest flow accumulation point occurred where river Kakum drains into the ocean.

Slope
The slope map showed high percentage values in the northernmost and southeastern portions of the metropolis, while the southwestern portion had the lowest slope value of zero (Figure 3b). Slope and elevation maps of the study area show a similar output since the slope is a derivative of elevation.

Flow Accumulation
Flow accumulation shows the cell in an output raster where all flowing water will gather downslope. A higher value is a measure of intense flow hence higher flood hazard and vice versa [40]. The categorization revealed flow accumulation values from 0 to 1500, being the lowest, to the highest which is over 12,000 (Figure 3c). The highest flow accumulation point occurred where river Kakum drains into the ocean. Figure 3d shows the stream density in the study area. We found areas with high stream density from the midsection to the coastal portions in the Metropolis. The highest values recorded for stream density were along the main channel of the largest river in CCM-i.e., River Kakum.  Figure 3d shows the stream density in the study area. We found areas with high stream density from the midsection to the coastal portions in the Metropolis. The highest values recorded for stream density were along the main channel of the largest river in CCM-i.e., River Kakum.

Distance from Streams
Our findings indicate a link between distance from streams and floods. Places close to rivers were identified to be more susceptible to floods and vice versa. The data show many of the flood-prone towns are within 500 m from river networks (Figure 4a).

Land Use/Land Cover (LULC) Classification
The largest LULC types covering more than 70% of the area were vegetated areas (both sparse (35.78%) and dense (34.82%)). While dense vegetation was predominant in the northernmost part of the metropolis, sparse vegetation was concentrated in the middle portion of the area ( Figure 4d). Water was the smallest land cover type occupying barely 1.4% of the total land cover. The southern and middle sectors of the metropolis had many built-up areas (28.0%) compared with the north.  3.1.6. Soil Drainage Capacity Figure 4b reveals the drainage capacity of soils in the study area. Soils with drainage capacity classified as "moderate" covered a large portion of the area, while only a small portion near a wetland at Abakam was classified as having "very poor" drainage.

Normalized Difference Vegetation Index (NDVI)
The CCM recorded higher NDVI values up north compared with the lowest values in the south dominated by built-up environments and water bodies (Figure 4c). The period that the image was captured affected the NDVI values because in January when the harmattan period has already set in, many plants had become distressed due to the dry conditions.

Land Use/Land Cover (LULC) Classification
The largest LULC types covering more than 70% of the area were vegetated areas (both sparse (35.78%) and dense (34.82%)). While dense vegetation was predominant in the northernmost part of the metropolis, sparse vegetation was concentrated in the middle portion of the area ( Figure 4d). Water was the smallest land cover type occupying barely 1.4% of the total land cover. The southern and middle sectors of the metropolis had many built-up areas (28.0%) compared with the north.
Validation of the LULC classification was performed using an error matrix method (see Table 2). The overall accuracy and kappa coefficient were 87.86% and 0.816, respectively. A kappa coefficient of 1 signifies a perfect agreement between a classified image and a reference one, while 0 indicates no agreement [42]; hence, our classification was almost a perfect one.

Flood Hazard Zones
The resulting FHI (Figure 5a) was obtained through a combination of the causal factors earlier highlighted in the text. The FHI ranged from 40 (the least score) to 140 (the highest score). A high score is an indication of areas with high susceptibility to flooding and vice versa. To make it easy to interpret and enhance end-user understanding, it was reclassified as a color map into four categories ranging from "no" to "high" flood hazard zones (Figure 5b). This was achieved using the natural break method [11].   Our analysis revealed that high flood hazard zones within the Cape Coast Metropolis (CCM) were concentrated around the coastal areas extending inland up to the middle zone of the study area (Figure 5b). Unlike the outcome in the south of CCM, many places in the hinterlands fall within no (18.18%) or low (33.80%) flood zones. High flood hazard zones cover about 21.94% of the study area. Some of the major characteristics of the hotspots include low slope and elevation values, as well as within few meters from streams. However, both high (21.94%) and medium (26.07%) flood hazard zones constitute 48% of the entire study area.

Flood Hazard Index Validation
We validated the resulting flood hazard index (Figure 5a) using historical flood records from the Central Regional Secretariat of the National Disaster Management Organization (NADMO) to ascertain its authenticity. Historical flood records are relied upon for verifying the accuracy of flood hazard maps generated using GIS [11,40]. Table 3 shows the list of communities in CCM classified as flood-prone. When we plotted these communities on the flood hazard map, as shown in Figure 5b, the outcome confirmed seventeen percent as correctly classified under a medium flood hazard zone with the rest accurately classified as belonging to high flood hazard territory. This proves the success of our technique in mapping potential flood areas. North Ola/Ameen Sangari 12,210 Source: Central Regional National Disaster Management Organization (NADMO) Secretariat, Cape Coast.

Discussion
Based on geospatial techniques, this study mapped flood hotspots within the Cape Coast Metropolis (CCM) using an index-based approach. To the best of our knowledge, this is the first study to employ an index-based approach to mapping flood hotspots for the study area. The finding that the midsection of the metropolis down to the coastal areas was within medium and high flood hazard zones should be a great concern for both flood management institutions and residents of the CCM. As many communities within the high flood intensity areas were located close to streams as observed in Kwapro, Ankaful, Amamoma, and Duakor (please see Figure 6), much consideration should be given to the relocation of residents close to these streams or building a defense block to prevent flood disaster in the future. Our suggestion corresponds with conclusions in some previous studies indicating that high flood risk zones are usually close to river bodies [11,41] and that relocation of residents can be a life-saving opportunity for people living in such high flood hazard areas [4,41].
Additionally, the dredging of streams can be an effective way of reducing flood hazards in the CCM as in the case of 2016 where concerned authorities sponsored the dredging of the Kakum River at Kwapro, leading to reducing flood cases in Kwapro (Figure 7).
The findings from this study further suggest that flood hazards in the CCM result from the interplay of several factors: elevation, slope, flow accumulation, stream density, and distance from streams; Normalized Difference Vegetation Index (NDVI), soil drainage capacity, and Land use/Land cover (LULC). Areas with low elevation, small slopes, high stream density, low NDVI values, and near the highest flow accumulation were found to be the main hotspots for flooding and vice versa, and these findings confirm the conclusions in previous studies elsewhere [6,13].
Water 2020, 12, x FOR PEER REVIEW 12 of 16 Figure 6. Dwellings located close to the Kakum River in the study area (Credit: Central Regional NADMO Secretariat, Cape Coast).
Additionally, the dredging of streams can be an effective way of reducing flood hazards in the CCM as in the case of 2016 where concerned authorities sponsored the dredging of the Kakum River at Kwapro, leading to reducing flood cases in Kwapro (Figure 7). The findings from this study further suggest that flood hazards in the CCM result from the interplay of several factors: elevation, slope, flow accumulation, stream density, and distance from streams; Normalized Difference Vegetation Index (NDVI), soil drainage capacity, and Land use/Land cover (LULC). Areas with low elevation, small slopes, high stream density, low NDVI values, and near the highest flow accumulation were found to be the main hotspots for flooding and vice versa, and these findings confirm the conclusions in previous studies elsewhere [6,13].  The findings from this study further suggest that flood hazards in the CCM result from the interplay of several factors: elevation, slope, flow accumulation, stream density, and distance from streams; Normalized Difference Vegetation Index (NDVI), soil drainage capacity, and Land use/Land cover (LULC). Areas with low elevation, small slopes, high stream density, low NDVI values, and near the highest flow accumulation were found to be the main hotspots for flooding and vice versa, and these findings confirm the conclusions in previous studies elsewhere [6,13].
Additionally, distance to stream networks played a crucial role in the flood hazard index; many of the communities classified under high flood territories were within a 500 m buffer of river bodies. This finding corresponds with the assumption that distance from river channels plays a critical role in flood inundation; therefore, river overflows should be considered since most rivers overflow their banks when water exceeds the carrying capacity of the river channel [40].
In line with the findings of Nyarko [23], land use played a contributory role in flooding in the Cape Coast Metropolis (CCM), and authorities in the CCM need to discourage deforestation and promote afforestation among residents of the locality. Moreover, it will be prudent for the authorities within CCM to devise strategies to protect wetlands from further encroachment as the resulting flood hazard map showed that most of the wetlands in the CCM are inhabited by people. Additionally, distance to stream networks played a crucial role in the flood hazard index; many of the communities classified under high flood territories were within a 500 m buffer of river bodies. This finding corresponds with the assumption that distance from river channels plays a critical role in flood inundation; therefore, river overflows should be considered since most rivers overflow their banks when water exceeds the carrying capacity of the river channel [40].
In line with the findings of Nyarko [23], land use played a contributory role in flooding in the Cape Coast Metropolis (CCM), and authorities in the CCM need to discourage deforestation and promote afforestation among residents of the locality. Moreover, it will be prudent for the authorities within CCM to devise strategies to protect wetlands from further encroachment as the resulting flood hazard map showed that most of the wetlands in the CCM are inhabited by people.

Policy Recommendations
Based on our findings, some recommendations are provided: i.
The Cape Coast Metropolitan Assembly (CCMA) should fund the dredging of rivers and streams as shown in Figure 7 every two years because the periodic dredging of river networks could salvage the situation. This would require huge financial allocation from the local assembly but would go a long way to reduce flood cases in the Metropolis. ii.
Considering the role of vegetation in reducing runoff, it will be prudent on the part of city authorities to plant trees along riverbanks to serve as a buffer. This should be done with the total participation of the local folks since their involvement would prove crucial in ensuring success. iii.
NADMO, CCMA, and Ghana Police should work in tandem to ensure strict compliance with buffer zone policies. Authorities must enforce such laws devoid of political interference or favor to deter residents from building close to riverbanks as shown in Figure 6. iv.
The Assembly should make the necessary effort to construct concrete banks to prevent sediments from filling the river channel. v.
Residents within the CCM especially those living in flood hazard zones should take education and early warning communication by disaster management organizations seriously to minimize the effects of floods should they occur. vi.
Moreover, residents in CCM should construct houses with recourse to building codes since many communities lay within high flood hazard zones.

Limitations of the Study
Although this study demonstrates novel contributions to the mapping of flood hazards using geospatial technology, some limitations should be noted. Firstly, the study excluded direct rainfall from the factors responsible for flooding. This also comes at the backdrop that the CCM has a fairly good drainage capacity as shown in our findings-the soils in CCM allow water to percolate and did not influence high runoff. Secondly, we did not capture the spilling of excess water from the Brimsu dam in the metropolis. Excess water from this dam is commonly spilled over when the dam reaches its peak capacity. This exercise is intended to protect the dam from collapse but ends up increasing flood risk in the neighboring communities, such as Ankaful. The spilling of the water was not included in the list as this is usually a sponsored venture and the quantification of the amount of water spilled is difficult to measure.

Conclusions
Our study explored the use of remote sensing and GIS in mapping flood hotspots within the Cape Coast Metropolis (CCM) using an index-based approach. The study has contributed to the knowledge that simple mapping techniques can be used to determine flood hotspots in areas where data and expertise are major concerns. This study has identified that elevation, slope, flow accumulation, stream density, NDVI, distance from streams, and LULC are key contributory factors to floods in the CCM, and preventive measures can be focused on these factors. It is also worth noting that based on our method, the soil in CCM has been found as a less contributory factor to flooding in the CCM.
To prevent future flood disasters, however, authorities and residents of the CCM can utilize the new findings based on our flood hazard map for the effective planning of the metropolis. We believe that the suggested recommendations if carefully implemented would save lives, properties, and historical monuments in the CCM. The findings and approach used in this study can also be useful for future studies aimed at mapping other flood hazard areas in Ghana and even beyond.