Hydrogeology and Hydrogeochemistry of the Lauria Mountains Northern Sector Groundwater Resources (Basilicata, Italy)

In this study, the hydrogeological characterization of the northern sector of the Lauria Mounts carbonate hydrostructure (southern Apennines, Basilicata region) has been carried out and the hydrochemical properties of different collected groundwater samples have been characterized. Several normal springs drain the hydrostructure, some of them characterized by high annual mean discharges. Groundwater samples were collected from different springs; many parameters such as pH, electrical conductivity, and total dissolved solids have been measured, and major (cations and anions) elements and stable isotopes have been analysed following standard test procedures. Other chemical characteristics were derived from the analysed quality parameters. The results elucidate that the main hydrogeochemical processes control the chemical content and assess the quality of the groundwater within the hydrostructure. The analyses highlight that the chemical compositions of groundwater are strongly influenced by the lithology, especially limestones and dolomitic limestones; they explain and confirm the hydrogeological setting of the system. The groundwater system displays light different geochemical signatures. The processes contributing to the concentrations of major ions depend primarily on carbonate dissolution. The analysis, in all studied groundwater samples, shows that the facies groundwater type is Ca–HCO3, bicarbonate is the dominant anion, and calcium is the dominant cation with appreciable magnesium concentrations. To identify the aquifer’s recharge areas, the environmental stable isotopes oxygen and hydrogen, deuterium, and O were analysed. The unaltered δO and δD signatures for the groundwater of the major springs allows identifying the recharge area of these emergencies at elevations ranging from 900m to 1000m (a.s.l.), pointing out the presence of deeper flow regime feeding of these springs. The groundwater sample isotopic characteristics of D and O suggest that most of the groundwater is recharged directly by infiltration in a high-permeability medium.


Introduction
In Mediterranean countries, water availability is a key factor for social and economic development. Recent changes in demography and lifestyle, agricultural land use and irrigation, tourism development, and climate change projections point to an increase in water scarcity and environmental problems with negative implications towards current and future sustainability [1][2][3]. The balance between water demand and water availability has reached critical levels in many regions; for this reason, a sustainable approach to water resource management is required. This is able to consider the interactions among climate change impacts and water scarcity, surface water and groundwater pollution, and water engineering and human systems including societal resilience and adaptations [4,5]. Groundwater represents important freshwater resources, used as drinking water supply for agricultural and industrial purposes. Climate variability and human change affect groundwater systems contributing to the decrease in water availability, both directly through modifications in aquifer recharge and indirectly through changes in groundwater quality and use [6]. According to the European and Italian Laws [7][8][9] on the protection of groundwater against deterioration and pollution, it is necessary to apply the most appropriate actions and procedures at the hydrogeological basin scale. These actions, finalized to define strategic sustainable management processes, must be able to prevent or contain the potential depletion and pollution of groundwater resources, also taking into account mitigation strategies for climate change impacts [10,11].
These precious resources, widely exploited, are not distributed at regional level and within each country, in a homogeneous way. The increase in groundwater exploitation poses a severe risk to the availability of the water resources, and the resulting resource scarcity is a major concern in most countries of the Mediterranean region.
Thus, in the last decades, many researches and studies have been planned, devoted to the assessment of the groundwater quantity, quality, and trends. The hydrogeological and groundwater hydrogeochemical assessments have a primary importance to define the hydrogeological conditions of the systems for implementing a sustainable integrated management of groundwater.
Hydrochemical characteristics and stable isotope compositions of groundwater provide important data regarding the water-rock interaction along flow paths, groundwater mixing due to different aquifers and hydrological processes, and area and timing of groundwater recharge for the comprehensive understanding of the hydrodynamic and hydrogeological setting [12][13][14][15]. In particular, stable isotopes of oxygen and hydrogen (δ18O and δD) are useful as space tracers of hydrological processes in aquifers, to determine the elevation and conditions of the recharge areas [16][17][18].
In the Basilicata region (southern Italy), extensive karst and fractured carbonate hydrostructures due to their peculiar geostructural features and hydrogeological conditions are characterized by the presence of a considerable number of springs of good quality.
In order to contribute to an improved understanding of the hydrogeological setting and the groundwater hydrogeochemical characterization of the still poorly understood carbonate hydrostructure of the Lauria Mountains northern sector (southern Apennines), a comprehensive methodological approach was applied. Investigations included geological and hydrogeological field studies to understand the specific characteristics of the hydrogeological system, groundwater sampling at the springs, and chemical analysis of the major elements and stable isotopes (δ18O and δD) to define the hydrogeochemical processes and the recharge area of the major springs.
This will help to provide insight and support in defining groundwater appropriate development actions and strategies, into comprehensive, sustainable, and efficient water resource assessment and management to overcome the future and increase freshwater demands, during the drought periods and in a changing climate scenario.

