Next Article in Journal
Quantitative Analysis of the Sub-Cloud Evaporation of Atmospheric Precipitation and Its Controlling Factors Calculated By D-Excess in an Inland River Basin of China
Next Article in Special Issue
Water Level Prediction of Emergency Groundwater Source and Its Impact on the Surrounding Environment in Nantong City, China
Previous Article in Journal
Big Data Analytics and Its Role to Support Groundwater Management in the Southern African Development Community
Previous Article in Special Issue
Wildfire Effects on Groundwater Quality from Springs Connected to Small Public Supply Systems in a Peri-Urban Forest Area (Braga Region, NW Portugal)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Urban Groundwater Processes and Anthropogenic Interactions (Porto Region, NW Portugal)

1
Laboratory of Cartography and Applied Geology (LABCARGA), Department of Geotechnical Engineering, School of Engineering (ISEP), Polytechnic of Porto, 4200-072 Porto, Portugal
2
GeoBioTec Centre (Georesources, Geotechnics, Geomaterials research group, University of Aveiro, 3810-193 Aveiro, Portugal
3
CITEUC, Department of Earth Sciences, Faculty of Sciences, University of Coimbra, 3030-790 Coimbra, Portugal
4
Centro de Recursos Naturais e Ambiente (CERENA), Instituto Superior Técnico, University of Lisbon, 1049-001 Lisbon, Portugal
5
Centro de Ciências e Tecnologias Nucleares (C2TN), Instituto Superior Técnico, University of Lisbon, 2695-066 Bobadela, Portugal
*
Author to whom correspondence should be addressed.
Submission received: 28 July 2020 / Revised: 30 September 2020 / Accepted: 3 October 2020 / Published: 9 October 2020

Abstract

:
Groundwater in fissured rocks is one of the most important reserves of available fresh water, and urbanization applies an extremely complex pressure which puts this natural resource at risk. Two-thirds of Portugal is composed of fissured aquifers. In this context, the Porto urban region is the second biggest metropolitan area in mainland Portugal. In this study, a multidisciplinary approach was developed, using hydrogeological GIS-based mapping and modeling, combining hydrogeochemical, isotopic, and hydrodynamical data. In addition, an urban infiltration potential index (IPI-Urban) was outlined with the combination of several thematic layers. Hydrogeochemical signatures are mainly Cl-Na to Cl-SO4-Na, being dependent on the geographic proximity of this region to the ocean, and on anthropogenic and agricultural contamination processes, namely fertilizers, sewage, as well as animal and human wastes. Isotopic signatures characterize a meteoric origin for groundwater, with shallow flow paths and short residence times. Pumping tests revealed a semi- to confined system, with low long-term well capacities (<1 L/s), low transmissivities (<4 m2/day), and low storage coefficients (<10−2). The IPI-Urban index showed a low groundwater infiltration potential, which was enhanced by urban hydraulic and sanitation features. This study assessed the major hydrogeological processes and their dynamics, therefore, contributing to a better knowledge of sustainable urban groundwater systems in fractured media.

1. Introduction

Urbanization is the foremost global phenomenon of our time, and groundwater from springs and wells has been a vital source of urban water supply since the first settlements [1]. Urban hydrogeology is a key scientific field, and an understanding of groundwater in urban areas has increased greatly in the last decades. Attention was given to the relationship between urban development and groundwater resources beginning in the 1950s and 1960s of the 20th century, when accelerated growth after the Second World War began to cause significant diversity of hydrological problems. Most of these problems were related to urban runoff and floods, therefore, urban hydrology established itself in an affirmative way [2,3]. Thus, urban groundwater rapidly emerged as a necessary and urgent branch of study. Particularly, after the 1990s of the 20th century, several scientific publications on issues addressing the subject of hydrogeology in urban areas are found, namely [2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22].
The impact of urban development on geological processes, in general, and groundwater systems, in particular, has been highlighted for a long time (e.g., [12,16,23,24,25,26,27,28,29,30,31,32,33,34,35,36]). It is remarkable that only one century ago, 20% of the world population lived in urban areas, whereas, presently, 54% of the world population lives in urban environments [37]. Anthropogenic activity is the major geomorphic agent affecting the Earth’s land surface, while urbanization and agriculture are, perhaps, the major processes currently affecting the land (e.g., [7,14,20,23,38,39,40,41] and references therein). The need for provision of safe water, sanitation, and drainage systems are key elements which are vital to the understanding and management of groundwater resources in urban environments. Groundwater problems associated with urbanization and urban sprawl include both quantity and quality issues, namely (e.g., [2,42]), subsidence, covering of shallow systems, increased recharge from leaky water ascribed to sewage systems and irrigation return flow, changes in subsurface secondary porosity, permeability from utility systems and other constructions, as well as the effects of imported water resources, saline intrusion, and contamination from point and nonpoint sources. In this context of global anthropogenic activity, groundwater is a key source of urban supply worldwide, and aquifer storage represents a key resource for achieving water-supply security under climate change that will contribute to aggravate these pressures, with decreasing rainfall, as well as longer and more frequent and extended drought periods (e.g., [1,43,44]). Additionally, the impact of the complex system of man-made infrastructures (e.g., sewer, storm sewers, pipes, trenches, tunnels, and other buried structures) and impervious surfaces and or cover areas must be highlighted as anthropogenic disturbances in groundwater media (e.g., [11,42,45,46,47]).
Integrated multidisciplinary approaches are often required to address the scientific issues involving groundwater resources. In this way, environmental isotopes such as 2H, 18O, 13C, and 3H, can be used as fingerprints to assess some of these groundwater problems and to provide information needed for the rational management of these water bodies (e.g., [48,49,50]). Furthermore, the delineation of groundwater potential zones (GWPZ’s) using remote sensing, geographic information systems (GIS), and geovisualization techniques are effective tools for groundwater exploration and groundwater recharge (e.g., [13,17,18,35,36,51,52]). This approach encompasses the integration of several factors, organized from different geodatabases and a multiparametric approach, which influence the occurrence of groundwater, namely lithology, tectonic lineaments, land use, slope, drainage, precipitation, and urban hydraulics. Nearly half of the world depends on groundwater as the main source of drinking water [53]. However, groundwater is an underexploited resource in many urban areas, due to its invisibility, inadequate management, scientific uncertainty, and political strategies that favor the use of surface waters (e.g., [54,55,56,57]).
The main objectives of this study in the Porto urban region (NW Portugal) were:
  • To develop an integrated geoenvironmental assessment of groundwater resources in urban environments using geotechnological capabilities, particularly GIS mapping and geovisualization techniques, and also extensive field and laboratory work;
  • To evaluate groundwater quality and groundwater flow paths, by combining hydrogeochemical and environmental isotopic data and to identify the leading processes responsible for groundwater disrepair;
  • To assess groundwater quantity, by means of pumping test data;
  • To delineate groundwater infiltration potential zones at a regional scale, using the innovative infiltration potential index for urban areas (IPI-Urban) that integrates several layers of information properly weighted and overlaid in a GIS platform; a tool which should improve our understanding of complex urban groundwater recharge processes in future investigations;
  • To refine the regional hydrogeological conceptual model, merging all the data, to improve the understanding of urban groundwater systems in the Porto urban region.

2. Study Area: Land Cover and Hydrogeological Background

The Porto urban region is a densely populated region which includes other cities such as Vila Nova de Gaia, Matosinhos, Maia, and Vila do Conde, in the Porto Metropolitan Area, with about 1.7 million inhabitants [58], (Figure 1). According to the land use map for mainland Portugal [59], most of the investigated region (93.3%) is occupied by the following three classes: urban and industrial areas (41.4%), agricultural areas (28.5%), and forest areas (23.4%). Urban areas are concentrated primarily in the municipality of Porto and in the surrounding municipalities. Agricultural spaces are located mainly in the northern part, especially in the municipality of Vila do Conde. Forest areas are mainly found in the NE, especially to the east of Vila do Conde and in the municipality of Trofa.
Considering agriculture and farming, the most important food supplies produced in this region are potatoes, vegetables, fruits, and maize, as well as cattle raising, namely bovines. According to [20], the two principal types of fertilization are organic which uses mostly manure, as well as sewage and inorganic. The inorganic fertilizers are mostly composed of nitrogen, phosphorous, potassium, and nitrogen which is the most prevalent. A sprinkler system is the leading irrigation system.
This region is drained by the following four major streams: the Ave River, in the north; rivers Onda and Leça in the centre; and Douro River, in the south. Located close to the Atlantic Ocean, this region has a temperate climate with dry and mild summers (Köppen climate classification Csb), with a mean annual precipitation around 1200 mm and a mean annual air temperature of 14 °C. There is a water deficit from June to September, mainly in July and August [60,61,62]. The availability of groundwater resources is reflected by these climatic conditions, along with geologic and morphotectonic features. The regional geological framework comprises a crystalline fissured basement of highly deformed and overthrust Late Proterozoic/Palaeozoic metasedimentary and granitic rocks. The post-Miocene sedimentary deposits frequently cover the bedrock [63,64,65] and references therein. The morphotectonic background includes a littoral platform characterized by a regular planation surface dipping gently to the west, ending around 100 m a.s.l [66]. The Ave, Leça, and Douro rivers have incised valleys that cross the flatness of this morphological surface, suggesting a regional tectonic control. The main regional hydrogeologic units [12] in this area are (Figure 2) the following: (i) sedimentary cover, particularly beach and dune sands, alluvia, sandy silts, and clays; (ii) metasedimentary rocks, namely micaschists, metagraywackes, and paragneisses; and (iii) granitic rocks, mainly two-mica granite, medium to coarse grained and biotitic granite, medium to fine grained.

3. Materials and Methods

The study area was assessed using various tools, as proposed by [67,68], namely geological and hydrogeological mapping, hydrogeochemical and isotopic techniques, and hydrodynamic characterization. Fieldwork campaigns were first carried out to identify major geologic features accountable for groundwater circulation pathways and to assess litho-structural heterogeneities. The field and laboratory data were evaluated in a GIS environment (ArcGis|ESRI, Redlands, CA, USA, and OCAD for Cartography Inc., Baar, Switzerland).
An urban hydrogeological field inventory was developed, which included more than 300 water points (springs, water mines, fountains, boreholes, and dug wells). Three fieldwork campaigns were carried out in May 2008, November 2008–January 2009, and May–June 2009. A total of 40 water points were selected for groundwater sampling and for chemical analyses consisting of 34 boreholes (mean depth of 117 m), 3 dug wells (mean depth of 11 m), 2 springs, and 1 fountain (connected to a water mine) (see Figure 2), most of them used for industrial and agricultural water supply. The inclusion of 2 dug wells and 2 springs in the fieldwork campaigns basically developed in the period 2005–2007 and were justified as follows: (i) to have dug wells represented in the hydrogeological unit “micaschists, metagraywackes and paragneisses”; (ii) to have springs represented (the 2 springs selected were two important springs of the historic Paranhos and Salgueiros water galleries in Porto City, which date back to the 16th Century).
Among the 40 samples, isotopic determinations of δ18O, δ2H, and 3H were performed in 14 boreholes, 2 dugwells, and 2 springs in one fieldwork campaign. All the inventory and sampling sites were georeferenced with a high-accuracy GPS (Trimble® GeoExplorer, Sunnyvale, CA, USA). In situ measurements included water temperature (°C), pH, electrical conductivity (µ Scm−1), dissolved oxygen (mg L−1), and redox potential (mV) using a multiparametric portable equipment (Hanna Instruments, HI9828). The hydrogeochemical analyses included major and minor element concentrations and were acquired at the “Instituto Nacional de Saúde Doutor Ricardo Jorge” (Porto, Portugal). Environmental isotopes (18O, 2H, and 3H) were measured at the ‘‘Instituto Tecnológico e Nuclear’’, meanwhile renamed “Campus Tecnológico e Nuclear/Instituto Superior Técnico (CTN/IST)” (Bobadela, Portugal). The δ18O and δ2H results were reported to the Vienna-Standard Mean Ocean Water (V-SMOW), and determinations were carried out using a mass spectrometer SIRA 10 VG-ISOGAS, applying the analytical methods described in [73,74], respectively. The tritium content was determined using the electrolytic enrichment and liquid scintillation counting method described by [75,76], using a PACKARD TRI-CARB 2000 CA/LL (see [77,78]). AquaChem 5.1 software was used for the hydrochemical interpretation.
The hydrodynamic interpretation of the groundwater systems was carried out based on pumping tests which were gathered from hydrogeological and geotechnical scientific and technical reports. Pumping data were analyzed using the Theis method, the Cooper–Jacob method, and the Theis recovery method [79,80,81,82], and the transmissivity and storage coefficient were evaluated. Moreover, the Logan approximation [83] was used to estimate transmissivity. The values of transmissivity, storage coefficient, and long-term well capacity were compiled from the hydrogeological reports and also contributed to the discussion.
The delineation of groundwater infiltration potential zones was assessed by the IPI-Urban index (details in [13,17,18,19]). The IPI-Urban is a weighted sum of the following eight factors, calculated using the analytical hierarchy process (AHP) [84]: land use, hydrogeological units, tectonic lineament density, drainage density, slope, and three other factors related to urban hydraulic and sanitation, sewer network density, stormwater network density and water supply network density.

4. Results and Discussion

4.1. Hydrogeochemical Approach

The physical and chemical parameters of groundwater samples collected during the fieldwork campaigns are summarized in the supplementary materials (Table S1). These groundwaters are slightly acidic with a median pH value of 6.2, and electrical conductivity (EC) values ranging from 95 to 1035 µS/cm, with a median value of 382 µS/cm, therefore, characterized as medium mineralized waters. These values are in good agreement with other data reported for this area [20,61,62]. The dissolved oxygen (DO) values in these groundwaters are in the range of 0.30–4.14 mg/L, with a median value of 1.75 mg/L. According to [85], these values are in accordance with those expected for most groundwaters (DO concentrations <5 mg/L, frequently <2 mg/L). The redox potential (Eh) measurements revealed that most of the studied groundwaters were in the range +0.1–+0.3 V, with a median of +0.15 V, which were values characteristic of most groundwaters [86] and within the water stability domain.
The major anions analyzed included bicarbonate (HCO3), sulphate (SO42−), chloride (Cl), and nitrate (NO3). The groundwaters showed the following global rating in terms of median concentrations (mg/L): HCO3 (45.0) > Cl (41.2) > NO3 (40.0) > SO42− (35.9). HCO3 has the highest standard deviation, while SO42− has the lowest, and 31% of the water samples have concentrations of NO3 > 50 mg/L. Considering the European guidelines for the use of groundwater for human consumption and the Portuguese legislation for irrigation activities [87], nearly one-third of these water samples exceed the parametric value for nitrate (50 mg/L). This trend corroborates the results presented by [20,61,62]. Moreover, [20] classified a peri-urban area of Vila do Conde, between Ave and Onda rivers, with a potential high to moderate risk of nitrate contamination, based on the SINTACS and IPNOA indexes.
The major cations studied were sodium (Na+), calcium (Ca2+), potassium (K+), and magnesium (Mg2+), and the global rating was the following in terms of median concentrations (mg/L): Na+ (37.0) > Ca2+ (13.0) > Mg2+ (7.9) > K+ (4.0). The highest and the lowest standard deviation belong to Na+ and Mg2+, respectively. Among the other cations, most of the samples have concentrations of aluminium (Al3+) and iron (Fe2+) below 50 µg/L and concentrations of ammonia (NH4) below 0.1 mg/L. Concerning the environmental isotopic approach, the median values for δ18O, δ2H, and 3H are, respectively, −4.42‰, −26.5‰ and 1.0 TU.
Considering the major anions and cations, we conclude that the hydrogeochemical signatures of these groundwaters are predominantly Cl-Na to Cl-SO4-Na, however, the following facies may occur, Cl-SO4-Na-Mg, Cl-SO4-Na-Ca, and Cl-SO4-Mg-Na (Figure 3a). Nevertheless, in the three fieldwork campaigns, five groundwater samples from boreholes drilled in diverse geological environments revealed a distinct hydrogeochemical facies in the bicarbonate domain: 110_17, HCO3-Ca (biotitic granite, medium to fine grained); 110_64, HCO3-Na to HCO3-Na-Mg (mica schists, metagraywackes, and paragneisses); 110_65, HCO3-Na (biotitic granite, medium to fine grained); 122_32, HCO3-Na-Ca (two-mica granite, medium to coarse grained); 122_34, HCO3-Ca-Mg (two-mica granite, medium to coarse grained). Moreover, two other borehole groundwater samples had HCO3-Cl-Na facies, i.e., 110_02 (geological contact between two-mica granite, medium to coarse grained and gneisses) and 110_05 (geological contact between two-mica granite, medium to coarse grained and migmatites). [11,60,61,62,71,88] achieved similar hydrogeochemical signatures for the same hydrogeological units.
In order to better understand the hydrochemical signatures of these waters and the main processes controlling water mineralization, several correlations were established. Figure 3b–d indicates the following:
  • The strongest linear correlations were obtained in the water samples belonging to the shallow groundwater systems, represented by dug wells, fountains and springs.
  • Cl and Mg seem to have the same source for the shallow and deep groundwater systems, while the relation Cl and Na seems to point to a unique source only for the shallow systems; however, the common source of Cl and Mg is not marine, and the sea spray seems to be partially responsible for the Cl and Na; moreover, the proximity of several data to the Mg/Na seawater line outlines a partially common marine origin for these parameters.
  • The isotopic signatures of the groundwaters characterize a meteoric origin, since most of the water samples are positioned very close to the Global Meteoric Water Line (GMWL), defined by [89] and later improved by [90,91,92], and similar to the isotopic composition of the precipitation water samples, from the Portuguese Isotopes in Precipitation Network [93]. From the isotopic point of view, the following two main groundwater clusters have been identified: Group I stands for groundwaters collected from dug wells and boreholes, which presents a more enriched isotopic composition, similar to the Porto precipitation [93], corresponding to normal deviations related with seasonal variations of δ18O and δ2H on precipitation; Group II is composed of the spring samples and presents more depleted δ18O and δ2H values, which may be attributed to the fact that they could be ascribed to random precipitation events, resulting into a direct infiltration of meteoric waters along the fractured granitic rocks;
  • The meteoric origin of the shallow groundwaters seems to be reinforced by the SO4-Ca facies, concerning the relationship with the precipitation data; therefore, the partial origin of SO4 should be atmospheric pollution, enriched in SO2 gases.
  • The common source of Cl and Mg must be anthropogenic, related to organic fertilizers, including sewage and livestock residues (liquid and semiliquid manure), and animal waste (e.g., bovines), in areas with a higher agricultural and/or livestock production activity, especially in the Vila do Conde, Trofa, and Northern Maia municipalities. Moreover, the good relationship between NO3 and Cl for the shallow groundwater systems shows their partially common source; in fact, NO3 is also an important constituent of fertilizers, either organic or synthetic, sewage, and animal and human wastes (e.g., [94]). The studied groundwaters do not present such a trend, indicating different sources for Cl and NO3. Several studies developed in this region have reached similar conclusions (e.g., [20,62,95,96,97]).
  • The relations among Cl and Mg, and Cl and NO3, may be ascribed to the urban and industrial sectors, particularly in the Porto, Matosinhos, and Vila Nova de Gaia municipalities, due to numerous groundwater potential contamination activities and their high density in some areas (cf. [17] for Porto and Vila Nova de Gaia urban areas), i.e., wastewater leakages, cesspools, and solid waste tanks contamination, hydrocarbons present from vehicle fuels and industrial processes, such as solvents and degreasing agents, namely trichloroethylene, which is one of the most common.
  • HCO3 and Ca have a reasonable correlation among the borehole water samples, and the correlation between HCO3 + NO3 and Na + Ca is good for the same water samples. This trend seems to indicate that these parameters have partially the same origin that probably should be ascribed to water–rock interaction, namely the hydrolysis of plagioclases presented in granitic rocks.
Therefore, the mineralization of these groundwaters seems to be controlled by the following three main sources/processes: (i) meteoric, including sea spray, contributing with the Cl and Na+ ions, and atmospheric pollution, particularly responsible by the presence of SO42−; (ii) anthropogenic contamination, particularly introducing Cl, Mg2+, and NO3; and (iii) water–rock interaction, adding HCO3, Na+, and Ca2+.
Since the altitude differential of the springs sampling sites is rather small, in the study region, (spring waters are the most representative of local precipitation), it was difficult to estimate the altitude of recharge areas ascribed to the groundwaters from boreholes (usually recharged far from the discharge). In fact, the correlation between the altitude of the spring sites and δ18O values was very low (r = −0.14).
Tritium concentrations, in the range 0.0–4.6 TU, are consistent with the precipitation 3H record measured at the Porto meteorological station, with a weight annual mean of 4.5 TU (see [100]). Therefore, for coastal areas, according to [49], the 3H values related to the studied groundwaters could correspond to a mixture of recent waters, with short and fast underground flow paths (<5–10 years) and relatively close to the surface, and submodern waters where recharge took place before the thermonuclear events in 1952.

4.2. Hydrodynamical Assessment

For this evaluation, the study area was divided into the following two key sectors (Figure 4a), each sector dominated by a hydrogeological unit: (i) Sector 1, in the geological contact of two-mica granite, medium to coarse grained and metasedimentary rocks, and also in the biotitic granite, medium to fine grained, with 16 boreholes; (ii) Sector 2, in the two-mica granite, medium to coarse grained, with 18 boreholes and five piezometers.
The median depth of the boreholes is 121 m in Sector 1, and higher in Sector 2, around 154 m (Figure 4b). For piezometers, the depth range is 15–34 m. The median value for the thickness of the weathering profile is null for both sectors, and it varies in the following ranges: 0.0–17.0 m in Sector 1, and 0.0–26 m in Sector 2. Concerning long-term well capacity, the minima, maxima, and median values are very similar in both sectors, i.e., 0.4 L/s, 3.5 L/s, and 1.1 L/s, respectively (Figure 4b). There was no correlation found between depth of boreholes and long-term well capacity, or between depth of screens and long-term well capacity.
For the assessment of transmissivity and storage coefficient values, we selected two boreholes and one piezometer, i.e., boreholes 110_02 and 122_09 in Sectors 1 and 2, respectively, and piezometer 122_23 in Sector 2 (Figure 5). Concerning borehole 110_02, the transmissivity values obtained during the pumping period were 21.8 m2/day and 34.4 m2/day, for the borehole and the observation well 110_03, respectively, while, for the recovery period, it was 24.3 m2/day. Additionally, the storage coefficient value was 3.1 × 10−4. Regarding borehole 122_09, the transmissivity values obtained during the pumping period were 2.8 m2/day and 6.7 m2/day for the borehole and the observation well 122_07, respectively, while, for the recovery period, it was 2.1 m2/day. Moreover, the storage coefficient value was 2.5 × 10−5. For piezometer 122_23, the transmissivity values obtained during the pumping period were 51.3 m2/day and 55 m2/day. Furthermore, the storage coefficient values were 3.1 × 10−2 and 3.9 × 10−3.
Due to the lack of well-documented aquifer conditions and complete pumping tests, a Logan approximation was achieved using specific capacity values for 69% of the boreholes in two sectors (Figure 6). According to [101], the Logan method can give a reasonable first estimate of transmissivity, provided that near equilibrium conditions are achieved when the maximum drawdown is measured. Therefore, Figure 6a shows good correlations among the Logan method and the other methods, particularly the Theis and the Cooper–Jacob methods. Additionally, transmissivity values are quite variable in both sectors, in the range of 0.1–26 m2/day. Nevertheless, the median values of transmissivity are very similar for the different methods and even for the two sectors, 1.8–4.0 m2/day (Sector 1) and 2.5–4.2 m2/day (Sector 2), (Figure 6b).
Considering the storage coefficient, values are also variable, in the range of 2.5 × 10−5–4 × 10−2, which characterizes confined to semi-confined groundwater conditions.
Very good correlations (R2 > 0.85) were found between long-term well capacities and transmissivities in both sectors. The highest values of long-term well capacity and transmissivity could be explained by the local presence of a very productive fracture, or to a direct connection with a nearby surface water body or a porous overlying aquifer.
All these conclusions confirm other investigations, performed by several authors, in fractured-rock aquifers in the Porto region and Portugal (e.g., [60,61,72,102,103] and references therein).
The hydrodynamic behavior of these aquifer systems meets the criteria mentioned in the literature for fractured-rock aquifers over the last four decades (e.g., [104,105,106,107,108,109,110,111,112,113]).

4.3. Urban Infiltration Potential Index (IPI-Urban)

The spatial distribution and explanation of six thematic layers (Figure 7) involved in the computation of the IPI-Urban (Figure 8), i.e., tectonic lineaments density, slope, drainage density, sewage network density, stormwater network density, and water supply system density are presented in this section.
Land use and hydrogeological units’ thematic layers were previously introduced in Section 2. Nevertheless, briefly, land use plays a crucial role in this region, and the following two classes dominate: urban and industrial areas (ca. 41%), and agricultural and forest areas (ca. 52%). This suggests that surfaces of low to very low permeability cover more than one-third of the region and must be responsible for the reduction in infiltration and, consequently, direct groundwater recharge. Hydrogeological units also play a key role in groundwater occurrence. Granitic rocks are the most representative lithology, namely the two-mica and biotitic granites (ca. 56%).
Concerning tectonic lineaments, most of the region (ca. 65%) has moderate to low densities (0.5–1.5 km/km2). The highest densities (>1.5 km/km2, ca. 20%) occur disseminated in the region, and the lowest densities (<0.5 km/km2, ca. 15%) are mostly concentrated in the west-northwestern (W-NW) area, where the sedimentary cover has its most representative domain.
As for slope, much of the region (ca. 81%) has very gentle to gentle gradients (<5°). Moderate and strong gradients (>5°) arise in almost 19% of the region and are mostly connected with the Ave, Leça and Douro river valleys.
The region has mostly (ca. 92%) a very low to low drainage density (<4 km/km2). The highest densities (4–8 km/km2) are closely related with the four major rivers and streams (Ave, Onda, Leça, and Douro rivers).
Regarding sewer and stormwater networks, the largest part of the region (ca. 58%) has low to moderate densities (5–15 km/km2). The highest densities (>15 km/km2) are concentrated (ca. 21%) in the most representative municipalities, namely Porto and Matosinhos.
Ultimately, the dominant water supply system classes (ca. 66%) are low to very Low (<12 km/km2). Most of the remaining parts of the region (ca. 25%) has moderate densities (12–18 km/km2) which are concentrated in Porto, Matosinhos, and Maia municipalities.
Regarding the IPI-Urban index, the overall scenario shows that in 96% of the urban region, the dominant classes are low (ca. 54%) and moderate (ca. 42%). The low class seems to be mostly connected with the two-mica and biotitic granites in areas where these hydrogeological units coexist with an agricultural and forest cover. However, the moderate class appears to be associated firstly with the sedimentary cover. Moreover, this class occurs in areas covered by an urban and industrial fabric. Although this may seem incoherent, the higher densities of the urban hydraulic and sanitation features (sewer, stormwater, and water supply networks) are found in these areas. Therefore, on a regional scale, the major role of these features is to enhance the natural low infiltration potential. The very low class occurs in approximately 4% of the region and is linked to the “Micaschist, metagreywacke and paragneiss” and “Micaschist, gneiss and migmatite” hydrogeological units, especially where these units are covered by urban and industrial areas, as well as the airport and harbor areas. The high class represents less than 1% of the region and corresponds to very confined areas, mainly between the Leça and Douro rivers and to the south of the Douro river.

5. Hydrogeological Conceptual Ground Models

The building of a conceptual hydrogeological model is presented as the first step in the entire process of modeling Earth systems. A multicriteria methodology supported in an accurate GIS mapping allows a perspective from geosciences for the paradigm of smart cities, namely through multidisciplinary, interdisciplinary, and transdisciplinary approaches (e.g., [35,36]).
In the study of crystalline aquifers, hydrogeomorphology has a fundamental role in the refinement of conceptual hydrogeological models, which aim at the sustainability of groundwater resources [17,18]. In this kind of system, there are many factors, some of them very complex, that control the infiltration and circulation of groundwater at the local/regional scale.
Figure 9 shows the hydrogeological conceptual models for the Porto urban region (Figure 9I), and the Porto urban area (Figure 9II). This conceptualization included the identification of different kinds of aquifer systems in an urban and pre-urban scenario. The main criteria considered were the hydrogeological units features, the morphostructure, the land use, the hydrogeochemical, isotopic and hydrodynamical characteristics, the infiltration potential, and the recharge processes (Table 1).
Regarding the Porto urban region, the following two perspectives are presented: a pre-urban context (Figure 9IA) and an urban context, with the main levels of land use (Figure 9IB), which is relevant to evaluate the impact on infiltration and direct recharge of these aquifer systems. For the Porto urban area, two viewpoints are shown, i.e., the main land use classes (Figure 9IIC) and the main IPI-Urban classes (Figure 9IID). For both views, several hydroclimatic and hydrodynamical parameters were added.
The analysis of both models allowed us to characterize the following three aquifer systems, i.e., superficial, intermediate, and deep, which are interconnected, sometimes in a discontinuous manner:
  • A superficial unit corresponds essentially to the sedimentary cover and the weathered/fractured zones of metasedimentary and granitic rocks that constitute a porous medium with hydraulic connection to the drainage system. The water table is close to the surface (<5 m). The more suitable exploitation structures are dug wells and springs associated, or not, to galleries, and the long-term well capacities are low (1 < Q < 2 L/s). Sedimentary cover can reach thicknesses of almost 30 m. In crystalline rocks, the weathering thickness is variable, and can reach values of 20–40 m locally (e.g., [116]), affecting transmissivity, which is generally low (<5 m2/day), and storage coefficient. The sedimentary cover constitutes an unconfined aquifer, while, in crystalline rocks, it corresponds to a semi-confined one. The hydrochemical facies is mostly Cl-Na. The infiltration potential is moderate to high and the groundwater recharge is direct, through infiltration of precipitation.
  • Intermediate aquifers constitute a fractured media, which may have a hydraulic connection to the drainage system. The more suitable exploitation structures are boreholes and the long-term well capacities are mostly very low (Q < 1 L/s). Transmissivity values are low (<5 m2/day) and these aquifers are semi-confined to confined. Groundwater has a short and shallow circuit with a Cl-Na to Cl-SO4-Na hydrochemical facies. The infiltration potential is moderate to low and the groundwater recharge occurs by leakage of the overlain levels or directly from the surface, namely by geological structures (e.g., geological contacts, faults, and veins), with favorable geo-hydraulic characteristics for the groundwater flow (e.g., major deep, open and not filled fractures, and intersections between tectonic lineaments).
  • Deep aquifers correspond to unweathered and massive crystalline bedrock, with closed fractures. They constitute a fissured media, where groundwater flow tends to have a weak regime with very low transmissivities and the hydraulic characteristics are confined.

6. Conclusions

Expanding our understanding of the behavior of fractured bedrock systems in urban environments will continue to be a research challenge requiring multidisciplinary approaches such as the one presented in this work.
The Porto urban region is the second largest metropolitan area in Portugal mainland, with a high population density, living in a land covered by urban and industrial areas (ca. 41%), mostly in the south part of the region, and agricultural and forest areas (ca. 52%), in the north.
The climatic conditions, the morphotectonic and geological background, mainly comprised of crystalline formations, namely granites and metasedimentary rocks, are responsible for the groundwater resources availability.
These groundwaters have median values for pH, electrical conductivity, dissolved oxygen, and redox potential of 6.2, 382 µS/cm, 1.75 mg/L, and +0.15 V, respectively. Considering major anions, the median concentrations for bicarbonate, sulphate, chloride, and nitrate are 45.0 mg/L, 35.9 mg/L, 41.2 mg/L, and 40.0 mg/L, respectively. In addition, for the major cations the median concentrations for sodium, calcium, potassium, and magnesium are 37.0 mg/L, 13.0 mg/L, 4.0 mg/L, and 7.9 mg/L, respectively. The hydrogeochemical facies are mainly Cl-Na to Cl-SO4-Na. Moreover, the median values for δ18O, δ2H, and 3H are −4.42‰, −26.5‰, and 1.0 TU, respectively. All these characteristics are mostly dependent on the proximity of the Porto urban region to the Atlantic Ocean, the anthropogenic and agricultural contamination processes, and, with a lesser impact, the water–rock interaction, since groundwater has fast, shallow flow paths, and short residence times.
This fractured environment is the key to several hydrodynamic characteristics, which characterize confined to semi-confined groundwater conditions, namely low to very low long-term well capacities, with a median value of 1 L/s; low transmissivities, with median values in the range 2–4 m2/day; and, generally, low storage coefficients, in the range 10−5–10−2.
The infiltration potential and the groundwater recharge are typically low in this fissured background. This low potential largely arises in zones where granitic rocks and agricultural/forest areas coexist. However, a moderate infiltration potential also occurs, naturally associated with the sedimentary cover, and artificially to the urban hydraulic and sanitation features.
The main conclusions summarized above contribute to outline recommendations for policymakers for the environmental protection, reasoning on the planning with nature and to the sustainable water resources management for the urban areas. In fact, integrated groundwater resources management is vital to evaluate water use, water availability, and to design balanced solutions in a changing society, environment, and climate.
Last but not least, in the current context of climate change, water shortage and groundwater depletion, along with a pandemic scenario, the authors support this impressive message by [117], “There is a need for the incorporation of groundwater in national IWRM plans and promotion of groundwater resource management in the programs of river basin agencies. Planned conjunctive management of groundwater and surface-water resources in many cases represents the best prospect for improving water-supply security for urban and irrigation use, and for sustainable resources.”

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4441/12/10/2797/s1, Table S1: Physical, chemical, and isotopic parameters for groundwater samples collected during the fieldwork campaigns.

Author Contributions

Conceptualization, M.J.A. and H.I.C.; methodology, M.J.A., L.F., J.M.M., P.M.C. and H.I.C; software, L.F. and M.J.A.; investigation, M.J.A., L.F., J.M.M., P.M.C., F.R., A.J.S.C.P. and H.I.C.; visualization, M.J.A., L.F. and H.I.C.; writing—original draft preparation, M.J.A., L.F., H.I.C, J.M.M., P.M.C., F.R., and A.J.S.C.P.; writing-review and editing, M.J.A., L.F., and H.I.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially financed by FEDER-EU COMPETE Funds and the Portuguese Foundation for the Science and Technology, FCT (UID/GEO/04035/2020, UID/Multi/00611/2020, and GroundUrban project POCI/CTE-GEX/59081/2004), and by the Labcarga|ISEP re-equipment program (IPP-ISEP|PAD’2007/08). The research was also funded by a doctoral scholarship from the Portuguese Foundation for Science and Technology (FCT) to L. Freitas (SFRH/BD/117927/2016). P.M. Carreira acknowledges the FCT support through the FCT-UIDB/04349/2020 project and J.M. Marques recognizes the FCT support through the UID/ECI/04028/2020 project.

Acknowledgments

This work was partially financed by FEDER-EU COMPETE Funds and the Portuguese Foundation for the Science and Technology, FCT (UID/GEO/04035/2020, UID/Multi/00611/2020 and GroundUrban project: POCI/CTE-GEX/59081/2004), and by the Labcarga|ISEP re-equipment program (IPP-ISEP|PAD’2007/08). The research was also funded by a doctoral scholarship from the Portuguese Foundation for Science and Technology (FCT) to L. Freitas (SFRH/BD/117927/2016) and, currently, by Project ReNATURE (Centro-01-0145-FEDER-000007-PT). P.M. Carreira acknowledge the FCT support through the FCT-UIDB/04349/2020 project and J.M. Marques recognize the FCT support through the UID/ECI/04028/2020 project. We are grateful to the anonymous reviewers for the comments that helped to improve the focus of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Foster, S.D.; Tyson, G. Resilient Cities & Groundwater; Strategic Overview Series; International Association of Hydrogeologists: Goring/Reading, UK, 2015. [Google Scholar]
  2. Howard, K.W.F.; Israfilov, R.G. Current Problems of Hydrogeology in Urban Areas, Urban. Agglomerates, and Industrial Centres; Kluwer: Dordrecht, The Netherlands, 2002. [Google Scholar]
  3. Howard, K.W.F. (Ed.) Urban. Groundwater. In Proceedings of the Meeting the Challenge: Selected Papers from the 32nd International Geological Congress (IGC), Florence, Italy, 20–28 August 2004; Taylor & Francis: London, UK, 2007. [Google Scholar]
  4. Massing, H.; Packman, J.; Zuidema, F.C. Hydrological Processes and Water Management in Urban Areas; IAHS Publication No. 198; IAHS: Wallingford, UK, 1990; p. 362. [Google Scholar]
  5. Johnson, S.P. The Earth Summit: The United Nations Conference on Environment and Development (UNCED); Graham & Trotman/Martinus Nijhoff: London, UK, 1993. [Google Scholar]
  6. Wilkinson, W.B. (Ed.) Groundwater Problems in Urban Areas; T. Telford: London, UK, 1994; p. 453. [Google Scholar]
  7. Chilton, J. (Ed.) Groundwater in the Urban Environment: Problems, Processes and Management. In Proceedings of the 27th Congress, International Association of Hydrogeologists, Nottingham, UK, 21–27 September 1997. [Google Scholar]
  8. Chilton, J. Groundwater in the Urban Environment: Selected City Profiles; A. A. Balkema: Rottterdam, The Netherlands, 1999; p. 342. [Google Scholar]
  9. Howard, K.W.F. Groundwater for Socio-Economic Development—The Role of Science; UNESCO IHP-VI Series on Groundwater, 9, Published as CD.; UNESCO: Paris, France, 2004; ISBN 92-9220-029-1. [Google Scholar]
  10. Tellam, J.H.; Rivett, M.O.; Israfilov, R.G. Urban. Groundwater Management and Sustainability; NATO Science Series, IV Earth and Environmental Sciences; Springer: Dordrecht, The Netherlands, 2006; Volume 74, p. 491. [Google Scholar]
  11. Afonso, M.J.; Chaminé, H.I.; Marques, J.M.; Carreira, P.M.; Guimarães, L.; Guilhermino, L.; Gomes, A.; Fonseca, P.E.; Pires, A.; Rocha, F. Environmental issues in urban groundwater systems: A multidisciplinary study of the Paranhos and Salgueiros spring waters, Porto (NW Portugal). Environ. Earth Sci. 2010, 61, 379–392. [Google Scholar] [CrossRef]
  12. Afonso, M.J.; Freitas, L.; Pereira, A.J.S.C.; Neves, L.J.P.F.; Guimarães, L.; Guilhermino, L.; Mayer, B.; Rocha, F.; Marques, J.M.; Chaminé, H.I. Environmental groundwater vulnerability assessment in urban water mines (Porto, NW Portugal). Water 2016, 8, 499. [Google Scholar] [CrossRef] [Green Version]
  13. Afonso, M.J.; Freitas, L.; Chaminé, H.I. Groundwater recharge in urban areas (Porto, NW Portugal): The role of GIS hydrogeology mapping. Sustain. Water Resour. Manag. 2019, 5, 203–216. [Google Scholar] [CrossRef]
  14. Sharp, J.M.; Hibbs, B.J. Special Issue on Hydrogeological Impacts of Urbanization. Environ. Eng. Geosci. 2012, 18, 111. [Google Scholar]
  15. Schirmer, M.; Leschik, S.; Musolff, A. Current research in urban hydrogeology—A review. Adv. Water Resour. 2013, 51, 280–291. [Google Scholar] [CrossRef]
  16. Freitas, L.; Afonso, M.J.; Devy-Vareta, N.; Marques, J.M.; Gomes, A.; Chaminé, H.I. Coupling hydrotoponymy and GIS cartography: A case study of hydrohistorical issues in urban groundwater systems, Porto NW Portugal. Geogr. Res. 2014, 52, 182–197. [Google Scholar] [CrossRef]
  17. Freitas, L.; Afonso, M.J.; Pereira, A.J.S.C.; Deleru-Matos, C.; Chaminé, H.I. Assessment of sustainability of groundwater in urban areas (Porto, NW Portugal): A GIS mapping approach to evaluate vulnerability, infiltration and recharge. Environ. Earth Sci. 2019, 78, 140. [Google Scholar] [CrossRef]
  18. Freitas, L.; Chaminé, H.I.; Pereira, A.J.S.C. Coupling groundwater GIS mapping and geovisualisation techniques in urban hydrogeomorphology: Focus on methodology. SN Appl. Sci. 2019, 1, 490. [Google Scholar] [CrossRef] [Green Version]
  19. Freitas, L.; Chaminé, H.I.; Afonso, M.J.; Meerkhan, H.; Abreu, T.; Trigo, J.F.; Pereira, A.J.S.C. Integrative groundwater studies in a small-scale urban area: Case study from the municipality of Penafiel (NW Portugal). Geosciences 2020, 10, 54. [Google Scholar] [CrossRef] [Green Version]
  20. Barroso, M.F.; Ramalhosa, M.J.; Olhero, A.; Antão, M.C.; Pina, M.F.; Guimarães, L.; Teixeira, J.; Afonso, M.J.; Delerue-Matos, C.; Chaminé, H.I. Assessment of groundwater contamination in an agricultural peri-urban area (NW Portugal): An integrated approach. Environ. Earth Sci. 2015, 73, 2881–2894. [Google Scholar] [CrossRef] [Green Version]
  21. Hibbs, B.J. Groundwater in urban areas. J. Contemp. Water Res. Educ. 2016, 159, 143. [Google Scholar] [CrossRef] [Green Version]
  22. Guimarães, L.; Guilhermino, L.; Afonso, M.J.; Marques, J.M.; Chaminé, H.I. Assessment of urban groundwater: Towards integrated hydrogeological and effects-based monitoring. Sustain. Water Resour. Manag. 2019, 5, 217–233. [Google Scholar] [CrossRef]
  23. Sherlock, R.L. Man as a Geological Agent; H.F. & G. Witherby: London, UK, 1922; p. 372. [Google Scholar]
  24. Legget, R.F. Cities and Geology; McGraw-Hill: New York, NY, USA, 1973; p. 579. [Google Scholar]
  25. Dunne, T.; Leopold, L.B. Water in Environmental Planning; W.H. Freeman Co.: San Francisco, CA, USA, 1978; p. 818. [Google Scholar]
  26. Utgard, R.O.; McKenzie, G.D.; Foley, D. Geology in the Urban. Environment; Burgess Publishing Company: Minneapolis, MI, USA, 1978; p. 355. [Google Scholar]
  27. Leveson, D. Geology and the Urban. Environment; Oxford University Press: Oxford, UK, 1980; p. 386. [Google Scholar]
  28. Zaporozec, A. (Ed.) Cities and water. Geo J. 1985, 11, 203–283. [Google Scholar]
  29. McCall, G.J.; Demulder, E.; Marker, B.R. Urban. Geoscience; AGID Special Publication Series; Taylor & Francis: Rotterdam, The Netherlands, 1996; Volume 20, p. 279. [Google Scholar]
  30. Gehrels, H.; Peters, N.E.; Hoehn, E.; Jensen, K.; Leibundgut, C.; Griffioen, J.; Webb, B.; Zaadnoordijk, W.J. Impact of Human Activity on Groundwater Dynamics; International Association of Hydrological Sciences, Publication no. 269; IAHS Press: Wallingford, UK, 2001; p. 369. [Google Scholar]
  31. Pokrajac, D. (Ed.) Groundwater in urban areas. Urban. Water 2001, 3, 171–237. [Google Scholar] [CrossRef]
  32. Bocanegra, E.; Hernández, M.; Usunoff, E. Groundwater and Human Development. International Association of Hydrogeologists Selected Papers; Taylor & Francis: London, UK, 2005; Volume 6, p. 262. [Google Scholar]
  33. Ehlen, J.; Haneberg, W.C.; Larson, R.A. Humans as Geologic Agents. Reviews in Engineering Geology; The Geological Society of America: Boulder, CO, USA, 2005; Volume 16, p. 158. [Google Scholar]
  34. Culshaw, M.G.; Reeves, H.J.; Jefferson, I.; Spink, T.W. Engineering Geology for Tomorrow’s Cities; Engineering Geology Special Publications; Geological Society: London, UK, 2009; Volume 22, p. 315. [Google Scholar]
  35. Chaminé, H.I.; Afonso, M.J.; Freitas, L. From historical hydrogeological inventory through GIS mapping to problem solving in urban groundwater systems. Eur. Geol. J. 2014, 38, 33–39. [Google Scholar]
  36. Chaminé, H.I.; Teixeira, J.; Freitas, L.; Pires, A.; Silva, R.S.; Pinho, T.; Mon- teiro, R.; Costa, A.L.; Abreu, T.; Trigo, J.F. From engineering geosciences mapping towards sustainable urban planning. Eur. Geol. J. 2016, 41, 16–25. [Google Scholar]
  37. UN-Habitat [United Nations Human Settlements Programme]. Urbanization and Development: Emerging Futures—World Cities Report 2016; United Nations Human Settlements Programme, World Urban Forum edition: Nairobi, Kenya, 2016; p. 264. [Google Scholar]
  38. Underwood, J.R. Anthropic rocks as a fourth basic class. Environ. Eng. Geosci. 2001, 7, 104–110. [Google Scholar] [CrossRef]
  39. Baker, L.A. (Ed.) The Water Environment of Cities; Springer Science & Business Media: London, UK, 2009; p. 307. [Google Scholar]
  40. Sharp, J.M. The impacts of urbanization on groundwater systems and recharge: Aqua Mundi. Aqua Mundi 2010, 1, 51–56. [Google Scholar] [CrossRef]
  41. Gogu, C.R.; Campbell, D.; de Beer, J. The Urban Subsurface—From Geoscience and Engineering to Spatial Planning and Management. Procedia Eng. 2017, 209, 224. [Google Scholar] [CrossRef]
  42. Hibbs, B.J.; Sharp, J.M. Hydrogeological impacts of urbanization. Environ. Eng. Geosci. 2012, 18, 3–24. [Google Scholar] [CrossRef]
  43. Taylor, R.G.; Scanlon, B.; Döll, P.; Rodell, M.; van Beek, R.; Wada, Y.; Longuevergne, L.; Leblanc, M.; Famiglietti, J.S.; Edmunds, M.; et al. Ground water and climate change. Nat. Clim. Chang. 2013, 3, 322–329. [Google Scholar] [CrossRef] [Green Version]
  44. Foster, S.D.; Tyson, G. Global Change & Groundwater; Strategic Overview Series; International Association of Hydrogeologists: Goring/Reading, UK, 2016. [Google Scholar]
  45. Wiles, T.J.; Sharp, J.M. The secondary permeability of impervious cover. Environ. Eng. Geosci. 2008, 14, 251–265. [Google Scholar] [CrossRef]
  46. Pujades, E.; De Simone, S.; Carrera, J.; Vázquez-Suñé, E.; Jurado, A. Settlements around pumping wells: Analysis of influential factors and a simple calculation procedure. J. Hydrol. 2017, 548, 225–236. [Google Scholar] [CrossRef]
  47. Lyu, H.M.; Shen, S.-L.; Yang, J.; Yin, Z.-Y. Inundation analysis of metro systems with the storm water management model incorporated into a geographical information system: A case study in Shanghai. Hydrol. Earth Syst. Sci. 2019, 23, 4293–4307. [Google Scholar] [CrossRef] [Green Version]
  48. IAEA [International Atomic Energy Agency]. Stable Isotope Hydrology. In Deuterium and Oxygen-18 in the Water Cycle; Technical Reports Series 210; IAEA: Vienna, Austria, 1981. [Google Scholar]
  49. Clark, I.D.; Fritz, P. Environmental Isotopes in Hydrogeology; CRC Press: Boca Raton, FL, USA; Lewis Publishers: Boca Raton, FL, USA, 1997; p. 328. [Google Scholar]
  50. Hoefs, J. Stable Isotope Geochemistry; Completely Revised. Updated and Enlarged Edition; Springer: Berlin, Germany, 1997. [Google Scholar]
  51. Şener, E.; Davraz, A.; Ozcelik, M. An integration of GIS and remote sensing in groundwater investigations: A case study in Burdur, Turkey. Hydrogeol. J. 2005, 13, 826–834. [Google Scholar] [CrossRef]
  52. Thapa, R.; Gupta, S.; Gupta, A.; Reddy, D.V.; Kaur, H. Use of geospatial technology for delineating groundwater potential zones with an emphasis on water-table analysis in Dwarka River basin, Birbhum, India. Hydrogeol. J. 2018, 26, 899–922. [Google Scholar] [CrossRef]
  53. Morris, B.L.; Litvak, R.G.; Ahmed, K.M. Urban groundwater protection and management: Lessons from developing cities in Bangladesh and Kyrghyztan. In Current Problems of Hydrogeology in Urban Areas, Urban Agglomerates and Industrial Centres; NATO Science Series, IV Earth and Environmental Sciences; Howard, K.W.F., Israfilov, R.G., Eds.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2002; Volume 8, pp. 77–102. [Google Scholar]
  54. Sharp, J.M. Ground-water supply issues in urban and urbanizing areas. In Groundwater in the Urban Environment: Problems, Process and Management; Chilton, J., Ed.; A. A. Balkema: Rotterdam, The Netherlands, 1997; pp. 67–74. [Google Scholar]
  55. Barrett, M.H.; Howard, A.G. Urban groundwater and sanitation: Developed and developing countries. In Current Problems of Hydrogeology in Urban Areas, Urban Agglomerates and Industrial Centres; NATO Science Series, IV Earth and Environmental Sciences; Howard, K.W.F., Israfilov, R.G., Eds.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2002; Volume 8, pp. 39–56. [Google Scholar]
  56. Foster, S.D.; Hirata, R.; Howard, K.W.F. Groundwater use in developing cities: Policy issues arising from current trends. Hydrogeol. J. 2011, 19, 271–274. [Google Scholar] [CrossRef]
  57. Elshall, A.S.; Arik, A.D.; El-Kadi, A.I.; Pierce, S.; Ye, M.; Burnett, K.M.; Wada, C.A.; Bremer, L.L.; Chun, G. Groundwater sustainability: A review of the interactions between science and policy. Environ. Res. Lett 2020, 15, 093004. [Google Scholar] [CrossRef]
  58. INE—Instituto Nacional de Estatística. Statistical Information about Portuguese Population: Porto City. 2011. Available online: http://www.ine.pt/ (accessed on 30 January 2019).
  59. Caetano, M.; Igreja, C.; Marcelino, F.; Costa, H. Estatísticas e Dinâmicas Territoriais Multiescala de Portugal Continental 1995–2007–2010 com Base na Carta de Uso e Ocupação do Solo (COS). Relatório Técnico. Direção-Geral do Território (DGT). 2017. Available online: http://www.dgterritorio.pt/ (accessed on 28 February 2019).
  60. Afonso, M.J. Hidrogeologia de rochas graníticas da região do Porto (NW de Portugal). Cad. Lab. Xeol. Laxe 2003, 28, 173–192. [Google Scholar]
  61. Afonso, M.J. Hidrogeologia e Hidrogeoquímica da Região Litoral Urbana do Porto, entre Vila do Conde e Vila Nova de Gaia (NW de Portugal): Implicações Geoambientais. Ph.D. Thesis, Instituto Superior Técnico da Universidade Técnica de Lisboa, Lisboa, Portugal, 2011. [Google Scholar]
  62. Afonso, M.J.; Chaminé, H.I.; Carvalho, J.M.; Marques, J.M.; Gomes, A.; Araújo, M.A.; Fonseca, P.E.; Teixeira, J.; Marques da Silva, M.A.; Rocha, F.T. Urban groundwater resources: A case study of Porto City in northwest Portugal. In Urban Groundwater: Meeting the Challenge. International Association of Hydrogeologists Selected Papers; Howard, K.W.F., Ed.; Taylor & Francis Group: London, UK, 2007; Volume 8, pp. 271–287. [Google Scholar]
  63. Pereira, E.; Ribeiro, A.; Carvalho, G.S.; Noronha, F.; Ferreira, N.; Monteiro, J.H. Carta Geológica de Portugal, escala 1/200000. Folha 1; Serviços Geológicos de Portugal: Lisboa, Portugal, 1989. [Google Scholar]
  64. Chaminé, H.I.; Gama Pereira, L.C.; Fonseca, P.E.; Noronha, F.; Lemos de Sousa, M.J. Tectonoestratigrafia da faixa de cisalhamento de Porto–Albergaria-a-Velha–Coimbra–Tomar, entre as Zonas Centro-Ibérica e de Ossa-Morena (Maciço Ibérico, W de Portugal). Cad. Lab. Xeol. Laxe 2003, 28, 37–78. [Google Scholar]
  65. Chaminé, H.I.; Afonso, M.J.; Robalo, P.M.; Rodrigues, P.; Cortez, C.; Mon- teiro Santos, F.A.; Plancha, J.P.; Fonseca, P.E.; Gomes, A.; Devy-Vareta, N.F.; et al. Urban speleology applied to groundwater and geo-engineering studies: Underground topographic surveying of the ancient Arca D’Água galleries catchworks (Porto, NW Portugal). Int. J. Speleol. 2010, 39, 1–14. [Google Scholar] [CrossRef] [Green Version]
  66. Araújo, M.A.; Gomes, A.; Chaminé, H.I.; Fonseca, P.E.; Gama Pereira, L.C.; Pinto de Jesus, A. Geomorfologia e geologia regional do sector de Porto-Espinho (W de Portugal): Implicações morfoestruturais na cobertura sedimentar Cenozóica. Cad. Lab. Xeol. Laxe 2003, 28, 79–105. [Google Scholar]
  67. Struckmeier, W.F.; Margat, J. Hydrogeological Maps: A Guide and a Standard Legend; International Association of Hydrogeologists: Hannover, Germany, 1995; p. 177. [Google Scholar]
  68. Assaad, F.A.; LaMoreaux, P.E.; Hughes, T.H.; Wangfang, Z.; Jordan, H. Field Methods for Geologists and Hydrogeologist; Springer: Berlin, Germany, 2004; p. 420. [Google Scholar]
  69. Carríngton da Costa, J.; Teixeira, C. Carta Geológica de Portugal na escala de 1/50000. Notícia Explicativa da Folha 9-C (Porto); Serviços Geológicos de Portugal: Lisboa, Portugal, 1957; p. 38. [Google Scholar]
  70. Teixeira, C.; Medeiros, A.C. Carta geológica de Portugal na escala 1:50000. Notícia explicativa da folha 9A-Póvoa de Varzim; Serviços Geológicos de Portugal: Lisboa, Portugal, 1965. [Google Scholar]
  71. Pedrosa, M.Y. Notícia explicativa da Carta Hidrogeológica de Portugal, à escala 1/200000. Folha 1; Instituto Geológico e Mineiro: Lisboa, Portugal, 1999; p. 70. [Google Scholar]
  72. Carvalho, J.M. Prospecção e Pesquisa de Recursos Hídricos Subterrâneos no Maciço Antigo Português: Linhas Metodológicas. Ph.D. Thesis, Universidade de Aveiro, Aveiro, Portugal, 2006; p. 292. [Google Scholar]
  73. Epstein, S.; Mayeda, T. Variations of 18O content of waters from natural sources. Geochim. Cosmochim. Acta 1953, 4, 213–224. [Google Scholar] [CrossRef]
  74. Friedman, I. Deuterium content of natural waters and other substances. Geochim. Cosmochim. Acta 1953, 4, 89–103. [Google Scholar] [CrossRef]
  75. IAEA [International Atomic Energy Agency]. Procedure and Technique Critique for Tritium Enrichment by Electrolysis at IAEA Laboratory; Technical Procedure No. 19; IAEA: Vienna, Austria, 1976. [Google Scholar]
  76. Lucas, L.L.; Unterweger, M.P. Comprehensive review and critical evaluation of the half-life of tritium. J. Res. Natl. Inst. Stand. Technol. 2000, 105, 541–549. [Google Scholar] [CrossRef]
  77. Carreira, P.M.; Marques, J.M.; Graça, R.C.; Aires-Barros, L. Radiocarbon application in dating “complex” hot and cold CO2-rich mineral water systems: A review of case studies ascribed to the northern Portugal. Appl. Geochem. 2008, 23, 2817–2828. [Google Scholar] [CrossRef]
  78. Carreira, P.M.; Marques, J.M.; Carvalho, M.R.; Capasso, G.; Grassa, F. Mantle-derived carbon in Hercynian granites Stable isotopes signatures C/He associations in the thermomineral waters N-Portugal. J. Volcanol. Geotherm. Res. 2010, 189, 49–56. [Google Scholar] [CrossRef]
  79. Theis, C.V. The relation between lowering of the piezometric surface and rate and duration of discharge of a well using ground-water storage. Trans. Am. Geophys. Union 1935, 16, 519–524. [Google Scholar] [CrossRef]
  80. Cooper, H.H.J.R.; Jacob, C.E. A generalized graphical method for evaluating formation constants and summarizing well-field history. Trans. Am. Geophys. Union 1946, 27, 526–534. [Google Scholar] [CrossRef]
  81. Kruseman, G.P.; de Ridder, N.A. Analysis and Evaluation of Pumping Test. Data, 2nd ed.; International Institute for Land Reclamation and Improvement: Wageningen, The Netherlands, 1990; Volume 47, p. 377. [Google Scholar]
  82. Sterrett, R.J. Groundwater and Wells, 3rd ed.; Johnson Screens, A Weatherford Company: New Brighton, MN, USA, 2007; p. 812. [Google Scholar]
  83. Logan, J. Estimating transmissibility from routine production tests of water wells. Ground Water 1964, 2, 35–37. [Google Scholar] [CrossRef]
  84. Malczewski, J.; Rinner, C. Multicriteria Decision Analysis in Geographic Information Science; Springer: New York, NY, USA, 2015. [Google Scholar]
  85. Custodio, E.; Llamas, M.R. Hidrología Subterránea; Segunda Edición Corregida; Ediciones Omega: Barcelona, Spain, 2001; p. 2350. [Google Scholar]
  86. Younger, P.L. Groundwater in the Environment: An Introduction; Blackwell Publishing: Hoboken, NJ, USA, 2007; p. 318. [Google Scholar]
  87. MA—Ministério do Ambiente. Decreto-Lei nº 236/98, de 1 de Agosto. Diário da República—I Série-A, Nº 176; Ministério do Ambiente: Lisbon, Portugal, 1998.
  88. Pedrosa, M.Y. Carta Hidrogeológica de Portugal, à escala 1/200000. Folha 1; Instituto Geológico e Mineiro: Lisboa, Portugal, 1998. [Google Scholar]
  89. Craig, H. Isotopic variations in meteoric waters. Science 1961, 133, 1702–1703. [Google Scholar] [CrossRef] [PubMed]
  90. Rozanski, K.; Araguás-Araguás, L.; Gonfiantini, R. Isotopic patterns in modern global precipitation. In Climate Change in Continental Isotopic Records, Geoph. Monog. Series; Swart, P.K., Ed.; AGU: Washington, DC, USA, 1993; pp. 1–36. [Google Scholar]
  91. Bowen, G.J.; Wilkinson, B. Spatial distribution of δ18O in meteoric precipitation. Geology 2002, 30, 315–318. [Google Scholar] [CrossRef]
  92. Terzer, S.; Wassenaar, L.I.; Araguás-Araguás, L.J.; Aggarwal, P.K. Global isoscapes for δ18O and δ2H in precipitation: Improved prediction using regionalized climatic regression models. Hydrol. Earth. Syst. Sci. 2013, 17, 4713–4728. [Google Scholar] [CrossRef]
  93. Carreira, P.M.; Araújo, M.F.; Nunes, D. Isotopic composition of rain and water vapour samples from Lisbon region: Characterization of monthly and daily events. In IAEA-TECDOC-1453 Isotopic Composition of Precipitation in the Mediterranean Basin in Relation to Air Circulation Patterns and Climate; IAEA: Vienna, Austria, 2005; pp. 141–155. [Google Scholar]
  94. Wakida, F.; Lerner, D. Non-agricultural sources of groundwater nitrate: A review and case study. Water Res. 2005, 39, 3–16. [Google Scholar] [CrossRef] [PubMed]
  95. Heitor, A.M.F. Contaminação das águas subterrâneas no Norte de Portugal. In Las Aguas Subterráneas en el Noroeste de la Península Ibérica; Textos de las Jornadas, Mesa Redonda y Comunicaciones; Samper, J., Leitão, T., Fernández, L., Ribeiro, L., Eds.; A Coruña. AIH-Grupo Español & APRH. ITGE: Madrid, Spain, 2000; pp. 295–308. [Google Scholar]
  96. Pedrosa, M.Y.; Brites, J.A.; Pereira, A.P. Carta das Fontes e do Risco de Contaminação da Região de Entre-Douro-e-Minho. Folha Sul, Escala 1/100000, Nota Explicativa; Instituto Geológico e Mineiro: Lisboa, Portugal, 2002. [Google Scholar]
  97. Correia, M.; Barroso, A.; Barroso, M.F.; Soares, D.; Oliveira, M.B.P.P.; Delerue-Matos, C. Contribution of different vegetable types to exogenous nitrate and nitrite exposure. Food Chem. 2010, 120, 960–966. [Google Scholar] [CrossRef] [Green Version]
  98. Begonha, A.; Sequeira Braga, M.A.; Gomes da Silva, F. A acção da água da chuva na meteorização de monumentos graníticos. In IV Congresso Nacional de Geologia. Memórias do Museu e Laboratório Mineralógico e Geológico da Faculdade de Ciências da Universidade do Porto; Borges, F.S., Marques, M., Eds.; MLMGFCUP: Porto, Portugal, 1995; Volume 4, pp. 177–181. [Google Scholar]
  99. Begonha, A. Meteorização do Granito e Deterioração da Pedra em Monumentos e Edifícios da Cidade do Porto; Colecção monografias, FEUP Edições: Porto, Portugal, 2001; p. 445. [Google Scholar]
  100. Carreira, P.M.; Valério, P.; Nunes, D.; Araújo, M.F. Temporal and seasonal variations of stable isotopes (2H and 18O) and tritium in precipitation over Portugal. In Proceedings of the Isotopes in Environmental Studies—Aquatic Forum, Monte Carlo, Monaco, 25–29 October 2004; IAEA: Vienna, Austria, 2006; pp. 370–373. [Google Scholar]
  101. Misstear, B.D.R. The value of simple equilibrium approximations for analysing pumping test data. Hydrogeol. J. 2001, 9, 125–126. [Google Scholar] [CrossRef]
  102. Afonso, M.J.; Carvalho, J.M.; Marques, J.M.; Chaminé, H.I. Hydrodynamic constraints of the Porto urban area crystalline bedrock (NW Portugal, Iberian Massif): Implications on groundwater resources. In Proceedings of the 7th Hellenic Hydrogeological Conference and 2nd MEM Workshop on Fissured Rocks Hydrology, Athens, Greece, 5–6 October 2005; Stournaras, G., Pavlopoulos, K., Bellos, T.h., Eds.; The Geological Society of Greece (Hellenic Committee of Hydrogeology): Athens, Greece, 2005; Volume 2, pp. 77–81. [Google Scholar]
  103. Carvalho, J.M.; Chaminé, H.I.; Afonso, M.J.; Espinha Marques, J.; Medeiros, A.; Garcia, S.; Gomes, A.; Teixeira, J.; Fonseca, P.E. Productivity and water cost in fissured-aquifers from the Iberian crystalline basement (Portugal): Hydrogeological constraints. In Water, Mining and Environment. Book Homage to Professor Rafael Fernández Rubio; López-Geta, J.A., et al., Eds.; Instituto Geológico y Minero de España: Madrid, Spain, 2005; pp. 193–207. [Google Scholar]
  104. Larsson, I. Groundwater in Hard Rocks. Studies and Reports in Hydrology; UNESCO: Paris, France, 1984; p. 234. [Google Scholar]
  105. Wright, E.P.; Burgess, W.G. The Hydrogeology of Crystalline Basement in Africa; Geological Society Special Publication, 68; GSL: London, UK, 1992; p. 264. [Google Scholar]
  106. Lloyd, J.W. Water Resources of Hard Rock Aquifers in Arid and Semi-arid Zones; Studies and Reports in Hydrology, 58; UNESCO: Paris, France, 1999; p. 284. [Google Scholar]
  107. Robins, N.S.; Misstear, B.D.R. Groundwater in the Celtic Regions: Studies in Hard Rock and Quaternary Hydrogeology; Geological Society of London: London, UK, 2000; p. 273. [Google Scholar]
  108. Stober, I.; Bucher, K. Hydrogeology of Crystalline Rocks; Water Science and Technology Library; Kluwer Academic Publishers: Dordrecht, The Netherlands, 2000; Volume 34, p. 284. [Google Scholar]
  109. Krásný, J.; Sharp, J.M. Hydrogeology of fractured rocks from particular fractures to regional approaches: State-of-the-art and future challenges. In Groundwater in Fractured Rocks. International Association of Hydrogeologists Selected Papers; Krásný, J., Sharp, J.M., Eds.; Taylor & Francis Group: London, UK, 2007; pp. 1–30. [Google Scholar]
  110. Singhal, B.B.S.; Gupta, R.P. Applied Hydrogeology of Fractured Rocks, 2nd ed.; Springer: Dordrecht, The Netherlands, 2010; p. 408. [Google Scholar]
  111. Gustafson, G. Hydrogeology for Rock Engineers; BeFo: Stockholm, Sweden, 2012. [Google Scholar]
  112. Sharp, J.M. Fractured Rock Hydrogeology; CRC Press: Boca Raton, FL, USA, 2014; p. 408. [Google Scholar]
  113. Ofterdinger, U.; Macdonald, A.M.; Comte, J.C.; Young, M.E. (Eds.) Groundwater in Fractured Bedrock Environments: Managing Catchment and Subsurface Resources; Geological Society: London, UK, 2019. [Google Scholar]
  114. Freitas, L. Avaliação Integrada de Recursos Hídricos em Áreas Urbanas: Aplicações Para a Sustentabilidade e o Ordenamento Territorial. Ph.D. Thesis, Faculdade de Ciências e Tecnologia da Universidade de Coimbra, Coimbra, Portugal, 2019; p. 425. [Google Scholar]
  115. Carvalho, J.M.; Espinha Marques, J.; Afonso, M.J.; Chaminé, H.I. Prospecção e pesquisa de recursos hidrominerais e de água de nascente no Maciço Antigo Português. Boletim e Minas 2007, 42, 161–196. [Google Scholar]
  116. Begonha, A.; Sequeira Braga, M.A. Weathering of the Oporto granite: Geotechnical and physical properties. Catena 2002, 49, 57–76. [Google Scholar] [CrossRef]
  117. IAH [International Association of Hydrogeologists]. The UN-SDGs for 2030: Essential Indicators for Groundwater. 2017. Available online: https://iah.org/wp-content/uploads/2017/04/IAH-Groundwater-SDG-6-Mar-2017.pdf (accessed on 30 June 2020).
Figure 1. Land use of the urban coastal region between the Municipalities of Vila do Conde and Vila Nova de Gaia (adapted from [59]).
Figure 1. Land use of the urban coastal region between the Municipalities of Vila do Conde and Vila Nova de Gaia (adapted from [59]).
Water 12 02797 g001
Figure 2. Regional hydrogeological setting of the Porto urban region (NW Portugal). Geological background updated from [63,65,69,70]. Hydrogeological framework adapted from [61,71,72]. Hydrogeological inventory and groundwater sampling points from [61].
Figure 2. Regional hydrogeological setting of the Porto urban region (NW Portugal). Geological background updated from [63,65,69,70]. Hydrogeological framework adapted from [61,71,72]. Hydrogeological inventory and groundwater sampling points from [61].
Water 12 02797 g002
Figure 3. Hydrochemical and isotopic relationships for groundwater samples collected during the three fieldwork campaigns. (a) Piper diagram; (b) Scatter diagrams for several hydrochemical parameters, R-squared value (R2) is indicated for each studied linear relationship, and seawater relations are those reported by [85]; (c) The average chemical composition of Porto precipitation, reported by [98] and [99]; (d) δ18O vs. δ2H signatures of the studied groundwater samples. The global meteoric water line (GMWL) (δ2H = 8δ18O + 10, [89]), the local (Porto meteorological station) meteoric water line (LMWL) (δ2H = (6.20 ± 0.33) and δ18O + (1.12 ± 1.16), r = 0.93, [93]), and the long-term mean isotopic composition of Porto precipitation (Prec. Porto) (δ2H = −26.9‰ and δ18O = 4.54 [93]) were plotted as reference.
Figure 3. Hydrochemical and isotopic relationships for groundwater samples collected during the three fieldwork campaigns. (a) Piper diagram; (b) Scatter diagrams for several hydrochemical parameters, R-squared value (R2) is indicated for each studied linear relationship, and seawater relations are those reported by [85]; (c) The average chemical composition of Porto precipitation, reported by [98] and [99]; (d) δ18O vs. δ2H signatures of the studied groundwater samples. The global meteoric water line (GMWL) (δ2H = 8δ18O + 10, [89]), the local (Porto meteorological station) meteoric water line (LMWL) (δ2H = (6.20 ± 0.33) and δ18O + (1.12 ± 1.16), r = 0.93, [93]), and the long-term mean isotopic composition of Porto precipitation (Prec. Porto) (δ2H = −26.9‰ and δ18O = 4.54 [93]) were plotted as reference.
Water 12 02797 g003aWater 12 02797 g003bWater 12 02797 g003c
Figure 4. Inventory of boreholes and piezometers used for the hydrodynamical assessment. (a) Study sites (Sectors 1 and 2) are shown on the right side of the figure; (b) Depth of boreholes and long-term well capacity for Sectors 1 and 2.
Figure 4. Inventory of boreholes and piezometers used for the hydrodynamical assessment. (a) Study sites (Sectors 1 and 2) are shown on the right side of the figure; (b) Depth of boreholes and long-term well capacity for Sectors 1 and 2.
Water 12 02797 g004
Figure 5. Interpretation of pumping tests. (I) 110_02, Cooper–Jacob method: (a) Pumping period data for 110_02; (b) Pumping period data for the observation well 110_03, located 370 m from 110_02; (c) Theis recovery method; (II) 122_09, Cooper–Jacob method: (a) Pumping period data for 122_09; (b) Pumping period data for the observation well 122_07, located 128 m from 122_09; (c) Theis recovery method. (III) 122_23, Cooper–Jacob method for two observation wells: (a) Located 16 m (b) Located 63 m.
Figure 5. Interpretation of pumping tests. (I) 110_02, Cooper–Jacob method: (a) Pumping period data for 110_02; (b) Pumping period data for the observation well 110_03, located 370 m from 110_02; (c) Theis recovery method; (II) 122_09, Cooper–Jacob method: (a) Pumping period data for 122_09; (b) Pumping period data for the observation well 122_07, located 128 m from 122_09; (c) Theis recovery method. (III) 122_23, Cooper–Jacob method for two observation wells: (a) Located 16 m (b) Located 63 m.
Water 12 02797 g005
Figure 6. Transmissivity values for Sectors 1 and 2. (a) Logan vs. other methods; (b) Box-and-whisker diagrams for Sectors 1 and 2.
Figure 6. Transmissivity values for Sectors 1 and 2. (a) Logan vs. other methods; (b) Box-and-whisker diagrams for Sectors 1 and 2.
Water 12 02797 g006
Figure 7. Tectonic lineaments density, slope, drainage density, sewage network density, stormwater network density, and water supply system density in the Porto urban region (NW Portugal).
Figure 7. Tectonic lineaments density, slope, drainage density, sewage network density, stormwater network density, and water supply system density in the Porto urban region (NW Portugal).
Water 12 02797 g007
Figure 8. Urban infiltration potential index (IPI-Urban) in the Porto urban region (NW Portugal).
Figure 8. Urban infiltration potential index (IPI-Urban) in the Porto urban region (NW Portugal).
Water 12 02797 g008
Figure 9. Urban hydrogeological conceptual models (revised and updated from [61,114]. (I) Porto urban region: (A) Pre-urban scenario; (B) Urban scenario). (II) Porto urban area: (C) Land use context; (D) IPI-Urban framework.
Figure 9. Urban hydrogeological conceptual models (revised and updated from [61,114]. (I) Porto urban region: (A) Pre-urban scenario; (B) Urban scenario). (II) Porto urban area: (C) Land use context; (D) IPI-Urban framework.
Water 12 02797 g009
Table 1. Regional hydrogeological units (cf. Figure 2) and related features in the Porto urban region (summarized and updated from [17,18,20,61,62,71,72,103,115]).
Table 1. Regional hydrogeological units (cf. Figure 2) and related features in the Porto urban region (summarized and updated from [17,18,20,61,62,71,72,103,115]).
Hydrogeological GroupsSedimentary CoverMetasedimentary RocksGranitic Rocks
Regional Hydrogeological Units (RHU)Beach and Dune Sands; AlluviaSandy Silts and ClaysMicaschists, Metagraywackes and ParagneissesTwo-mica Granite, Medium to Coarse GrainedBiotitic Granite, Medium to Fine GrainedMicaschists, Gneisses and MigmatitesAplite-pegmatite
and Quartz Veins
HYDROGEOLOGICAL FEATURESThickness (m)<1210–20not applicable
Weathering
profile
low (m)not applicable10–20 5–10<5
high (m) 20–4020–40
silty and/or clayeyX XXX
sandy XXXX
Connectivity to the drainage systemwithXX
possible XXXXX
Type of media
flow
porousXX
fissured XXXXX
Hydrochemical faciesCl-Na to NO3-NaCl-Na to Cl-SO4-Na
Environmental isotopesδ18O (‰)not determined−6.5 to −3.5not determined
δ2H (‰)−40 to −20
3H (TU)< 5not determined
Hydrodynamic parameterslong-term well capacity,
Q (L/s)
very low
(Q ≤ 1)
XXXXX
Low
(1 < Q <2)
X X
Transmissivity (T, m2/day)15–20< 11–30.5–21–3
Storage coefficient (S)10−1–10−210−2–10−310−3–10−5
Aquifer confinementunconfinedaquitardsemi-confined to confined
More suitable exploitation structuresdug-wells, galleries, and springsXX X
boreholes XXXXX
Direct groundwater recharge (%)25–3020–2510–155–1015–20
Urban infiltration
potential index
(IPI-Urban)
highX
moderateXX XX
low XXXXXX
very low X X

Share and Cite

MDPI and ACS Style

Afonso, M.J.; Freitas, L.; Marques, J.M.; Carreira, P.M.; Pereira, A.J.S.C.; Rocha, F.; I. Chaminé, H. Urban Groundwater Processes and Anthropogenic Interactions (Porto Region, NW Portugal). Water 2020, 12, 2797. https://0-doi-org.brum.beds.ac.uk/10.3390/w12102797

AMA Style

Afonso MJ, Freitas L, Marques JM, Carreira PM, Pereira AJSC, Rocha F, I. Chaminé H. Urban Groundwater Processes and Anthropogenic Interactions (Porto Region, NW Portugal). Water. 2020; 12(10):2797. https://0-doi-org.brum.beds.ac.uk/10.3390/w12102797

Chicago/Turabian Style

Afonso, Maria José, Liliana Freitas, José Manuel Marques, Paula M. Carreira, Alcides J.S.C. Pereira, Fernando Rocha, and Helder I. Chaminé. 2020. "Urban Groundwater Processes and Anthropogenic Interactions (Porto Region, NW Portugal)" Water 12, no. 10: 2797. https://0-doi-org.brum.beds.ac.uk/10.3390/w12102797

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop