ALBERT

Description: Transient numerical simulations of coupled fluid flow and heat transport processes were conducted, to investigate (i) the impact of hydraulic anisotropy on a complex hydrogeological system and (ii) the anomalous geothermal gradient in the Upper Cretaceous (A7/B2) aquifer of the Lower Yarmouk Gorge (LYG). The N-S directed geological profile starts on the basalt-covered Golan Heights, crosses the LYG and ends at the Jordanian Ajloun Plateau. Heated fresh groundwaters ascend within the LYG from the confined A7/B2 limestone aquifer through artesian Meizar wells and artesian Muhkeibeh well field. A hydrogeological study based on high frequency water-table measurements at the Meizar wells suggested strong impact of abstraction at the Mukheibeh wells on the hydraulic head distribution north of the LYG (Shalev et al., 2015). Contrastingly, hydrochemical investigations conducted in the area concluded that recharge areas of the A7/B2 aquifer are: (i) the foothills of Mountain Hermon, (ii) the Ajloun Plateau and (iii) the Syrian Hauran Plateau, indicating the presence of a zone of high-hydraulic anisotropy along the main LYG axis. Due to this still debatable hydraulic feature, flow along the LYG principle axis is enhanced whilst flow perpendicular to it is constrained (Siebert et al., 2014). In agreement, transient simulations based on a NW-SE profile supported further the hypothesis of a structural feature existent (Magri et al., 2015). The modeled profile cross-cuts the heterogeneous zone and hence the heterogeneity was implemented as an impermeable zone by employing the Equivalent Porous Media approach. Initial 2D models managed to successfully reproduce the natural hydraulic head and temperature distributions. In subsequent simulations, by implementing Meizar abstraction rates, results revealed that mixed convection explains the anomalous temperature gradient in the area as temperature patterns of these simulations are in accordance with a temperature-depth borehole profile. Most significantly, the existence of a zone of anisotropic hydraulic conductivity at the LYG, as suggested above, is supported. Sensitivity analysis of the heterogeneous permeability zone (ranging between 1e-7 m/s and 1e-10 m/s) simultaneously displayed hydraulic connectivity and prevented fluid flux occurrence between the northern and southern LYG flanks. Within the studied aquifer system, aquifer topography directs groundwater flow towards the LYG (N-S and S-N). Once groundwaters reach the zone of hydraulic anisotropy the direction is diverted westwards, towards the Lower Jordan Valley. A final finding pointed to the inadequate size of the A7/B2 drainage basin south of the LYG as the Mukheibeh well field requires a larger drainage area.

Description: The Lower Yarmouk Gorge (LYG) extends on the eastern margin of the Jordan Rift Valley (JRV) adjacent to the Kinneret basin which makes part in the chain of pull-apart basins along the Dead Sea Transform (DST). The LYG is bounded to the south by the Ajloun Plateau (northern Jordan) and to the north by the Golan Heights (Israel). It acts as the outflow of the Yarmouk drainage basin into the Jordan River a few kilometers south of Lake Tiberias. Although topographically the Golan Heights and the Ajloun seem to be different provinces separated by the LYG, the northern Ajloun and southern Golan form the southern flank of a major synclinal structure. Morphologically, the LYG resembles the Sheikh-Ali strike-slip fault and other NE-SW striking faults related to the DST. However, the existence of faulting along that line is debatable. As the gorge serves also as state boundary, no seismic lines cross the LYG. Quaternary landslides, mostly on the southern flank of the LYG cover possible surface evidence of faulting although Quaternary basalts located at the gorge path may indicate possible vertical pathways. Moreover, hydrological studies (Siebert et al., 2014, Goretzki et al., 2016) show that permeability anisotropy along the LYG line allows heated groundwater to emerge along the gorge with temperatures rising up to 60C. The presented study uses well data from northern Jordan and southern Golan Heights as well as seismic data from the southern Golan Heights to bridge over the information gap. Based on the available information we present our hypothesis on the geology of the LYG aiming to contribute to the discussion regarding active tectonics at the Lower Yarmouk Gorge. Our preliminary results show that a major fault along the gorge path is dubious however, perpendicular faulting is more likely to occur. We expect the results of current research to contribute to the understanding of the local complex geohydrological system. Furthermore, results are expected to deepen our understanding of pull-apart basin tectonics effect on the evolution of marginal structures.

