Monitoring and Management of Karstic Coastal Groundwater in a Changing Environment (Southern Italy): A Review of a Regional Experience

Similar documents
Monitoring and management of coastal karstic aquifers

The piezometric stress in the coastal aquifers of a karstic region, Apulia, Italy

In press. Fidelibus M. D. and Tulipano L.

Genetic Algorithm based Simulation Optimization Approach to Seawater Intrusion

POSSIBILITIES FOR GROUNDWATER DEVELOPMENT FOR THE CITY OF JAKARTA, INDONESIA

Potential effects evaluation of dewatering an underground mine on surface water and groundwater located in a rural area

TSUNAMI DISASTER AND GROUNDWATER RESOURCES: THE SOUTHEASTERN COAST OF INDIA. Bhanu Neupane Regional Hydrologist, UNESCO

J.Carneiro, A. Correia Geophysical Centre of Évora, Évora, Portugal. A. Rimi Scientific Institute, Rabat, Morocco

SALINE CONTAMINATION PHENOMENA IN THE KARSTIC AND FISSURED CARBONATIC AQUIFER OF THE SALENTINE PENINSULA (SOUTHERN ITALY)At'ffi THEIR EVOLUTION * * *

Numerical Modeling of Groundwater Flow in Karst Aquifer, Makeng Mining Area

SALINISATION OF THE NORTHERN COASTEL AREA OF THE NETHERLANDS DUE TO LAND SUBSIDENCE AND SEA LEVEL RISE

Water Resources on PEI: an overview and brief discussion of challenges

The effect of Hilti / Salahi Recharge Dam in Batineh Area, Sultanate of Oman on the Quality of Groundwater Using GIS

Piezometric Variations of Mikkes Groundwater Basin (MOROCCO)

Estimation of Natural Groundwater Recharge in Qatar Using GIS

Climate Change and Sea Level Rise

Modeling the Managed Aquifer Recharge for Groundwater Salinity Management in the Sokh River Basin

POLLUTION POTENTIAL OF GROUND- WATER RESOURCES IN ANTALYA CITY

AQUA th September 2015, ROME

Analysis of long-term ( ) annual discharges of the karst spring Fontaine de Vaucluse (France)

Groundwater. Groundwater Movement. Groundwater Movement Recharge: the infiltration of water into any subsurface formation.

Geology 627, Hydrogeology Review questions for final exam h t 1/ 2

GLY 155 Introduction to Physical Geology, W. Altermann. Grotzinger Jordan. Understanding Earth. Sixth Edition

CHAPTER 2. Objectives of Groundwater Modelling

B. The use of isotopes in groundwater hydrology

STRAWMAN OUTLINE March 21, 2008 ISWS/ISGS REPORT ON THE OPPORTUNITIES AND CHALLENGES OF MEETING WATER DEMAND IN NORTH-EAST ILLINOIS

Climate Change and the Eastern Mediterranean Precipitation Regime

Indicator Fact Sheet (WQ01c) Water exploitation index

HYDROLOGICAL IMPACTS OF IRAQI BADUSH DAM ON GROUNDWATER

Numerical Simulation and Prediction of Groundwater Flow in a Coastal Aquifer of Southern India

Simulation of horizontal well performance using Visual MODFLOW

DECLARATION OF MALTA

Groundwater modelling to assess the effect of interceptor drainage and lining

1 INTRODUCTION 2 MOTIVATION FOR THE PRESENT EXPERIMENT. Dynamics of the Salt-Freshwater Mixing Zone in Ocean Beaches. N. Cartwright 1 and P.

Water Resources Management: Vulnerability of Coastal Aquifers to Climate Change & Human Effects.

Well Hydraulics. The time required to reach steady state depends on S(torativity) T(ransmissivity) BC(boundary conditions) and Q(pumping rate).

Issue paper: Aquifer Water Balance

Our Eastern Shore Groundwater Part IV Groundwater Quality on the Eastern Shore: How safe is our groundwater and are there ways we can protect it?

WATER AND THE HYDROLOGIC CYCLE

Environmental Indicator: Groundwater in British Columbia

Steady Flow in Confined Aquifer

CHAPTER 13 OUTLINE The Hydrologic Cycle and Groundwater. Hydrologic cycle. Hydrologic cycle cont.

1. Introduction. Keywords Groundwater, Vulbnerability, Aquifer, Aquitard, Vadose zone. Alsharifa Hind Mohammad

Isotope hydrology: Stable isotopes in water resources assessment IAEA International Atomic Energy Agency

MEASUREMENT OF GROUNDWATER TEMPERATURE FOR UNDERGROUND THERMAL UTILIZATION

Mitigation of Seawater Intrusion using Subsurface Intakes of Desalination Plants

Design of a passive hydraulic containment system using FEFLOW modelling

Salt Water Intrusion in Coastal Aquifers March, 2015 L. Pilgrim, P.Geo.

Alternative Futures for the City of La Paz, Mexico

22 Tubewell Drainage Systems

TORO ENERGY LAKE MAITLAND DEWATERING AND WATER BALANCE REVIEW

Keywords: MODFLOW; water availability; water resources management

Aquifer Science Staff, January 2007

Prediction for Natural Recharging In Langat Basin and Ukm Campus as Case Study

Pilot Study. Adaptation in water and coastal areas in Puglia, Italy

Groundwater Bodies in the UK

NOTE ON WELL SITING AT CORNER OF CASTLE VALLEY DRIVE AND SHAFER LANE FOR TOWN OF CASTLE VALLEY, GRAND COUNTY, UTAH

Chapter 2: Aquifers and groundwater

Using groundwater flow simulation of the Chania Plain area to propose a proper irrigation plan

Schlumberger WATER SERVICES. Schlumberger Water Services. Solutions for Water Exploration and Sustainable Supply

Numerical Groundwater Model for the Kaweah Delta Water Conservation District

Numerical Simulation of Variabledensity Groundwater Flow: Submarine Groundwater Discharge

Groundwater Monitoring Protocols for Seawater Intrusion Example of Challenges and Experiences in a Coastal Groundwater Basin

Northern Arizona Hydrogeology

Groundwater Flow Evaluation and Spatial Geochemical Analysis of the Queen City Aquifer, Texas

Monitoring and modelling catchment water quantity and quality

ANALYSIS OF RAINFALLINFILTRATION

Capão Xavier Mine Water Drainage Management (Minas Gerais, Brazil)

Hydrology and Water Management. Dr. Mujahid Khan, UET Peshawar

Groundwater Balance Study in the High Barind, Bangladesh. A.H.M.Selim Reza 1, Quamrul Hasan Mazumder 1 and Mushfique Ahmed 1

Water Framework Directive. Groundwater Monitoring Programme. Site Information. Ballygar

Susan P. Abano Engineer IV Policy and Program Division

Prepared For: Town of Castle Valley, Utah

Capture Zone Analyses For Pump and Treat Systems. Internet Seminar Version: September 4, 2008

ENVIRONMENTAL HYDROGEOLOGY

Comparison between Neuman (1975) and Jacob (1946) application for analysing pumping test data of unconfined aquifer

GROUNDWATER SALINIZATION IN VERSILIA (ITALY)

Uncertainty and Urban Water Recharge for Managing Groundwater Availability: A Case Study and Methods Development for Karst Aquifers ABSTRACT

Post-closure mining pit lake water quality modelling David Dettrick 1, Brian Barnett 2, Damon Grace 3, Hiro Toki 4

Detailed technical design for acid mine water treatment in Novo Brdo Mine, Artana / KOSOVO

CASE STUDY OF THE GROUNDWATER RESOURCES USAGE FOR A COASTAL ARID REGION IN THE AL BATINAH, OMAN

ESTIMATION OF CLIMATE CHANGE IMPACT ON WATER RESOURCES BY USING BILAN WATER BALANCE MODEL

(,,,) = ( )exp ( + C(x,y,z,t) = the concentration of the contaminant at location x, y, z from the source at time t.

Predicting Unmet Irrigation Demands due to Climate Change An integrated Approach in WEAP

SOUTHERN WATER USE CAUTION AREA REVIEW

HYDROGEOLOGY OF THE HUMBOLDT RIVER BASIN, IMPACTS OF OPEN-PIT MINE DEWATERING AND PIT LAKE FORMATION

STRATEGIES FOR CHARACTERIZING SUBSURFACE RELEASES OF GASOLINE CONTAINING MTBE

Importance of irrigation return flow on the groundwater budget of a rural basin in India

UNESCO EOLSS. The various methods available for the treating and disposing of treated sewerage in natural receiving water bodies are outlined.

Subject: Peer review of Duntroon Quarry hydrogeologic modelling: Report #1

UNIT HYDROGRAPH AND EFFECTIVE RAINFALL S INFLUENCE OVER THE STORM RUNOFF HYDROGRAPH

MEMORANDUM. RAI Responses Related to East Lake Road Wellfield Drawdown Analysis, WUP No SDI Project No. PCF-180.

NAMIBIA: Assessing Recharge to the Southeast Kalahari Artesian Basin

What is a representative pan factor value for the eastern Pilbara? Aditya Jha

Downloaded 02/05/14 to Redistribution subject to SEG license or copyright; see Terms of Use at

CHAPTER 4 DISCUSSION. Total Flow. For the study watershed, the water budget equation is defined as inflow equals

East Maui Watershed Partnership Adapted from Utah State University and University of Wisconsin Ground Water Project Ages 7 th -Adult

Illinois State Water Survey Division

Coffin Bay Groundwater Lens A. Resource Assessment

Drilling for Natural Gas in the Marcellus and Utica Shales: Environmental Regulatory Basics

Transcription:

water Review Monitoring and Management of Karstic Coastal Groundwater in a Changing Environment (Southern Italy): A Review of a Regional Experience Maurizio Polemio Italian National Research Council-Research Institute for Geo-Hydrological Protection (CNR-IRPI), Via Amendola 122/I, Bari 70125, Italy; m.polemio@ba.irpi.cnr.it; Tel.: +39-80-5929584 Academic Editors: Robert Puls and Robert Powell Received: 21 December 2015; Accepted: 6 April 2016; Published: 13 April 2016 Abstract: The population concentration in coastal areas and the increase of groundwater discharge in combination with the peculiarities of karstic coastal aquifers constitute a huge worldwide problem, which is particularly relevant for coastal aquifers of the Mediterranean basin. This paper offers a review of scientific activities realized to pursue the optimal utilization of Apulian coastal groundwater. Apulia, with a coastline extending for over 800 km, is the Italian region with the largest coastal karst aquifers. Apulian aquifers have suffered both in terms of water quality and quantity. Some regional regulations were implemented from the 1970s with the purpose of controlling the number of wells, well locations, and well discharge. The practical effects of these management criteria, the temporal and spatial trend of recharge, groundwater quality, and seawater intrusion effects are discussed based on long-term monitoring. The efficacy of existing management tools and the development of predictive scenarios to identify the best way to reconcile irrigation and demands for high-quality drinking water have been pursued in a selected area. The Salento peninsula was selected as the Apulian aquifer portion exposed to the highest risk of quality degradation due to seawater intrusion. The capability of large-scale numerical models in groundwater management was tested, particularly for achieving forecast scenarios to evaluate the impacts of climate change on groundwater resources. The results show qualitative and quantitative groundwater trends from 1930 to 2060 and emphasize the substantial decrease of the piezometric level and a serious worsening of groundwater salinization due to seawater intrusion. Keywords: monitoring; groundwater management; karstic aquifer; seawater intrusion; numerical modeling; climate change 1. Introduction In focusing on high-quality groundwater resources, which are larger than superficial water resources but less renewable, the main risks of degradation can be categorized into two types (Figure 1). The first type, quality degradation, is partly caused by intrinsic aquifer characteristics, which determine the aquifer vulnerability, and also by the existence of pollution sources due to anthropization as well as sources of salinization, which, though natural, can lead to the phenomenon of seawater intrusion [1]. Quality degradation causes the reduction of high-quality groundwater. The primary effect of aquifer overexploitation, i.e., the piezometric lowering, which is mainly due to increasing water demand, droughts, and/or climate change, is amplified in coastal aquifers in terms of quality degradation due to seawater intrusion and consequent salinization, bringing together quality and quantity degradation. As a result of seawater intrusion, wells could be abandoned and/or large volumes of fresh groundwater could be salinized to the point that it becomes dangerous for a number of uses, causing reduction of both groundwater quality and availability [2]. Water 2016, 8, 148; doi:10.3390/w8040148 www.mdpi.com/journal/water

Water 2016, 8, 148 2 of 16 Karstic rock outcrops cover approximately 12% of the world s surface. Karstic aquifers are highly vulnerable to contamination and anthropogenic modifications [3]. Water 2016, 8, 148 2 of 16 Karstic aquifers are a key global source of water resources; approximately 20% 25% of the global population depends on karstic aquifer groundwater [4]. In addition, 29% of European2 of 16 Water 2016, 8, 148 Karstic rock outcrops cover approximately 12% of the world s surface. Karstic aquifers are groundwater is obtained from karst aquifers; the majority of these aquifers are coastal aquifers [5,6]. highly vulnerable to contamination and anthropogenic modifications [3]. Karstic aquifers are a key global source of water resources; approximately 20% 25% of the global population depends on karstic aquifer groundwater [4]. In addition, 29% of European groundwater is obtained from karst aquifers; the majority of these aquifers are coastal aquifers [5,6]. Figure 1. Schematicofofthe themain main risks risks of degradation. Figure 1. Schematic ofgroundwater groundwater degradation. Groundwater discharge increased over the second half of the 20th century; as a consequence, a Karstic rock outcrops cover approximately 12% of the Karsticinaquifers are highly great number of aquifers are currently overexploited. Thisworld s overusesurface. is concentrated coastal areas, vulnerable to contamination and anthropogenic modifications [3]. where increasing population, growth of urban areas, and increases of irrigation and industrial Karstic aquifers are aare keyoccurring global source demands and tourism [1,6 8].of water resources; approximately 20% 25% of the global Most of thefigure aquifers of such as of Spain, France,groundwater Portugal, population depends oncoastal karstic aquifer [4]. In addition, 29% European is 1. Schematic groundwater of Mediterranean the main risks ofcountries, groundwater degradation. Slovenia, Greece, Albania, Turkey, Egypt, Morocco, Tunisia, Syria, Israel, and Italy, obtained fromcroatia, karst aquifers; the majority of Cyprus, these aquifers are coastal aquifers [5,6]. are at risk of overexploitation and seawater intrusion (Figure 2). Groundwater discharge increased over the second half of the 20th century; as a consequence, a This overuse overuse is is concentrated concentrated in in coastal coastal areas, areas, great number of aquifers are currently overexploited. This where increasing population, growth of urban areas, and increases of irrigation and industrial demands increasing population, growth of urban areas, and increases of irrigation and industrial and tourism occurring [1,6 8]. [1,6 8]. demands andare tourism are occurring aquifers of Mediterranean countries, such as Spain, Slovenia, Most of ofthe thecoastal coastal aquifers of Mediterranean countries, such as France, Spain, Portugal, France, Portugal, Croatia, Greece, Albania, Cyprus, Egypt, Morocco, Syria, Israel,Syria, and Italy, areand at risk of Slovenia, Croatia, Greece,Turkey, Albania, Turkey, Cyprus, Egypt,Tunisia, Morocco, Tunisia, Israel, Italy, overexploitation and seawaterand intrusion (Figure 2). (Figure 2). are at risk of overexploitation seawater intrusion Figure 2. The main Mediterranean coastal aquifers highly affected by groundwater quality degradation due to seawater intrusion and anthropogenic contamination. In addition to water demand problems, any modification of the quality and quantity of coastal karstic groundwater outflow can severely affect the hydrological and ecological equilibrium of a coastal water body, including highly vulnerable wetlands [9,10]. Figure 2.2.The The main Mediterranean coastal aquifers highly by affected by groundwater quality main Mediterranean coastal aquifers highly affected groundwater quality degradation due to seawater andintrusion anthropogenic contamination. degradation dueintrusion to seawater and anthropogenic contamination. In In addition addition to to water water demand demand problems, problems, any any modification modification of of the the quality quality and and quantity quantity of of coastal coastal karstic groundwater outflow can severely affect the hydrological and ecological equilibrium karstic groundwater outflow can severely affect the hydrological and ecological equilibrium of of aa coastal water body, including highly vulnerable wetlands [9,10]. coastal water body, including highly vulnerable wetlands [9,10].

Water 2016, 8, 148 3 of 16 Water 2016, 8, 148 3 of 16 In previous In previous decades, decades, climate climate variations and droughts in in the the Mediterranean Mediterranean contributed contributed to the to the increasing increasing overexploitation trend, trend, with with negative negative consequences in terms in terms of seawater of seawater intrusion intrusion [6,7,11 14]. These trends [6,7,11 14]. werethese alsotrends observed were in also Italy, observed where in seawater Italy, where intrusion seawater is intrusion the main is the cause main of cause groundwater of quality groundwater degradationquality in coastal degradation karst aquifers, in coastal the karst largest aquifers, of which the largest are located of which inare thelocated Apulian in region the [15] Apulian region [15] (Figure 3). Climate change, particularly with regard to changes in temperature, (Figure 3). Climate change, particularly with regard to changes in temperature, precipitation, sea precipitation, sea level, and sea water salinity, is a key factor for future decades, potentially level, and worsening sea water the risks salinity, of groundwater is a key factor degradation, for future as in the decades, case of potentially coastal Mediterranean worsening aquifers the risks of groundwater [16,17]. degradation, as in the case of coastal Mediterranean aquifers [16,17]. Figure 3. Apulian simplified geological scheme and hydrogeological structures (after [18,19]). Figure 3. Apulian simplified geological scheme and hydrogeological structures (after [18,19]). (A) Apulia together with the other Italian regions; (B) geological scheme and hydrogeological (A) Apulia structure together map (light with blue the other is used Italian for surface regions; water); (B) (C) geological schematic section scheme (note: andnot hydrogeological to scale; faults are structure map (light plotted blue with is used red lines; forfresh surface and water); salt groundwater (C) schematic are schematically section (note: distinguished not to scale; by light faults blue and are plotted with redlight lines; violet fresh shading, and salt respectively). groundwater are schematically distinguished by light blue and light violet shading, respectively). In Apulia, the widespread karst makes this region extremely poor in surface water but rich in groundwater resources, which have proven to be of strategic importance for economic and social development, as they can meet the high water demand of activities such as agriculture and tourism. If the variations of natural conditions that determine the recharge, boundary conditions, and growing water demand are not properly considered for management purposes, this will greatly contribute to the degradation of the quality and quantity of Apulian water resources; this trend still continues [12,20 23]. The implementation of new management strategies of groundwater resources has become extremely important [15,23,24]. These strategies have to be validated with the use of numerical models of key processes [2,16]. In Apulia, the widespread karst makes this region extremely poor in surface water but rich in groundwater resources, which have proven to be of strategic importance for economic and social development, as they can meet the high water demand of activities such as agriculture and tourism. If the variations of natural conditions that determine the recharge, boundary conditions, and growing water demand are not properly considered for management purposes, this will greatly contribute to the degradation of the quality and quantity of Apulian water resources; this trend still continues [12,20 23]. The implementation of new management strategies of groundwater resources has become extremely important [15,23,24]. These strategies have to be validated with the use of numerical models of key processes [2,16]. The scientific activity performed in recent years by the Hydrogeology Group of CNR IRPI (http://hydrogeology.ba.irpi.cnr.it/) aims to develop methodologies and knowledge for the management and safeguard of groundwater resources. The rapid national socioeconomic growth,

Water 2016, 8, 148 4 of 16 which has occurred over the past few decades, has deeply affected the main hydrogeological systems, thereby leading to different risk conditions, especially for the largest karst coastal aquifers, located in the Apulian region. The approach and methodologies described in this paper, which are based on Apulian aquifers as a case study, include a recharge and discharge trend analysis, a regional monitoring network, multi-parameter well logging for rapid groundwater quality classification, a discussion of regulation of groundwater utilization, a salinity threshold and multi-temporal spatial analysis, and large-scale numerical modeling. 2. Hydrogeological Features of the Main Aquifers Bounded by the Appenninic chain on the western side, Apulia is bounded for more than 800 km by coastline on the Adriatic and Ionian Seas (Figure 3). The Apulian carbonate platform includes large limestone and dolomite Mesozoic outcrops, which constitute the physiographic units of Gargano, Murgia, and Salento. Outcrops of detrital organogenic soils and rocks (Tertiary and Quaternary) can be found in the topographic troughs that partially overlap the abovementioned platform, as in the case of the Tavoliere physiographic unit [18,19]. Four hydrogeological structures (HSs) can be distinguished in the Apulian region: Tavoliere, Gargano, Murgia, and Salento, the last three of which are karstic. The Apulia region hosts the largest coastal karstic Italian aquifers. The Tavoliere HS is mainly characterized by a porous shallow aquifer, whose groundwater flow is limited by a clay bottom a few hundred meters thick, both lying on top of the carbonate platform. Groundwater flows in phreatic conditions upstream, far from the sea, and is confined downstream, close to the coast. Inland, the piezometric head is approximately 300 masl; only along a very narrow coastal strip is the aquifer sufficiently deep to allow seawater intrusion. Salinization is a negligible phenomenon, both from a scientific and practical point of view (also due to the local availability of surface water resources), for this aquifer. The Tavoliere HS hosts aquifers of secondary relevance, whereas the karstic HSs of Gargano, Murgia, and Salento create large carbonate aquifers that are hundreds of meters deep and are deeply influenced by karstic processes [19]. These aquifers exhibit varying degrees of fracturing and permeability [18]. Gargano is a high horst that strikes NE SW in the shape of a promontory extending into the Adriatic Sea (peak altitude 1065 masl). The carbonate rock types of Gargano allow it to be divided into two distinct parts: the eastern sector consists of frontal reef depositions and transition depositions, and the central-western sector is dominated by reef facies and, to a lesser degree, by back reef depositions. Due to the mountainous features of Gargano, the level of anthropogenic activity is very low in terms of the degree of groundwater utilization [16,24]. The Murgia HS, which corresponds to a plateau (peak altitude 680 masl) and has a large asymmetric horst caged by two direct fault systems (striking NW SE and NE SW), and Salento HS, an almost flat area (peak altitude 150 180 masl) bounded by two seas, together constitute a continuum for groundwater flow. These HSs are also divided by a morphological-structural feature called the Messapian Threshold, which is multiples of kilometers wide and covers an area that extends from sea to sea (Figure 3) [18]. Moving from Murgia to Salento HS, the hydrogeological parameters improve (in the latter, higher hydraulic conductivity and storage coefficients are observed), the depth to water decreases considerably (from hundreds to tens of meters), the mean piezometric head decreases, the mean distance from the sea decreases, the level of groundwater utilization increases, and, on the whole, the risk of quality degradation as a result of seawater intrusion increases [12,15,22]. Recharge of both of these HSs occurs mainly from rainfall infiltration, but a large amount of the Murgia recharge discharges to the sea by flowing across the Salento HS [18,19,25]. Focusing on groundwater flow and the piezometric surface, a hydrogeological boundary divides the Salento HS into two portions: the northern portion, in which the main inflow is due to recharge from the Murgia HS and discharge to the Ionian and Adriatic Seas, together with direct rainfall infiltration, and the

Water 2016, Water 8, 2016, 148, 148 5 of 165 of 16 into two portions: the northern portion, in which the main inflow is due to recharge from the Murgia southern HS and portion discharge (SS-HS), to the inionian whichand the Adriatic rechargeseas, is only together due towith direct direct rainfall infiltration. infiltration, This and southern portion southern of portion Salento(SS-HS), area experiences in which greater the recharge effects is of only seawater due to intrusion direct rainfall [22]. infiltration. This southern portion of the Salento area experiences greater effects of seawater intrusion [22]. 3. Recharge and Discharge Trends 3. Recharge and Discharge Trends Monthly time series of rainfall and temperature from 1916 (temperature from 1925) to 2005, based Monthly time series of rainfall and temperature from 1916 (temperature from 1925) to 2005, on data published by the National and Regional Hydrological Service, were analyzed. Mean annual based on data published by the National and Regional Hydrological Service, were analyzed. Mean rainfall varies from 430 to 1170 mm in the Apulia region, with an average value of 644 mm [11]. Over annual rainfall varies from 430 to 1170 mm in the Apulia region, with an average value of 644 mm the entire study period, the Apulian annual rainfall trend is equal to 0.80 mm/year [26]. [11]. Over the entire study period, the Apulian annual rainfall trend is equal to 0.80 mm/year [26]. Mean Mean annual annual net net rainfall rainfall varies varies from from 52 52 to 675 mm mmin inthe theregion, region, with with a regional a regional average average of of 146 mm. 146 mm. For For each each time time series, the net rainfall trend trend was was strongly strongly negative; negative; the absolute the absolute value of value net of net rainfall trend got worse more than than rainfall rainfall trend. trend. The net The rainfall net rainfall trend ranged trend from ranged 0.23 from to 3.52 0.23 to 3.52mm/year. In the In whole the whole period, period, the decrease the decrease of net rainfall of netamounts rainfall to amounts 42.2% of to the 42.2% mean of net the rainfall. mean net rainfall. The The five-year five-year moving moving average average of annual of annual net rainfall net rainfall was assessed was assessed (Figure (Figure 4). Negative 4). Negative trends are trends are evident for all spatial domains; an an anomalous anomalous sequence sequence of negative of negative values values is evident is evident since 1980. since 1980. Figure Figure 4. Normalized 4. deviation of of the five-year moving averageof of annual net net rainfall rainfall (continuous lines) lines) and and linear linear trend trend (dashed lines) assessed for the main Apulian physiographic units. units. A unique A unique color color was was used used for each for each physiographic unit unit to to plot the time seriesand andthe the trend trend line. line. A few deep wells existed in Apulia, located in Tavoliere, in 1936 [27]. The level of Apulian Agroundwater few deep utilization wells existed was very in Apulia, low (zero located in large inareas) Tavoliere, until the in end 1936 of [27]. the Second The level World ofwar. Apulian groundwater A continuously utilization increasing was trend very low of groundwater (zero in large abstraction, areas) until starting thein end the of second the Second half of the World 1950s, War. A continuously was observed. increasing The current trend of well groundwater number is not abstraction, known due starting to the high in the percentage second half of unauthorized of the 1950s, was observed. wells; The however, current this well number number is approximated is not known to be due many to the tens high of thousands. percentage The of continued unauthorized effect of wells; however, recharge this variability number isand approximated increasing exploitation to be manycan tens be of assessed thousands. in terms Theof continued groundwater effect quantity of recharge variability degradation and increasing by measurement exploitation of the piezometric can be assessed head and in terms coastal of spring groundwater discharge quantity [28]. degradation A regional piezometric study was realized for the period 1965 2010, using monthly data, with by measurement of the piezometric head and coastal spring discharge [28]. gaps of 30 selected well time series (Figure 5) [29]. A regional piezometric study was realized for the period 1965 2010, using monthly data, with gaps of 30 selected well time series (Figure 5) [29]. The statistical reliability of the detected trends was tested with the Mann Kendall test at a confidence level of 0.05; the piezometric trend was quantified with the Angular Coefficient (AC) of the straight-line regression. The piezometric trend was generally downward (negative AC) because there was a widespread tendency, albeit in some cases a very slow one, towards a piezometric decline (Table 1). These results were confirmed by measurements of coastal spring outflow. In the case of the Fiume Grande Spring (Figure 5), data with gaps are available from 1926 to 2010: on the basis of a mean value of 0.574 m 3 /s, the trend was equal to 0.00443 m 3 /year [29].

Water 2016, 8, 148 6 of 16 Water 2016, 8, 148 6 of 16 Figure Figure 5. Location 5. Location map map of the of the selected wells and springs (x (xand and y y coordinates were were plotted plotted as meters). as meters). The statistical reliability of the detected trends was tested with the Mann Kendall test at a Table 1. Statistical values of the selected piezometric time series (AC: angular coefficient of straight confidence level of 0.05; the piezometric trend was quantified with the Angular Coefficient (AC) of linethe trend). straight-line regression. The piezometric trend was generally downward (negative AC) because there was a widespread tendency, albeit in some cases a very slow one, towards a piezometric decline (Table 1). These results Period were confirmed by measurements Value of coastal (masl) spring outflow. In AC the Well No. case of the Fiume Grande from Spring (Figure 5), todata with gaps Min are available Mean from 1926 Max to 2010: on (m/y) the basis of a mean value of 0.574 m 3 /s, the trend was equal to 0.00443 m 3 /year [29]. 1 September 1973 December 2009 0.81 1.45 1.64 0.0007 2 Table 1. Statistical September values 1973 of the selected December piezometric 2009 time series 3.81(AC: angular 5.98 coefficient 6.58 of straight 0.0045 3 line trend). December 1965 December 2009 1.61 0.82 2.02 0.0017 4 June 1975 December 2009 4.95 29.60 63.54 0.1004 8 June 1975 Period December 2009 42.79 Value 47.22 (masl) 54.36 AC 0.0039 Well No. 9 May 1975 from Decemberto 2009 31.87 Min Mean 44.64 Max 50.99 (m/y) 0.0083 19 1 September September 1973 1973 September December 2003 2009 1.51 0.81 1.45 1.99 1.642.20 0.0007 0.0120 20 2 October September 1973 1973 September December 2003 2009 1.95 3.81 5.98 2.30 6.586.02 0.0045 0.0300 21 3 October December 1973 1965 September December 2003 2009 0.45 1.61 0.82 1.42 2.021.79 0.0017 0.0132 22 4 September June 1973 1975 September December 2003 2009 0.99 4.95 29.60 2.07 63.54 2.40 0.1004 0.0192 23 8 September June 1973 1975 September December 2003 2009 24.55 42.79 47.22 27.45 54.36 28.18 0.0039 0.0492 32 9 OctoberMay 1975 1975 February December 20102009 1.23 31.87 44.64 4.78 50.99 9.27 0.0083 0.0035 38 19 September September 1973 1973 September September 2003 2003 0.11 1.51 1.99 1.68 2.20 2.93 0.0120 0.0180 39 September 1973 February 2010 0.65 0.13 0.92 0.0003 20 October 1973 September 2003 1.95 2.30 6.02 0.0300 41 October 1973 January 2010 6.40 1.42 1.85 0.0044 21 October 1973 September 2003 0.45 1.42 1.79 0.0132 43 September 1973 January 2010 1.10 2.21 2.94 0.0035 22 September 1973 September 2003 0.99 2.07 2.40 0.0192 44 September 1973 March 2010 1.44 3.26 3.73 0.0041 23 September 1973 September 2003 24.55 27.45 28.18 0.0492 45 October 1973 March 2010 1.01 2.12 2.41 0.0023 32 October 1975 February 2010 1.23 4.78 9.27 0.0035 46 October 1973 March 2010 3.38 4.63 5.17 0.0026 38 September 1973 September 2003 0.11 1.68 2.93 0.0180 47 January 1978 January 2010 0.18 1.11 1.99 0.003 48 39 July September 1968 1973 January February 2010 2010 1.17 0.65 0.13 2.24 0.92 2.61 0.0003 0.002 49 41 JuneOctober 1975 1973 November January 1996 2010 4.46 6.40 1.42 4.84 1.855.58 0.0044 0.023 50 43 October September 1973 1973 February January 2010 2010 2.21 1.10 2.21 3.19 2.943.62 0.0035 0.001 51 44 July September 1975 1973 February March 2010 2010 1.52 1.44 3.26 2.80 3.733.56 0.0041 0.003 52 September 1973 February 2010 7.40 14.72 16.03 0.019 54 March 1973 February 2010 2.82 4.87 6.56 0.003 55 September 1973 February 2010 1.13 2.28 2.68 0.002 56 September 1973 November 2009 2.36 6.36 8.70 0.012 119 May 1975 March 2010 17.94 23.50 38.31 0.010 125 January 1965 February 2010 1.00 1.99 3.24 0.001

Water 2016, 8, 148 7 of 16 The most probable spatial piezometric trend of each HS ranges from a low decrease, as in the case of Gargano (well 32), to a moderate decrease, as in the case of Salento, and to a high decrease rate, as in the case of Murgia and Tavoliere. The widespread trend towards piezometric lowering increases the risk of salt pollution by seawater intrusion, especially in Salento, which is higher than in the rest of the region, due to peculiarities of the boundaries and the piezometric range (the Salento historical maximum piezometric head of approximately 4 m asl is orders of magnitude lower than maxima of Murgia and Gargano). 4. Groundwater Regional Monitoring Network To characterize the modification in terms of groundwater availability and quality due to salinization and/or human-related pollution, a regionally based continuously operating hydrogeological monitoring network should be considered an affordable option [12,23]. In 1998, the Apulian monitoring system used 118 wells, selected from 120 available wells, for monitoring purposes (the well map is shown in Figure 6), some of which were hundreds of meters deep [30]. Different types of wells were equipped with electrical sensors along the vertical axis and connected to a GIS system permitting piezometric (P, measured in each well), temperature (T), electric conductivity (EC), total dissolved solid (TDS, estimated on the basis of T and EC measurements), Eh, and dissolved oxygen Water (DO) 2016, measurements. 8, 148 8 of 16 Figure 6. Figure Monitoring 6. Monitoring wells wells and and log log types. types. Log Log type and selected examples of of wells wells and and logs logs (from (from A A to F in Figures to F in 6Figures and 7): 6 (1) and type 7): (1) A; type (2) A; type (2) type B; (3) B; type (3) type C; C; (4) (4) type D; (5) typee; E; (6) (6) type type F. F. 5. Multi-Parameter Logging for Rapid Groundwater Quality Classification Three types of wells were used. Piezometric wells, in which only a P sensor was installed, The multi-parameter logs are measurements of some groundwater parameters, generally T, EC, accounted TDS, for ph, 78 Eh, wells and (66% DO, along of the the total); water three columns couples of the (at wells. different The rapid depths) groundwater of T and quality EC sensors (permitting classification TDS assessment) method is based wereon installed detection in with 12multi-parameter wells (10%), i.e., logs the of typical coastal parameter wells, located trends along the coast that and are generally due to natural very hydrogeological deep (the well conditions bottomof is the below studied the aquifer. transition zone between fresh and saline groundwater); The method T, was EC, tested ph, Eh, in 120 and Apulian DO sensors wells over were one installed year, with insurveys quality conducted control every wells, three accounting for 24% months of the total, (Figure which 6) [31]. were It was located observed where that these the withdrawal types of trends rate are and rather groundwater recurrent in space contamination and time (Figure 7). Six typical trends of the multi-parameter logs were distinguished: the inner or hazards were higher. recharge area (type A), in which pure fresh groundwater, typical of a carbonate aquifer, flows; the Oncoastal the basis strip of trend the(type monitoring B), characterized network by a data, salinity high knee differences trend due were to seawater assessed intrusion; in terms the of the seawater Tavoliere intrusion shallow in the aquifer karstic (type structures. C), in which The the lateral low groundwater seawaterquality intrusion is exhibited in Murgia by the isph, bounded Eh, in a narrow and coastal DO logs; strip, the usually transition lesszone thanbetween 3 km in the width Murgia and 14 the kmtavoliere at maximum; (type D); apart the transition from this zone strip, the groundwater between salinity the Tavoliere is lowand further Gargano inland (type (less E), than in which 0.5 g/l), the trends withchange a depth from ofone 500 900 aquifer mto the below sea other; and the transition zone between the Murgia and the Salento (type F), in which the effect of the level. In Salento, only narrow areas, located inland, exhibit pure fresh groundwater. leakage from inland Murgia to Salento corresponds to a steady and widespread decrease of the piezometric head and a significant increase in permeability. It is reasonable that an extensive use of the suggested method can be easily applied to the preliminary detection of the hydrogeological conditions. Some examples include preliminary and rapid assessment (which can be further developed with more expensive hydrogeological studies) of the natural groundwater quality; if aquifer boundaries are uncertain, the preliminary recognition of

Water 2016, 8, 148 8 of 16 5. Multi-Parameter Logging for Rapid Groundwater Quality Classification The multi-parameter logs are measurements of some groundwater parameters, generally T, EC, TDS, ph, Eh, and DO, along the water columns of the wells. The rapid groundwater quality classification method is based on detection with multi-parameter logs of typical parameter trends that are due to natural hydrogeological conditions of the studied aquifer. The method was tested in 120 Apulian wells over one year, with surveys conducted every three months (Figure 6) [31]. It was observed that these types of trends are rather recurrent in space and time (Figure 7). Six typical trends of the multi-parameter logs were distinguished: the inner or recharge area (type A), in which pure fresh groundwater, typical of a carbonate aquifer, flows; the coastal strip trend (type B), characterized by a salinity knee trend due to seawater intrusion; the Tavoliere shallow aquifer (type C), in which the low groundwater quality is exhibited by the ph, Eh, and DO logs; the transition zone between the Murgia and the Tavoliere (type D); the transition zone between the Tavoliere and the Gargano (type E), in which the trends change from one aquifer to the other; and the transition zone between the Murgia and the Salento (type F), in which the effect of the leakage from inland Murgia to Salento corresponds to a steady and widespread decrease of the piezometric head and a significant increase in permeability. Water 2016, 8, 148 9 of 16 Figure 7. Examples Figure 7. Examples of multi-parameter of log log types after [31]. [31]. Parameters: (1) Temperature (1) Temperature (T, C); (T, C); (2) TDS (g/l); (3) ph; (4) dissolved oxygen (DO and (5) oxidation-reduction potential (Eh, mv). (2) TDS (g/l); (3) ph; (4) dissolved oxygen (DO mg/l); and (5) oxidation-reduction potential (Eh, mv). Log type and the selected example well (the well locations are in Figure 6): (A) Inner or recharge Log type andarea; the (B) selected coastal strip example or sea well water (the intrusion wellarea; locations (C) Tavoliere are inhs; Figure (D) transition 6): (A) zone Inner between or recharge area; (B) coastal strip Murgia or sea and water Tavoliere intrusion HSs; (E) area; transition (C) Tavoliere zone between HS; Tavoliere (D) transition and Gargano zone between HSs; and Murgia and (F) transition zone between Murgia and Salento HSs. Tavoliere HSs; (E) transition zone between Tavoliere and Gargano HSs; and (F) transition zone between Murgia and 6. Regulation Salento HSs. of Groundwater Utilization Until 1984, laws regulated groundwater extraction in Apulia (and roughly in Italy) mainly from an administrative point of view, without considering the effects of seawater intrusion and the application of hydrogeological management criteria. In 1984, the regional authorities implemented a safeguarding and decontamination plan called PRA or PLAN1, which determined the quality zonation of Apulian groundwater and the regulation of groundwater utilization as a function of the risk of groundwater degradation, mainly considering seawater intrusion effects [32]. As an example, the drilling of new pumping wells was forbidden where the salinity was too high due to the mixing with saline groundwater (no discharge zone, called NOD1 zone hereafter, Figure 8A).

Water 2016, 8, 148 9 of 16 It is reasonable that an extensive use of the suggested method can be easily applied to the preliminary detection of the hydrogeological conditions. Some examples include preliminary and rapid assessment (which can be further developed with more expensive hydrogeological studies) of the natural groundwater quality; if aquifer boundaries are uncertain, the preliminary recognition of a tapped aquifer; and the detection of anthropogenic effects on quality, i.e., the occurrence of human-related pollution or seawater intrusion. Human-related pollution is highlighted by an anomalous trend of one or generally more than one of the following parameters: T, ph, Eh, and DO. The effect of seawater intrusion (upconing or lateral intrusion) requires repeated logs over time, focusing on EC or TDS variations; in this way, overexploitation effects can be detected by observing the vertical modification of the transition zone between fresh and saline groundwater [18]. This method may prove particularly useful for practical purposes because the execution of logs is quite simple, quick, and comes at a low cost. 6. Regulation of Groundwater Utilization Until 1984, laws regulated groundwater extraction in Apulia (and roughly in Italy) mainly from an administrative point of view, without considering the effects of seawater intrusion and the application of hydrogeological management criteria. In 1984, the regional authorities implemented a safeguarding and decontamination plan called PRA or PLAN1, which determined the quality zonation of Apulian groundwater and the regulation of groundwater utilization as a function of the risk of groundwater degradation, mainly considering seawater intrusion effects [32]. As an example, the drilling of new pumping wells was forbidden where the salinity was too high due to the mixing with Water saline2016, groundwater 8, 148 (no discharge zone, called NOD1 zone hereafter, Figure 8A). 10 of 16 Figure 8. Groundwater quality zonationfor for regulation (after [27]). [27]). (A) (A) As As regulated in in 1984 1984 (PRA (PRA or PLAN1, or PLAN1, no no more more applicable from from 2009): (1) (1) zone without discharge restrictions; (2) low-quality groundwater with salt contamination (NOD1 zone); and (3) high-quality groundwater (safeguarded for drinking use); (B) As regulated in 2009 (PTA or PLAN2, applicable from 2009): (1) zone without discharge restrictions of seawater of seawater intrusion intrusion effects; (2) effects; low-quality (2) low-quality groundwater groundwater with salt contamination with salt (new contamination fresh groundwater (new fresh abstraction groundwater is not abstraction allowed, is NOD2 not allowed, zone); NOD2 (3) qualitative zone); (3) and qualitative quantitative and protection quantitative zone protection (QQP zone, regulated (QQP zone, newregulated pumping new wells); pumping and (4) quantitative wells); and (4) protection quantitative zone (QP protection zone, new zone abstraction (QP zone, new permits abstraction suspended). permits suspended). The areas where the water quality was the lowest were located along the coast, as an effect of seawater intrusion. The boring of new abstraction wells was regulated where groundwater quality was suitable for drinking utilization (Figure 8A). This was roughly the case for the main recharge areas of the Apulian aquifers. In 2009, the Water Protection Plan (PTA or PLAN2), still applicable, replaced PLAN1, and considers other phenomena in addition to seawater intrusion [33]. It distinguishes zones of quantitative protection (QP zone), as in the example of the shallow Tavoliere aquifer, the most overexploited Apulian aquifer, in which new abstraction permissions were suspended (Figure 8B)

Water 2016, 8, 148 10 of 16 was suitable for drinking utilization (Figure 8A). This was roughly the case for the main recharge areas of the Apulian aquifers. In 2009, the Water Protection Plan (PTA or PLAN2), still applicable, replaced PLAN1, and considers other phenomena in addition to seawater intrusion [33]. It distinguishes zones of quantitative protection (QP zone), as in the example of the shallow Tavoliere aquifer, the most overexploited Apulian aquifer, in which new abstraction permissions were suspended (Figure 8B) [34]. Focusing on seawater intrusion, zones affected by seawater intrusion (hereinafter referred to as the NOD2 zone, to be distinguished from the NOD1 zone of the previous PLAN1), located along the coast, and zones defined as qualitative and quantitative protection zones (QQP zone), located between recharge or inland areas, are distinguished from each other. In the NOD2 zone, the issue of permits for constructing new wells and for abstraction of fresh groundwater is suspended, even if the probability of finding fresh groundwater in this zone is low. It is only possible to drill for new pumping wells in the case of using deep and saline groundwater for specific purposes, but a specific hydrogeological study and assessment are required. The authorization renewal of pre-existing wells is only possible on the condition of two criteria: the well bottom (m below sea level) should be less than 30, 25, or 20 times the piezometric head (m above sea level) in the cases of Gargano, Murgia, and Salento, respectively; additionally, the maximum abstraction rate should cause a drawdown (m) of less than 50% (Murgia) or 30% (Gargano and Salento) of the piezometric head (m asl). In the QQP zone, new authorizations are again regulated considering the well bottom (a condition quite similar to the NOD2 zone but declined if groundwater is locally confined) and the maximum abstraction allowed. It requires a drawdown (m) of less than 60% (Murgia) or 30% (Salento) of the piezometric head (masl); TDS should be less than 1 g/l, and chloride concentrations should be less than 500 mg/l. Both plans do not explain the criteria defined to plot NOD1, NOD2, and QQP zones, but it seems the high level of salinity and pumping well density were considered as key criteria. Neither plan considered the spatial and temporal modification of seawater intrusion effects. The PLAN1 approach failed to protect groundwater quality and availability. It is too soon to assess the effects of PLAN2, but the described criteria seem to be overly linked to a single-well criterion; for example, the Ghyben Herzberg principle of a sharp interface does not consider the overlapping effects of nearby wells and the whole amount of abstraction compared to the variable recharge. 7. Salinity Threshold and the Multi-Temporal Spatial Analysis Method A very simple criterion to detect the absence or presence of seawater intrusion quality degradation was defined for the Apulian karstic aquifers [23]: the threshold criterion. The key concept was to define a boundary or threshold between pure fresh groundwater and any type of mixing between fresh and saline groundwater. On the basis of a hydrogeochemical classification of approximately 500 groundwater laboratory analyses (samples taken from the karstic HSs), groundwater samples of the hydrochemical parameters (water types), i.e., Ca 2+ -HCO 3, Ca 2+ -Mg 2+ -HCO 3, and Mg 2+ -Ca 2+ -HCO 3, were considered to be representative of pure fresh karstic groundwater, as an effect of chemical modification on infiltrated rain water due to the characteristics of the topsoil and the karstic aquifer. The simplest grouping of samples is as follows: pure fresh groundwater (group F) and the remaining samples (group S), consisting of pure fresh groundwater mixed with variable percentages of seawater (all the pure fresh groundwater samples came from wells located in the inner portion of Murgia and Salento). On the basis of salinity statistical analysis of group F, the upper salinity limit (TDS) or threshold value of pure fresh groundwater was determined as equal to 0.5 g/l. With this threshold, the analysis of the salinity spatial evolution over time was based on low frequency and high-density salinity data used for plotting the threshold contour line (TCL) of 1981, 1989, 1997, and 2003, using all available data, including those for private discharge wells [23,27].

(group S), consisting of pure fresh groundwater mixed with variable percentages of seawater (all the pure fresh groundwater samples came from wells located in the inner portion of Murgia and Salento). On the basis of salinity statistical analysis of group F, the upper salinity limit (TDS) or threshold value of pure fresh groundwater was determined as equal to 0.5 g/l. With this threshold, the analysis of the salinity spatial evolution over time was based on low Water 2016, 8, 148 11 of 16 frequency and high-density salinity data used for plotting the threshold contour line (TCL) of 1981, 1989, 1997, and 2003, using all available data, including those for private discharge wells [23,27]. Multi-temporal spatial analysis was was then then performed performed on the basis on the of TCLs basis overlapping, of TCLs overlapping, considering the considering time variations the time from variations 1981 tofrom 20031981 (Figure to 2003 9; the (Figure 2003 TCL 9; the was 2003 not TCL plotted was not to simplify plotted to the simplify figure reading). the figure Three reading). zones Three of the zones contour of the line contour were distinguished. line were distinguished. Figure 9. 9. Threshold (0.5 (0.5 g/l) g/l) and and salinity salinity vulnerability map: map: (1) 1997 (1) contour; 1997 contour; (2) 1989(2) contour; 1989 contour; (3) 1981 contour; (3) 1981 contour; (4) inland(4) zone inland with zone salinity with lower salinity than lower thethan threshold; the threshold; (5) coastal (5) zone coastal with zone salinity with salinity greater than greater thethan threshold; the threshold; and (6) intermediate and (6) intermediate zone (spanned zone (spanned by the threshold by the contour threshold line contour over time). line over time). In the first zone or inland zone, salinity remained permanently below the threshold, indicating a low vulnerability In the first zone to seawater or inland intrusion zone, salinity [18]. TCLs remained were always permanently downward below orthe onthreshold, the boundary indicating of this zone a low (between vulnerability the zone to seawater and theintrusion coast). In[18]. the second TCLs were or coastal always zone, downward salinity or was on always the boundary above the of threshold due to seawater intrusion, showing a high vulnerability to seawater intrusion; TCLs were always upward of this zone or on the boundary of this zone. In the third zone or intermediate zone, each point was spanned by TCLs over time. The largest inland movement of TCL (worst salinization effects) was observed as an effect of huge drought periods, as in the case of the 1989 TCL, whereas a salinity or quality improvement was an effect of rainy periods, as in the case of the generalized seaward movement of the 1997 TCL [11,26,27]. The first zone corresponds to a wide portion of inland Murgia and a narrow strip in the middle of the Salento peninsula; it had not been involved in seawater intrusion. The groundwater salinization can be considered a long-standing phenomenon in the second zone. In the third zone, salinity is highly sensitive to climate, the water cycle, exploitation variability and, predominantly, human management of groundwater resources. This zonation was preliminarily validated using all available previous data, similar to how the analyzed groundwater samples were used to calculate the threshold value. Data of sufficient density to plot contour lines are not available after 2004, but available scattered data confirm the validity of the zonation approach. Based on these results, some basic management criteria could be proposed. As an example, discharge for drinking purposes could be increased in the inland zone considering a limitation rule based on the annual recharge (in any case, almost less than total recharge), operating on the basis of both the yearly recharge values and the statistical analysis of low recharge values. The discharge in the intermediate zone could be managed using two criteria: the discharge of the intermediate zone should be reduced or forced to be considerably lower in the residual part (not

Water 2016, 8, 148 12 of 16 used in the inland zone) of the annual recharge and validated by monitoring salinity effects with the purpose of ensuring a buffer zone around the inland zone, in which the salinity should be maintained close to the threshold. The discharge in the coastal zone should be authorized only for specific uses for which brackish or saline water is acceptable, as long as it does not cause negative economic or environmental effects. Saline water discharge should be encouraged if it entails positive effects on the pure fresh groundwater resources. In both cases, the effects of discharge in terms of degradation risks should be assessed in detail. The discharge plan should be defined taking into account different discharge and recharge scenarios. In the intermediate and coastal zones, the discharge wells (number, location, geometry, and construction techniques) should be optimized on the basis of local hydrogeological characteristics. Each one of these rather simple proposals are still neglected by regulations. 8. Large-Scale Numerical Model Approach In the previous section, basic management criteria were proposed on the basis of the TCL approach. The selected management criteria, which could be at the core of a future groundwater protection plan, should be defined by the management authority to determine measureable groundwater quality and quantity objectives (milestones) and conduct pre-validation and ongoing post-validation surveys. The pre-validation should be realized before the plan is applicable using numerical modeling to test the feasibility of the planned objectives based on future scenarios; the post-validation should be based on monitoring to check the ongoing feasibility of the planned objectives without plan adjustment. Pre-validation and post-validation were not realized in the case of PLAN1 and PLAN2 by management authorities. Pre-validation in the case of large, deep, coastal, and karstic aquifers, which represent a very complex hydrogeological problem, is uncommon. This research activity started with the purpose of demonstrating the feasibility of conducting pre-validation of the Apulian HSs by testing the future effects of current management criteria. For this purpose, SS-HS, the area undergoing the most significant degradation due to seawater intrusion, i.e., southern Salento, was selected. An equivalent porous continuous medium (EPCM) was selected to be applied to SS-HS, as described by [16]. The numerical codes used were MODFLOW [35] and SEAWAT [36] (the former was basically used in the initial phases for steady-state calculations and initial calibrations; the latter, which is a computer program designed to simulate three-dimensional variable-density groundwater flow, was used in the rest of the activities). The active domain of the study area (active cells) covered approximately 2300 km 2 with 45,925 cells. Vertically, the area was divided into 12 layers, from 214 to 350 m asl, to allow an accurate lithological and hydrogeological discretization. The thickness and geometry of the layers were defined based on 3D knowledge of the hydrogeological complexes. The natural flow and salinity conditions were calculated using a steady state model, calibrated with mean historical data of the 1930s, a period in which the level of discharge was so low that it was considered negligible; these steady-state results could be considered as natural groundwater conditions and are used as a reference to assess the following modifications. Piezometric and salinity variations from 1930 to 2060 were simulated under three past scenarios (up to 1999) and three future scenarios of climate change to test different hypotheses regarding discharge for different types of utilization, and changes in sea level and salinity. The model was validated using surveyed piezometric and salinity data of different decades. The past scenario results confirmed the piezometric lowering and the salinization worsening not only due to the lateral seawater intrusion but mainly due to upconing caused by groups of drinking wells [2,37]. Based on the reliable A1B scenario model [13], rainfall was predicted to decrease: equal to 3.9%, 5.9%, and 9% of the 1960 1980 mean value for the periods of 2000 2020 (future scenario F1), 2021 2040 (F2), and 2041 2060 (F3), respectively. Mean temperature increases, ranging from 0.9 C

Water 2016, 8, 148 13 of 16 for F1 to 2.4 C for F3, were predicted. Based on these results, the amount of annual recharge was calculated. As groundwater discharge cannot be predicted, a steady discharge level for all uses at the level of 2000 was hypothesized for each future scenario, which is almost equal to the discharge values observed from the implementation of PLAN2 (2009). This hypothesis is approximate and will likely underestimate the future well discharge if adaptation measures are not implemented, as is currently the case. Continuing the application of current management criteria, which do not consider the physical characteristics of seawater intrusion and the need of adaption measures to regulate the discharge (to be almost lower than the recharge), could cause effects similar to those observed in Figure 10 (the displayed layer was selected because it includes the transition zone between pure fresh groundwater and saline water over a large portion of the aquifer, far from the coast), in which the results of the final Water 2016, 8, 148 14 of 16 year of the future simulation are summarized. Figure 10. SS-HS piezometric (A, m); and salinity (B, mg/l) maps for a selected layer (layer 50 to Figure 10. SS-HS piezometric (A, m); and salinity (B, mg/l) maps for a selected layer (layer 50 to 65 masl) for the year 2060, assessed according to the hypothesis of an unchanged discharge level 65 masl) for the year 2060, assessed according to the hypothesis of an unchanged discharge level observed at the beginning of the PTA application [33]; piezometric head and salinity increase with observed at the beginning of the PTA application [33]; piezometric head and salinity increase with respect to natural or steady state values, based on data from the 1930s (after [16]). The light blue area respect to natural or steady state values, based on data from the 1930s (after [16]). The light blue area denotes the shallow aquifer. denotes the shallow aquifer. 9. Conclusions A dramatic piezometric decrease of more than 2.5 m was predicted in 2060 compared to steady-state or natural conditions (Figure 10A). The predicted 2060 salinity shows a huge increase from steady-state Decades of conditions, surveys and greater studies than on large 5000 mg/l karstic incoastal wide areas aquifers (Figure were 10B). summarized This increase to offer wasan predicted overview along of efforts the coastal to pursue zone, particularly the sustainable the utilization western Ionian of the area, main but natural also far resources from theof coast. a karstic In the coastal latter case, Italian the region. increase and the consequent dramatic effects on drinking water resources are due to the enlargement The main and results worsening show the ofprogressive the upconing degradation phenomenon, of Apulian which can groundwater be observedresources, around groups both in ofterms drinking of quality wells. and quantity. The methods and approaches described were designed and tested to be An proposed evident to piezometric support decision drop was makers confirmed to adopt for decisions the past; useful a similar for the dramatic purpose drop of maximizing appears to be the likely groundwater in the future. utilization The lateral to intrusion ensure the andsteady upconing availability effects ofand seawater quality intrusion of groundwater were non-negligible resources. infuture the past research and willdevelopments be considerable will in the be finalized future. All on phenomena the proposal considered of management here, including criteria, sea including level and those sea described salinity, exhibited in this paper, not negligible which are effects optimized on coastal using groundwater. the pre-validation approach, together with other potential criteria based on adaptation strategies focused on agricultural uses, the greatest type of water demand in this Mediterranean region. Acknowledgments: I would thank the whole Hydrogeological Group of CNR-IRPI, which contributed to realize the described researches, with peculiar reference to Andrea Romanazzi and Livia Zuffianò. Conflicts of Interest: The authors declare no conflict of interest. References 1. Polemio, M.; Dragone, V.; Romanazzi, A. La risorsa idrica. Sfruttamento, depauperamento dei serbatoi sotterranei e utilizzo razionale nel caso della Calabria. In L acqua in Calabria: Risorsa o Problema?; Dramis, F.,

Water 2016, 8, 148 14 of 16 9. Conclusions Decades of surveys and studies on large karstic coastal aquifers were summarized to offer an overview of efforts to pursue the sustainable utilization of the main natural resources of a karstic coastal Italian region. The main results show the progressive degradation of Apulian groundwater resources, both in terms of quality and quantity. The methods and approaches described were designed and tested to be proposed to support decision makers to adopt decisions useful for the purpose of maximizing the groundwater utilization to ensure the steady availability and quality of groundwater resources. Future research developments will be finalized on the proposal of management criteria, including those described in this paper, which are optimized using the pre-validation approach, together with other potential criteria based on adaptation strategies focused on agricultural uses, the greatest type of water demand in this Mediterranean region. Acknowledgments: I would thank the whole Hydrogeological Group of CNR-IRPI, which contributed to realize the described researches, with peculiar reference to Andrea Romanazzi and Livia Zuffianò. Conflicts of Interest: The authors declare no conflict of interest. References 1. Polemio, M.; Dragone, V.; Romanazzi, A. La risorsa idrica. Sfruttamento, depauperamento dei serbatoi sotterranei e utilizzo razionale nel caso della Calabria. In L acqua in Calabria: Risorsa o Problema?; Dramis, F., Mottana, A., Eds.; Aracne: Roma, Italy, 2013; pp. 2 29. 2. Van Camp, M.; Mtoni, Y.; Mjemah, I.C.; Bakundukize, C.; Walraevens, K. Investigating seawater intrusion due to groundwater pumping with schematic model simulations: The example of the Dar es Salaam coastal aquifer in Tanzania. J. Afr. Earth Sci. 2014, 96, 71 78. [CrossRef] 3. Polemio, M.; Casarano, D.; Limoni, P.P. Karstic aquifer vulnerability assessment methods and results at a test site (Apulia, Southern Italy). Nat. Hazards Earth Syst. Sci. 2009, 9, 1461 1470. [CrossRef] 4. Ford, D.; Williams, P. Karst Hydrogeology and Geomorphology; John Wiley & Sons Ltd.: Chichester, UK, 2007. 5. Biondic, B., Bakalowicz, M., Eds.; Hydrogeological Aspects of Groundwater Protection in Karstic Areas; Cost Action 65; European Commission: Luxembourg, 1975. 6. Tulipano, L.; Fidelibus, M.D.; Panagopoulos, A. Groundwater Management of Coastal Karst Aquifers; COST Action 621, Report EUR 21366; European Commission: Brussels, Belgium, 2005. 7. Polemio, M. Degradation risk owing to contamination and overdraft for Apulian groundwater resources (southern Italy). In Proceedings of the Water Resources Management in a Vulnerable Environment for Sustainable Development, Perugia, Italy, 23 25 November 1998. 8. Taniguchi, M.; Dausman, A.; Howard, K.; Polemio, M.; Lakshmanan, E. Trends and Sustainability of Groundwater in Highly Stressed Aquifers; IAHS Publications: Wallingford, UK, 2009; Volume 329, p. 312. 9. Pisinaras, V.; Petalas, C.; Tsihrintzia, V.A.; Zagana, E. A groundwater flow model for water resources management in the Ismarida plain, North Greece. Environ. Model. Assess. 2007, 12, 75 89. [CrossRef] 10. Zuffianò, L.E.; Basso, A.; Casarano, D.; Dragone, V.; Limoni, P.P.; Romanazzi, A.; Santaloia, F.; Polemio, M. Coastal hydrogeological system of Mar Piccolo (Taranto, Italy). Environ. Sci. Pollut. Res. 2015. [CrossRef] [PubMed] 11. Cotecchia, V.; Casarano, D.; Polemio, M. Piovosità e siccità in Italia meridionale tra il 1821 ed il 2001. L Acqua 2003, 2, 99 106. 12. Polemio, M. Seawater intrusion and groundwater quality in the Southern Italy region of Apulia: A multi-methodological approach to the protection. In Progress in Surface and Subsurface Water Studies at the Plot and Small Basin Scale; Maraga, F., Arattano, M., Eds.; UNESCO-IHP: Paris, France, 2005; Volume 77, pp. 171 178. 13. Giorgi, F.; Lionello, P. Climate change projections for the Mediterranean region. Glob. Planet. Change 2008, 63, 90 104. [CrossRef]

Water 2016, 8, 148 15 of 16 14. Candela, L.; von Igel, W.; Elorza, F.J.; Aronica, G. Impact assessment of combined climate and management scenarios on groundwater resources and associated wetland (Majorca, Spain). J. Hydrol. 2009, 376, 510 527. [CrossRef] 15. Cotecchia, V.; Polemio, M. L inquinamento e il sovrasfruttamento delle risorse idriche sotterranee pugliesi. In Proceedings of the VI Workshop del Progetto Strategico Clima, Ambiente e Territorio nel Mezzogiorno, Taormina, Italy, 13 15 December 1995. 16. Romanazzi, A.; Gentile, F.; Polemio, M. Modelling and management of a Mediterranean karstic coastal aquifer under the effects of seawater intrusion and climate change. Environ. Earth Sci. 2015, 74, 115 128. [CrossRef] 17. IPCC. Summary for Policymakers Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change; IPCC: Paris, France, 2007. 18. Cotecchia, V.; Grassi, D.; Polemio, M. Carbonate aquifers in Apulia and seawater intrusion. G. Geol. Appl. 2005, 1, 219 231. 19. Maggiore, M.; Pagliarulo, P. Circolazione idrica ed equilibri idrogeologici negli acquiferi della Puglia. Geol. Territ. 2004, 1, 13 35. 20. Polemio, M.; Limoni, P.P. L evoluzione dell inquinamento salino delle acque sotterranee della Murgia e del Salento. Mem. Della Soc. Geol. Ital. 2001, 56, 327 331. 21. Barrocu, G. Seawater intrusion in coastal aquifers of Italy. Coastal aquifers intrusion technology: Mediterranean Countries. Minist. Cienc. Tecnol. 2003, 2, 207 223. 22. Polemio, M.; Dragone, V.; Limoni, P.P. Salt contamination of Apulian aquifers: Spatial and time trend. In Proceedings of the 19th SWIM & 3rd SWICA Jointed Meeting, Cagliari, Italy, 23 27 June 2008; pp. 115 121. 23. Polemio, M.; Dragone, V.; Limoni, P.P. Monitoring and methods to analyse the groundwater quality degradation risk in coastal karst aquifers (Apulia, Southern Italy). Environ. Earth Sci. 2009, 58, 299 312. 24. Polemio, M.; Di Cagno, M.; Virga, R. Le acque sotterranee del Gargano: Risorse idriche integrative e di emergenza. Acque Sotter. 2000, 68, 41 58. 25. Doglioni, A.; Simeone, V. Data-Driven Modeling of the Dynamic Response of a Large Deep Karst Aquifer. Procedia Eng. 2014, 89, 1254 1259. [CrossRef] 26. Polemio, M.; Casarano, D. Rainfall and drought in southern Italy (1821 2001). In The Basis of Civilization Water Science?; Rodda, J.C., Ubertini, L., Eds.; IAHS: Rome, Italy, 2004; Volume 286, pp. 217 227. 27. Polemio, M.; Casarano, D.; Limoni, P.P. Apulian coastal aquifers and management criteria. In Proceedings of the 21 SWIM Meeting, Azores, Portugal, 16 20 May 2010; pp. 203 206. 28. Polemio, M.; Dragone, V.; Limoni, P.P. The piezometric stress in the coastal aquifers of a karstic region, Apulia, Italy. In Sustainability of Groundwater in Highly Stressed Aquifers; Taniguchi, M., Dausman, A., Howard, K., Polemio, M., Lakshmanan, E., Eds.; IAHS: Hyderabad, India, 2009; Volume 329, pp. 138 144. 29. Polemio, M.; Dragone, V.; Limoni, P.P. La disponibilità di acque sotterranee in Puglia negli ultimi 80 anni. In Le Modificazioni Climatiche e i Rischi Naturali: Bari; Polemio, M., Ed.; CNR IRPI: Bari, Italy, 2011; pp. 201 204. 30. Cotecchia, V.; Polemio, M. Apulian groundwater (Southern Italy) salt pollution monitoring network. Flem. J. Nat. Sci. 1999, 79, 197 204. 31. Cotecchia, V.; Limoni, P.P.; Polemio, M. Identification of typical chemical and physical conditions in Apulian groundwater (southern Italy) through well multi-parameter logs. In Proceedings of the XXXIX IAH Congress Hydrogeology and Land Use Management, Bratislava, Slovenská, 6 10 September 1999; Fendeková, M., Fendek, M., Eds.; pp. 353 358. 32. Apulia Region. Piano Regionale Risanamento delle Acque. Boll. Uff. Regione Puglia, Bari, Italy, 1984, 57. 33. Apulia Region. Piano di Tutela delle Acque (PTA); Regional Water Protection Plan (PTA); Water Protection Service: Bari, Italy, 2009. (In Italian) 34. Polemio, M.; Dragone, V.; Limoni, P.P. Effetti antropici e naturali sul degrado quantitativo delle acque sotterranee del Tavoliere. In Proceedings of the III Conv. Naz. Sulla Protezione e Gestione delle Acque Sotterranee per il III Millennio 1999, Parma, Italy, 13 15 October 1999. 35. McDonald, M.G.; Harbaugh, A.W. A Modular Three-Dimensional Finite-Difference Ground-Water Flow Model; Techniques of Water-Resources Investigations of the United States Geological Survey; US GPO: Washington, DC, USA, 1988; p. 6.

Water 2016, 8, 148 16 of 16 36. Guo, W.; Langevin, C.D. User s Guide to SEAWAT: A Computer Program for Simulation of Three-Dimensional Variable-Density Ground-Water Flow; Open-File Report 01-4340; U.S. Geological Survey: Miami, FL, USA, 2002. 37. Cotecchia, V.; Daurù, M.; Limoni, P.P.; Mitolo, D.; Polemio, M. La valutazione della vulnerabilità integrata degli acquiferi. La sperimentazione nell area campione di Corigliano in Salento. Acque Sotter. 2002, 77, 9 20. 2016 by the author; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC-BY) license (http://creativecommons.org/licenses/by/4.0/).