Groundwater Natural Contamination by Toluene in Beja and Faro Districts, Portugal

The focus area is located along Beja and Faro districts (southwestern of Portugal), where the Baixo Alentejo Flysch Group exists, is composed of the Mertola, Mira, and Brejeira Formations. The aim of this article is to characterize the groundwater contamination by toluene, one of the petroleum products that may originate from natural petroleum systems and have been introduced into the environment by natural seepage. For this study, twenty soil samples and thirty-one water samples were collected, in 2016, from artesian wells, boreholes, and springs, and analyzed in a GC-FID system. Sampling was performed using headspace vials, with distilled water and a bactericide to inhibit any bacterial activities. The results showed that 93.5% of water samples range between 1000 and 6000 μg/L, and 55% are higher than 3000 μg/L of toluene concentration. In soil samples, amounts of toluene are less than 1000 μg/L in general. These levels are much higher than the World Health Organization health-based guideline for drinkable water (700 μg/L) and higher than those usually found in other places, being only comparable with values from severe industrial contamination. Contaminated groundwater represents a serious concern for human health, and residents should be informed to avoid drinking contaminated ground water or using it for other domestic proposes.


Introduction
The classic Medical Geology approach considers the inorganic elements and substances with more relevance because of their positive or negative effects on human and animal health. However, the natural organic substances must be considered as responsible for many troubles that affect living beings as well [1,2]. Among the 160 EPA (United States Environmental Protection Agency) priority pollutant volatile and semi-volatile organic compounds, only 20 are petroleum-type hydrocarbons [3]. Bandeira de Mello & Miller [4] emphasized the importance of natural organic compounds, for the Medical Geology. Many of these substances, such as toluene, are present on the "ATSDR's Substance Priority List [5]" as the most significant potential threat to human health due to their known or suspicious toxicity and potential for human exposure.
The injection of petroleum hydrocarbons into groundwater system is the object of study for many different branches of science such as biology, hydrogeology, environmental geology, environmental Seepage along faults can be more an effusive process, rather than microseepage [32]. However, pressure gradients cause hydrocarbon fluids to migrate to areas of lower pressure. These migrations may be in directions other than vertical. Natural overburden pressure is a mechanism that forces gas to the surface [33].
Generally, according to Etiope [31], the identification and characterization of hydrocarbon seeps are important for the main reasons: -Seeps can be indicators of natural petroleum systems; -Seeps indicate the occurrence of a fault. According to Etiope [26], gas, in fact, follows preferential pathways of movement (less resistance paths), as determined by fractures and faults; -Seeps can represent a geo-hazard for societal community and industry; -Seeps are natural sources of greenhouse gas.

Objectives
This paper aims to detect and characterize the presence of hydrocarbon fluids, mainly toluene, in soil and groundwater levels of Beja and Faro Districts (Portugal). Following the health-based guideline from World Health Organization, it aims to alert the Portuguese authorities about the potential consequences for the local people. It also aims to discuss the possible sources of these Adapted from [27].
The rate and volume of hydrocarbon seepage to near-surface sediments controls the surface expression of geological, geophysical, biological, and geochemical features, which can range from macro to microseepage [28]. According to Etiope [26], the surface expressions of natural gas seepage can be classified on the basis of spatial dimension, visibility, and fluid typology.
Macroseeps (or seeps) are "channeled" flows of gas, typically related to fault systems. Gas flux is expressed in terms of mass/time (e.g., kg/day or t/year). This surface macroseeps (visible gas vents or oil leaks from the soil or rock outcrops) are generally an indication of a fault in an active petroleum seepage system [29] belonging to a total petroleum system [30,31].
Microseepage is the slow, continual, or episodic loss of light hydrocarbons from gas-oil-prone sedimentary basins where thermal degradation of indigenous organic matter occurs [32]. It is essentially the pervasive, diffuse exhalation of gas from soil, independent of the presence of macroseeps [26].
Seepage along faults can be more an effusive process, rather than microseepage [32]. However, pressure gradients cause hydrocarbon fluids to migrate to areas of lower pressure. These migrations may be in directions other than vertical. Natural overburden pressure is a mechanism that forces gas to the surface [33].
Generally, according to Etiope [31], the identification and characterization of hydrocarbon seeps are important for the main reasons: -Seeps can be indicators of natural petroleum systems; -Seeps indicate the occurrence of a fault. According to Etiope [26], gas, in fact, follows preferential pathways of movement (less resistance paths), as determined by fractures and faults; -Seeps can represent a geo-hazard for societal community and industry; -Seeps are natural sources of greenhouse gas.

