Geothermal renewable energy prospects of the African continent using GIS

Geothermal energy potential is identified through the deep investigation of geological, geothermal, and geophysical information over a small area; however, this task is both expensive and complex. Geographic Information Systems (GIS) can integrate different types of data (or thematic layers) over large regions and use them to identify zones of high geothermal potential. In this study, we use GIS as a tool to construct the first regional-scale geothermal potential map of Africa from different datasets. The key objective of this study is to estimate the geothermally promising areas within Africa by integrating geological thematic layers (rock units and faults), geophysical layers (heat flow derived from aeromagnetic data and seismicity), and geothermal layers (hot springs and volcanoes) within the GIS database. A weighted overlay technique within the GIS environment is applied to these data to generate the geothermal potential map. The result shows 14 regions with a high geothermal favorability index. The geothermal potential map of Africa is useful for targeting and exploring new geothermal renewable energy sites and can reduce exploration costs and pinpoint investigation areas during preliminary geothermal studies. The first geothermal potential map of Africa was constructed using GIS. A weighted overlay method was used to integrate various geo-datasets. The geothermal favorability map of Africa detected 14 high-potential zones. The first geothermal potential map of Africa was constructed using GIS. A weighted overlay method was used to integrate various geo-datasets. The geothermal favorability map of Africa detected 14 high-potential zones.

strategies to satisfy their energy demands by using renewable energy resources. Thus, geothermally active areas of Africa have great potential to provide vital energy resources to the continent. The utilization of renewable energy in general (particularly geothermal energy) has allowed African countries to shift their policies towards energy independence and reduce costly energy imports. Geothermal energy is considered to be clean, virtually inexhaustible, and with much less negative environmental impact compared to hydrocarbons.
At present, East African economies are using available geothermal resources to generate an estimated 630 MW of power annually. Kenya is the African leader for operational geothermal power plants, with its geothermal power production totaling more than 40% of the country's total electricity production. Kenya launched an investigation into geothermal energy during the late 1970s. As indicated by the Geothermal Council Resource, Kenya's progress in geothermal applications is within the top ten ranks in the world. The distribution of volcanoes (obtained from the USGS Worldwide GIS global coverage database), geysers (©OpenStreetMap), hot springs (Ait Ouali et al., 2019;Waring, 1983;Saibi et al., 2006), and installed geothermal power plants throughout the African continent. The map also depicts the structural elements and tectonic framework of the African continent, which were obtained from a geological map of Africa at a scale of 1:10 million (Thiéblemont et al., 2016) There are several stages of scientific investigation preceding the development of any geothermal project. The exploration stages include: geological field surveys (basic geology and structural geology); in situ measurements of geothermal properties (temperature, pH, and flow); geophysical surveys (gravity, magnetics, and magnetotellurics); hydrochemistry of naturally heated waters; and numerical modeling. These approaches usually require deep well data to constrain the models, which is often limited in availability and extremely costly to acquire. To overcome these limitations, in this study, we instead use a weighted overlay method based on GIS, which helps assess regions for geothermal energy potential by integrating different types of geothermal data and applying weighting to different factors connected to geothermal energy output. This approach has been used by many authors, who have demonstrated its applicability in geothermal fields around the world. For example, in Mexico, geothermal zones were evaluated by integrating geological and geophysical data using a GIS-based approach at Los Azufres (Prol-Ledesma 2020). In Indonesia, a GIS model was developed for mapping the potential of geothermal resources in West Java (Carranza et al. 2008). In Japan, the geothermal positivity map for Iwate and Akita in Northern Japan was created utilizing a GIS-based weighted overlay model applied to constrain geologically, thermally, and geochemically suitable areas (Noorollahi et al. 2007). In Turkey, the weighted overlay technique was used to map the geothermal favorability of Western Anatolia using geothermal and earthquake epicenter datasets, tectonic lineament maps, and magnetic anomaly thematic layers (TL) of evidence (Tüfekci et al. 2010). In addition, in the Akarcay Basin (Turkey), a GIS-based Minor Civil Division method was applied using surface temperature, density, and land surface proximity to current geothermal geological fields, geothermal surface manifestations, and active faults (Yalcin and Gul 2017). In Italy, a data integration tool using GIS was described to build-up a map of the geothermal favorability of Sicily, southern Italy (Trumpy et al. 2015). More recently, in Egypt, a GIS-based weighted overlay technique was applied to produce geothermal favorability map using different datasets including tectonic lineament data, temperature logs, and geophysical potential field data (Abdel Zaher et al. 2018;Abuzied et al. 2020;Abdel-Fattah et al. 2021).
The direct investigation and exploitation of geothermal resources involve potentially enormous costs due to the expense of drilling; therefore, prospecting for geothermal resources usually instead relies on indirect investigation through geological and geophysical examinations, which give a generalized perspective of the subsurface conditions. This study demonstrates the utilization of a GIS-based approach in geothermal investigation by integrating constraints from different datasets such as the distribution of geothermal manifestations (volcanoes, hot springs, and geysers), rock types, tectonics, and geophysical data (airborne magnetic and seismic activity data). These datasets are combined into several TL of evidence utilizing a GIS environment, which is used to map the geothermal potential of the African continent, showing the priority zones for possible geothermal exploitation. The input evidence layers in the GIS model are rock types, faults, surface manifestations, seismic activity, and heat flow map derived from the aeromagnetic map of Africa.
This study is, to our knowledge, the first to explore the geothermal potential of the African continent using a GIS-based model. This research will help future geothermal investigation studies by better targeting of study areas. The uniqueness and novelty of this approach is that it enables the geothermal potential to be assessed over a broad geographic area using relevant data derived from a GIS-based model, as opposed to costly and spatially limited borehole-based investigation techniques. Furthermore, the developed GIS model is flexible and can be adapted to accept future changes by integrating new geoscientific data as they becomes available and applying simple reclassifications to produce a potentiality map for any part of Africa. Site selection of geothermal resources is certainly not a simple task-it is complex and involves a diverse criterion for appropriate selection, however, our proposed model represents a powerful tool to aid in site suitability prediction.

