The Geological Characteristics of the Vadose Zone Influence the Impact of Treated Wastewater on the Groundwater Quality (SCA.Re.S. Project 2019–2020)

This study evaluated whether some chemical and microbial contaminants in treated sewage effluents from two wastewater treatment plants (WWTP) reached the groundwater when they drained through a fractured karst vadose zone (WWTP-K) and a porous vadose zone (WWTP-P). Forty-five samples of sewage water (SW), treated water (TW), and monitoring well (MW), collected from WWTP-P (24) and WWTP-K (21), were analyzed for a range of microbiological and chemical properties. The E. coli and Salmonella counts were below the limits outlined in the Legislative Decree 152/06 in effluents from both types of WWTP. Enteric viruses were found in 37.5% and 12.5% of the SW and TW from WWTP-P, respectively. The percentages of Pepper mild mottle virus isolated were higher in TW (62.5% in WWTP-P, 85.7% in WWTP-K) than in SW and MW. The residual concentrations of contaminants of emerging concern (CEC) of each drug category were higher in the MW downstream of WWTP-K than of WWTP-P. Our results showed that the porous vadose zone was more effective at reducing the contaminant loads than the fractured karst one, especially the CEC, in the effluent. The legislation should include other parameters to minimize the risks from treated effluent that is discharged to soil.


Introduction
Worldwide, groundwater represents the main source of water for drinking purposes and agricultural and livestock needs [1]. In arid and semiarid areas, groundwater is used to meet the needs of the population where surface water resources are scarce or overexploited. Table 1. Chemical and physical parameters, reported as mean values (± standard deviation), measured in the sewage water (SW), treated water (TW), and monitoring well (MW) samples in the porous aquifer (WWTP-P), and karst fractured aquifer (WWTP-K).

Parameters
Unit of Measure All the values of biochemical oxygen demand (BOD 5 ) and chemical oxygen demand (COD) for the treated wastewater (TW) were below the limits allowed by Lgs. D. 152/06. This leads us to affirm that, in general, the wastewater treatment process in both plants proved to be effective against organic substances. The parameters relating to the nutrients contained in TW showed that the total nitrogen (TN) was sufficiently reduced, while total phosphorus (TP) was slightly higher than the limit (2 mg/L) at both the treatment plants.
Overall, the values of all parameters decreased significantly (p < 0.05), from sewage water (SW) to the monitoring wells (MW). However, it should be noted that the value for COD, an indicator of organic pollution, was always negative in the MW of WWTP-P, while the mean value was of 2.5 ± 3.1 mg/L in WWTP-K with a significant difference (p < 0.05) between the two plants.

Bacteria
The results from tests for the bacterial parameters in the WWTP effluents, expressed as mean concentrations (± SD), are presented in Table 2. The mean loads of E. coli were always below the limit allowed by Italian Lgs. D. 152/06 (5000 cfu/100 mL). E. coli and enterococci were detected in the effluents from both WWTP-P (893 ± 1138) and WWTP-K (6 ± 5) but were not detected in the respective MWs.
C. perfringens was detected in TW and MW (in one of eight samples) from WWTP-P, and was detected in TW from WWTP-K. This microorganism results in being more resistant in the environment and to disinfection treatments than coliforms and enterococci [32]. It is the main etiological agent of myonecrosis of connective tissues [33] and may also be responsible for food poisoning and diarrhea [34]. Table 2. Average concentrations (±SD) of bacterial parameters in the treated water (TW) and monitoring well (MW) samples of the porous (WWTP-P) and the karst fractured (WWTP-K) aquifers. P. aeruginosa was higher in the MW samples than in the TW samples at both the porous and karst sites. It is a free-living bacterium that mainly occurs in soil and natural waters (lakes, rivers, and groundwater), therefore the concentration of P. aeruginosa increased during the filtration process both in porous and karst fractured aquifer. It is one of the most important opportunistic pathogens of humans involved in a variety of infections, including respiratory and urinary tract infections, wound and soft tissue infections [35].
Salmonella spp. (results not shown) was not detected either at WWTP-K or WWTP-P. The transport of bacterial cells through a porous soil is influenced by several mechanisms such as physical straining, which in turn depends on the pore size, the bacterial cell size, and the hydraulic load rate, as well as adsorption to porous soils, which, on the other hand, is affected by the content of organic matter, the development of biofilm and electrostatic attraction due to ion strength of the solution or electrostatic charges of celland particle surfaces [36,37]. Furthermore, the concentration of viable cells in groundwater may depend on the inactivation rate of microorganisms during the percolation phase which is influenced by abiotic and biotic factors such as moisture content, pH, temperature, organic matter, bacterial species, predation, and antagonistic symbiosis among microorganisms [38,39].