Study Area
The morphostructure of the northern sector of the Lauria Mounts, geographically located along the Calabrian-Lucanian regional border, in the southern-western part of the Basilicata region, is oriented in the W-E direction from Lauria to Castelluccio villages. It is bordered by the watershed of the Sinni River to the N, by the Torrente Peschiera to the NE, and by the urban areas of Castelluccio Superiore and Castelluccio Inferiore to the E. The Mercure River basin bounds the carbonate massif in the SE sector, the Valico Prestieri and Fosso Mancosa in the southern sector, Serra La Nocara to the SW by the Lauria urban area, and the lower part of the Torrente Caffaro tributary of the Noce River to the W. The morphostructure includes the reliefs of La Spina (1652 m) and Zaccana (1580 m) Mounts, Castello Starsia (1387 m), Serra Rotonda (1285 m), and Serra Tornesiello (1185 m) ( Figure 1).
The geographic location of the study area, immediately inside the imposing Apennine Chain, which rises along the Tyrrhenian coast, has a strong impact on climate factors. The nearness to the sea confers its beneficial effect to the climate which results to be substantially mild, despite the elevations of the territory, with a very high average annual precipitation of about 1550 mm. The annual typical Mediterranean rainfall regime involves very rainy winters and dry summers. The maximum monthly precipitation value is recorded in December, the minimum in July. The regime of temperature is essentially characterized by the mean maximum value of about 23°C in summer that always occurs in the months of July-August and the mean minimum value of about 6.5°C in winter. The average annual temperature is instead about 14°C. The hydrographic network is mainly present on the carbonate rocks that form the central structure and on the alluvial deposits characterized by good permeability and low slope gradients. The development of minor streams is controlled by the state of fracturing of the carbonate and is distributed around several directions.
The straight segments of streams (e.g., Torrente Caffaro, Fosso Mancosa) are strongly related to the main fault systems that evidently have conditioned the hydrographic network development. In the SE edge of Serra Rotonda, a small lake, Lago La Rotonda, with a seasonal regime is located in a small endorheic basin; it has a water-extended surface in rainy months and dry in the summer. The lake is generated in correspondence of a karst polje, in which, at the bottom, detritical and palustrine deposits are present. The morphostructure hosts an important carbonate hydrostructure, drained by many springs some of which are characterized by huge groundwater discharges.

