Wildﬁre E ﬀ ects on Groundwater Quality from Springs Connected to Small Public Supply Systems in a Peri-Urban Forest Area (Braga Region, NW Portugal)

: Peri-urban areas are territories that combine urban and rural features, being particularly vulnerable to wildﬁre due to the contact between human infrastructures and dense vegetation. Wildﬁres may cause considerable direct and indirect e ﬀ ects on the local water cycle, but the inﬂuence on groundwater quality is still poorly understood. The aim of this study was to characterize the chemistry of several springs connected to small public supply systems in a peri-urban area, following a large wildﬁre that took place in October 2017. Groundwater samples were collected in four springs that emerged within burned forests, while control samples were from one spring located in an unburned area. Sampling took place from October 2017 until September 2018, starting 15 days after the wildﬁre occurrence, to evaluate the inﬂuence of the time after ﬁre and the e ﬀ ect of precipitation events on groundwater composition. Groundwater samples collected in burned areas presented increased content of sulfate, ﬂuoride and nitrogen and variability in pH values. Iron, manganese and chromium contents also increased during the sampling period. Post-ﬁre concentrations of polycyclic aromatic hydrocarbons (PAHs), mainly the carcinogenic ones, increased especially after intense winter and spring rain events, but the levels did not exceed the guideline values for drinking water.


Introduction
Wildfires are one of the main obstacles to the sustainability of forests and related ecosystems. The resulting devastation may cause severe economic and social costs, with the loss of lives and infrastructures, as well as the disturbance of the provision of goods and services, including ecosystem services, together with environmental damages such as the loss of carbon sequestration [1][2][3].
The urbanization of forest areas constitutes a new fire risk scenario. Over the years, throughout the world, the contact areas between human infrastructures and wild vegetation increased (the so-called peri-urban or wildland-urban interfaces) [4][5][6]. In Portugal, these areas result mainly from the Poisoning may occur from continued consumption of water containing high concentrations of copper (Cu), with gastrointestinal symptoms. Arsenic (As) and cromium (Cr) (specially hexavalent Cr) may be carcinogenic, while aluminium (Al), lead (Pb) and mercury (Hg) are toxic when consumed in sufficient quantities for prolonged periods. Following wildfire, increased exports of nitrogen (N) and phosphorous (P), in various forms, can also be problematic for managers of water supply catchments. High concentrations of nitrates (NO 3 − ) and nitrites (NO 2 − ) also present a potential risk to human health, primarily through reduction of NO 3 − to NO 2 − , which may affect oxygen transport in red blood cells, while high concentrations of ammonium NH 3 + /NH 4 + may corrode copper pipes and fittings.
N and P are limiting nutrients for growth of aquatic plants, algae and cyanobacteria in water bodies. Eutrophication increases the risk of potentially toxic blooms, with implications for human health, aesthetic problems (taste, odor and color), and aquatic ecosystem function disturbance [14,17]. Polycyclic aromatic hydrocarbons (PAHs), such as benzo[a]pyrene, are known for their potential teratogenicity, carcinogenic and mutagenic properties, explained by the formation of adducts between the DNA bases and epoxides derived from hydrocarbons after an oxidizing process in the liver [18]. In addition to cancer, long-term exposition can also cause chronic bronchitis, skin problems and allergies. Some compounds are also classified as potential endocrine disruptors because they have estrogenic activity, generally coupled with a high potential for bioaccumulation [15,16,19,20]. The United States Environmental Protection Agency (USEPA) designated several PAHs as priority pollutants that should be regulated due to their high toxicity and adverse effects [21]. PAHs were also designated as priority hazardous substances by the European Commission, in Directive on Environmental Quality Standards (Directive 2008/105/EC) [22].
The effects of fire on surface water are more evident that those on groundwater. Recharge rates, net infiltration and water balance can change after a wildfire, leading to difficulties in maintaining the supply of potable groundwater to populations [23]. Nowadays, despite the widespread access to tap water provided by the Portuguese municipalities, many people still prefer to drink groundwater from public fountains supplied by nearby springs, as it is associated with high purity and pleasant organoleptic characteristics. However, the impact of wildfires on the composition of groundwater from nearby springs is unknown.
The goal of this study was to identify the impact of the wildfire that occurred on October 2017, on the chemical characteristics of water from springs connected to small public supply systems used for human consumption in the peri-urban area of the city of Braga, in NW Portugal. The study included the analyses of major ions and trace elements, namely PAHs, which were carried out during one year after the wildfire.

