A multi-environmental tracer study to determine groundwater residence times and recharge in a structurally complex multi-aquifer system

Despite being the main drinking water resource for over five million people, the water balance of the Eastern Mountain Aquifer system on the western side of the Dead Sea is poorly understood. The regional aquifer consists of 15 fractured and karstified limestone — aquifers of Cretaceous age and can be separated in Cenomanian aquifer (upper aquifer) and Albian aquifer (lower aquifer). Both aquifers are exposed along the mountain ridge around Jerusalem, which is the main recharge area. From here, the recharged groundwater flows in a highly karstified aquifer system towards the east, to discharge in springs in the Lower Jordan Valley and Dead Sea region. We investigated the Eastern Mountain Aquifer system on groundwater flow, groundwater age and potential mixtures, and groundwater recharge. We combined 36Cl/Cl, tritium and 20 the anthropogenic gases SF6, CFC-12 and CFC-11, CFC-113 as “dating” tracers to estimate the young water components inside the Eastern Mountain Aquifer system. By application of lumped parameter models, we verified young groundwater components from the last 10 to 30 years and an admixture of a groundwater component older than about 70 years. Concentrations of nitrate, Simazine® (Pesticide), Acesulfame K® (artificial sweetener) and Naproxen® (drug) in the groundwater were further indications of infiltration during the last 30 years. The combination of multiple environmental 25 tracers and lumped parameter modelling helped to understand the groundwater age distribution and to estimate recharge despite scarce data in this very complex hydrogeological setting. Our groundwater recharge rates support groundwater management of this politically difficult area and can be used to inform and calibrate ongoing groundwater flow models. https://doi.org/10.5194/hess-2019-451 Preprint. Discussion started: 9 September 2019 c © Author(s) 2019. CC BY 4.0 License.


