Assessing Climate Change and Land ‐ Use Impacts on Drinking Water Resources in Karstic Catchments (Southern Croatia)

: The Mediterranean freshwater resources, mostly represented by groundwater, are under increasing pressure due to natural and anthropogenic factors. In this study, we investigated possible negative effects of climate change and land ‐ use practices on water quality and availability from five springs in the karstic catchments in southern Croatia. The investigated springs are used in the regional public water supply system. Firstly, we employed hydrogeochemical field and laboratory analyses to detect possible traces of anthropogenic activity originating from specific land use. Additionally, we performed hydrological and climate modeling to detect changes in the air temperature, precipitation, and runoff. In particular, we used three regional climate models (Aladin, RegCM3, and Promes). The results estimated an increase in the mean annual air temperature, changes in the precipitation patterns, and reductions in runoff in the study area. Hydrochemical analyses showed standard ion concentrations for karst groundwaters, elevated sulfates due to evaporite deposits in the hinterland, surprisingly low nitrate levels which disproved expected agricultural pollution, and high microbiological activity. Significant water losses are expected in the near future which require immediate attention in order to develop adaptation strategies that focus on sustainable utilization and resilience of freshwater resources. This paper was based on the Interreg Central Europe PROLINE CE project research in the South Dalmatia.


Introduction
Sustainable management of fresh water resources is one of the key environmental challenges of the 21st century. Degradation of its quality or unsustainable management of its quantity can cause environmental and social issues, leading to great consequences. Water quality is affected by various natural and/or anthropogenic factors, and their relative influence changes depending on the scale of observation and the investigated time span [1]. Among anthropogenic factors, land-use practices and changes are often the main drivers of the degradation of water quality [2]. Land use describes the function and processes that take place on the land, while land-use change describes the processes by which natural land is transformed by human activity [3]. The last two centuries have been characterized by rapid industrialization and global population growth that intensified negative land-use change and unsustainable practices, such as: (i) the expansion of artificial surfaces and agricultural areas on the account of forest and wetland areas, (ii) poor waste management, and (iii) unmonitored use of plant protection products. Unsustainable and improper land practices could be reflected in: (i) the degradation of phosphates, chemical residues, and insoluble mineral particles which become infiltrated into the groundwater, surface water, and soil ( Figure 1). A global trend of fertilizer consumption increase over the years has changed the nutrient composition of runoff and leachate resulting in groundwater, surface water, and soil quality deterioration [21][22][23]. The diagram in Figure 1 shows the various stresses that impact water resources. In the coming years, adverse climate change influence will only continue to grow and more severe floods, heat waves, and consequential water losses are to be expected. Land misuse through extensive clear cut of the forested areas which obstruct evapotranspiration, urban impervious surfaces that hinder water infiltration, pesticide and fertilizer overuse, and industrial leachate discharge have a detrimental effect on the quality of water. Additionally, the lack of sanitary protection zone implementation on any jurisdictional level poses a grave issue for the proper resource preservation. Therefore, it is essential to adjust to the changing climate by lessening the harmful anthropogenic effect and adapting a more sustainable approach.
This study aims to provide a hydrogeological and hydrological characterization of the regional karstic groundwater resources in southern Croatia, and to assess the impacts of climate change and land-use practices on its quality and quantity. To detect the pollution originating from various types of land use (e.g., agricultural, urban), we conducted hydrogeochemical analyses and in situ measurements on groundwater samples from five springs in southern Croatia. Therefore, it is necessary to demonstrate how the hydrochemical parameters could be used as indicators of catchment processes [24]. Furthermore, we employed hydrological modeling to determine the impacts of imminent climate change on water resource quantities. The investigations were carried out within the scope of the PROLINE-CE project, funded by the transnational cooperation program Interreg Central Europe.

Geographical Setting
The study area is located in the south of Croatia, in the region of Dalmatia (Figure 2). It is a relatively narrow belt on the eastern shore of the Adriatic Sea, characterized by a well-indented coastline and 79 islands parallel or subparallel to the coast. The study area was divided into two parts: (i) Imotski polje in the northern part, and (ii) Prud, Klokun, and Mandina mlinica spring catchments in the southern part ( Figure 2). Imotski polje was singled out as an individual area due to the characteristics of karst polje, i.e., a closed limestone depression (plateau), usually with springs on one side of the polje and ponors (steep-sided swallow-holes) on the other side [19]. Prud, Klokun, and Mandina mlinica were clustered together due to their proximity and the overlapping catchment areas. The climate in the northern part is characterized by the transition between the Mediterranean and sub-Mediterranean climate, with cold, humid winters, and very dry, warm summers. The climate in the southern part is the Mediterranean climate with mild, humid winters, and hot, dry summers. Mean annual air temperature in the period from 1961 to 1990 was 12.2 °C and 12.7 °C, and mean annual precipitation was 1402 mm and 1424 mm in Imotski polje and catchments of Prud, Klokun, and Mandina mlinica, respectively. The vegetation cover is the typical Mediterranean, with dominant holm oak, Aleppo pine, Dalmatian black pine, maquis shrubs, and garrigue [25,26].

