Journal cover Journal topic
Hydrology and Earth System Sciences An interactive open-access journal of the European Geosciences Union
Journal topic
Hydrol. Earth Syst. Sci., 23, 763-771, 2019
https://doi.org/10.5194/hess-23-763-2019
Hydrol. Earth Syst. Sci., 23, 763-771, 2019
https://doi.org/10.5194/hess-23-763-2019

Research article 11 Feb 2019

Research article | 11 Feb 2019

# Faulting patterns in the Lower Yarmouk Gorge potentially influence groundwater flow paths

Faulting patterns in the Lower Yarmouk Gorge
Nimrod Inbar1,2,3, Eliahu Rosenthal1, Fabien Magri4,5, Marwan Alraggad6, Peter Möller7, Akiva Flexer1, Joseph Guttman8, and Christian Siebert9 Nimrod Inbar et al.
• 1Porter School of the Environment and Earth Sciences, Tel Aviv University, Tel Aviv, Israel
• 2Department of Geophysics and Space Sciences, Eastern R&D, Ariel, Israel
• 3Department of Physics, Ariel University, Ariel, Israel
• 4UFZ Helmholtz Centre for Environmental Research, Dept. Environmental Informatics, Leipzig, Germany
• 5Freie Universität Berlin, Chair of Hydrogeology, Berlin, Germany
• 6Water, Energy and Environment Center, University Of Jordan, Amman, Jordan
• 7Helmholtz Centre Potsdam German Research Centre for Geosciences – GFZ, Section 3.4, Telegrafenberg, Potsdam, Germany
• 8Mekorot, the National Water Company, Hydrology Department, Tel Aviv, Israel
• 9UFZ Helmholtz Centre for Environmental Research, Dept. Catchment Hydrology, Halle, Germany
Abstract

Recent studies investigating groundwater parameters, e.g., heads, chemical composition, and heat transfer, argued that groundwater flow paths in the Lower Yarmouk Gorge (LYG) area are controlled by geological features such as faults or dikes. However, the nature of such features, as well as their exact locations, were so far unknown. In the present paper, we propose a new fault pattern in the LYG area by compiling and revising geological and geophysical data from the study area, including borehole information, geological map cross sections, and seismic data from the southern Golan Heights and northern Ajloun mountains. The presented pattern is composed of strike–slip and thrust faults, which are associated with the Dead Sea transform system and with the Kinnarot pull-apart basin. Compressional and tensional structures developed in different places, forming a series of fault blocks probably causing a non-uniform spatial hydraulic connection between them. This study provides a coarse fault-block model and improved structural constraints that serve as fundamental input for future hydrogeological modeling which is a suggested solution for an enigmatic hydrological situation concerning three riparian states (Syria, Jordan, and Israel) in a water-scarce region. In areas of water scarcity and transboundary water resources, transient 3-D flow simulations of the resource are the most appropriate solution to understand reservoir behavior. This is an important tool for the development of management strategies. However, those models must be based on realistic geometry, including structural features. The study at the LYG is intended to show the importance of such kinds of structural investigations for providing the necessary database in geologically stressed areas without sufficient data. Furthermore, during the hydrogeological investigation, a mismatch with results of pull-apart basin rim fault evolution studies was discovered. We argue that this mismatch may result from the settings at the eastern rim of the basin as the en-echelon changes from pull-apart basins (Dead Sea, Kinnarot, Hula) to a push-up ridge (Hermon).

1 Introduction

