A First Step towards Developing a Decision Support System Based on the Integration of Environmental Monitoring Activities for Regional Water Resource Protection

: The topic of diffuse pollution is of particular interest from technical, scientiﬁc, and administrative management points of view. Diffuse pollution is deﬁned as the contamination or chemical, physical, or biological alterations of environmental matrices caused by diffuse sources and not attributable to a single origin. In this study, various sources of diffuse pollution such as nitrates, pesticides, metals, and plastics were analysed. This was aimed at the implementation of a decision support system able to represent the state of environmental matrices degradation, with particular attention to water resources, and to make decisions evaluating similar environmental contexts. The potential of the developed system makes it possible to identify areas with the same environmental characteristics, referring to the various activities that create diffuse pollution and areas with the same pressure values on the environmental matrices. The system provides the political decision-maker with greater awareness of the environmental state, thus enabling him to apply more accurate land management policies. The created system, based on open-source software, which can be implemented with additional available data sources, is characterised by a data processing workﬂow that provides output information at the municipal level, so that it can be managed both by mayors and regional managers who are able to share the same information with all.


Introduction
The Puglia region (south-eastern Italia peninsula) is ruled by a flat and hilly territory, with a coastline 800 kilometres long between the Adriatic and Ionian Seas. The territory is mainly karstic [1]. It is characterised by a wealth of water resources, mainly underground [2]. The water resources to be safeguarded in Puglia include mainly aquifers, rivers, many lakes, and artificial basins. As for surface water resources, among the main streams are the Ofanto, one of the longest rivers in southern Italy, the Cervaro, the Fortore, the Bradano, and the Basento. Among the lakes stands out Lesina Lake, located in the coastal plain of Tavoliere, characterised by an important bird and fish fauna. The territory of the Puglia region and the environmental matrices are constantly under pressure due to numerous factors. The results of regional monitoring of underground rivers and aquifers indicate the presence of several threats [3,4].
There is a need to implement specific actions to combat the depletion of water resources. Biological, physical, and chemical factors can compromise a water body compared to its beneficial use, generating a condition of contamination or degradation. The extent of pollution depends on the pollutant, the type of water body, its location, and the benefits it brings. Natural events, able to compromise the balance of water ecosystems, can occur but in most cases, the pollution arises from anthropogenic activities.
Quarries can potentially become lakes when certain conditions are met, such as precipitation, interaction with groundwater, or with human intervention. Quarries under the water table imply an irreversible change in land use, as they form artificial lakes. Quarry lakes are particularly numerous in the floodplains, as in northern Italy in the Po Valley. Such systems create an area of high vulnerability to contaminants that could reach the aquifer [29]. Quarries flooded by groundwater can be used as hydropower storage systems. Large masses of water have an impact on adjacent aquifers, especially for high hydraulic conductivity and porosity rock values [30].
Illegal landfills have a significant impact on the environment and public health, especially if there is asbestos [31]. Waste can contain hazardous chemicals and produce toxic gases that can cause respiratory problems and other long-term diseases. Furthermore, illegal landfills can contaminate groundwater and compromise air quality. To address the problem of illegal landfills, we need vigilant policies on the part of the authorities and the introduction of tougher sanctions for those who violate the legislation.
Furthermore, the impact of urbanisation on the quality of surface water compromises its use for irrigation purposes to the point that it can no longer be used for this purpose. Many indices such as Water Quality Index (WQI) and the Nemerow Pollution Index (NPI) [32,33] can be used to monitor this aspect.
This work aims to treat all these aspects together to enable the decision-maker to make the best possible choices, based on the integration of data from various types of monitoring at the regional level, not only to avoid the possible pollution of water resources but also to improve the state of chemical and biological quality [16,34]. To achieve this goal, a data processing workflow has been developed to support decisions and increase the effectiveness of the analysis. This provides support for all those who need to make strategic decisions for issues that cannot be solved easily. This software falls into the category of decision support system (DSS). The terms indicate attention to decision-making processes based on the support of computer technologies through human-machine interaction [35]. The main function of a DSS is to extract in a short and versatile way useful information for the decision-making processes, coming from a significant amount of data. The DSS created does not require long series of data inputs, powerful computers, and advanced know-how to make it work as could happen with other very sophisticated tools [36].