Description: Volcanic rocks covering vast areas in central north Jordan and in central and southern Syria erupted during 6 different phases starting in Miocene and continuing - with major interruptions - into the Holocene. The petrological composition of the different flows of the Harrat ash Shaam Basalt complex is quite homogeneous with the major minerals: Plagioclase, K-feldspar, clinopyroxene, amphibole, biotite, olivine, magnetite, limonite, goethite, pyrite and chalcopyrite. The oldest basalts cover Cretaceous and Paleogene sediments, which at that time formed the land surface of drainage basins. The basaltic aquifer contains groundwater with a wide range of salinities. They represent a continuous sequence of increasingly mineralized groundwater originating from precipitation over Jebel Druz flowing radially into all directions, in coincidence with the topographic slopes. Along the flow-paths halite and gypsum are dissolved. Ca2+ not only depends on gypsum dissolution but also increases proportionally to Mg. This may suggest that the combination of Ca2+, Mg2+ and sulfate is a saline endmember fluid originating from the underlying carbonate formations of the basalt. Mixing with recharge water could explain the chemical composition of the various types of water. The signature of dissolved gypsum and halite indicates dissolution of evaporites that might have formed by evaporation either before the basalt covered the area or due to the hot basalts heating up the underlying carbonates and their enclosed fluids. Evaporation of water precipitated evaporites. Ca and Mg halides are hygroscopic, thus they are only present in solution. Such saline water, however, has not affected the low saline groundwater because their increase in Ca depends neither on the increase of Mg2+ nor of SO2�� 4 . This leaves the formation of clay minerals as the probably sink for Na. Inverse modelling applying PHREEQC with phreeq.dat database reveals that the mineralization of groundwater increases due to dissolution of increasing amounts of halite and gypsum which are mass-wise, the most important reactants. Concurrently, albite increasingly precipitates. Montmorillonite, gibbsite and calcite form, whereas kaolinite is consumed. Sulfides are oxidized. D and 18O of well and springs fit an evaporation line rooted on the Ajloun MWL. Hydrochemically, there are two sources of salts: Mixing with a saline endmember brine and/or dissolution of evaporites. Near Jebel Druz, dissolution of evaporites dominates, whereas mixing with a saline endmember and formation of clay minerals occur at greater distances.

Description: The Tiberias Basin (TB) is located within the Jordan–Dead Sea Transform and is bordered to the west by the Lower Galilee (Israel), where Pliocene basalts cover an area of 35 km2. Hydrochemical analyses highlight that two types of brines exist around Lake Tiberias (LT) (Mandel 1965; Möller et al. 2009): (1) Along the eastern side of LT, brine is characterized by Mg/Ca〉1, which resulted from evaporation of seawater during the Late Miocene, whereas (2) along the western side of the lake, brine is characterized by Mg/Ca〈1, possibly formed out of the Mg-rich brine by dolomitization of limestones (Möller et al. 2012). Dolomitization of limestones occurs at temperatures of at least 100 C. We suppose that basalts which erupted through numerous fissures, forming nowadays sills within the Cretaceous and Eocene limestones, sufficiently heated the formations, which build up the Lower Galilee, west of LT. As a result, the Cenomanian Formations, where the original brine is mostly buried, were only sufficiently heated in the eastern Galilee. In this study, we try to estimate to which extent and through which mechanisms fissure eruptions have induced heated brine to flow within the limestone aquifers. 2D simulations of coupled heat transport and fluid flow show that different aspects control the heat transport in the limestone aquifer. Preliminary results indicate that conductive heat transport generates sharp temperature fronts that extent 30 meters after 5 years of continuous magmatic intrusion from the fissures, ruling out heat conduction as a major mechanism for dolomitization. By contrast, convective cells in the Turonian and Cenomanian aquiferous formations have the potential to develop at different scales that depend on (a) hydraulic conductivity and porosity of the aquifer, (b) the orientation of the regional flow and (c) the topography of the aquifer. As a result of convective flow, brines surrounding a single fissure intrusion are heated more than 100 C up to a width of 2 km in both directions. The thermally induced flow velocities are in the range of 2 m/year. In simulations with multiple intrusions, the thermal plume can stretch over 5.5 km. The simulations indicate that magmatic induced advective/convective heating may have generated temperature conditions favorable for dolomitization, which in turn may explain the existence of two different brines that are found around the LT.