Figure 1Regional location map on top of a digital elevation model, presenting the faults at the eastern rim of the Dead Sea Transform (DST), Sheikh Ali Fault (SAF), Nov Fault Zone (NFZ), wells belonging to Meizar (MI) and Mukheibeh (MU), and the location of the Hammat Gader (HG) springs. The locations of seismic lines (DS-3543 and GP-3662; Fig. 2) and of a geotectonic cross section along the LYG (A–A'; Fig. 3) are indicated as well.

The Lower Yarmouk Gorge (LYG) is a prominent geomorphological feature located in the triangle between Israel, Syria, and Jordan, east of the Dead Sea Transform (DST) between the Sheikh Ali Fault to the north (in the central Golan Heights, henceforth called GH) and the Zarka Fault to the south, in northwestern Jordan (Fig. 1). Structurally, the LYG is located along the southern extremity of the Golan syncline bounded by Mt. Hermon in the north and by the Ajloun Dome in the south (Mor, 1986). The hinge line of the Golan syncline is considered to be located along the Sheikh Ali tensional fault zone, a few kilometers north of the LYG (Meiler, 2011; Shulman et al., 2004).

The LYG drains the natural flow of the 6833 km2 large surface catchment of the Yarmouk River basin (Fig. 1) and flows into the Jordan River south of the Sea of Galilee. Additionally, the LYG receives groundwater flows from the Ajloun Dome (El-Naser, 1991), from Mt. Hermon through deep aquiferous formations in the subsurface of the Golan Heights (Gvirtzman et al., 1997), and from the Hauran Plateau (Siebert et al., 2014), situated to the NE.

Figure 2Seismic lines at southern Golan. (a) Seismic line DS-3545 showing a flower structure reinterpreted following Meiler (2011). Vertical scale in meters. (b) Seismic line GP-3662 showing the thrust fault reinterpreted following Bruner and Dekel (1989). Vertical scale in two-way time.

Figure 3Geotectonic cross section along the Lower Yarmouk Gorge, showing the interpretation of seismic line data and borehole information, which result in the given faults and depth position of geological formations.

Through the springs of Hammat Gader and Mukheibeh and the wells of Mukheibeh (henceforth abbreviated as MU) and Meizar (henceforth abbreviated as MI), all located along the LYG, groundwater emerges with different hydraulic heads, chemical compositions, and temperatures (Siebert et al., 2014). Most studies of the springs emerging on both sides of the gorge were carried out between 1970 and the 1990s, considered mainly their geochemical and hydrological characteristics, and concluded that their outflow is a mixture of shallow, cold, and fresh groundwater with a deep water body of higher salinity and temperature (Arad and Bein, 1986; Arad et al., 1986; Baijjali et al., 1997; Eckstein and Simmonsi, 1977; El-Naser, 1991; Mazor et al., 1973, 1980). However, the mechanisms and pathways allowing the deep hot brine to ascend remained obscure. Deep-reaching faults below the area of the LYG would be promising structural features controlling such hydrogeological systems.

Recent studies attempted to explain the mechanisms responsible for the outflow of hot springs and high-pressure groundwater identified in wells located along the LYG. Roded et al. (2013) applied a conceptual structural model of the Golan syncline to numerically simulate regional groundwater flow and heat transfer, introducing a 5 km wide zone of enhanced vertical permeability below the LYG. Based on geochemical evidence, Siebert et al. (2014) suggested that the groundwaters in the LYG are fed by water originating from (a) Mt. Hermon, (b) northern Jordan, and (c) the Hauran Plateau in Syria, and do not mix arbitrarily, which suggests the existence of a zone of hydraulic anisotropy along the LYG. In that weakness zone, brines likely heated by buried dikes may ascend and feed the Hammat Gader springs and some boreholes in the central LYG.

Numerical modeling of a 2-D section crossing the LYG inferred the existence of a fault in the gorge and revealed a complex groundwater flow pattern in the study area (Magri et al., 2015). This was followed by a 3-D model of a hypothetical fault tracing along the LYG and suggested a mechanism of heat-driven convection cells (Magri et al., 2016). Consequently, the anomalous heat flow was studied by a method of inverse problem, suggesting that the thermal constraints of the system require one of two scenarios, (a) relatively permeable and continuous faults, cutting through the entire geological section, reaching at depth the Triassic beds, or (b) local fractures interconnected by a highly permeable Cretaceous aquifer (Goretzki et al., 2016). Tzoufka et al. (2018) were able to substantiate the existence of a zone of hydraulic anisotropy along the LYG, possessing high hydraulic conductivity along its course, while S  N oriented groundwater migration across the plane is impeded. However, the above-mentioned studies are based on simplified assumptions of fault locations and orientation.

2 Stratigraphy

The exposed stratigraphy of the Golan Heights reveals mostly Pliocene to Quaternary basalts (Mor, 1986; Heimann et al., 1996; Dafny et al., 2003). In the north, the entire sequence from Quaternary basalts to Jurassic limestone is exposed close to Mt. Hermon (Hirsch, 1996; Picard and Hirsch, 1987). In the central part of the Golan Heights and along its southern margins Eocene to Miocene sediments are exposed in the wadis flowing towards the DST (Michelson, 1979). Middle Cretaceous rocks crop out in the Ajloun Dome, revealing a Campanian formation which does not exist either west of the DST (Flexer, 1964) or in the northern part of the GH. This formation is known by the Jordanian nomenclature as the Amman Silicified Limestone (ASL) or – by its hydrological term – as the B2 Aquifer (Andrews, 1992; El-Naser, 1991).

By its distinctive lithology of silicified limestone and chert, the Campanian ASL formation is easily identified in outcrops and boreholes. It overlies the Santonian (B1) limestone and is overlain by a thick Maastrichtian (B3) layer of marl. The transition to the underlying Santonian (B1) is somehow difficult to identify because the B1 unit is a thin layer (30–50 m) of limestone without cherts. However, the occurrence of dolomite or of dolomitic limestone clearly defines the Turonian (A7). Contrarily to the hydrological characteristics of the Senonian aquiclude in central and northern Israel, the ASL (B2) in northern Jordan makes up the upper part of the most essential aquifer in Jordan. The MI-1, MI-2, and MI-3 boreholes drilled in southern Golan up to a distance of 6 km north of the LYG reveal a lithostratigraphic sequence similar to the Ajloun, confirming that the units continue across the gorge.

South of the Ajloun Dome, in Wadi Zarka (Fig. 1), the exposed stratigraphy reveals the lithostratigraphic sequence down to Jurassic beds. The full section down to the Precambrian basement occurs in several deep boreholes drilled in northern Jordan and within the eastern escarpment of the lower Jordan Valley (Abu-Saad and Andrews, 1993). In the southern GH, MI-1 is the northernmost borehole drilled to Turonian beds providing a complete section of the overlaying lithology.

3 Tectonics

The Kinnarot basin is a link in the chain of pull-apart basins scattered along the DST formed by the left lateral movement along the Sinai–Arabian Plate boundary that started during the Early to Middle Miocene (Garfunkel, 1981; Gvirtzman and Steinberg, 2012). The transform itself was located on the eastern side of the basin by seismic interpretation showing the compressional structure of the Tel-Qatzir elevated block (Inbar, 2012). To the north, Mt. Hermon manifests a shift in the en-echelon arrangement from left- to right-stepping resulting in a restraining geometry (Weinberger et al., 2009) which causes uplifting of the Lebanon and Anti-Lebanon mountains (Beydoun, 1977) as well as the deepening of the Golan syncline. The LYG, the southern GH, and the northern Ajloun, located along the eastern rim of the Kinnarot basin, are considered to be subject to the regional forces applied by this active plate boundary.

Models of pull-apart basins with different scales and from various places around the world (e.g., Corti and Dooley, 2015; Sagy and Hamiel, 2017; Smit et al., 2008, 2010; Sugan et al., 2014) address mostly the shape of their boundaries and internal processes. However, these models predict fault evolution at the rims as well. Such a scenario as the collision of the Kinnarot basin's eastern rim with Mt. Hermon has not been modeled to the best of our knowledge. The conspicuous Sheikh Ali fault in the northeastern corner of the basin (Fig. 1) is a clear indication of this different tectonic setting. Therefore, the Kinnarot basin seems to provide a unique opportunity to explore the effect of such a setting on fault evolution at pull-apart basin rims.

Along the eastern side of the DST, strike–slip faults are known to branch out and penetrate inland, directed NE–SW and E–W, across its rims (Fig. 1; Andrews, 1992; Shulman, 2004). Though the LYG is a prominent morphological feature east of the DST, its structure was never explored by continuous land geophysics (e.g., seismic, electric) due to its function as an international border. Consequently, attempting to bridge the gap in data resulted only in theoretical and qualitative models.

Considering the extreme thickness difference between Jurassic sections in Mt. Hermon and in the Ajloun (2000 and 400 m, respectively) and by comparing seismic lines from the southern GH and from the Coastal Plain of Israel, Shulman et al. (2004) suggested the existence of a major fault tracing along the LYG. However, based on interpolation between the results of the interpretation of seismic data in the southern GH and deep borehole data in Jordan, it was suggested that the thickness of the Jurassic rock sequence increases gradually from the south northwards, and therefore a major fault in the gorge area should not be inferred (Meiler, 2011). The present study shows that although there is no evidence for large vertical displacements at the surface, seismic and borehole data indicate strike–slip faults crossing the LYG, forming fault blocks. Therefore, such faults should be taken into consideration when discussing groundwater hydrology.

Table 1Campanian chert-bearing formation (B2) – thickness and elevation in wells adjacent to the LYG.

* Source of information: Margane and Hobler (1994) available via DAISY (2017). ${}^{**}$ Source of information: well log. + Partial thickness, not fully penetrated unit.

4 Data and methods

The current study relies on reviewing, compiling, and evaluating available geological and geophysical data from southern Golan and northern Ajloun. Most of the seismic data were collected at the Golan Heights during the 1980s (Shulman and May, 1989), published later as time-migrated data analyses (Shulman et al., 2004), and reprocessed by Meiler (2011) using PSDM (pre-stack depth migration) to analyze regional aspects. In the current study, we aimed to refine Meiler's faulting interpretation of PSDM line DS-3545 (Figs. 1 and 2) to shift the precision for a less broad perspective (i.e., Meiler neglected the westernmost fault in the regional context, which might be relevant at a local scale). Additionally, seismic data from southern Golan (Mevo Hamma; Fig. 1) (Shulman et al., 2004; Bruner and Dekel, 1989; Meiler, 2011) were reinterpreted to considerably improve the structural information in the vicinity of the Lower Yarmouk Gorge. Based on these data, in combination with lithological data from boreholes (Table 1), a new geological cross section was constructed along the LYG, running SW–NE (Fig. 3). Faults from surface mapping (Michelson, 1972), reinterpreted seismic lines, and geological profiles (Sahawneh, 2011) were considered in order to generate an aerial view of faults and of fault-block patterns in the study area based on the tectonic concept of pull-apart basin rims. Such models predict evolution of different types of faults at the margins of the main basin mainly according to its maturity, size, and symmetry (obliqueness) (Wu et al., 2009; Rahe et al., 1998; Sugan et al., 2014; Smit et al., 2008, 2010). Finally, the previously mapped faults were traced together with those identified during the present study. It is important to note that the fault lines drawn on aerial view maps are a representation of an actual near-surface fault zone which converges to a single deep root as illustrated by seismic interpretation and the geological cross section.

5 Results and discussion

## 5.1 Reinterpretation of seismic data

The DS-3545 ENE–WSW trending and 17 km long seismic line (Fig. 2a) runs parallel to the LYG, about 6 km north of it. Reinterpretation of the PSDM conducted by Meiler (2011) shows an additional flower-structure fault at its southwestern part, indicating a set of strike–slip faults crossing that line. Although little to no horizontal displacement is visible in that seismic section, folds between the flower-structure fault branches clearly indicate lateral displacement. Moreover, the deep roots of the traced strike–slip faults suggest that it is related to significant regional tectonics. Its proximity to the DST combined with its effect on shallow lithology may advocate its connection to the DST tectonics.

Seismic line GP-3662 (Fig. 2b) trends NE–SW and is 4.5 km long. Its southwestern end is located near MI-2. The time-migrated line described by Bruner and Dekel (1989) was re-evaluated considering the results from the reinterpretation of seismic line DS-3545 and MI-1 and MI-2 stratigraphy. A near-top Turonian marker was adjusted to fit its depth in the nearby Meizar wells and DS-3545 seismic line by depth-to-time conversion using rms velocities of the original GP-3662 time-migrated section. The near-top Turonian reflector (Fig. 2b, green markers) is interrupted by a thrust fault indicating compressional stress along the line. The different dip angles of shallower reflectors (Fig. 2b, yellow markers) at the southwestern side of the line as well as folded reflectors in the upper part of the section support that idea.

Surface and shallow geological data in the southern GH and northern Ajloun indicate possible surface and shallow fault patterns (Sahawneh, 2011; El-Naser, 1991; Michelson, 1979). All these data suggest the existence in the study area of short faults of limited vertical displacement. However, deep seismic data collected in the southern GH (Shulman et al., 2004; Meiler, 2011) revealed that normal faults with minor to no vertical displacement at the surface and shallow subsurface may indicate deep strike–slip faults (Fig. 2a).

## 5.2 Geological section along the LYG

Borehole information from the Meizar and Mukheibeh wells, drilled along the LYG close to the area of the Hammat Gader springs (Figs. 1 and 4a), provides a unique opportunity to explore the complex faulting pattern along major parts of the gorge. Based on the lithological interpretation of well logs (IHS, Margane and Hobler, 1994) that were drilled along the LYG, a geological cross section was constructed and validated by the lithological description of all other well sections in the area (Fig. 4b). At least the uppermost part of the B2–A7 aquifer, the Campanian B2 horizon, is disclosed in all wells of the study area. Hence, borehole data about its elevation and thickness provide good coverage in the study area and are the key to the current work (Table 1).

Figure 4Results of the current study in a map local frame similar to the one used in Fig. 1. The proposed faults are in the area of the Lower Yarmouk Gorge, east of the Kinnarot pull-apart basin, Mevo Hamma Fault (MHF), and Lower Yarmouk Fault (LYF). It is suggested that those faults result from the settings at the eastern rim of the DST at the Golan Heights where the en-echelon changes from pull-apart basins to a push-up ridge at Mt. Hermon (Fig. 1).

The base of B2 displays large elevation changes over small distances, e.g., in borehole MU-6 located between MU-JRV1 and MU-8 (Fig. 4a), at distances of 1.2 and 0.6 km, respectively. Although the B2 horizon appears at similar elevations in boreholes MU-JRV1 (−408 m a.m.s.l.) and MU-8 (−403 m a.m.s.l.), it was encountered considerably deeper in well MU-6 (−480 m a.m.s.l.). Even more distinct, between boreholes MU-7A, MU-4, and MU-2, the elevation of base B2 varies over a horizontal distance of 40 m to an extent of 130 m. These differences result from faulting (Fig. 3).

Considering places where the B2 formation was fully penetrated (Table 1) only, the average thickness of the chert-bearing B2 layer is 189±14 m. The MU-8 well contains an exceptionally thick Senonian sequence exceeding 290 m. Though the Senonian sequence in the region is well known for its thickness variations (Rosenthal et al., 2000a, b; Rosenthal, 1972), such variations over short distances are exceptional and require a different explanation. It is therefore suggested that the unit is either strongly tilted or that thickness was “doubled” in the drilling due to crossing a thrust fault.

In the study area, the thickness of Turonian beds (A7) is fairly uniform and ranges between 300 and 350 m. However, in borehole MU-JRV1, a section of about 700 m consisting of two repeating sequences of marly limestone and dolomitic limestone is regarded as a Turonian unit “doubled” by thrust faulting (Fig. 3).

## 5.3 Tracking fault paths

Following the seismic interpretation, it is suggested that faults detected on the surface or at shallow depths in the study area are likely to indicate deep-seated strike–slip faults. Based on that hypothesis, faults shown on geological cross sections in the western Ajloun area (Sahawneh, 2011) were considered the southern extension of the Lower Yarmouk Fault (LYF). At the southernmost end, the fault branches out from the DST and continues to the NE. The fault crosses the LYG west of the Hammat Gader–Meizar–Mukheibeh area (Fig. 1). North of the gorge, the fault turns further eastwards following the outcropping fault lineaments mapped by Michelson (1979). Another possibility is a northward continuation of the fault joining with the Nov Fault Zone (NFZ) (Shulman et al., 2004).

Another SE–NW strike–slip fault was drawn according to previous interpretation and current reinterpretation of seismic data. Although parts of this fault were previously mapped as normal faults, our revised interpretation suggests that it is most likely a strike–slip fault (Mevo-Hama Fault, MHF). In between these two newly suggested strike–slip faults (MHF and LYF), there is a thrust fault which clearly stands out in two seismic lines, GP-3661 and GP-3662.

6 Summary and conclusions

A comprehensive reinterpretation and compilation of available geological and geophysical data across the LYG at a local scale are presented. Former seismic analysis of the deep structure in the Golan Heights (Meiler, 2011; Shulman et al., 2004; Shulman and May, 1989) used deep borehole data from Jordan and Syria to study regional-scale structures. However, in these studies, all rock units between the Turonian and Pliocene basalts were presented as one undivided unit, and some of the faults at the southernmost DS-3545 line were neglected. The shallow seismic survey conducted in the Mevo Hamma area (Bruner and Dekel, 1989) was originally interpreted simultaneously with preliminary results of Shulman and May (1989) and with no access to data from northern Ajloun. Since then, the Mevo Hamma data have not been used by any other study done in the area.

As the Mukheibeh well field was further developed in later years, more data of the Upper Cretaceous have been made available alongside hydrological observations. The rejuvenation of hydrological studies of the LYG started with conceptual models (Roded et al., 2013; Siebert et al., 2014), which speculated about potential reasons for local hydraulic anisotropies in order to explain the complex hydrological system manifested by groundwater pressure, temperature, and chemical variations. This was followed by a series of numerical 2-D models (Goretzki et al., 2016; Magri et al., 2015; Tzoufka et al., 2018), which used a combination of previous works to construct their structural model. The results achieved by those studies simulate groundwater flow and heat transfer over a fault in the geological sequence of the area. However, the transition of a 2-D model into semi-3-D (Magri et al., 2016) emphasized the problem of numerically simulating groundwater dynamic processes using a speculated structural model and the need for a unified structural solution in a complex setting like the LYG.

While tectonic models of pull-apart basins have been developed for different scales and various places around the world (e.g., Corti and Dooley, 2015; Sagy and Hamiel, 2017; Sugan et al., 2014), including models for the Dead Sea basin itself (Smit et al., 2008, 2010), they mostly address the shape of basin boundaries and internal processes. The unique situation at the eastern rim of the Kinnarot pull-apart basin includes the collision of the Golan Heights at Mt. Hermon. To the best of our knowledge, such a scenario has not yet been modeled. Hence, the suggested fault system does not entirely comply with the findings of previous research on the surrounding tectonics of pull-apart basins performed in other tectonic settings. It is suggested that this mismatch may result from the settings at the eastern rim of the basin where the en-echelon changes from a chain of pull-apart basins to a push-up ridge.

Using available geological and geophysical data (Fig. 1), from the southern Golan and northern Ajloun, a new fault pattern has been delineated across the Lower Yarmouk Gorge. It includes a series composed of strike–slip and thrust faults, which may be associated with the regional Dead Sea transform system and with the local Kinnarot pull-apart basin. It seems that these compressional and tensional structures have been developed to form a series of fault blocks, causing a non-uniform spatial hydraulic connection between them.

The main motivation for the work presented here was to provide fundamental structural and geological information, which may have hydrological implications. It is common that structural features do affect groundwater movement. To develop management strategies for a reservoir, which is particularly important in areas of water scarcity and transboundary water resources, transient 3-D flow simulations of the resource are the most appropriate solution to understand reservoir behavior. However, they must be based on realistic geometry, including structural features. Though formulated for the LYG, the study is intended to show the importance of such studies in providing the necessary database in geologically stressed areas without sufficient data.

Data availability
Data availability.

The current study relies on reviewing, compiling, and evaluating available geological and geophysical data from southern Golan and northern Ajloun. The data have been collected from the cited sources (papers, reports, publicly open data bases).

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Special issue statement
Special issue statement.

Acknowledgements
Acknowledgements.

This research was supported by the German Science Foundation (DFG, grant MA4450/2) within the special program to support peaceful development in the Middle East. We greatly acknowledge the DFG for their continuous support. We thank Eyal Shalev and Miki Meiler for sharing seismic data of the southern Golan Heights. We thank Armin Margane and two anonymous reviewers for their very helpful comments, which significantly improved the quality of the paper.

Edited by: Efrat Morin
Reviewed by: Armin Margane and two anonymous referees

References

Abu-Saad, L. and Andrews, I. J.: A database of stratigraphic information from deep boreholes in Jordan, Subsurface Geology Bulletin, edited by: NRA, A., Natural Resources Authority (NRA), Amman, 6, 182 pp., 1993.

Andrews, I. J.: Cretaceous and Paleogene lithostratigraphy in the subsurface of Jordan, Subsurface Geology Bulletin, edited by: NRA, A., Natural Resources Authority (NRA), Amman, Jordan, 5, 60 pp., 1992.

Arad, A. and Bein, A.: Saline versus fresh water contribution to the thermal waters of the northern Jordan Rift Valley, Israel, J. Hydrol., 83, 49–66, 1986.

Arad, A., Kafri, U., Halicz, L., and Brenner, I.: Genetic identification of the saline origins of groundwaters in Israel by means of minor elements, Chem. Geol., 54, 251–270, https://doi.org/10.1016/0009-2541(86)90140-3, 1986.

Baijjali, W., Clark, I. D., and Fritz, P.: The artesian thermal groundwaters of northern Jordan: insights into their recharge history and age, J. Hydrol., 192, 355–382, https://doi.org/10.1016/S0022-1694(96)03082-X, 1997.

Beydoun, Z.: Petroleum prospects of Lebanon: reevaluation, AAPG Bull., 61, 43–64, 1977.

Bruner, I. and Dekel, I.: High resolution seismic reflection survey in the southern Golan Heights, The Institute for Petroleum Research and Geophysics, Holon, Israel, Rep. 120/1939/88, 1989.

Corti, G. and Dooley, T. P.: Lithospheric-scale centrifuge models of pull-apart basins, Tectonophysics, 664, 154–163, https://doi.org/10.1016/j.tecto.2015.09.004, 2015.

Dafny, E., Gvirtzman, H., Burg, A., and Fleischer, L.: The hydrogeology of the Golan basalt aquifer, Israel, Israel J. Earth. Sci., 52, 139–153, 2003.

Eckstein, Y. and Simmonsi, G.: Measurement and interpretation of terrestrial heat flow in Israel, Geothermics, 6, 117–142, https://doi.org/10.1016/0375-6505(77)90023-2, 1977.

El-Naser, H.: Groundwater resources of deep aquifer system in NW-Jordan: hydrogeological and hydrogeochemical quasi 3-dimensional modelling, Uni. Würzburg, Würzburg, 144 pp., 1991.

Flexer, A.: The paleogeography of the Senonian and Maestrichtian in northern Israel, PhD, Hebrew University, Jerusalem, 157 pp., 1964.

Garfunkel, Z.: Internal structure of the Dead Sea leaky transform (rift) in relation to plate kinematics, Tectonophysics, 80, 81–108, 1981.

Goretzki, N., Inbar, N., Kühn, M., Möller, P., Rosenthal, E., Schneider, M., Siebert, C., Raggad, M., and Magri, F.: Inverse Problem to Constrain Hydraulic and Thermal Parameters Inducing Anomalous Heat Flow in the Lower Yarmouk Gorge, Enrgy. Proced., 97, 419–426, https://doi.org/10.1016/j.egypro.2016.10.038, 2016.

Gvirtzman, H., Garven, G., and Gvirtzman, G.: Thermal anomalies associated with forced and free ground-water convection in the Dead Sea rift valley, Geol. Soc. Am Bull., 109, 1167–1176, https://doi.org/10.1130/0016-7606(1997)109<1167:taawfa>2.3.co;2, 1997.

Gvirtzman, Z. and Steinberg, J.: Inland jump of the Arabian northwest plate boundary from the Levantcontinental margin to the Dead Sea Transform, Tectonics, 31, TC4003, https://doi.org/10.1029/2011TC002994, 2012.

Heimann, A., Steinitz, G., Mor, D., and Shaliv, G.: The Cover Basalt Formation, its age and its regional and tectonic setting: implication from K-Ar and Ar40 Ar39 geochronology, Israel J. Earth Sci., 45, 55–72, 1996.

Hirsch, F.: Geology of the southeastern slopes of Mount Hermon, Current Research, Geological Survey of Israel, 10, 22–27, 1996.

Inbar, N.: The evaporitic subsurface body of Kinnarot Basin – Stratigraphy, Structure, Geohydrology, Tel Aviv University, Tel Aviv, 131 pp., 2012.

Magri, F., Inbar, N., Siebert, C., Rosenthal, E., Guttman, J., and Möller, P.: Transient simulations of large-scale hydrogeological processes causing temperature and salinity anomalies in the Tiberias Basin, J. Hydrol., 520, 342–355, https://doi.org/10.1016/j.jhydrol.2014.11.055, 2015.

Magri, F., Möller, S., Inbar, N., Möller, P., Raggad, M., Rödiger, T., Rosenthal, E., and Siebert, C.: 2-D and 3-D coexisting modes of thermal convection in fractured hydrothermal systems – Implications for transboundary flow in the Lower Yarmouk Gorge, Mar. Petrol. Geol., 78, 750–758, https://doi.org/10.1016/j.marpetgeo.2016.10.002, 2016.

Margane, A. and Hobler, M.: Groundwater Resources of Northern Jordan, Vol. 3: Structural Features of the Main Hydrogeological Units in Northern Jordan – Technical Cooperation Project “Advisory Services to the Water Authority of Jordan”, BGR & WAJ, Amman, BGR archive no. 118702:1-3, 30 app., 38 ann., 57 pp., 1994.

Mazor, E., Kaufman, A., and Carmi, I.: Hammat Gader (Israel): Geochemistry of a mixed thermal spring complex, J. Hydrol., 18, 289–303, https://doi.org/10.1016/0022-1694(73)90053-X, 1973.

Mazor, E., Levitte, D., Truesdell, A. H., Healy, J., and Nissenbaum, A.: Mixing models and ionic geothermometers applied to warm (up to 60 C) springs: Jordan Rift Valley, Israel, J. Hydrol., 45, 1–19, https://doi.org/10.1016/0022-1694(80)90002-5, 1980.

Meiler, M.: The Deep Geological Structure of the Golan Heights and the Evolution of the Adjacent Dead Sea Fault System, Tel Aviv University, Tel Aviv, 2011.

Michelson, H.: The hydrogeology of the southern Golan Heights, Tahal, Tel Aviv, HR/72/037, Tahal, Tel Aviv, Israel, 89, 1972.

Michelson, H.: The Geology and Paleogeography of the Golan Heights, PhD, Tel Aviv University, Tel Aviv, 163 pp., 1979.

Mor, D.: The volcanism of the Golan Heights, Geological Survey of Israel, Jerusalem, Israel, GSI/5/86, 159, 1986.

Picard, L. and Hirsch, F.: The Jurassic stratigraphy in Israel and the adjacent countries, The Israel Academy of Sciences and Humanities, Jerusalem, Israel, ISBN: 965-208-076-4, 106 pp., 1987.

Rahe, B., Ferrill, D. A., and Morris, A. P.: Physical analog modeling of pull-apart basin evolution, Tectonophysics, 285, 21–40, https://doi.org/10.1016/S0040-1951(97)00193-5, 1998.

Roded, R., Shalev, E., and Katoshevski, D.: Basal heat-flow and hydrothermal regime at the Golan–Ajloun hydrological basins, J. Hydrol., 476, 200–211, https://doi.org/10.1016/j.jhydrol.2012.10.035, 2013.

Rosenthal, E.: Hydrogeology and Hydrogeochemistry of the Bet Shean and Harod Valleys, Geological Survey of Israel, Jerusalem, Israel, Rep. H/2/1972, 1972.

Rosenthal, E., Weinberger, G., Almogi-Labin, A., and Flexer, A.: Late Cretaceous-Early Tertiary Development of Depositional Basins in Samaria as a Reflection of Eastern Mediterranean Tectonic Evolution, AAPG Bull., 84, 997–1014, https://doi.org/10.1306/A9673B94-1738-11D7-8645000102C1865D, 2000a.

Rosenthal, E., Weinberger, G., Almogi-Labin, A., and Flexer, A.: Late Cretaceous – Early Tertiary paleostructures in Samaria, Israel, Geological Survey Current Research, Geological Survey of Israel, Jerusalem, Israel, 12, 159–167, 2000b.

Sagy, A. and Hamiel, Y.: Three-dimensional structure and evolution of an asymmetric pull-apart basin, Int. J. Earth Sci., 106, 1619–1630, https://doi.org/10.1007/s00531-016-1374-4, 2017.

Sahawneh, J.: Structural Control of Hydrology, Hydrogeology and Hydrochemistry along the Eastern Escarpment of the Jordan Rift Valley, JORDAN, PhD, Karlsruher Institut für Technologie (KIT), 284 pp., 2011.

Shulman, H. and May, P. R.: Seismic survey in the Golan Heights and its results, Isr. Geol. Soc. An. Meet., Ramot, Israel Geological Society, Jerusalem, Israel, 155–156, 1989.

Shulman, H., Reshef, M., and Ben-Avraham, Z.: The structure of the Golan Heights and its tectonic linkage to the Dead Sea Transform and the Palmyrides folding, Israel J. Earth Sci., 53, 225–237, 2004.

Siebert, C., Möller, P., Geyer, S., Kraushaar, S., Dulski, P., Guttman, J., Subah, A., and Rödiger, T.: Thermal waters in the Lower Yarmouk Gorge and their relation to surrounding aquifers, Chem. Erde-Geochem., 74, 425–441, https://doi.org/10.1016/j.chemer.2014.04.002, 2014.

Smit, J., Brun, J. P., Fort, X., Cloetingh, S., and Ben-Avraham, Z.: Salt tectonics in pull-apart basins with application to the Dead Sea Basin, Tectonophysics, 449, 1–16, 2008.

Smit, J., Brun, J. P., Cloetingh, S., and Ben-Avraham, Z.: The rift-like structure and asymmetry of the Dead Sea Fault, Earth Planet. Sc. Lett., 290, 74–82, 2010.

Sugan, M., Wu, J., and McClay, K.: 3-D analogue modelling of transtensional pull-apart basins: comparison with the Cinarcik basin, Sea of Marmara, Turkey, B. Geofis. Teor. Appl., 55, 699-716, https://doi.org/10.4430/bgta0129, 2014.

Tzoufka, K., Magri, F., Rödiger, T., Inbar, N., Shalev, E., Möller, P., Raggad, M., Rosenthal, E., and Siebert, C.: The effect of hydraulic anisotropies on intensely exploited groundwater resources: the numerical evaluation of a hydrothermal transboundary aquifer system in the Middle East, Hydrogeol. J., 26, 2875–2890, https://doi.org/10.1007/s10040-018-1827-x, 2018.

Weinberger, R., Gross, M. R., and Sneh, A.: Evolving deformation along a transform plate boundary: Example from the Dead Sea Fault in northern Israel, Tectonics, 28, TC5005, https://doi.org/10.1029/2008TC002316, 2009.

Wu, J. E., McClay, K., Whitehouse, P., and Dooley, T.: 4-D analogue modelling of transtensional pull-apart basins, Mar. Petrol. Geol., 26, 1608–1623, https://doi.org/10.1016/j.marpetgeo.2008.06.007, 2009.