Methodology
GIS has an extremely diverse range of applications in fields of study related to spatial information. In particular, in the context of this study, GIS-based resource evaluation is a significant and generally economical tool for distinguishing areas of interest for geothermal investigation. GIS allows different layers of spatial data to be examined and provides tools for investigating and displaying the interrelationships between TL. Recently, GIS-based multi-criteria decision analysis (MCDA) techniques have been widely utilized to produce geothermal suitability maps (e.g., Tüfekci et al. 2010;Noorollahi et al. 2007;Trumpy et al. 2015;Abdel Zaher et al. 2018;Abuzied et al. 2020). MCDA methods involve weighted and fuzzy overlay techniques. Fuzzy logic overlay analyses follow the same general steps as weighted overlay analysis, but with additional emphasis on certain steps. For example, in the fuzzy overlay method, the input is transformed (i.e., scaled) and the resulting values define the possibility of membership in multiple sets, whereas, the weighted overlay method uses a relative preference scale and weights each input accordingly (Esri 2016). The weighted overlay technique avoids the risk of disregarding potential sites because of sharp limits between spatial classes; furthermore, it gives clearer outcomes and provides reliable assessments from unsatisfactory to reasonable sites (Aliki and Hatzichristos 2019). The weighted overlay model allows various issues to be considered that are affected by multiple factors with different weights applied to each of them. The outcomes are expected to give a more in-depth assessment and provide more detailed information, as the weighting method ranks areas based on their appropriateness.
In this study, a GIS-based MCDA technique and a weighted overlay method are used to integrate various datasets, including geological, geophysical, and surface manifestations of geothermal activity. The data were processed digitally to constitute five criteria layers: distance to fault lines, geologic units, distance to surface manifestations (i.e., hot springs, volcanoes, and geysers), distance to seismic epicenters, and the value of heat flow, as derived from aeromagnetic data. An analytic hierarchy process (AHP) method was applied to weight each criterion. This is followed by combining each standardized, weighted layer, using the GIS to gain an overall geothermal potential map. The results were verified using comparison with the locations of known geothermal wells. Figure 2 shows a workflow chart of the input data and the process of defining the geothermal potentiality map of the African continent.

Criteria standardization
Standardization of criteria is the first step in preparing the raster layers. In this approach, the user-selected layers are assigned a common numeric range by classifying their attributes into discrete classes. In this paper, criteria layers were attributed to different classes on a scale from 1 to 9, based on the decision-maker's experiencebased suitability evaluation. The lower values represent less favorable areas, while higher values represent areas of greater geothermal potential.