Enteric Viruses
The results from testing for viruses in the samples from both plants by nested PCR are shown in Table 3. In WWTP-P, enteric viruses were detected in 37.5% of SW samples and 12.5% of TW samples but were not detected in MW. NoV-GII was detected in SW (12.5%) and TW (12.5%), while NoV-GI (12.5%) and EV (12.5%) were only detected in SW. Table 3. Percentage of virus positive samples in the sewage water (SW), treated wastewater (TW), and monitoring well (MW) samples from the porous (WWTP-P) and the karst-fractured (WWTP-K) aquifers, determined by PCR assay. AdV In WWTP-K, enteric viruses were found in 85.7% of SW samples, 57.1% of TW samples, and 14.3% of MW samples. Multiple virus families were detected in 85.7% of SW samples (up to three virus species in a single sample) and 28.5% of TW samples (up to two virus species in a single sample). NoV-GII and AdV were the most commonly detected viruses in SW (71.4%) and TW (42.9%). AdV was also found in MW (14.3%), while EV (57.1%) and NoV-GI (28.6%) were only detected in SW.
HEV, HAV, and RoV were not detected in any samples. Pepper mild mottle virus (PMMoV) was isolated in all samples of both WWTPs except for MW of WWTP-P and was isolated at a higher percentage in TW (62.5% in WWTP-P, 85.7% in WWTP-K) than in the other samples.
Six groups of enteric viruses were selected from the wide range of viruses excreted in human waste for this study because of their epidemiological significance as waterborne pathogens. Of these, only three, AdV, NoV, and EV, were detected in the analyzed water samples. It is important to note that molecular methods are only capable of detecting viral genomes, and do not provide information about infectivity, and thus positive samples do not necessarily indicate an actual threat to human health.
Viruses were present in the WWTP-P monitoring wells, perhaps because karst aquifers, unlike fissured and porous aquifers, are characterized by holes and channels through which water can flow rapidly, thereby transporting contaminants over great distances with little attenuation [27]. The degree of virus attachment is affected by several factors, including viral surface properties, groundwater quality, sediment surface charges, and soil moisture [40]. Further research will be carried out at the lab scale as part of the SCA.RE.S project to examine the processes of adsorption and/or inactivation of viruses through the unsaturated zone before reaching the groundwater.

Contaminants of Emerging Concern
Using targeted screening (Section 3.6), 35 CECs and their degradation products were identified and quantified in the TW and MW samples. All the substances were grouped into pharmacological categories, namely antimicrobials, anticonvulsants, nonsteroidal antiinflammatory drugs, beta-adrenoceptor blocking agents, UV filters, antipsychotic drugs, antihistaminic drugs, antidiabetic drugs, and X-ray contrast media.
The average concentrations of the CECs detected in the TW and MW samples of the WWTPs are listed in Table 4, along with the residual concentrations (RC) of each contaminant, expressed as the percentage of the detected concentration in MW compared with its concentration in TW. The average of the sum of the substances grouped by pharmacological category (in bold) and the respective RCs are also listed.   The concentrations of the detected CECs in the effluents of the two WWTPs did not differ significantly (p > 0.05) (Figure 1), perhaps because the two WWTPs use similar conventional processes to treat sewage. Moreover, studies have shown that the quantified contaminants are not efficiently removed during conventional wastewater treatment and are usually detected in secondary wastewater effluents at concentrations ranging from µg L −1 to ng L −1 [24].
Previous studies have reported that non-steroidal anti-inflammatory drugs (NSAIDs), b-blockers, antimicrobials, and anticonvulsants are widely used by the population and are widespread in WWTP influents [41,42]. Carbamazepine and lamotrigine were the main representatives of the anticonvulsant category, and they were detected at almost the same concentrations in the TW samples. Of the antimicrobials, fluconazole and levofloxacin were detected at the highest concentrations.
Venlafaxine, a widely used antipsychotic, and fluconazole were included in the Third Watch List of the Water Directive Framework through the Decision (EU) 2020/1161. In the last EU Decision, there was no final conclusion about these substances, so it is important to monitor for their presence in the aquatic environment, and to have information to support a decision about their possible inclusion in the List of Priority Substances. Pathogens 2022, 11, x FOR PEER REVIEW 8 of 15