Geological Setting of the Lauria Mounts
The Lauria Mountains consist of a series of ridges, coincident with morphostructures with a N120°SE pattern, located to the west of the Mercure Basin, which separates them from the neighbouring mountains of the Pollino Chain and represents the link between the Campanian-Lucanian Apennine and the domains belonging to the crystalline-metamorphic sedimentary rocks of the Calabrian-Peloritan Arc.
The Lauria Mountains, characterized by a significant geological and structural complexity, represent a high structure consisting of the Meso-Cenozoic calcareous-dolomitic succession, referable to the Alburno-Cervati-Pollino Unit 2 Geofluids [19], constituting the monoclines of the Calabrian-Lucanian border, confined by the Miocene argillite-marly flysch.
In particular, in the Lauria Mountains it is possible to distinguish the Lauria-Castelluccio northern ridge, consisting of the main reliefs of Castello Starsia, Serra Rotonda, La Spina, and Zaccana, where carbonates, from Triassic to Cretaceous, such as dolomites, dolomitic limestones, and limestones constitute the central structure, while terrigenous sediments are present along the lower part of the northern slopes. Quaternary sediments organized in different continental depositional environments fill the southern sector. This morphostructure is separated from the southern Trecchina-Laino ridge by a tectonic line with a transcurrent trend oriented N110°-120°that has cut the existing plicative structures. Cretaceous limestones of the Serra San Filippo and Serra La Nocara reliefs and Jurassic-Cretaceous limestones of the Rossino and Fossino Mounts constitute the southern morphostructure.
The transcurrent deformation, active from the Upper Pliocene, whose faults are subsequently reactivated in the extensional system [22,23], produces a complex structural pattern [24]. Several geological and structural studies of the entire area [21,23,25] made it possible to place a time constraint to the transcurrent activity, which is still present during the Sicilian, while the extensional tectonic, still active, develops starting from the passage between the Lower and Middle Pleistocene.
The tectonic structures oriented approximately N120°a re extremely important, as responsible for the genesis and evolution of many quaternary basins, including those present in the Calabrian-Lucanian border [22].
The morphostructure of the northern sector of the Lauria Mountains, oriented in the WNW-SSE direction, is located to the north-west sector of the Mercure Basin, which separates the carbonate structure from the neighbouring mountains of the Pollino Chain [26].
The study area highlights different fault systems, with predominantly Apennine faults orientation N120 ± 10°and distributed fault systems oriented N-S and N30°-50°, generally in association with thrusts and folds [26] (Figure 2).
The carbonate morphostructure is bordered by highangle transcurrent and extensional faults that generated a series of depressions filled by Quaternary continental deposits and by Miocene terrigenous deposits [25]. Carbonate successions are bounded to the north by an important thrust front with N-NE vergence, produced by the translational tectonics that led Liguridi and Sicilidi-Affinity Units to overlap on Carbonate Units of the Apennine Platform and on Lagonegresi Units [27].
The presence of the fault with NNW-SSE direction, which plunges to WSW with a fault throw of about 400 m, divides the limestone-dolomitic area in two zones, the first one to the west and the second to the east of the mentioned fault.
This fault starts from Lago La Rotonda and, passing between Castello Starsia and the right side of Serra Tornesiello, reaches Campo del Galdo, in correspondence  In the study area are different distinguishable stratigraphic-tectonic units that constitute the structure of the Lucanian Apennines [19] (Figure 2).
Carbonate Units of the Campanian-Lucanian Platform, represented by the calcareous and calcareous-dolomitic Mesozoic succession [32], overlay the Lagonegro Units. Carbonate Units constitute the entire structure of the Lauria Mountains. Carbonate successions belong to the tectonic units of the Maddalena Mounts and Foraporta Mount and to the Alburno-Cervati-Pollino Unit and comprehend, from the bottom to the top, red shales and cherty limestone of the Middle Triassic, followed by a powerful succession of dolomites (Upper Triassic) and limestones, with subordinate marls [32].
Liguridi Units, present in the area surrounding the northern morphostructure of the Lauria Mountains and extensively to the north-east sector of Castelluccio Superiore, include the turbidite succession of the Liguride Complex (Cretaceous-Oligocene) [33][34][35].
The aforementioned deposits appear to be those geometrically higher than the south-Apennine Chain and mainly outcrop along the western border of the morphostructure, at N and NE of Serra Rotonda.
Quaternary deposits include different clastic deposits such as gravelly-sandy and subordinately sandy-loamed deposits, as well as clayey and marly continental sediments of the Noce and Mercure fluviolacustrine basins. These deposits are present mostly in the SE part of the study area attributable to different accumulation phases and intercalations of paleosols formed by brown silty clays with skeletons of different sizes, which testify to repeated periods of stasis in alluvial sedimentation. Sandy-gravelly deposits, colluvial deposits, and detrital cones of the Upper Pleistocene-Holocene outcropping almost everywhere along the southern slopes of the hydrostructure form a series of alluvial fan deposits that are interdigitated to the middle-Pleistocene lacustrine deposits of the Mercure Basin [22].

