Hellenic Karst Aquifers Vulnerability Approach Using Factor Analysis : The Example of the Louros Karst Aquifers

The karst aquifers of the Hellenic territory display high vulnerability that is related to the existence of the epikarst zone. If present, the epikarst zone regulates the vertical movement of the groundwater, which results in the chemical pollutants staying longer within it, thus facilitating the inactivation of the microbial loads. The Hellenic karst is well developed; it displays high hydraulic conductivity, relatively large storage, and the characteristics of an important water reservoir. Although its porosity was calculated to be between 0.005–8%, the statistical processing of the springs’ discharge time series points to an increased storage. The increased storage was attributed to the presence of the epikarst zone, which regulates the water movement to the underlying formations, and finally to its discharge from the karst system through springs. In these areas, the influence of nitrogen compounds on the water quality is evident. The presence of cyanobacteria and toxic substances at concentrations that are too high for fresh, flowing waters whose temperature does not exceed 17 ◦C advocates the role of nutrients, even at concentrations that do not exceed the maximum permissible limits of drinking water. In mountainous areas, where agriculture is limited, the presence of these compounds is attributed to heavily polluting land uses. The impact of such uses became apparent from the application of factor analysis, which not only highlighted their contribution to the overall groundwater quality, but also determined their geographical distribution and origin. From the above, it becomes evident that there is a need for protection of the karst aquifers, and investigation of their vulnerability mechanisms. The Hellenic karst, which is expected to cover the drinking water needs in the near future, proves to be very vulnerable to any kind of pollution.


Introduction
In previous years, the drinking water and irrigation needs were covered by the easily accessible aquifers of granular geological formations.The relative isolation of carbonated aquifers in mountainous areas has contributed to maintaining the view of their good water quality characteristics.However, the water quality degradation of the granular aquifers, which was mainly due to seawater intrusion caused by over-pumping, has turned the interest to the study of carbonate formations, and highlighted the impact of anthropogenic activities on water pollution.Nowadays, the eutrophic conditions and the certain substances that have been observed in lakes' water are due to their pollution by human activities.
Although the large water quantities that are commonly circulating through the karst formations help in the dilution of the pollutants, the statistical processing of the chemical analyses reveals the domination of the nitrogen compounds, and thus the presence of eutrophication.
Karst aquifers are characterized by high groundwater velocities and low storativity.These two properties contribute to the immediate restoration of the quality of the polluted aquifers, provided that the pollution sources are eliminated.However, recent research has shown that the presence of the epikarst zone has a positive effect on the function of these aquifers, as it slows down the vertical movement of the water therein, and creates storage conditions.However, precisely due to its presence, the view of rapid restoration needs to be revised.
In the present study, aiming at the qualitative investigation of the vulnerability of karst aquifers: (a) the basic characteristics of the Hellenic karst are given where similarities among the karst formations in relation to their hydraulic properties and the presence of the epikarst zone can be noticed; (b) the role of the epikarstic zone as the regulator of the vertical movement of water is highlighted; and (c) Louros Basin, which is a typical example of one of the most significant and representative karst systems in the country, in terms of both its geological structure and its land uses, is discussed.
In order to investigate the origin of pollution, hydrochemical data has been processed with the use of factor analysis.The proposed methodology can be useful to study other karst aquifers that have similar characteristics.

