Hydrogeologic Investigations in Western Crete by Using Isotopic Analyses and GIS Techniques ()
1. Introduction
During the last decades groundwater systems have been thoroughly investigated to acquire the necessary data for undertaking protection measures against pollution pressures, especially in aquifers destined for public water supply [1] . Nevertheless, in many cases groundwater resources cannot be managed efficiently since the required information is either inadequate or absent which leads to gaps in understanding and designing the appropriate restoration practices.
Isotopic techniques can offer useful information in groundwater investigation efforts [2] -[5] and therefore can facilitate optimum water management decisions [6] [7] . Important hydrogeological information can be provided from stable isotopes such as the groundwater recharge zones, the interrelations between groundwater and surface waters, detection of the water velocity, age and the direction of groundwater flow. Also, basic hydraulic characteristics of aquifers such as transmissivity, porosity and dispersivity can be determined with the contribution of isotopes [2] . This information is valuable for the qualitative and quantitative assessment and management of groundwater resources, especially in aquifers where monitoring efforts of classic hydraulic parameters are limited or missing [8] -[10] . Moreover, environmental isotopes provide also useful information in modeling ground- water systems [11] . According to Koniger [12] , isotope studies led to a better understanding of the movement and mixing of the recharge procedures in the Upper Rhine valley aquifer in an area where classic techniques did not provide adequate information about the recharge and origin of groundwater.
Moreover, heavy stable isotopes such as oxygen-18 and deuterium and the radioactive isotopes like carbon-14 and Tritium are involved in groundwater investigations [13] . According to Gibson [14] , the isotopes (1H, 2H, 16O, 18O) are widely used as tracers in the hydrogeological research. Groundwater studies have been also based on isotopic tracers for identifying interactions between the groundwater bodies and the surrounding springs [14] . Specifically, stable isotopes can give excellent information about the water circulation while the radioactive isotopes can identify the residence time if the groundwater is not contaminated by other water bodies [13] .
Although isotopic measurements can give information about the movement of water and recharge in groundwater systems, they can also be used for identifying others useful factors. Oxygen isotopes can be used for the separation of groundwater from the runoff water in surface water bodies [15] . Changes in isotopic composition are caused by chemical interactions related to ground weathering, by evaporation processes and by mixing of water. Isotopic deviations from the precipitation to the groundwater composition show an extended residence time and a mixing of different sources of water [15] .
There are many different approaches to estimate recharge zones, including multicriteria analysis of hydrogeological properties, remote sensing, numerical modeling and geochemical investigations [16] -[18] . Nevertheless, there are significant limitations in these approaches, including the intense monitoring efforts required to acquire the necessary information, the time and costs involved as well as the lack of physical interpretation in some of the approaches. On the contrary, several research efforts have been conducted successfully, at a global scale, to delineate the recharge zones of important aquifers by using Oxygen isotopes and hydrogeologic data [18] -[20] . For this purpose, the isotopic signatures of groundwater and rainfall are usually used that reflect the discharge and recharge characteristics of the study catchment [21] . Particularly, Clark and Fritz [15] state that: “the groundwater isotopicsignature reflects a long term average of rain events while the average precipitation signature over a rain event has a discrete concentration that falls within the range of the mean annual concentration”.
The purpose of this study is to identify the relationships between the main groundwater bodies and the springs of Western Crete, estimate the recharge zones of each group of springs while propose management and protection measures for their long term sustainable use.
2. Methodology
2.1. Study Area
The study area covers the main springs of Western Crete which are used for domestic and irrigation purposes (Figure 1, Table 1). The total water demand in the study area (Prefecture of Chania) is 130 × 106 m3 on an annual basis and approximately 85% of this amount is allocated to agriculture and 13% to domestic uses [22] . The most important springs in terms of discharge and water use are Agyia (domestic supply in the city of Chania), Stylos (large scale bottling factory), and Meskla (Table 1). No regular monitoring of the springs’ discharge occurs in the particular area but some indicative measurements from the Water Resources Department of the Regional Authority of Crete illustrate Agya, Almyros, Armenoi and Stylos as the most important springs in terms of discharge (from 63 to 80 × 106 m3, Table 1). The geological substrate of the study area (Figure 1) comprises calcareous rocks (the White Mountains massif, approx. 42% of the study area), Flysch formations (western part, approx.: 27% of the area) and Neogene and Quaternary sediments (lowland, south part, approx.: 31% of the study area).All the springs under study are located within calcareous areas (limestones) with the exception of the
Figure 1. Study area and sampling sites in Western Crete.
Table 1. Location and characteristics of the studied springs.
lowland springs (Koleni, Stylos and Armenoi) that are found within alluvial deposits, overlaying usually calcareous formations (Figure 2).
2.2. Sampling Campaign and Analyses
A series of water samples for isotopic analysis has been collected over the period 2005-2007, from 15 different
Figure 2. Main geological formations and the average δ18O values.
sites (11 springs and 4 rain samples, Figure 2) at a seasonal basis (10 sampling campaigns). The samples were collected in plastic bottles of 100 ml, following the zero-head air principle and transferred to Demokritus Research Center for isotopic analysis in δ18O and Tritium (TU). The equipment used for the analyses of δ18O samples is a SIRA Series III sotope Ratio Mass Spectrometer while for the Tritium a Packard-Tricarb LSA Liquid scintillator has been used. The Electrical Conductivity (in μS/cm) of the water samples were measured in-situ with a portable Hanna Inc. Instrument. Moreover, a water quality campaign in July 2006, took place in the same sites to identify the qualitative status of the springs and propose specific management measures.
The Oxygen isotopic composition of the water samples were the main focus of this effort, which is defined as delta (δ) values, expressed in the per mille deviation from the standard isotopic composition known as SMOW (Standard Mean Ocean Water, [15] [23] . The isotopic ratio of 18O/16O is defined in the following equation:
‰
If the value of δ is negative then the sample is “depleted” in the light O isotope (low amount of evaporation and high recharge altitude) while if it is positive the sample is “enriched” in comparison to the standard which means lower recharge altitude and more exposed to the evaporation effect.
All the isotopic values together with the additionally collected data (Tritium, EC, Water quality parameters) underwent a thorough statistical analysis including descriptive statistical estimations, hierarchical clustering and regression analysis. Particularly, the δO18 values have been used to identify the groups of springs that have similar behavior and may be affected from the same aquifer. For this process the hierarchical clustering algorithm as well as a box plot analysis in SPSS software has been used. A pair-wise regression analysis followed between the stable isotope concentrations and the physicochemical parameters in order to identify common relationships and trends. Moreover, a correlation matrix has been also estimated to identify links between the major ions concentrations and the stable isotopes levels in all the springs. These results have been combined to interpret the interrelationships between sampling points and identify the common recharge areas. The rainfall samples collected from 4 different locations in Western Crete on January 2006 (single rainfall event), covering the local altitude natural variations, were used to estimate the relationship between altitude and δ18O values by using the least squared equation approach. The resulted equation has been then used to estimate the recharge altitude of the sampling sites. The average recharge zone altitude for each sampling site was estimated by using the following the following equation:
where y = the estimated average recharge zone altitude (m), a = the slope of the least squared regression line connecting rainfall δ18O values and altitude, x = the sampling points average δ18O value and b = the offset in y axis of the least squared regression line mentioned above. The particular recharge zone altitude estimations is an approximation of the reality since the available data to justify the accuracy of these estimations are limited but combined with other sources of information such as the geological and topographic substrates can provide the basis for the recharge zones identification.
The clustering exercise used the hierarchical clustering algorithm in R software to identify group of sites with similar behaviour as far as δ18O values are concerned. Regression analysis between δ18O and TU, Electrical Conductivity and Ca were used to validate the aforementioned groups. Then, the particularly outcome was combined with the recharge zone altitude estimations and the main geological properties of the area to identify common recharge areas and the potential limits of the respective aquifers.
A water quality survey was attempted in all the sites at July 2006 and the acquired water samples were analysed in the hydrochemistry laboratory of the Hellenic Centre for Marine Research. The examined parameters included: SiO2, Ca2+, Mg2+, Na+, , and. The correlation matrix of these chemical analyses, including the same dates’ δ18O and Tritium values has been estimated to depict the potential relationships between them while an assessment of the water quality status followed to reassure the classification of springs in groups based on the δ18O values.
The potential interactions of the examined springs became apparent through the aforementioned approach and a qualitative assessment followed in all of the studied aquifers to present specific management and protection measures.
3. Results
3.1. δ18O Fluctuations
The results indicated that the monitored springs illustrate a large variation in their δ18O values both in seasonal and interannual basis (Table 2). Particularly, the highest average δ18O value was detected in the Agya spring with the value of −8.15‰ followed by the Armenoi spring (−8.10‰). The lowest δ18O value was observed in the Koleni spring (−4.81‰) followed by the spring named Argyroupoli with the value of −4.83‰ (Table 2). Furthermore, Argyroupoli and Almyros had a low standard deviation (0.46 and 0.47 respectively) indicating that the intrannual fluctuations are limited while Ag. Eirini and Fournes 2 had a high standard deviation (0.77 and 0.67 respectively) value showing that the δ18O values throughout the year were fluctuating significantly. Also, the skewness was calculated (Table 2) which describes the degree of asymmetry of δ18O values’ distribution around its mean [24] . The majority of the springs had a positive value of skewness and therefore most of their δ18O values were lower than their means. Fournes 2 was the only spring which its skewness was negative showing that most of its values were higher than the mean. Moreover, Ag. Eirini and Meskla (followed by Stylos and Armenoi) that present high StD and Skewness are probably affected more from snowmelt than the other sites. This is because snowmelt in high altitude has very low δ18O values while rainfall in moderate and low altitudes present relatively high δ18O values, especially in Crete. Thus, the aforementioned sampling sites undergo the impact of extremely low δ18O values more frequently than the other sites.
The seasonal fluctuation of δ18O values (Figure 2 and Figure 3) indicates that the highest values were observed in October 2006 and July 2006, while the lowest ones in January 2006 and November 2005. The springs such as Koleni, Almyros, Argyroupoli, Meskla, Agiya, followed more closely this particular seasonal pattern than the rest of the sampling sites. Armenoi and Stylos illustrated their highest values in April and November 2006 while their lowest values were the same of the other sites. This behavior could be related to the main water supply mechanism that changes throughout the year. Particularly, for the sampling sites that have recharge zones affected by snowmelt mostly the spring and summer periods are the ones with the highest δ18O values (Armenoi
Table 2. Descriptive statistics for the δ18O values in the studied sampling sites.
Figure 3. The δ18O values fluctuations in the 11 springs from November 2005 to May 2007.
and Stylos) due to the characteristic δ18O values of snowmelt and the extended travel time that can fluctuate from 2 - 5 months. The sites that are both affected by snowmelt and rainfall have higher seasonal fluctuations but always in the lower levels of δ18O values (Meskla, Agya, Ag. Eirini) while the ones that are located in lower altitudes and are mostly affected by rainfall have significant seasonal fluctuations at a lower δ18O level (Almyros, Koleni, Argyroupoli).
The similarity index of the cluster analysis indicated that Stylos and Armenoi are very close not only in absolute values but also in their seasonal and interannual fluctuations. The same stands for Ag. Eirini and Sempronas 2 while Meskla and Fournes 2 illustrate the most dissimilar behavior in relation to the all other sites.
Thus, according to the seasonal fluctuations and the cluster analysis two main group of springs can be observed from the data; the first one includes Agiya, Amenoi, Stylos and Meskla while the second one Koleni, Argyroupoli, Almyros and Potistiri (Figure 4).
Figure 4. The cluster analysis of the δO18 values for the springs.
The Box and Whisker plot of all the data (Figure 5) shows the substantially different behavior of the various sampling sites as far as δ18O is concerned. The majority of springs did not have a normal distribution apart from the springs named Ag. Eirini and Fournes 2 which presented a distribution close to normal. The highest median of δΟ18 values appeared in Stulos spring and the lowest one in Koleni (Figure 5). The maximum range of the δ18Ο values was found in Fournes 2 while the minimum one in Almyros. The clustering groups of sampling sites described above are also accredited by the Box and Whisker plot since Stulos and Armenoi have very similar distribution of values while the same stands for Almyros, Argyroupoli and Potistiri, as well as for Agyia and Meskla and to a lesser extend for Ag. Eirini and Sempronas (Figure 5).
3.2. Correlation between δ18O and Other Parameters
The δ18O values of the sampling sites have been correlated with other measured parameters such as Electrical Conductivity (EC), Tritium and Ca to identify potential relationships between these parameters and distinguish group of sites with similarities in their behavior.
Regarding EC, Almyros spring had the greatest values, ranging from 2909 μS/cm to 2180 μS/cm which was expected due to its proximity to the sea and the dominant geological substrate that includes limestones with fluctuating amounts of gypsum (Figure 6). The average conductivity in the rest sampling sites fluctuated between 633.40 μS/cm and 252.20 μS/cm which are typical values for karstic groundwater with relatively limited residence time.
The least squared regression between EC and δ18O indicates a quite low correlation coefficient (0.364) between the two parameters which is mainly due to Almyros site that appears as an outlier in the analysis. According to their scatterplot (Figure 6), there is a positive linear relationship between conductivity and δ18O which groups the sites in a similar way as the previously mentioned analyses.
The Tritium values show a slightly positive linear relationship with the δ18O values but again the correlation coefficient is low (0.1531) showing that this relationship is not statistically significant (Figure 7). The greatest range of Tritium values was observed in Fournes 2 with the maximum and minimum ones to be 16.7 TU and 4.2 TU respectively (Figure 8). The average Tritium values of the rest of the springs fluctuated between 5.66 TU and 2.76 TU while Sempronas 2 had the lowest average Tritium value (2.79 TU). In terms of groundwater age, based on the Tritium concentrations, Ag. Eirini and Sempronas have relatively low values which indicates a mixture of recent and old (prior to 1952) water [15] while the rest of the sampling sites have Tritium values from 4 to 10.5 TU which indicates modern water. The regression analysis between Tritium and δ18O values (Figure 7) indicates a clear grouping of the samping sites (group 1: Stylos, Armenoi, Agyia, Meskla, group 2: Almyros, Potistiri, Koleni and Argyroupoli, group 3: Ag. Eirini and Sempronas).
The relationship between Ca and δ18O values in the particular sampling sites were weak but reassured some of the aforementioned groups (Figure 9). Almyrosand Koleni springs had the highest Ca values (107 mg/l and 61mg/l respectively) while Sempronas 2 and Ag. Eirini had the lowest ones (20 mg/l and 35 mg/l respectively). Stulos, Armenoi and Agyia presented very similar Ca values while Argyroupoli and Koleni had almost identical
Figure 5. The boxplots of the δΟ18 values in the different springs.
Figure 6. The scatterplot of the conductivity against δΟ18.
Figure 7. The scatterplot of the Tritium against δO18.
Figure 8. The temporal variations of the Tritium values in the study area.
Figure 9. The scatterplot of the Ca against δO18.
concentrations (Figure 9). The different concentrations of Ca in the studied springs express the high variation of the geological formations and the groundwater residence times in the particular area. Springs with relatively low Ca concentrations such as Sempronas and Ag. Eirini indicate probably aquifers with limited calcareous content and/or short travel times. Proximity in the sea also affects the Ca concentrations significantly, such as the case of Almyros spring.
3.3. Recharge Zones
The rainfall sampling sites and their δ18O values have been used to estimate the relationship between δ18O and recharge altitude in the study area (Figure 10). The coefficient of determination for the identified relationship was very high (0.84) and therefore the associated regression equation could be used to estimate the potential, average recharge zone altitude for all the sampling sites (Table 3).
Stulos spring had the greatest average recharge altitude (1040 m) followed by Armenoi and Agiya springs
Figure 10. Regression analysis between δ18O and recharge altitude.
Table 3. Estimated recharge altitude for all the sampling sites and dates.
(1020 m). Koleni and Argyroupoli had the lowest average recharge altitude (460 and 480 m respectively), followed by Almyros (560 m) and Potistiri (570 m, Table 3). The potential seasonal fluctuations in the recharge zone altitudes of each site indicate that Ag. Eirini and Meskla have the highest values (range of 610 and 620 m respectively), while Fournes and Potistiri have the lower variations (350 and 380 m respectively, Table 3).
The scatterplot of recharge altitude against δ18O values (Figure 11), illustrates again the strong negative linear relationship between these parameters while re-validates the grouping effects mentioned before in the particular study.
By combining the sampling sites clustering and topographic information together with their recharge zones altitude and the dominant geological conditions, an effort has been made to delineate the recharge zone boundaries for each site as well as the common aquifers shared by each group of sites (Figure 12). These recharge zones present only the potential areas of water supply for each group of springs, based on the fluctuations of the average recharge altitudes for each sample and the dominant geological substrate. Nevertheless, more work needs to be done to validate and possibly amend these areas by using more isotopic data and reference values to optimize and/or accredit the relationship between recharge altitude and isotopic values.
4. Water Quality Results
The main pollution pressures in the study area are agrochemicals and untreated urban waste from villages and
Figure 11. The scatterplot of the recharge altitude against δ18O.
Figure 12. Potential recharge zones of the main springs of Western Crete.
touristic infrastructure. Therefore, if the pollution pressures are exceeding the carrying capacity of the hydrogeological systems in terms of natural attenuation, elevated levels of NO3, SO4 and Na are expected. However, the water quality analyses (Table 4) indicates that there is no evidence of significant pollution impacts in the particular springs with the exception of Potistiri site where NO3 reaches 15 mg/l. In the rest of the study sites the NO3 concentrations were less than 4 mg/l and SO4 levels were below 35 mg/l with the exception of Almyros which is influenced by the sea (brackish waters).
The correlation matrix (Table 5) indicates that Ca, Na, Mg, SO4 and Conductivity have very high correlation coefficients which is due to their common origin in this area (either from rocks weathering or sea intrusion/ spray effect) while the rest of the parameters do not illustrate significant relationships amongst them. Intermediate correlation coefficient have been estimated between HCO3 and δ18O (0.422) as well as between NO3 and δ18O (0.395, Table 5). This indicates that in many cases, the more positive the δ18O, the highest the HCO3 and NO3 values in groundwater. This is an expected trend that has to do with the groundwater residence times and the density of the pollution sources that is higher in the lowlands (where most of the agricultural land and touristic areas are located).
The variations in the water quality parameters accredit the grouping of sites occurred based on the δ18O values (Table 4). Armenoi and Stylos illustrate common Ca (40 and 44 mg/l respectively) and Mg (6.7 and 6.8 mg/l respectively) values and the same stands for Argyroupoli and Potistiri (Ca: 57 and 56 mg/l respectively while Mg: 10 and 9.8 mg/l). Almyros is clearly influenced by the sea (Na: 778 mg/l and SO4: 313 mg/l) due to its close proximity to the shore and therefore the water quality parameters cannot provide appropriate information to justify its prior grouping based on the δ18O values. Moreover, Ag. Eirini and Sempronas 2 present very similar Mg (5 and 3 mg/l respectively) and SiO2 (5.66 and 5.6 mg/l respectively) levels but have some differences in
Table 4. Water quality analysis and isotopic levels for the study sites.
Table 5. Correlation matrix for water quality and isotopic parameters for the study sites.
the Ca concentrations (35 and 20 mg/l respectively, Table 4). This could be due to the slightly different travel time of the water discharges in these 2 springs.
5. Discussion and Conclusions
The use of stable isotopes as tracers for the study of groundwater movement and interactions between different aquifers has both advantages and disadvantages. The most important advantage is that the isotopes concentration in groundwater cannot be influenced by the potential reactions between rocks and the groundwater, under normal temperatures. There is also a relationship between the isotopes of oxygen and altitude, latitude and season which allows for the identification of potential recharge zones for springs. Areas with higher latitude are characterized by more negative δ18O values than the sites with lower latitude [15] .
Nevertheless, the use of isotopes for the study of hydrogeological groundwater characteristics has also disadvantages. The expensive cost of the isotopic analyses is a reason for the limited use of them in the groundwater research. A great number of samples must be collected and a significant number of isotopic analyses are required for a complete study of a groundwater system [25] .
Particularly, the present isotope study in Western Crete led to a better understanding of the movement and mixing procedures of the different aquifers and their relationships with the main springs of the area. Springs that are fed from common recharge zones have been identified and thus areas of special concern for groundwater management and restoration have been illustrated. Based on the potential recharge zones altitude, the examined springs have been classified in 4 different groups:
1st group: Stylos, Armenoi, Agyia, Meskla-recharge alt.: ~700 - 1300 m (Meskla presents a bit lower rescharge altitude (500 - 1100 m).
2nd group: Almyros, Potistiri, Argyroupoli,-recharge alt.: ~250 m - 750 m (Argyroupoli has a bit lower value in recharge alt. and its UT is very low which means mixing with relatively old groundwater)
3rd group: Ag. Eirini, Sempronas 2-recharge alt. ~400 - 950 m
4th group: Koleni has probably a local recharge area which agrees with the clastic sediments found south of the spring.
The western part of Lefka Ori mountain aquifer (karstic) is discharged through Ag. Eirini, Sempronas, Meskla, Agyia (this spring is probably recharged by the central limestone formation of Lefka Ori which is in higher altitude than the surrounding limestone and dolomitic formations). Stylos and Armenoi discharge the central part of Lefka Ori (limestones and dolomites) while Almyros, Argyroupoli and Potistiri discharge the eastern part of the aforementioned mountain which has relatively lower altitudes (limestone and dolomitic formations). Based on the aforementioned clustering of the studied springs, the fluctuations in the altitude of the potential recharge zones and the dominant hydrogeological conditions of the study area, the potential boundaries of the springs’ aquifers have been determined. Particularly, Potistiri and Argyroupoli seem to share a common aquifer while Almyros has its own local recharge zone in the eastern part of the studied area, partially affected by the seawater. Ag. Eirini and Sempronas have also a common aquifer in the center of the study area while Meskla and Agyia partially share a high altitude aquifer in the northern part of the White Mountains. Stylos and Armenoi have also a quite extensive common aquifer in the northeast part of the White Mountains while Koleni and Fournes seem to have smaller, lower altitude and local recharge zones.
The water quality data available accredited the initial grouping of springs based on δ18O values and indicated significant variations in the pollution pressures of the studied springs. Potistiri NO3 values indicate the existence of substantial pollution pressures, possibly from agricultural activities, at a local level while the rest of the sites have quite low NO3 levels (less than 4 mg/l). The springs with the highest recharge altitude (Armenoi, Stylos and Meskla) illustrate the lowest SO4 values (4 to 5 mg/l) while there is a general trend for SO4 to increase as recharge altitude decreases.
Therefore, the springs that have large recharge zones in the higher altitudes and thus higher discharges and good quality waters (Stylos, Armenoi, Meskla, Agyia) should be intensively protected. This could be achieved by avoiding the installation of hazardous activities in the particular significant recharge zones and by imposing specific restrictions in the amount of water abstractions in these areas. The area covering the Potistiri recharge zone should be examined in detail for pollution pressures and immediate restoration measures should be initiated since the current levels of NO3 are significantly elevated and alarming. Nevertheless, further investigations should be done to calibrate and validate the estimated recharge zones of the main springs in the broader area of Chania. The sampling campaign should include more seasons and many more rainfall isotopic analyses in order to accredit the relationship between the recharge altitude and the δ18O values of the springs. However, the present approach in the particular study area provided important information for the main aquifers’ and springs’ properties which can be used for the better management and protection of these, important for humans and the environment aquifers.
Acknowledgements
Special gratitude to Dr. Nikolaos Zouridakis, the responsible researcher of the isotopic hydrology laboratory of the Demokritos Research Center where most of the isotopic analyses took place.
NOTES
*Corresponding author.