Objectives
This paper aims to detect and characterize the presence of hydrocarbon fluids, mainly toluene, in soil and groundwater levels of Beja and Faro Districts (Portugal). Following the health-based guideline from World Health Organization, it aims to alert the Portuguese authorities about the potential consequences for the local people. It also aims to discuss the possible sources of these organic compound seeps and compare the results with other available information from other countries.

Study Area
The focus area is situated in a wide region along Beja and Faro districts (south of Portugal- Figure 2), which include the Baixo Alentejo Flysch Group (BAFG), composed of Mira and Brejeira Formations [34,35]. The SPZ (South Portuguese Zone) is located in the southernmost sector of the pre-Mesozoic Iberian Variscan Chain, which is composed by upper Paleozoic shales and greywackes (interbedded), of low grade metamorphism.
Geosciences 2018, 8, 9 4 of 18 organic compound seeps and compare the results with other available information from other countries.

Study Area
The focus area is situated in a wide region along Beja and Faro districts (south of Portugal- Figure 2), which include the Baixo Alentejo Flysch Group (BAFG), composed of Mira and Brejeira Formations [34,35]. The SPZ (South Portuguese Zone) is located in the southernmost sector of the pre-Mesozoic Iberian Variscan Chain, which is composed by upper Paleozoic shales and greywackes (interbedded), of low grade metamorphism. All these formations contain significant thicknesses of black shales (especially the southernmost formation, Brejeira). The SPZ has several wide areas with significant average total organic carbon content (TOC) and thermal maturation values within hydrocarbons generation zones [36][37][38][39]. The BAFG includes also one more formation (fm.), Mértola. However, because of the highest thermal maturation, conjugated with lowest TOC content among SPZ formations, Mértola Fm. was not included in the surface geochemical sampling [39]. All these formations contain significant thicknesses of black shales (especially the southernmost formation, Brejeira). The SPZ has several wide areas with significant average total organic carbon content (TOC) and thermal maturation values within hydrocarbons generation zones [36][37][38][39]. The BAFG includes also one more formation (fm.), Mértola. However, because of the highest thermal maturation, conjugated with lowest TOC content among SPZ formations, Mértola Fm. was not included in the surface geochemical sampling [39].

Europe
From the SPZ hydrogeological point of view, the water is high mineralized, hard, with high concentration of chlorides, sodium, etc. These characteristics are unexpected from water flowing through low reactive materials like shales and greywacke. The causes of these chemistry aspects are still unclear. The spatial distribution of the main species concentrations is not homogeneous, verifying a clear trend for high concentrations at east, close to Guadiana Basin. The water from this lithologic group is not recommended for human consumption, being almost all MRL (maximum recommended value) overpassed, except for hardness and pH [40].

Surface Geochemical Prospecting
The procedures for surface geochemical sampling campaign begin with preliminary planning actions, such as area selection (target), geological base maps, satellite imagery, etc., as well as choosing the most suitable months for field operations [41]. Once the research area has been established, the efforts are focused on potential hydrocarbon systems (conventional or unconventional), thus defining the positioning of the sampling points. After these adjustments, the most appropriate sampling methodologies and respective standard procedures are defined ( Figure 3).
Geosciences 2018, 8, 9 5 of 18 From the SPZ hydrogeological point of view, the water is high mineralized, hard, with high concentration of chlorides, sodium, etc. These characteristics are unexpected from water flowing through low reactive materials like shales and greywacke. The causes of these chemistry aspects are still unclear. The spatial distribution of the main species concentrations is not homogeneous, verifying a clear trend for high concentrations at east, close to Guadiana Basin. The water from this lithologic group is not recommended for human consumption, being almost all MRL (maximum recommended value) overpassed, except for hardness and pH [40].