Methodology
Along the Louros watershed, five samplings were carried out between October 2008 and May 2010 from springs, boreholes, and surface waters.In the present study, the data from the 27-28 February 2010 sampling-which was considered to be the most representative-was elaborated using factor analysis.The samples were collected from 65 sampling points according to the United States (US) Environmental Protection Agency (EPA) (1976) procedures [1].Seawater and rainwater samples were also collected [2,3].
The water temperature, pH, dissolved oxygen, electrical conductivity, and redox potential were measured in situ using a Hana ® HI 9828 portable instrument.Additionally, alkalinity was measured at the site using a Hach ® Digital Titrator.
Three different water samples in acid-washed, well-rinsed polyethylene bottles for subsequent analyses were collected from each sampling site.The first 1-L volume contained a bulk water sample for the anion analysis.The second 0.1-L volume contained water filtered through a Whatman 0.45-µm cellulose membrane and acidified with 0.5 mL of ultrapure HNO 3 for the determination of metals.The third contained bulk sample for the isotopic analyses.
All of the chemical analyses were performed in the Laboratory of Hydrogeology, University of Patras within hours from collection.Anions (NO 3 − , NO 2 − , SO 4 2− , and F − ), P, NH 4 + , and SiO 2 concentrations were measured with a Hach ® DR 4000 spectrophotometer.The Cl − content was determined using AgNO 3 0.1 N titration.
Major cation (Ca 2+ , K + , Mg 2+ , Na + ) concentrations were determined using a GBC ® Avanta P flame atomic absorption spectrophotometer.The charge balance error for all of the chemical analyses is within a highly acceptable range (±5%).Chemical calculations such as the saturation indices were performed using PHREEQC version 3.0 [4].
Along the Louros Basin, there is a large number of livestock, poultry farms, and fish farms, as well as inhabited areas.In order to detect the possible pollution of the groundwater resulting from these uses, samplings and microbiological analyses were also carried out in February 2010, which is the month with the highest water levels, on 46 surface and groundwater samples.The samples were stored in sterile glass bottles and analyzed for microorganisms, total coliforms, Escherichia coli, Enterococcus, sulfite-reducing clostridia, and the total number of microorganisms at 22 • C and 37 • C in the Public Health Laboratory at the Department of Medicine of the University of Patras.

The Hellenic Karst
A large part of the Greek territory, exceeding 35%, is predominantly covered by Mesozoic carbonate formations.The Hellenic karst systems were developed during the Alpine orogeny that took place between Jurassic and Miocene [5].Lateral movements in the upper Eocene (e.g., the Pindos napple, the Pelagonian napple), intense tectonism at the end of the upper Miocene, as well as the eustatic movements of the Mediterranean Sea during the glacial and interglacial periods created the necessary conditions for the karstification of the carbonate formations of the different geotectonic zones.
According to [6], in the Hellenic karst, the following types can be distinguished: a.
The "paleokarst" that was developed during the upper Cretaceous and upper-lower Miocene under warm-wet climates (Savanna type), which is not currently evolving and is considered as inactive, fossilized.b.
The "younger karst" that was developed during Pleistocene-Holocene between the interglacial periods, under climatic conditions that are usually warmer and wetter than today.The main karstification process of the Greek territory took place during the same period, and was also accompanied by the creation of new karst forms, and sometimes the reactivation of older ones that often continue to evolve, even to this day.c.
The "recent karst", which is created under the current climatic conditions, and is characterized by surface karst features, e.g., karrens.
The large karst basins, which were already formed in the Miocene, constitute a well-developed paleokarst contrary to the next generation of smaller-sized basins that began to form during Plio-Pleistocene and lower Pleistocene [7].
The majority of the available data for the Hellenic karst originates from the IGME (Institute of Geology & Mineral Exploration) surveys.The most well-known studies focus on the central Peloponnese and Western Greece, where the carbonate formations of the Ionian, Gavrovo-Tripolis, and Olonos-Pindos zones are mainly developed.Additionally, individual studies provide important information for other areas, for example [8] for the Pre-Apulian zone, [9,10] for the Parnassos-Giona zone, [11] for the Pelagonian and Attico-Cycladic zones, [12] for the Koziakas Mountain, [13] for the crystalline limestones of the Rhodope massif, and others [14].
The main development of the karst, in its holokastric form, is observed in the lithological homogeneous thick bedded limestones of the Cretaceous-Senonian age of the Ionian, Gavrovo-Tripolis, Parnassos-Giona, and Pelagonian zones (Figure 1).The thin-bedded limestones, particularly due to the presence of radiolarian cherts such as those of the Ionian zone (Vigla), Pindos zone, etc., exhibit a lower karstification degree.According to [15], they form karst systems that are classified as merokarst.
In the western and central Greece and in the Peloponnese, where holokarst formations are found in the Ionian, Parnassos-Giona, Pelagonian, and Gavrovo-Tripolis Zones, respectively, all of the exokarst and endokarst features such as dolines, poljes, and caves are present.
In merokarst, the movement of water is carried out through cracks and fractures.In the case of the coexistence of two karst structures, the strong differentiation in hydraulic conductivity forces the latter to behave as relatively impervious to the former, e.g., Senonian and Vigla limestones of the Ionian zone [17].
In the central Peloponnese, the lithological variations of these structures, as well as their stratigraphic position-the merokarst limestones of the Pindos zone overlap with the holokarst formations of Gavrovo-Tripolis zone-lead to the development of sinkholes in the thick-bedded, and lithologically homogeneous limestones of the Gavrovo-Tripolis zone [18].