Introduction
About 20 % of the Earth's land surface is covered by carbonate karst or sulfate aquifers and serves as the primary water resource for at least 25 % of the world's population. In addition, about a fifth of the world's karst systems are in (semi-)arid areas, whose water scarcity is aggravated by strong population growth (Ford and Williams, 2007). Karst systems represent abundant but highly variable water resources whose extremely heterogeneous and anisotropic flow behaviour prevents exact predictions regarding mass transport and the usable water quantities (Bakalowicz, 2005). Nevertheless, the water balance meets the basic requirement for sustainable management and protection of any water resource.
Environmental tracers play an important role in sustainable water management strategies because they allow for estimating the groundwater age distribution with depth, and together with simple lumped parameter models they can be used to quantify groundwater infiltration rates (Vogel, 1967;Solomon et al., 1995;IAEA, 2006). Especially in karst aquifers, wide ranges of residence times are observable due to the strongly heterogeneous hydraulic system, allowing water to rapidly flow through conduits and fractures and very slowly flow through the small pores of the matrix. That leads to large heterogeneities in the groundwater age distribution, hence requiring the application of multiple tracers to constrain the age distribution, especially in karst. We define young groundwater as having measurable concentrations of anthropogenic tracers, thus a mixing component recharged after about 1950, while in old groundwater these tracers are not detectable (e.g. Plummer et al., 1993;Cook and Herczeg, 2000). The atmospheric tracer CFC-11 (chlorofluorocarbon), CFC-12, CFC-113, SF 6 (sulfur hexafluoride) and 36 Cl / Cl and tritium from bomb tests or anthropogenic organic trace pollutants like pesticides, sweeteners or drugs are increasingly used as tracers of young groundwater (IAEA 2006). Gas tracers like CFCs and SF 6 move through the unsaturated zone primary by diffusion, leading to a time lag at the water table compared to the atmosphere . A time lag is also possible for the water bound tracers tritium and 36 Cl, since the advection through the unsaturated zone may take decades (Suckow et al., 1993;Lin and Wei, 2006), and in infiltration areas dominated by sand or clay, water-bound tracers are generally slower than gas tracers (Solomon et al., 1992;. This can be very different in karst systems, where preferential flow in karst "tubes" allows for a fast recharge to the groundwater table, and fluctuations of groundwater level may allow for further gas exchange thereafter. The deconvolution of measured tracer concentrations into recharge rates therefore needs modelling. If the age of water would be known as function of depth, any flow model could be directly constrained, and the recharge rates could be deduced. However, the "idealized groundwater age", which is often understood as the time span an imaginary water parcel needs between infiltrating at the groundwater surface and being sampled at a well or spring (Suckow, 2014a) is not directly measurable. In addition, groundwater mixes both along its natural flow through the aquifer and during sampling in the well. Therefore, simple lumped parameter models (LPMs) are used to interpret the measured tracer concentrations as mean residence times (MRTs) via a convolution integral, which in combination with the underlying assumptions on the flow system allow for deducing recharge rates. We applied the piston flow model (PM), the dispersion model (DM) and the partial exponential model (PEM) to approximate the age distribution in our groundwater samples.
In our study, the Eastern Mountain Aquifer (EMA) system in the western Dead Sea catchment is the pivotal water resource for some million people in the West Bank and Israel. An unequal distribution of borehole information leads to poor and limited data for studying that aquifer system. Previous studies considered age dating tracers to quantify water movement and flow velocity within the EMA and associated aquifers. Paul et al. (1986) and Yechieli et al. (1996) studied 36 Cl / Cl to detect very old groundwater brines in the Dead Sea area. In contrast young-age dating tracers such as tritium or anthropogenic trace gases (CFCs and SF 6 ) were used to quantify the duration of water flow from recharge areas to the springs (Lange, 2011). Environmental tracer investigations of the main Cretaceous aquifers (Upper Cenomanian and Albian) in the western Dead Sea catchment attempted to quantify the duration of water flow from the recharge area to the springs in the mountain region uphill of the Dead Sea coast. The young-age dating tracers demonstrated a large young groundwater fraction with a mean residence time of less than 30 years in the springs of the mountain region and fast connections to the recharge area. All previous studies together show large heterogeneities in the groundwater age distribution (Avrahamov et al., 2018).
In this study we combine for the first time in this area bomb-derived 36 Cl, anthropogenic organic trace substances and environmental tracers like tritium, CFCs, SF 6 in combination with lumped parameter models to interpret the distribution of these tracers to quantify recharge. Małoszewski and Zuber (1982 have shown that LPMs are a useful tool for interpreting tracer data obtained at separate sampling sites when it is not possible to use distributed parameter models, as the latter require more detailed and often unavailable knowledge about distributed parameters for the investigated system. In detail this work aims (i) to validate young rainwater input and short groundwater travel times via karst conduits, related to rapid flow paths from the recharge area; (ii) to quantify the time lag of gas tracers in the unsaturated zone; (iii) to quantify groundwater mixing of groundwater components with different ages via lumped parameter models; and (vi) to estimate groundwater recharge and support calculations of future groundwater resource development.

Study area
The study area, which represents the western surface drainage basin of the Dead Sea, is embedded in a region that is morphologically and geologically dominated by the tectonic processes associated with the Jordan-Dead Sea rift, active since the late Oligocene (Garfunkel et al., 1981;Rosenfeld and Hirsch, 2005). The western rift fault separates the Cretaceous aquifer formations that form the graben shoulder from the deeply subsided graben and its Quaternary filling. In addition, rift tectonics induced a series of faults within the western graben shoulder, resulting in down-faulted blocks, which find their surface expression in a strong morphologi- Figure 1. The location of the study site including geological information based on Begin (1974), Mor and Burg (2000), Raz (1986), Roth (1973), Shachnai (2000), Avni (2011), andRoth (2012). Universal Transverse Mercator (UTM) coordinates: UTM WGS84 Z36N. cal gradient. Within less than 25 km, the land surface drops from +800 m m.s.l. in the west to sea level at the rift margin, with a terminal step to −430 m m.s.l. at the Dead Sea in the east (Fig. 1).
The semi-arid-to-arid Mediterranean climate leads to precipitation during the winter season, with a strong decline from west to east due to which the study area can be divided into three hydrological zones: (i) the recharge area in the upland that receives annually up to 580 mm of precipitation; (ii) the transition zone occupying the hillsides of the upland down to the rift margin, receiving 100-400 mm yr −1 and (iii) the major discharge area of groundwater in the lower Jordan Valley-Dead Sea area, receiving less than 100 mm yr −1 of precipitation.
The graben shoulder hosts a thick aquifer system that is mainly built of fractured and layered karstified Upper Cetaceous limestone and dolostone, which are overlain by a Senonian chalky aquitard. A marly aquiclude (Lower Cenomanian) divides the system into a lower aquifer (Albian) and upper aquifer (Upper Cenomanian). Impervious beds in the upper aquifer permit the development of a perched and locally important aquifer (Turonian), being built of homogenous and fissured limestones and holding springs, which emerge in deeply incised valleys within the transition zone (e.g. Wadi Qilt springs; Fig. 2). While all aquifer units are recharged in the mountain area, the natural discharge of the two regional aquifers occurs through springs at the base of the graben shoulders, where the groundwater leaves the aquifers and approaches the prevalently impervious Quaternary graben filling.
Subsequently, groundwater emerges along the shore in spring clusters, forming ecologically important oases such as Ein Feshkha. Where groundwater can percolate into the Quaternary sediments; they intensely dissolve the evaporite minerals contained within (halite, anhydrite and aragonite) and get saline on their flow path to the lake shore. In addition, the chemical and isotopic composition of approaching fresh groundwaters is systematically modified by an admixing of brines, already in the vicinity of the major rift fault (Katz and Kolodny, 1989;Stein et al., 1997;Ghanem, 1999;Yechieli, 2000;Klein-BenDavid et al., 2004;Khayat et al., 2006a, b;Möller et al., 2007;Siebert et al., 2014;Starinsky and Katz, 2014). Groundwater recharge rates are controlled by climate conditions and have been investigated earlier to force regional groundwater models (Guttman et al., 2004;Yellin-Dror et al., 2008;Gräbe et al., 2013;Schmidt et al., 2014), which allowed for a detailed insight into the regional groundwater flow dynamics (Laronne Ben-Itzhak and Gvirtzman, 2005;Sachse, 2017).
Human groundwater abstraction takes place mainly in the mountain ridge inside the recharge area and along the transition to the Jordan Valley. This unequal distribution of the  et al., 1990), SF 6 and CFCs (Plummer et al., 2006) and decay-corrected Tritium in rain water of Beit Dagan and Vienna (IAEA/WMO, 2019).
sampling possibilities led to a data scarcity for the entire central area of the aquifer, which is also visible in Fig. 1.

Fundamentals of the method
Age distributions of young groundwater can be characterized by applying anthropogenic trace gases like CFC-11, CFC-12, CFC-113 and SF 6 in lumped parameter models if an input function is available. For the last 6 decades, that function has been derived for gas tracers from (i) their known concentrations in the atmosphere (Fig. 2), (ii) the observation that they are well mixed in the atmosphere and (iii) the observation that their solubility at the temperature of recharge is known from Henry's law (Plummer and Busenberg 2000). The large-scale production of CFC-11 and CFC-12 (both used as cooling fluid) started in the early 1940s, while the production of CFC-113 only started in the 1960s. Inevitably they leaked into the environment, with atmospheric concentrations rising until the 1990s, when their moratorium took effect. Another industrial gas SF 6 , widely used as an electrical insulator, has been detectable in the atmosphere since 1960s with still exponentially rising concentrations. Atmospheric concentrations of SF 6 and the CFCs for the period 1953-2006 ( Fig. 2) are derived from Plummer et al. (2006). Since gas solubility in recharging precipitation depends on temperature and atmospheric pressure, the average air temperature of the winter season in the study area (15 • C) was used for the former and an altitude of 700 m m.s.l., which is the average altitude of the infiltration area, was used for the latter. The rain water salinity was set to 0 ‰ since rain water is low mineralized.
Tritium ( 3 H) is the naturally occurring isotope of hydrogen and is mainly produced by fast secondary neutrons from cosmic radiation. It decays to 3 He with a half-life of 12.32 ± 0.02 years (Lucas and Unterweger, 2000), making 3 H usable for groundwater age dating in a time frame of < 40 years (e.g. Schlosser et al., 1988;Solomon et al., 1992;Cook and Solomon, 1997;Sültenfuß and Massmann, 2004), including quantifying changes of measured 3 He / 4 He ratios in groundwater. At the study area, observations of tritium in precipitation (Global Network of Isotopes in Precipitation stations Beersheba, Beit Dagan and Tirat Yael; GNIP; IAEA/WMO, 2019) were available only for 1960-2001, while input data are required until 2014. We therefore applied tritium data from the Vienna station, Austria (IAEA/WMO, 2019), which are adjusted to the longitudinal and latitudinal difference by a factor of 0.4 to match the stations in Israel (Fig. 2). For the synoptic tracer plots the decay correction is set to 31 October 2013. The pre-bomb input value for tritium was set to a mean tritium concentration of 3 TU obtained from the GNIP database (IAEA/WMO, 2019). 36 Cl is produced naturally via cosmic-ray and solar protons induced by nuclear reactions of argon in the atmosphere and 35 Cl in marine aerosols (Alvarado et al., 2005). However, comparable to 3 H, the atmospheric concentration of "bomb" 36 Cl peaked during the 1950s as an effect of nuclear weapon tests and was washed out from the atmosphere by the end of the 1960s. The 36 Cl bomb peak precedes the tritium peak by half a decade. The 36 Cl input curve for our study area ( Fig. 2) was obtained from Iceland ice core measurements from Synal et al. (1990), which were corrected to the location of the study area by applying a latitudinal correction with a factor of 3, according to Heikkilä et al. (2009), who modelled 36 Cl fallout for different latitudes. Natural "pre-bomb" concentrations of 36 Cl / Cl were assumed to be 10 −14 , which is based on 36 Cl / Cl in rainwater, sampled during the winter of 2014-2015 with an average value of 36 Cl / Cl = 8×10 −15 .
Anthropogenic organic trace pollutants in groundwater are associated with nutrition, medication, or agricultural and industrial development. During the last decades, artificial sweeteners have played an important role as a surrogate in the nutrition industry. Particularly acesulfame K (ACE-K) has been used since the 1990s and is stable against wastewater treatment (WWT) processes, making it an ideal tracer for domestic wastewater. A second chemical marker for human intake is naproxen (NAP), widely applied as an antiinflammatory drug since the 1980s. Though NAP is partly eliminated during WWT and may be adsorbed along flow paths to sediments (Chefetz et al., 2008;Yu et al., 2008;Teijón et al., 2013), it occurs in effluents of sewage plants and increasingly in natural waters (Arany et al., 2013). Contrasting that, ACE-K is hydrophilic (e.g. Buerge and Poiger, 2011) and inert against degradation, which makes it a valuable substance to trace transport from the recharge towards the discharge area. In addition to the urban indicators, pesticide traces in groundwater were used to identify the agricultural contributions to the water resources. Simazine, which has been available since the 1950s, is one of the most applied herbicides and absorbs into soil where it may be eliminated though bacterial degradation. However, the use of simazine has been phased out in Israel since 2012-2014 (Berman et al., 2014). Nitrate is another indicator for anthropogenic input, whose distribution in groundwater originates from the nitrification of NH 4 during WWT and fertilizer excess in agriculture.

Sampling and analytical methods
In the study area, 22 groundwater samples were taken from springs and active wells (sampling locations in Fig. 1 and Table 1) after reaching stable conditions for temperature, electrical conductivity and pH. Major anion samples were filled into high-density polyethylene (HDPE) bottles after passing a 0.45 µm cellulose acetate filter. Samples for cations were acidified using HNO 3 . Samples for 36 Cl analyses are filled into 500 mL HDPE bottles, which have been specifically precleaned using ultrapure HNO 3 . Samples for 36 Cl analyses were acidified with HNO 3 . Following the methodology de-scribed in Oster et al. (1996), sampling for CFCs and SF 6 was performed using glass bottles fully submerged in tins, filled with sampling water. Tritium samples were collected in 500 mL HDPE bottles. Organic trace elements were sampled in pre-cleaned, methanol-flushed 1000 mL brown glass bottles with the use of a glass microfiber filter (0.7 µm).
Major element analyses in water samples were performed at the Helmholtz Centre for Environmental Research (UFZ) by using matrix-adjusted inductively coupled plasma optical emission spectroscopy (ICP-AES, Spectro Arcos) for Na + , K + , Ca 2+ , Mg 2+ and Sr 2+ and by using ion chromatography (ICS-2000, Dionex) for Cl − , Br − and SO 2− 4 . Bicarbonate was determined in situ by titration. Analyses of 36 Cl were carried out in the Helmholtz-Zentrum Dresden-Rossendorf at the accelerator mass spectrometry (AMS) facility DREsden AMS (DREAMS) (Akhmadaliev et al., 2013). The main preparation steps consist of (i) the precipitation of chloride by adding an AgNO 3 solution (10 %) and subsequent dissolution of the AgCl in NH 4aq ; (ii) the separation of chlorides from sulfates by co-precipitation of BaSO 4 with BaCO 3 (CO 2 from air) using a saturated BaNO 3 solution and filtration through a syringe filter made of polyvinylidene fluoride (pore size: 0.45 µm); and (iii) the re-precipitation of AgCl in the filtrate with nitric acid. 36 Cl is measured with AMS relative to the stable Cl isotopes, 35 Cl and 37 Cl (Pavetich et al., 2014;Rugel et al., 2016), and is given as the ratio of 36 Cl/( 37 Cl+ 35 Cl) (termed 36 Cl / Cl in the text). All data are normalized to the standards of SM-Cl-12 and SM-Cl-13 (Merchel et al., 2011).
Tritium sample preparation and measurement were conducted by the isotope hydrology group at UFZ, following the preparation steps of Trettin et al. (2002) by first enriching a 400 mL water sample electrolytically and measuring it via liquid scintillation counting with a detection limit of 0.5 TU.
CFCs and SF 6 were analysed in the Spurenstofflabor Dr. Harald Oster (Wachenheim, Germany) by gas chromatography (Bullister and Weiss, 1988;Oster et al., 1996). Measurements of organic trace elements were conducted with high-pressure liquid chromatography mass spectrometry (HPLC-MS) at the UFZ. To transport the samples to the lab, organic components were stabilized by solid phase extraction (SPE) on cartridges containing a polar-modified polystyrene-divinylbenzene copolymer (Chromabond Easy, Macherey-Nagel). At the lab, analytes were eluted with methanol and measured with HPLC-MS/MS (Agilent 1000, Agilent Technologies, Germany; coupled with an API2000 mass spectrometer, AB Sciex, Germany). Limits of detections are 1, 0.6 and 0.3 µg L −1 for naproxen, acesulfame K and simazine, respectively. Analytical results are given in Table 1.

Lumped parameter model "LUMPY"
LPMs are pre-defined analytical solutions of simplified flow systems. They describe the tracer output mathematically with  a convolution integral that combines the tracer input history, weighed with the age distribution valid for the flow system in question (Maloszewski andZuber, 1982, 2002).
Mathematically all lumped parameter models for steadystate flow systems with a time-variable tracer input are convolution integrals (Eq. 1).
where t is the calendar time, t is the transit time of the tracer, C out is the output concentration, C in is the input concentration and g(t ) is the weighting function or system response function. All weighting functions of all models are normalized, following Eq. (2).
The mean residence time is the main fitting parameter, while some models (e.g. dispersion model) require additional parameters like the (i) Péclet number, (ii) top and bottom of the screened section or (iii) saturated thickness of the aquifer. Input data are (i) regional atmospheric tracer input curves and (ii) selected hydrogeological characteristics such as the infiltration temperature and elevation. In this study the LPM code LUMPY (Suckow, 2012) was used to implement the convolution integral for flow systems that can be described by the piston flow model, the exponential model (including the partial exponential model) and the dispersion model.
The PM describes the movement of a water parcel along a defined flow path from the aquifer surface towards the spring or well filter, neglecting any mixing, dispersion or diffusion. The DM characterizes transport influenced by dispersion and advective flow. The relative magnitude of both is expressed as the Péclet number Pe (Eq. 3).
where l is the flow length of the system under consideration, ν is the velocity and D is the dispersion constant (Huysmans and Dassargues, 2005). In our model approach, best fits were obtained by applying a Péclet number of 30, characterizing a dominantly advective transport. The PEM is related to the exponential model (EM). Based on homogenous infiltration into a homogeneous aquifer as in Vogel (1967), the PEM describes mixing of those flow lines reaching the filter screen of a well. In the special case where the filter screen extends over the whole thickness of the aquifer the PEM is equivalent to the exponential model, and the MRT (mean residence time of the sample) is the only fitting parameter. The mathematical equation of DM, EM and PEM are described in Maloszewski andZuber (1982, 2002) and Jurgens et al. (2016), respectively.

Parameterization and model setup of LUMPY
To parameterize the unsaturated (vadose) zone in the recharge area, characteristic wells (Samia 2 and Herodion 1) were used ( Table 2).

Influence of a thick unsaturated zone
Gas tracers like CFCs and SF 6 predominantly pass the unsaturated (vadose) zone in the gas phase, posing certain problems for their interpretation. In unsaturated zones of less than 5m thickness, the gas composition of soil air resembles that of the atmosphere Engesgaard et al., 2004). However, a time lag may occur for the diffusive transport of CFCs and SF 6 through thick unsaturated zones of porous aquifers . This time lag is a function of the tracer diffusion coefficients, tracer solubility in water and moisture content (Weeks et al., 1982;. A time lag always results in a gas tracer age older than the time of groundwater recharge. In fractured (or karstic) aquifers, however, the time lag may be much shorter (e.g. Darling and Gooddy, 2007;Gooddy et al., 2006) resulting in ages of CFCs or SF 6 measurements obtained from groundwater which effectively represent residence time of groundwater since recharge approached the groundwater table, without a time lag and as if the tracer were transported within the saturated zone only. The present study aims to estimate travel times in the thick unsaturated zone Plummer et al., 2006) by estimating the time lag as the relative difference in mean residence time between gas and water-bound tracers, applying both, water-bound tracers ( 3 H and 36 Cl) and gas tracers (CFCs and SF 6 ). In particular, the wells Samia 2 and Herodion 1 of the Upper Cenomanian are selected in order to con-sider wells close to the recharge area with a thick unsaturated zone of 200-240 m (Table 2). They are in the eastern part of the recharge zone (Fig. 1), where the limestone is intensely fissured vertically and partly karstified.

