Integrative Groundwater Studies in a Small-Scale Urban Area: Case Study from the Municipality of Penafiel (NW Portugal)

: The occurrence of groundwater in urban hard rock areas is mostly controlled by secondary permeability caused by the fracturing degree and weathering grade. Vulnerability GIS-mapping, monitoring, and infiltration/recharge of water resources, and the delineation of the environmental protective background are key issues in evaluating, planning, managing, and decision-making for urban water systems. A small-scale urban area in Penafiel City was used as a case study. Historical and current scenarios were compared using hydrogeological GIS-based modelling. Water sources and potential contamination sources were mapped around the study area. Groundwater is mainly acidic, with a low mineralization and a Cl-Na to Cl-SO 4 -Na or Cl-Ca to Cl-SO 4 -Ca hydrogeochemical facies. Considering potential contamination sources, the moderate and the moderate-high classes dominate the area. The combined approach of the vulnerability indexes suggests that most of the area has a low to moderate vulnerability to contamination. The Infiltration Potential Index in urban areas (IPI-urban) indicates the prevalence of a low to moderate infiltration. Groundwater recharge evaluation indicates values < 2.3 mm/year and <60 mm/year for summer and winter scenarios, respectively. Finally, the investigation addresses insights for decision-makers for the appropriate planning and sustainable management of groundwater resources in small-scale urban areas.