Hydrogeological Framework
The city of Braga (41 • 32 N; 8 • 25 W) is located in NW Portugal (Figure 1), in the Minho Province. The municipality has a resident population of 181,494 inhabitants (in 2011) [24], representing the seventh largest municipality in Portugal (by population) whit an area of 183.40 km 2 . The city is surrounded by peri-urban areas consisting of agroforestry systems, especially in the higher and steeper terrain (Figures 1 and 2). To the southeast, the urban area is bordered by a mountainous ridge ( Figure 1) with altitude above 500 m at the summits of Santa Marta (562 m), Monte Frio (548 m) and Sameiro (572 m). To the north, this region falls into the Cávado river catchment, while to the south it falls into the Ave river catchment. Braga region is located in the following geomorphological units [27]: Central System (1st level unit) and Entre Douro e Minho Open Valleys and Hills and Atlantic Front of NW Iberia Mountains (2nd level units). The region is also located in the Central-Iberian Zone of the Iberian Massif [28]. The dominant regional geological units are granitic rocks and metasedimentary  The wildfire started at 13:00h on 15 October and was controlled within in a few hours. However, a strong reactivation projected the fire in several directions, with a very significant speed of propagation, and by 17:17h the authorities requested additional means to protect the houses, since the fire had two very long and intense fronts. Only by 10:09h, on 16 October, the wildfire was finally controlled [30].

Water Sampling
Groundwater samples were collected in springs that supply non-treated water used for human consumption, located in near Braga city ( Figure 1): i) Santa Maria Madalena Fountain (S1); ii) Santa Marta de Leão Fountain (S2); iii) Tanque de Dadim Fountain (S3); iv) Depósitos Spring (S4); v) Monte de Dadim Spring (S5). The recharge areas of S1 to S4 springs were affected by the wildfire, while S5 corresponds to a spring located in an unburned area, which was used as control.
All springs are characterized by relatively shallow water circulation. In fact, water is abstracted by means of hand dug galleries which, in this region, are typically associated with springs connected to the upper saturated zone and, sometimes, also to interflow. Moreover, springs S1 and S2 are located at hilltops and, therefore, have very short circulation paths. Finally, The Braga climate has Atlantic features, with a mean value of annual precipitation around 1449 mm at the local climatic station (located at 190 m of altitude). The temporal distribution pattern of precipitation is similar to the one usually observed in NW Portugal: December is the rainiest month (about 220 mm) and July is the driest month (about 22 mm). The mean annual air temperature is around 15 • C, with maximum in July and August (21 • C) and minimum in January (9 • C). The Braga Köppen-Geiger climate classification is Csb, which is dominant in Northwestern Iberia and corresponds to warm temperate (C), with dry and warm summer (sb) [25,26].
Braga region is located in the following geomorphological units [27]: Central System (1st level unit) and Entre Douro e Minho Open Valleys and Hills and Atlantic Front of NW Iberia Mountains (2nd level units). The region is also located in the Central-Iberian Zone of the Iberian Massif [28]. The dominant regional geological units are granitic rocks and metasedimentary rocks; sedimentary cover areas are residual. Therefore, the major processes in the local groundwater cycle take place mostly in fractured circulation media and, to a minor extent, in porous media.
The unsaturated zone depth often reaches more than 20 m at hilltops and less than 10 m at valley bottoms. In addition, the structure of the unsaturated zone encompasses a soil cover (mainly Leptosols, Regosols and Cambisols, with an umbric A horizon) overlying granite or metasedimentary rock.