Well construction, aquifer data and the calculation of recharge rates
The convolution integral of the partial exponential model can be further constrained by the well construction data like (saturated) depth to the top of the screen, the screen length and saturated aquifer thickness (Jurgens et al., 2016). These data were taken from the construction logs of the investigated production wells (Table 2). Based on the estimated MRT of the applied lumped parameter models DM, PM and PEM, recharge rates can be estimated. However, the formulas to apply differ slightly between the different models. The following Eq. (4) allows for calculating the recharge rate R for the PEM (Vogel, 1967;Jurgens et al., 2016): where φ is the porosity of the aquifer, L is the saturated aquifer thickness and the MRT is valid for the whole aquifer (which is an output in LUMPY derived from the fitted MRT for the sample and using the well parameters Z 1 , Z 2 and L in Table 2). As for the DM and PM, L is the distance from the groundwater surface to the depth of the centre of the screened section of the well. Depth to the centre of the water-filled screen section from groundwater level is 23 m for Samia 2 and 53.5 m for Herodion 1 ( Table 2). As for the DM, not the MRT but the peak time, is used, which is also an output of LUMPY, being calculated by applying Eq. (5) (Suckow, 2014b).

Peak Time =
MRT P e 9 + P 2 e − 3 Tracer transport in karst areas is influenced by double porosity effects, particularly a retardation of the tracer can be assumed due to diffusive loss into the adjacent limestone of the fissures. A correction of calculated MRTs and recharge rates may be possible applying a retardation factor (Eq. 6) (Maloszewski et al., 2004;Purtschert et al., 2013).
where ∅tot ist the total porosity and ∅eff is the effective porosity. An assumed porosity in the carbonate karst aquifers may vary between < 2 % (representing the open fissures and solution pipes) to 20 % (total pore space in the aquifer rock). Groundwater recharge was also estimated using the chloride mass balance (CMB) method, which was successfully applied elsewhere (Eriksson and Khunakasem, 1969;Allison and Hughes, 1978;Wood and Sanford, 1995;Purtschert et al., 2013;Crosbie et al., 2018) and which assumes Cl input to groundwater originates from Cl concentration in precipitation, which becomes enriched due to evaporation only (Eq. 7).
The formula consists of the mean annual precipitation P and the chloride concentrations in precipitation Cl P and groundwater Cl GW in mg L −1 . The mean annual precipitation in the recharge area is about 550 mm. The average long-term Cl content in rain water is ca. 5 mg L −1 (Herut et al., 1992), while it reaches 28 and 35-38 mg L −1 in the groundwaters of wells Herodion 1 and Samia 2, respectively.