Introduction
"Water is the blood of a city; of your good or poor quality depends on much of the health status of a population" (Bourbon e Noronha, 1885) Sustainable urban groundwater systems are increasingly considered to be significant in global development issues such as management, protection, supply, security, and services (e.g., References [1][2][3][4][5][6][7][8][9][10][11][12][13]). Urban population growth, in combination with better-quality living standards, lead to increased water use and demand. In addition, climate change and environmental subjects address several issues related to the sustainable availability of water resources.
An approach to monitoring the availability and sustainability of water resources need to focus on all features of ecosystems and necessarily requires a transdisciplinary methodology that encompasses socio-economic and cultural perspectives, technical-scientific solutions based on a design and nature, and an ethical approach and attention to societal dynamics. The Sustainable Development Goals (SDG) of the United Nation's 2030 Development Agenda highlight the significant role of sustainable urban groundwater management and interactions with sanitation issues [14][15][16]. Reference [17] addresses the need to use a methodological approach to analyze groundwater interlinkages with SDG targets, supported by the clear identification of evidence-based and logical interlinkages, including in urban areas.
Reference [18] points out that the concept of urban water security addresses the provision of all different water system services, considering the overall well-being in terms of social equity and environmental sustainability, and also includes hazards, risks, uncertainties, and resilience. In addition, Reference [13] defines urban water security as an adequate supply of clean freshwater to support humans and ecosystems and calls for an improved understanding and the better management of the quantity, quality, and societal nexus of water resources. Reference [19] also highlights a vital topic focused on access to good-quality water as a basic human right. Reference [20] emphasizes urban groundwater as a core resource to human health, and well-being, decisive for sustainable development, including the rural in areas [21].
Groundwater represents a major portion of the available freshwater resources and is widely exploited around the world [22,23]. Sustainable urban water resources must be managed without causing adverse environmental impacts and urban planning must be designed to be sustainable with nature. In addition, urban water resources are variable and even asymmetric on the regional scale, with natural limitations to infiltration and/or recharge, because water is a finite resource and is an impervious material to surface block infiltration.
Hydrological conceptualization is a decisive step in urban groundwater site model development. It provides a systematic, internally consistent overview of system boundaries, properties, and processes, bridging the gap between a reliable hydrogeological characterization and groundwater modelling (e.g., References [24][25][26]).
The useful concepts of urban groundwater systems are supported by earth-based models and numerical modelling (probabilistic, deterministic, and/or stochastic) to outline prediction scenarios [24,25] and are sometimes incorporated in a multi-model approach [27]. They should be built into a robust and calibrated framework on a back-analysis and must include the intrinsic geological ground variability and uncertainty of earth-based systems, as well as a multi-model and risk assessment in a multi-hazard environment approach (see References [25,26], and references therein). In addition, Geographical Information Systems (GIS) and related technologies provide an accurate tool to improve the data analysis of geospatial information about urban hydrological systems. GIS mapping and GIS-enabled site modelling of coupled urban groundwater data are based on a multi-analysis approach based on reliable and valid high-quality data sources (e.g., References [24][25][26][28][29][30][31][32]).
The main goal of this research is to outline an integrative urban water system assessment by using hydrogeology, hydrogeomorphology, hydraulic, and GIS-based mapping studies, including a historical urban water supply understanding of a small-scale urban area (Penafiel city, NW Portugal). Figure 1 shows a methodological flowchart that allows for a hydrogeological conceptual site model in a small-scale urban area (the site investigated here is ca. 5.1 km 2 ). The study corresponds to a multidisciplinary approach using GIS technology based on a multicriteria methodology. The research is grounded on basic data collection in order to create a geo-database to organize and analyze the spatial data. The cartographic data were collected for the characterization of the study area relative to its topography, land use, geology, morphotectonics, hydroclimatology, net recharge, and urban hydraulics and sanitation. The inventory techniques allowed us to obtain the hydrogeological features, the surface and shallow groundwater potential contamination activities, and hydrotoponymical survey. The different fieldwork inventories (e.g., hydrogeological features, potential contamination sources and hydrotoponymical features) were performed using a highaccuracy GPS device (Trimble GeoExplorer) and field inventory datasheets were produced. The potential contamination activities were classified in categories by their sources of origin (details in References [33][34][35]): urbanization, industry, agriculture, water mismanagement, and miscellaneous. Groundwater vulnerability assessment cannot be directly measured in the field. In order to obtain groundwater vulnerability maps, GIS technology is often used. It is a powerful tool to analyze, process, and combine spatial layers. Given the complexity of urban areas, the application of various indexes allowed for the comparing and improving of vulnerability analyses [36]. The aim of using these various indexes was to carry out an integrated evaluation between the different methods, allowing for an assessment of the vulnerability of urban groundwater. In fact, for the urban groundwater vulnerability, an integrated assessment was performed using several indexes, including DRASTIC-A standardized system for evaluating groundwater pollution potential [37], GODS-Groundwater Quality Protection [38], DRASTIC-fm-A standardized system for evaluating groundwater pollution potential for fracture media [39], SINTACS-Assessing and mapping groundwater vulnerability to contamination [40], and SI-Susceptibility Index [41]. The hydrogeological background was the basis of the vulnerability assessment. In the vulnerability maps, the international color code for DRASTIC was used. Finally, the urban groundwater vulnerability assessment was achieved using an integrated evaluation of all the methods.