Description: Lake Kinneret (Lake Tiberias or Sea of Galilee) is the most important freshwater reservoir in the Northern Jordan Valley. Simulations that couple fluid flow, heat and mass transport are built to understand the mechanisms responsible for the salinization of this important resource. Here the effects of permeability distribution on 2D and 3D convective patterns are compared. 2D simulations indicate that thermal brine in Haon and some springs in the Yamourk Gorge (YG) are the result of mixed convection, i.e. the interaction between the regional flow from the bordering heights and thermally-driven flow (Magri et al., 2014). Calibration of the calculated temperature profiles suggests that the faults in Haon and the YG provides paths for ascending hot waters, whereas the fault in the Golan recirculates water between 1 and 2 km depths. At higher depths, faults induce 2D layered convection in the surrounding units. The 2D assumption for a faulted basin can oversimplify the system, and the conclusions might not be fully correct. The 3D results also point to mixed convection as the main mechanism for the thermal anomalies. However, in 3D the convective structures are more complex allowing for longer flow paths and residence times. In the fault planes, hydrothermal convection develops in a finger regime enhancing inflow and outflow of heat in the system. Hot springs can form locally at the surface along the fault trace. By contrast, the layered cells extending from the faults into the surrounding sediments are preserved and are similar to those simulated in 2D. The results are consistent with the theory from Zhao et al. (2003), which predicts that 2D and 3D patterns have the same probability to develop given the permeability and temperature ranges encountered in geothermal fields. The 3D approach has to be preferred to the 2D in order to capture all patterns of convective flow, particularly in the case of planar high permeability regions such as faults.

Description: Numerical investigations of 2D and 3D modes of large-scale convection in faulted aquifers are presented with the aim to infer possible transport mechanisms supporting the ormation of thermal springs through fractured aquicludes. The transient finite elements models are based on dealized structural features that can characterize many hydrothermal systems. The sensitivity analysis of the fault permeability showed that faults cross-cutting the main regional flow direction allow groundwater to be driven laterally by convective forces within the fault planes. Therein thermal waters can either discharge along the fault traces or exit the fault through adjacent permeable aquifers. In the latter case, the resulting flow is helicoidally and transient. The location and the spacing between discharge areas can migrate with time, is not strictly constrained to the damage zones and reflects the wavelength of the multicellular regime in the fault zone. An illustrative example based on simplified structural data of the Lower Yarmouk Gorge (LYG) is presented. The numerical calculations indicate that crossing flow paths result from the coexistence of fault convection, developing for example along NE-SW oriented faults within the Gorge, and additional flow fields. The latter are induced either by topography NeS gradients, e.g. perpendicular to the major axe of the Gorge, or by local thermal convection in permeable aquifers below the Eocene aquiclude. Sensitivity analysis of fault hydraulic conductivity (K) and the analytical solutions based on viscousdependent Rayleigh theory show that K values between 2.3e_7 m/s and 9.3ee 7 m/s (i.e. 7 m/yr and 30 m/yr, respectively) favor coexisting transport processes. The uprising thermal plumes spread over several hundred meters forming clusters of springs, in agreement with observation, and which temperature fall within the measured ranges, i.e. 20 _C_60 _C. To some extent the models also reproduced the transient behavior of the spring temperature. Owing to the idealized nature of the presented models, the numerical results and the associated analytical solution can be applied to study the onset of thermal convection and resulting flow patterns of any fractured hydrothermal basin.

