Environmental Research Communications LETTER • OPEN ACCESS Influence of land use changes on submarine groundwater discharge To cite this article: Martí Rufí-Salís et al 2019 Environ. Res. Commun. 1 031005 View the article online for updates and enhancements. This content was downloaded from IP address 158.109.161.66 on 02/05/2019 at 08:12 Environ. Res. Commun. 1 (2019) 031005 https://doi.org/10.1088/2515-7620/ab1695 LETTER Influence of land use changes on submarine groundwater discharge OPEN ACCESS Martí Rufí-Salís1,2 , Jordi Garcia-Orellana1,3 , Gerard Cantero4, Jordi Castillo5, Almudena Hierro6, RECEIVED 14 November 2018 Joan Rieradevall1,2 and Joan Bach7 1 Institut de Ciència i Tecnologia Ambientals (ICTA), Z Building, Universitat Autònoma de Barcelona (UAB), Campus UAB, 08193 REVISED 2 April 2019 Bellaterra, Barcelona, Spain 2 Department of Chemical, Biological and Environmental Engineering, Q Building, Universitat Autònoma de Barcelona (UAB), Campus ACCEPTED FOR PUBLICATION 5 April 2019 UAB, 08193 Bellaterra, Barcelona, Spain 3 Department of Physics, Facultat de Ciències, Universitat Autònoma de Barcelona (UAB), Campus UAB, 08193 Bellaterra, Barcelona, PUBLISHED Spain 30 April 2019 4 Fundació Mar, Llibertat Street 1, 17200 Palafrugell, Girona, Spain 5 Monteditorial, Carme Street 37, 08500 Vic, Barcelona, Spain Original content from this 6 Institut de Técniques Energètiques, Universitat Politècnica de Catalunya, 08028 Barcelona, Spain work may be used under 7 External Geodynamics and Hydrogeology Group, Universitat Autònoma de Barcelona (UAB), Campus UAB, 08193 Bellaterra, the terms of the Creative Commons Attribution 3.0 Barcelona, Spain licence. E-mail: Jordi.Garcia@uab.cat Any further distribution of this work must maintain Keywords: land use, submarine groundwater discharge, groundwater, urbanization, biogeochemical cycles attribution to the author(s) and the title of Supplementary material for this article is available online the work, journal citation and DOI. Abstract The 20th century has been characterized by an exponential population growth, with a high density in coastal zones. The aim of this work is to study the impact of land use changes on the hydrological cycle, and possible consequences in marine environment. The study has focused on the relationship between coastal urbanisation and submarine groundwater discharge (SGD) that may cause an alteration in the biogeochemical cycles of marine ecosystems. An analysis of land use changes, historical salinities and nutrients data in the Maresme (Barcelona) coastal zone have revealed that the percentage of urban and forest zone land use has increased since 1990, corresponding to a decrease of nutrients in coastal waters through SGD. These impacts may cause changes in coastal biogeochemical cycles, like the decrease in chlorophyll concentrations. Results denote the correlation between land use, coastal aquifer dynamics and the effects of SGD in coastal environments. 1. Introduction Urbanization of coastal areas worldwide has rapidly accelerated during the second part of the 20th and the first decade of the 21st centuries, and is one of the main causes of environmental changes (Grimm et al 2008). San Francisco, New York, Shanghai, Cape Town, Rio de Janeiro and Barcelona are remarkable examples of these rapidly growing coastal mega-cities. This urbanization increase has led to a cascade of environmental changes, as for example, by impacting groundwater resources (Foster 2001, Michael et al 2005). Groundwater resources, in terms of aquifer recharge and water geochemical properties, are also affected by rural flight or spread of forest areas. The relationship between land uses changes, characterised by an urban area expansion, and the quality of groundwater resources has been extensively studied (McFarlane 1984, Tellam 1995, Cronin et al 2003, Taniguchi et al 2007). The majority of the studies have reported worsening groundwater quality in response to increased urbanization. There are several anthropogenic impacts that may alter the characteristic complexity of coastal groundwater systems, such as the construction of new neighbourhoods, facilities (e.g. harbours, airports), expansion of roads, railroads, malls, etc. One of the most important characteristics of this urban expansion with respect to coastal aquifer dynamics is the large-scale change in surface permeability, driven by an increase in paved-surfaces. Decreased infiltration can affect both water quality and quantity of aquifers systems, © 2019 The Author(s). Published by IOP Publishing Ltd Environ. Res. Commun. 1 (2019) 031005 assessed through specific parameters: piezometric levels, conductivity or nutrients concentrations (e.g. nitrates or silicates). Submarine Groundwater Discharge (SGD) is an important component of the hydrological cycle in coastal aquifers and is a source of nutrients and micronutrients to coastal areas (Rengarajan and Sarma 2015, Rodellas et al 2015, Trezzi et al 2016). SGD is generally composed of a mixture between freshwater and circulated seawater, the latter of which is either recirculated through intertidal and subtidal sediments or derived from seawater intrusion which has infiltrated into depleted coastal aquifers (Burnett et al 2003). Therefore, any alteration in coastal aquifers due to changes in land use could trigger an impact in SGD flows, and hence its impact on coastal biogeochemical cycles. For example, along the Kona coast of Hawaii, SGD-driven nitrate and silicate loads were found to vary along shoreline segments with respect to population density, the percentage of bare land and the presence or absence of a golf course. On Maui, SGD-driven nitrogen loads are greatest for sugarcane and pineapple farming land use areas (Bishop et al 2017). A Long Island, NY embayment experiences high SGD-driven nitrogen loadings along shoreline segments that are classified as high-development land use (Young et al 2015). Considering these premises, the aim of this work is to study the interrelation between land use changes, the alteration of coastal aquifers dynamics and the consequent alteration of submarine groundwater discharge (SGD) flow and nutrients fluxes to coastal ocean. 2. Materials and methods 2.1. Study area Maresme county has an area of 398 km2 and is located in the western Mediterranean Sea (figure S1(A) is available online at stacks.iop.org/ERC/1/031005/mmedia), in the region of Catalonia (figure S1(B)). The county is situated along the sea, with 50 km of coastline and a width between 5 and 15 km. The climatology of this region is characterized by Western Mediterranean conditions, and the region experiences high temperatures in summer (mean ∼23 °C) and mild temperatures in winter (mean ∼10 °C). The annual precipitation of the area ranges from 500 to 800 mm and is mainly related to intense and short storms during the fall and spring season. The urban extension of Maresme county has significantly increased (Mateu 1985, Benaiges and Corbera 2007) during the second part of the 20th century, in part due to its proximity to the metropolis of Barcelona, with a population of more than 3,000,000 inhabitants. The present study is focused on the central area (CA) of the Maresme region, including the municipalities of Argentona, Cabrera de Mar and Mataró (county capital). This decision was based on the importance of the Argentona ephemeral stream (i.e. Riera d’Argentona) where surface water flows are restricted to the rainy seasons, particularly during major rainfall events (e.g., >50 mm in a few hours) (Cerdà-Domènech et al 2017). Relevant data from the northern area (NA) (with the municipalities of Arenys de Munt, Arenys de Mar and Sant Vicenç de Montalt) and southern area (SA) (with the municipalities of Alella, Teià and El Masnou) were included in this study, when necessary, to help with the discussion. The Maresme region is composed of a series of coastal aquifers, which can be assessed according to: (1) land use, (2) aquifer dynamics and (3) SGD and its influence on the coastal environment. 2.2. Land use temporal quantification Land use was assessed by the digitalization of orthophotos from the Cartographic and Geological Institute of Catalonia, by using software Miramon 7.1. The reference years selected were 1956 (first cartographic data available), 1990, 2000 and 2010 (figures S2 to S5, respectively). Types of land use considered were Urban, Forest and Agriculture. 2.3. Hydrology data acquisition A summary of data used in the present study is shown in table S1. Data sources are described in Content 1 of the Supporting Information. The location of the sampling points is shown in figures S6 to S9. 2.4. Hydrology data assumptions There are several techniques and analytical methods that have been used as tracers of SGD (e.g. salinity, Rn, Ra-isotopes, heat, seepage meters, etc). Dissolved inorganic silicate (DSi) data was selected as the main approximation of SGD because unlike nitrate or phosphate, DSi is not a fundamental macronutrient in the Mediterranean Sea and consequentially is less biologically-sensitive in the coastal marine environment (Hwang et al 2005, Garcia-Solsona et al 2010, Null et al 2014). Chlorophyll data was selected as an approximation for phytoplankton biomass coastal productivity (Eisner et al 2016). Recharge estimation was obtained according to the values stated at table 1. Recharge model used is described in the Supplementary Information. 2 Environ. Res. Commun. 1 (2019) 031005 Figure 1. Land use relative percentages (forest, agriculture and urban) from 1956 to 2010 in the central area (CA) of the Maresme region. Table 1. Parameters used in the recharge model. Parameter Forest Agriculture Urban Field capacity (volumetric 0.25 0.25 — contents) Initial humidity (volumetric 0.25 0.25 — contents) Root thickness (m) 0.5 0.15 — Wilting point (volumetric 0.1 0.1 — contents) Lamination value (mm) 100 100 0 Net reservoir (mm) 75 22.5 0 Initial reservoir (mm) 75 22.5 0 3. Results and discussion 3.1. Land use Land use variation over the last 60 years is shown in figure 1. It can be observed that in 1956 there was a predominant use of land for agricultural purposes (71%) with respect to forests and urban areas (25% and 3%, respectively). In the succeeding years following 1956, the percentage of urban areas increased from 3% (1956) to 30% (2010). This urbanization process was greatest during the period 1990–2000 (from 12% to 29%). During this time, agricultural areas, historically known to be the most important in the Catalonia region and dedicated to growing potatoes (until 1968) and strawberries and ornamental plants (from 1965) (Badosa 2013), were abandoned and as a consequence, the percentage of forest mass increased. Regarding to the spatial distribution of land uses changes, urban areas tend to be located in the flattest zone of the region (closer to the coast), while forests are more likely to be located near the mountain range (Serra de Marina) (figures S2 to S5). 3.2. Aquifers Piezometric levels of the aquifer did not follow a clear temporal trend (figure 2). It is difficult to establish a simple correlation between piezometric levels and a single factor, such as precipitation or land use. This complexity could be caused by the influence of other factors such as the supply of an external water network. Specifically, the local drinking water company (ATLL) has supplied an average of 21.8 hm3 of external water to the Maresme county yearly since 2005. As a consequence, groundwater pumping and withdrawal in the region is minor, thereby recharging the aquifer. In addition, another factor that could decrease water pumping is the drop of agricultural activities in the analysed time period, reducing the overall water needed for crop irrigation. In terms of conductivity, different values and trends are observed, depending on the distance of any individual well to coastline (figure 3). Conductivity values in piezometers closest to the coastline were higher (up to 2200 μS·cm−1 in some cases) than inland piezometers as a result of coastal aquifer-seawater interactions (figures 3 and S8). Nitrate concentrations generally decreased in the majority of wells analysed over time. However, a notable peak was observed during 2006–2008, corresponding to a drought in Catalonia during those three years (March and Saurí 2013) (figure 2). 3 Environ. Res. Commun. 1 (2019) 031005 Figure 2. Evolution of nitrate and piezometric levels in inland wells from different areas of the Maresme county (Central Area (AC), Northern Area (AN), and Southern Area (AS)), and the evolution of aquifer recharge and rainfall. Piezometric units (m) refer to meters above sea level. Figure 3. Evolution of conductivity at 20 °C in inland wells from different areas: Central Area (AC) and Northern Area (AN) of the Maresme region. Data in parenthesis refer to distance to coastline. 3.3. Seawater quality In the central part of the studied area (CA), an increase of salinity in coastal waters can be observed, from ∼34 to 38 psu from 1994 year to 2014 year (figure S10). Silicate concentrations decreased over the same time period (figure S11), from 0.31 mg·L−1 (1995) to 0.08 mg·L−1 (2012), while chlorophyll concentrations (figure S12) 4 Environ. Res. Commun. 1 (2019) 031005 Figure 4. Evolution of nitrate, silicate and chlorophyll concentrations in seawater and the temporal change in land use. decreased from 3.86 μg·L−1 (1994) to 1.33 μg·L−1 (2012). Both trends are observed throughout the study period in the seawater adjacent to the Argentona ephemeral stream (CA). In other locations of the Maresme region, silicate concentrations and salinity were nearly constant over the studied period (figures S11 and S12). These parameters denote that there were not visible changes in the SGD flow in these zones during the period under study. 3.4. Interrelation between land use, aquifers and submarine groundwater discharge Threshold precipitation (necessary precipitation for runoff) decreased with increasing urban land use (378 ha·y−1 from 1990 to 2000), leading to a ground waterproofing process. This urbanization process triggered an increase of urban water (i.e. rainwater fallen on urban land) during the 1956–2010 period, from 5 to 49 hm3. Concurrently, agricultural land also decreased (−523 ha·y−1 from 1990 to 2000) (figure 1). Additional factors influenced the decrease in recharge, including the annual variation in rainfall and water uptake and evapotranspiration of the forest mass due to its increase in surface area (50 to 70 hm3 in the 1956–2010 period). The drop of agriculture activities (from 3.5% to 1.6% in GDP contribution of the primary sector in the 1991–2006 period (Tirado 2007)), related to the reduction of groundwater pumping from aquifers, could also influence groundwater recharge due to a pressure decrease from groundwater extraction. Nitrates concentrations decreased in most of the wells closest to the shoreline (figure 2), (e.g. well 08009–0060 from 118.8 mg·L−1 in 1996 to 77.1 mg·L−1 in 2015) where agricultural and urban land uses have historically been located (figures S2 and S3). Despite the influence of the 2006–2008 drought on both nitrates and conductivity, the latest was also found to decrease overtime. The decrease in nitrate concentration and conductivity could be linked to the decrease in agricultural land use (20,750 to 6,895 ha in the 1956–2010 period), related to the extensive use of fertilizers and water demand for agricultural purposes, respectively. There was a strong relationship between the decrease in precipitation and piezometric levels, and the increase in nitrate concentrations, during the general drought of 2006–2008 (yellow shade in figure 2). Moreover, the influence of the sea is clearly observed on the groundwater conductivity (figure 4). Higher average conductivity values (up to 2000 μS·cm−1) and wider fluctuations denote a higher sensitivity of wells closest to the sea. In specific times, there were piezometric levels (in sample points close to the coastline) that were below sea level (figure 2); such difference could drive saltwater intrusion into the coastal aquifer, as indicated by higher conductivity values in nearby wells. Focusing on seawater, figure 4 shows weak correlations between nutrient concentrations (nitrate and silicate) and biotic activity (chlorophyll concentration), such as silicates (1.4 mg·L−1) and chlorophyll (28.3 μg·L−1) in 1995. However, other factors like ocean currents (Raimbault et al 1993) or nutrient uptake by biological activity may affect these correlations, as in 1998 (peaks of nitrate and silicate). The largest variations in both nutrients and chlorophyll (3.5 μg·L−1 on average) were located within the period 1990–2000, where land use showed the largest variation. Land use variation got stabilised in 2000, corresponding to fairly constant nutrient (e.g. silicates 0.10 mg·L−1) and chlorophyll (1.5 μg·L−1) concentrations. This suggests that there was a relationship between land use changes and the alteration of the biochemical cycles in the coastal sea. These observations are similar to those showed by other studies where Ra-isotopes were used as a SGD tracer (Tang et al 2015). 5 Environ. Res. Commun. 1 (2019) 031005 Figure 5. Correlation between salinity and silicate concentrations in different sample points and different periods in the Maresme county. Circle indicates low silicate concentrations and high salinity values. During the studied period of time, there was an inverse correlation between silicate concentration and salinity in the Argentona ephemeral stream (Mataró) before 2000, for example in 1994 (0.27 mg·L−1; 35.4 psu) and in 1998 (0.19 mg·L−1; 36.5 psu) (figure 5). This inverse correlation shows the existence of important SGD flows (considering silicate as SGD tracer). Furthermore, as seen in other sampling sites such as Arenys de Mar or El Masnou, also in the Maresme county, the correlation between silicate concentration and salinity was not as significant as in the Argentona ephemeral stream, suggesting that SGD was less important in these areas. In addition, the silicate-salinity correlation in the Argentona ephemeral stream faded with time, reaching a salinity of ∼38 psu with negligible silicate by year 2000. Therefore, this drastic fading of the relationship between salinity and silicates, the stabilization of salinity after 2000 and the no observation of any significant peak of silicates, nitrates and chlorophyll would indicate a decrease in the constant influence of the SGD in the area rather than a decrease in the influence of the ephemeral stream. Therefore, land impermeability due to the construction of new infrastructures such as malls, motorways or airports and changes in land use from agriculture to urban or forest areas around large coastal cities, such as Barcelona, may alter land-ocean interactions and reduce, for example, allochthonous nutrient fluxes to coastal water triggered by SGD. These changes and the pollution of water and sediments due to anthropogenic activities may considerably alter coastal ecosystems, reducing primary production and biodiversity, especially in semi- arid regions as the oligotrophic Mediterranean Sea, where runoff is limited. 4. Conclusions There was a clear relationship between urbanization over areas occupied by unconfined coastal aquifers and SGD in Maresme (Barcelona). The exponential population growth of the last century has created a high population-density in coastal areas. Urbanization has therefore triggered changes in the hydrological cycle by altering groundwater recharge through the paving of large surface coastal areas, the abandonment of agricultural areas and an increase of forest cover. In the present example, it was not possible to establish quantitative relationships in terms of historical piezometric levels, salinity and nutrient concentrations in groundwater (DIN) and coastal seawater (DSi) with respect to land use changes; however, we have established a qualitative pattern. There was no clear trend in piezometric levels and land use, according to the location of the piezometers; this may be partially driven by the abandonment of agricultural areas, improvements in agricultural irrigation or external water supply from a drinking water company. Moreover, the increase in forest cover could have caused a decrease in the water table in some specific areas. Nitrate concentrations decreased over 2000–2010 time interval, probably due to the abandonment of agricultural areas in the flatter part of the region and improvements of agriculural techniques. These combined alterations to the hydrological cycle and groundwater quality altered submarine groundwater discharge (SGD) and its associated nutrient flux to the coastal ocean. Analysis of historical salinity and nutrient concentrations in coastal waters showed a steady reduction in the central area (CA) of the study area. Land use driven changes to groundwater quantity and quality may have caused changes in some biogeochemical cycles of the marine environment through SGD, as qualitatively observed by decreasing chlorophyll concentrations in coastal surface waters (from 7.7 μg·L−1 in 6 Environ. Res. Commun. 1 (2019) 031005 1995 to 1.3 μg·L−1 in 2012). Further research is needed to extrapolate the results found in the present study to worldwide locations where the process of urbanization affects coastal zones. This aspect is of great importance in terms of groundwater resources management and for biotic activity in coastal seawater. Acknowledgments This work was partially funded by the project CGL2016-77122-C2-1-R/2-R of the Spanish Government. The authors are grateful to the Universitat Autònoma de Barcelona for awarding a research scholarship (PIF-UAB 2017) to M Rufí-Salís. The authors would like to thank Agència Catalana de l’Aigua (ACA) for the data provided for this research. Jordi Garcia-Orellana wants to thank the support of the Generalitat de Catalunya to MERS (2018 SGR-1588) and J Tamborski. This work is contributing to the ICTA ‘Unit of Excellence’ through the ‘María de Maeztu’ program (MinECo, MDM2015-0552). All data used in this study is available upon request to the correspoding author. ORCID iDs Martí Rufí-Salís https://orcid.org/0000-0003-3696-1033 Jordi Garcia-Orellana https://orcid.org/0000-0002-0543-2641 References Badosa R 2013 Evolució de l’agricultura de la comarca del Maresme en els darrers cent anys 16 201–9 Online:http://revistes.iec.cat/index. php/DA/ Benaiges N and Corbera J 2007 Imatges d’una transformació. El Maresme 60 anys després L’Atzavara 15 133–46 Bishop J M, Glenn C R, Amato D W and Dulai H 2017 Effect of land use and groundwater flow path on submarine groundwater discharge nutrient flux J. Hydrol. Reg. Stud. 11 194–218 Online:https://sciencedirect.com/science/article/pii/S2214581815001135 Burnett W C, Bokuniewicz H, Huettel M, Moore W S and Taniguchi M 2003 Groundwater and pore water inputs to the coastal zone Biogeochemistry 66 3–33 Online:http://link.springer.com/10.1023/B:BIOG.0000006066.21240.53 Cerdà-Domènech M, Rodellas V, Folch A and Garcia-Orellana J 2017 Constraining the temporal variations of Ra isotopes and Rn in the groundwater end-member: implications for derived SGD estimates Sci. Total Environ. 595 849–57 Online:http://sciencedirect.com/ science/article/pii/S0048969717305107 Cronin A A, Taylor R G, Powell K L, Barrett M H, Trowsdale S A and Lerner D N 2003 Temporal variations in the depth-specific hydrochemistry and sewage-related microbiology of an urban sandstone aquifer, Nottingham, United Kingdom Hydrogeol. J. 11 205–16 Online:http://link.springer.com/10.1007/s10040-002-0246-0 Eisner L B, Gann J C, Ladd C, D Cieciel K and Mordy C W 2016 Late summer/early fall phytoplankton biomass (chlorophyll a) in the eastern Bering Sea: spatial and temporal variations and factors affecting chlorophyll a concentrations Deep Sea Res. Part II Top. Stud. Oceanogr. 134 100–14 Online:http://linkinghub.elsevier.com/retrieve/pii/S0967064515002386 Foster S S 2001 The interdependence of groundwater and urbanisation in rapidly developing cities Urban Water 3 185–92 Online:http:// linkinghub.elsevier.com/retrieve/pii/S1462075801000437 Garcia-Solsona E, Garcia-Orellana J, Masqué P, Garcés E, Radakovitch O, Mayer A, Estradé S and Basterretxea G 2010 An assessment of karstic submarine groundwater and associated nutrient discharge to a Mediterranean coastal area (Balearic Islands, Spain) using radium isotopes Biogeochemistry 97 211–29 Online:http://link.springer.com/10.1007/s10533-009-9368-y Grimm N B, Faeth S H, Golubiewski N E, Redman C L, Wu J, Bai X and Briggs J M 2008 Global change and the ecology of cities Science (80-.). 319 756–60 Online:http://science.sciencemag.org/content/319/5864/756 Hwang D-W, Kim G, Lee Y-W and Yang H-S 2005 Estimating submarine inputs of groundwater and nutrients to a coastal bay using radium isotopes Mar. Chem. 96 61–71 Online:http://sciencedirect.com/science/article/pii/S0304420304003123 March H and Saurí D 2013 La sequera del 2007-2008 a la ciutat de Barcelona: gènesi, gestió i visions discordants1 Treballs la Soc. Catalana Geogr. 76 289–96 Online:http://academia.edu/download/34961689/TSCG_MArch_Sauri2013.pdf Mateu J 1985 17–45 Quaderns agraris. (ICEA) Online:http://raco.cat/index.php/QuadernsAgraris/article/view/240661 McFarlane D 1984 Effects of urbanization on groundwater quantity and quality in Perth, Western Australia (https://doi.org/10.4225/23/ 5ab48fe537f9e) Michael H A, Mulligan A E and Harvey C F 2005 Seasonal oscillations in water exchange between aquifers and the coastal ocean Nature 436 1145–8 Online:http://nature.com/articles/nature03935 Null K A, Knee K L, Crook E D, de Sieyes N R, Rebolledo-Vieyra M, Hernández-Terrones L and Paytan A 2014 Composition and fluxes of submarine groundwater along the Caribbean coast of the Yucatan Peninsula Cont. Shelf Res. 77 38–50 Online:http://sciencedirect. com/science/article/pii/S0278434314000247 Raimbault P, Coste B, Boulhadid M and Boudjellal B 1993 Origin of high phytoplankton concentration in deep chlorophyll maximum (DCM) in a frontal region of the Southwestern Mediterranean Sea (algerian current) Deep Sea Res. Part I Oceanogr. Res. Pap. 40 791–804 Online:http://linkinghub.elsevier.com/retrieve/pii/096706379390072B Rengarajan R and Sarma V V S S 2015 Submarine groundwater discharge and nutrient addition to the coastal zone of the Godavari estuary Mar. Chem. 172 57–69 Online:http://sciencedirect.com/science/article/pii/S0304420315000626 Rodellas V, Garcia-Orellana J, Masqué P and Font-Muñoz J S 2015 The influence of sediment sources on radium-derived estimates of submarine groundwater discharge Mar. Chem. 171 107–17 Online:https://sciencedirect.com/science/article/pii/ S0304420315000389 Tang G, Yi L, Liu L and Cheng X 2015 Factors influencing the distribution of 223Ra and 224Ra in the coastal waters off Tanggu and Qikou in Bohai bay Cont. Shelf Res. 109 177–87 Online:http://sciencedirect.com/science/article/pii/S027843431530056X 7 Environ. Res. Commun. 1 (2019) 031005 Taniguchi M, Uemura T and Jago-on K 2007 Combined effects of urbanization and global warming on subsurface temperature in four Asian Cities Vadose Zo. J 6 591 Online:https://soils.org/publications/vzj/abstracts/6/3/591 Tellam J 1995 Urban groundwater pollution in the Birmingham Triassic Sandstone aquifer Proc 4th Annu. IBC Conf. Groundw. Pollution Online:https://scholar.google.es/scholar?q=Urban+groundwater+pollution+in+the+Birmingham+Triassic+sandstone +aquifer&btnG=&hl=ca&as_sdt=0%2C5 Tirado D 2007 Estudi Socioeconòmic de la comarca del Maresme 2007 Online:http://cambradigital.com/c/document_library/get_file? uuid=d8d98fa2-9797-483a-af9e-3b904db55847&groupId=1533402 Trezzi G, Garcia-Orellana J, Rodellas V, Santos-Echeandia J, Tovar-Sánchez A, Garcia-Solsona E and Masqué P 2016 Submarine groundwater discharge: a significant source of dissolved trace metals to the North Western Mediterranean Sea Mar. Chem. 186 90–100 Online:https://sciencedirect.com/science/article/pii/S030442031630113X Young C, Tamborski J and Bokuniewicz H 2015 Embayment scale assessment of submarine groundwater discharge nutrient loading and associated land use Estuar. Coast. Shelf Sci. 158 20–30 Online:https://sciencedirect.com/science/article/pii/S0272771415000530#! 8