Geomorphological and Structural Factors' Role on the
Aquifer System. The geostructural and geomorphological peculiarities of the morphostructure strongly affect the complex hydrogeological setting of the hydrostructure, the groundwater flow directions, and the presence of the springs. The carried-out hydrogeological and hydrogeochemical surveys highlight that the geometrical configuration of the carbonate hydrostructure of the Lauria Mounts northern sector is structured according to different carbonate aquifers. Major faults and karst or fracturing conditions define the hydrogeological watersheds, to places with underground water interchanges, with different hydrogeological and hydrodynamic characteristics, and organize in the carbonate rocks' distinct groundwater flow patterns.
Based on the distinct geolithological, structural features and relative permeability of the stratigraphic-tectonic units, different hydrogeological complexes have been distinguished in the study area [40], each of them characterized by common spatial arrangement, prevalent permeability, and a degree of relative permeability restrained in a narrow range [41].
Limestones, dolomitic limestones, and dolomites constitute the Carbonate Complex that host the most important aquifers drained by the springs with major discharge amounts. It can be divided into three subcomplexes depending on the presence or absence of the karst process and forms (dolines and cavities) and the density of rock fracturing state; the hydraulic conductivity ranges from 10 -4 m/s to 10 -2 m/s [42], and the relative permeability varies from medium-high to very high, representing the favoured complex for the infiltration.
The Conglomerate Complex, mainly outcropping in the south-eastern part of the investigated area, consists of conglomerates and breccias. Characterized by a medium-high degree of relative permeability, it presents a discrete groundwater flow fed by the carbonate complexes. This complex stores part of the groundwater amount and allows the emergence of smaller springs.
The Calcarenite Complex, made of a gray well-layered calcarenite succession with frequent intercalations of marls, has a subordinate role in the local hydrogeological framework, which presents a medium degree of relative permeability due to the presence of the rock fractures.
The Detritical Complex, constituted by heterogeneous conglomerates and gravels with different cementation degrees, includes the detritical and clastic deposits that constitute the alluvial fans and the debris cones. Its hydraulic conductivity ranges from 10 -7 m/s to 10 -5 m/s with a medium degree of relative permeability [42].
The Fluviolacustrine Complex, characterized by a medium-low degree of relative permeability, comprehends clays, marly clays, layered clayey marls with sandy layers, and polygenic conglomerates rich in sandy matrix. It is largely made up by the quaternary deposits of the Mercure Basin that borders the carbonate hydrostructure along its southern boundary.
The Clayey-Marly Complex includes pelagic and chaotic sequences of varicoloured clays, marls, and marly limestone, belonging to the Sicilidi-Affinity and Lagonegresi Units and quartzites; clays and marls referable to the Frido Unit, in which even ophiolitic elements can be distinguished [43][44][45]; and arenaceous-pelitic alternation belonging to the Flysch of Albidona. This complex has a low degree of relative permeability, borders the carbonate hydrostructure, and plays an important hydrogeological role representing the limit of permeability in some parts of the system; the hydraulic conductivity is <10 -9 m/s [42].
The geometry of the carbonate hydrostructure and the underground watershed are mainly defined by the structural elements present throughout the entire hydrostructure, such as faults and thrust, tectonic deformation bands, and karst and fracturing state of rocks.
The assessment of the karst phenomena in the area was based mainly on the recognition of the surface karst forms; the entity of the underground karstic development is evinced considering the geomorphological context and the spring's emergency areas that leads to the consideration of the presence of evolved fractures and karstic conduits. The presence of the karst morphologies, on the dolomitic limestones, suggests a local prevalence of the calcareous component on the dolomite. They represent the surface trace of a developed underground karst network.
This specific hydrogeological setting bordered laterally by major faults and characterized by subsurface boundaries, due to tectonic elements, lithological arrangement, and development of karst network, influences the different groundwater flow directions and the emergence of numerous springs, distinguishing two aquifers with different hydrogeological and hydrodynamic characteristics ( Figure 3).
The Lauria and La Spina-Zaccana aquifers, both characterized by the presence of important springs, located in the western and eastern parts of the carbonate massif, at altitudes between 315 m and 504 m a.s.l., represent strategic 5 Geofluids water resources for the area and for the Basilicata region. Several other emergencies with smaller groundwater amounts are present along the border of the carbonate hydrostructure ( Figure 3).
The groundwater flow, conditioned by different factors such as tectonic arrangement, rock fracturing degree, karst activity, and presence of the basal dolomite complex at the south-eastern part of the aquifer system, essentially occurs in two preferential directions.
The prevailing NE-SW and NW-SE directions are oriented towards the areas where the major springs emerge, most of which are not adequately exploited or even uncollected.