Study Area
The Puglia region is located in southeast Italy with a geographical extension of about 19,345 km 2 . The region is characterised mainly by coastal plains, hills, and plateaus. The mountains are not very high; the highest point of the region is Mount Cornacchia at 1168 m. The geomorphology of the region is extremely varied. The western part of Puglia is characterised by flat and fertile coastal plains, while the eastern part is dominated by rocky mountains and hills. The area of Alta Murgia is the most important mountain range, which extends near Bari ( Figure 1).
Agriculture is one of the most important sectors of the economy of Puglia. The region is particularly known for its olive trees and olive oil production [37]. Wheat, corn, potato, fruit, and vegetables are important products in Apulian agriculture. The crops are widespread everywhere but mainly in the Tavoliere area to the north. Products of animal origin are another important economic resource in Puglia. Cheese production is widespread in the region. Beef, pork, and sheep are also important in the zootechnics of Puglia, mainly in Murgia dei Trulli. Very important are also the cultivation of mussels in the areas of the Mar Piccolo and Laguna di Lesina [37].
There are many natural parks of environmental importance, including the Gargano National Park, the Reserve of Bosco Incoronata, the Reserve Saline di Margherita di Savoia, and the Regional Nature Reserve Bosco delle Pianelle, just to name a few. These nature parks are characterised by a wide range of natural habitats [38]. Thanks to its cultural and natural heritage, the region of Puglia is a very popular destination for tourists. National Park, the Reserve of Bosco Incoronata, the Reserve Saline di Margherita di Savoia, and the Regional Nature Reserve Bosco delle Pianelle, just to name a few. These nature parks are characterised by a wide range of natural habitats [38]. Thanks to its cultural and natural heritage, the region of Puglia is a very popular destination for tourists. Finally, Puglia includes over 200 municipalities, among these the most populous being Bari, Brindisi, Taranto, Lecce, and Foggia.

Pesticide Analysis Methods
Several methods for the analysis of organic micropollutants, including pesticides, have been optimised over the last decade. Considering the wide variety of plant protection products on the market, refinement of traditional analysis techniques, including gascromatography (GC) and high-performance liquid chromatography (HPLC) with different types of detectors, has been observed over the years [39].
European regulations [40] have imposed very restrictive legal limits in the determination of active phytosanitary substances in water, to protect living beings and the environment. For this, it is often necessary to precede the instrumental analysis with techniques of preparation and concentration of the sample (liquid-liquid extraction, solidphase extraction, solid-phase-microextraction) [9].

Nitrate Analysis Methods
Nitrate ions can be quantified using laboratory techniques, e.g., ion chromatography detection with the Metrohm 930 compact flex IC [41]. Nitrate detection can, also, be performed with in situ sensor-based methods. In the first case, it is an indirect determination method which allows high sensitivity and a good detection limit (LOD). However, it has Finally, Puglia includes over 200 municipalities, among these the most populous being Bari, Brindisi, Taranto, Lecce, and Foggia.

Pesticide Analysis Methods
Several methods for the analysis of organic micropollutants, including pesticides, have been optimised over the last decade. Considering the wide variety of plant protection products on the market, refinement of traditional analysis techniques, including gas-cromatography (GC) and high-performance liquid chromatography (HPLC) with different types of detectors, has been observed over the years [39].
European regulations [40] have imposed very restrictive legal limits in the determination of active phytosanitary substances in water, to protect living beings and the environment. For this, it is often necessary to precede the instrumental analysis with techniques of preparation and concentration of the sample (liquid-liquid extraction, solid-phase extraction, solid-phase-microextraction) [9].

Nitrate Analysis Methods
Nitrate ions can be quantified using laboratory techniques, e.g., ion chromatography detection with the Metrohm 930 compact flex IC [41]. Nitrate detection can, also, be performed with in situ sensor-based methods. In the first case, it is an indirect determination method which allows high sensitivity and a good detection limit (LOD). However, it has the disadvantage of expensive instrumentation. In the second case, we refer to direct methods, such as those involving the use of biosensors. The latter is a viable alternative for in situ monitoring of nitrate ions with good data repeatability [42].