Wildfire Description
Between 13 and 18 October of 2017, more than 7900 wildfires affected Northwestern Iberia. During this month, 15,440 wildfires were active, 33 of which were of important size, which spread quickly due to the drought and to particular meteorological conditions: strong winds caused by the Ophelia hurricane that swept the western coast of the Iberian Peninsula and unusual air temperatures, above 30 • C. September 2017 was the driest month in 87 years [29]. Around 81% of the territory was under severe drought and 7.4% in extreme drought. In the district of Braga, 23 wildfires were registered on 14 October and 43 wildfires on 15 October. The wildfire of greatest impact on the study region started in Leitões, Guimarães, and quickly reached Braga (particularly the mountain area between Santa Marta and Sameiro, see Figure 2), consuming over 1200 ha of forest [29].
The wildfire started at 13:00h on 15 October and was controlled within in a few hours. However, a strong reactivation projected the fire in several directions, with a very significant speed of propagation, and by 17:17h the authorities requested additional means to protect the houses, since the fire had two very long and intense fronts. Only by 10:09h, on 16 October, the wildfire was finally controlled [30].

Water Sampling
Groundwater samples were collected in springs that supply non-treated water used for human consumption, located in near Braga city ( Figure 1): (i) Santa Maria Madalena Fountain (S1); (ii) Santa Marta de Leão Fountain (S2); (iii) Tanque de Dadim Fountain (S3); (iv) Depósitos Spring (S4); (v) Monte de Dadim Spring (S5). The recharge areas of S1 to S4 springs were affected by the wildfire, while S5 corresponds to a spring located in an unburned area, which was used as control.
All springs are characterized by relatively shallow water circulation. In fact, water is abstracted by means of hand dug galleries which, in this region, are typically associated with springs connected to the upper saturated zone and, sometimes, also to interflow. Moreover, springs S1 and S2 are located at hilltops and, therefore, have very short circulation paths. Finally, S3, S4, and S5 springs are all located at slopes and abstract water circulating in the granite weathering mantle.
The site selection criteria were the existence of permanent flow throughout the year, the location of the recharge areas regarding the burnt region, and the sampling feasibility. These sampling points were also chosen to be preserved as much as possible from other anthropogenic impacts, reducing the risk of water contamination by sources of pollution other than wildfires (Figures 1 and 2 and Table 1). The control spring analytical results were used to provide a framework of reference for contaminant concentrations on burned areas. The discharge flow in all springs was under 0.5 L/s. Samples collected at S1, S2, S3 and S4 were designated as Burnt samples (BR), and the ones collected at S5 as unburnt samples (NB). The sampling plan included five campaigns, carried out during one year, from October 2017 until September 2018, starting 15 days after the wildfire occurrence (before the first post-fire rain event), and included 25 samples.
Samples were collected according to ISO 5667-3:2003(E) Water quality-Sampling-Part 3: guidance on the preservation and handling of water samples, during non-storm conditions, before and after the first overland flow event and aquifer recharge following the fire.
The influence of wildfire on water quality constituents was investigated. The physicochemical analyses of key water quality constituents included the major ions, organic matter and heavy metals, which may be mobilized according to fire intensity, post-fire precipitation, geological and geomorphological conditions and vegetation cover of each site. Compounds that can serve as specific markers, indicative of wildfire, such as the high molecular weight PAHs, were also monitored. Water pH, electrical conductivity (EC) and temperature were measured in situ during sampling.
The temporal distribution of precipitation and climate data (measured at Braga meteorological station) during the study period, provided the knowledge of the influence of rainfall on the input of pollutants into groundwater [31].