Hydrogeological Features of the Lauria and La
Spina-Zaccana Aquifers. The Lauria Aquifer is characterized by an evident stratigraphic-structural complexity. The carbonate structure is affected by a dense network of high angle direct faults, which develop the geometry of en-echelon from Lago La Rotonda to Campo del Galdo. This direct fault system shows a vertical fault throw of at least 400 meters and gets down the entire western block of the Lauria Ridge with respect to the La Spina-Zaccana Ridge.
The most evident consequence of this remarkable dislocation is represented by the difference of the lithological outcrops between the two blocks: in fact, while the limestones outcropping mainly characterizes the lowered western block of the Lauria side, the eastern block is largely constituted by the outcrops of the deepest dolomitic limestones and dolomites.
In correspondence of this fault system, the distribution of the layered carbonate succession defines an anticlinal fold with the axial zone localized between Castello Starsia and La Spina Mount; it can be considered the tectonic structure that presumably corresponds to the hydrogeological watershed between the two aquifers.
Argillite-marly flysch and clayey-marly flysch bound, to the north and to the west, the carbonate aquifer, consisting limestones and dolomitic limestones. Argillite and marly limestones outcrop to the south, and the most markedly dolomitic terms emerge to the east characterized by an intense fracturing state. The groundwater pathways, influenced by the presence of different faults, has the preferential flow direction oriented NE-SW. It starts from the hydrogeological watershed with the La Spina-Zaccana aquifer and is directed towards the main emergencies located at W of the hydrostructure. The main springs, located in the municipality of Lauria, such as Caffaro-Mandarino, Caffaro, and Arena Bianca, reach a total capacity of about 845 L s -1 [46].
The three different important emergencies of Caffaro-Mandarino with the flow rate of about 0.81 m 3 s -1 [46], located on the edge of the carbonate structure, is tamponed by the clayey-marly flysch formations of the Sicilidi-Affinity Unit. The group of Caffaro-Mandarino springs represents  6 Geofluids the delivery points of the deep aquifer groundwater flow. Nowadays, its discharge is partly used for hydroelectric purposes. The Caffaro Spring did not naturally emerge, but the derivation of the draining tunnel was realized for a hydroelectrical system; the flow rate of the Caffaro is equal to 25 L s -1 [46]. Some quantities of water flow in the detritical deposits or feed the minor springs: Alagia, Martino, and Montepesco, placing a different altitude, refer to shallow underground water circuits that characterize the variable regime of the small flow rates. Those positioned on average from 645 to 1120 m a.s.l. are located in two distinct areas enclosed between the northern slopes of Castello Starsia and the southern side of Serra Tornesiello. Starsia and Menavoli springs with 4 L s -1 flow rates are used for local needs (Figure 3). The presence of some springs at the northern edge of the hydrostructure is due to the underground water transfers of the carbonate rocks towards the flysch deposits in the most permeable areas and towards the detrital cover.
The Although dolomitic limestone and dolomite are considered part of a single succession in original stratigraphic continuity, a band of intense cataclastic deformation is evident at their contact, which is therefore to be considered of tectonic type. It has already been recognized in the study area [20,47] and interpreted as a low-angle tectonic contact with a predominantly extensional or transtensive kinematics of a younger on older [48]. The carbonate aquifer is bounded almost in all directions by exclusively tectonic contacts of different nature and kinematics and stratigraphic places.
The north-western boundary of the hydrostructure is located in the anticlinal hinge zone of the Castello Starsia-La Spina Mount. At its north-eastern edge, however, the aquifer is bounded by a thrust, to places lowered by some normal faults, which superimposes the argillite-marly formations on the carbonates. Along the southern edge, the La Spina-Zaccana substructure is bounded by an important fault, indicated as Castelluccio Fault, with high inclination, oriented in the N120°direction and submerged towards SSO. It is formed in the Pleistocene with a left transcurrent kinematics subsequently reactivated as a direct fault in the upper Pleistocene.
The high and steep slope that characterizes the southern side of the La Spina-Zaccana ridge morphologically highlights the high rate of activity of this fault and its considerable vertical throw. Within its limits, the La Spina-Zaccana morphostructure does not seem to be crossed by major faults, except for the NS subvertical direct fault present in the southern slope that intersects the Castelluccio Fault and which lowers the eastern sector.
The Castelluccio Fault acts as the morphostructural southern limit of the structure and from a hydrogeological point of view represents the superimposed permeability threshold of the San Giovanni spring front.
The La Spina-Zaccana carbonate aquifer represents significant groundwater resources, well fed and drained in large part from the San Giovanni spring; minor emergencies of less potentiality are scattered throughout the south-eastern sector of the structure (Figure 3).
Along the entire southern slope of the La Spina-Zaccana Ridge, the Castelluccio Fault places in tectonic contact the carbonate rocks of the aquifer with the Pleistocene fluviolacustrine deposits of the Mercure Basin. The geometry of the aquifer, its main hydrogeological characteristics, the underground water circulation, and the groundwater emergencies are closely connected to the geological, structural, and karst arrangement. The recharge area, essentially located at high altitudes of the hydrostructure, includes the main reliefs of the study area, towards the W-E direction.
The groundwater flow, predominantly, in the NW-SE direction, feeds the San Giovanni springs group, Acqua del Lavatoio and Pantanello springs along the eastern edge of Monte Zaccana. These springs have a total flow rate of about 580 L s -1 [46]. The San Giovanni spring front (at least 7 normal permanent springs) (mean discharge of 450 L s -1 ) [46], one of the most important groundwater resources of the Lucanian Apennines, is located along a steep slope at an altitude of between 480 m and 504 m. Groundwater is mainly, if necessary, exploited for drinking purposes and irrigation use. Other water emergencies, such as Peschiera, Salice Sopra, Salice Sotto, and Sorgituro, drain limited surficial aquifers characterized by shallow groundwater flow paths.

