Assessment of Hydrogeochemistry and Environmental Isotopes in Karst Springs of Makook Anticline, Kurdistan Region, Iraq

Karst springs of the Makook anticline were investigated to get more insight into the hydrogeologic, hydraulic, and hydrodynamic behavior of this complex karst system. Eight springs were studied in terms of hydrogeochemical constituents as well as δ 18 O and δ 2 The regional meteoric line plots between the global and Mediterranean meteoric water lines. The majority of the spring samples plot between the local and Mediterranean meteoric water lines implying influences by Mediterranean air masses, rapid infiltration of rainfall through the karst system and a short residence time (shallow karst aquifers). The correlation of d-excess and saturation index of halite was used to evaluate the evaporation process in karst waters of the area. The temporal variation in isotopic composition of karst springs was observed making it possible to distinguish between different origins for springs and indicate a possible connection between the aquifers and depending on this, the karst system can be classified into three aquifers: Behkme aquifer, Kometan aquifer, and Shiranish aquifer.


Introduction
Karst springs often discharge substantial amounts of water, which draws the interest of water resources researchers especially in arid and semi-arid areas facing water scarcity.Different methods and techniques were used by researchers to understand the response of the karst system to recharge events.Thermal patterns of karst springs were used by [1] as a key base for identifying two types of systems: thermally ineffective springs corresponding to rapid flow of event waters and effective thermal systems characterized by small fractures and slower water movements leading to equilibration of recharge and aquifer water's temperature.Electrical conductivity of karst waters may represent a good indicator to outline the flow system, whether it has a conduit or diffuse flow (retain piston flow).A piston flow phenomenon is the displacement of more mineralized water by intense rainfall recharge water toward the spring [2,3].The hydrodynamic and hydrogeochemical characteristics together can outline the behavior of karst springs [2].
Karst springs of the Makook anticline represent the main source of water for drinking, agriculture, and other uses in the areas of Balisan and Akoyan, Kurdistan Region, Iraq.The population of villages and towns in the part of Ranya district in which the studied area is located comprised more than 20,000 inhabitants.Water demand in the area has increased sharply in the last 20 years because of the increase in population and their needs.The mountainous area with poorly developed infrastructure creates certain limitations for studying the area.Few works have been done in the field of isotopes in the area in recent years.Only two local meteoric water lines (LMWL) were determined in nearby areas.The LMWL for Haji Omaran area was extrapolated by Mawlood [4] from isotopic data of Halab station in Syria and Adana station in Turkey, and the meteoric water line of Basara area was published by Hamamin and Ali [5].During 2000-2003 geologic and hydrologic investigations were carried out by a team from the Food and Agriculture Organization (FAO), including groundwater monitoring of Sarwchawa spring and wells in the area.Al Manmi [6] did an environmental isotopic study of surface and groundwater in Ranya area, which is close to the studied area and includes some of the springs considered in this work.According to [6], water of the Ranya area is of good quality with low salinity and thus suitable for drinking and irrigation purposes.Some possible recharge sources were determined depending on isotopic fingerprints.Karst aquifers in Makook system are mentioned as one thick Bekhme aquifer, which, according to [6,7], may contain Bekhme, Qamchuqa, and Kometan Formations.Bekhme aquifer consists of highly fissured, well-karstified carbonate rocks containing many caves and channels.
The objective of this study is to explore the applications of environmental isotopes of δ 18 O and δ 2 H in the field of karst springs.This work includes an attempt to build a local water line for the area, which had not done before.The local meteoric water line was correlated to the regional lines to investigate the source and determine the residence time of this water till it reaches the springs.This study also addresses the relation between hydrodynamic and hydrogeochemical characteristics of the karst springs.

Geography
The area of interest is located in the Kurdistan region in the northeast of Iraq, mainly in Sulaimani governorate, Ranya district (Figure 1).The area is situated within latitudes 36°10ʹ-36°35ʹ north and longitudes 44°30ʹ-44°50ʹ east in the elevation range of 500 to ˃2000 meter (m) above sea level (a.s.l).
Makook Anticline Mountain lies in a northwest-southeast direction and is surrounded by Shawre valley in the northeast and Balisan valley in the southwest.The area of interest comprises about 400 km 2 and is part of the Dokan lake catchment.The studied area is characterized by humid to moist climate.Rainy season starts from October and mainly ends in May (with minor showers in June, July, and September).The mean annual climatic parameters are presented in Figure 2 for the period 1980 to 2012.The highest precipitation rate, highest relative humidity, and lowest temperature are observed in February.The mean annual rainfall for the mentioned period is 570.5 mm and the average annual temperature is 20 °C.The mean annual evaporation (pan) is about 2000 mm and the mean relative humidity is about 47% [8].