Criteria weight assignment
An important step includes the assignment of criteria weighting coefficients. The analytic hierarchy process (AHP) is the most frequently used method for weighting coefficient calculation (Bana e Costa et al. 2016; Meng et al. 2021;Saaty 1977). Saaty's (1977) method was applied to compare the evaluation index and to study its importance ( Table 1). The established model helps to estimate the weighting for each index. A matrix is constructed according to Saaty (1994) which helps to determine the importance of the layer's elements using the following equation: where x ij , (i, j = 1, 2, 3,…, n) is the ratio of the importance of criteria i to that of criteria j.
(1) The geometric mean method (Perzina and Ramíka 2014) was used to calculate the weights of each individual criterion based on the values in the pairwise comparison matrix. The consistency index (CI) of the evaluation matrix is defined as follows: where n is the order of the pairwise comparison matrix and λ max is the highest eigenvalue of the matrix. Then the consistency ratio (CR) is calculated via the following equation: where RI is represented by average CI values gathered from a randomly generated pairwise comparison matrix. The use of CR allows one to evaluate inconsistency regarding the judgments within the comparison matrix. The suggested value of the CR should be no higher than 0.1 (Saaty 1988); if the CR value is less than 0.1, this indicates that the judgments were consistent.
The MCDA GIS method model is developed with on overlapping multiple criteria layers after multiplying the standardized layers by their corresponding weights. The adding of the processed grid layers leads to a composite value that represents the geothermal potential map for the entire study area.

Geological and structural TL
The first TL of evidence is obtained from the lithological and structural map of Africa, which provides evidence of the genesis of geothermal phenomena. The distribution of the Phanerozoic and Precambrian rocks was digitized from the geological map of Africa (Fig. 3) at a scale of 1:10 million (Thiéblemont et al. 2016). Young volcanoes or volcanic rocks of Pleistocene-Recent age may act as a heat source, while the presence of faults and fractures creates fluid pathways (Poux and Suemnicht 2012). Pleistocene-age volcanic rocks were also integrated into the geological and structural TL (Persits et al. 1997), as they have a strong correlation with the geothermal manifestations of northern Algeria and eastern Africa.
(3) CR = CI/RI, Table 1 The fundamental scale of the analytic hierarchy process (AHP) (Saaty, 1994) Intensity of importance on an absolute scale Definition If activity "i" has one of the above numbers assigned to it when compared with activity "j", then "j" has the reciprocal value when compared with "i". i and j are subscripts of "x" in Eq. (1) Regionally continuous basement exposures exist in the eastern part of the African continent along the East African Rift System (EARS) and the western coast of the Red Sea, the eastern coast of the Atlantic Ocean, in the western craton of the African continent, and East Madagascar. The crystalline basement rocks of Africa contain three significant associations: (1) the granite-gneiss-greenstone and other high-grade metamorphic assemblages which make up the Archaean cratons; (2) strongly deformed mobile belts (mostly of Proterozoic age); and (3) anorogenic intrusions and extrusive products, which include rift-related Phanerozoic magmatic and volcanic rocks (Key 1992). South Africa also has an abundance of granites and gneisses of different ages and chemical compositions, a large portion of which show anomalous radioactivity and above-average heat flow. Despite the influence of basement rock lithology, the presence of intensive zones of rock fracturing and faulting is the main control on the activity of geothermal fluids. These faults are mostly associated with recent tectonic events in the African Plate, in particular, Cenozoic faulting along tensional fractures of the EARS.
Faults and fractures are important in controlling geothermal potential since hydrothermal fluids migrate most effectively through faults (Hanano 2000). Faults on the geological map of Africa (Thiéblemont et al. 2016;Meghraoui et al. 2016) were collected and digitized (Fig. 1). Much of the geothermal activity in Africa is concentrated in the east of the continent around the EARS, where continental rifting is associated with geothermal systems with magmatic and volcanic heat sources. The EARS is one of the most major active extensional tectonic regimes on Earth; the area is  Thiéblemont et al., 2016). The geological map has been simplified to cover the major rock units of clastic sedimentary, carbonate, and Precambrian rocks in accordance with their thermal conductivity properties characterized by prominent geothermal potential spatially associated with the Quaternary volcanoes and geysers located along the rift axis.