Results
The 36 Cl / Cl ratios in precipitation, which fall in the recharge area, are assumed to have been stable since the 1980s, as indicated by rainwater samples collected during the winter of 2014-2015, which show 36 Cl / Cl ratios of 1.3 × 10 −14 to 5.5 × 10 −14 , resembling results from the early 1980s (Herut et al., 1992). Contrastingly, the tritium concentration in precipitation continuously declined to about 4-6 TU today (IAEA/WMO, 2019). Fitting the hydrogeological trichotomy of the study area well, analytical results (Table 2) resemble the regional situation and group according to the individual aquifers (Fig. 3).

Recharge area
In the wells Samia 2 and Herodion 1, representative for the recharge area of the Upper Cenomanian aquifer, tritium concentrations of 2.1 and 1.2 TU as well as 36 Cl / Cl values of 1.41 × 10 −13 and 1.57 × 10 −13 , respectively, are observable. Further south, in the well Al Rehiyya, groundwater in the aquifer show lower 36 Cl / Cl and 3 H values of 1.73 × 10 −14 and < 0.5 TU, respectively. The gas tracer concentration in that part of the aquifer is low but detectable and showed SF 6 = 1 ± 0.2 fmol L −1 , CFC-11 = 2.3 ± 0.3 pmol L −1 , CFC-12 = 1.2 ± 0.1 pmol L −1 and CFC-113 = 0.19±0.05 pmol L −1 in Samia 2 and comparable values in Herodion 1 ( Table 1). The well Jerusalem 1, which represents the Lower Cenomanian aquifer extract groundwater, showed 3 H concentrations (1.7 TU) as low as those observable in the upper aquifer, but it had much higher 36 Cl / Cl values (8.9 × 10 −13 ).