Hydrogeological and Hydrological Settings
The region of the study area belongs to the geotectonic unit of External Dinarides, which are a part of the Dinaric mountain system [27]. Mesozoic carbonate rocks that are characterized by very deep and irregular karstification prevail. As a result of the dissolution of carbonate rock mass by water, specific karst morphological forms, such as karrens, dolines, caves, pits, and ponors are formed. The study area is characterized by an extremely complex geological setting that results in its heterogeneous hydrogeological setting. Permeable deposits are primarily represented by Cretaceous limestones and dolomites in which karst aquifers were formed. Impermeable flysch deposits, located along the coast and to the north of Imotski polje, represent the regional hydrogeological barrier and play a significant role in which they prevent major groundwater outflows to the sea, and divert groundwater flow in a general direction of NW-SE towards the southern part of the study area [18] (Figure 2).
Imotski polje is dominantly made of Cretaceous limestones. The infiltration of meteoric water is quick due to the high porosity/permeability of the formations as a result of the intense rock fracturing and the superimposed karstification. Quaternary alluvial deposits (sand, gravel, loam mixed with terra rossa and humus) up to 100 m thick, represent a low permeability lateral barrier to the groundwater flow ( Figure 2). The largest and most significant spring in Imotski polje is Opačac, which is formed at the contact between the permeable carbonate formations and the impermeable Quaternary cover. It is the main spring feeding the Vrljika River, which is the main recipient in the area. The mean annual discharge measured at Opačac from 1989 to 2001 was 6.62 m 3 /s [28]. The spring is utilized for the water supply of Imotski and the surrounding hinterland (approximately 15,000 inhabitants) and the abstraction rates are up to 0.2 m 3 /s [29]. The major part of the Opačac catchment is located in the neighboring Bosnia and Herzegovina. Significant changes in the hydrological regime of Imotski polje were induced by the construction of the accumulation Ričice. Despite the water regulation infrastructure, a large part of Imotski polje is still susceptible to periodical flooding. Along its northwestern edge, the world-renowned karst hydrogeological phenomena of Crveno and Modro lakes are located. These lakes, hosted in deep pits, testify to the complex processes that shape the karstic environment and represent the outstanding site of natural heritage, which is under evaluation for obtaining the status of UNESCO Global Geopark.
Prud, Klokun, and Mandina mlinica are the most important springs of the southern part of the study area. These springs have a vast recharge area that extends deeply into the neighboring Bosnia and Herzegovina and it is mostly made of karstified carbonate rocks. The area is characterized by numerous karst features such as estavelles, ponors, poljes, and sinking rivers. Unlike in the Imotski polje, the surface hydrographic network is poor and is usually associated with rivers in karst poljes. The only exception is Baćina lakes, a morphological depression where seven lakes occur with a total surface area of 1.38 km 2 . Baćina lakes are important from the aspect of flood mitigation, as they receive excess flood waters from the Jezero polje to the north and discharge them into the sea near Ploče to the south through two hydrotechnical tunnels.
Prud is located in the easternmost boundary of the study area. It is a spring of the river Norin, a right tributary of Neretva, the largest river in the region. The minimum and maximum discharge of Prud ranges from 2.5 m 3 /s to 20 m 3 /s [30,31]. Abstraction rates for public water supply at Prud are up to 0.382 m 3 /s [29], and the water is distributed to the cities of Metković and Opuzen, and through the submarine pipeline to the islands of Korčula, Mljet, Lastovo, and to the Pelješac peninsula. Occasionally, Prud exhibits increased chloride concentrations as a result of the hydraulic contact and mixing with the Neretva River, in which the seawater intrusion is very pronounced due to natural and anthropogenic factors (e.g., hydropower plants). Additionally, the construction of irrigation canals has increased the vulnerability to seawater intrusion even further [31,32].
Klokun is located in the southern part of the study area. Its average discharge in the period from 1999 to 2013 was 1.4 m 3 /s, and the minimum was 0.13 m 3 /s [31]. Abstraction rates for the public water supply of Ploče city and Gradac municipality are up to 0.170 m 3 /s [29]. Klokun is not susceptible to seawater intrusions due to its distance from the sea, but it is prone to turbidity during heavy rain.
Mandina mlinica is a cluster of coastal springs situated along the Žrnovnica bay and several vruljas (submarine karst springs) that are visible only during high discharge. The majority of springs are brackish due to the contact with seawater. The spring with the lowest chlorinity was utilized in the past for the public water supply of Gradac municipality, however, it was discontinued due to salinization issues in the summer months and the introduction of Klokun into the water supply system. In the summer of 2015, the spring was put to use to fill the water carrier ships that were used to supply the Hvar Island due to the water shortage. The total discharge from the Mandina mlinica springs ranges from 0.5 to 4 m 3 /s, and they are prone to increased turbidity during heavy rain. Tracer test has shown a direct connection of the discharged water in the Žrnovnica bay with waters from the Jezero polje [33]. In this study, we investigated two springs from the site, i.e., Mandina mlinica spring (MMS) and Mandina mlinica abstraction site (MMA).

Land Use and Practices
The study area is dominantly covered by forest and seminatural areas, i.e., broadleaved forest, transitional woodland-shrub, sclerophyllous vegetation, and natural grassland ( Figure 3). Agricultural areas, with arable land, permanent crops, and pastures are mostly concentrated in karst poljes, characterized by fertile soils, gentle relief, and water availability. Despite the high irrigation potential, Croatia is amongst the lowestirrigated countries in the EU, as less than 1% of agricultural land and less than 10% of irrigable land is irrigated [34]. In the study area, irrigation is usually performed by rudimentary canals from the rivers in karst poljes. Artificial areas (i.e., urban areas) do not cover a significant portion of the land. The majority of the population is located in the cities of Imotski, Vrgorac, Ploče, and Metković. A smaller percentage lives in the scattered rural or coastal environment. Regarding industry and transport, the most significant infrastructure is the state highway that passes through the entire study area, the port of Ploče, and several smaller quarries. Wetland areas are represented by Prološko blato in Imotski polje, and inland marshes in the southern part of the study area. These areas are important from the aspect of flood mitigation, as they can retain floodwaters. Regarding the legal protection of water resources, Opačac and Prud springs are protected by the drinking water protection zones (DWPZ in further text). In Croatia, DWPZs for fractures and karstified aquifers are divided into four zones: (i) zone of strict protection regime and surveillance, (ii) zone of strict limitation and surveillance, (iii) zone of limitation and surveillance, and (iv) zone of limitation. Each zone has its characteristic limitations on permitted activities and practices, which are defined by the Croatian Republic's Ordinance on the conditions for the establishment of sanitary protection zones.
In Imotski polje, the dominant agricultural crop is grape for wine production, followed by various fruits and vegetables. Wine cultivars do not require extensive irrigation and are very resistant to floods. In Croatia, vineyards constitute only 3.6% of the arable land, and approximately 32% of the total pesticide use is related to their production [35]. In particular, inorganic copper or sulfur-based fungicides are most commonly used. Furthermore, grape production is considered very intensive from the aspect of fertilizer consumption. Romić et al. [35] reported that the annual consumption of fertilizers in Croatia related to grape production was the second-highest, just after corn production, which is the dominant crop in Croatia. The most commonly used fertilizers are nitrogen-based, followed by potassium and phosphorous-based fertilizers.
In the southern part of the study area, in the catchments of Prud, Klokun, and Mandina mlinica, fruit production is dominant, followed by olive groves and grapes, which also cover a significant portion of agricultural land. The hot spot of the fruit production is Jezero polje to the north of Baćina lakes and Klokun spring, which is characterized by intensive production of strawberries, grapes, and citrus fruits. Animal husbandry constitutes a significant proportion of agricultural activity, and the dominant livestock is represented by cattle, cows, sheep, and goats. Due to the scattered production, grazing is not intensive.
In the recent decades, significant land-use changes did not occur in the study area, except for the increase in natural forest succession due to depopulation and abandonment of previously arable land. A minor change was evidenced in the reduction of wetland areas in Imotski polje due to the expansion of agricultural areas. Input data source is the European Environment Agency & Copernicus LAND Service Corine Land Cover [36].

Identification of Risks for Water Resources
From the aspect of sustainable use of water resources and their protection, several risks and conflicts have been identified. Agricultural practices in the study area can be considered conventional and unsustainable, with little or no consideration for over-use of plant protection products (e.g., pesticides, herbicides) or fertilizers, over-tilling, and planting of monocultures. Furthermore, despite the structural flood defense measures, karst poljes are still exposed to periodical flooding, causing deterioration of both surface water and groundwater quality due to transport of hazardous chemicals, enhanced erosion, and turbidity, as well as social and economic consequences. From the urban aspect, the high pressure is exacerbated by the insufficient extent of the sewage network. In Imotski polje, only Imotski and just south of the city are connected to the sewage system, while most of the population (approximately 80%) relies on septic pits resulting in the potential transport of contaminants and pathogens through a highly fractured karst network. Furthermore, the existing system is outdated, resulting in leaks and possible contamination of soil or groundwater. The system collects both domestic wastewater and meteoric water, and outflows occur during heavy rains. Collected wastewater is purified in a treatment facility before being discharged to surface water recipients.
Similar issues are present in the southern part of the study area, and the only wastewater treatment facility was constructed in Vrgorac. Coastal settlements and cities (e.g., Ploče), discharge wastewater into the sea, with minor or without treatment. The most significant waste disposal sites are Kozjačić in Imotski, Ajdanovac in Vrgorac, and Lovornik in Ploče, and these sites lack the proper barrier liner systems, gas venting, and leachate collection systems. Furthermore, there are numerous speleological objects (caves, pits, swallow holes) in the study area where the waste is illegally dumped. These sites have a high potential for pollution of soil, groundwater, and surface water, as waste is disposed of in barren karst terrain [37].
Despite the current abundance of water resources, climate change projections for Croatia and the whole Mediterranean region indicate an increase in air temperature, changes in precipitation regimes, increase in evapotranspiration, and intensification of floods and droughts in the near future (e.g., [5,9,12,13,[38][39][40][41]), which could have a negative effect on the water availability in the study area. This problem is aggravated by the extreme losses in the water supply system, which are high above the national average of 42%, i.e., up to 80% in the water supply of Imotski, and up to 70% in Prud and Klokun water supply systems [29].

Hydrogeochemical Investigations
As carbonate aquifers belong to the most abundant groundwater sources, their solubility contributes significantly to the higher mineral content of water. In addition to the high chemical weathering rates, the heterogeneous nature of karst and its dependability on hydrological conditions could be reflected in its hydrochemical footprint. To understand the climate change impact on water resources, it is important to study the natural geochemical conditions which occur in the aquifer [42]. The main cations in groundwater originate dominantly from the chemical weathering of rocks, whereas the anions may derive from non-lithologic sources such as agricultural runoff, wastewater, and sewage seepage [43]. Piper [44] used a trilinear diagram for water classification with two separate ternary plots representing cations and anions. The plots are projected into a diamond field used to represent the water's composition with respect to both ionic groups, illustrating the dominant hydrochemical facies. Furthermore, the major ion content of spring water may indicate water residence time since groundwater is derived from infiltrated precipitation and thus, poor in dissolved mineral content. However, the water-rock interaction during a period of time enriches the groundwater in major ion concentrations [45].
Periodical monthly hydrogeochemical field and laboratory investigations of the springs analyzed in this study were conducted from August 2019 to September 2020. Field investigation included in situ monitoring of pH, temperature, and electrical conductivity (EC) by WTW multi-parameter probe. Furthermore, alkalinity, i.e., the concentration of bicarbonate ions, was determined by volumetric titration with 1.6N H2SO4 to pH 4.5. Spring water samples were collected in 200 mL polyethylene bottles that prevent evaporation and were kept at the temperature closest to the temperature of the spring during the sampling procedure. Laboratory investigations were performed in the hydrochemical laboratory of the Department of Hydrogeology and Engineering Geology at the Croatian Geological Survey. Principal ion composition (HCO3 − , Cl − , SO4 2− , Ca 2+ ; Mg 2+ ; Na + ; K + ; NO3 − ) was analyzed by ion chromatography on DIONEX ICS-6000 DP (Thermo Fischer Scientific Inc., Waltham, Massachusetts, USA) from spring water samples of Opačac, Klokun, Prud, and Mandina mlinica spring and abstraction site. This method is best suited for non-metal ion analysis in aqueous solutions [46]. Moreover, three major springs from the study area (Opačac, Klokun, and Prud) are used for the local water supply and are under the surveillance of the water utilities. The authors of this paper obtained microbiological data which included total organic carbon (TOC), aerobic, total and fecal coliform, streptococci, Escherichia coli, and inorganic nitrate content. The microbiological analysis was conducted in the laboratory of the Croatian Institute for Public Health.

Hydrological and Climate Modeling
Hydrological modeling in the study area was conducted based on the size and the water balance elements of the analyzed catchments. We employed the empirical expressions by Turc [47] and Langbein [48], which are most commonly used in the Croatian karst [49] to estimate the runoff based on the spatial distribution of the mean annual precipitation and air temperature.
Turc's formula expresses runoff deficit (D) as a function of precipitation (P) and a temperature factor (L): where D is the runoff deficit as a function of precipitation P and a temperature factor L, given as: where T represents the increasing air temperature. The Langbein formula (1962) is based on a unique relationship between precipitation-temperature factor ratio (P/K) and runoff-temperature factor ratio (Q/K), i.e., it uses precipitation and air temperature as the key attributes when estimating surface runoff. Air temperature (T) is incorporated through a temperature factor (K), which increases exponentially as air temperature increases. Langbein's expression states: The spatial distribution of the analyzed parameters was facilitated by Geographical information systems, in a similar way as it was done in the work of Horvat and Rubinić [50]. The empirical models were calibrated by the hydrological observations in the period from 1961 to 2011. As a referent station, the climatological station Opuzen was used due to having a long time series of historical measurements. The spatial distribution of the mean annual air temperature and precipitation from 1961 to 1990, in a form of a 1000 m resolution raster, was prepared by the Croatian Meteorological and Hydrological Service (DHMZ). The inflow in the analyzed catchments, i.e., the precipitation that infiltrated the karst catchment, was determined by the same empirical expressions and was based on the spatial distribution of the mean annual precipitation and air temperature, as well as the size of the catchment. The validation of the empirical model was based on the comparison of the modeling results and the measured data. Based on the historical meteorological and climatological observations, the synthetic series of data for the period from 2012 to 2050 were generated. By combining the data on the mean annual air temperature and annual precipitation, generated by regional climate models (RCMs), with estimated catchment runoff and effective precipitation, synthetic series of runoff in the analyzed catchments were generated.
Climate change, i.e., the mean annual air temperature, precipitation, and mean annual runoff, up to 2050, have been projected using regional climate models (RCMs), namely, Aladin [51], RegCM3 [52], and Promes [53]. The RCMs were forced by the observed concentration of greenhouse gases from 1951 to 2000, and from 2001 onwards, the IPCC A1B scenario was applied [54]. The initial and boundary conditions for each RCM were obtained from different global climate models (GCMs): (i) Aladin was forced by the Arpege GCM, (ii) RegCM3 by the ECHAM5 GCM, and (iii) Promes by the HadCM3Q GCM.
Furthermore, we have used the E-OBS gridded datasets [55] to determine the change of climate indicators that could have a significant impact on the availability of freshwater resources. In particular, we observed the changes in consecutive dry days (CDD), consecutive wet days (CWD), and the maximum annual 1-day precipitation (Rx1D) in the period from 2021 to 2050 and from 2071 to 2100 versus the period from 1971 to 2000. In the scope of the project PROLINE-CE, the climate predictions of these indicators were obtained and downscaled to the study area by an ensemble of high-resolution EURO-CORDEX regional climate models.

Hydrochemical Analysis of the Study Area
The minimum, average, and maximum values of ion concentration, pH, and electrical conductivity (EC) are shown in Table 1. The average values of the major ions were used to construct a Piper diagram ( Figure  4). At Opačac, Prud, and Klokun springs, the distribution of major ions is characteristic for limestone-dominant aquifers, with Ca 2+ as the most abundant cation, followed by Mg 2+ , Na + , and K + . Mandina mlinica spring site is abundant in Na + and Cldue to sea proximity; Mg 2+ and K + are elevated for the same reason, whereas Mandina mlinica abstraction site shows typical chemistry for karst groundwaters. The observed values of pH ranged from 7.28 to 8.14, and they were within the standard range of the freshwater (i.e., from 6.5 to 8.5). Values of electrical conductivity (EC) at Prud, Klokun, and MMA were elevated but well within the drinking water threshold (<2500 μS/cm) [56], as a result of marine influence due to their position close to the coast. At MMS, even the minimum values of EC and Cl − evidenced total mixing of spring water with the seawater. Contrary values observed at the two adjacent springs at MMS and MMA demonstrated the heterogeneity of karst systems in terms of preferential flows and differences in freshwater/seawater interface. Furthermore, the dissolution of sulfate-rich minerals (i.e., gypsum, anhydrite) from evaporite deposits could also be reflected in the elevated concentration of sulfates (SO4 2− ) and subsequently, EC. Evaporite deposits are common in the central and southern Adriatic basins, and Bosnia and Herzegovina [57]. As the majority of Prud's catchment is located in the neighboring country, evaporite deposits in the underground were the most likely cause of increased sulfate concentration. Additionally, the anthropogenic activity (industry, agriculture, domestic wastewater) could also add up to the high concentration of SO4 2− . Sulfates are abundant in seawater at 2.712 g/L in 35 Practical salinity units (PSU) [58], therefore, a marine origin is also plausible. Perry et al. [59] have given a strong argument for using the SO4 2− /Cl − ration as a seawater mixing indicator. The average value for seawater is 10.3. Adhering to that, the results confirm the authorsʹ assumption regarding seawater mixing in Prud and Klokun (Table 1).
Relatively high bicarbonate (HCO3 − ) concentrations confirmed that the groundwater flow occurs in a carbonate aquifer. Nitrate concentrations were generally low, i.e., up to 5.35 mg/L at all analyzed springs, and they did not exceed the standard drinking water threshold (50 mg/L) [56]. Although there are many possible sources of nitrate in the groundwater, the majority of them are linked to human activities. However, intensive agriculture and husbandry in the study area were not reflected in high nitrate concentrations in the groundwater. The most probable causes of a low concentration of NO3 − were: (i) high discharges at Opačac, Prud, and Klokun (i.e., on average, 4.6, 5.7, and 1.4 m 3 /s) resulting in a high dilution of nitrates; (ii) preferential and piston groundwater flow paths through a highly fractured and karstified underground, and (iii) complex groundwater recharge and mixing patterns in the overlying catchments. The reliable interpretation is constrained by the lack of data on nitrate concentrations from surface waters and soil. In particular, rivers in karst poljes are the main recipients and receive a surplus of nitrogen that is not taken up by plants, and a high concentration of nitrates is expected, particularly during several spraying seasons throughout the year. This was supported by low values of potassium at all analyzed locations except in MMS, where high values of K − originated from the seawater. A similar study was conducted in China [60] where the authors confirmed karst hydrogeochemical parameters such as pH, conductivity, HCO3 − , Ca 2+ , SO4 2− , and Mg 2+ as good indicators of various land use and other environmental changes.  Figure 4). Opačac, Klokun, Prud, and MMA belong to the calcium type, i.e., Ca-HCO3 and mixed Ca-HCO3 and Ca-Mg-Cl facies. MMS belongs to the Na-Cl type due to its coastal location. Since the study area consists of carbonate deposits, i.e., limestone and dolomite, they add a significant amount of Ca 2+ and Mg 2+ to the water.
The bacteria, total organic carbon, and inorganic nitrogen content are shown in Table 2. Assessing the given microbiological data, what immediately stands to attention is the high aerobic bacteria content (6000 units/mL) at Prud. Total coliforms follow suit with 1800 units in 100 mL, as well as fecal bacteria. This reflects the poor purification capacity of karst aquifers due to the thin or non-existent overlying soil and, generally, the absence of fine-grained materials in the aquifer (i.e., clay, sand) [61]. Total organic carbon values and inorganic nitrogen content show a slight correspondence with bacterial colonies, especially for Opačac and Prud springs. Pronounced fractions of fecal bacteria are derived from animal manure coupled with higher values of inorganic nitrogen building a connection with agricultural runoff (fertilizers) [62]. Groundwater pollution generated by intense nitrogen fertilization is a common issue in numerous agricultural communities [62][63][64].
In general, the water is of good quality, with slightly elevated sulfates most likely on account of evaporite deposits from the neighboring Bosnia and Herzegovina. The microbiological data is high, but the water is being chlorinated prior to distribution to end-users and thus properly managed. Even though the study supposed that the continual monthly hydrochemical monitoring over a period of hydrological year is sufficient, the obtained results showed otherwise. The limitations of the hydrochemical study were the sparse data set and lack of continued daily monitoring, therefore, a more intricate groundwater monitoring network is needed for a conclusive study, over a longer period of time [64].

Assessment of Climate Change Impacts
Hydrological modeling and assessment of climate change impacts on water resources were performed for Imotski polje and the catchments of Prud, Klokun, and Mandina, respectively. The measured data on mean annual air temperature and precipitation from 1961 to 2012, as well as the data on the spatial distribution of the runoff, and their generated values for the period until 2050, were used in the analyses. Statistics of the mean annual air temperature and the values generated by Aladin, RegCM3, and Promes are shown in Table 3. Table 3. Statistics of the mean annual air temperature (minimum, average, and maximum) in the study area during 30-year reference periods: (i) from 1961 to 1990 (measured data), (ii) from 1991 to 2020 (partly measured and partly generated data), and (iii) from 2021 to 2050 (generated data). In comparison to the 30-year reference period (from 1961 to 1990), all models evidenced an increase in the average values of mean annual air temperature until 2050. In Imotski polje, Promes evidenced the highest increase in the average value of mean annual air temperature, up to 1.4 °C, followed by Aladin with 0.5 °C, and RegCM3 evidenced a negligible increase of 0.1 °C. Slightly lower values were observed in the catchments of Prud, Klokun, and Mandina mlinica, where increases of 0.9 °C, 0.6 °C, and 0.1 °C were indicated by Promes, Aladin, and RegCM3, respectively. In both parts of the study area, the maximum values of the mean annual air temperature tend to increase more significantly until 2050. For Imotski polje, Promes showed the highest increase, up to 2 °C, followed by Aladin with 1.3 °C, and RegCM3 with 0.5 °C. In the catchments of Prud, Klokun, and Mandina mlinica, slightly lower increases of 1.4 °C, 1.3 °C, and 0.8 °C were evidenced by the models. Regarding the minimum values of the mean annual air temperature, mixed results were obtained, whereas RegCM3 evidenced a decrease of −0.4 °C in both parts of the study area, while Promes and Aladin showed a slight increase.