Description: The Lower Yarmouk Gorge (LYG) lies on the eastern margin of the lower Jordan Rift Valley (JRV), bounded to the south by the Ajlun and to the north by the Golan Heights. It allows the outflow of the Yarmouk drainage basin and flow into the Jordan River, a few kilometers south of Lake Tiberias. The main aquifer system of the LYG is built mostly of Cretaceous sandstones and carbonates confined by Maastrichtian aquiclude. Fissures allow hydraulic connections between the major water-bearing formations from Quaternary to Upper Cretaceous age. It is supposed that the gorge acts as the mixing zone of two crossing flow pathways: N-S from the Hermon Mountains and from the Ajlun Dome, and E-W from Jebel al Arab Mountain in Syria (also known as Huran Plateau or Yarmouk drainage basin). As a result, several springs can be found within the gorge. These are characterized by widespread temperatures (20 – 60 C) which indicate that, beside the complex regional flow, also ascending thermal waters control the hydrologic behavior of the LYG. Previous simulations based on a conceptual simplified 3D model (Magri et al., 2016) showed that crossing flow paths result from the coexistence of convection, that can develop for example along NE-SW oriented faults within the gorge or in permeable aquifers below Maastrichtian aquiclude, and additional flow fields that are induced by the N-S topographic gradients. Here we present the first 3D hydrogeological model of the entire LYG that includes structural features based on actual logs and interpreted seismic lines from both Israeli and Jordanian territories. The model distinguishes seven units from upper Eocene to the Lower Triassic, accounting for major aquifers, aquicludes and deep-cutting faults. Recharges are implemented based on the numerical representation developed by Shentsis (1990) that considers relationships between mean annual rain and topographic elevation. The model reveals that topography-driven N-S and E-W flows strongly control the location of discharge areas while the anomalous spring temperature is not necessarily linked to the presence of fault convection. Local permeability anisotropy due to aquifers folding or facies changes are features sufficient for the rising of hot fluids.

Description: In the present study, inverse problems (IP) are applied with FEPEST® in FEFLOW® to find parameter distributions of hydraulic and thermal conductivity that lead to the observed thermal anomalies in the Lower Yarmouk Gorge (LYG) and a more accurate fit of borehole data. Results indicate reduced hydraulic conductivities in shallower parts of the system than previously estimated. The absence of a fault under the Lower Yarmouk Gorge in deeper sediments can be superseded by presence of interconnecting local fractures. The contribution of the conductive thermal regime turns out to be negligible.

Description: Previous investigations postulate the presence of a zone of high hydraulic anisotropy parallel to the principle axis of the Lower Yarmouk Gorge (LYG) in the Middle East. Driven by mixed convection, heated fresh groundwaters ascend within the gorge from confined Cretaceous units through artesian wells and Hammat Gader springs. Two-dimensional transient numerical simulations of coupled fluid flow and heat transport processes are used to investigate the impact of (1) a zone of hydraulic anisotropy and (2) abstraction on hydraulic heads and temperature profiles in the shallow aquifers. The models successfully reproduce hydraulic head distribution pre- and post-groundwater abstraction; dominance of conductive and advective heat transport processes is also shown. The models further support the existence of a structural feature along the principle axis of the gorge, which hydraulically connects groundwaters in both flanks, while cross flow of groundwaters is prevented. That implies a subsurface anisotropic zone, which lets the gorge act as a complex conduit-barrier system where adjacent N–S and S–N flow-fields confluence and get drained towards the Jordan Rift. The present numerical investigations support the hypothesis that, most likely, structural features that represent physical anisotropies control the hydrothermal system of the LYG. Furthermore, the study provides an example of numerical investigation of a complex transboundary aquifer system, with emphasis on existent anisotropies, structural ambivalence and restricted field accessibility.