Transition zone
Groundwater, emerging from the perched Turonian aquifer (Ein Auja) and from the Upper Cenomanian aquifer (Mizpe Jericho 2), shows a low 36 Cl / Cl ratio of 6.35 × 10 −14 and 7.99 × 10 −14 and 3 H concentrations of 2.9 and 1.4 TU, respectively. The well Mizpe Jericho 2 stands out due to its high CFC contents, reaching values of CFC-11 = 25 pmol L −1 and CFC-113 = 4 pmol L −1 , higher than possible in equilibrium with the atmosphere (4.5 and 0.5 pmol L −1 , respectively). Further south, the Arugot spring discharges in the upper Arugot Valley, about 300 m above the Dead Sea from the Upper Cenomanian aquifer, closely to the brim of the graben flank but within the transition zone. The emerging groundwater is high in SF 6 (2.1 fmol L −1 ), CFC-11 (3.8 pmol L −1 ) and CFC-12 (1.8 pmol L −1 ) but contains very low 3 H (0.5 TU), suggesting a well-developed karst network which allows for sufficient gas exchange of older groundwater with recent atmosphere along its flow path. Groundwater in the Albian aquifer (wells PWA 3, Herodion 4 and Bani Naim 3) has a low concentration regarding CFCs, SF 6 and 3 H (< 1 TU) and shows low 36 Cl / Cl ratios of 1.9-4.9 × 10 −14 . An exception is well Azaria 3, at the edge between the recharge and transition zones, which shows a higher 36 Cl / Cl value (9.59 × 10 −14 ) than recent precipitation, much less tritium (0.5 TU) and high concentrations of CFC-12 (5.2 pmol L −1 ), which are again, well above values in equilibrium with the atmosphere (max 2.3 pmol L −1 ).

Discharge area
Groundwater from the Upper Cenomanian aquifer (Jericho 2) shows similarly low 36 Cl / Cl ratios of 2.04 × 10 −14 and 2 TU similar to groundwater upstream in the transition zone. In groundwaters pumped in Jericho 4 and 5 from the Albian aquifer, 36 Cl / Cl and 3 H contents are even lower: 1.03-1.15 × 10 −14 and < 0.6 TU, respectively. These low 36 Cl / Cl values might result from admixing brines, which are abundant within the rift. It becomes evident in the groundwater of Ein Feshkha, where interstitial brines ( 36 Cl / Cl = 1.08 × 10 −14 ) hosted in the interstitial space of the Quater-nary sediment, get leached by approaching fresh groundwaters and cause the lowest observable 36 Cl / Cl values of 1.01-1.51 × 10 −14 .
In general, the following patterns can be extracted from the hydrochemical and tracer data: i. spring water of the Upper Cenomanian (Ein Auja) contains groundwater infiltrated after the era of atmospheric bomb testing, indicated by tritium and 36 Cl / Cl values comparable with recent rain water; ii. groundwater of the Upper Cenomanian close to the recharge area shows 36 Cl / Cl ratios higher than recent rain combined with low tritium contents, possibly referring to admixtures of water from the early fission bomb testing but before thermonuclear devices; iii. along the flow path and in the lower aquifer 36 Cl / Cl values are shifted to lower ratios due to admixture of saline water, with the end-member being brine similar to the Dead Sea probably admixed as pore water from earlier higher sea-levels; iv  ing no recent recharge, although SF 6 and CFC values are high, as recent gas exchange in karst structures may have reset the SF 6 and CFC "clocks".
Anthropogenic organic trace substances simazine, NAP and ACE-K are detectable in trace concentrations in nearly all sampled wells of the upper and lower aquifer (Fig. 4), indicating input younger than 40 years. Combining nitrate contents with concentrations of the herbicide simazine in the sampled groundwaters allows for distinguishing between the origin of NO − 3 , either from domestic wastewater or from agriculture (Fig. 4a). All groundwaters show positive corre-lations between nitrate and simazine, while there is no clear systematic trend observable that is specific for one of the aquifers nor a specific region. Studying the samples according to aquifers suggests that groundwater in the Upper Cenomanian aquifer is more strongly influenced by wastewater inflow than groundwaters in the Albian aquifer. Groundwater from Mizpe Jericho 2, with high NO − 3 (43 mg L −1 ) values at comparably low simazine concentrations clearly underlines the contribution of wastewater. Contrastingly, groundwaters from Jericho 2 and Jerusalem 1 may refer to a higher contribution from agriculture.
Looking at ACE-K and NAP as pure wastewater indicators (Fig. 4b), two distinctly opposite trends are distinguishable. Trend 1 is characterized by very low NAP concentrations and high concentrations of ACE-K as observable in the wells Jericho 4 (56 ngL˙− 1 ) and Jericho 2 (136 ng L −1 ). The opposite trend is found in Ein Feshkha spring D, Jerusalem 1 and Samia 2, with high NAP concentrations (170, 168 and 137 ng L −1 , respectively) at low ACE-K concentrations (< 20 ng L −1 ). If one excludes Ein Feshkha, NAP concentrations decrease from the recharge area (Jerusalem 1 and Samia 2) to the discharge area (Jericho 2 and 4), most probably due to adsorption onto the aquifer matrix along the flow path. As for Ein Feshkha D, the NAP contamination here is comparably high in the recharge area and much larger than the NAP contents in Ein Feshkha C and Enot Zukim. This suggests a significantly shorter residence time of the contaminant in the aquifer, which requires a source much closer to the spring, probably a leakage in the treated wastewater (TWW) pipeline, passing the area just upstream of the spring. The only groundwater without any anthropogenic contamination is the Arugot spring.