Geologic Setting
Rocks of Mesozoic and Cenozoic age can be found in the region of interest (Figure 3).Two units of Jurassic rocks are cropping out, representing older rocks: (I) Sarki and Sehkaniyan formation, which are composed mainly of limestone, dark dolomite, and shale [9]; and (II) Sargelu formation (Bajocian-Bathonian), Naokelekan formation (Late Oxfordian-Early Kimmeridgian), Barsarin formation (Kimmeridgian-Early Tithonian), and Chia Gara formation (Late Tithonoan), which include bedded to massive dolomite, limestone, bituminous limestone, and marl.The Cretaceous unit comprises Balambo and Sarmord formation (marly and bedded limestone, bedded dolomite, and marl), Qamchuqa formation (mainly massive limestone and dolomite), Bekhme formation and locally Kometan formation (well bedded to massive limestone), Shiranish formation (well bedded limestone and marl), and Tanjero formation (sandstone, claystone, and conglomerate) [10].Kometan formation does not appear in the geological map (Figure 3), while it is more regional map and the existence of Kometan formation was in time of map creation still suspicious.Near Sarwchawa town, Kometan formation changed laterally to Bekhme formation through a transitional zone [11], and both formations are possibly overlain by Shiranish formation.Jurassic and Cretaceous rocks are occasionally overlain by Quaternary fan deposits (boulders, gravel, and fine clastics) and a thin layer of soil (mainly in the valleys, because of soil erosion in that mountainous area).
Makook Anticline represents a double plunging anticline (NW-SE) within in the parallel trend of Zagros folded structures between Ranya and Palawan anticlines and Sahwre and Balisan synclines (in between) in the NE and SW, respectively (Figure 3).The area is characterized by tectonic distortion, especially in the northwest part of Makook Anticline.

Figure 3.
Geological map of the area (modified from [12]).Numbers 1 to 8 refer to the springs where sampling was conducted.

Hydrogeology
Generally, all springs are within the Makook Anticline Mountain, except Chewa spring, which is located in the Balisan valley.Two main Karst aquifers are recognized in the region of interest.First, Bekhme karst aquifer supplies major springs along Makook anticline, including Zewa, Chewa, Bla, and Gullan springs (Figure 3).Bekhme aquifer is distinguished as a well-karstified, thick, semi-confined aquifer underlain by Sarmord formation (aquiclude).The second one is Kometan formation, which represents the source of Sarwchawa and Betwata springs.Kometan formation is a well-karstified, highly fissured aquifer with confined to semi-confined conditions and overlain by Shiranish formation [7].The aquifer is composed of carbonate rocks and contains a large amount of groundwater, varying in space and time.Strong karstification cycles took place in the aquifer system during the Paleocene-Miocene period [13].It is assumed that karstification happens recently.One of the springs (Shkarta) drains poorly developed karst (with small fissures).This spring is located in the Shiranish formation (marl and marly limestone), which was inconsistently recognized as a fissured aquitard (locally in the marly limestone part) by [6].The existence of Qamchuqa formation aquifer in studied area is possible, as mentioned by [6], but no field evidence has been observed.