Volcanic and geothermal surface manifestations TL
All available surface data related to geothermal activity were collected to characterize the geothermal systems. These data comprise all geothermal manifestations (i.e., volcanoes, hot springs, and geysers) (Fig. 1). The hot spring locations in Africa were gathered from diverse sources (Ait Ouali et al. 2019;Waring 1983;Saibi et al. 2006) and were merged in the GIS model. These surface manifestations were used as the fourth TL of evidence in the GIS model. Volcanic elements such as craters, calderas, and active or young volcanoes are also direct indicators for the existence of an underground heat source. The south and north of Africa have not been considered until now as feasible areas for geothermal energy production due to their tectonic stability (Dhansay et al. 2014;Enerdata 2013); however, some countries in the north and south of Africa (e.g., Egypt, north Algeria, and South Africa) have surface thermal manifestations, such as hot springs of various temperatures.
The United Nations Environment Program and the Infrastructure Consortium for Africa estimate a geothermal potential capacity of more than 20 GW of geothermal energy across Eastern Africa (Teklemariam 2018), which has encouraged countries, including Comoros, Eritrea, Djibouti, Rwanda, Uganda, and Tanzania, to conduct preliminary investigations for geothermal resources. Ethiopia has a future plan to reach 1 GW production from geothermal energy by 2021. Additionally, Uganda, Burundi, and Zambia are working to establish new small-scale geothermal power plants (Hafner et al. 2018). Table 2 summarizes the geothermal activities, temperature gradients, and heat flows of different African countries, compiled from the previous literature.