Discussion
The patterns observed for the different measured substances deliver a heterogenous picture of the study area. Starting at the top, a shallow perched aquifer system with short residence times is indicated for the Ein Auja spring based on a high 36 Cl / Cl ratio (6.35 × 10 −14 ), tritium content (2.9 TU) and low mineralization like recent precipitation.
In the recharge area of the Cenomanian and Albian aquifers, clear indications for a contribution of recharge enriched with 36 Cl from nuclear-bomb tests is observable in freshwater wells (Cl content <67 mg L −1 ) Jerusalem 1, Herodion 1 and Samia 2, which show 36 Cl / Cl values of 1.4-8.9 × 10 −13 , which are much larger than in precipitation from after 1980 ( 36 Cl / Cl =< 5.5 × 10 −14 ). Furthermore, gas tracer results in Herodion 1 and Samia 2 imply that the aquifer contains groundwater recharged after 1960. Organic pollutants concretize it even more. High simazine content in Jerusalem 1 and high loads of NAP in Samia 2 and Jerusalem 1 indicate significant contamination through a young water fraction in both aquifers, the Upper Cenomanian and the Albian.
Within the transition zone, the wells Al Azaria 3 (Albian) and Mizpe Jericho 2 (Upper Cenomanian) show 36 Cl / Cl values still larger than in precipitation from the last four decades, indicating at least similar or even older recharge periods than of groundwater in the recharge zone. However, high concentrations of CFC-11 and ACE-K in Mizpe Jericho 2 and high concentrations of CFC-12 and NAP in Al Zaria 3 also show a significant contribution of the young (waste-)water fraction to these wells. Contrastingly, anthropogenic gas tracers are close to or even below the limit of detection in groundwater samples from Albian aquifer at the wells PWA 3, Herodion 4 and Bani Naim 3, indicating no freshwater input from the last decades and travel times longer than 70 years. A similar figure results from taking these three wells and forming a N-S transect through the Albian aquifer in the mountain range: their low 36 Cl / Cl ratios, which are well within the range of recent precipitation, decrease from north to south, while chloride remains stable. Since the 3 H content is below the limit of detection (< 0.5 TU), the observed groundwaters are considered to be mainly pre-bomb water and hence older than 6-7 decades.
In the southern part of the study area, survey stations are even rarer. However, the well Al Rehiyya shows very low 36 Cl / Cl and low 3 H values, indicating either unmixed prebomb water or the dilution of recent precipitation with much older water possibly originating from aquifer parts even further southwards. The latter is promoted by results from numerical flow modelling, which indicate an SW-NE-directed upstream flow (Laronne Ben-Itzhak and Gvirtzman, 2005; Gräbe et al., 2013). Negligible concentrations of NAP, ACE-K and simazine in Arugot spring water, which emerges at the furthest end of the southern transition zone, show no anthropogenic contamination by wastewater. This leads to the assumption that open karst conduits allow for efficient exchange with the atmosphere along the subsurface flow of the spring water, resulting in CFC-11 and CFC-12 content in equilibrium with the atmosphere.
Groundwaters from the Albian (Jericho wells 4 and 5) and from Upper Cenomanian (Jericho 2) in the discharge area are characterized by 36 Cl / Cl < 2 × 10 −14 values. While the former are free of tritium, Jericho 2 shows 2 TU, indicating at least a fraction of younger water. Since the well Jericho 2 is drilled directly at the mouth of the Qilt Valley, we assume that the young water fraction may reach the well through rapid infiltration through fractures and karst conduits within the valley. The young water fraction in Jericho 2 is also proven by a high simazine and ACE-K concentration in the groundwater of that well.
Low 36 Cl / Cl ratios in the wells drilled into the Quaternary section east of Jericho (ca. 2 × 10 −14 ) and in the Ein Feshkha springs (ca. 1 × 10 −14 ) refer to the enrichment of fresh groundwater with different sources of salinity, showing low 36 Cl / Cl ratios. The brackish wells east of Jericho are known to draw groundwater from the Upper Cenomanian aquifer (Khayat et al., 2006a), which passes the graben fault and enters into the sediment body, where it may leach abundant evaporitic minerals. It also contains water from agricultural irrigation and wastewater as confirmed by high NO 3 concentrations, the presence of simazine and occasionally remarkable ACE-K content.
Contrastingly, Dead Sea brines, which get released from the sedimentary body, mix in Ein Feshkha with approaching fresh groundwaters from the graben shoulder. The low time, lake level was at −390 m m.s.l. and higher (Hassan and Klein, 2002). At the time, Dead Sea brine infiltrated into the shallow lake bed, where nowadays the observed springs emerge, at an elevation of −392 to −395 m m.s.l. However, the Ein Feshkha springs C and D and Enot Zukim receive young water fractions, since NAP and ACE-K contents in the former are as high as in groundwater in the recharge zone (Jerusalem 1) and in the latter at least enhanced.
The next section will exemplify a more detailed modelling of mean residence times with transient tracers where the data allow for this approach. This was possible only on two wells. The wells of the Albian aquifer in the transition zone are not interpreted using LPM due to contamination with CFCs (e.g. Al Azaria 3 and Table 1) or very low concentrations of anthropogenic gases (e.g. Bani Naim 3). If the measured concentrations of the anthropogenic trace gases CFCs and SF 6 are much higher than expected from the solubility equilibrium with the atmosphere, this is regarded as contamination (e.g. Mizpe Jericho 2). The tritium concentrations of 1.5 and 2 TU in Mizpe Jericho 2 and Jericho 2, respectively, indicate a reasonable fraction of recent rain. Due to the continuous outflow of the springs along the coastline and chemical mixing patterns of the spring water, it is assumed that there is also a connection to the water resources of the Lower Cretaceous aquifer, which could provide water with longer residence times. All springs in the discharge area (e.g. Ein Feshkha) have values less than 0.5 TU and a 36 Cl / Cl value of < 1.5 × 10 −14 . However, the organic tracers ACE-K, simazine and NAP are also detectable also in these springs. Their values in Ein Feshkha are between 5 to 30 times above the detection limit (Table 1 and Sect. 2.2). Since the input function of the organic tracers, cannot be quantified, the fraction of young water, representing that concentration cannot be calculated. This discrepancy (no tritium and no anthropogenic 36 Cl but including organic tracers) can be understood by quantifying the detection limit of "young water" that tritium allows: recent rain in the area has 4 TU, so a value of < 0.5 TU is equivalent to less than 12 % of recent rain. This implies the young end-member in the mixture must have a concentration of ACE-K, simazine and NAP > 10 times higher than that measured in the groundwater samples. A more precise quantification of the young water fraction in these springs is only possible with either a lower detection limit for tritium (e.g. determination of 3 H via ingrowth with a detection limit of 0.005 TU) (Bayer et al., 1989;Beyerle et al., 2000;Sültenfuß et al., 2009) or a precise quantification of the input of ACE-K, simazine and NAP. Both are beyond the scope of the present study.