Quarry Mapping
A standard methodology was used to study mining activities in the Puglia region. The software used was Quantum GIS, a free and open-source geographic information system application [43]. The reference system used was the World Geodetic System 1984, EPSG:4326, of the global geodetic coordinate system. The base maps used for the identification of the sites were the orthophotos of the Puglia region. The time range used was from 2006 to 2019 (years 2006, 2010, 2013, 2016, 2019). For the analysis of the extractive zones of the territory, the database of the authorised quarries of the Puglia region was fundamental (provided by administrative offices). To facilitate the mapping of the regional area, a 3 km × 5 km grid was created. In this way, it was possible to analyse individually, at a more detailed scale, portions of territory for more representative results. A remapping of the previously mapped quarries was carried out, based on the most recent orthophotos, with vector typing tools in QGIS. New areas with clear evidence of mining activity have been mapped. The project has been implemented with various themes including the municipalities of Puglia, the official list of natural protected areas (EUAP), land use, and urbanised and other themes. This framing allowed us to have a full view of the area to be monitored and to obtain the perimeters of unauthorised quarries through the application of change detection techniques [44]. Through the activities of geoprocessing and overlay mapping, we have identified areas where the presence of excavation activities outside the known authorisation limits is evident. The areas thus identified have been saved in a geodatabase. It is important to note that the use of this methodology can provide a rough estimate of the extractive activities. Further investigations and verifications are needed for a complete and accurate assessment, as a field finding of remote data collected.

Further Data Sources and Metrics Implementation for the Evaluation of the Most Compromised Areas
We used additional information from regional databases by aggregating and normalising municipal data. Below, we explain the type and reasons for using them: • cultivated hectares on municipal extension as a source of pesticides and plastics. The data come from the Corine Land Cover Project [45]; code 2 was used, i.e., agricultural areas. There are currently no perimeters of areas vulnerable to pesticides, so agricultural areas are a useful alternative source of information; • perimeter delimitation of Nitrate Vulnerable Zones (NVZ) as a source of nitrate pollution of established agricultural origin [12]; • number of sites with waste and illegal dumps discovered and their extension as a source of pollution by metals, other pollutants, and plastics (confidential dataset); • number of quarries and quarry extensions as a source of pollutants in the subsoil [4]; • the perimeter of protected areas, values of the intrinsic vulnerability of aquifers, and results of analytical monitoring of pesticides in Apulian water bodies as elements for discussion [9].
These data make it possible to carry out a qualitative assessment of the pressures on the water resource, also providing the cognitive elements for the identification of the protection and recovery measures aimed at guaranteeing the achievement or maintenance of the environmental objectives envisaged in the Water Framework Directive (WFD) [46]. Pressure assessment can be used to evaluate whether a specific element that varies causes an increase in specific pressure and therefore whether it is possible to intervene by defining guide criteria and identifying any prescriptions and measures to be implemented to reduce this pressure.

Cluster Analysis
To be able to compare the data from the different databases, we proceeded to normalise the values from 0 to 1 by dividing the value of each by the maximum value. In the case of sites subject to waste dumping and illegal dumping and quarrying, we normalised the coefficients of the average area extents by dividing the total extents of the sites by the number of them. With the normalised data, a cluster analysis was performed to identify subgroups of observations within a dataset. The aim was to identify as many different groups as possible that correspond to different observations. The method applied was k-means clustering, an unsupervised, machine-learning method for classifying subgroups. This method classifies records in such a way that objects within the same cluster are as similar as possible (high intra-class similarity), whereas objects in different clusters are as different as possible (low inter-class similarity). In k-means clustering, each cluster is represented by its centroid, which corresponds to the average of the points assigned to the cluster. Therefore, the total variation within clusters must be as low as possible. The algorithm for detection is from Hartigan-Wong (1979) [47], which defines the total variation within the cluster as the sum of the squared Euclidean distances between the elements and the corresponding centroid, according to the following: with: • V as total intra-cluster variation; • x i as single record belonging to the k cluster; • mk as the average value of the points assigned to the k cluster.
Each observation is therefore assigned to cluster k; thus, the sum of the squares of the distances of the observations from the centroids is minimal. Before running the algorithm, the number of clusters that will be generated must be identified and specified. This decision is usually made by the data analyst. To determine the optimal number of clusters for our purposes, we used the Elbow method [48]. Applying this method, we calculated the clustering algorithm for different values of k; for example, varying k from 1 to 10 clusters. For each value of k, the sum of the squares of the intra-cluster distances is calculated. The trend curve is drawn and the points where it breaks are observed.

