2. Geodynamic and Hydrogeological Context
The westward propagation of the Gulf of Aden at a velocity of 30 mm/year gave rise to the emergence of the Asal Rift [
2]. The opening of the Asal Rift segment dates from the Miocene and its geology is characterized by the abundance of volcanic rocks formed in an environment of extension and opening. This emerged rift segment is 15 km long and 14 km wide. The opening velocity is approximately 16 mm/year [
6] in the direction of N°45° E. The structure of this rift is characterized by a set of fracture networks and normal faults in the N130°E direction between Asal Lake and Ghoubbet Bay (
Figure 1). Having a thin thickness and basaltic lavas at the surface in the form of basalt series, the land crust of the Asal Rift is considered transitional, i.e., an intermediate between the continental and nascent oceanic crust [
4]. The first effusive event linked to the birth of the rift dates from 853,000 to 315,000 years ago [
7]. Effusive magmatic activity continued in the northern part of the Asal Rift for 315,000 years and the southern part has been under seawater for 326,000 years, with the formation of hyaloclastites [
7]. The activity of the central Fieale Caldera controlled the evolution of the rift between 326,000 and 100,000 years ago. Indeed, it allowed the injection of large volumes of basaltic lava into the interior floor that masked the previous faults [
8]. Over the last 50,000 years, volcanic activity has decreased and the successive basaltic flow formations that make up Fieale Caldera have gradually been shifted by normal faults [
7]. The modern rift structure began about 40,000 to 30,000 years ago with the development of the faults with outer margins H and δ [
9]. From this period, the Fieale Caldera gradually extinguished and volcanic activity continued along the interior floor with small volcanic edifices and eruptive fissures [
7]. With the most recent magmato-tectonic event dated in November 1978, two earthquakes of magnitude 5 and 5.3 were recorded in the Asal-Ghoubbet Rift and, following this event, 0.7 m of vertical subsidence and 2 m of horizontal extension in the direction N °40 E occurred in the rift [
9]. A one-week eruption accompanied by basaltic fissures to the northwest of the rift gave rise in November 1978 to the axial volcanic chain and then to the eruption of the Ardoukoba volcano (
Figure 1). At present, most of the rift deformation is concentrated in its north-eastern part and around the edifice of the Fieale Caldera [
4,
5,
6].
Thermal springs, fumaroles, and the existence of alteration on the surface of the Asal Rift are indicators of volcanic activity and the presence of potential geothermal resources at depth. Hot springs are abundant around the Asal Lake, especially in the eastern part, with temperatures ranging from 30 to 90 °C (
Figure 1). Between the Asal Lake and the Ghoubbet Sea, only fumaroles are observed at the surface because the hydrostatic level is below the topographic surface.
Chemical analyses of gases and determination of δD and δ
18O water stable isotopes of steam condensate have been performed for some of these fumaroles [
10,
11]. The results allowed estimating the source temperature at a depth higher than 230 °C, and even than 300 °C, in the studied Asal Rift. The isotopic data from the steam condensate suggest the existence of primary steam (likely separated at high temperature) originating from a deep heavy brine in this area [
11]. The existence of this primary steam is in good agreement with the presence of a dominant vapor zone in the conceptual geothermal model proposed by BRGM [
10,
11,
12].
Owing to the high salinity of Asal Lake and according to geochemical analyses, about 90% of the recharge of Asal Lake is believed to be due to infiltration and evaporation of Ghoubbet seawater [
3], depleted in sulfate because of CaSO
4 precipitation [
13]. From a geochemical study of water sampled from the rift thermal springs, the authors of [
3] showed that the chemical and Sr isotope composition of these waters can be explained by variable seawater interaction with basaltic rocks at different temperatures and evaporation processes. Among thermal waters, those of the Manda group (
Figure 1) indicate the lowest seawater–basalt interaction and evaporation degree, suggesting that their source temperature is not high and their circulation rate is relatively significant between the Ghoubbet Sea and the Asal Lake. The thermal waters of the Eadkorar group, located in the northeastern part of the lake (
Figure 1), constituted by a mixing of seawater, Asal Lake water, and meteoric water, as suggested by the δD and δ
18O water stable isotopes [
14,
15], indicate the highest seawater–basalt interaction grade and seem to be the only direct leaks of a deep fluid discharged from a geothermal reservoir at 200–210 °C, as estimated by chemical geothermometry [
3]. The Korilii and Kalou thermal waters (
Figure 1) can result from a mixing between the deep geothermal water at 260–270 °C with seawater and Asal Lake water. The stable isotopes suggest the Korilii water is also diluted by meteoric waters [
14,
15]. The Eounda Alifitta waters, located in the north-eastern most part of the lake (
Figure 1), indicate relatively low basalt–seawater interaction degrees and significant contributions of Asal Lake water. A contribution of meteoric water is also suggested by the water stable isotopes [
15].
The geothermal wells A1, A3, and A6 are within a radius of 300 m (
Figure 1). These wells have confirmed the existence of an intermediate reservoir located between 300 and 700 m depth, with a maximum temperature of 185 °C, and another potentially deep reservoir located between 1200 and 2000 m depth, with a maximum temperature of 260–280 °C (
Figure 2). However, production tests proved to be ineffective for a continuous commercial geothermal energy production over time for two reasons:
The authors of [
3] showed that the deepest geothermal fluid from the wells is constituted by a mixing of seawater and Asal Lake water, interacting with basalt rocks at 260–270 °C. The water stable isotopes also suggest a contribution of meteoric water [
15], similar to that observed in the Korilii, Eadkorar, and Eounda Alifitta thermal waters, which could mainly be the meteoric water that transits from Asal local rift to Sakalol-Harralol depression, located at the north-west of Asal Lake, or/and the meteoric water coming from the deep circulating regional aquifers [
11,
12,
13,
14]. The low boron concentration observed in these deep waters suggests the existence of a steam phase in the geothermal system [
3]. The existence of such a zone with dominant vapor is in agreement with the conclusions relative to the gas and steam condensate data drawn by [
11] and the conceptual geothermal model proposed by BRGM [
10,
11,
12]. The equilibrium temperature calculated from all gaseous reactive species associated with these well waters, except H
2S, is very close to that measured at the bottom of the geothermal wells Asal 3 and Asal 6. This temperature corresponds to that of the identified geothermal reservoir and varies between 250 and 270 °C. According to [
3], the geothermal waters located at a depth between 300 and 700 m correspond to a mixing between the deepest geothermal fluid and seawater. The circulation of hydrothermal fluids in the northern part of the rift is believed to happen deeper than in its southern part [
17].
Assuming the highly saline water of Asal Lake contributes to the recharge of the geothermal reservoir, this contribution is relatively recent as the water of Asal Lake was fresh and its salinity was recently acquired [
3,
4,
5,
6,
7,
8,
9]. This constrains the fluid circulation time in the deep hydrothermal system located between 1 and 3 km, having high enthalpy. The water level of Asal Lake 9 to 6 thousand years ago was 160 m above sea level. Over the past 5 to 6 thousand years, the level has declined by 310 m and the lake level is currently 150 m below sea level [
9]. This level has remained constant for the last centuries. Therefore, this suggests that, despite the intense evaporation from the arid climate of the area and the infiltration of water from the lake, this evaporation of lake water is compensated by a constant recharge mainly supplied by the Ghoubbet Sea.
The Asal 2 and Asal 3 wells are believed to penetrate the potential deep reservoir unit, but Asal 2 well was non-productive, although a temperature of 230 °C was measured at a depth of 923 m, while the Asal 3 and Asal 6 (located with a distance of about 300 m) were producers with permeable zones at a depth of 1030 m (
Figure 2). The fluid produced was hyper saline and essentially composed of liquid at 260° C [
16].
5. Discussion
The directions of faults and fractures are parallel to the direction of believed groundwater flow (
Figure 1 and
Figure 4). Unfortunately, the available MT data do not cover the entire area of the emerged Asal Rift. A lack of data in the area of the normal faults ε and δ to the north-east of the rift limits the interpretation (
Figure 1). However, previous geochemical analyses conducted by [
36] in littoral thermal springs could help to place our interpretations in a logical context and may explain the origin and the evolution of thermal waters that are not far at the eastern flank of the normal faults ε and δ located at the north-east of the rift. These authors conclude that a mixture of hot seawater with cold seawater emerges through the fractures. Moreover, the presence of deep hot seawater-derived geothermal fluid with a temperature of 210 °C was confirmed [
37]. Despite the lack of data in this zone, our interpretation of pronounced fluid circulation at the north-east of the rift can be justifiable. The comparison between the
87Sr/
86Sr ratio and strontium concentration of the brine waters from the Red Sea bottom and Asal Lake have indicated that Asal Lake is a better location to study how hydrothermal fluids of seawater origin may evolve in the Red Sea rift [
38].
The descriptive lithology of the wells drilled in the Asal Rift consists mainly of a succession of basaltic series, which were deposited over geological time and have different ages [
39]. Even if the abundance of shale and/or graphite is known to cause a decrease in electrical resistivity [
21], no significant quantity of the latter is found in the wells drilled, and more precisely between 1 and 2 km depth [
39]. Within the Earth’s crust, the dominant electrical conduction mechanism is electrolytic [
21]. Therefore, low electrical resistivity in the context of the Asal Rift is generally interpreted as an indicator of the presence of a fluid. At depths of 1200 and 1400 m (
Figure 4), groundwater seems to be present mainly between the H fault and the δ fault according to the following evidence. In this zone, observed fractures were caused by the tectonic and volcanic activity associated with the eruption of the Ardoukoba volcano that occurred in 1978 (
Figure 1). The low electrical resistivity suggests the presence of water and hence more important flow to the southeast between the H and F faults, and to the north-east of Fieale Caldera. The existence of dykes at a depth near the rift axis (between the F and β faults) has been hypothesized based on the deformations observed at the surface [
9]. In this case, rift extension and subsidence are controlled by dyke injection. Moreover, the stretching model of Asal Rift is likely accommodated by the magma activity [
40]. This extension can result in open cracks where groundwater flows easily. The proximity of the sea and the coexistence of recent open faults and fractures is a fair evidence for such flow illustrated with the simulations. The electrical resistivity model, especially images obtained at a depth between 1.2 and 2 km (
Figure 4), indicates a possible hydrothermal system in the north-east part of the Asal Rift and with more pronounced regional groundwater flow in this area rather than in the south-west zone (
Figure 4). Considering the contrast in electrical resistivity between the north-east and south-west parts of the Asal Rift, it seems likely that the hanging wall of the H fault acts as a hydrogeological barrier separating the Asal Rift between two zones. A first conductive zone, whose high electrical conductivity is due to the presence of groundwater, as well as a second less conductive zone with an increase of the electrical resistivity gradient towards the south of the rift, imply lower porosity and/or less groundwater in pore space and fractures. Only the major J fault and the minor fault passing through well Asal 4 exist in the south zone of the Asal Rift. In addition, according to the bathymetry map of the region, the J fault is not connected to the Ghoubbet sea. The absence of the open cracks or fractures observed with the November 1978 eruption in the southern zone of the Asal Rift is another argument that limits the possibility of possible fluid flow or water infiltration in this sector, while the presence of activated faults in 1978 in the northern part is dominant (
Figure 1).
Wells A5, F1, and F2 were drilled in the north-east zone of the H fault, but only the temperature profile of well A5 was available for this study (
Figure 2). This temperature profile shows a first convection cell from 300 to 500 m and a second convection cell from 500 to 1200 m, and then a constant geothermal gradient that can be related to heat conduction between 1.2 and 1.8 km. This information, therefore, correlates with the conductive zone to the north-east of fault H. Furthermore, the geothermal gradient estimated between 1.2 and 1.8 km in well Asal 5 is 18.7 °C/100 m. In addition, well Asal 5 is located in a relatively conductive zone with an electrical resistivity greater than 20 Ω m, which can corroborate the absence of groundwater that can have a high salinity associated with electrical resistivity less than or equal to a value of 10 Ω m [
21]. This Asal 5 well can be close to a heat source where heat is transferred by deep conduction, resulting in a strong conductive geothermal gradient at a depth between 1.2 and 1.8 km. The absence of groundwater circulation between 1.2 and 2 km in this well is corroborated by both its temperature profile and the electrical resistivity model presented in this study. The latter also shows the Fieale Caldera in close proximity to well Asal 5 has a particular electrical resistivity range of 20–30 Ω m. It is highly probable that this resistivity is characteristic of the nature of the rock in this Caldera zone, which can be hot, but less permeable, limiting the circulation of fluid. Another important fact is that hydrothermal mineral assemblage is in good agreement with the measured temperatures in all geothermal wells of Asal Rift, except for well Asal 5 [
39], where a chlorite-epidote mineral assemblage zone that normally should exist at a high temperature was observed at a depth between 500 and 1200 m, where an important temperature inversion was measured (
Figure 2). This information is helpful allows to assume that a pre-existing high temperature reservoir can have existed before the infiltration of cold sea water, and the temperature inversion may likely be related to a recent infiltration of seawater.
Geodetic measurements indicate that the opening of the Asal–Ghoubbet Rift has a velocity of 16 mm/year and an opening direction perpendicular to the rift axis, highlighting the asymmetrical behavior of the Asal Rift [
6]. These authors studied the horizontal and vertical deformation of the rift in directions parallel and perpendicular to the rift axis. They consider that the north-eastern part of the Asal Rift is the site of significant continuous deformation, whereas for its south-western part, the deformation is weak and practically non-existent. These results corroborate the presence of active normal faults in the north-eastern part of the rift and the existing seismic activity [
4]. The rate of deformation along the axis perpendicular to the rift axis increases towards the north-east of the rift [
6]. This is a surprising correlation between the increase in this deformation rate and the increase in electrical conductivity highlighted by the model developed in this work pointing in the same direction (
Figure 4).
The regional groundwater flow affected by sea water was believed to be the main source of recharge to Asal Lake [
3], but our simulation results show the groundwater flow is directed toward both the Asal Lake and the Ghoubbet Sea, with a maximum point on the SE side diving flow in the two directions and a pronounced flow to Asal Lake. This can limit the pontifical depth of penetration of sea water intrusions near the surface. In this case, recharge of Asal Lake can occur through aquifer (C1), potentially containing fossil seawater that is likely fed by a deep hydrothermal circulation. In all the different simulation scenarios, it is important to highlight the evidence of groundwater flow into Asal Lake and hydrothermal circulation beneath Asal Lake and the Ghoubbet sea.
The superposition of three permeable zones with the predominance of flow towards Asal Lake and the Ghoubbet sea are important elements to advance the hypothesis of the existence of a large hydrothermal system under the rift that would be compartmentalized into these three permeable zones, which can be individualized under the tectonic activity and the injection of a deep heat source since the formation of the rift over the geological time. The subsurface flow at a depth between 300 and 700 m (
Figure 7) can be influenced by the Ghoubbet sea and the local rift topography, while the intermediate flow at a depth between 1200 and 2000 m (
Figure 7) is interpreted as the result of the development of a hydrothermal system, and the deep flow can be linked to a deep hydrothermal circulation feed by a heat source originating from a deep magmatic system (
Figure 7). The presence of the cap rock layer limits possible communications between the second permeable layer, identified here as a potential deep reservoir, and the shallow aquifer C1 (
Figure 5 and
Figure 7), and facilitates the development of a hydrothermal reservoir in the formation, where this cap rock layer covers and helps to contain thermal energy in the potential deep reservoir.
Similar simulation results obtained with scenario 5 (
Figure 7) support the idea that regional flow variation is not strongly influenced by the permeability anisotropy in layer C1 of the conceptual model (
Figure 5). In other words, the regional flow is mainly controlled by tectonics that affect deep structures and allow the deep heat source to be directed to shallower geological strata and formations. Simulations with an isotropic permeability of 10
−15 m
2 inside the deep reservoir and equal to 10
−17 m
2 inside the heat source were conducted for scenario 3 to indirectly evaluate the vertical dependence of permeability, which can decrease with depth. The development of hot convection cells within this unit was less pronounced than in scenarios 1 and 2 (
Figure 6 and
Figure 7). Scenarios with permeability dependent temperature (scenarios 1, 2, and 5) seem to be more realistic than scenario 3 with isotropic permeability. However, the results from scenario 3 are not representative of conceptual hydrothermal models proposed in previous studies [
41,
42], where hydrothermal activity is concentrated in areas with elevated topography. Hydrothermal activity in permeability dependent temperature scenarios better correlates with the central zone of high topography (
Figure 6 and
Figure 7). In scenarios 1, 2, and 5, the lateral extension of the hydrothermal reservoir with the isotherm 210 °C becomes important compared with the results of scenarios 3 and 4 (
Figure 6 and
Figure 7). This pattern is in good agreement with a previous study conducted in a tectonically active rift-ridge zone in Iceland [
43]. The authors of this study confirm that emplacement and geometry of the upflow zones are mainly controlled by the location and the permeability values of fault zones. The geometry, location, and permeability values of fault zones remain the same in all the different scenarios for this study, but we deduce from our results that the permeability dependent temperature rock formulated in both heat source and potential deep reservoir can play a major role to reproduce a representative dynamic heat transfer in hydrothermal systems. Permeability dependent temperature formulated solely for the heat source unit may not produce significant variation when the potential deep reservoir is set with a constant permeability like the simulated conditions of scenario 4. In addition, the results are approximately the same as the results of scenario 3, where isotropic permeability was set for both the heat source and potential deep reservoir units. In such cases, we hypothesize that permeability dependent temperature would be an appropriate formulation in a hydrothermal system as in the Asal Rift, when it is adopted for both the heat source and potential deep reservoir or, uniquely, for the potential deep reservoir like in scenario 2.
The comparison between the temperature profile measured in the Asal 6 well and the temperature simulated appears more realistic with scenarios 1, 2, and 5 compared with other scenarios (
Figure 8). This indicates that the conditions used are representative of the rift and the associated thermal anomaly can be related to the injection of a magma forming a heat source and the development of vertical conduits qualified as active faults to convey this internal thermal energy. In scenario 4 (
Figure 8), the system at the end of the simulation began to cooled down, as the temperature at 100,000 years is less than the temperature measured in well Asal 6. It can be reasonable to interpret this as a consequence of an extinguished deep hydrothermal fluid circulation inside the heat source layer, where permeability could be higher than that assigned in scenario 4. However, the shape of the simulated temperature profile at the location of well Asal 6 is maintained and has the same appearance as the measured temperature in well Asal 6. Furthermore, the same reservoir was recognized at a depth between 1050 and 1300 m in three wells, Asal 1, 3, and 6, and the geologic formation hosting this reservoir is called Dalha basal, which has an average age of 4 to 9 My, where temperatures measured in this reservoir range from 260 to 280 °C [
33] and are close to the simulated temperatures in this study. In addition, permeability anisotropy considered in scenario 5 inside the C1 layer does not lead to cooling of the system and the simulated thermal state at well Asal 6 is close to the measured one (
Figure 8). It can be assumed that the thermal anomaly caused by this near-surface groundwater flow has less effect at the global scale of the system. Nevertheless, the assumption of pronounced horizontal flow near Asal Lake in scenario 5 is corroborated by our interpreted conceptual model based on electrical resistivity inversion (
Figure 5), where the thickness of the shallow aquifer called C1 increases toward Asal Lake. As results of MT inversion and numerical simulations showed, this study is focused for the potential deep reservoir between 1 and 3 km. The developed 2D model simulate a regional groundwater flow that gives a capture of fluid circulation in a regional scale. Then, the lack of a better correlation between measured and simulated temperature at Asal 6 well for the shallow part of the system above 600 m of depth (
Figure 8) may be related to a localized fractures or convections cells. Another possible interpretation would be resulted by the difference of density fluid between the upper part of the system and the lower part. May be if a higher density fluid (as saline fluid) was considered in the shallow part than the deeper part, the shallow part of temperature profile measured at Asal 6 well could be well represented.
Temperatures measured in the thermal springs emerging on the periphery of Asal Lake range from 30 to 90 °C [
3], and the results of the different scenarios show approximate emergence temperatures that corroborate with those measured in the thermal springs (
Figure 7). Based on the chemical and isotopic characteristic of water in the Asal region, three main types of water may co-exist in the subsurface: seawater intrusion from Ghoubbet sea, Asal Lake water, and meteoric infiltration probably originating from the meteoric water that transits from Asal Rift to Sakalol-Harralol depression or/and the meteoric water coming from the deep circulating regional aquifer [
14]. The low temperature emergence observed in the hot springs of group Manda [
3], their chemical composition of which is close to that of seawater, and the high rate of flow show that seawater circulation through the faults and fractures is the main source of their recharge [
15]. This information can be related to the emplacement of the hot springs group Manda that are close to the H fault and are aligned more or less to this H fault and with two other major faults located south of the H fault (
Figure 1). Moreover, in this case, the information corroborates our hypothesis that seawater intrusion occurs where faults/fractures exist (
Figure 1 and
Figure 4).
Dissolved salts in water can change the phase conditions of liquid–vapor system and the phase relation of saline water approximated by the H
2O-NaCl system is not the same as those of pure water considered in this study. In conditions of an H
2O-NaCl system such as that of the Asal Rift, the temperature and pressure of liquid–vapor co-existence are above the critical values of pure water [
44]. Despite that the fluid simulated from HYDROTHERM code is pure water [
19], for saline geothermal fluid, the emplacement depth of intrusion controls significantly whether phase separation is dominantly carried out by boiling or by condensation [
35]. Obviously, in this case, numerical simulations considering a saline geothermal fluid are expected to reveal the presence of a vapor phase, as proposed by the BRGM conceptual model [
17]. This model is in good agreement with the gas chemical results of the Asal fumaroles and the water stable isotope data of their steam condensates [
11], as well as the decrease of boron concentrations in the well geothermal fluids evidenced by geochemical analyses [
3].
The total flow rate of geothermal fluid produced in well Asal 3 showed a considerable decrease of production rate and a decrease of bottom hole pressure, which were related to sulphide and silica deposits in the well with a total dissolved solids (TDS) of 116,000 ppm [
33], suggesting a high potential of permeability reduction.
The δD and δ
18O water stable isotope data for the geothermal fluids of the Asal wells and for some hot springs (Korilii, Eadkorar, and Eounda Alifitta groups;
Figure 1) indicated contributions of meteoric water [
3]. This meteoric contribution probably originates from the meteoric water that transits from Asal Rift to Sakalol-Harralol depression or/and the meteoric water coming from the deep circulating regional aquifer [
14]. These authors showed, in the area of Sakalol-Harralol, the presence of hot springs, with high TDS Na-Cl waters, aligned in NW-SE along the main faults like in the Asal Rift. The high altitude zones of rift margins at the north and south do not seem to be potential recharge sources. However, a more representative numerical model can incorporate these zones and be tested to confirm or invalidate these hypotheses. As suggested by [
14], the existence of a deep regional aquifer that may have a wide extension in Djibouti and where the Asal geothermal water has apparently common features in terms of chemical and water-isotope compositions is a good argument to highlight the need for a numerical multiphase model developed at a large regional scale. Coastal aquifers in Djibouti hosted geothermal water that could be more evolved in terms of water–rock interaction between deep fossil water and sea water intrusions [
38].