Methodology
Eight karst springs (Sarwchawa, Shkarta, Betwata, Zewa, Chewa, Bla, Qala Saida, and Gullan) were monitored and sampled during six periods (September 2011, December 2011, April 2012, June 2012, September 2012 and November 2012).Precipitation samples (two snow and eight rain samples) were collected between November 2011 and April 2012 from three locations near Ranya town, Qala Saida and Gullan springs.Precipitation samples were collected by rain collector non-systematically depending on precipitation events.Field parameters (pH, conductivity and water temperature) were measured on-site (springs outlet).Redox potential (Eh) was measured via a WinLab redox meter (Windaus Labortechnik, Germany) and WTW SenTix-ORP electrode (WTW: a Xylem Group brand, Germany) and checked by standard redox buffer (pH = 7 and 220 millivolts at 25 °C).The parameters pH, electrical conductivity (EC), dissolved oxygen (DO) and water temperature (T) were measured on-site using a multi-parameter WTW model 3430, WTW pH electrode SenTix-940, WTW TetraCon-925 conductivity electrode and WTW FDO-925 optical dissolved oxygen sensor.Three-point calibration was carried out for the pH electrode via technical WTW buffers (pH = 4.01, pH = 7 and pH = 10).The EC electrode was checked with WTW standard control solution; settings were chosen so that the output of the electrode was converted to a water temperature of 25 °C.The optical sensor for dissolved oxygen measurement was checked regularly using a water vapor saturated air vessel.A certified mercury-thermometer was employed for air temperature measurement.Alkalinity was determined immediately after sampling by means titration and converted to hydrogen-carbonate ( HCO ) according to [14].Because no continuous discharge measurements of the springs were available, different methods (velocity-area using current meter and volumetric readings) were employed to record the discharge of the springs during the monitoring period from September 2011 to November 2012.Due to difficulties in the field, flow measurements and field parameters were measured only six times (September 2011, December 2011, April 2012, June 2012, September 2012 and November 2012) accompanying the sampling processes (Table S1).
All samples were filtered by 0.2 µm filters and stored at 4 °C until analysis in the laboratory.Samples analyzed by ICP-MS were acidified with an ultra-pure 1:1 nitric acid (30% HNO3) mixture after filtration in the field.Batches of 48 samples were analyzed in the laboratories of the Chair of Hydrogeology, Technische Universität Bergakademie Freiberg, Germany.Lab work included determination of Si and Sr ions with an ICP-MS XSeries-2 model (Thermo Scientific, Bremen, NS, Germany).Major cations Ca 2+ , Mg 2+ , Na + , K + were determined by ion chromatography (IC) using an 850 Professional IC Metrohm (Switzerland) and Metrosep C4-150 column with 2 mM dipicolinic acid eluent.Major and minor anions F¯, Cl¯, PO , SO , Br¯, and NO were determined by a Metrohm Compact IC Pro 881 (Switzerland) and Metrosep A sup 15-150 column with 3 mM NaHCO3 and 3.5 mM Na2CO3 eluent.The reproducibility of IC and ICP-MS determination were less than 2% and 5%, respectively.Stable isotopes of oxygen (δ 18 O) and deuterium (δ 2 H) were determined by means of an LGR liquid-water isotope analyzer-DLT 100 (ABB Company, Mountain View, CA, USA) with a precision of <0.1‰ for δ 18 O and <0.3‰ for δ 2 H.
Non-parametric two-tailed correlation analysis (Kendall and Spearman correlations) was performed for the hydrogeochemical and field parameters in addition to basic statistical treatment using SPSS software package.

Results and Discussion
Hydrologic and hydrochemical data of the Makook Karst springs are shown in the electronic Tables S1 and S2.The hydrochemical data were assessed by PHREEQC [15].The relative analytical error (E%) of the most samples was within an acceptable range of ±2%, except Shkarta and Zewa springs, which have a higher deviation in April 2012 and September 2012, respectively, but generally the accuracy is acceptable (Table S2).

Discharge and Recharge Estimations
From the discharge of springs, mean annual precipitation, catchment area, and average chloride concentration in precipitation and groundwater, different hydrologic parameters were estimated.Maximum discharge of karst springs has an important impact on hydrological modeling [16].The maximum estimated discharge (Dme) for an area is defined as: where P is mean annual precipitation in mm and Ac is catchment area in square meter (m 2 ).The mean annual precipitation is 570.5 mm (1.81 × 10 −5 L/m 2 /s) and the studied area is 4 × 10 8 m 2 ; therefore Dme of the studied karst area is 7240 L/s.The recharge in the area of interest was estimated using two methods: chloride mass-balance (CMB) and empirical Kessler method.The CMB method has often been used for recharge estimation [17,18] and is applicable in semi-arid regions when additional sources of chloride are absence [19,20].The precipitation in the area of interest was sampled 10 times during hydrological year (2011-2012).The temporal variation of chloride seems to be rather low (Table 1).Thus it is assumed that this average can be used for recharge calculation.Most of the chloride in the karst water is from precipitation (about 0.08 meq/L).Only 0.7% (0.02 meq/L Cl) of the chloride in karst waters stems from carbonate rocks, except in Shkarta spring where the chloride concentration is about 0.2 meq/L by the interaction of water with clay minerals of Shiranish aquifer (Tables 1 and S2).Most of the recharge areas of the springs are not arable land (except some areas related to Sarwchawa and Shkarta springs).Therefore fertilizers as additional source of chloride were excluded.Surface runoff in the recharge area of the springs is not significant and can thus be omitted from the recharge estimation.Recharge (R) in the area of interest was estimated by the CMB method: where Clp is average reference chloride concentration in precipitation and Clsp is average chloride concentration in spring water (both in mg/L).The re-calculation of Cl in precipitation to reference Cl is to consider the impact of evapotranspiration on the mass-balance estimation.The average reference Cl of precipitation is 2.52 mg/L and the average Cl in studied springs is equal to 2.9 mg/L, thus the estimated recharge is 496 mm/year.The recharge estimation by Kessler's method [21] was originally adapted for karst aquifers in semi-arid regions, and is based on the pluviometric data of last four months of the year.The estimated recharge in the area of interest by Kessler method is equal to 525 mm/year (the methodology of recharge calculation is described in [21]).The recharge rate according to Kessler's estimation is 92% and of CMB is 87%.The difference between the two estimations is normal (being different approaches were used), but the latter is more reasonable depending on the field observation, and conclusions of [22] in karst systems of other semi-arid region.These estimations are convenient with the calculated effective recharge (RE) that defined as: where RE is the effective recharge (the amount of water that reach the groundwater), Dsp is the total average discharge of the springs (4037 L/s), and Dme is maximum estimated discharge.The effective recharge in the studied area is 55.8%, and depending on these estimations it can be assumed that minority of the precipitation will loss by evaporation and transpiration.