Surface Geochemical Prospecting
The procedures for surface geochemical sampling campaign begin with preliminary planning actions, such as area selection (target), geological base maps, satellite imagery, etc., as well as choosing the most suitable months for field operations [41]. Once the research area has been established, the efforts are focused on potential hydrocarbon systems (conventional or unconventional), thus defining the positioning of the sampling points. After these adjustments, the most appropriate sampling methodologies and respective standard procedures are defined ( Figure 3). According to Abrams [28], the light hydrocarbons require special handling and containers to limit volatile loss and prevent post-sampling microbial alteration. Several methods are used for direct hydrocarbon sampling such as probe and headspace for free gases, blender and disrupter for occluded gases, and acid extraction for adsorbed gases and adsorbent fibers [41,42].
In general, the research methodologies in surface geochemistry aim to quantify hydrocarbon movement in a medium, on a time scale [43]. Under standard conditions for temperature and pressure (25 °C and 1 atm), light hydrocarbons (methane [CH4], ethane [C2H6], propane [C3H8], isobutane [i-C4H10], and n-butane [n-C4H10]) are present, on the surface, as a gas. Hydrocarbons of slightly larger molecular weights such as n-pentane (n-C5) and n-hexane (n-C6) are present, on the surface, in the liquid phase. However, they can be found as gas (vapor) just with a small increase of the temperature.
The sampled volume goes on to (1) review of research investigations and case history studies of every gas that might act as an indicator or pathfinder; and (2) elaborate on the models that explain the observed results [44].
The study sampling zone was selected based on geological features such as faults, folds, and thrusts that could work as pathways for the hydrocarbon migration to groundwater ( Figure 1). For this study, the collection of all samples was performed using IsoJar ® (Isotech Laboratories Inc., Champaign, IL, USA) ( Figure 4B)  According to Abrams [28], the light hydrocarbons require special handling and containers to limit volatile loss and prevent post-sampling microbial alteration. Several methods are used for direct hydrocarbon sampling such as probe and headspace for free gases, blender and disrupter for occluded gases, and acid extraction for adsorbed gases and adsorbent fibers [41,42].
In general, the research methodologies in surface geochemistry aim to quantify hydrocarbon movement in a medium, on a time scale [43]. Under standard conditions for temperature and pressure (25 • C and 1 atm), light hydrocarbons (methane [CH 4 ], ethane [C 2 H 6 ], propane [C 3 H 8 ], isobutane [i-C 4 H 10 ], and n-butane [n-C 4 H 10 ]) are present, on the surface, as a gas. Hydrocarbons of slightly larger molecular weights such as n-pentane (n-C5) and n-hexane (n-C6) are present, on the surface, in the liquid phase. However, they can be found as gas (vapor) just with a small increase of the temperature.
The sampled volume goes on to (1) review of research investigations and case history studies of every gas that might act as an indicator or pathfinder; and (2) elaborate on the models that explain the observed results [44].
The study sampling zone was selected based on geological features such as faults, folds, and thrusts that could work as pathways for the hydrocarbon migration to groundwater (Figure 1). For this study, the collection of all samples was performed using IsoJar ® (Isotech Laboratories Inc., Champaign, IL, USA) ( Figure 4B), with distilled water (just for soil/rock samples; water samples do not need to be filled with distilled water) and a bactericide added to inhibit any bacterial activities post sampling. Twenty soil samples were collected for this study, using a drilling machine ( Figure 4D), a metal tube, and a hammer ( Figure 4A,C). Thirty-one water samples were collected from artesian wells, boreholes, and springs ( Figure 2).
Geosciences 2018, 8,9 6 of 18 not need to be filled with distilled water) and a bactericide added to inhibit any bacterial activities post sampling. Twenty soil samples were collected for this study, using a drilling machine ( Figure  4D), a metal tube, and a hammer ( Figure 4A,C). Thirty-one water samples were collected from artesian wells, boreholes, and springs ( Figure 2).