Mean Annual Air
The linear regressions evidenced increasing trends of 2.3 °C/100 years by Promes, followed by Aladin with 0.9 °C/100 years, for Imotski polje. Similar results were obtained for Prud, Klokun, and Mandina mlinica catchments, where Promes evidenced an increase of 1.4 °C/100 years, and Aladin 1 °C/100 years. In both parts of the study area, RegCM3 provided a negligible increase in the mean annual air temperature until 2050. Due to the long tradition of using RegCM3 in Croatia for climate predictions, the models often display conservative projections, small warm biases in air temperature, and general overestimation of precipitation [9]. Furthermore, the increasing trends in the mean annual air temperature obtained by Promes and Aladin fit well into regional studies and attest to the hypothesis of faster warming trends in the Mediterranean region. In particular, Branković et al. [9] described warming rates in the Adriatic region and showed an increase in the mean annual air temperature of 0.07-0.22 °C/decade for the period from 1951 to 2010, and 0.29-0.71 °C/decade for the period from 1981 to 2010. Furthermore, they projected the most significant increase for the summer and autumn periods, with an increase of up to 5.5 °C until 2100. Similarly, a higher intensity of warming during the warmer parts of the year (i.e., spring and summer) was evidenced by Bonacci et al. [5]. Moreover, the observed results are in accordance with Berkeley Earth [11], grid-based datasets of global air temperature. They described the current mean rate of air temperature change of 2.94 ± 0.26, 3.54 ± 0.33, and 3.35 ± 0.40 °C/100 years for Europe, Croatia, and Split, respectively.
Statistics of the mean annual precipitation and the values generated by Aladin, RegCM3, and Promes are shown in Table 4. Table 4. Statistics of the mean annual precipitation (minimum, average, and maximum) in the study area during 30-year reference periods: (i) from 1961 to 1990 (measured data), (ii) from 1991 to 2020 (partly measured and partly generated data), and (iii) from 2021 to 2050 (generated data). All models pointed out the stable or slightly decreasing trends in precipitation until 2050. In particular, Aladin evidenced a decrease of approximately 57 mm/100 years, RegCM3 of 61 mm/100 years, and Promes by 129 mm/100 years in Imotski polje. Furthermore, a more significant decrease in the maximum values of the mean annual precipitation was evidenced by all models, i.e., by approximately 100-250 mm (Table 4). Promes indicated the most significant decrease, while Aladin and RegCM3 showed similar predictions. Conversely, the minimum values of the mean annual precipitation are expected to increase by 200-250 mm until 2050. A nearly identical situation was observed in the catchments of Prud, Klokun, and Mandina mlinica, where models indicated slightly decreasing trends in precipitation of 57 mm/100 years (Aladin), 62 mm/100 years (RegCM3), and 131 mm/100 years (Promes). Additionally, the behavior of the minimum and maximum values of the mean annual precipitation is similar to the one observed in Imotski polje. These results are in accordance with most climate models for the Mediterranean region which predict a gradual decrease in precipitation and an increase in variability. However, these model-based changes in precipitation are usually characterized by high uncertainty [5,9,12].

Precipitation (mm) Imotski
Based on the available data on the mean annual air temperature and precipitation, and their spatial distribution in the period from 1961 to 1990, the mean annual runoff for the same 30-year reference periods were calculated by Turc and Langbein empirical expressions. Table 5 shows the statistics (minimum, average, maximum) of the time series of the calculated mean annual runoff in the period from 1961 to 1990. Generated time series of the mean annual runoff in the periods from 1991 to 2020, and from 2021 to 2050 according to Aladin, RegCM3, and Promes are shown in Table 6. Furthermore, the time series of the historical and generated values of the mean annual runoff in Imotski polje and catchments of Prud, Klokun, and Mandina mlinica are shown in Figure 5. The results evidenced that the average values of the mean annual runoff in Imotski polje will decrease until 2050. In particular, Aladin and RegCM3 evidenced a decrease of 10% and 8%, respectively, while Promes indicated a 15% decrease when compared to historical runoff (from 1961 to 1990, Tables 4 and 5). The maximum values of the mean annual runoff showed a more significant decrease, namely by 10%, 8%, and 22% according to Aladin, RegCM3, and Promes, respectively. All three models indicated an extreme increase of the minimum mean annual runoff, namely by 44%, 100%, and 88%. In the catchments of Prud, Klokun, and Mandina mlinica, all three models evidenced an almost identical decrease in the average values of the mean annual runoff as in Imotski polje. Projected changes in the maximum values of the mean annual runoff are similar to those projected for Imotski polje, with a foreseen decrease up to 10%, 12%, and 22%. Regarding the minimum values of the mean annual runoff, the same increase as in Imotski polje was evidenced by all three models.
Furthermore, from 1971 to 2000, E-OBS datasets evidenced 34 to 40 CDD in Imotski polje and catchments of Prud, Klokun, and Mandina mlinica. EURO-CORDEX RCMs predicted an increase in CDD by 3 to 8 days from 2021 to 2050, for the 75th percentile under RCP4.5 and RCP8.5 scenarios. In the period from 2071 to 2100, CDD is expected to increase by 8 to 13 days in coastal parts of the study area, and by 3 to 8 days in the hinterlands, where the continental effect of the climate is more pronounced. This increase in CDD could be reflected in the prolonged droughts, especially during the warmer season (i.e., spring and summer) when the hydrological minimum occurs, posing a threat to agricultural production, which is characterized by a low percentage of irrigated areas and a relatively high percentage of production that requires irrigation (e.g., fruits). Regarding the CWD, which is far greater in mountainous regions, E-OBS indicated 9 to 11 CWD in the period from 1971 to 2000 in the study area. Models showed a mixed-signal for the period from 2021 to 2050, i.e., for the 25th percentile, both scenarios predicted a decrease by 1 to 4 days, while stable trends were provided for the 50th and 75th percentile. The same result was projected for the period from 2071 to 2100. Furthermore, the median value (50th percentile) under RCP8.5 predicted reductions in CWD mostly across Italy and Balkan countries, which could have a significant impact on the duration and intensity of the flood waves in the study area. Concerning Rx1D, climate models indicated a clear increasing signal on a short time horizon (2021-2050), and according to the median value, the maximum annual daily rainfall could increase from 5 to 15 mm. The same result was obtained for the 75th percentile. On a long time horizon (2071 to 2100), the most significant increases in Rx1D were detected for the 75th percentile, from 15 to 25 mm in lowland parts of the study area, and 25 to 35 mm in mountainous areas.
Since the pressure on water resources will continue to grow, owing to adverse climate change effects (prolonged drought, consequently, greater seawater intrusion), observation and extended analyses of the endangered water resources is an obligation.