Laboratory Analyses
Analyses were performed according to procedures outlined in Standard Methods for the Examination of Water and Wastewater 23 rd edition and in Le Rodier-L'analyse de l'eau 10 e édition. The laboratory has been accredited under ISO/IEC 17025 since 2007. Precision and accuracy were calculated for all analytical methods with values <10%. Uncertainties were also calculated with results varying from 2% to 10%.
Water turbidity was measured in a Hach 2100N Laboratory Turbidity Meter. Electrical conductivity (EC) and pH were determined in a Crison MultiMeter MM 41. Total alkalinity and total hardness were analyzed by titration and reported as milligrams per liter of calcium carbonate (mgL −1 CaCO 3 ). Color, PO 4 2− and total phosphorus, expressed as P ( Atomic Absorption Spectrometer (Varian Inc., Palo Alto, CA, USA) and in a Continuous Segmented Flow Instrument (San-Plus Skalar, Skalar Analytical, Breda, The Netherlands), respectively. PAHs were analyzed by dispersive liquid-liquid microextraction coupled to gas chromatography/mass spectrometry (DLLME-GC/MS) methodology in a Shimadzu GCMS-QP2010 gas chromatograph mass spectrometer equipped with an auto injector AOC5000 (Shimadzu Corporation, Kyoto, Japan), according the procedure described in Borges et al. [32]. Analytical standards were supplied by Sigma-Aldrich (Steinheim, Germany) and Merck (Darmstadt, Germany). The reference standard mixture containing the 15  Methanol, dichloromethane and acetonitrile were organic trace analysis grade SupraSolv and were supplied by Merck (Darmstadt, Germany). Ultrapure water was highly purified by a Milli-Q gradient system (18.2 mΩ/cm) from Millipore (Milford, MA, USA).

Statistical Studies
Data on chemical concentration was analyzed throughout time, location, and throughout time in burned areas. All dependent variables from every analyzed parameter were tested for distribution of the residuals with the Shapiro-Wilk's test. Chemical concentrations were studied using a one-way analysis of variance (ANOVA), if normal distribution of the residuals was confirmed. Welch correction was applied when the homogeneity of variances was not verified. Whenever statistical significances were found, Tukey's test or the Tamhane's test post-hoc tests were applied for mean comparison, depending on variances assumption or not.
If normal distribution of the residuals was not found, parameters analyzed were studied using a Kruskal-Wallis test. Whenever statistical significances were found, Dunn's post-hoc test was applied for median comparison.
All analyses were performed at 5% significance level, using XLSTAT for Windows version 2014.5 (Addinsoft, Paris, France).

Results and Discussion
The effects of wildfires on the catchment's hydrologic responses were noted worldwide, although the specific impacts are unpredictable in terms of both the magnitude of potential effects and the persistence of the influence. So, there is no clear pattern for wildfire effects on water bodies, which may not be affected or, on the other hand, experience fire-related changes that can range from aesthetic concerns (taste or appearance) to potential toxicity or carcinogenicity with prolonged exposure, as well as environmental damages [14,33].
The concentration of major and trace elements monitored in groundwater collected at the peri-urban area of Braga, after the wildfire of October 2017, shown significant differences between samples for several compounds. However, for others, a clear tendency that could be attributed to the wildfire was not observed, or the variation over time was similar to that observed in control samples.
Post-fire analytical results are shown in Table 2 and corresponding statistical analysis on Supplementary Table S1.
The analyses of the results show a slight decrease followed by an increase in pH values in S1 and S2, and an increase in pH values in S3 and S4 during the sampling period, wherein 18 January and April 2018 were grouped, whereas the remaining sampling times were independent of each other (p = 0.006). This is in accordance with studies that described wood ash as being alkaline and rich in carbonates and metal oxides [14]. However, ash chemical composition is highly variable as it reflects the type of vegetation and the part of the plant burned, as well as the soil type and combustion conditions, with different implications on water quality parameters [14,34]. No statistical differences (p > 0.050) were found for turbidity, color and silica. The electrical conductivity, which is a water quality indicator for estimating the amount of mineralization and total dissolved solids, remained also constant during the study period, with no significant variations in BR samples (p > 0.050).
Wildfires often induce quantitative and qualitative changes in soil organic matter, sometimes with significant losses due to the partial or total removal of the litter layer and, possibly, some organics from the upper few centimeters of mineral soil [35]. However, in this study, the total organic carbon concentrations appeared to be unaffected by fire, with mean values of 0.37 ± 0.15 mgL −1 and 0.33 ± 0.17 mgL −1 in BR and NB samples, respectively, and the chemical oxygen demand increased immediately after the fire, reaching a peak during spring, but then slowly declined as vegetation re-established. Nevertheless, these variations could not be attributed to the wildfire, as they occurred in the NB samples too (p > 0.050).   Few compositional changes regarding major ions were observed. The Piper diagram presented in Figure 3 highlights the main changes: all sampling points reveal a minor shift from the first campaign (October 2017, at the end of the dry season) to the third campaign (April 2018, after the intense March precipitation events). The cations content shift results mainly from a decrease in calcium and an increase in sodium, while the anions content shift results from an increase in chloride and sulfate. However, changes in calcium and chloride ions observed in the burnt area are similar to those observed in the unburnt area and, thus, may not be connected to the wildfire. A decrease in bicarbonate is also observed in S1 to S4 points, but not in S5 point. Bicarbonate is originated by water-rock interaction and depends on the nature of the groundwater flowpaths that, in this case, are poorly known. with significant losses due to the partial or total removal of the litter layer and, possibly, some organics from the upper few centimeters of mineral soil [35]. However, in this study, the total organic carbon concentrations appeared to be unaffected by fire, with mean values of 0.37 ± 0.15 mgL −1 and 0.33 ± 0.17 mgL −1 in BR and NB samples, respectively, and the chemical oxygen demand increased immediately after the fire, reaching a peak during spring, but then slowly declined as vegetation reestablished. Nevertheless, these variations could not be attributed to the wildfire, as they occurred in the NB samples too (p > 0.050).
Few compositional changes regarding major ions were observed. The Piper diagram presented in Figure 3 highlights the main changes: all sampling points reveal a minor shift from the first campaign (October 2017, at the end of the dry season) to the third campaign (April 2018, after the intense March precipitation events). The cations content shift results mainly from a decrease in calcium and an increase in sodium, while the anions content shift results from an increase in chloride and sulfate. However, changes in calcium and chloride ions observed in the burnt area are similar to those observed in the unburnt area and, thus, may not be connected to the wildfire. A decrease in bicarbonate is also observed in S1 to S4 points, but not in S5 point. Bicarbonate is originated by waterrock interaction and depends on the nature of the groundwater flowpaths that, in this case, are poorly known. Regarding the sampling date, ANOVA analyses of major ions also reveal that chloride, sodium, potassium, calcium and magnesium concentrations did not vary significantly (p > 0.050). In contrast, Regarding the sampling date, ANOVA analyses of major ions also reveal that chloride, sodium, potassium, calcium and magnesium concentrations did not vary significantly (p > 0.050). In contrast, the increase of sulfate in groundwater from the burnt areas was statistically significant (p = 0.034) and is probably due to the oxidation of sulfur in soil organic matter after the wildfire. Similar results were already observed in an earlier study carried out in Caramulo region, in Central Portugal [36] and reported in the literature by other authors [14].
The levels of fluoride have also risen (p = 0.003) in May by approximately two orders of magnitude regarding the control point (S5), decreasing to levels similar to the initial ones in the last campaign.
Nutrient export from burnt soils usually increases after wildfires, and this process may affect groundwater composition. However, wildfire effects on stream exports of total phosphorus (TP) and total nitrogen (TN) vary significantly [14,37,38]. In the present study, it was observed a small decline (despite this, with no statistical significance) of phosphate, but an increase of TP with multiple change of 1.2 to 4.2 times the initial values. The wildfire effect on nitrogen should be examined with caution, because small differences regarding pre-fire values or control samples could be very important [39,40]. During the first six months after fire, the concentration of nitrogen species increased significantly (p < 0.001), especially after precipitation, compared with concentrations in the reference samples. Combined concentrations of post-fire nitrite, nitrate and ammonium, which are dissolved forms, varied from 0.52 mgL −1 in control to 6.28 mgL −1 (mean values) in S4 sample. The land cover at the S4 recharge area consists only of forest, without other pollution sources besides the wildfire. The higher nitrate content observed in S4 should be a result of the wildfire and is in agreement with the thicker layer of ash observed in the recharge area of this spring. Several factors could explain the increase of nitrogen exports in post-fire situations. On one hand, there is a lower plant demand and the nitrogen mineralization is stimulated (due to changes in pH and electrolytes). On the other hand, the nitrate form is mobile in soil-water systems and leaches through soil into catchments and drainages after heavy overland flow events that follow the wildfire [39,41]. In addition, because the wildfire removes forest cover and litter, rain interception decreases and nutrient transport via infiltration may increase. As well as nitrate, ammonium loading may increase as it is volatilized during fire and can dissolve into water. This compound may be retained in soil in its exchangeable form and subsequently be leached.
Hazardous chemicals, such as metals, were also monitored in Braga, as wildfires may influence the concentration of trace metals differently, with potentially harmful effects to human health and the environment. In this study no significant differences between samples collected in NB and BR locations were observed for cadmium (Cd), arsenic (As), lead (Pb), nickel (Ni), copper (Cu) and zinc (Zn). In contrast, the mean levels of iron (Fe), manganese (Mn) and chromium (Cr) were about 1.3, 7.0 and 2.8 times above the NB sample levels, respectively. Substantial post-fire increases in total iron and total manganese have been reported in the literature, indicating an added influx of these metals as part of an increase in particulates [42]. Even at values below the criteria established for aquatic systems, results must be considered as iron and manganese are related to aesthetic issues of the water (taste and color), and chromium, namely the hexavalent form, is carcinogenic. Concentrations of metals in groundwater samples are displayed in Table 3. Table 3. Descriptive statistics of selected trace elements in groundwater samples after the wildfire.

PAH Analyses
This study also included the analyses of the 15 PAHs designated as priority hazardous pollutants. Thirteen hydrocarbons were found during the sampling period. The most abundant compounds in samples collected in burnt areas were Ant (26%), Acy (17%), BaA (15%) and BaP (14%), while in control samples Acy was the dominant compound (33%), followed by BaA (21%) and Ant (18%) (Figure 4). The sum of total concentrations, as well as the variety of PAHs, has increased throughout the year in all sampling points, with maximum values in May, after the wet season precipitation events. Mean values ranged from negligible amounts in S5 sample in January to 0.029 µgL −1 in S3 in May. Comparing the initial to the last campaign, total values increased by factors of 1.3 to 2.2 times the initial concentrations. Furthermore, in control samples, only light PAHs, with three to four rings (Acy, Ace, Flu, Phe, Ant, Flt, Pyr, BaA, and Chr), were detected, in contrast to the samples of burnt areas, where a different profile was observed throughout the study period. Heavy PAHs, with five (BbF, BkF, BaP, and DahA) and six rings (InP, BghiP), have increased until May followed by a decrease in September, which probably indicates a natural remediation concerning these compounds. The compositional pattern of PAHs by ring size for the water samples is shown in Figure 5.   Significant differences for some PAHs were observed if the date of sampling was considered. Results from Shapiro-Wilk's test shown that PAHs do not have a normal distribution of the residuals. Kruskal-Wallis analysis for the sum of PAHs revealed that October 2017 and May 2018 sample results were significantly different. (Table S1).
Regarding the carcinogenic PAHs (BaA, Chr, BbF, BkF, BaP, Ind and DahA), the concentrations were also significantly higher in May 2018 (p = 0.022), at levels around 0.010 µgL −1 (0.0098 to 0.013 µgL −1 ) in BR samples, representing 60.4%, 52.9%, 36.2% and 62.6% of total PAHs in S1, S2, S3 and S4,    Significant differences for some PAHs were observed if the date of sampling was considered. Results from Shapiro-Wilk's test shown that PAHs do not have a normal distribution of the residuals. Kruskal-Wallis analysis for the sum of PAHs revealed that October 2017 and May 2018 sample results were significantly different. (Table S1).
Regarding the carcinogenic PAHs (BaA, Chr, BbF, BkF, BaP, Ind and DahA), the concentrations were also significantly higher in May 2018 (p = 0.022), at levels around 0.010 µgL −1 (0.0098 to 0.013 µgL −1 ) in BR samples, representing 60.4%, 52.9%, 36.2% and 62.6% of total PAHs in S1, S2, S3 and S4, Significant differences for some PAHs were observed if the date of sampling was considered. Results from Shapiro-Wilk's test shown that PAHs do not have a normal distribution of the residuals.
Kruskal-Wallis analysis for the sum of PAHs revealed that October 2017 and May 2018 sample results were significantly different. (Table S1).
Regarding the carcinogenic PAHs (BaA, Chr, BbF, BkF, BaP, Ind and DahA), the concentrations were also significantly higher in May 2018 (p = 0.022), at levels around 0.010 µgL −1 (0.0098 to 0.013 µgL −1 ) in BR samples, representing 60.4%, 52.9%, 36.2% and 62.6% of total PAHs in S1, S2, S3 and S4, respectively. Approximately half of these values correspond to BaP, with concentrations ranging from 0.0046 µgL −1 in S4 to 0.0077 µgL −1 in S1. In S3, concentrations of BaP increased in April 2018 and remained approximately constant up to September 2018, with a mean value of 0.0052 ± 0.0006 µgL −1 . In S5 samples (NB), BaP was not detected during the study period ( Figure 6). In S5 samples (NB), BaP was not detected during the study period ( Figure 6). Groundwater quality standards for PAHs (referred to as "threshold values") were established by several European Member States taking into account identified risks [43]. Ranges of threshold values throughout Europe varied from 0.005 to 0.03 µgL −1 for BaP, 0.01 to 0.1 for BghiP and BbF, and 0.05 to 0.1 µgL −1 for BkF. Regarding these values, our results reveal several nonconformities for BaP, namely in samples collected in April and May in S1, S2 and S3. BaP, is the most extensively studied carcinogenic PAH, classified by IARC as a Group 1 or a known human carcinogen [44], and is in the top ten priority pollutants designated by the Agency for Toxic Substances and Disease Registry (ATSDR) in 2015 [45]. BaP is also considered a powerful endocrine disruptor compound [46]. Despite the parametric value for drinking water established by the European Union Council Directive 98/83/EC for BaP has not been exceeded (0.010 µgL −1 ) [47], attention is needed especially regarding mixture cancer potency, as individual PAHs occur as part of environmental mixtures and cumulative risk assessment should be considered [48]. The established limit of 0.100 µgL -1 for the sum of concentrations of BbF, BkF, BghiP and Ind was also not exceeded.
Regarding PAHs sources, they may have a pyrogenic origin, linked to processes of incomplete combustion of organic matter (in particular, vegetation or fossil fuels) or a petrogenic origin, resulting from the transformation of organic matter that occurs in geological materials [49]. Another source of PAH is the activity of plants, algae/phytoplankton, and microorganisms (biogenic PAHs) [50,51]. Diagnostic ratios based on PAHs physical and chemical properties and stability against photolysis can be used to distinguish between petrogenic and pyrogenic origins. Congener ratios of (Ant/Ant + Phe) > 0.1 and (Phe/Ant) < 10 indicate pyrogenic sources. The (BaA/BaA + Chr) and (Flt/Flt + Pyr) ratios are used to further distinguish combustion sources: high ratios (>0.35 and >0.50, respectively) indicate grass, wood or coal combustion, intermediate ratios (0.20-0.35 and 0.40-0.50, respectively) indicate liquid fossil fuel combustion or mixed petrogenic and pyrogenic origin and low ratios (<0.20 and <0.40) usually imply petrogenic sources [52,53].
The ratios of (Ant/Ant + Phe) and (Phe/Ant) were calculated with results between 0.68-0.82 and between 0.23-0.50, respectively, for samples in which PAHs were detected. The application of Groundwater quality standards for PAHs (referred to as "threshold values") were established by several European Member States taking into account identified risks [43]. Ranges of threshold values throughout Europe varied from 0.005 to 0.03 µgL −1 for BaP, 0.01 to 0.1 for BghiP and BbF, and 0.05 to 0.1 µgL −1 for BkF. Regarding these values, our results reveal several nonconformities for BaP, namely in samples collected in April and May in S1, S2 and S3. BaP, is the most extensively studied carcinogenic PAH, classified by IARC as a Group 1 or a known human carcinogen [44], and is in the top ten priority pollutants designated by the Agency for Toxic Substances and Disease Registry (ATSDR) in 2015 [45]. BaP is also considered a powerful endocrine disruptor compound [46]. Despite the parametric value for drinking water established by the European Union Council Directive 98/83/EC for BaP has not been exceeded (0.010 µgL −1 ) [47], attention is needed especially regarding mixture cancer potency, as individual PAHs occur as part of environmental mixtures and cumulative risk assessment should be considered [48]. The established limit of 0.100 µgL −1 for the sum of concentrations of BbF, BkF, BghiP and Ind was also not exceeded.
Regarding PAHs sources, they may have a pyrogenic origin, linked to processes of incomplete combustion of organic matter (in particular, vegetation or fossil fuels) or a petrogenic origin, resulting from the transformation of organic matter that occurs in geological materials [49]. Another source of PAH is the activity of plants, algae/phytoplankton, and microorganisms (biogenic PAHs) [50,51]. Diagnostic ratios based on PAHs physical and chemical properties and stability against photolysis can be used to distinguish between petrogenic and pyrogenic origins. Congener ratios of (Ant/Ant + Phe) > 0.1 and (Phe/Ant) < 10 indicate pyrogenic sources. The (BaA/BaA + Chr) and (Flt/Flt + Pyr) ratios are used to further distinguish combustion sources: high ratios (>0.35 and >0.50, respectively) indicate grass, wood or coal combustion, intermediate ratios (0.20-0.35 and 0.40-0.50, respectively) indicate liquid fossil fuel combustion or mixed petrogenic and pyrogenic origin and low ratios (<0.20 and <0.40) usually imply petrogenic sources [52,53].
The ratios of (Ant/Ant + Phe) and (Phe/Ant) were calculated with results between 0.68-0.82 and between 0.23-0.50, respectively, for samples in which PAHs were detected. The application of (BaA/BaA + Chr) and (Flt/Flt+Pyr) ratios was also performed, with values > 0.35 and >0.50, also suggesting a pyrogenic source of PAHs (Figure 7).

Conclusions
This study identified the main impacts of a large forest wildfire on groundwater quality of springs connected to small public supply systems in a peri-urban area. Results pointed out that the best practice for assessing wildfire hydrochemical effects is to start monitoring programs immediately after the wildfire event, and proceed with sampling campaigns for at least 12 months, since several parameters show considerable variation through time. It was also found that extreme events, like intense precipitation, were much more important for groundwater contamination than long-term average changes.
An increase in several parameters, such as sulfate, fluoride, phosphorous, nitrogen compounds, iron, manganese and chromium was observed in Braga peri-urban aquifers. Concerning PAHs, the results reflected the fire impact mainly through the profile of the compounds that appear in the BR springs, which differ significantly from the control spring. Six months after the wildfire, and after the first intense rain event, carcinogenic PAHs, including BaP, began to be detected in considerable concentrations, corroborating the idea of it being difficult to predict the long-term impacts of wildfires on groundwater quality.
Although the connection between groundwater depletion and destructive wildfires might seem tenuous at first glance, and the parametric values for drinking water established by international guidelines have not been exceeded, the results clearly demonstrate the vulnerability of aquifers to wildfires, especially for PAHs, which constitutes an issue yet poorly understood in terms of both the magnitude and persistence.
More information is needed on appropriate monitoring strategies in order to identify a standard set of trace and major compounds to be analyzed and establish protocols to effectively assess water quality, which is essential for developing sustainable water resources management practices.
Supplementary Materials: The following are available online at www.mdpi.com/xxx/s1, Table S1: Values for statically analyses trends in the different sampling time for burned areas.

Conclusions
This study identified the main impacts of a large forest wildfire on groundwater quality of springs connected to small public supply systems in a peri-urban area. Results pointed out that the best practice for assessing wildfire hydrochemical effects is to start monitoring programs immediately after the wildfire event, and proceed with sampling campaigns for at least 12 months, since several parameters show considerable variation through time. It was also found that extreme events, like intense precipitation, were much more important for groundwater contamination than long-term average changes.
An increase in several parameters, such as sulfate, fluoride, phosphorous, nitrogen compounds, iron, manganese and chromium was observed in Braga peri-urban aquifers. Concerning PAHs, the results reflected the fire impact mainly through the profile of the compounds that appear in the BR springs, which differ significantly from the control spring. Six months after the wildfire, and after the first intense rain event, carcinogenic PAHs, including BaP, began to be detected in considerable concentrations, corroborating the idea of it being difficult to predict the long-term impacts of wildfires on groundwater quality.
Although the connection between groundwater depletion and destructive wildfires might seem tenuous at first glance, and the parametric values for drinking water established by international guidelines have not been exceeded, the results clearly demonstrate the vulnerability of aquifers to wildfires, especially for PAHs, which constitutes an issue yet poorly understood in terms of both the magnitude and persistence.
More information is needed on appropriate monitoring strategies in order to identify a standard set of trace and major compounds to be analyzed and establish protocols to effectively assess water quality, which is essential for developing sustainable water resources management practices.