Analytical Methodology
For the analysis of the samples, a gas chromatographer (Agilent 7890B GC; Agilent Technologies, Santa Clara, CA, USA) equipped with a 30 m 0.32 mm Gaspro column and a 50 m 0.32 mm Porabond Q column (both from Agilent) was used. The Gaspro column is coupled to a flame ionization detector and used for analysis of hydrocarbons. The samples were injected into a heated injection valve on the GC using a gastight syringe. The detector responses were calibrated against a reference gas (Scott Specialty Gases Netherlands BV, Breda, The Netherlands) several times a day. Some samples need to be modified for the specific requirement of a particular analytical technique [45]. Most consumer products and geological samples (as hydrocarbons) are composed of a wide variety of compounds that differ in molecular weight, polarity, and volatility [46].
A gas extract is ideally suited for analysis by GC, and this combination is called ''headspace-gas chromatography'' (HS-GC sampling), the fastest and cleanest method for analyzing volatile organic compounds. A headspace sample is normally prepared in a vial containing the sample, the dilution solvent, a matrix modifier like a bactericide, and the headspace ( Figure 5) [45,46].

Analytical Methodology
For the analysis of the samples, a gas chromatographer (Agilent 7890B GC; Agilent Technologies, Santa Clara, CA, USA) equipped with a 30 m 0.32 mm Gaspro column and a 50 m 0.32 mm Porabond Q column (both from Agilent) was used. The Gaspro column is coupled to a flame ionization detector and used for analysis of hydrocarbons. The samples were injected into a heated injection valve on the GC using a gastight syringe. The detector responses were calibrated against a reference gas (Scott Specialty Gases Netherlands BV, Breda, The Netherlands) several times a day. Some samples need to be modified for the specific requirement of a particular analytical technique [45]. Most consumer products and geological samples (as hydrocarbons) are composed of a wide variety of compounds that differ in molecular weight, polarity, and volatility [46].
A gas extract is ideally suited for analysis by GC, and this combination is called "headspace-gas chromatography" (HS-GC sampling), the fastest and cleanest method for analyzing volatile organic compounds. A headspace sample is normally prepared in a vial containing the sample, the dilution solvent, a matrix modifier like a bactericide, and the headspace ( Figure 5) [45,46].
Volatile components from complex sample mixtures can be extracted from non-volatile sample components and isolated in the headspace or vapor portion of a sample vial. An aliquot of the vapor in the headspace is delivered to a GC system for separation of all of the volatile components [46].
In order to achieve the best performance when using HS-GC, careful attention should be used in sample preparation and instrument setup. Key issues to address when setting up HS-GC systems include minimizing system dead volume, maintaining inert sample flow paths, and achieving efficient sample transfer [46]. products and geological samples (as hydrocarbons) are composed of a wide variety of compounds that differ in molecular weight, polarity, and volatility [46].
A gas extract is ideally suited for analysis by GC, and this combination is called ''headspace-gas chromatography'' (HS-GC sampling), the fastest and cleanest method for analyzing volatile organic compounds. A headspace sample is normally prepared in a vial containing the sample, the dilution solvent, a matrix modifier like a bactericide, and the headspace ( Figure 5) [45,46]. Figure 5. Illustrative scheme of a headspace vial Isojar ® , filled with soil/water sample, distilled water, and bactericide. Water samples do not need to be filled with distilled water. Adapted from [42].
distilled water + bactericide soil/water sample Figure 5. Illustrative scheme of a headspace vial Isojar ® , filled with soil/water sample, distilled water, and bactericide. Water samples do not need to be filled with distilled water. Adapted from [42]. The adopted procedures for taking cuttings samples in IsoJar ® , following Isotech [42] recommendations, are presented: 1.
Put the sample in jar (approximately 1 cup); 2.
Only add distilled water to reach the line on the label; 3.
Do not fill the jar with water as the headspace gap is needed to allow gas to desorb into the gap; 4.
Leave about a 1-inch gap between the water and top of the jar; 5.
Add 10 drops of the dilute bactericide, Benzalkonium Chloride, to the jar; 6.
Screw the lid on as tight as possible; 7.
Tape the lid to keep it tight and from vibrating loose during shipment. Tape in same direction (clockwise) that the lid is screwed on.
According to Abrams [28], the fifth point is used to prevent in-situ microbial activity. The amount of bactericide required to prevent microbial activity has not been rigorously examined. To the best of the author's knowledge, no comprehensive studies have been conducted to evaluate the effectiveness of these anti-bacterial agents, how much is required, or if the anti-bacterial agents impact the geochemistry results.
In order to discard the possibility of the expressive amounts of any hydrocarbon coming from the bactericide, the Isotech lab ran a test and kindly shared the results for this study (after sitting 24 h with bactericide the water concentration in C 6+ does not changed) and, with the same purpose, the Isolab b.v. injected the used bactericide in the GC-FID. Isolab b.v. also ran a test injecting pure toluene in the GC, in order to confirm the compound vapor pick.