Study Scenarios
This project was carried out in two WWTPs of the Salento Peninsula (Southern part of Apulia Region, Italy) characterized by vadose zones of different compositions: one in an area mainly comprised of fractured and karst limestone and the other in an area mainly comprised of porous carbonate rock.
The WWTPs located nearby the town of Carpignano Salentino (WWTP-P), and the town of Soleto-Galatina (WWTP-K) discharge effluent into infiltration pond systems from where it infiltrates through the vadose zone, thereby exploiting its natural capacity for filtering microbiological and chemical contaminants [46,47].
The infiltration ponds of the Carpignano plant have been dug into the carbonate porous rock of the Pliocene, known as the Sabbie di Uggiano Formation. The vadose zone below the ponds, about 60 m thick, is constituted by a porous rock that belongs to the Calcareniti of Andrano and Pietra Leccese Formations, both of Miocene age. The groundwater lies in a layered aquifer system, characterized by a fault system-oriented NW-SE that separates into two distinct groundwaters: the shallow one hosted in the Calcarenite of Andrano formation, while the deeper one in the cretaceous fractured limestone belongs to the Calcare of Altamura formation, locally known as the Dolomie of Galatina formation, and stratigraphically located below the Pietra Leccese ( Figure 2).
Instead, the bottom of the infiltration ponds of Soleto-Galatina plant corresponds to the top of the fractured and karstic limestone belonging to Calcare di Altamura Formation of Cretaceous age. Here, the vadose zone, about 57 m thick, is constituted by limestone, dolomitic limestone, and dolomite that are considerably diaclased, karstified, and often faulted. These rocks have a high degree of permeability. The groundwater lies in the deep carbonate aquifer and flows mainly in the fractures and within secondary karst conduits. Olmesartan, irbesartan, and flecainide (b-blockers) were found at concentrations above 1 µg/L. In a recent study, flecainide was detected in almost all the effluents tested [43]. Ladhari et al. [44] found that, because of their low degradation rate, they were often released into the environment and detected in different water sources.
Quantitatively, compounds in the beta-adrenoceptor blocking agent category, detected at average concentrations of 11.62 ± 1.14 µg/L in WWTP-P and 9.44 ± 2.98 µg/L in WWTP-K, were the most abundant in TW. Antimicrobials, with average concentrations of 1.53 ± 0.69 µg/L and 2.46 ± 1.27 µg/L in WWTP-P and WWTP-K, respectively, ranked second.
The average concentrations of 2-Phenyl-5-benzimidazolesulfonic acid (UV filter) were high in both WWTP-P and WWTP-K. The earlier study of WWTP-K showed that the concentrations of 2-Phenyl-5-benzimidazolesulfonic acid varied by season, as it was detected at different concentrations in summer and autumn.
The overall RC (%) results clearly highlighted that CECs behaved differently and had different fates in the porous and karst soils. The RCs of each pharmacological category were higher in the MW from WWTP-K than WWTP-P. These results probably reflect the hydrogeological characteristics of the areas where the plants are located.
WWPT-K is in a central-western area of Salento with karst fractured subsoil. This means that the groundwater is highly vulnerable, and contaminants can pass easily from the surface to the deep aquifer [45]. The groundwater velocity is usually high in karst aquifers that do not allow sufficient time for groundwater self-purification [11]. WWPT-P, on the other hand, is in the central-eastern part of the Salento Peninsula. The vadose zone is characterized by the presence of porous layers that could retain contaminants contained by the treated wastewater as it percolates through the subsoil. In the soil, CECs can undergo different processes, such as adsorption-desorption, transport, and biotic and/or abiotic transformations. The interactions between CECs and soil particles are also influenced by the properties of the soil (mainly the soil organic matter content, pH, clay content, and clay type) and the CEC molecular properties (such as lipophilicity, size, fractions of H-bond donors/acceptors, and charge) [42].