Seismic activity TL
The seismic database of Africa was obtained from numerous sources, including the United States Geological Survey (USGS) (https:// earth quake. usgs. gov/ earth quakes), the International Seismological Center (ISC) (http:// www. isc. ac. uk) and the Egyptian Seismological Center (ENSN 2014) (Fig. 4). "Seismic events" or epicenters refer to the existence of active faults which, as noted above, play an important role in geothermal systems by providing the permeability pathways required to bring waters heated at depth to close to the surface.
The distance to seismic activity (epicenters) was utilized as the fifth TL of evidence to map the geothermal potential of Africa. The African plate experiences frequent seismic events, the most recent of which are principally situated along rift zones, active volcanic fields, thrust and fold mountain belts, and offshore transform faults. Some regions in Africa are relatively aseismic, whereas, others have experienced disastrous seismic events. The most conspicuous tectonically active structures are located in the east of Africa, where the intracontinental EARS is located (Kebede and Kulhanek 1991; Ring 2014). This rift system is the source of most earthquakes in the African continent and extends from the Afar Triangle through Ethiopia to Mozambique; these areas have been volcanically and magmatically active since rifting began in the Miocene (Wafula 2011). In contrast, the western and central regions of Africa show generally subdued seismicity.

Heat flow TL calculated from magnetic data
The aeromagnetic map of Africa was used to calculate the regional Curie Point Depth (CPD) and geothermal structure in the form of the heat flow map of Africa using the spectral centroid technique. Aeromagnetic information was acquired from the various datasets that were combined under the "African Magnetic Mapping Project" (AMMP) (Getech 1992;Green et al. 1992). The CPD, which is the depth at which ferromagnetic properties of minerals change to paramagnetic due to rising temperature (> 500 °C), was estimated using the centroid technique which involves spectral analysis of the aeromagnetic data (using a window size of 1000 km 2 ). The centroid technique depends on the appraisal of separated magnetic anomalies and the statistics of magnetic groups (Spector and Grant 1970;Bhattacharyya and Leu 1977). The depth to the top of the magnetic source (Z t ) is computed from the slope of the longest wavelength portion of the spectrum, and the depth to the centroid (Z o ) is calculated from the spectrum divided by the wavenumber "|k|". The base of the magnetic source (Z b ), assumed to be the CPD, can then be obtained from the equation (Okubo et al. 1985). The temperature gradient was calculated using the basal depth Z b and the Curie point temperature of 580 °C: dT/dz = 580 °C/Z b (Maden 2010). The heat flow map of Africa (Fig. 5) was then computed using the formula: q = λ (dT/dZ) = λ (580 °C /Z b ), where λ is the coefficient of thermal conductivity (assumed to be 2.0 W/m °C, which is representative of Precambrian lower crust (Seipold 1992)). The calculated heat flow values from borehole bottom hole temperatures for 23 African countries were compared with the values derived from CPDs in order to confirm the results (Fig. 5). This comparison showed an overall coefficient of correlation of 0.7. The highest heat flows are found over the eastern part of the African continent and relate to the Afar Triangle and EARS, which comprise highly volcanic areas. High heat flows continue northwards along the Red Sea rift, but also southwards to the Kalahari craton. The northwestern margin of Africa bordering the Atlas Mountains is also a broad area of high heat flow. Individual intracontinental hotspots correspond to plume activity in the Tibesti (Chad) and Darfur (Sudan) regions.

Results and discussion
The compilation of a geothermal potential map of Africa comprises two main stages. The first stage is the identification of the main sources of evidence that predict geothermal potentiality; the second stage is then the preparation of thematic maps for each of these pieces of evidence. The five prepared thematic maps (i.e., geothermal surface manifestations, geological maps, faults, seismic activity map, and heat flow map) were used as ArcGIS model inputs. All TLs were inserted into the ArcGIS geodatabase and interpolated using the natural neighboring technique which is a weighted-average interpolation approach that generates natural neighborhood zones around each point in the dataset. The TLs were then reclassified to a common suitability scale (from 1 to 9) where higher rates indicate greater geothermal potentiality. The geothermal potential depends upon the threshold values of every parameter, as dictated by the AHP-derived weighting values.
The geothermal surface manifestation layers include a "distance to volcanic activity" layer and a "distance to hot springs and geysers" layer. Each of these layers is separated into two classes. In areas close to these features, a value of 9 was assigned; the greater the distance from the geothermal surface manifestation, the less the influence, and, beyond a threshold distance, a smaller value of 1 was applied. Noorollahi et al. (2007) showed that nearly all (97%) of geothermal wells are situated within 4 km of hot springs in Northern Japan. However, to reduce the chance of disregarding potentially promising geothermal zones, we have assumed that the geothermal surface manifestations recorded are restricted to a region of 5 by 5 km. Thus, distances less than 5 km (a similar distance to that used by Yousefi et al. (2010) for hot springs in Iran) are considered as close to a geothermal feature, resulting in the assignment of the highest value of 9 for this parameter.
The geological maps include lithological data and "distance to faults" data as two separate GIS layers. For the lithological layer, sedimentary rocks generally have low thermal conductivity values compared to volcanic and basement rocks, therefore, they are assigned a lower value on the potentiality scale. Volcanic and basement rocks have higher thermal conductivity and correspondingly higher geothermal potential (assigned as class 8 for basement rocks and class 9 for volcanic rocks). In geothermal areas, fractures and faults are crucial because they control the subsurface movement of fluids. In previous studies, Noorollahi et al. (2007) found that around 95% of geothermal wells are situated within 6 km of regional faults in Japan and Yousefi et al. (2010) used the same distance to develop a geothermal resource map of Iran. Accordingly, similar to the geothermal surface manifestations TL above, the "distance to faults" data layer is quantified as a scale of two classes-locations less than 6000 m from faults are assigned the highest value (9).
The geophysical data comprise the "distance to seismic activity" layer and the "heat flow calculated from the CPD" layer. The "distance to seismic activity" layer is treated the same way as the distance to faults layer; earthquake foci that have depths less than 10 km are classified into two classes. As the distance from the epicenter of the earthquake increases (> 6000 m), the influence on geothermal activity decreases, and, thus, a lower class is identified. Regarding the higher heat flow layer value, the more geothermally promising an area is, and higher potentiality classes are assigned accordingly.
After standardization of each layer to a scale from 1 to 9, the five criteria layers can be compared (Fig. 6). According to the previously determined evaluation ranking of the indices, the judgment pairwise matrix was attained (see "Methodology" Section) and the weight of each layer was estimated (Table 3). The calculated consistency ratio (CR) is equal to 0.1, therefore, the pairwise comparison matrix is consistent. Fig. 6 Standardization of the five criteria layers based on decision-maker's experience-based suitability evaluation (from 1 to 9), where higher rates imply a better potential for geothermal energy: A geologic units, classified into three major types: sedimentary rocks (lower rank), basement, and volcanic rocks (higher rank); B distance to fault lines, where sites less than 6000 m from fault lines are ranked top; C distance to surface manifestations, with locations less than 5 km ranking highest; D distance to epicenters, where sites less than 6000 m from epicenters are ranked top and E heat flow, in which a greater value correlates to a higher class of potentiality All five TLs were combined in the ArcGIS environment using the weighted overlay tool that multiplies several layers by their corresponding given weight values and adds them together. The geothermal surface manifestations and seismic activity are the principal indications of areas of geothermal potential because these are based on pre-existing data that are directly indicative of geothermal activity. The surface manifestations, including volcanoes, hot springs, and geysers, were assigned the maximum influence value of 38.8%. Faults and seismic activity were each assigned an influence of 20.7% per layer, as these features also directly relate to the occurrence of geothermal areas. The geology layer was assigned a weight of 10.1%, and the heat flow TLs are assigned an influence of 9.7% for each layer. The heat flow had the lowest influence as this parameter was based on data estimated using geophysical techniques. The weighted overlay of the various TL evidence produced an overall geothermal potential map of Africa with a range of values. The final geothermal potentiality map was classified using the Jenks natural breaks method which reduces the variance within classes and maximizes the differences between classes. Higher values indicate zones of high geothermal potential, whereas lower values suggest low potential zones (Fig. 7). The resultant geothermal potential map is a guide to the priority for geothermal exploration in the future.
To validate our findings, we used known geothermal wells in different countries to reliably verify the resultant geothermal potential map. Wells with high calculated heat flow values correspond to regions of high-ranked geothermal potential in our resulting map. Furthermore, our findings were compared to those of earlier studies conducted in various regions of the African continent. Limberger et al. 2018 used global heat flow data to discuss surface heat flow; their main areas of geothermal prospectivity in Africa (i.e., northwestern and eastern Africa) are broadly in good agreement with our study. A key difference between their study approach and ours is that we have also integrated constraints from surface thermal manifestations to confirm and constrain the most prospective areas.
Furthermore, Macgregor (2020) investigated the distribution of the geothermal gradient and heat flow across the African plate. His paper showed direct heat flow and geothermal gradient data from deep wells and are very localized and cannot be extrapolated to other regions that have no direct data in Africa. The current study shows geothermal heat flow data for the African continent and has the potential to help reduce exploration costs and time and, instead, allow future researchers to focus on high-potential zones. In addition, the GIS model can be further modified and updated to accommodate any new geoscientific data as they become available.

Conclusions
In total, five TL of evidence: (1) geological TL, (2) structural TL, (3) volcanic and geothermal manifestation TL, (4) seismic activity TL, and (5) heat flow TL were used to map the geothermal potential of the African continent. The different TL of evidence were assigned weights and incorporated within a GIS model to map the overall potential of geothermal resources across Africa. Note that determining the weight value Fig. 7 Geothermal potential map of Africa derived from the weighting of overlays, each representing analyses of thematic layers. The map identifies 14 areas with high geothermal potential that might be exploited for more geothermal investigation and development and criteria for each TL forms a crucial step in this approach and exerts a major influence on the resulting map. This was performed by taking into account the opinions of geothermal experts, existing literature, and the characteristics of the study area. The resulting geothermal potential map of Africa shows 14 regions of high geothermal potential, which could be used for further geothermal reconnaissance, investigation, and development. These 14 high-potential geothermal regions are located mainly in the northernmost of Africa, the Gulf of Aqaba and Gulf of Suez in Egypt, the Red Sea access, Liberia-Ivory Coast, Djibouti, the Main Ethiopian Rift, the Kenyan Rift, the Lake of Albert Rift (DRC-Uganda), the Lake Kivu Rift, northern Malawi, central Zambia, Namibia, and the Botswana-South Africa border.
The majority of Africa's current geothermal exploitation activities are related to the EARS zones that formed during the Neogene to Quaternary through normal faulting. The rift stretches from the intersection of Afar/Red Sea triple-junction through Mozambique, Malawi, Tanzania. Kenya, and Ethiopia. Investigations of geothermal reconnaissance in Eastern Africa showed a capability of 15,000 MW (Kombe and Muguthu 2019). Kenya alone has operational geothermal plants, where high enthalpy fields are situated inside the axis of the Kenya rift. However, there are future plans to raise the number of geothermal establishments in Eastern Africa by over 4000 MW through the next decade, and at present, Rwanda, Ethiopia, Djibouti, and Tanzania are making significant progress in harvesting their untapped resources of geothermal energy. In contrast, at present, geothermal activities in the North African nations are currently principally used for direct applications. This is the first study to produce a geothermal map of Africa on this scale. The approach will allow future geoscientific data to be included to update the model and the geothermal potential map has important implications for future geothermal exploration activity on the African continent.