Lumped parameter modelling: delay of gas tracer in the unsaturated zone
The detailed lumped parameter modelling of tritium, 36 Cl / Cl and the gas tracers was performed only for the wells Samia 2 and Herodion 1 since only for these there is a con- sistent data set for all tracers (Table 1). Earlier interpretation of groundwater hydraulics based on groundwater level measurements determined a very fast transfer velocity of the water phase through karst holes to the groundwater table (Jabreen et al., 2018). LPM was therefore done in several steps. First a mean residence time in agreement with measured values for tritium and the bomb spike of 36 Cl / Cl values was derived. Then a delay for the gas tracers CFCs and SF 6 was derived assuming a simple piston flow transport to describe the residence time of these tracers in the unsaturated zone. Once the delay created an agreement between the water-bound tracer tritium and 36 Cl / Cl values with the gas tracers CFC-11, CFC-12, CFC-113 and SF 6 , the mean residence time in the saturated zone was investigated using the PM, DM and PEM to describe flow in the saturated zone. From the MRT in the saturated zone, groundwater recharge was estimated in a final step and compared with results from the CMB method and earlier numerical groundwater models. Results for tritium and 36 Cl / Cl values fit to MRT of about 10 years in the saturated zone using the PM, the dispersion model results for a 20-year MRT and for the partial exponential model 16-and 20-year MRTs are estimated (Table 4). These estimated MRT values of the water-bound tracers then allowed for determining the gas delay in the unsaturated zone. For every applied model (PM, PEM and DM) the gas delay is estimated using the calculated model curve of the water-bound tracers and using the delay as a parameter to fit the gas tracer measurements of Samia 2 and Herodion 1. This resulted in specific gas delays for every concentration of the CFCs and SF 6 , exemplified for CFC-12 in Fig. 5 and provided for all gas tracers in Table 3. Estimates for every model result in different gas delays are related to the MRT of water-bound tracers in the wells Samia 2 and Herodion 1 ( Table 3). The delays of the CFCs are comparable to each other and differ from 14 to 24 years for Samia 2 (Fig. 5) and from 14 to 30 years for Herodion 1 over all calculated models. The gas delay in Samia 2 is lower than in Herodion 1, which is reasonable regarding the thicker unsaturated zone in Herodion 1 (Table 2). Ideally the transport of gas tracers in the subsurface would result in similar delays for all gas tracers . The delay determined for SF 6 was significantly lower than for the CFCs (Table 3), indicating higher concentrations of SF 6 as compared to the expectation from the CFC model results. This could be a result of biological degradation of the CFCs in the upper part of the unsaturated karst zone or from considerable excess air, which influences SF 6 much more than the CFCs. A decision between these two processes is possible determining excess air independently measuring the concentrations of all noble gases, which were not available in the present study.