Sampling and Analytical
Methods. Groundwater samples were collected from 15 springs, some of them partially used for drinking, irrigation, and hydroelectrical purposes, during a single field trip in May 2017. During the rainy season, another sampling was carried out only for isotopic analyses (November 2017).
The sampling sites, located mainly in the western and eastern sectors of the hydrostructure, were selected based on the complex hydrogeological geometry and hydrodynamic characteristics of the aquifers (Figure 3).
Groundwater springs were sampled to define the main hydrogeological and hydrogeochemical processes controlling the chemical properties. Physicochemical parameters including groundwater temperature, pH, and electrical conductivity were measured in situ. One filtered aliquot (using 0.45 μm cellulose acetate filters) and one filtered and acidified (with HNO 3 ) aliquot were collected and stored in lowdensity polyethylene bottles for laboratory analysis of major elements. The major elements were determined by ion chromatography, IC Metrohm 850 Professional. Cations (Na + , K + , Ca 2+ , and Mg 2+ ) were separated by a Metrosep C 4 250/4.0 column using 3.0 mM HNO 3 as eluent and a flow ) by a Metrosep A supp7 250/40 column using 3.6 mM Na 2 CO 3 as eluent at a flow rate of 0.7 mL/min. Alkalinity was determined by titration with 0.1 N HCl. The analytical precision was always better than 5% for each analysed species, and the ion balance (i.e., the difference between the cations and anions sum) was always within 3%.
The geochemical software PHREEQC 3.0 [49,50] was used to calculate aqueous speciation of the investigated elements. Saturation index, defined as SI = log IAP/Kt , IAP being the ion activity product of the mineral-water reaction and Kt the thermodynamic equilibrium constant at the measured temperature, were calculated by ion activities.
The field chemical composition of the analytical results and the saturation index (SI) values for calcite and dolomite are reported are listed in Tables 1 and 2. Validation procedures characterize the analytical performance of laboratory tests in order to understand its capability and limitations. The limit of detection and the limit of quantitation (LOD and LOQ) are the parameters used to describe the analytical sensitivity of an analytical procedure.
The detection limit, LOD, represents the smallest concentration of an analyte in a sample that can be detected with reliable statistical certainty under the stated conditions of the where X b1 is the mean concentration of the blank and S b1 is the standard deviation of the blank. The LOQ is the lowest analyte concentration, determined with acceptable precision and accuracy under the stated conditions of the test, and corresponding to the sample blank value plus ten standard deviations as shown in the following equation: where X b1 is the mean concentration and S b1 is the standard deviation of the blank. The weakness is that there is no objective evidence to prove that a low concentration of analyte will indeed produce a signal distinguishable from a blank (zero concentration) sample [51]. Results of the detection and quantification limits of the analyzed elements are expressed in Table 3.
To evaluate the interrelationships affecting the distribution of the major chemical elements, an R-mode factor analysis was performed.
Factor analysis provides a way to reduce the data to an easily interpretable form, by identifying chemical signatures expressed through statistically significant factors determined by the analysis.
R-mode factor analysis was undertaken by the Statgraphics Centurion XVI package [52], and factor component matrix is extracted (Table 4).
This matrix, showing the impact coefficient of the two extracted components, evidences the link between these components and the groups of variables.
The 18 O/ 16 O and 2 H/ 1 H ratios were determined by isotope ratio mass spectrometry (IRMS). Oxygen isotope measurements were carried out by using a Gas Bench peripheral coupled with a Thermo Delta V mass spectrometer. A TC-EA peripheral interfaced by means of a ConFlo IV with a Thermo Delta XP mass spectrometer was used for hydrogen isotopes. Results are expressed in terms of d-values (i.e., δD and δ 18 O) with respect to the Vienna Standard Mean Ocean Water (VSMOW) as follows: where R sample is the stable-isotope ratio ( 2 H/ 1 H or 18 O/ 16 O) of the sample and R VSMOV is the ratio for VSMOW. The analytical error for IRMS was ±0.1‰ for δ 18 O and better than ±1‰ for δD. The isotopic data of the two sampling periods are reported in Table 5.

Results and Discussion
5.1. Chemical Composition of Groundwater. Groundwater samples are characterized by low electrical conductivity    The hydrogeochemical facies of groundwater, determined using the Piper diagram [54], shows that the investigated groundwater has a homogenous distribution with bicarbonate alkaline-earth composition (Figure 4).
To assess the role of the dissolution mechanisms of carbonate compounds the relationship of HCO 3 to Ca 2+ Mg 2+ has been verified. The observation points tend to fall above y = x (dissolution of carbonates), indicating that Ca 2+ and Mg 2+ ions would originate from calcite and dolomite mineral dissolution in the carbonate rock areas ( Figure 5).
In consideration of the geological context in which the investigated springs occur, it is to be considered that this correlation is mainly due to the dissolution of the calcium carbonate and in a minor way of the magnesium carbonate, dominant in these carbonate deposits.
The carbon dioxide of the soil reacts with water and carbonate rocks dissolving calcite and dolomite according to the following reactions: In Figures 6(a) and 6(b), the saturation indices (SI) for calcite and dolomite, expressed in logarithmic form, in relation to Mg/Ca ratio are shown. Most groundwater samples are saturated or nearly saturated with respect to calcite (Figure 6(a)). On the contrary, almost all of groundwater samples are undersaturated with respect to dolomite (Figure 6(b)).
The dolomite dissolution and concurrent precipitation of calcite maintain water-mineral equilibrium (dedolomitazion) and increase Mg/Ca ratios along flow paths. The scatter plot of Na + versus Cl - (Figure 7) indicates that groundwater are well correlated (r 2 = 0 75) and show a Na/Cl ratio higher than seawater (Na/Cl = 0 86) and halite dissolution (Na/Cl = 1). The Na + /Cl − > 1 relationship may be  , to exclude the contribution of ions from carbonates and silicates. Na + concentration is corrected by Cl − to exclude Na + from atmospheric deposition [55].
Our data do not show a good correlation suggesting that the ionic exchange process is not possible. Therefore, the relative Na-excess is probably due to interaction of water with alkaline feldspar belonging to clastic rocks of the Frido Unit.
Factor analysis allowed verifying the important interdependences among the chemical elements. The factor component matrix (Table 4) determined by factor analysis explains the relationship between the components and the variables.
The analysis shows that component 1 is related to the sodium, sulphate, and chloride variables while electrical conductivity, calcium, potassium, and bicarbonates are linked to component 2.
The physical interpretation of the two factors is graphically evidenced in the two-component plot of Figure 9, where two main groups of variables related to physicochemical data of groundwater samples are shown and, even if they both fall within the positive value field, are anticorrelated.
In component 1, it is clear that the ions of Na + , SO 4 2-, and Cl -, even if characterized by low concentrations, have a high contribution in terms of correlation, indicating the interaction of water with silicate minerals of the arenaceous-pelitic flysch deposits. The loading of Na and Cl and their interrelation, can derive also, to a lesser extent, from the contribution of the sea spray provided by the high rainfall amounts of the study area, due to the orographic conditions and the nearness of the hydrostructure to the sea. The high permeability of the carbonate formations due to the intense fracturing state and karst processes promotes the infiltration of the precipitations into the groundwater system.
The results of component 2 indicate that the high loading factors of EC, Ca 2+ , HCO 3 -, and K + can be associated with the dissolution of carbonate into the aquifer system, which increases the concentrations of these ions. Carbonate dissolution is the most important process affecting the ion concentrations of the groundwater of the samples. The loading factor of Mg 2+ represents the processes related to the dissolution of the dolomites and dolomitic limestones.  )] vs. (Na + -Cl − ). The investigated water samples are not correlated highlighting that the ion exchange process does not take place. 12 Geofluids Moreover, in the carbonate environments, in the presence of the arenaceous and marly flysch deposits, the effect of water-flysch interactions acts as the dissolution mechanisms of the silicates involving the release of detectable quantity of Ca 2+ , HCO 3 -, Mg 2+ , K + , Na + , and H 4 SiO 4 [56,57]. The factor scores of K + and NO 3 -, characterized by very low concentrations, do not have any discernible pattern, and the nitrate ion has no significant lithologic source; the correlation degree suggest that it may be associated mainly with zootechnical activity or surface runoff of fertilizer used for agricultural purposes. 18 O Groundwater Data. Groundwater shows isotopic compositions ranging between -7.5‰ and -8.6‰ for oxygen and -42‰ and -54‰ for deuterium.