The Role and Presence of the Epikarstic Zone in General, and in the Hellenic Karst
The importance of the epikarst, as a zone of particular characteristics, has emerged from studies of various researchers since the 1970s.The identification of aquatic fauna at the cave's ceiling pre-required the existence of a saturated zone in between the ground and the ceiling of the caves.
Mangin [20,21] studied karst springs' hydrographs, and concluded that there should be some kind of water storage occurring on the upper part of the unsaturated zone.Bakalowicz et al. [22] observed variations in the hydrochemical and isotopic composition of groundwater before it reaches the aquifer, which was attributed to water storage on the upper part of the unsaturated zone.Kiraly and Morel [23], in order to obtain a typical hydrograph form, during simulations of karst aquifers, highlighted the need for water storage on the upper part of the unsaturated zone.The speleologists [24,25], during their investigations in central Asian mountains, concluded that the presence of water in the unsaturated zone after a long period of drought can only be attributed to water storage on the upper part of this zone.This resulted in the presence of a zone near the soil surface that behaves as a recharge zone for the underlying karst system.Mangin [20,21] named this zone as "epikarst", and pointed to the presence of a perched aquifer within it.
The great significance of epikarst, in relation to the propagation of pollutants in karst aquifers, was highlighted by White and White [26] for its role in the dissemination of NAPLs.The epikarst often acts as a reservoir that delays the karst recharge [21,[27][28][29][30][31].
The presence of the epikarst zone was investigated by many researchers through the statistical processing of springs' discharge data.The correlogram describes the memory of the karst aquifers."Memory time" is the time needed for the original conditions to no longer affect the functioning of the aquifer, that is, the aquifer may "cease to remember the original conditions".It gives information on the degree of karstification as well as on the water reserves [32-34], since long memory corresponds to high porosity and small development, and vice versa.
The research on the karst formations that is known to have taken place in Greece does not refer to the epikarst zone.The published hydrogeological data mainly concern the determination of hydraulic properties using the standard methodology of pumping tests, wherefrom very low storage values were obtained.On the contrary, the exploration of different karst formations throughout the country (Figure 2), and the statistical analysis of the springs' discharge data, pointed out that in every case, there seems to be a factor that regulates the groundwater discharge that would correspond to increased porosity.In particular, the following results were derived: Panagopoulos and Lambrakis [35] comparatively studied the characteristics of hydrographs of the daily springs' discharge values of these karst aquifers, and found that they display long memory (lag time >50 days), which corresponds to increased porosity and direct aquifers' response to the rainfall, as well as instant drainage and high velocity of groundwater, i.e., increased hydraulic conductivity.Similar conclusions were also drawn for the Almyros aquifer by [36], who associated the increased memory to an increased porosity of the formations.However, since the measured porosity is very small, the increased memory could be attributed to the presence of the epikarst zone.
(b) The karst system of the Louros Basin consists of carbonate formations of the Ionian zone.From the springs' discharge data for the Eocene limestones of the area, a long memory emerged just as in the previous case (lag time between 70-100 days) corresponding to increased porosity, which is explained by the presence of the epikarst zone [2,37].It was also found that the very good drainage of the aquifers corresponds to a developed karst system.(c) The karst system hosted in the marbles of Aggitis presented similar results regarding the autocorrelation and cross-correlation functions, and thus similar characteristics with the above-mentioned systems [13].The lag time of 50 days for the autocorrelation function was attributed to the system's long memory, corresponding to an increased storage, which could explain the huge water reserves that were calculated by the recharge area of the springs.Given the very low porosity that was measured in the area, the increased memory was attributed to the presence of snow that accumulates on mountain tops and can play a regulatory role in the discharge of the springs' groundwater.According to [14], the increased memory could also be attributed to the presence of the epikarst zone.
Geosciences 2018, 8, x FOR PEER REVIEW 2 of 14 huge water reserves that were calculated by the recharge area of the springs.Given the very low porosity that was measured in the area, the increased memory was attributed to the presence of snow that accumulates on mountain tops and can play a regulatory role in the discharge of the springs' groundwater.According to [14], the increased memory could also be attributed to the presence of the epikarst zone.

Geological and Hydrogeological Setting of Louros Basin
The study area receives increased precipitation rates throughout most of the year (mean annual precipitation around 1700 mm).The extent of the basin is relatively limited, and its morphology does not display extensive flattened areas and exokarstic forms.However, the occurrence of karst features and the presence of large areas covered by terra rossa are common.
It is hosted in the formations of the Ionian zone [39].The oldest ones consist of evaporitic layers of early to middle Triassic age (Figure 3), locally overlain by breccias [40].
These formations are overlain by a thick sequence of carbonate and clastic rocks reflecting a continuous sedimentation from late Triassic to upper Eocene [41].The Pantokrator and Upper Senonian limestones are considered to be permeable, and their storage coefficient values range between 2-3% [42,43].Clastic limestones of Eocene and Vigla limestones are considered to be of low permeability.

Geological and Hydrogeological Setting of Louros Basin
The study area receives increased precipitation rates throughout most of the year (mean annual precipitation around 1700 mm).The extent of the basin is relatively limited, and its morphology does not display extensive flattened areas and exokarstic forms.However, the occurrence of karst features and the presence of large areas covered by terra rossa are common.
It is hosted in the formations of the Ionian zone [39].The oldest ones consist of evaporitic layers of early to middle Triassic age (Figure 3), locally overlain by breccias [40].
These formations are overlain by a thick sequence of carbonate and clastic rocks reflecting a continuous sedimentation from late Triassic to upper Eocene [41].The Pantokrator and Upper Senonian limestones are considered to be permeable, and their storage coefficient values range between 2-3% [42,43].Clastic limestones of Eocene and Vigla limestones are considered to be of low permeability.The latter formation [17], when it underlies limestones of high permeability such as the Senonian ones, constitutes the base level of the local aquifer, under which karstification is very limited or absent.Dolomites, shales, Ammonitico Rosso, and flysch formations are considered to be of very low permeability.Finally, the Neocene and Quaternary deposits also display low permeability compared to the carbonate rocks.The major tectonic features of the region are folds with their axes trending SW-NE and NNW-SSE.The broader area is characterized by a continuous succession of synclines and anticlines with a general direction of NNW to SSE.Such synclines, filled with flysch, limit the area to the east and west and control the movement of water.
Geosciences 2018, 8, x FOR PEER REVIEW 3 of 14 The latter formation [17], when it underlies limestones of high permeability such as the Senonian ones, constitutes the base level of the local aquifer, under which karstification is very limited or absent.Dolomites, shales, Ammonitico Rosso, and flysch formations are considered to be of very low permeability.Finally, the Neocene and Quaternary deposits also display low permeability compared to the carbonate rocks.The major tectonic features of the region are folds with their axes trending SW-NE and NNW-SSE.The broader area is characterized by a continuous succession of synclines and anticlines with a general direction of NNW to SSE.Such synclines, filled with flysch, limit the area to the east and west and control the movement of water.Strike-slip fault systems display the same directions while younger ones, normal to the previous ones, create a tiered structure and contribute to the formation of individual karst units.Finally, diapirism resulted in the slip surfaces of gypsum, which facilitated the over-thrusting, contributing to the complex tectonic structure, and affecting groundwater quality.
The most important aquifers of the area are hosted in the Pantokrator and Upper Senonian limestones.They display high permeability due to intense karstification and fracture porosity.Individual karst units, which are more or less interconnected, form a karst system that is drained by Strike-slip fault systems display the same directions while younger ones, normal to the previous ones, create a tiered structure and contribute to the formation of individual karst units.Finally, diapirism resulted in the slip surfaces of gypsum, which facilitated the over-thrusting, contributing to the complex tectonic structure, and affecting groundwater quality.
The most important aquifers of the area are hosted in the Pantokrator and Upper Senonian limestones.They display high permeability due to intense karstification and fracture porosity.
Individual karst units, which are more or less interconnected, form a karst system that is drained by the Louros River.These units are discharged by a number of karst springs that are distributed at three distinct levels where the three aquifers are developed [45].

The Pollution Sources and Origin of Groundwater in Louros Basin
Groundwater samples of the Louros Basin fall under the category of oligometalic carbonate waters of Ca-HCO 3 type.They are freshwaters that are affected by the diapirism of the evaporites in the central part of the basin, while they are affected in the southern part and the lowlands by the presence of hypothermal water [46].The electrical conductivity (EC) ranges between 147-923 µS/cm.The higher values correspond to the samples from the southern part of the basin, while the average value above that part is 598 µS/cm.Nitrate concentrations range between two and 11 mg/L.Katsanou et al. [46] using the chemical components and isotopic signatures of groundwater, reported that the Louros karst system is recharged by rainwater, while it is isolated from adjacent ones, including the Ioannina Basin, beyond the northern margins of the Louros Basin.Thus, both nutrient and microbial loads are produced within the basin, and are attributed to human activities.
In the autumn months, when the temperature is higher, there is a thick layer of cyanobacteria on the surface of Lake Ziros.The water surface of the lake is a continuation of the surface of the aquifer hosted in the carbonated rocks of the wider area, and is being discharged through springs along its margins.
According to [47], cyanobacterial blooms is a phenomenon that occurs in Lake Ziros due to the presence of the Planktothrix rubescens species.At the same time, microcystin concentrations were determined to be the highest that were ever identified in Greek ecosystems.The level of the microcystins that were measured in cyanobacterial pellets were also found to be among the highest in the world.According to [48], microcystins are toxic substances that could cause serious damage to human health through drinking water, while studies at the University of Oregon showed that they grow in stagnant fresh waters under high temperature and nutrient supply [49].
Although the density of the sampling points is not considered satisfactory due to objective difficulties, i.e., the lack of sampling points in mountainous karst areas, microbial determinations revealed the presence of microbial load in the majority of samples.From the distribution of microbial parameters (Figure 4), it appears that non-pathogenic microorganisms are found throughout the basin with significant concentrations as a result of point and not diffuse sources.

The Application of R-Mode Factor Analysis in Louros Basin
The application of factor analysis on chemical data can play a key role in hydrochemical research.Several methods such as principal component analysis, cluster analysis, discriminant analysis, and factor analysis are available [50].The R-mode factor analysis is applied to investigate the relationships between a large number of hydrochemical parameters (15) and samples (65), aiming to create a much smaller number of factors that will express without loss-if possible, or with the least possible loss of statistical information, the groundwater quality [51].The theoretical part is omitted, while SPSS 24.0 [52] was used to process the data.In the chemical analyses of the groundwater samples that were collected from the area, both Rmode and Q-mode factor analyses were applied.For the R-mode, the first four factors were selected on the basis of the criterion that the values of the coefficients are greater than or equal to one (Table 1, column 1).The overall variation amounted to 76.8% of the total statistical information (Table 1, column 3).
The first factor, which accounts for about 40% of the total variance, presents high loadings on major ions, HCO3 − , K + , Na + , Mg 2+ , NH3, P, and Cl − (Table 2), which suggests that the chemical composition of water is mainly formed by these ions.This relationship is valid for the majority of In the chemical analyses of the groundwater samples that were collected from the area, both R-mode and Q-mode factor analyses were applied.For the R-mode, the first four factors were selected on the basis of the criterion that the values of the coefficients are greater than or equal to one (Table 1, column 1).The overall variation amounted to 76.8% of the total statistical information (Table 1, column 3).
The first factor, which accounts for about 40% of the total variance, presents high loadings on major ions, HCO 3 − , K + , Na + , Mg 2+ , NH 3 , P, and Cl − (Table 2), which suggests that the chemical composition of water is mainly formed by these ions.This relationship is valid for the majority of natural waters [50].The presence of phosphorus, which is not only linked to human activities such as NH 3 , but might also be related to the presence of phosphorites in the area, is highlighted by the first factor.This factor indicates the significant contribution of P to the groundwater's chemical composition, which increases the nutrients' content, and leads to the development of eutrophication.
The second factor, with a total variation of 15%, exhibits high loadings for Ca 2+ and SO 4 2 ions, clearly indicating the effect of gypsum on the groundwater quality of the area [46].The third factor highlights the relationship between NO 2 − concentrations and temperature.The highest NO 2 − concentrations are abundant in the central part of the area (Figure 5a).In this part, the water level is very close to the ground surface, and the groundwater is influenced by the atmospheric temperature.
The fourth factor links the NO 3 − concentrations to pH.The highest concentrations of nitrates are related to the inhabited areas (Figure 5b).The last two factors, with a cumulative variance of about 20% of the total variance, demonstrate the effect of nitrogen compounds on the chemical composition of groundwater, as well as the presence of pollution.
The application of the Q-mode factor analysis led to the extraction of the factors' scores that in turn led to the mapping of the sampling positions and the respective contribution of each of them to the presence of these factors.Figure 5 shows the distribution maps of the third and the fourth factor that are related to pollution processes.Figure 5a presents the distribution of the third factor, which has to do with the presence of nitrites.It appears that the higher factor scores occur in the central and eastern segments of the basin where most agglomerations and human activities are concentrated.Figure 5b is completely analogous to the previous one.

Discussion and Conclusions
The inner continental karst basins are promising water reservoirs for the near future.They are mainly known for: (a) their inaccessibility due to their steep relief; (b) their isolation from the effect of the sea, since they are mostly remote from coastal areas; (c) their cold climate and their high precipitation rates; (d) the large springs that they host, and (e) the limited human activities that they accumulate.
The Hellenic karst aquifer systems present both similarities and differences due to their geological structure and karstification degree.The similarities include good drainage conditions, low storage, and long memory, which corresponds to a regulating factor of the groundwater discharge of these formations that could be attributed to the presence of an epikarst zone.
Although the combination of large quantities of water ensures the dilution of pollutants and the inaccessibility of the area limits farming, and hence the use of fertilizers, it appears that land uses such as the inhabited areas and the poultry farms are sufficient to create widespread pressures and the degradation of water quality.
The concentrations of nitrates determined in surface and groundwater samples by [46] did not exceed 25 mg/L, but along with other nutrients, they are considered to be responsible for the presence of cyanobacteria.
Although karst formations are ultimately susceptible to pollution, their susceptibility is very difficult to calculate.Chemical pollutants are favored by the presence of the epikarst zone because during polluting episodes, they accumulate within it, and are gradually released during the first rains of the next wet period to the underlying karst aquifers.Microbial loads propagate directly during wet periods mainly through surface drainage.Therefore, the need to protect these aquifers, and study their pollution susceptibility, becomes imperative.

Discussion and Conclusions
The inner continental karst basins are promising water reservoirs for the near future.They are mainly known for: (a) their inaccessibility due to their steep relief; (b) their isolation from the effect of the sea, since they are mostly remote from coastal areas; (c) their cold climate and their high precipitation rates; (d) the large springs that they host, and (e) the limited human activities that they accumulate.
The Hellenic karst aquifer systems present both similarities and differences due to their geological structure and karstification degree.The similarities include good drainage conditions, low storage, and long memory, which corresponds to a regulating factor of the groundwater discharge of these formations that could be attributed to the presence of an epikarst zone.
Although the combination of large quantities of water ensures the dilution of pollutants and the inaccessibility of the area limits farming, and hence the use of fertilizers, it appears that land uses such as the inhabited areas and the poultry farms are sufficient to create widespread pressures and the degradation of water quality.
The concentrations of nitrates determined in surface and groundwater samples by [46] did not exceed 25 mg/L, but along with other nutrients, they are considered to be responsible for the presence of cyanobacteria.
Although karst formations are ultimately susceptible to pollution, their susceptibility is very difficult to calculate.Chemical pollutants are favored by the presence of the epikarst zone because during polluting episodes, they accumulate within it, and are gradually released during the first rains of the next wet period to the underlying karst aquifers.Microbial loads propagate directly during wet periods mainly through surface drainage.Therefore, the need to protect these aquifers, and study their pollution susceptibility, becomes imperative.
Funding: This research received no external funding.
(a) The coastal karst aquifers of Trifylia consist of Eocene carbonate formations, while the corresponding areas of Almyros spring in Heraklion area consist of the Triassic carbonate rocks of the Ionian zone and Cretaceous limestones of the Gavrovo-Tripolis zone.

Figure 4 .
Figure 4.The pie charts of THE microbial load along Louros Basin.

Figure 4 .
Figure 4.The pie charts of THE microbial load along Louros Basin.

Figure 5 .
Figure 5. Geographical distribution of a. the third, and b. the fourth factor in Louros Basin.

Figure 5 .
Figure 5. Geographical distribution of a. the third, and b. the fourth factor in Louros Basin.

Table 1 .
Rotation Sums of Squared Loadings.