Study Scenarios
This project was carried out in two WWTPs of the Salento Peninsula (Southern part of Apulia Region, Italy) characterized by vadose zones of different compositions: one in an area mainly comprised of fractured and karst limestone and the other in an area mainly comprised of porous carbonate rock.
The WWTPs located nearby the town of Carpignano Salentino (WWTP-P), and the town of Soleto-Galatina (WWTP-K) discharge effluent into infiltration pond systems from where it infiltrates through the vadose zone, thereby exploiting its natural capacity for filtering microbiological and chemical contaminants [46,47].
The infiltration ponds of the Carpignano plant have been dug into the carbonate porous rock of the Pliocene, known as the Sabbie di Uggiano Formation. The vadose zone below the ponds, about 60 m thick, is constituted by a porous rock that belongs to the Calcareniti of Andrano and Pietra Leccese Formations, both of Miocene age. The groundwater lies in a layered aquifer system, characterized by a fault system-oriented NW-SE that separates into two distinct groundwaters: the shallow one hosted in the Calcarenite of Andrano formation, while the deeper one in the cretaceous fractured limestone belongs to the Calcare of Altamura formation, locally known as the Dolomie of Galatina formation, and stratigraphically located below the Pietra Leccese ( Figure 2).

Sampling
Both plants perform primary treatments, such as screening and sand separation, followed by secondary treatment such as denitrification, oxidation, and secondary sedimentation. Before being released onto the soil through a system of infiltration trenches, the effluents undergo a disinfection treatment (tertiary treatment), which is done by chlorination (sodium hypochlorite) in the WWTP-P, and by chlorination (sodium hypochlorite) and UV rays in the WWTP-K (Figure 3). Instead, the bottom of the infiltration ponds of Soleto-Galatina plant corresponds to the top of the fractured and karstic limestone belonging to Calcare di Altamura Formation of Cretaceous age. Here, the vadose zone, about 57 m thick, is constituted by limestone, dolomitic limestone, and dolomite that are considerably diaclased, karstified, and often faulted. These rocks have a high degree of permeability. The groundwater lies in the deep carbonate aquifer and flows mainly in the fractures and within secondary karst conduits.

Sampling
Both plants perform primary treatments, such as screening and sand separation, followed by secondary treatment such as denitrification, oxidation, and secondary sedimentation. Before being released onto the soil through a system of infiltration trenches, the effluents undergo a disinfection treatment (tertiary treatment), which is done by chlorination (sodium hypochlorite) in the WWTP-P, and by chlorination (sodium hypochlorite) and UV rays in the WWTP-K (Figure 3).

Sampling
Both plants perform primary treatments, such as screening and sand separation, followed by secondary treatment such as denitrification, oxidation, and secondary sedimentation. Before being released onto the soil through a system of infiltration trenches, the effluents undergo a disinfection treatment (tertiary treatment), which is done by chlorination (sodium hypochlorite) in the WWTP-P, and by chlorination (sodium hypochlorite) and UV rays in the WWTP-K (Figure 3). The quality of the wastewater before and after the secondary treatment and of the groundwater from wells close to the treatment plants was assessed from the chemical, physical, and microbiological properties. The quality of the wastewater before and after the secondary treatment and of the groundwater from wells close to the treatment plants was assessed from the chemical, physical, and microbiological properties.
Samples were collected monthly from May 2019-December 2020. On the collection day, samples were collected between 8:00 and 11:00 in calm atmospheric conditions, with no rain.
A total of 45 water samples (24 from WWTP-P and 21 from WWTP-K) were collected, comprising sewage water (SW) (8 from WWTP-P and 7 from WWTP-K), treated water (TW) (8 from WWTP-P and 7 from WWTP-K), and water from a monitoring well (MW) (8 from WWTP-P and 7 from WWTP-K). The monitoring wells were located no more than 500 m downstream from each WWTP and were used routinely to monitor the plant activity by local authorities.
Microbiological, chemical, and physical parameters, including some beyond those required by Italian law [22], were determined in the TW and MW samples. Samples of water for bacterial property analysis (2 L) were collected in sterile glass containers by using an automatic sampler. Samples of different volumes (2.5 L samples of SW, 40 L of TW, and 1000 L of MW) were also collected and analyzed for virus detection. Samples of TW and MW were collected and analyzed for chemical and physical properties (1 L) and CECs (200 mL) determination.