δD and δ
The δD and δ 18 O relationship is plotted in Figure 10. All groundwaters plot between the Global Meteoric Water Line (GMWL) [58] and the Mediterranean Meteoric Water Line (MMWL) [59]. This indicated that groundwaters have a meteoric origin and secondary processes occurring after the precipitation (i.e., evaporation, evapotranspiration, and/or water-rock interaction) do not significantly change the isotope ratios of groundwater.
Most of the groundwater springs sampled in November 2017 show an enrichment in O and D heavier isotopes than do data measured in May 2017. The difference in δ 18 O isotopic composition between the samplings is superior to the analytical error (0.1‰ and 1‰ for δ 18 O and δD values, respectively). This is possible because in these springs, the contribution by seasonal rainfall is present, and this may be due to an aquifer with low water capacity and relatively short and shallow hydrogeological circuits through the sedimentary units. Only a limited number of groundwater (P10, P11, and P12 springs), sampled in November 2017, show isotopic values similar to data recorded in May 2017. These springs are not influenced by the seasonal variation of the meteoric recharge. Therefore, under these conditions we may assume that these springs drain the deep aquifers characterized by deep groundwater flow pathways, where the temporal variations of the isotopic ratios are usually within the measurement error. In addition, these groundwaters, showing mean oxygen isotope ratios between -8.4‰ and -8.6‰, are similar to the springs, located at high altitudes >900 meters above sea level (a.s.l.), and pertaining to relatively short and shallow hydrogeological circuits. This finding allows estimating that the average altitude of the recharge areas for the major springs characterized by homogenous groundwater isotopic data ranges from 900 m to 1000 m a.s.l.

Conclusions
This study focused on the hydrogeological setting and groundwater hydrogeochemical assessment of the carbonate hydrostructure of the Lauria Mountains northern sector (southern Apennines). The complex hydrogeological environment and the hydrodynamic characteristics of the aquifer system reflect the geostructural peculiarities of the carbonate hydrostructure. According to the stratigraphic and tectonic setting, deep groundwater flow pathways feed the major springs. The presence of these important springs plays a considerable socioeconomic role for the interest of the entire territory in terms of water availability.
Results highlight that hydrogeochemical and stable isotopic investigations are useful to defining the hydrogeological conceptual model of the aquifer system and representing the essential tool to characterize groundwater chemical assessment. Hydrogeochemical analysis elucidates that the chemical compositions of groundwater depend from the lithology and from the hydrodynamic characteristics of the systems. It suggests that carbonate dissolution is the controlling factor for groundwater chemical properties. The factor 13 Geofluids analysis has demonstrated a greater usefulness in interpreting the hydrogeochemical data relating to water-rock interaction processes with carbonate rocks and silicate minerals.
The main aim of this study was to achieve a deeper understanding of the hydrogeological system and groundwater hydrogeochemical characteristics of the carbonate hydrostructure, using a hydrogeological and hydrogeochemical integrated approach, even in hydrogeological systems with limited data availability, finalized at identifying available groundwater resources. At present, it is possible to argue that considerable groundwater amounts can be used for supplementary and substitutive uses, or to compensate for water emergencies, through sustainable and opportune management actions and strategies.
The current changes in hydrogeology should be faced by improved focus on monitoring techniques and exploitation of information; therefore, the relevance of this carbonate hydrostructure requires further and more detailed knowledge. The findings indicate that effort should be directed, in the future, to the monitoring of the springs discharge, to the detailed isotopic study of rainfall and groundwater, for a better understanding of the preferential groundwater flow directions, of recharge time, taking into account the impact of climate change, in the perspective of protecting these resources from depletion and contamination.

Data Availability
The data used to support the findings of this study are included within the article. 14 Geofluids