Conclusions
Climate models pointed out an increase in the mean annual air temperature, changes in the precipitation patterns, and the reductions in runoff in the study area. Coupled with changes in CDD, CWD, and Rx1D, a significant negative effect could be induced on the environment and its processes. Despite the current abundance of water resources and the low ratio of utilization of available resources for public water supply, climate models pointed out the possibility of significant losses of water resources in the near future. Additionally, this problem is aggravated by the high uncertainty and lack of knowledge regarding the groundwater reserves in karst aquifers. This shortcoming could be partly solved by establishing a groundwater level monitoring network, which is practically nonexistent in Dalmatia. The available data on groundwater levels would greatly contribute to the determination of the aquifer's response to hydrological extremes (e.g., floods or droughts), over-abstraction for potable or irrigation purposes, as well as its recharge processes. Currently, agricultural demand for fresh water is not exhibiting high pressure on water resources, but the existing spatial plans foresee a construction of a large irrigation system in Imotski polje and Jezero polje, hence the production could shift to cultures that require more extensive irrigation (e.g., fruits). Additionally, conservational use of water resources in the future should focus on increasing the efficiency of existing irrigation (e.g., sprinkler or drip systems), prevention of illegal connections to water supply systems, promotion of cultures that require little or no irrigation (e.g., vineyards), reduction of extreme losses in the water supply systems, and fostering alternative sources of water (e.g., rainwater harvesting).
Since the study area belongs to the rich agricultural zone of South Dalmatia, the authors of this paper expected a definite hydrochemical trace of the intense production. Hydrochemical data presented in this study are insufficient to quantify the land-use impact on water quality. Agriculture is the dominant activity of the region and many bad practices were observed throughout the PROLINE-CE project implementation ranging from excessive use of pesticides, application of restricted pesticides acquired from the neighboring Bosnia and Herzegovina, bad soil remediation after erosion and flood events (e.g., prescribed burning), and improvised irrigation. Poor land-use management greatly amplifies the pressure upon the environment, i.e., the water resources and ecosystems are its most endangered component. Furthermore, the nature of the karstic aquifer should also be kept in mind, i.e., its low retention and low purification capacity, rapid infiltration from surface to groundwater, and generally swift underground flow that makes it vulnerable to pollution with a widespread reach [61]. Additionally, climate change will affect the frequency and intensity of weather extremes, i.e., floods and droughts. Current flood mitigation infrastructure cannot cope with the intensity of high periodical flood waves, especially in the karst poljes. A combination of structural (i.e., new retention area) and holistic approaches should be fostered to provide a balance between the further development of the area and adequate flood mitigation. A holistic approach to flood mitigation includes prevention of construction in flood-prone areas, preservation of natural flood retention land (e.g., pasture, wetland, and forest), improvement of hydrological and meteorological forecasting, and promotion of green, permeable infiltration surfaces in the urban environment.
All of the above mentioned should raise consciousness to sustainable agriculture and improved water resource management as it affects us directly. Results of this study could be used by decision-makers, planners, water suppliers, and other major water users (i.e., agronomists, forest industry, energy production) to improve the sustainable utilization of groundwater resources and to increase the resilience of local communities and natural ecosystems in an uncertain climatic future. The authors of this study suggest continued water quality monitoring in the agricultural areas to anticipate adverse effects linked to anthropogenic and climate change impact. Awareness-raising events should also be a keymeasure in order to encourage more environmentally conscious actions. Decision-makers ought to speed up the implementation of drinking water protection zones that will restrict agricultural activity in the immediate vicinity of a spring used for water supply, as well as foster more green approaches (e.g., organic farming).