Detection of Bacterial Indicators
Bacterial indicators of fecal contamination (Escherichia coli, Enterococci, Salmonella spp., Pseudomonas aeruginosa, Clostridium perfringens) were assessed in the TW and MW water samples.
To detect Salmonella spp., a 1000 mL aliquot of each sample was filtered through a cellulose nitrate membrane filter (0.45-µm pore size) that was successively transferred into buffered peptone water (BPW) (Biolife Italiana srl, Milan, Italy) and incubated for 18-24 h at 36 ± 1 • C. A 0.1-mL aliquot of culture was inoculated into 10 mL of Rappaport Vassiliadis broth (Microbiol and C. s.n.c, Uta, Italy) and incubated at 41.5 • C for 24 + 24 h.
The broth was then streaked on xylose lysine deoxycholate agar plates (Biolife Italiana srl) and Hektoen Enteric Agar (Merck, Darmstadt, Germany). After incubating for 24 h at 36 ± 1 • C, the colonies with a typical morphology were sub-cultivated on Tryptic Soy Agar plates (Biolife Italiana srl) at 36 ± 1 • C and biochemically confirmed using the API 20E test (Biomèrieux, Marcy l'Etoile, France). Finally, colony typing was carried out using specific serological tests, as described in the APAT CNR IRSA manual 7080 [50].
To detect P. aeruginosa, 250 mL of each sample was filtered through a cellulose ester membrane filter (0.45-µm pore size). The membrane was placed onto a plate containing Pseudomonas selective agar supplemented with cetrimide (0.20 g) and nalidixic acid (15 mg) (Microbiol, Cagliari, Italy) and incubated at 36 ± 2 • C for 44 ± 4 h. Blue-green pyocyanin producing colonies were confirmed to be P. aeruginosa [51].
To detect C. perfringens, 100 ml of each sample was pre-treated at 75 ± 5 • C for 15 ± 1 min in a water bath and then filtered through a cellulose nitrate membrane filter (0.45-µm pore size). The membranes were placed on Tryptose Sulfite Cycloserine Agar (Biolife Italiana srl, Milan, Italy) and the plates were incubated at 44 • C for 24 h under anaerobic conditions (GasPak EZ Gas Generating Pouch Systems -BD Diagnostics, 7 Loveton Circle, Sparks, Maryland, USA). Black colonies were considered spores of sulfite-reducing clostridia [52].

Detection of Viruses
Molecular methods were used to test samples of SW, TW, and MW for different enteric viruses, namely adenovirus (AdV), norovirus genogroups I and II (NoV-I and NoV-II), enterovirus (EV), hepatitis A virus (HAV), hepatitis E virus (HEV), and rotavirus (RoV). The samples were also tested for the Pepper mild mottle virus (PMMoV).
Before molecular analysis, composite samples over a 24 h period were collected from the WWTP influent (SW) post the inlet screens (250 mL) and then concentrated using the two-phase (polyethylene glycol,PEG-dextran) separation method, in line with the method in the WHO guidelines for environmental surveillance of poliovirus circulation [53]. Samples from TW (40 L) and MW (1000 L) were filtrated by Nanoceram electropositive cartridges (Argonide Corporation, Sanford, FL, USA) and concentrated by the virus adsorptionelution (VIRADEL) technique, as described by Iaconelli et al., 2017 [54], and Montagna et al., 2020 [9]. Viral nucleic acid was extracted from concentrated samples after chloroform treatment (5 mL) using a semi-automated nucleic acid extraction platform (NucliSENS MiniMag, bioMerieux, Marcy l'Etoile, France), following the manufacturer's instructions. Extracted RNA was stored at −80 • C until the molecular analysis.
Enteric viruses and PMMoV were detected by reverse transcription-nested polymerase chain reaction (PCR), following the method described by Bonanno Ferraro et al., 2021 [55].
The PCR products were loaded onto 2% agarose gel containing GelRed staining (Biotium, Fremont, CA, USA). The amplified products were purified using a Montage PCRm96 Microwell Filter Plate (Millipore, Billerica, MA, USA) and subjected to Sanger sequencing on both strands (Bio-Fab Research, Rome, Italy). Consensus sequences from each sample were compared with those available in the GenBank database using BLAST [56].

Chemical and Physical Parameters
The samples were analyzed for the physical and chemical parameters required by Lgs D. 152/06 [22]. The temperature, pH, and conductivity were measured in situ with a multi-parameter probe (WTW MultiLine P4) with the potentiometric [57] and conductivity [58] methods. Total suspended solids (TSS) were determined by filtering samples through glass fiber filters [59], the biochemical oxygen demand (BOD5) was determined by the respirometric method [60], and the chemical oxygen demand (COD) was detected using the sealed tube method from the International Standard Organization [61]. The total nitrogen (TN) concentrations were measured by UV spectrometry after oxidative digestion with sodium persulfate using an actuator that operates in a coordinated analytical sequence [62]. The nitrate (NO 3 ) concentrations were determined from the dissolved anions by liquid phase ion chromatography [63], and the TP concentrations were determined with the ammonium antimony-phospho-molybdate colorimetric method after sequential analysis [64].

Contaminants of Emerging Concern
The TW samples from WWTP-P and the MW samples from the whole project period were examined for CECs. The detection was performed with a liquid chromate graphy system (Ultimate 3000; Thermo Fisher Scientific, Waltham, MA, USA) interfaced with a high-resolution mass spectrometer (TripleTOF® 5600+; AB Sciex, 500 Old Connecticut Path, Framingham, MA 01701, USA). The UPLC-QTOF/MS/MS conditions implemented for the detection and quantification of CECs followed a well-established method already described in the previous article [9].
Briefly, each sample was filtered through a 0.20 µm regenerated cellulose filter and spiked with an internal standard, i.e., CBZ D10, at a level of 10 ng/mL, and then was injected (2000 µL) using an optimized online solid phase extraction method (online SPE) that consisted of a pre-concentration column and an analytical column for analyte separation.
The data were processed with AB Sciex software. Each sample was screened for the presence of target compounds by comparing the retention time and the MS/MS fragmentation pattern of the detected compounds with those of the standard compounds. A reference calibration curve was inserted into the 0.01-5 µg/L range (0.01, 0.05, 0.1, 0.5, 1, and 5 µg/L) to quantify all the compounds identified in the samples.

Statistical Analysis
The results from the chemical and microbiological analysis were entered into a Microsoft Excel database and statistically processed using MedCalc Software version 12.3 (MedCalc Software bvba, Ostend, Belgium). The values of the microbiological parameters that were not detected were set at half (0.5 CFU/100 mL) of the detection limit (<1 CFU/100 mL), as described in Lorimer and Kiermer, 2007 [65].
The arithmetic mean, standard deviation, maximum, and the median were calculated for each group of quantitative variables. The D'Agostino-Pearson normality test was used to determine if the data sets conformed to a normal distribution, and the Tukey test was used to detect outliers. The RC of each contaminant of emerging concern in the MW was calculated as a percentage of the concentration of the same contaminant found in the TW.
The Student's t-test was used to determine if there were significant differences between the RCs of the substances grouped by pharmacological category at the different monitoring points.

Conclusions
The findings from the present investigation, combined with the findings from the earlier study, confirmed that, based on the chemical, physical, and microbiological measurements, the quality of the WWTP effluent discharged into the subsoil met the requirements of the current Italian Legislation. However, microbiological parameters beyond those required by the legislation, including viruses, were detected. The results demonstrate that parameters beyond Salmonella and E. coli should be included in the legislation that governs the discharge of effluent onto soil or into subsoil.
Of the karst subsoil and the porous soil, the porous one was better at retaining the chemical and microbial contaminants contained in the effluent from the WWTP-P, including emerging contaminants, probably because of its self-filtering capacity. This meant that effluent that had passed through the porous vadose zone was less likely to contaminate the groundwater than effluent that passed through the fractured and karst one. Thus, the risk of groundwater contamination is greater under a karst fractured vadose zone than under a porous vadose zone. In addition, the results confirm that pollutants are transported and diffuse more rapidly through a karst fractured domain than through a non-karst zone, without sufficient time to be retarded by chemical reactions or other attenuation mechanisms.
To date, few studies have considered how factors such as the precipitation rate, moisture, and seasonality influence the transport of effluents containing microbes through a non-karst vadose zone to the groundwater. A sound understanding of these processes would be beneficial for protecting the groundwater to ensure it was suitable for irrigation and drinking purposes. Therefore, further lab-and field-scale studies will be carried out in the SCA.RE.S project to investigate (1) the viability of enteric viruses in the WWTP effluents and (2) the variability of the groundwater quality recharged by effluent passing through vadose zones of different depths and lithology. Funding: This research was funded by the Apulia Regional Government, Bari, Italy (DGR 1346 del 04/08/2021).

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable. Data Availability Statement: All data generated or analyzed during this study are included in this published article.