Springs Hydro-Chemograph Analysis
Relating to the hydro-chemographs, different interpretations for the studied springs can be concluded (Figure 4).Flow rates in the springs (discharge) of the study area vary according to the capacity of the aquifer and degree of karstification in the spring aquifer (Figure 4a and Table S1).The Sarwchawa spring shows higher flow rates (2370 and 4630 L/s in dry and rainy seasons, respectively) and represents the most developed karst spring in the Makook karst system.The Shkarta spring (Shiranish formation) represents the less developed karst or non-karst spring with a low flow rate of 0.75 L/s, originating from fissured marly limestone [6].Sarwchawa, Shkarta and Qala Saida springs show a direct response of increasing discharge with respect to October 2011/April 2012 precipitation events (Figure 4a,b).In Shkarta spring, the response of discharge to the recharge events is less compared to the majority of the springs.As results of low flow rate in marly rocks, the residence time then will be high [23].Therefore, low flow rates in micro-fissures of Shiranish marly limestone can be the reason of high residence time and gradually varied discharge confirms this hypothesis.Betwata, Zewa, Chewa, Bla and Gullan springs show a nearly two months delay (maximum peak discharge appear in June 2012).The hydrograph of Betwata spring shows two peaks in response to the October-December 2011 and April-June 2012 rainfall events.A similar pattern can be noticed in Chewa spring's pattern but with lesser response to the October/December 2011 raining event.A gentle rise in Zewa spring's hydrograph was recorded and a steep decrease of discharge was recorded, which refers to a low storage capacity in the karst aquifer that feeds Zewa spring.Decrease in Bla's hydrograph in the December 2011 rainfall event was observed, which makes this spring different from all other studied springs.
The pH of the springs ranges between 7.1 and 8.4 (average 7.5) showing neutral to slightly alkaline water (Figure 4b and Table S1).The variability in pH of the springs is caused by changes in flow paths, water discharge, and residence times within the aquifer [24].Two of the springs (Zewa and Gullan) show decrease in pH during the December 2011 precipitation event.The similarity in pH patterns refers to the similarity in origin of water (same aquifer).The other springs show increasing of pH within the same recharge event (December 2011 precipitation).The addition of the newly recharged precipitation (December 2011-April 2012) caused an obvious decrease in pH of all the springs, and this was repeated in the November 2012 precipitation event.The pH value of precipitation collected during the current study indicates slightly acid conditions (Table 1).According to the nature of the karst aquifer, different impacts of this acidic recharged water are expected.Therefore, temporal variation in pH value can be used as a tracer in karst waters.
The redox potential of the springs showed partially oxidizing waters (323-478 mV, with an average of 383 mV).Chemistry of the karst waters in the unsaturated zone is modified by local flow variations and changes in redox condition [25].Analysis of the redox graphs reveals a falling pattern of Eh for all springs (Figure 4c).Variation of Eh graphs is controlled by changes in karst water temperature, and this relation is more significant in Shkarta, Chewa, Bla and Gullan springs.The increasing of Eh in Shkarta, Betwata, Chewa and Bla springs can be interpreted as a result of recharge with water containing more oxygen [26] in a shallower depth compared with the first group of springs.
Electrical conductivity ranges between 265 and 1117 μS/cm outlining differences in recharge rate and seasonal effects between wet and dry periods (Figure 4d).Generally, conductivity decreases in wet seasons, reflecting dilution by precipitation water.But this situation became completely the opposite (except Zewa and Gullan springs) during the November 2012 precipitation event when EC increases drastically during dilution by recharge.The dilution processes followed by a piston flow effect assumed to be the reason of the EC increase during the recharge period [3].The water temperature of the springs ranges between 6 °C and 23.4 °C (Table S1).The water temperature pattern in Shkarta spring was like that of Sarwchawa spring, but with higher amplitude in Shkarta (water temperature variation: 17.7 °C), which was recorded as the maximum among the springs studied (Table S1).A lateral groundwater influx from one aquifer to another with presence of clay and marl layers in the vicinity of springs leads to mixing and causes a decrease in water temperature [27].This scenario is the reason for the drastic decrease of spring water temperature (as obvious in the thermal pattern of Shkarta spring) in the rainfall period of October-November 2012 (Figure 4e).According to the stratigraphic superposition, a lateral influx of water from Kometan aquifer to Shiranish aquifer assumed to be the reason of the October-November 2012 temperature decrease in the Shkarta spring.The water temperature pattern in Betwata is the same as in the two former springs, but with nearly two-month time lag until recovery.The time lag in temperature recovery in karst springs results from low transit times, low development in karstification, and diffuse flow pattern [3].The variability of temperature in karst springs is considered a tracer for recharge and geometry of the aquifer [1].Slow movement of recharge water through Sarwchawa and Gullan spring aquifer results in thermal equilibrium between recharge and discharge water and indicates long residence time and rather deep karst.The enrichment of Sarwchawa spring water with Ca, Mg and SO4 indicates higher residence time, more intense water-rock interaction, and slower movement of recharged water through the aquifer.Gullan spring shows smoother thermal fluctuation at lower temperature (8 °C-14.1 °C) compared to Sarwchawa (11 °C-18.7 °C), which may be caused by differences in altitude of the springs outlet, elevation of recharge area, flow type, and depth of the aquifer (lower altitude, deeper aquifer and conduit flow for Sarwchawa spring).There is some evidence that convective heat transport of infiltrating water (recharge) is responsible for the observed temperature patterns [28].Unfortunately the data given are not sufficient to calculate the recharge by inverse modeling of the temperature distribution in the aquifer over time.
Temporal variation of major cations (Ca 2+ and Mg 2+ ) and anions (SO and HCO ) between wet and dry periods correlates to corresponding precipitation events (September 2011 to September 2012) and EC as illustrated in Figure 4.The calcium chemographs of spring water mainly show decreasing calcium concentration during recharge periods, except Zewa spring (Figure 4f), which shows the reverse (increasing Ca concentration with increasing recharge).Similar decreasing of Ca through dilution by low mineralized recharge water was reported by Vesper and White [29] for karst springs in Kentucky, USA.A zigzag-like pattern in Ca chemographs of Chewa and Bla springs was observed and indicates the similarity in recharges and flow systems for both springs.A time lag was enhanced by the most noticeable changes in hydrogeochemistry after the recharge events.This time lag suggests the presence of piston phenomena [29].
Chemograph of Si shows low variation, implying that silica in karst waters is conservative (Figure 4g).Shkarta spring shows higher Si concentration compared to the other springs, indicating presence of silicate mineral in this aquifer (Shiranish marly limestone).
Sulfate concentrations show the same variation as Ca corresponding to the dilution with recharge water, except in Zewa and Gullan springs, which showed the opposite patterns (Figure 4h).The sulfate concentration is higher in Sarwchawa spring.Sulfate-rich waters may form by oxidation of pyrite in carbonate rocks [30].Oxidation of pyrite concretions [9,31,32] and gypsum impurities along the flow path of Kometan aquifer leads to formation of this sulfate rich water.The direct response of sulfate to recharge is interpreted as presence of conduit type of flow and is reported from different karst springs worldwide [29,[33][34][35].
A general trend of increasing PCO2 with increasing discharge was observed in the karst springs.This reflects the decrease of diffuse atmospheric CO2 and increase of CO2 in soil [36].The pattern of PCO2 in karst springs is consistence with sulfate patterns rather than the Ca ones (Figure 4j).Increasing of PCO2 is accompanied by increasing of Ca, sulfate, EC, T and Eh.

Water-Rock Interaction and Geochemical Processes
The variation in hydrochemistry of water can be used as an indicator to understand water-rock interactions and geochemical processes during flow in the aquifers [37].In order to identify the hydrochemical type of the spring waters different techniques were used.Karst waters were classified depending on the relative abundance of Ca 2+ , Mg 2+ , Na + , and K + as major cations and HCO3 − , SO4 2− , Cl − , and NO3 − as major anions (in meq/L).According to the abundance, all springs in all sampling periods were classified as Ca-HCO3 water except Sarwchawa spring, which was classified as Ca-HCO3-SO4 water.The excess of Ca and SO4 in Sarwchawa spring relates to the longer residence time.A Chadha diagram classification [38] was also used for classifying the spring water, and all springs show Ca-Mg-dominant-HCO3 type (5th-sub field), which indicates that alkaline earths and weak anions exceed both alkali metals and strong anions.
Different hydrogeochemical relations indicate different sources of geologic material for the karst waters of the area (Figure 5).The significant relations of Na-Cl and Na-(NO3 + PO4) refers to different aquifers and sources of geochemical materials in addition to meteoric water (Figure 5a,b).Up to 0.017 meq/L of sodium stems from rainwater, while 0.093 of Na stems from limestone and dolomite of Kometan and Bekhme formation.Excess Na in karst water stems from clay minerals and feldspars [39].Therefore clay minerals of Shiranish marly limestone seem to be the source of excess Na in water of Shkarta spring.According to this variety of sources of geochemical materials, three different aquifers were recognized in the karst system of Makook anticline.(1) Springs feed by Shiranish aquifer (marly limestone) are characterized by high Na, Cl, NO3 and PO4; (2) springs feed by Kometan aquifer (limestone) shows medium content of Na, Cl, NO3 and PO4; and (3) springs feeds by Bekhme aquifer (dolomite) characterized by lower Na, Cl, NO3 and PO4 compared to the other aquifers (Figure 5a,b).
Calculated saturation indices (Table S3 and Figure 5c) indicate that Zewa, Chewa and Qala Saida springs are over-saturated (SI ≥ 0.05) with calcite while Sarwchawa, Shkarta, Betwata and Bla springs are in equilibrium (−0.05 ≤ SI ≤ 0.05).De-dolomitization processes are in progress in waters in equilibrium with calcite and under-saturated with dolomite [40].Only Gullan spring was under-saturated (SI ≤ −0.05) with respect to calcite, which suggests shorter flow path and residence time in the aquifer compared to the other springs (Figure 5c).Regarding the SIDolomite, only Qala saida spring was in equilibrium, and the rest were under-saturated with respect to dolomite.All springs were under-saturated with respect to gypsum, fluorite, and halite.Dissolution is assumed to be incongruent if a mineral dissolves while another one is oversaturated and thus precipitate.Obviously, incongruent dissolution occurs in waters of Sarwchawa, Shkarta, Betwata and Bla springs.SIGypsum is significantly correlated with SICelestine (Figure 5d) indicating presence of gypsum and celestine (SrSO4) impurities in karst rock, especially in springs drained from limestone.The sources of Sr are controlled by the stratigraphy [41] as well as the lithology.Thus, limestone of Kometan formation provides the karst water with more Sr than the two other aquifers (Figure 5d).  1.
A nonparametric Kruskal-Wallis-H test was performed to test the significance of the above grouping.The data were taken randomly and independently of each other.The results confirm that the three groups of spring water are different by means of p-value ≤ 0.05 for Ca, Mg, Na, K, Cl, SO4, HCO3, F, NO3, Si and Sr, except PO4, which shows higher p-values.

Development of a Local Meteoric Line
Oxygen (δ 18 O) and hydrogen (δ 2 H) isotope compositions for the precipitation and spring samples are given in Tables 1 and 2, respectively.The isotopic composition of precipitation in the studied area ranged between −12.28 to −3.35‰ for δ 18 O and −82.3 to −1.95‰ for δ 2 H.The regression equation of the LMWL is δ 2 H‰ = 20‰ + 7.9 × δ 18 O (R 2 = 0.95, p ≤ 0.001, and n = 10).In relation to the Global Meteoric Water Line (GMWL), given by equation δ 2 H = 10‰ + 8 × δ 18 O [42], and the Mediterranean Meteoric Water Line (MMWL), given by equation δ 2 H‰ = 22‰ + 8 × δ 18 O [43], the LMWL is located between these two (Figure 6), and it is closer to MMWL.The LMWL of the studied area is very close to that of Haji Omaran area (50 km northeast of the studied area, 36°40ʹ00ʺN, 45°03ʹ00ʺE and with higher altitude) δ 2 H‰ = 20‰ + 8 × δ 18 O [4], but far away from that of Basara area (100 km southeast of the studied area, 35°20ʹ30ʺN, 44°57ʹ40ʺE and with lower altitude) δ 2 H‰ =14.4‰ + 7.7 × δ 18 O [5].The d-excess higher than 20‰ in precipitation results from arid vapor sources [44].The vapor sources of some precipitation samples (R-2, R-3, R-7, and R-9) are arid (low humidity), because the d-excess is higher than 20‰.The evaporation of precipitation can reduce the d-excess [44].Therefore, the d-excess in samples R-1, R-4, R-5, R-6, R-8 and R-10 was reduced as a result of water loss by evaporation (Table 1).A continental effect is a possible reason for the high d-excess in some precipitation samples (R-3, R-7 and R-9), because they were higher than that of the MMWL.Another possible reason for this high d-excess in R-7 and R-9 is the origin and nascence condition [45], because they represent snow samples.

Factors Affecting Stable Isotopes in Precipitation
The air temperature is correlated with d-excess (d-excess ‰ = 31.9‰− 1.4‰/°C × TAir (°C), R 2 = 0.67, p = 0.03, n = 10).This correlation may emphasize the seasonal effects and thermodynamic conditions in the atmosphere [45,46] (Figure 7a) and may also be underlined by the increase of evaporation with temperature (Evaporation (mm) = 0.26 mm/°C × TAir (°C) -0.44 mm, R 2 = 0.6, p = 0.05, n = 10).The amount effect (Figures 8 and 9b) shows a negative correlation (higher rain amount shows a lower isotopic signature), which was reported worldwide except in polar areas [47].StD = Analytical standard deviation; NA = Not available; Spring codes can be used to trace the locations of the springs in Figure 1.

Isotopic Pattern of the Karst Springs
The spring water samples have isotopic compositions between −8.29‰ and −6.14‰ for δ 18 O and −45.02‰ to −31.89‰ for δ 2 H (Table 2).Two regression lines for the karst springs were developed (Figure 6), one for all spring samples, excluding those of November 2012 period (δ 2 H‰ = 2.3‰ + 5.56 × δ 18 O, R 2 = 0.95, p ˂ 0.001 and n = 40), and the other for November 2012 samples (δ 2 H‰ = 5.56 × δ 18 O − 1.36‰, R 2 = 0.99, p ˂ 0.001 and n = 8).The majority of the samples plot between LMWL and MMWL.From this it can be concluded that the precipitation that recharges the karst aquifers originates from Mediterranean vapor masses, reflecting rapid infiltration of rainfall through the karst system and a short transit time (the karst aquifer is shallow, therefore, short residence time does not lead to reliable change in isotopic composition of the recharged water).Shkarta spring samples are isotopically enriched compared to the other springs.In argillaceous carbonates, evaporation leads to enrichment of isotopes in the pore water [48].The isotopic enrichment of Shkarta spring relates to slow infiltration because of low permeability in the marly limestone of Shiranish aquifer and thus increased evaporation.The first regression line of the springs (green dashed line in Figure 6) is assumed to be a local evaporation line (LEL), which may occur due to the influence of Shkarta spring.The relation between δ 2 H and chloride concentrations of the springs (δ 2 H‰ = 21.83‰+ 0.473 × Cl, R 2 = 0.9, and p ˂ 0.001) supports the evaporation trend in karst springs of the area (Figure 8a).Minerals like halite (NaCl) react fast with water attaining equilibrium within a short time compared to the residence time of groundwater [49].Saturation of halite during evaporation of aqueous solution is accompanied by δ 18 O enrichment [50].A significant correlation between the d-excess and SIHalite is notice (Figure 8b).This plot is interesting and supports the idea of evaporation, because it shows the same evaporation trend of δ 2 H-Cl.According to this correlation (Figure 8a), increasing of SIHalite corresponds to decrease in d-excess.The correlation between SIHalite and d-excess is regarded to the effect of Mediterranean origin of the precipitation in the area of interest.Therefore, the saturation index of halite could represent a tool to evaluate the evaporation process in environmental isotopes of karst waters.The trend of evaporation is increasing from Zewa spring water, which is less enriched with deuterium and chloride, in comparison to Shkarta, which is the most enriched spring (Figure 8).

Altitude Effect
The isotopic composition of the groundwater in relation to altitude represents an indicator for locating the groundwater recharge area [51,52].The plots of δ 18 O and δ 2 H in the spring water of the area versus the corresponding altitude of the springs are presented in Figure 9.The altitude effect on groundwater of the area is expressed by the linear correlation of δ 18 O and δ 2 H with altitude: where A is altitude in meter, R 2 is equal to 0.5 and 0.4, p-value is equal to 0.04 and 0.08 and n is 48 for δ 18 O and δ 2 H, respectively.The local isotopic gradient (LIGO&D) shows depletion of stable isotopes by about −0.17‰ and −0.8‰ per 100 m elevation for δ 18 O and δ 2 H, respectively (Figure 9).The most depleted springs are Zewa, Bla, and Gullan, while Shkarta is the least depleted one (lower elevation).

Temporal Variation
The temporal variation in the isotopic composition of the springs is obvious in all springs (Figure 10).This fluctuation shows more or less the different fingerprints of dry and wet periods.The variation in isotopic composition of karst springs indicates different origins of recharged water and possible connection between the aquifers [53].The gap between Shkarta and other springs (Figure 10) indicates a connection between Bekhme and Kometan aquifers.Different isotopic signatures for the springs were recognized, separating them into three aquifers: Shiranish aquifer, represented by Shkarta spring (the most enriched isotopes), Kometan aquifer (Sarwchawa and Betwata), mainly represents intermediate isotopic signatures and Behkme aquifer (Zewa, Chewa, Bla, and Gullan springs), generally represented by depleted isotopic signatures if compared with the above two aquifers.The δ 2 H is more sensitive to the temporal variation than the oxygen isotope one, because the mean difference ranges from 0.17 to 0.61‰ for δ 18 O (about eight times the analytical standard deviation of δ 18 O, 0.07‰) and from 3.15 to 4.5‰ for δ 2 H (about 11 times, 0.4‰).Depletion in δ 2 H isotopic signature of the springs was noticed in samples of November 2012 (Figure 7) and this may be due to the nearly lowest recorded temperature of the spring waters during this sampling event and the amount of winter precipitation.This depletion-temperature dependency was also indicated in Figure 7 (lower water temperature characterized by lower δ 18 O).Depletion was also recorded in Chewa spring for both δ 18 O and δ 2 H April 2012, which probably occurred by the effect of precipitation amount of rainy season (January to April).

Conclusions
The main objective of this study was to understand the isotopic and hydrogeochemical relations in a complex karst system of Makook anticline.EC, pH, Eh, and water temperature temporal variations can be used as tracer for residence time, aquifer storage capacity and depth of the karst aquifers.Hydrogeochemical classification and isotopic composition can be used as indicators for the mean residence time in karst aquifers.The sources of water constituents and driving forces of hydrogeochemical processes are magnified by interrelations of karst water constituents.The sources of the geochemical properties were also indicated by these relations.
The local meteoric water line was derived for the area and is plotting between the global and Mediterranean meteoric water line (however, much closer to the Mediterranean MWL).From the d-excess in precipitation, the vapor source was found to be an arid, low humidity, and evaporation-dominant Mediterranean type.Temperature and amount effects are obvious from the isotopic behavior of the precipitation of the area.From the isotopic plot of springs it can be concluded that the recharged precipitation is of Mediterranean origin and that it infiltrates rapidly.Enrichment in stable isotopes gives rise to the hydraulic and hydrodynamic properties of karst aquifers.The local evaporation line is obvious from the oxygen-deuterium plot from Zewa to Shkarta spring and this trend was confirmed by the chloride-deuterium and d-excess-SIHalite relations.It can be concluded from the temporal variation of stable isotopes that connections between the aquifers exist.

Figure 1 .
Figure 1.Digital elevation model (DEM) showing the area of interest in Kurdistan Region, Iraq.Numbers 1 to 8 refer to the springs where sampling was conducted.

Figure 2 .
Figure 2. Climatological data of the Dokan Meteorological Station (690 m a.s.l) located in the investigation area.Averages for the period 1980-2012 are shown [8].

Figure 6 .
Figure 6.Meteoric water line (LMWL) and springs isotopic line for the studied area.Blue line is the global meteoric water line (GMWL) from [42]; black line is the Mediterranean meteoric water line (MMWL) derived by [43] and the red line is the local meteoric water line (LMWL) from the present study.

Figure 7 .
Figure 7. Relation between meteorological and isotopic parameters.(a) Plot of d-excess and temperature in precipitation and springs water, data of springs are represents mean d-excess and water temperature of six sampling period for each spring; (b) Plot of event precipitation amount and δ 2 H in precipitation.

Figure 8 .
Figure 8. Illustration of evaporation trend from springs of Bekhme aquifer to Shiranish aquifer spring (Shkarta spring).(a) δ 2 H-chloride plot; (b) d-excess and SIHalite plot; values of δ 2 H, d-excess, SIHalite and chloride represent mean values of six samples for each spring: September 2011to September 2012.

Figure 9 .
Figure 9. Relation between altitude and δ 18 O and δ 2 H in the studied area.(a) Plot of δ 18 O in spring water and altitude of the spring outlet and (b) for δ 2 H. LIGO is local isotopic gradient of δ 18 O and LIGH is for δ 2 H.

Figure 10 .
Figure 10.Temporal variation in isotopic composition of the karst springs.(a) Variation of δ 18 O during the period September 2011 and November 2012; (b) Variation of δ 2 H during the period September 2011 and November 2012.

Table 2 .
Isotopic composition of spring water samples.