Materials and Methods
The definition of the areas with the greatest infiltration potential was established on a quantitative basis and a multicriteria analysis was performed in a GIS environment. Evaluation of recharge and discharge used the Infiltration Potential Index in Urban Areas (IPI-Urban), for details see References [34,35]. The factors used in the IPI-Urban were based, revised, and updated from the key bibliographical sources (e.g., References [34,35,[42][43][44][45][46][47], and their relative weight/score for each factor was calculated using the analytical hierarchy process (AHP). The AHP technique allowed for the comparison of factors to each other, assigning a weighting to each factor to describe the importance of each of those factors that contributed to the overall goal (e.g., References [48][49][50][51]). For the hierarchy of factors, the knowledge gathered in the fieldwork was fundamental. The IPI-Urban was a weighted sum of eight factors: hydrogeology, tectonic lineaments density, land use, drainage density, slope, sewer network density, stormwater network density, and water supply network density. Each factor was spatially represented and analyzed in a GIS environment. The application of the IPI-Urban allowed us to obtain a map (regular grid of 4 × 4 m 2 ) representing the combination of all factors, ranging from 0 to 100. The areas where the IPI-Urban index showed higher values, corresponding to the more favorable conditions for water infiltration, according to all factors.
For the Penafiel urban area, two scenarios were developed, covering both summertime and wintertime. Table 1 presents the relative weight of the explained factors for each scenario. Since the study area is hot and dry during the summer, the stormwater network density was not considered for that scenario.
In order to calculate the urban recharge (mm/year), an analysis of the amount of rainfall was performed. That allowed us to conclude that the period between June and August was predominantly dry. According to regional hydrogeological key studies [47,52], performed in fissured media and surrounding the study area, an initial recharge rate of 8% was considered. The yearly average for dry months (35.4 mm/year) was used to calculate the urban recharge in the summer scenario, while the average for the remaining months (wet period) was used for the winter scenario (1299.3 mm/year).

General Setting and Hydroclimatology
Penafiel municipality occupies the interflow of the rivers Sousa-Tâmega-Douro in North Portugal. The site investigation is situated in the Northern district of Porto and it is called the capital of the Tâmega Sub-region. The studied area-the small-scale urban area of Penafiel-has 8761 inhabitants [53], with an area of 5.1 km 2 , which corresponds to a low population density (1717 people/km 2 ).
For this study, the urban area was delimited using the criteria of land use as well as geomorphological features. The study area is mainly located in a planned surface with an NE-SW orientation. Figure 2 presents the geographical and geomorphological setting of the Penafiel urban area. The region is mainly constituted by soils, essentially granitic, suitable for intensive agricultural purposes. The urban area is formed by granitic rocks and they outcrop in some parts of the urban area, particularly in the NE sector (Santa Marta) [54].
The urban area of Penafiel is characterized by a temperate climate with Atlantic influence, with the existence of a hot and dry period that corresponds to the summer months-more precisely, Csb in the Köppen climate classification [55,56]. Climate series from the Portuguese Water Resources Information System (SNIRH) were used to obtain a local characterization of the climate in the Penafiel region. The average annual rainfall is 1335 mm. The average monthly precipitation values of above 150 mm are found between November and March and the driest months correspond to July and August, with values rarely exceeding 30 mm monthly. The average annual temperature is 14.1 °C, with a minimum value in January (7.9 °C) and a maximum value in July (21.1 °C).

Urban Hydro-Historical Milestones
Landamnn [63] made impressive remarks in terms of the landscape, geology and water resources about Penafiel urban area and surroundings during the Napoleonic invasions of the Iberian Peninsula: "Hence we obtain a view of Penafiel on the summit of a third chain, two leagues distant, and from which we are separated by a beautiful valley richly cultivated, abundantly wooded, refreshed by numerous springs, the sources of small rivulets, and interspersed with substantial villas and cottages.
[…] The town of Penafiel is situated on a mountain, which extends from the river Souza to a small rivulet at two leagues beyond Penafiel: the town is not very considerable.
[…], the streets are very irregular, tolerably well paved, but narrow, and few of the houses exceed two stories in height.
Soon after our departure, in looking towards the north and west, we have a very grand view: Penafiel is on the left, whence towards the right extends the valley we traversed from Balthar, and beyond this horizon is bounded by the second great chain of mountains from Porto, where Balthar is seated near the summit; the fore-ground is composed of the declivities of the mountain on which we stand, forming extensive ravines, where scattered pine-trees, and large masses of soft granite, many of them nearly spherical, are seen on the surface in great profusion." References [64] and [65] clearly report the extremely bad sanitary conditions of the Penafiel urban area, circa 1928-1934. In fact, there was no urban sewage system, the urban stormwater arrangements were deficient and smelly water was often fed into the water supply system (details in References [64,65]). Only in the late 1930s and early 1940s were some key hydrological and hydraulic studies made about the Penafiel area [66,67].
The recorded urban water systems conditions are very similar to the concerns of other Portuguese municipalities in the 19th century (e.g., Lisbon, Coimbra, Figueira da Foz, Porto, Guimarães, and Braga) related to initiating comprehensive groundwater studies for urban water supply (e.g., References [57,60,62,[68][69][70][71][72][73][74][75][76], and references therein). In addition, References [77] and [78] reported several cholera epidemic outbreaks in Portugal, showing poorer hygienic conditions than most Portuguese cities. More recently, some urban groundwater studies and investigations were performed in Penafiel, aimed at the use of the urban water resources in the irrigation of green spaces and public fountains [26,79,80]. Figure 3 outlines the evolution of the key historical landmarks concerning the water supply and sanitation in the Penafiel urban area, as well as the overall historical framework in Portugal (e.g., References [63,81]).

Urban Hydraulics and Penafiel's Water Supply
Water supply systems should make use of the best available water sources, whether groundwater or surface water. Priority should be given to solutions based on existing water sources. By using the water from the existing underground mining galleries in the city of Penafiel, the efficient use of available water resources is promoted, contributing to sustainable management.
Water mine galleries capture fissured aquifers occurring in hard rocks, and are dug on hillslopes, where the topography intersects with the water table, almost always enabling gravity water supplies. This type of catchment can have significant flow fluctuations throughout the year, making them sometimes less reliable water sources.
The hydraulic study of Reference [67] refers to three water mines in the upper part of the city of Penafiel that fed public fountains, drinking fountains, and community wash-houses: (i) Parque Mine, (ii) Ponte Mine, and (iii) Santa Marta (or Chãos) Mine.
With population growth and increasing per capita consumption, the pressure on the water resources grew and there was a need to explore other sources. This led to the opening of galleries and water mines in Perafita (the Duas Igrejas site), whose water was sent to a reservoir in Penafiel (Alto do Sameiro), today wedged between the Penafidelis restaurant and the Pena hotel. Still, due to the lack of water in hot and dry years, in the 1970s water supply reinforcement works were carried out from the Sousa River to Penafiel, and an additional water reservoir was built in the Alto do Castro (Santa Marta). Later, a water catchment in the Tâmega River was constructed and this is now supplying the public water distribution network of Penafiel.
The old system was therefore abandoned and has deteriorated due to lack of use and maintenance. It should be noted that, in addition to the fact that there is no inventory of the old system, it was found, locally, that a large part of the network was destroyed when the public water supply system was built and also during urban regeneration works. To promote sustainable urban management, the feasibility of recovering the old water supply system was assessed, considering the available mine flows and the conservation status of existing pipelines.
It is noted that the main fountains and public green spaces of the city are located between the topographic elevations 300 and 350 m. Of the three mines, only the Santa Marta mine gallery can satisfy this range of elevations, as it is located at a higher topographic elevation, approximately 352 m. This provides the necessary pressure for the city's garden irrigation network. This management of water resources is also environmentally advantageous, as water from the mine does not carry harmful substances to green spaces.
An assessment of the flow rates for the current needs of the water supply network to public fountains and public green spaces reveals that the existing mine cannot meet all contemporary needs. Accordingly, a study was conducted on the use of the mine flow to partially meet these needs. A rational use of the flow rates was proposed, leading to the re-use of these waters and minimizing the costs associated with the public supply network consumption.
In the proposed supply system, the water is transported by gravity through buried pipelines. Initially, a connection between the Santa Marta Mine and a reservoir to be located at the top of Sameiro's Garden is proposed. The proposed solution replicates the old system that was deactivated a few years ago but still exists. It is a buried reservoir, built in stone masonry, from which the gardens and fountains were supplied.
The proposed network starts from a new reservoir to be located at the site of the old one, and its scheme is outlined in Figure 4. The design network begins by establishing the connection between the Santa Marta Mine (MNA) and the reservoir located at the top of Sameiro's Garden (RSV).
The need for water storage through the construction of a reservoir is due to hydraulic, economic, and safety reasons. This storage is ensured by the reservoir, which simultaneously performs the following functions: (i) to serve as a regularization device, compensating for consumption fluctuations; (ii) to set up water reserves to ensure distribution in the event of voluntary or accidental disruption of the upstream network; (iii) and to balance the pressures in the distribution network.
From the reservoir, three branches develop: the first pipeline conducts water to the gardens (Sameiro Garden-JD1, Praça da República Garden-JD2 and Calvário Garden-JD3); the second pipeline directs the water to a set of fountains, such as Largo da Ajuda (CH3); and the third pipeline leads the water to the fountain located on the steps of the Sanctuary of Nossa Senhora da Piedade (CH2) and Santos Passos (CH1). This kind of "cascade" layout solution seeks to achieve better water rationing.
The layout of the pipelines was defined, considering economic optimization criteria, passing on roads and paths, and assessing the costs of possible pavement replacements and interference with residents, natural obstacles, and other infrastructure, seeking to preserve the historical area.

Results and Discussion
A hydrogeological inventory was performed, and a total of 44 water sources were mapping (Figure 5a and 5d). Most of the water sources are fountains/springs and water tanks, and for nearly 50% of them, the origin of water is natural. The remaining water sources are supplied by the municipal network.
Two fieldwork campaigns were performed and in situ determinations included temperature (°C), pH, and electrical conductivity (µScm −1 ). Water samples were collected in each campaign for hydrogeochemical and radioisotope analyses. The hydrogeochemical analyses included major element concentrations, as well as nitrates. To evaluate the impact of radioisotopes on water quality, radon gas ( 222 Rn) was analyzed. The laboratory groundwater results were analyzed with AquaChem software (version 5.1., Waterloo Hydrogeologic, Waterloo, Canada).
The groundwater has an acidic character (median pH = 5.5) with relatively low mineralization (electrical conductivities < 300 µScm −1 ) and three water sources are nitrate enriched: ws1 (65 mgL −1 ), ws4 (94 mgL −1 ), and ws6 (68 mgL −1 ). The hydrogeochemical facies for most of the groundwater are Cl-Na to Cl-Ca and Cl-SO4-Na to Cl-SO4-Ca. Concerning radon gas ( 222 Rn), the concentrations are in the range 35-549 BqL −1 and the median value is 359 BqL −1 (Figure 5c). Considering the relatively low mineralization and the hydrogeochemical facies, the groundwater appears to have a shallow flow path, with short periods of time between recharge and discharge areas. Moreover, sulphates and nitrates are likely to have an anthropogenic origin, most probably related to the proximity of agricultural fields and public gardens, where fertilizers and pesticides are used. These results are in good agreement with other data reported for urban areas in granitic terrains (see e.g., References [34,52,[82][83][84].
Here, 222 Rn concentrations in groundwater are within the variation interval observed in Portugal for granitic aquifers [84,85].
An inventory of potential groundwater contaminant activities is based on thorough and systematic work. In urban areas, it is fundamental to understand whether there is contamination, and what its origin is (natural or anthropogenic). It is also an important contribution to a future analysis of the vulnerability of water resources to contamination. It should also be treated as an indispensable part for the development of a groundwater protection strategy [33].
In the studied area, 70 groundwater potential contamination sources were identified, and a field datasheet was created for each one. The groundwater potential contamination sources were inventoried and classified according to the degree of contamination (adapted from Reference [38]). The degree of contamination was divided into four classes, namely: low, moderate, moderate to high, and high. Figure 6 shows the location of the potential contamination sources, as well as the degree of contamination, relative to the hydrogeological units. From an analysis of the figure, it can be concluded that, in the urban area of Penafiel, the high contamination class represents 37.1%, while the moderate and moderate to high classes represent 55.7% of the potential contamination sources. The low class represents 7.1%. Due to the potential for contamination, the nature of contaminations sources may contribute to the locally increasing vulnerability to contamination. In order to carry out the urban vulnerability assessment in the Penafiel urban area, several indexes were applied, namely the GODS, DRASTIC, DRASTIC-fm, SINTACTS, and SI (Figure 7).
The GODS index (Figure 7a) reflects the hydrogeological units, which can be justified by the smaller number of parameters used in its calculation. Most of the area presents negligible vulnerability (57.5%), and the remaining area corresponds to the low vulnerability class ( Table 2). The class with the highest representation occurs in HGU 2 (monzonite and porphyritic granite, medium grained- Figure 8). Figure 7b represents the DRASTIC index. The very low vulnerability class comprises 52.7% of the area, while the remaining area presents a low vulnerability. The relation between the hydrogeological units and the vulnerability classes ( Figure 8) allows us to conclude that the very low vulnerability class mostly occurs in HGU 2 and the low vulnerability in HGU 4 (porphyritic granite, coarse grained). The DRASTIC-fm index, which is an adaptation of the DRASTIC method by the addition of the parameter fm (fracture media), was also applied According to DRASTIC-fm (Figure 7c), the urban area presents low vulnerability (69%) overall. The predominance of this class of vulnerability occurs in both hydrogeological units. The study area, according to the SINTACS index (Figure 7d), is classified predominantly as having high vulnerability (77.9%). Finally, the application of the SI index allows the classification of the area into four classes of vulnerability (Figure 7e). Three classes of vulnerability are clearly distinguished: moderate (50.8%), low-moderate (24.7%), and low (20.5%). These correspond to 95.9% of the entire study area. In the calculation of this index, the LU (Land Use) parameter is added, and it can be concluded that its relationship with the hydrogeological units is very low.
In the study area, according to the various indexes applied, the vulnerability classes have different levels. It can be concluded that the increase of parameters used, as well as their weighting, contributes to a more accurate characterization of the investigation site. Finally, the integrative approach of intrinsic vulnerability is essential for sustainable water resource management in urban areas.   The application of the IPI-Urban allows for the definition of areas where infiltration is most favorable, according to a set of factors [34,35,47]. As such, the factors considered in the IPI-Urban are represented in Figures 8 and 9. Regarding hydrogeology (Figure 9a), HGU 4 (porphyritic granite, coarse grained) and HGU 2 (monzonite and porphyritic granite, medium grained) stand out. These hydrogeological units make up 99.6% of the total urban area. As can be seen in Figure 9b, the low slope class (55%) stands out, followed by the very low slope class (26.6%). Regarding land use (Figure  9c), the "urban fabric" class stands out, corresponding to 58.9% of the study area. In the tectonic lineaments' density (Figure 9d), the most representative classes are low and very low, corresponding to 83.6% of the urban area. Regarding the drainage density (Figure 10a), it was found that the most representative class is the high class, representing 75.8% of the area under study. The remaining parameters relate to urban hydraulics and sanitation. With regard to the sewer network density (Figure 10b), the most representative classes are low and moderate, corresponding to 96.1% of the urban area. In turn, the stormwater network density (Figure 10c) is mostly low to moderate (96.1%). Regarding the last factor considered, the water supply network density (Figure 10d), two classes that represent 75% of this urban area can be highlighted, namely the moderate and the high-density classes. Figure 11 presents the IPI-Urban of the Penafiel area, and both the summer (Figure 11a) and winter scenarios (Figure 11b). The analysis of the figure, together with Table 3, allows us to conclude that, in this urban area and in both scenarios, the low to moderate infiltration class predominates. In the summer scenario, the low to moderate class corresponds to 71% of the area, with 82.9% for the winter scenario. The evaluation of the urban recharge was performed according to the methodology described above. Additionally, here, the two scenarios (summer and winter) are presented regarding recharge in the urban area of Penafiel.  In the summer scenario ( Figure 11c and Table 3), the recharge classes are relatively evenly balanced, with the most representative being the class < 1.3 mm/year (36.1%), followed by the 1.5-2.3 mm/year (35.6%) class. The class least represented is the class 1.3-1.5 mm/year (28.3%). Recharge increases substantially in the winter scenario (Figure 11d and Table 3). Thus, the class 50-60 mm/year stands out, representing 46.6% of the urban area. It is followed by the 40-50 mm/year class, representing 30.8% of the area. The recharge class 60-70 mm/year covers 20.8% of this urban area.    Figure 12 is the conceptual hydrogeological model of the urban area of Penafiel and its surroundings. The analysis of the figure shows that the recharge is done by meteoric waters, mainly in the flattened areas between 350 and 385 m (Santa Marta), and in areas occupied by forests. Precipitation values vary between 1300 and 1350 mm/year, with infiltration rates of around 3%, reaching values in the range of 40-60 mm/year. The water table follows the topography and is practically at the surface, between 5 and 10 m depth. The occurrence of springs is generally associated with tectonic nodes, with flow rates generally below 0.4 L/s. Underground flow occurs at shallow depths of up to 50-60 m, along the weathered and fractured zones of the granitic rocks; the deepest flow occurs along both the main tectonic axes with adequate connectivity and the veins, namely aplite and aplite-pegmatite.

Urban Hydrogeological Site Conceptual Model
In the weathered zone (ranging 10-20 m), unconfined aquifer systems are found, with dominant interstitial permeability, hydraulic conductivity of over 0.1 m/day, and an effective porosity of around 5%. In deep fractured areas, zones of discontinuous circulation will occur, with hydraulic conductivity of less than 0.01 m/day, which may be compared to semi-confined aquifers.
The inventoried groundwater is short cycle, poorly mineralized water, mainly installed in the most altered and fractured area of the granite massif. At greater depths, there will be waters with longer residence time in the reservoir and higher mineralization.

Conclusions
The study of the ancient water supply system of the Penafiel urban area (NW Portugal) was performed to seek out the feasibility to recover the old water system, mainly springs, public fountains and community wash places. That approach was based on a comprehensive urban groundwater characterization and evaluation, including the understanding of the historical milestones of the urban water supply in the municipality of Penafiel (NW Portugal). In addition, a detailed survey, inventory, and inspection of the still accessible springs, water mines, and water facilities conservation state was performed.
According to hydro-historical studies (e.g., References [66,67]) and the results of current investigations, a water yield reference value of 34.5 m³/day (≈ 0.40 L/s) is obtained. Groundwater flow is mainly related with the Santa Marta spring. In fact, this value will partly guarantee a sustainable water supply to public fountains and urban green spaces. The connection to the public water supply should remain, which can be used whenever the water mine proves to be insufficient in order to ensure the reliable supply of the quantity and quality of water required. In order to take advantage of the water supply availability for irrigation of the green spaces, an alternating mode was suggested, i.e., during the day the supply of the fountains and at late evening switching of the flow to the irrigation of green spaces. In addition, this urban groundwater resources management is also environmentally balanced because of the water mine quality. According to the average irrigation requirements obtained by hydroclimatological analysis, a daily average of 13.3 m³ is considered acceptable to meet the needs of the grass-covered areas of the main city gardens (e.g., Sameiro Garden, Praça da República Garden, and Calvário Garden).
The proposed design and restoration of the ancient water supply system oriented to the irrigation of the main gardens and for the supply of some fountains will promote environmental protection, a design with nature and sustainable water resource management for the Penafiel city. A comprehensive approach is focused on the adaptation to climate changes and is suitable for reasoning planning and management of groundwater resources in a small-scale urban area.
Author Contributions: L.F. and H.I.C. designed the research. L.F., H.M., T.A., J.F.T., M.J.A., A.J.S.C.P. and H.I.C. gave input on the topics of fields survey and inventory, urban hydrogeology and hydrogeomorphology, hydrohistorical milestones, and GIS mapping of the study site. A.J.S.C.P. contributed with radioisotope hydrology. T.A. and J.F.T. wrote the urban hydraulics section. All authors contributed to the data analysis, interpretation and discussed results. L.F. and H.I.C. wrote the manuscript with contributions of all authors. All authors have read and agreed to the published version of the manuscript.