Lumped parameter modelling: mean residence times in the saturated zone
In the next step, the determined gas delays of the gas tracers (Table 3) in the unsaturated zone were used to estimate the MRT of the water-bound tracers and the gas tracers in the saturated zone. Only the gas delay allowed for a fit to the measurement combination of the gas tracer concentrations and the water tracer results of Samia 2 in the time range of 10 to 15 years for the MRT for the piston flow model, 17 to 28 years for the MRT for the dispersion model and 26 to 46 years for the partial exponential model (Table 4). Comparing the different model approaches, estimates of the saturated MRT systematically increase in the sequence of PM to DM to PEM in both wells, which is a result of the increasing amount of dispersion and mixing in these models. In contrast to Samia 2, it was not possible to fit the results of Herodion 1 without an admixture of tracer-free old water. With and without delay the lumped parameter curves do not fit the measurements of the lower tritium concentra- Figure 6. Modelling results for groundwater from the Herodion 1 well, applying PM (red), PEM (black) and DM (blue) and (a) tritium vs. CFC-12 and (b) CFC-12 vs. SF 6 . The calculated delays of CFC-12 are 30, 20, and 28 years according to the models PM, DM and PEM, respectively (ref. Table 3). tion in Herodion 1 (1.2 TU), whereas they do for Samia 2. The well Herodion 1 is situated in the recharge zone, but it is also in the lower part of the Upper Cenomanian aquifer. The model results therefore indicate an admixture of at least one older water component ascending from the Albian aquifer and diluting the younger tritium concentrations. This can be described with a mixing line of two end-members, which fits the tritium measurement of Herodion 1 (Fig. 6).

Recharge rates and dual porosity
The mean residence times resulting from the lumped parameter models PM, DM and PEM in Table 4 were used to calculate recharge rates according to Eqs. (4) and (5). Within a single LPM the MRT have uncertainties in the range of 50 %, and this uncertainty increases to approximately a factor of 2 if all models are regarded as equally probable (Table 4). Therefore, a major uncertainty for the calculation of recharge results from the unknown porosity of the karst aquifer, which is estimated to be somewhere between 2 % and 20 %, creating an additional uncertainty of a factor of 10. It is therefore useful to consider whether this uncertainty can be further restricted.
Recharges rates from chloride mass balance highly depend on the value used for chloride in precipitation since the chloride concentration in the Cenomanian aquifer for youngest groundwater (tritium > 1.5 TU) varies only by 22 % between 28 and 44 mg L −1 ( Table 1). The Cl − values measured in precipitation during the present study vary between 5 and 30 mg L −1 resulting in a spread between 62 and 375 mm yr −1 for the recharge derived from CMB according to Eq. (7), assuming an average precipitation amount of 550 mm yr −1 . This is comparable to recharge rates used in numerical flow modelling studies of the eastern groundwater catchment at the Dead Sea, which estimated 124-292 mm yr −1 (Guttman, 2000) and 100-300 mm yr −1 (Gräbe et al., 2013) in the mountainous recharge area. Tracer-derived recharge values compare best with the recharge estimates from CMB and numerical modelling if a porosity of 10 % is assumed ( Table 5). The tracer-derived values for 20 % porosity are in the range of 40 % to 100 % of the mean annual precipitation of Jerusalem, and especially the latter value seems unrealistically high. The highest values originate from the PM and the PEM. The PM neglects mixing and retardation along the flow path completely, which of course is the least probable assumption. The PEM, however, considers mixing in the well, but it relies on the flow system being extremely homogeneous, as in the "Vogel" aquifer (Vogel, 1967;Jurgens et al., 2016). A few preferential flow paths can easily disturb this picture towards smaller MRT values and thus larger recharge values. An indication of why the tracer-derived recharge values using 2 % porosity are unrealistically low comes from the hydrogeology of the karst aquifer itself. The 2 % porosity can be attributed to structural fractures and macropores in the limestone, which account for the rapid reaction of the groundwater surface to precipitation events (Jabreen et al., 2018). However, according to Eq. (5) the tracers are not expected to indicate this flow, but they would be retarded by diffusion of the tracer into and out of the rock matrix (Maloszewski et al., 2004;Purtschert et al., 2013;Suckow et al., 2019). This means the tracers never "see" only the 2 % fracture and macropore space of the karst but also at least a part of the total rock matrix with its total porosity of 20 %. If the effective porosity in the karst system is 2 % and the total porosity is 20 %, the retardation factor according to Eq. (5) would be 10. Applied to the lumped parameter models, this retardation means that the calculated MRT values are too old, and the calculated recharge rates with 2 % porosity are too low.

Conclusions
The present study derived groundwater flow patterns, mixing end-members, transport times and recharge estimates in the upper and lower Cretaceous aquifers of the EMA between the recharge area around the central mountain ridge and the discharge zones close to the lower Jordan Valley and the Dead Sea, shared by Israel and Palestine. This was possible, despite a very low number of measurements and a complicated karst hydrogeological setting, using a powerful combination of multiple lines of evidence from hydrogeology, hydrochemistry, anthropogenic organic trace substances and classical environmental age-dating tracers like tritium, CFCs, SF 6 and 36 Cl / Cl values. The present study is one of presently only a handful demonstrating the useful application of atmospheric bomb-test-derived 36 Cl to study groundwater movement. The lower Cretaceous aquifer was found to be basically free of tritium and other anthropogenic environmental tracers like CFCs, SF 6 and bomb-derived 36 Cl and, therefore, it has groundwater transport times larger than 50 years. Groundwater in the upper Cretaceous aquifer contains anthropogenic trace substances, indicating time scales of groundwater flow of a few decades and shows clear indications of preferential flow paths as expected for karstified groundwater systems. In two springs the combination of several environmental tracers ( 3 H, bomb-derived 36 Cl, CFC-11, CFC-12, CFC-113 and SF 6 ) allowed for estimates of the mean residence time in the saturated zone, gas transfer delays in the unsaturated zone and mixing ratios with older groundwater. These estimates were confirmed by anthropogenic substances from agriculture (nitrate and pesticides like simazine), pain killers (naproxen) and sweeteners (acesulfame K). Calculated recharge rates based on MRT estimates compare well with chloride mass balance and numerical flow modelling. The multi-tracer methodology presented here is applicable in other data-sparse areas with complex hydrogeology (karst or fractured) with or without anthropogenic influence.
Data availability. The data on which the study is based are given in the paper (Table 1).
Author contributions. CW, TR, BM and CS outlined and conducted field work. SM, SP and GR conducted AMS measurements. CW, ASu and CS analysed the data. BM, TR, UM, CM, ASa and SMW were involved in discussion process. CW, ASu and CS wrote the paper. All authors revised the paper and approved the final version.
Competing interests. The authors declare that they have no conflict of interest.
Special issue statement. This article is part of the special issue "Environmental changes and hazards in the Dead Sea region (NHESS/ACP/HESS/SE inter-journal SI)". It does not belong to a conference.