Results and Discussion
As shown in Figure 6, toluene was detected in all analyzed water samples with 93.5% of the values between 1000 and 6000 µg/L and 55% higher than 3000 µg/L. Figure 7 exemplifies three chromatograms of contaminated water samples (W06, W17, and W23) and the chromatogram of pure toluene. The amounts of toluene in groundwater samples are significantly higher than the values in soil samples, which in general are less than 1000 µg/L. Taking into account that the source of this contamination is the natural petroleum systems, more toluene will likely be found in the groundwater than in the surface (soil). In order to evaluate the spatial distribution of toluene in groundwater, a thematic map was designed (Figure 8), conjugating geological information (faults and formations) with two of concentration values: lower and higher than 1000 µg/L, for water and soil samples.
It is important to note that these amounts were obtained from gas (headspace) phase. This means that, following the partition coefficient of toluene in water-air systems at 20 • C (K = 4.39), the concentration of compound in the original sample (liquid phase), before analysis, is at least 4.39 times bigger [45][46][47]. According to Tipler [47], compounds with a high value for partition coefficients (>1) will favor the liquid phase, whereas compounds with low coefficients (<1) will favor the headspace phase.
The highest level reported in USA groundwater, in 1983, was 1.4 µg/L [2]. Toluene has been reported in industrial effluents in Canada. Raw process effluent from petroleum refineries in Ontario contained an average concentration of toluene of 0.6 µg/L. The average daily concentration for the refinery with the greatest discharge of toluene was 2.1 µg/L (maximum of 17.1 µg/L); the average daily loading was 0.05 kg/day in the process effluent. Still, in Ontario, samples taken in a contaminated shallow aquifer at a depth of 6 m beside an existing industrial chemical waste disposal lagoon, concentrations of toluene above 3900 µg/L were reported [10]. In Europe, maximum toluene levels of 1100 µg/L were reported in groundwater contaminated by point emissions [48].
In a total of more than 800 water samples taken across Canada from 1985 to 1988, concentrations of toluene in only six samples were greater than 0.5 µg/L. These included one surface water sample (0.9 µg/L), one drinking water sample (0.6 µg/L), two groundwater samples (0.6 and 3.9 µg/L), and two samples of undiluted sewage treatment plant effluent (31 and 32 µg/L) [10].
Individual gasoline hydrocarbons were detected infrequently in aquifers in USA. Only one gasoline hydrocarbon, toluene, was detected in more than 1% of aquifer samples, according to Zogorski et al. [17], and in approximately 1% of all groundwater-derived public drinking-water systems in the USA, toluene levels are above 0.5 µg/L [2].
Concentrations of 0.8 µg/L and 1.9 µg/L have been reported in the Rhine River in Germany and Switzerland, respectively [49]. Concentrations in the Morava River in Slovakia range from a winter maximum of 0.58 µg/L to a summer maximum of 3.49 µg/L [2]; in Spain, levels as high as 22 µg/L have been detected at the mouth of the Besos River [50]. In coastal waters, maximum levels of 1 µg/L have been found [2].
Geosciences 2018, 8, 9 8 of 18 lagoon, concentrations of toluene above 3900 μg/L were reported [10]. In Europe, maximum toluene levels of 1100 μg/L were reported in groundwater contaminated by point emissions [48]. In a total of more than 800 water samples taken across Canada from 1985 to 1988, concentrations of toluene in only six samples were greater than 0.5 μg/L. These included one surface water sample (0.9 μg/L), one drinking water sample (0.6 μg/L), two groundwater samples (0.6 and 3.9 μg/L), and two samples of undiluted sewage treatment plant effluent (31 and 32 μg/L) [10].
Individual gasoline hydrocarbons were detected infrequently in aquifers in USA. Only one gasoline hydrocarbon, toluene, was detected in more than 1% of aquifer samples, according to Zogorski et al. [17], and in approximately 1% of all groundwater-derived public drinking-water systems in the USA, toluene levels are above 0.5 μg/L [2].
Concentrations of 0.8 μg/L and 1.9 μg/L have been reported in the Rhine River in Germany and Water samples compounds between 4 and 12 carbon atoms per molecule, including toluene) originate thermogenically and are therefore associated with petroleum systems [41]. Following the natural contamination point of view, Collins [52] believes that natural concentrations of toluene and benzene in water are hydrochemical indexes that indicate the presence of petroleum system in subsurface. Accumulations of oil and gas can give rise to the presence of a number of low molecular weight hydrocarbons (petroleum is the main source of benzene and toluene in oilfield waters) that have l According to Seneshen et al. [53], there are strong toluene anomalies in soils over the northwest part of Lisbon field (San Juan County, UT, USA) that are upwind of hydrocarbon production. Natural organic contamination of groundwater with BTEX (benzene, toluene, ethyl-benzene, and xylene) are abundant next to oil and gas fields in the Gulf of Mexico [54,55].
No significant benzene levels were detected above GC-FID detection limit, but this absence can be explained by the biodegradation and oxidation in water, which can be more active in benzene than in toluene. According to Prince & Walter [16], it is reasonable to expect that the biodegradation of a broad range of compounds will be enhanced in the presence of toluene. Gülensoy & Alvarez [56] reported enhanced aerobic biodegradation of benzene, p-xylene, and naphthalene by several pseudomonads when toluene was present, and a similar phenomenon occurred under anaerobic conditions [57]. According to the ARH Alentejo [51], a measurement of the toluene levels was carried out in the industrial effluent from the Sines refinery, directed for treatment at Ribeira de Moinhos sewage treatment. The results shown toluene levels below 1 µg/L (Figure 2-the refinery is located 62 km from the nearest sampling point and 83 km from the farthest sampling point, on the NW of the sampling zone).
From the health impact point of view and following the health-based guideline (700 µg/L) suggested by the World Health Organization [13], toluene levels of the water for human consumption, in the study area, are of concern. Following the recommended WHO value, considering that the lower measured values (<1000 µg/L) are inserted into this limit, and the values from gas phase are the same as in the original sample (toluene favor the liquid phase in water-air headspace system), the groundwater in the sampling zone is contaminated by toluene (Figures 6 and 8). The analysis of the values mentioned above highlight the degree of contamination reported in this study. The detected levels are much higher than is usually found in groundwater, being only comparable with severe industrial contamination [10,15].
Light hydrocarbons (C 1 to C 4 ) very often originate from both thermogenic and biogenic processes. Normally, the hydrocarbons in the gasoline range (that consist of volatile organic compounds between 4 and 12 carbon atoms per molecule, including toluene) originate thermogenically and are therefore associated with petroleum systems [41]. Following the natural contamination point of view, Collins [52] believes that natural concentrations of toluene and benzene in water are hydrochemical indexes that indicate the presence of petroleum system in subsurface. Accumulations of oil and gas can give rise to the presence of a number of low molecular weight hydrocarbons (petroleum is the main source of benzene and toluene in oilfield waters) that have low odor thresholds in drinking-water [13,52].
According to Seneshen et al. [53], there are strong toluene anomalies in soils over the northwest part of Lisbon field (San Juan County, UT, USA) that are upwind of hydrocarbon production. Natural organic contamination of groundwater with BTEX (benzene, toluene, ethyl-benzene, and xylene) are abundant next to oil and gas fields in the Gulf of Mexico [54,55].
No significant benzene levels were detected above GC-FID detection limit, but this absence can be explained by the biodegradation and oxidation in water, which can be more active in benzene than in toluene. According to Prince & Walter [16], it is reasonable to expect that the biodegradation of a broad range of compounds will be enhanced in the presence of toluene. Gülensoy & Alvarez [56] reported enhanced aerobic biodegradation of benzene, p-xylene, and naphthalene by several pseudomonads when toluene was present, and a similar phenomenon occurred under anaerobic conditions [57].
It is also important to note the presence of other two hydrocarbons in soil-water samples in the study area. The significant presence of 2-methylpentane ( Figure 9 exemplifies three chromatograms of contaminated soil samples, SPZ-07, SPZ-15, and SPZ-16) in all soil samples, that, according to Seneshen et al. [53], is anomalous in two adjacent samples over a distance of 150 m and in one isolated sample over a gas field in the Lisbon field (San Juan County, UT, USA). Two water samples (W33 and W34) from two different boreholes (a few meters apart) contain worrying levels of methane (CH 4 ). According to Eltschlager et al. [33], recommended action levels for CH 4 dissolved in water suggest that when the level is between 10,000 and 28,000 µg/L, monitoring is required, and immediate action is needed above 28,000 µg/L. The sample W33 has a concentration of 12,000 µg/L of CH 4 , and W34 contains 70,000 µg/L (in the headspace phase), presenting a serious risk of borehole explosion, as the water is fully saturated with methane and the air is now flammable.
According to WHO [13], the maximum contaminant level (MCL), or the health-based guideline, for toluene in drinkable water is 700 µg/L (Table 1). It is readily absorbed completely from the gastrointestinal tract after ingestion and rapidly distributed in the body, with a preference for adipose tissue then for the kidneys, liver, and brain [1,13]. Toluene is rapidly metabolized and, following conjugation, excreted predominantly in urine [13]. Therefore, potential health effects from long-term exposure above the toluene MCL are in nervous systems, brain, kidney, and/or liver, with problems with fatigue and drowsiness being the most obvious symptoms [1,58].
The existence of the vertical or near-vertical migration of hydrocarbons from petroleum systems has been controversial for many years. Arguments against it have included the following. (1) If reservoirs leaked enough to be detected, they would have been depleted a long time ago; (2) If seepage did occur, it would be deflected by groundwater movements, and anomalies would be displaced laterally for relatively long distances [27]. According to Hale [40], the ways that gases move through subsurface natural environment are more diverse and less well characterized than mechanisms of solid and liquid phases dispersion. However, in the gas phase percolation, diffusion and mass flow are both well-established mechanisms and each has a role to play in gas transport in the subsurface. Their relative contributions to gas migration at different positions in the crust are less certain. The existence of the vertical or near-vertical migration of hydrocarbons from petroleum systems has been controversial for many years. Arguments against it have included the following. (1) If reservoirs leaked enough to be detected, they would have been depleted a long time ago; (2) If seepage did occur, it would be deflected by groundwater movements, and anomalies would be displaced laterally for relatively long distances [27]. According to Hale [40], the ways that gases move through subsurface natural environment are more diverse and less well characterized than mechanisms of solid and liquid phases dispersion. However, in the gas phase percolation, diffusion and mass flow are both well-established mechanisms and each has a role to play in gas transport in the subsurface. Their relative contributions to gas migration at different positions in the crust are less certain.
A more controversial mechanism, applicable to gas transport through groundwater, is gas streaming [44]. Buoyant hydrocarbon microbubbles rising through the water-filled network of fractures, joints, and bedding planes provide a logical mechanism for vertical migration of hydrocarbons from reservoirs to the surface and commonly occupy the pore voids in rocks, overburden, and soil [27,44]. Kristiansson & Malmqvist [50,59] and Malmqvist & Kristiansson [60,61] hypothesize that, in the zone of saturation, pressure gradients and pressure shocks cause oversaturation, leading to the formation of gas bubbles. These stream upward at a comparatively rapid rate until they reach the water table (any hydrocarbon that dissolves in groundwater could, given the appropriate conditions, migrate by streaming) and mix with the soil air. The resulting mixture is then driven slowly further upward by the pressure gradient caused by the bubble stream [44].