Results
In the Puglia region, the mining activity represents about 10% compared to the national one; this depends mainly on the presence in the territory of economically important carbonate rocks for industry and construction.
There are three main extraction basins in the area: Apricena-Poggio Imperiale, in the northwestern sector of the region characterised by Mesozoic carbonate rocks; Trani-Bisceglie-Ruvo-Fasano, in the central sector with Cretaceous limestones; and Salento with the famous "Pietra Leccese", in the south-western sector characterised by marly calcarenite. Extraction in the regional territory, in authorised sites, and production is mainly represented by limestone, clay, sand, and gravel [49]. The annual quantities (m 3 ) by type of material extracted are 140,000 sand and gravel, 375,000 ornamental stones, 4,280,000 limestone, and 39,000 clay.
The mapping activity led to the identification of 7700 geometric polygons which identify quarries of various sizes in the regional territory, of which 3900 are or have been authorised while 271 have never been authorised. These sites are characterised by the fact that the concession limits have been exceeded with new excavation fronts or new quarries have been built in areas where it is not possible to do it. Therefore, this map was carried out by identifying those sites outside the authorised perimeters and which, sometimes, may also have disappeared, over time, because they were filled in (the nature of the materials is unknown to us). An example of a record identified through change detection is the one shown in Figure 2, in which the start of a quarry activity outside the authorised limits, in the period from 2006 to 2019 and the subsequent filling, is evident over the years.
Hydrology 2023, 10, x FOR PEER REVIEW 7 of 18 materials is unknown to us). An example of a record identified through change detection is the one shown in Figure 2, in which the start of a quarry activity outside the authorised limits, in the period from 2006 to 2019 and the subsequent filling, is evident over the years. Another sad phenomenon concerns the sites affected by uncontrolled waste abandonment. In these sites of various extensions, from punctiform sites to those of a few hectares, continuous abandonments of waste of various kinds occur over time because it is difficult for the authorities to keep this phenomenon under control with the men and means at their disposal. The ascertained data of these illegal dumps indicate the presence of about 2000 sites throughout the territory (some examples are shown in Supplementary Materials File S1). Figure 3 shows the aggregate data of the number of sites currently active and already reclaimed by the province to represent the distribution of the phenomenon in the regional territory. In the whole regional territory, the cultivated hectares are 1,419,638, while the hectares identified as NVZ are 452,905. These areas are more concentrated in the Tavoliere delle Another sad phenomenon concerns the sites affected by uncontrolled waste abandonment. In these sites of various extensions, from punctiform sites to those of a few hectares, continuous abandonments of waste of various kinds occur over time because it is difficult for the authorities to keep this phenomenon under control with the men and means at their disposal. The ascertained data of these illegal dumps indicate the presence of about 2000 sites throughout the territory (some examples are shown in Supplementary Materials File S1). Figure 3 shows the aggregate data of the number of sites currently active and already reclaimed by the province to represent the distribution of the phenomenon in the regional territory.
materials is unknown to us). An example of a record identified through change detection is the one shown in Figure 2, in which the start of a quarry activity outside the authorised limits, in the period from 2006 to 2019 and the subsequent filling, is evident over the years. Another sad phenomenon concerns the sites affected by uncontrolled waste abandonment. In these sites of various extensions, from punctiform sites to those of a few hectares, continuous abandonments of waste of various kinds occur over time because it is difficult for the authorities to keep this phenomenon under control with the men and means at their disposal. The ascertained data of these illegal dumps indicate the presence of about 2000 sites throughout the territory (some examples are shown in Supplementary Materials File S1). Figure 3 shows the aggregate data of the number of sites currently active and already reclaimed by the province to represent the distribution of the phenomenon in the regional territory. In the whole regional territory, the cultivated hectares are 1,419,638, while the hectares identified as NVZ are 452,905. These areas are more concentrated in the Tavoliere delle In the whole regional territory, the cultivated hectares are 1,419,638, while the hectares identified as NVZ are 452,905. These areas are more concentrated in the Tavoliere delle Puglie area in the north of the region. Figure 4 shows a map with the cultivated areas of the Puglia region and with the areas identified as NVZ. Puglie area in the north of the region. Figure 4 shows a map with the cultivated areas of the Puglia region and with the areas identified as NVZ. The data also made it possible to obtain a qualitative assessment of the pressures on water bodies. The pressures considered are those envisaged by the WFD [50]. This analysis allows us to identify a certain level of pressure on the water resources, which we have classified as low, medium, and high ( Figure 5).
This aspect is important as the monitoring programs are closely related to the pressures on water bodies for the selection of the quality elements to be monitored, i.e., the chemical, biological, and hydro-morphological ones, and therefore act as a useful decision-making tool to support the DSS. The data also made it possible to obtain a qualitative assessment of the pressures on water bodies. The pressures considered are those envisaged by the WFD [50]. This analysis allows us to identify a certain level of pressure on the water resources, which we have classified as low, medium, and high ( Figure 5). Regarding the results of the cluster analysis, we first report those relating to the identification of the number of clusters ( Figure 6). Looking at the curve obtained, the number  This aspect is important as the monitoring programs are closely related to the pressures on water bodies for the selection of the quality elements to be monitored, i.e., the chemical, biological, and hydro-morphological ones, and therefore act as a useful decision-making tool to support the DSS.
Regarding the results of the cluster analysis, we first report those relating to the identification of the number of clusters ( Figure 6). Looking at the curve obtained, the number of clusters to choose would range from two to five. Considering the territorial variability and the different databases used, by applying the value of two clusters, we would not give good emphasis to some aspects and processes taking place in the territory because we enclose too many variables in just two dimensions, while giving a value of five implies the creation of a number of clusters equal to the number of variables, and therefore, it would not make much sense. The value of four seems to us a good compromise to represent the different variables involved. Regarding the results of the cluster analysis, we first report those relating to the identification of the number of clusters ( Figure 6). Looking at the curve obtained, the number of clusters to choose would range from two to five. Considering the territorial variability and the different databases used, by applying the value of two clusters, we would not give good emphasis to some aspects and processes taking place in the territory because we enclose too many variables in just two dimensions, while giving a value of five implies the creation of a number of clusters equal to the number of variables, and therefore, it would not make much sense. The value of four seems to us a good compromise to represent the different variables involved. The results of the cluster analysis are shown in Figure 7. On the one hand, there is a certain overlapping of the clusters which represents an environmental continuity due to the situations present in different municipalities (this is normal), and on the other hand, each cluster is characterised by territorial peculiarities which allow us to easily identify the environmental contexts influencing the quality status of the environmental matrices, in particular the water resource. As previously explained, the clusters show a certain degree of overlap between them, which is expected since they are phenomena that do not The results of the cluster analysis are shown in Figure 7. On the one hand, there is a certain overlapping of the clusters which represents an environmental continuity due to the situations present in different municipalities (this is normal), and on the other hand, each cluster is characterised by territorial peculiarities which allow us to easily identify the environmental contexts influencing the quality status of the environmental matrices, in particular the water resource. As previously explained, the clusters show a certain degree of overlap between them, which is expected since they are phenomena that do not have a clear solution of continuity. The clusters also show the presence of some points at the borders and are very detached from others and the centroid.  Subsequently, the number of the cluster was associated with the relative municipality and this allowed the representation of these values through the GIS system with a thematic map (Figure 8). Subsequently, the number of the cluster was associated with the relative municipality and this allowed the representation of these values through the GIS system with a thematic map (Figure 8). Subsequently, the number of the cluster was associated with the relative municipality and this allowed the representation of these values through the GIS system with a thematic map (Figure 8). Finally, the values of the means of the normalised data relating to the clusters were compared to identify the driving forces that can explain the classification (Figure 9). Finally, the values of the means of the normalised data relating to the clusters were compared to identify the driving forces that can explain the classification (Figure 9). Non-overlapping values, identified by points, can be very spread out as in the case of clusters 1, 2, and 3 or can show some degree of linear continuity as in the case of cluster 4. Cluster 3 also shows a certain bidirectional variability, and some records completely cross the centre of the graph, positioning themselves both at the top and at the bottom, ranging from municipalities with intense tourist activity to others with a high presence of natural areas. Non-overlapping values, identified by points, can be very spread out as in the case of clusters 1, 2, and 3 or can show some degree of linear continuity as in the case of cluster 4. Cluster 3 also shows a certain bidirectional variability, and some records completely cross the centre of the graph, positioning themselves both at the top and at the bottom, ranging from municipalities with intense tourist activity to others with a high presence of natural areas.
Clusters 1 and 2 are almost superimposable, which does not happen for some records, i.e., for municipalities characterised by intense mining activity. The normalised data of variable A appear to be significant only for the subdivision of clusters 1, 2, and 3 from cluster 4 and this confirms the previous deduction. Cluster 4 mainly includes municipalities with less agricultural activity and a greater zootechnical and tourist vocation. This cluster includes the municipalities belonging to various protected areas such as the Gargano National Park, the Valle d'Itria, and the sub-Appennino Dauno, known for their tourist vocation and production of local dairy and meat products.
The data produced by variable E identify cluster 1, while it is evident that variable B identifies cluster 2. Cluster 3 represents the most common state in which municipalities are, unfortunately, characterised by the highest value of parameter C (this is not a good sign for environmental protection). Only in cluster 4, these values are slightly lower, probably because there is greater attention to the territory to avoid a perception by tourists of abandoned and dirty territory.