Conclusions
The toluene contamination levels, in the water samples from Beja and Faro Districts (1000-6000 μg/L), are much higher than is usually found in other places (even considering just the headspace phase concentration), being only comparable with values from severe industrial contamination. Possible anthropogenic/industrial contamination is practically discarded because the sampling area is isolated, sparsely populated, and with few gas stations. Furthermore, the presence of 2-   [13].

Occurrence
Concentration of a few micrograms per litre have been found in surface water, groundwater and drinking-water;point emissions can lead to higher concentrations in groundwater (up to 1000 µg/L);It may also penetrate plastic pipes from contaminated soil TDI 223 µg/kg body weight based on a LOAEL of 312 mg/kg body weight per day for marginal hepatotoxic effects observed in a 13-week gavage study in mice, adjusting for daily dosing and using na uncertainty factor for 1000 (100 for interspecies and intraspecies variattion and 10 for the short suration of the study and use of a LOAEL instead of a NOAEL) Limit of Detection 0.13 µg/L by GC with FID;6 ppb by GC-MS A more controversial mechanism, applicable to gas transport through groundwater, is gas streaming [44]. Buoyant hydrocarbon microbubbles rising through the water-filled network of fractures, joints, and bedding planes provide a logical mechanism for vertical migration of hydrocarbons from reservoirs to the surface and commonly occupy the pore voids in rocks, overburden, and soil [27,44]. Kristiansson & Malmqvist [50,59] and Malmqvist & Kristiansson [60,61] hypothesize that, in the zone of saturation, pressure gradients and pressure shocks cause over-saturation, leading to the formation of gas bubbles. These stream upward at a comparatively rapid rate until they reach the water table (any hydrocarbon that dissolves in groundwater could, given the appropriate conditions, migrate by streaming) and mix with the soil air. The resulting mixture is then driven slowly further upward by the pressure gradient caused by the bubble stream [44].

Conclusions
The toluene contamination levels, in the water samples from Beja and Faro Districts (1000-6000 µg/L), are much higher than is usually found in other places (even considering just the headspace phase concentration), being only comparable with values from severe industrial contamination. Possible anthropogenic/industrial contamination is practically discarded because the sampling area is isolated, sparsely populated, and with few gas stations. Furthermore, the presence of 2-methylpentane in all soil samples, and high concentration of methane in water samples from boreholes, are an indication that the natural unconventional petroleum system in SPZ is probably responsible for the groundwater and soil contamination by hydrocarbons.
Given that the geographical dispersion of soil and groundwater contamination is not yet known, the results presented in this paper are relevant for the Portuguese and Spanish Governments (SPZ formations are inserted in both territories). According to the MCL for toluene in drinkable water (700 µg/L), determined by WHO [13], contaminated groundwater represents a serious concern to the health of the residents, who should be informed to avoid drinking it or using it for other domestic purposes. Artesian wells with high methane levels represent a serious potential for borehole explosion, and immediate action is necessary to avoid this consequence.