Discussion
Aquatic ecosystems play a vital role in ecological balance and human well-being. However, pollution resulting from human activities is causing a significant depletion of the water resource and ecosystem services provided by these aquatic habitats [51][52][53][54][55]. Ecosystem services, i.e., the benefits humans derive from ecosystems, include, for example, the provision of drinking water, climate regulation, air purification, flood control, fishing, and tourism [56].
Pollution of water resources, vital for various ecosystem services, is mainly caused by human activities, such as industry, intensive agriculture, the use of fertilizers and pesticides, and incorrect or illegal waste disposal. These pollutants are released directly or indirectly into surface or groundwater, altering water quality and compromising aquatic ecosystems. Toxic chemicals, such as heavy metals and pesticides, can directly harm aquatic plants and animals, reducing their diversity and abundance [57][58][59]. This impact can interrupt food chains, influence ecological balances, and compromise the survival of species that play a key role in ecosystems.
Discharges of pollutants can lead to the destruction of aquatic vegetation, such as algae beds or macrophytes, which provide shelter, food, and breeding sites for many aquatic species; habitat loss can lead to declining populations of fish, waterfowl, and other organisms dependent on these habitats. Furthermore, the accumulation of harmful chemicals can make the water unsafe and unsuitable for human and animal use. Fisheries and aquaculture can be greatly affected by declining fish populations, while compromised water quality can make safe water supply and wastewater treatment difficult. Pollution can also negatively affect aquatic tourism, reducing the attractiveness of these destinations.
The problem of the diffusion of pollutants and above all their infiltration can be tackled and explored through data analysis with models that evaluate the information coming from the characteristics of the matrices and weather data to provide risk values for water resources.
Some hydrogeological models are useful to evaluate the state of water resources: hydrogeological models, water balance models, contaminant transport models, geomorphological models, and groundwater monitoring models.
A hydrogeological model is a mathematical and computational approach that integrates geological, hydrological, and hydrogeological data. This allows us to simulate and understand the flow of groundwater in a certain area. Water balance models [60] are used to quantify the balance between water inlets and outlets in a specific area. Contaminant transport models [61] can be used to study the dispersion and migration of contaminants into groundwater. Geomorphology models [62] can be used to understand the formation and evolution of karst systems. Groundwater monitoring models use real-time monitoring data to assess changes in groundwater quality and quantity. These models can be used to detect and monitor the effects of mining activities on water resources, including the influence of quarries on the surrounding aquifers. In addition, modelling approaches may vary depending on the specific geological and hydrogeological characteristics of the region of interest.
Hydrogeological models make it possible to predict the transport of pollutants in the aquifers over time and to predict possible contaminations in time and space, as well as to better identify critical situations through vulnerability assessment. In our case study, we compared the results of the cluster analysis with the intrinsic vulnerability values of the aquifers created according to the DRASTIC model [12,63]. (Figure 10). Some significant aspects emerge when comparing the results of our classification relating to the normalised qualitative assessment of pressures on water resources with the latest data available by the Regional Agency for Environmental Protection (ARPA) Puglia for surface water bodies [64] and groundwater [65].
For surface water bodies, the picture is very critical in the areas of Tavoliere and Ofanto, as well as some areas west of the Tavoliere Salentino all identified in cluster 2, while other areas that today are subject to numerous pressure factors were in 2018 in a stationary state from an ecological point of view but worsening from a chemical point of view, such as the areas to the east and close to Mar Piccolo in the Arco Jonico Tarantino. From an ecological status point of view, the water bodies which are improving are almost all coastal marine (pp. 55-56, [66]).
For groundwater bodies, the "poor" chemical status of the water resource in the areas of the Tavoliere and Ofanto is confirmed, as well as in the areas of Arco Jonico Tarantino west of Mar Piccolo of Taranto. Patchy "poor" chemical status is also established in areas of the Tavoliere Salentino and in almost all the municipalities belonging to cluster 1. Water bodies that are in a "good" qualitative status are those of the Alta Murgia and in the southern area of the Tavoliere Salentino and central area of the Salento delle Serre (p. 48, [67]).
It should be noted that the recently published data refer to the 2016-2018 control pe- Some significant aspects emerge when comparing the results of our classification relating to the normalised qualitative assessment of pressures on water resources with the latest data available by the Regional Agency for Environmental Protection (ARPA) Puglia for surface water bodies [64] and groundwater [65].
For surface water bodies, the picture is very critical in the areas of Tavoliere and Ofanto, as well as some areas west of the Tavoliere Salentino all identified in cluster 2, while other areas that today are subject to numerous pressure factors were in 2018 in a stationary state from an ecological point of view but worsening from a chemical point of view, such as the areas to the east and close to Mar Piccolo in the Arco Jonico Tarantino. From an ecological status point of view, the water bodies which are improving are almost all coastal marine (pp. 55-56, [66]).
For groundwater bodies, the "poor" chemical status of the water resource in the areas of the Tavoliere and Ofanto is confirmed, as well as in the areas of Arco Jonico Tarantino west of Mar Piccolo of Taranto. Patchy "poor" chemical status is also established in areas of the Tavoliere Salentino and in almost all the municipalities belonging to cluster 1. Water bodies that are in a "good" qualitative status are those of the Alta Murgia and in the southern area of the Tavoliere Salentino and central area of the Salento delle Serre (p. 48, [67]).
It should be noted that the recently published data refer to the 2016-2018 control period. The results reported by ARPA Puglia, if compared with those reported in the present methodology, indicate that the situation is worsening in the areas Tavoliere, Ofanto, Tavoliere Salentino, Arco Jonico Tarantino, and Puglia Centrale, while in other areas, there is a stationary state, except for coastal marine bodies in which the situation, from an ecological point of view, is improving.
Of course, this is symptomatic of the fact that some policies are succeeding while others are not. For example, the huge investments in the field of water treatment are bringing results while in the agricultural field, there is still a lot to do, as well as in the field of environmental controls on pollution from waste abandonment. In this regard, a vast regional action is about to be launched on the installation of hidden micro cameras in critical places whose images will be processed semi-automatically through machine learning.
Comparing the data of the results of the cluster analysis ( Figure 8) with the values of intrinsic vulnerability ( Figure 10) of the Puglia aquifers, it emerges that the areas of cluster 2 in the north of Puglia are located on less vulnerable aquifers, while those in the centre and south are on a more vulnerable aquifer: this implies that the risks of pollution and depletion of the groundwater deriving from agricultural activities are much greater in the latter areas. The municipalities belonging to cluster 1 seem to have the groundwater at greater risk especially if located in Puglia centrale, in the Alta Murgia, and partly in the Salento areas. These municipalities are precisely those characterised by a relatively greater number of quarries. In the areas present in cluster 4, there seem to be more critical issues for the municipalities of the Murgia dei trulli; in fact, compared to those of the Gargano, they are located on a more vulnerable aquifer with the risk of contamination by nitrates of zootechnical origin. As far as cluster 3 is concerned, the areas most at risk are those with higher vulnerability values: therefore, Puglia centrale with the southern Murgia and all of Salento.
To increase the levels of knowledge deriving from the information provided by these databases, the use of decision support systems induces considerable advantages: it is possible to have access to large amounts of data in real time, allowing decision-makers to make informed and timely decisions. For example, decision-makers could monitor the effects of the application of a monitoring program, assess whether there are environmental improvements, and better allocate resources for land management.
Precisely for this reason, we have preferred to carry out these elaborations based on the extensions of the municipal areas. Using the extensions of the river basins would create areas with an administrative management conflict, leading to a deadlock of improvement actions.
By automating data processing operations to obtain reports through programming scripts, we can leave much more time for decision-makers to focus on the most important decisions to be made.
The system also aims at minimizing the time required for drafting reports on the progress of the regional plans and programs of the individual different monitoring activities. In addition, decision support systems can use advanced algorithms and data analysis techniques to provide highly accurate and detailed information for site-specific purposes [68]. For example, it is possible to identify the factor limiting the growth in the well-being of the population of a given territory, the motivation why an area is not attractive from a tourist point of view based on elaborations on requests determined by particular needs.
On the other hand, decision support systems also have disadvantages. First of all, dependence on machines could lead decision-makers to neglect their intuition and personal experience. Furthermore, these systems are based on data and algorithms, and if there are errors in the latter, decisions could be compromised [69]. Sometimes decision support systems can require significant investments in terms of infrastructure, software, and qualified personnel. Furthermore, the process of implementing such systems can be very complex and require extensive staff training to be dedicated. Last but not least, we also mention the security risks: the collection and maintenance of sensitive data by decision support systems could expose it to security risks, such as data breaches.

Conclusions
Decision support systems can be very useful for administrators who wish to make informed and timely decisions on a regional scale.
In this case study, we reported how to implement a DSS based on the processing of environmental monitoring data at the municipal level. The created workflow is deliberately simple but this does not mean that it is unable to provide useful information on the status of some areas. The DSS can be applied on a regional scale to protect water resources, which are often subject to both punctual and widespread pollution events. The effort was to use available data coming from complex monitoring programs simply and immediately to offer valid information to decision-makers. The operational workflow represents a first step in monitoring the impacts on the environment and in particular on the water resource of some activities, both legal and illegal, on a regional scale.
Environmental pollution from pesticides, nitrates, heavy metals, and microplastics poses a threat to human health and the environment. To minimize these problems, there is a need for regulated management and the promotion of sustainable practices, such as the use of natural fertilisers, the reduction of pesticide use, and the reuse and recycling of waste. These practices could be adopted according to the specific territorial needs highlighted by this DSS.
As with any technology, however, there are also drawbacks to consider. To maximize the benefits of DSS, they should be implemented judiciously to ensure that decision-makers maintain a holistic view of the information and situation as they make important decisions.

Conflicts of Interest:
The authors declare no conflict of interest.