Next Article in Journal
Evaluation of Humanitarian Supply Chain Resilience in Flood Disaster
Next Article in Special Issue
Proposing the Optimum Withdrawing Scenarios to Provide the Western Coastal Area of Port Said, Egypt, with Sufficient Groundwater with Less Salinity
Previous Article in Journal
Stochastic Modeling of Hydroclimatic Processes Using Vine Copulas
Previous Article in Special Issue
Application of 2H and 18O Isotopes for Tracing Municipal Solid Waste Landfill Contamination of Groundwater: Two Italian Case Histories
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Irrigation with Coal Mining Effluents: Sustainability and Water Quality Considerations (São Pedro da Cova, North Portugal)

1
National Institute of Health Doutor Ricardo Jorge, Department of Environmental Health, University of Porto, 4050-083 Porto, Portugal
2
LAQV/REQUIMTE, Associated Laboratory for Green Chemistry of the Network of Chemistry and Technology, University of Porto, 4050-083 Porto, Portugal
3
Institute of Earth Sciences, Porto Pole, University of Porto, 4169-007 Porto, Portugal
4
Department of Geosciences, Environment and Spatial Planning, University of Porto, 4169-007 Porto, Portugal
5
Department of Earth Sciences, University of Coimbra, 3030-790 Coimbra, Portugal
*
Author to whom correspondence should be addressed.
Water 2021, 13(16), 2157; https://doi.org/10.3390/w13162157
Submission received: 21 June 2021 / Revised: 27 July 2021 / Accepted: 31 July 2021 / Published: 5 August 2021

Abstract

:
Two water effluents that drain from the abandoned coal mine of São Pedro da Cova (NW Portugal) were characterized in terms of their physic-chemical properties and suitability for irrigation purposes. Samples were also collected in a local surface stream, upstream and downstream from the mine drainage points, also used for irrigation by local farmers. Water samples were analyzed for major and minor ions and for trace element concentrations. Sampling campaigns started in 2017 and ended in 2019 and there were 46 water quality parameters tested. There were also proposed all-inclusive indices (the Water Quality Index and the Contamination Index, and also the Trace Element Toxicity Index) based on specific groups of 18 and 17 physic-chemical parameters, respectively, to achieve adequate monitoring requirements for mine effluents and surface water from coalfield. From the physical and chemical aspects of mine water it is inferred that the mine is not producing acid mine drainage. The coal mine water is of medium to high salinity, having almost neutral pH and a high thermal stability during the year, which is a distinguishing feature of the effluents. When compared to international irrigation water quality standards, as Food and Agriculture Organization of the United Nations admissible concentrations, the impacted waters are unsuitable for irrigation. The major outliers to the guidelines were iron, manganese, potassium, magnesium and bicarbonates, being also detected carcinogenic polycyclic aromatic hydrocarbons. Cost-effective ways of monitoring water quality parameters are needed to help control and manage the impact of coal mine effluents that should be treated before releasing into a ditch system that could be then used by local farmers to irrigate their crops.

1. Introduction

Irrigation is fundamental for agriculture but policies that push towards a restrained use of water are not popular among farmers, who are also not prepared to respond to drastic increases in water costs, which could decrease the economic profitability of their activities. Therefore, the use of unconventional free water sources, such as mine wastewaters, is regarded as a possible choice for irrigation. Effluents from coal mines are often considered severe and persistent forms of pollution, with environmental impact not only throughout the mine’s life cycle, but also long after the end of mining activities. Some deleterious impacts on the environment include the disruption of hydrological pathways, contamination of surface and groundwater, depression of the water table, soil contamination and loss of biodiversity [1,2]. The concentration of contaminants in coal mine drainage waters vary greatly and depends on a series of geological, hydrological and mining conditions, which are different from mine to mine. Therefore, the effluents can be alkaline, acidic, ferruginous, highly saline, or even clean [3]. Frequently, coal mine drainage is acid metal-rich waters with high concentrations of iron (Fe), copper (Cu), manganese (Mn) and nickel (Ni), formed during water–rock interaction involving sulfur-bearing minerals, such as pyrite (FeS2). These processes may cause red, orange, or yellow sediments with negative impacts on ecosystems and water resources. Additionally, mine effluents often contain high levels of total dissolved and suspended solids. The dissolved cations include mainly calcium (Ca), magnesium (Mg), sodium (Na) and potassium (K); the major anions are sulfate (SO4), chloride (Cl), fluoride (F), nitrate (NO4), bicarbonate (HCO3) or carbonate (CO3) [4].
The Douro Coalfield (NW Portugal) represents the most important coal-bearing deposit in Portugal [5,6,7,8] with 53 km length and width between 30 and 250 m (Figure 1).
São Pedro da Cova mining area was one of the principal centers of mining activity in Portugal (Figure 2a,b), with great economic and technological impacts, and a cultural significance from the end of the 17th century (1795), until the 20th century (1970) [9].
The São Pedro da Cova coal mine is an abandoned mine located in a peri-urban area, with a landscape consisting of a mosaic of urban, industrial, agricultural and forest areas (Figure 2a,b). The mine is located very close to a densely populated zone with several basic facilities including schools, a healthcare center and a leisure center.
The coal mining effluents from two mine drainage galleries are discharged around 1 km to SE of the mine, producing an ocher-colored sediment that is continuously accumulated in local watercourses (Silveirinhos stream and Ferreira river). Water from Silveirinhos stream is collected downstream from the mine drainage discharges and is locally used for agriculture irrigation. Local farmers have been using this polluted water for decades to produce a wide variety of products, including green vegetables, corn, barley and fruits (Figure 2c–g).
This study aimed at assessing the impact of coal mining effluents on the quality of irrigation water, by means of water quality individual parameters and indices.
This approach intended to contribute to the development of a specific evaluation methodology for coal mining regions in terms of i) the suitability of this type of water for irrigation purposes and ii) the environmental impacts resulting from the uncontrolled applications of these unconventional water resources in irrigated areas.

2. Assessment of Irrigation Water Quality: A Short State of the Art

Mine effluents are used for irrigation purposes worldwide, as they constitute an easily accessible and inexpensive source of water.
The quality of irrigation water is considered a key factor for safe food production [11]. However, when good quality water is scarce, water of marginal value is often considered for use in agriculture [12]. The prolonged and uncontrolled use of polluted water could result in reduced crop yields, deterioration of soil properties and severe environmental and health damages, requiring more complex management practices and more stringent monitoring procedures [13]. Therefore, the agricultural use of water resources in areas affected by mining activities requires not only a baseline water quality data, but also continuous monitoring.
There are several quality standards and guidelines for irrigation water proposed by various countries and organizations, combining conservative and liberal approaches. However, despite being useful, they are not always satisfactory due to the wide variability of hydrological and hydrogeological settings [14,15,16,17].
In 1985, the Food and Agriculture Organization of the United Nations (FAO) published a document concerning water quality for agriculture, which was reprinted in 1994, presenting a set of guidelines modified to give more practical procedures for evaluating and managing water quality-related problems, emphasizing long-term effects [18]. These general water quality classification guidelines help identify potential crop production problems associated with the use of conventional water sources, and are equally applicable to evaluate wastewaters for irrigation purposes in terms of their chemical constituents. As wastewater effluents may contain a number of toxic compounds, FAO also presented threshold levels for some selected trace elements.
FAO’s guidelines have been widely incorporated into national regulations all over the world and the following criteria have been considered the most relevant in defining quality [13,15,18]:
(i) Salinity hazard: the concentration of soluble salts in irrigation water, estimated in terms of electrical conductivity. Salinity has been deemed as the most important factor of irrigation water quality because high salinity in soil can create a hostile environment for the crop to absorb nutrients and cause specific ion toxicity;
(ii) Infiltration and permeability problems: the two most common water quality factors which influence the normal infiltration rate are water salinity (total quantity of salts in the water) and the sodium content relative to the calcium and magnesium content;
(iii) Specific toxicity hazard: certain ions and metals can accumulate in sensitive crops in concentrations high enough to cause damage and reduce yields. The ions of primary concern are boron, sodium and chloride;
(iv) Miscellaneous effects: these include high nitrogen concentrations in water, which supplies nitrogen to the crop and may cause excessive vegetative growth, lodging and delay crop maturity; unsightly deposits on fruit or leaves due to overhead sprinkler irrigation using water with high bicarbonate, gypsum, or iron contents; various abnormalities often associated with an unusual water pH.
In order to better classify water based on its specific characteristics and possible uses, a number of models of water quality, named Water Quality Indices (WQI), have also been developed since the 60s [19]. WQI are simplified representations of complex realities, used to assess the suitability of water for certain purposes based on specific characteristics. The use of WQI allows the representation of a large number of parameters in a single numerical value, facilitating the operational management of water resources and their allocation for different uses [20]. However, for this number to accurately represent the reality of a water body, the correct selection of environmental quality parameters is essential [21].
The first modern WQI was proposed by Horton (1965) [22] and was followed by numerous studies in this field. In recent years, many modifications have been considered in the WQI concept and several indices have been proposed and used by governmental agencies and researchers [19,20,23].
Misaghi et al. (2017) [24] introduced the first systematic WQI for irrigation purposes. However, this index considers a limited set of parameters for estimating water quality and does not take into account all potential impacting properties that could be critical, especially regarding wastewaters or other “marginal” quality waters.
The selection of variables is of major importance in calculating WQI as they should be independent and the most relevant ones, in order to define water quality and detect water quality deterioration. Several indices use parameters selected according to the opinion of specific expert panels, and the consequence is that the final evaluation can be highly subjective and variable. Other authors propose that the selection of parameters should be made according to the water management objectives, the location of the studied waters and the sampling periodicity [20].
To assess the impact of mining activities on irrigation water, some indices seemed to be more appropriate than others to be used as additional tools, as the Weight Arithmetic Water Quality Index method (WQIA), proposed by Cude (2001) [25], which is based on Horton’s principles but has been modified by introducing the normative values of the major factors of the water [20,25], the Contamination Index (CI), developed by Backman et al. (1998) [26,27,28] and the Trace Element Toxicity Index (TETI), by Ali et al. (2017) which is based on contaminant hazard intensity [28].
Furthermore, despite the usefulness and importance of normative guidelines, the effect of unusual or special water constituents is not always considered. An example is the contamination by polycyclic aromatic hydrocarbons (PAH), which are persistent, semi-volatile organic pollutants that can result from the oxidation and self-combustion of mine wastes and should be analyzed in mining surroundings due to their genotoxic, mutagenic and carcinogenic properties [29]. Although there are many PAHs, scientists and regulators have focused on 16 compounds that have been identified as priority control pollutants by the U.S. Environmental Protection Agency (USEPA). PAHs in underground mining environments may dissolve in mine water and eventually pollute the groundwater system. In addition, mine effluents can bring PAHs to the surface environment, polluting surface water and soils as well [30].

3. Materials and Methods

3.1. Study Area

The region of São Pedro da Cova (N 41°09′; W 8°30′) is situated in NW Portugal (Figure 1), in the city of Porto peri-urban region and has a resident population of around 16,500 inhabitants [31] and an area of 16 km2, where the territory comprises residential, industrial, and agroforestry areas (Figure 2a,b,g).
The region is located in the Central-Iberian Zone of the Iberian Massif [32]. The regional geological units consist of metasedimentary rocks (Figure 1) with minor sedimentary cover areas. For that reason, the prevailing groundwater circulation media are fractured.
The study region is located along the western limb of the Valongo Anticline, a regional megastructure [33], which created mountainous landforms reaching 367 and 385 m of altitude at the Santa Justa and Pias summits, respectively, on its western and eastern flanks. The regional morphology is dominated by two hill alignments that originated from differential erosion and are crosscut by Douro River. This structure controlled significantly the regional drainage network, that is part of the Douro river. In São Pedro da Cova it is possible to identify different metasedimentary formations with ages between the Precambrian and/or Cambrian, Ordovician, Silurian, Devonian and Carboniferous.
The old mine of São Pedro da Cova, where exploitation of anthracite A occurred for nearly 200 years, is located in one of the multiple coal deposits hosted in Douro Coalfield (from Upper Pennsylvanian), that represents the most significant Portuguese coal-bearing deposit. This deposit is elongated along NW–SE, presenting an approximate length of 53 km and variable width, ranging from 30 to 250 m. The sedimentary sequence comprises a basal breccia, followed by fossiliferous shales, siltstones and sandstones, along with interlayered conglomerates and coal seams [10].
The regional climate is Atlantic, the normal annual precipitation is around 1254 mm (with 195 mm in December and 18 mm in July) and the normal annual air temperature is around 15°C (with 20°C in July and August and 9°C in January). The Köppen-Geiger climate classification is Csb: warm temperate, with dry and warm summers [34,35].

3.2. Water Sampling

Groundwater samples were collected from two mine drainage galleries (G1 and G2). Additionally, surface water from Silveirinhos stream was sampled in two points, one located upstream (SS-U) and the other downstream (SS-D) from the mine galleries discharge (Figure 2d). Nine sampling campaigns were carried out in two periods, from April 2017 until December 2017 (Apr, Jun, Sep, Dec 2017), and from November 2018 until December 2019 (Nov 2018; Feb, May, Sep, Dec 2019), in a total of 34 samples. In the campaigns of October 2017 and September 2019 it was not possible to collect water from Silveirinhos stream upstream from the drainage galleries because, during the dry season, this part of the stream does not flow.
The sampling points were chosen in order to be preserved from other anthropogenic impacts besides coal mining, reducing the risk of water contamination by other pollution sources. Sampling was carried out according to standard methods: ISO 5667-3:2018 (E) Water quality—Sampling—Part 3 [36]. Samples were collected in glass or polyethylene bottles, were stored at low temperature (<5 °C) in the dark, and delivered to the laboratory within 5 h. The samples were collected with as little agitation or disturbance as possible. Special preservatives were required for some parameters. In this case, care was taken not to flush any preservative out of the bottle during sampling. Conditions such as strong winds or heavy rain were avoided during sampling.
Water samples were analyzed for a range of physical and chemical constituents in the laboratory, while temperature, pH and EC were measured in situ at the moment of sampling (Figure 2e), using a multiparametric meter from Hanna Instruments, model HI-991300, Woonsocket, RI, USA.
This set of sampling points provided a monitoring network for investigating the impact of mine drainage on environment, namely on irrigation water chemistry, including its seasonal variation.

3.3. Laboratory Analysis

Analyses were performed according to procedures outlined in Standard Methods for the Examination of Water and Wastewater 23rd edition [37] and in Le Rodier—L’analyse de l’eau 10e édition [38]. 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 (Hach Lange, Düsseldorf, Germany). Electrical conductivity (EC) and pH were determined in a Crison MultiMeter MM 41 (Hach Lange Spain, S.L.U., Barcelona, Spain). Total alkalinity, carbonates (CO32−) and bicarbonates (HCO3), were analyzed by titration. Phosphate (PO42−) was analyzed in a Shimadzu UV-1601 Spectrophotometer (Shimadzu Corporation, Kyoto, Japan). Chemical oxygen demand (COD) was evaluated in a Hach DR 2800 Spectrophotometer (Hach Company, Loveland, CO, USA). Major inorganic ions (Na+, K+, Mg2+, Ca2+, Li+, Cl, NO3, F and SO42−) were analyzed by ion chromatography (DionexTM system DX-120/ICS-1000, Dionex Corporation, Sunnyvale, CA, USA). Total organic carbon (TOC) was analyzed in a Shimadzu TOC-V (TOC-ASI-V, Shimadzu Corporation, Kyoto, Japan), heavy metals (Cr, Mn, Ni, Cu, Zn, As, Cd and Pb) and other components, such as Al, Fe, NO2, NH4+ and SiO2 were analyzed in a Varian AA240 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. 2018 [39].
Analytical standards were supplied by Sigma–Aldrich (Steinheim, Germany) and Merck (Darmstadt, Germany). The reference standard mixture containing the 16 EPA PAHs (naphthalene, Nap; acenaphthylene, Acy; acenaphthene, Ace; fluorene, Flu; phenanthrene, Phe; anthracene, Ant; fluoranthene, Flt; pyrene, Pyr; benz[a]anthracene, BaA; chrysene, Chr; benzo[b]fluoranthene, BbF; benzo[k]fluoranthene BkF; benzo[a]pyrene, BaP; dibenz[a,h]anthracene, DahA; benzo[ghi]perylene, BghiP; and indeno[1,2,3-cd]pyrene, Ind) was purchased from Sigma-Aldrich (Steinheim, Germany).
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).

3.4. Irrigation Water Quality Parameters

Water quality evaluation is necessary to assess the suitability of water to serve a specific purpose, and to determine appropriate treatments or precautions, if necessary. However, monitoring all parameters involved in a water source could be time-consuming and expensive. Therefore, reducing the subjectivity and the effective cost for assessing water quality is a great challenge.
This study focuses on the parameters adopted by FAO guidelines 29 [18] and on a set of quantitative assessment ratios which included the widely applied Sodium Adsorption Ratio (SAR), the Total Hardness (TH), the Residual Sodium Carbonate (RSC) and the Permeability Index (PI) (Table 1).
Moreover, in order to meet the hydrogeological and hydrogeochemical specificity of coal mining effluents, the water quality indices WQIA, CI and TETI were considered of significant importance and, therefore, were also calculated.
The water quality index based on the weighted arithmetic method (WQIA) was amended to be specific for irrigation, being adjusted taking into consideration the FAO recommendations, that is, the weights were defined as functions of the standards proposed in this guideline. For computation of the WQIA index, 18 water quality parameters were used, namely, the EC to estimate the salinity hazard; 3 elements with specific ion toxicity (Na+, Cl and B); 3 elements with miscellaneous effects (NO3, HCO3 and pH) and 11 trace elements with Recommended Maximum Concentration Values (Al3+, As, Cd2+, Pb2+, Cu2+, Cr3+, Fe2+, Mn2+, Ni2+, Zn2+ and F).
WQIA was calculated by using the following equation [20,43,44]:
WQI A = i = 1 n Q i W i / i = 1 n W i
The quality rating scale (Qi) for each parameter for a total of n water quality parameters is calculated by using this expression:
Qi = 100 V i V 0 / S i V 0
where, Vi is the actual value of the ith water quality parameter obtained from laboratory analysis, V0 is the ideal value of that water quality parameter obtained from standard Tables (V0 = 0, except for pH = 7.0) and Si is the recommended standard value of ith parameter.
The relative unit weight (Wi) for each water quality parameter is calculated by adopting the following formula:
W i = K / S i
where, K is the proportionality constant and can also be calculated by using the following equation:
K = 1 i = 1 n 1 / S i
The proposed index ranges from 0 to 100 and plain descriptions for index data were developed in order to provide a qualitative description of the index outcome [44,45]. The calculation of WQIA following the ‘weighted arithmetic index method’ involves the estimation of ‘unit weight’, assigned to each physic-chemical parameter considered for the calculation. By assigning unit-weights, all the concerned parameters of different units and dimensions are transformed to a common scale.
Weightage of each parameter means its relative importance in the overall water quality, and it depends on the permissible limits. Those parameters which have low permissible limits and can influence water quality to a large extent allocate high weighting, while parameters having high permissible limits are less harmful to the water quality and allocate low weighting.
Table 2 shows the irrigation water quality standards and the unit weights assigned to each parameter used for calculating the WQIA index. Maximum weights were assigned to cadmium (0.7135), arsenic and chromium (0.07135) and to copper, manganese and nickel (0.03568), thus suggesting the key significance of these trace elements in water quality assessment and their considerable impact on the index.
The CI was also calculated as a sum of the contamination factors of individual components (the 18 water quality parameters chosen for WQIA calculation, analyzed in the nine sampling campaigns), some of these exceeding the trigger values recommended by FAO [18]. The CI is determined by the following formula [26]:
CI = i = 1 n C Ai C Ni 1  
where CAi and CNi represent the analytical value and upper permissible concentration of the ith component, respectively. Note that CNi is taken as maximum allowable concentration.
Based on the CI index, values less than 1, 1–3, and more than 3 indicate low, medium and high levels of contamination, respectively [26].
TETI [28] was calculated based on the contaminants hazard intensity. The hazard intensity, or total score, of each parameter was determined according to the Toxicological Profiles of the Priority List of Hazardous Substances prepared by the Agency for Toxic Substances and Disease Registry (ATSDR), the Division of Toxicology and Environmental Medicine, Atlanta, GA, USA [46]. The ATSDR prioritization of substances is based on a combination of their frequency, toxicity, and potential for human exposure.
The concentration of each trace element detected was multiplied by its total score, and products were added to calculate TETI. The proposed TETI only considers toxic elements in water and is calculated by using the expression:
TETI = i = 1 n C i   × TS i
where Ci is the concentration of each individual trace element and TSi is its Total Score (ATSDR). This index clearly represents the impact of mine activities on the aquatic environment, where the lower TETI value represents better water quality.
In addition to the conventional parameters of water quality, studies on organic pollutants as PAHs are also very important, as several mining activities such as coal mining, processing or storage of coal provide the basic conditions for the generation and release of these compounds.
As PAHs can also become a source of pollution after the abandonment of coal mines, the 16 priority hydrocarbons were analyzed in São Pedro da Cova samples in order to investigate contamination levels and distribution [47].

3.5. Statistical Analysis

Data obtained for different parameters 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, the analyzed parameters were studied using a Kruskal–Wallis test. Whenever statistical significances were found, Dunn’s post-hoc test was applied for median comparison. All statistical analyses were performed at 5% significance level using R version 4.0.2 (R Project for Statistical Computing).

4. Results and Discussion

4.1. Water chemistry

The values of the physical and chemical parameters used to evaluate water quality from the two mine drainage galleries (G1 and G2, groundwater) and from Silveirinhos stream, collected upstream (SS-U) and downstream (SS-D) from the mine effluents discharge points are reported in Table 3. Results were compared to FAO guidelines in order to assess the water suitability for irrigation. The mining effluents G1 and G2 correspond to groundwater which circulates in the exploited rock massif as well as in the mine galleries. The SS-U water corresponds to surface water without mining influence and the SS-D water originates from the mixture of SS-U, G1 and G2 waters.
The water samples from Silveirinhos stream collected upstream from the galleries (SS-U) showed results within FAO’s permissible limits for irrigation purposes, with the exception of pH values, which were slightly below 6.5 due to the geological features of the catchment. Water samples from the mine drainage galleries G1 and G2 reflect the geochemical system of the coal seams and overlying strata. These mining effluents are not as acidic as one might expect, being neutral to slightly acidic, with pH values close to FAO’s inferior permissible limit.
Acidity in coal mine waters results mainly from the dissolution of oxidized pyritic materials associated with coal, during mining operations, which explains the existence of iron and sulphate in the water. The pH values of samples G1 and G2, concomitantly with their high levels of iron and sulphate, suggest the existence of an underground neutralization process. The scarcity of pyrites in particular layers and the predominance of carbonate minerals constitute the most common explanation for neutral or alkaline mine drainages [48]. However, the calcareous materials in this region are very scarce and do not constitute a reasonable justification. In São Pedro da Cova, a plausible origin of the neutralization process could be the mixture of acid groundwater, circulating in the shallow rock massif along the mine galleries and wells, with alkaline thermomineral water following deeper circulation paths, possibly through major faults. This hypothesis is corroborated by the relatively high fluoride content in G1 and G2 samples (a hydrogeochemical signature of thermomineral circulation), when compared to SS-U samples, and by the water temperature measured in situ during the nine sampling campaigns. Data analysis shows that the temperatures measured in mine effluents were higher than the average annual air temperature for the study area (15 °C), ranging from 9 °C in January to 20 °C in July. G1 and G2 temperatures remained constant during the study, with mean values of 18.9 ± 0.5 and 19.1 ± 0.4 °C, respectively, which are significantly higher than the mean water temperature in SS-U (15.4 ± 4.7 °C). SS-U water temperature is seasonal, and it closely follows changes in air temperature as can be seen in Figure 3. In SS-D the influence of the mining drainages is notorious, with a mean value of temperature of 17.3 ± 1.8 °C.
The number of particles in water can be expressed by turbidity. With the exception of the SS-U samples that recorded a turbidity median value of 1.1 NTU (Table 3), all the impacted waters have values that were far above the limits proposed by the US EPA of 2 NTU for directly consumed crops and unrestricted irrigation [49], and by Spain that recommends levels lower than 10 NTU for vegetable irrigation water [15].
High levels of turbidity can affect the performance of irrigation facilities, causing the clogging of the equipment, can lower the hydraulic conductivity of the soil, pollute the soil surface, and lead to aesthetic impairment of the water and of the vegetables produced. In addition, irrigating vegetables with turbid water could affect the quality of the products since microorganisms, such as parasites, bacteria and viruses can be attached to the solid particles and contaminate the crops [18].
The EC is also a significant parameter in determining the suitability of water for irrigation use, as it affects water salinity, which subsequently affects the productivity and yield of crops. The EC levels obtained in this study were all below the 3.0 dS/m permissible limit set by FAO for irrigation water [18]. However, the levels recorded in G1 and G2 were significantly higher than those recorded in Silveirinhos stream upstream from mine drainage galleries, these waters having a slight to moderate restriction on use.
Regarding major ions, the levels in SS-U were within FAO permissible limits, and significantly lower than in the impacted waters. Coal mining pollution significantly increases mineralization as a result of a greater water–rock interaction. Mg2+ and K+ were above the usual range for irrigation water, being the abundance of the ions as follows: SO42− > Mg2+ > Ca2+ > HCO3 > Na+ > Cl > K+ > CO32−.
The hydrogeochemical effect of coal mining, in terms of hydrogeochemical facies and major ion content, is illustrated by means of a Piper diagram (Figure 4) and a Stiff diagram (Figure 5). Surface water without mining influence (SS-U) has an intermediate SO4/Cl-Na/Mg classification while mine drainage waters (G1 and G2) as well as the water collected downstream from the mining effluents discharge (SS-D) have a hydrogeochemical SO4-Mg facies.
Sulfate is relatively common in water and has no major impact on soil other than contributing to the total salt content. Despite being within the usual range for irrigation (0–20 meq/L), concentrations of sulfate in G1 and G2 mine effluents, and in SS-D, were 26 and 17 times higher, respectively, than in SS-U, highlighting the mining influence. Chloride also contributes to the salinity of soils. It is necessary for plant growth in small amounts, but in high concentrations can inhibit plant growth or be toxic to some plants. In the studied waters, chloride levels were low, although G1 and G2 recorded two times the concentration of the upstream water samples.
Regarding the bicarbonate levels, G1 and G2 samples exhibited more than 25 times, and the SS-D samples more than 14 times, the concentration found in the SS-U samples, confirming the mining impact. High levels of bicarbonates can be directly toxic to some plant species. Levels greater than 1.5 meq/L are sufficient to cause concern. Concentrations of bicarbonates greater than 3.3 meq/L may pose a severe potential hazard. Bicarbonate reacts with calcium forming deposits of calcium carbonate and render calcium unavailable. Bicarbonate is also toxic to roots and reduces shoot growth, the uptake of phosphorus and of several micronutrients [13,15,18].
Calcium and magnesium are essential plant nutrients that occur naturally in water through the weathering of geological materials that contain these elements. However, in high concentrations, they are associated to soil aggregation and friability, being important qualitative criteria in the assessment of irrigation water quality. The average concentrations of magnesium in G1 and G2 mine waters (5.43 and 5.10 meq/L, respectively) were above FAO guidelines (0–5 meq/L). Calcium concentrations varied from 0.19 in SS-U to 4.12 meq/L in G2, far below the limit of 20 meq/L of FAO.
These results correlate with the TH values that were also calculated to categorize the water samples considering their calcium and magnesium contents, using the formula shown in Table 1. TH calculated mean values for G1, G2 and SS-D samples were 476, 461 and 308 mg/L, respectively, indicating very hard waters that can be considered harmful and unsuitable for irrigation use. In contrast, the SS-U samples are classified as soft, with a mean value of TH of 22 mg/L.
The RSC index was also calculated as it is an important parameter for irrigation used to indicate the alkalinity hazard for soil. RSC compares the relative concentrations of bicarbonate and carbonate ions with the concentrations of calcium and magnesium. The average RSC values were negative for all samples (<1.25 meq/L), indicating that according to this index waters are safe for irrigation.
Sodium and potassium also occur naturally in groundwater and surface water due to normal water–rock interaction. Among the soluble constituents of irrigation water, sodium is considered the most hazardous. High concentrations of sodium are undesirable because it adsorbs on to the soil cation exchange particles, causing deflocculation and pore sealing, decreasing soil permeability. In the studied water samples, sodium concentrations were low, with mean values varying from 0.32 in SS-U to 1.03 meq/L in G1 and G2. Otherwise, the mean concentrations of potassium in G1 and G2 were 5.2 and 6.7 mg/L, respectively, about three times higher than FAO’s limit of 2 mg/L [18]. In SS-D the mean concentration was 3.1 mg/L, also above the recommended limit.
Decades of research on the effect of irrigation water quality on soil permeability have established that the decreasing order of negative impacts of the four major cations follows the sequence: Na > K > Mg > Ca, although the current guidelines are still based on SAR and assume that potassium and magnesium pose no hazard. However, recent studies demonstrated that the negative effects of high K and Mg concentrations on soil permeability are substantial and that they should be taken into account through a new irrigation water quality parameter, the Cation Ratio of Structural Stability (CROSS) that can be directly incorporated into existing irrigation water quality guidelines by replacing SAR [50]. CROSS quantifies both the differing effects of Na and K as dispersing cations diminishing soil permeability and the differing effects of Mg and Ca as flocculating cations enhancing soil permeability. The interpretative guidelines for irrigation water quality involving SAR and CROSS are similar [51].
As mine waters are non-conventional irrigation waters and results revealed high levels of K and Mg, the CROSS ratio was also calculated by the following formula:
CROSS = Na + + 0.56   K +   Ca 2 + + 0.6 Mg 2 + / 2
where the concentrations of ions (Na, K, Ca, and Mg) are expressed in mmol/L. CROSS results were similar to SAR values, varying from 1.1 in SS-D to 0.8 in G1, G2 and SS-U. According to Richards (1954) [40] and FAO guidelines [18], all samples fall in the excellent section for irrigation, which is in accordance with RSC values. Concerning the permeability index (PI) results, only the SS-U samples were classified as of good quality, with more than 75% of maximum permeability. The impacted waters, although classified in Class II, presented values close to 25%, which can be considered unsuitable for irrigation (Table 3).
Thus, considering the characteristics discussed above, it can be concluded that the SS-U water is excellent for irrigation, but the mutual balancing of cations and anions leads to contradictory classifications of the aptness of G1, G2 and SS-D water samples for irrigation purposes. According to SAR, CROSS, and RSC values waters can be classified as good/excellent for irrigation, whilst according to TH and PI they are considered unsuitable.
Plant nutrient concentrations (nitrate, nitrite, ammonium and phosphate) and organic matter content (TOC and COD) were low in all water samples during the entire study period. Nitrate and phosphate are essential plant nutrients, but when in excessive amounts can cause water quality problems and accelerate eutrophication, altering the density and types of aquatic plants found in affected water bodies, promoting their degradation.
Boron is essential for the normal plant growth, but its occurrence in toxic concentrations makes it necessary to consider this element in assessing the water quality. Boron mean values in water samples ranged from 0.07 mg/L in SS-U to 0.26 mg/L in G2. Values are within acceptable threshold, not included in the restriction categories of the FAO classification.
Finally, several trace elements, mainly metals, were also analyzed as they are necessary for crop growth but when in high doses can cause serious environmental and health hazards. Their quantitative determination has shown that the waters affected by coal mining activities have higher metal content, especially iron, manganese, aluminum, nickel and arsenic, which may cause various health hazards such as cancer and environmental pollution. Only iron and manganese exceeded the FAO standards, with values far above the recommended concentrations, but the levels of aluminum almost duplicated in the impacted waters, and nickel and arsenic increased significantly in G1, G2 and SS-D, in a proportion of 7, 3, 3 times and 21, 48, 11 times, respectively, in comparison with the values in water samples collected upstream from the discharges.
Iron was the most abundant metal detected in the mine wastewater samples and in the samples collected downstream from the mine galleries. The median concentration of iron ranged from 0.12 mg/L in SS-U to 52.28 and 84.66 mg/L in G1 and G2, respectively, and 18.24 mg/L in SS-D. Iron can be a complex water quality problem, which not only affects plant growth, as it can compete with other needed micro-nutrients, but also can clog irrigation equipment. For micro-irrigation systems, iron levels need to be below 0.3 mg/L to prevent clogging. Above 1.0 mg/L, iron may lead to rust stains and discoloration on foliage plants in overhead irrigation systems, and above 5 mg/L iron is toxic to plant tissues.
Manganese presents many of the same issues as iron in irrigation water. It can clog irrigation equipment and cause foliar staining. The recommended drinking water standard for manganese is 0.05 mg/L, which is also the level where black staining and irrigation clogging may occur. Concentrations above 2.0 mg/L can be directly toxic to some plant species. In this study the mean concentration of manganese ranged from 0.06 mg/L in SS-U to 4.96 mg/L in G1.

4.2. Irrigation Water Quality Indices

As the results suggest that no unique parameter can sufficiently describe water quality, thus, the chemical status of the water samples was also assessed by using Water Pollution Indices. Indices were calculated covering a wide range of variables that were gathered in a single numerical value, allowing a simplified representation of a complex reality and also the evaluation of historical trends. The most significant parameters for the water quality evaluation were selected according to FAO’s guidelines, and for TETI according to the Toxicological Profiles of the Priority List of Hazardous Substances of ATSDR, in order to proceed with the calculation of the indices in a robust but simple way, using the smallest number of analytes. Indices calculations were also performed including other parameters, but no differences were found in the results and in the outcomes of the evaluation.
The results of the WQIA and the CI indices are shown in Table 4. From WQIA values waters were classified into five categories: excellent, good, poor, very poor and unsuitable, and CI values indicate a low, medium or high level of contamination.
The relation of precipitation events with the WQIA values is represented in Figure 6, with data highlighted in different colors according to the classification for irrigation purposes (Legend Table 4).
Results comparing upstream and downstream sampling sites point out the impact of mining effluents on surface water, being SS-U classified as excellent (WQIA values ranging from 1.4 to 3.2, with a mean value of 2.0) and SS-D as poor (WQIA values ranging from 16.2 to 105.4, with a mean value of 56.2). Analysis also revealed that G1 and G2 were, as expected, the two most polluted waters, reported as very poor, with values ranging from 66.2–129.8 and 69.4–140.1, respectively. Out of the 18 parameters considered for this study, iron and manganese were the two deciding parameters, followed by arsenic and EC, which exhibit the maximum influence (Qi x Wi) in the WQIA calculations.
Figure 6 shows fluctuations in WQIA values in both study periods: from April 2017 to December 2017 and from November 2018 to December 2019. In the first period, the effect of draught in water quality is clear: the highest WQIA values correspond to the driest months due to the lack of mixture of mine drainage with recently infiltrated precipitation. In the second period, a similar trend is observed in September 2019, with G1, G2 and SS-D presenting similar values. It was not possible to collect samples on SS-U because there was no streamflow due the drought conditions. The lower WQIA values from this period were observed in February and in December 2019 as a result of a dilution effect due the infiltration of precipitation in the previous months.
The results of the computed CI index are, in general, comparable with the WQIA values. The CI results for impacted waters exceeded the value of 3, ranging from a mean value of 10.1 in SS-D to 29.4 in G1 and 35.8 in G2, which indicates a high degree of pollution due mainly to iron, manganese and bicarbonate content, which exceed the limits of FAO guidelines. The SS-U samples have their computed CI values below 1, reflecting the absence of coal mining influence.
TETI results, based on the elemental toxicological impact, are shown in Table 5.
TETI indicates that manganese and potassium had the highest impact on the toxicological profiles of the polluted waters (G1, G2 and SS-D) with a total score above 2000, according to the ATSDR assessment, followed by ammonium, nitrate and aluminum [43]. Regarding SS-U water samples, nitrate was the most important constituent in terms of water quality concerns, followed by potassium, aluminum and zinc. G1 had the highest index value followed by G2 and SS-D, which clearly demonstrates the impact of mine activities on the water environment, where higher TETI values represent lower quality.
The evaluation of these three selected indices (WQIA, CI and TETI) highlights the coal mine inputs of metals and other pollutants in the study area. Although the three indices specify similar levels of contamination, their outcomes regarding the most important constituents are not uniform.
For example, for TETI calculation, potassium had a high impact with a total score 5–11 times higher than in SS-U. However, it is not considered in FAO or any international guideline for restriction on use purposes, hence it is not accounted for either the WQIA or CI indices. The same occurs regarding ammonia.
These findings clearly highlight the limitations of each index and of the international water quality guidelines that are, firstly, non-standardized between different countries and, secondly, do not provide guidelines for a number of pollutants of importance for specific matrices as is the case of coal mine effluents.

4.3. PAH Analyses

Regarding PAHs, analyses were performed in six campaigns (April 2017, September 2017, December 2017, November 2018, February 2019 and May 2019). According to the ring numbers, PAHs can be classified into three classes: 2–3 rings, 4 rings, and 5–6 rings composition, which represent low, medium, and high molecular weight hydrocarbons, respectively.
PAHs were detected in water samples at very low concentrations, with prevalence of low molecular weight compounds, with average concentration percentages that varied from 62% to 80%, which indicate a petrogenic origin consistent with the water circulation through the coal bearing rocks in the mine [53] (Figure 7).
The highest average concentration was detected in G1 (42.5 ng/L). SS-U and SS-D showed similar values of 35.8 and 34.1 ng/L, respectively.
The carcinogenic PAHs (BaA, Chr, BbF, BkF, BaP, InP, and DahA) were detected mostly in G1 and G2, accounting for 28.8% and 27.6% of the total average concentration of 16-PAHs, in contrast to 5.6% and 6.7% found in SS-U and SS-D, being the main components in mine effluents the BaA and the BkF. The carcinogenic PAHs are high-ring PAHs, i.e., 4–6 ring, which is related to the degree of coal maturation [47].
PAHs are of high environmental and human health concern, as they are toxic and persistent in the environment, susceptible to long-range atmospheric transport, and able to bioaccumulate. The uptake of PAHs by plants is important when considering their transfer from soils and water into the food chain. Recent studies also demonstrate that PAHs-metal co-contamination also alters PAHs uptake, attributing to the metal–soil or metal–root interactions [54].
As food chain is the most important pathway for pollutants’ entry into the human body, the uptake of carcinogenic PAHs and heavy metals through soil-to-root system and their translocation/accumulation in plant tissues is very important, particularly for food crops cultivated on non-treated wastewater-irrigated soils, as is the case.

5. Conclusions

This study points out the suitability of the coal mining effluents and the polluted surface water for irrigation purposes in São Pedro da Cova abandoned coalfield. The results allow proposing a cost-effective assessment methodology adjusted to specific problems of the water, minimize pollution of natural watercourses and soils and increase the potential of use of these effluents.
The evaluation of the use of mining effluents for irrigation can be an issue as specific water quality guidelines or legislation does not exist, and several dangerous chemicals are not included in routine water quality assays. In an attempt to standardize decision-making regarding irrigation with mining effluents, the criteria and data have been combined in user-friendly indices, which could assist in the practical implementation of mining effluents irrigation control plans as part of optimal mine-water management and reuse strategy. Eighteen parameters including selected anions, cations and trace elements were chosen after an exhaustive analysis and were used for water quality indices calculation. Source specific water quality indices are a very helpful tool to represent water quality in a simple and understandable manner, minimizing the data volume to a great extent and simplifying the expression of water quality status, giving efficiently the overall water quality of a specific area and for a specific use.
Results revealed that samples associated with mining activities have unacceptable index values due mainly to the high concentration of iron, manganese, bicarbonates, magnesium and potassium, and are not suitable for irrigation. WQIA indicated an overall Poor/Very Poor quality status, which is in accordance with the level of contamination (CI) of these waters. The TETI index, which only reflects the elemental toxicological impact of the waters, not considering other fundamental water quality parameters, indicated that the impacted water samples had higher index values, and manganese had the highest impact on the toxicological profiles, also showing the influence of coal activities on surface water quality. Risks to human health also arise from water pollution by organic substances such as PAHs with several carcinogenic compounds detected in G1 and G2. Levels in the water collected upstream from mine drainage points are within acceptable range for irrigation.
Scientific community and local authorities of mine areas are committed to mitigating the effects of past actions through the development of better management strategies for reducing environmental and health impacts in the mining area. The understanding of the mine water chemistry is fundamental for the design of an effective treatment system. Simple passive treatments could be an option for the contaminated drainage at this abandoned mine site according to Robert Hedin el al. 2013 (Effective Passive Treatment of Coal Mine Drainage. Robert Hedin, Ted Weaver, Neil Wolf, George Watzalaf. Paper presented at the 35th Annual National Association of Abandoned Mine Land Programs Conference, 2013). If properly designed, constructed and maintained, passive systems provide highly reliable treatment at a fraction of the cost of active alternatives.
As the overall quality of groundwater from mine drainage galleries revealed contamination, it is very important to raise awareness for rapid intervention in the area, since the mine is located near a population center and social infrastructures, as well as to mitigate the pollution on adjacent agricultural lands. Despite the evident deleterious impact for local communities, the health effects of mine drainage remain neglected in research and policy arenas, mainly because of the lack of documented evidence. Such communities suffer the effects of mine drainage principally through a perpetual risk posed by water pollution. This project intended to study the impact of mine drainage contamination towards an investment in temporary to long-term solutions, in order to reduce the risks caused by mining externalities.

Author Contributions

Conceptualization, C.M.; methodology, C.M., A.M., J.E.M.; validation, C.M., J.E.M.; formal analysis, C.M., A.M., J.E.M.; investigation, C.M., J.E.M., V.M., P.S., J.R., J.R.R., A.M. and D.F.; resources, J.E.M., C.M., D.F.; writing—original draft preparation, C.M., J.E.M., A.M., P.S., J.R.; writing—review and editing, C.M., J.E.M., A.M.; supervision, C.M., J.E.M., D.F.; project administration, D.F.; funding acquisition, D.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the project CoalMine—Coal mining wastes: assessment, monitoring and reclamation of environmental impacts through remote sensing and geostatistical analysis—financed by the Portuguese Science and Technology Foundation, FCT, call AAC nº 02/SAICT/2017 (POCI-01-0145-FEDER-030138) and framed within the activities of the ICT (Ref. UIDB/04683/2020).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data supporting the findings of this study are available within the article.

Acknowledgments

The authors acknowledge the Mining Museum of São Pedro da Cova for the permission to use the photograph from Figure 1a. This work received support and help from the UID/QUI/50006/2020 with funding from FCT/MCTES through national funds.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mahato, M.K.; Singh, A.K.; Singh, G.; Mishra, L.P. Impacts of coal mine water and Damodar River water irrigation on soil and maize (Zea mays L.) in a coalfield area of Damodar Valley, India. Arch. Agric. Environ. Sci. 2017, 2, 293–297. [Google Scholar] [CrossRef]
  2. Bott, T.L.; Jackson, J.K.; McTammany, M.E.; Newbold, J.D.; Rier, S.T.; Sweeney, B.W.; Battle, J.M. Abandoned coal mine drainage and its remediation: Impacts on stream ecosystem structure and function. Ecol. Appl. 2012, 22, 2144–2163. [Google Scholar] [CrossRef] [PubMed]
  3. Johnston, D.; Potter, H.; Jones, C.; Rolley, S.; Watson, I.; Pritchard, J. Abandoned Mines and the Water Environment; Environment Agency: Bristol, UK, 2008; p. 31. [Google Scholar]
  4. Tiwary, R.K. Environmental Impact of Coal Mining on Water Regime and Its Management. Water Air Soil Pollut. 2001, 132, 185–199. [Google Scholar] [CrossRef]
  5. Lemos de Sousa, M.J.; Wagner, R.H. General Description of the Terrestrial Carboniferous Basins in Portugal and History of Investigations. In The Carboniferous of Portugal. Memórias dos Serviços Geológicos de Portugal; Lemos de Sousa, M.J., Oliveira, J.T., Eds.; Direção Geral de Geologia e Minas: Lisboa, Portugal, 1983; Volume 29, pp. 117–126. [Google Scholar]
  6. Wagner, R.H.; Lemos de Sousa, M.J. The Carboniferous Megafloras of Portugal—A revision of identifications and discussion of stratigraphic ages. In The Carboniferous of Portugal. Memórias dos Serviços Geológicos de Portugal; Lemos de Sousa, M.J., Oliveira, J.T., Eds.; Direção Geral de Geologia e Minas: Lisboa, Portugal, 1983; Volume 29, pp. 127–152. [Google Scholar]
  7. Eagar, R.M.C. The non marine bivalve fauna of the Stephanian C of North Portugal. In The Carboniferous of Portugal. Memórias dos Serviços Geológicos de Portugal; Sousa, M.J.L., Oliveira, J.T., Eds.; Direção Geral de Geologia e Minas: Lisboa, Portugal, 1983; Volume 29, pp. 179–185. [Google Scholar]
  8. Fernandes, J.P.; Pinto de Jesus, A.; Teixeira, F.; Sousa, M.J.L. First palynological results in the Douro Carboniferous Basin (NW of Portugal). In Proceedings of XIII Jornadas de Paleontologia “Fósiles de Galicia” y V Reunión Internacional Proyecto 351 PICG “Paleozoico Inferior del Noroeste de Gondowana”; Sociedade Espanola de Paleontologia: Coruna, Spain, 1997; pp. 176–179. [Google Scholar]
  9. Custódio, J. Museum of Coal & Mines of Pejão; Museum Programme: Castelo de Paiva, Portugal, 2004. [Google Scholar]
  10. Pinto de Jesus, A. Sedimentary and tectonic evolution of the Douro Coalfield Basin (lower Stephanian C, NW Portugal). Cad. Lab. Xeolóxico Laxe 2003, 28, 107–125. [Google Scholar]
  11. HLPE. Water for food security and nutrition. In A Report by the High Level Panel of Experts on Food Security and Nutrition of the Committee on World Food Security; Food and Agriculture Organization of the United Nations: Rome, Italy, 2015. [Google Scholar]
  12. Pescod, M.B. Wastewater Treatment and Use in Agriculture—FAO Irrigation and Drainage Paper 47; FAO: Rome, Italy, 1992. [Google Scholar]
  13. Al-Hwaiti, M.S.; Brumsack, H.J.; Schnetger, B. Suitability assessment of phosphate mine waste water for agricultural irrigation: An example from Eshidiya Mines, South Jordan. Environ. Earth Sci. 2016, 75, 276. [Google Scholar] [CrossRef]
  14. Müller, K.; Cornel, P. Setting water quality criteria for agricultural water reuse purposes. J. Water Reuse Desalin. 2016, 7, 121–135. [Google Scholar] [CrossRef] [Green Version]
  15. Jeong, H.; Kim, H.; Jang, T. Irrigation Water Quality Standards for Indirect Wastewater Reuse in Agriculture: A Contribution toward Sustainable Wastewater Reuse in South Korea. Water 2016, 8, 169. [Google Scholar] [CrossRef] [Green Version]
  16. Alobaidy, A.; Al-Sameraiy, M.; Kadhem, A.; Majeed, A. Evaluation of treated municipal wastewater quality for irrigation. J. Environ. Prot. 2010, 1, 216–225. [Google Scholar] [CrossRef] [Green Version]
  17. Abdul-Rahman, S.; Saoud, I.P.; Owaied, M.K.; Holail, H.; Farajalla, N.; Haidar, M.; Ghanawi, J. Improving Water Use Efficiency in Semi-Arid Regions through Integrated Aquaculture/Agriculture. J. Appl. Aquac. 2011, 23, 212–230. [Google Scholar] [CrossRef]
  18. Ayers, R.S.; Westcot, D.W. Water quality for agriculture. In FAO Irrigation and Drainage Paper Rev. 1; FAO: Rome, Italy, 1994. [Google Scholar]
  19. Alves, M.T.R.; Teresa, F.B.; Nabout, J.C. A global scientific literature of research on water quality indices: Trends, biases and future directions. Acta Limnol. Bras. 2014, 26, 245–253. [Google Scholar] [CrossRef] [Green Version]
  20. Kachroud, M.; Trolard, F.; Kefi, M.; Jebari, S.; Bourrié, G. Water Quality Indices: Challenges and Application Limits in the Literature. Water 2019, 11, 361. [Google Scholar] [CrossRef] [Green Version]
  21. Poonam, T.; Tanushree, B.; Sukalyan, C. Water Quality Indices- Important Tools for Water Quality Assessment: A Review. Int. J. Adv. Chem. 2015, 1, 15–29. [Google Scholar]
  22. Horton, R.K. An index number system for rating water quality. J. Water Pollut. Control Fed. 1965, 37, 300–306. [Google Scholar]
  23. Abbasi, T.; Abbasi, S. Water Quality Indices; Elsevier: Amsterdam, The Netherlands, 2012. [Google Scholar]
  24. Misaghi, F.; Delgosha, F.; Razzaghmanesh, M.; Myers, B. Introducing a water quality index for assessing water for irrigation purposes: A case study of the Ghezel Ozan River. Sci. Total Environ. 2017, 589, 107–116. [Google Scholar] [CrossRef]
  25. Cude, C.G. Oregon water quality index a tool for evaluating water quality management effectiveness. J. Am. Water Resour. Assoc. 2001, 37, 125–137. [Google Scholar] [CrossRef]
  26. Backman, B.; Bodiš, D.; Lahermo, P.; Rapant, S.; Tarvainen, T. Application of a groundwater contamination index in Finland and Slovakia. Environ. Geol. 1998, 36, 55–64. [Google Scholar] [CrossRef]
  27. Biswas, P.K.; Uddin, N.; Alam, S.; -Us-Sakib, T.; Sultana, S.; Ahmed, T. Evaluation of Heavy Metal Pollution Indices in Irrigation and Drinking Water Systems of Barapukuria Coal Mine Area, Bangladesh. Am. J. Water Res. 2017, 5, 146–151. [Google Scholar]
  28. Ali, A.; Strezov, V.; Davies, P.; Wright, I. Environmental impact of coal mining and coal seam gas production on surface water quality in the Sydney basin, Australia. Environ. Monit. Assess. 2017, 189, 408. [Google Scholar] [CrossRef]
  29. EFSA. Polycyclic Aromatic Hydrocarbons in Food—Scientific Opinion of the Panel on Contaminants in the Food Chain. EFSA J. 2008, 6, 724. [Google Scholar] [CrossRef]
  30. Rodríguez-Eugenio, N.; McLaughlin, M.; Pennock, D. Soil Pollution, A Hidden Reality; FAO: Rome, Italy, 2018. [Google Scholar]
  31. INE. Censos 2011 Resultados Definitivos—Portugal; Instituto Nacional de Estatística, I.P.: Lisbon, Portugal, 2012. [Google Scholar]
  32. Julivert, M.; Fontboté, J.M.; Ribeiro, A.; Conde, L. Mapa Tectónico de la Península Ibérica y Baleares; IGME—Instituto Geológico y Minero de España: Madrid, Spain, 1974. [Google Scholar]
  33. Medeiros, A.; Pereira, E.; Moreira, A. Notícia Explicativa da Folha 9-D (Penafiel) da Carta Geológica de Portugal à Escala 1:50 000; Serviços Geológicos de Portugal: Lisbon, Portugal, 1980. [Google Scholar]
  34. Peel, M.C.; Finlayson, B.L.; McMahon, T.A. Updated world map of the Köppen-Geiger climate classification. Hydrol. Earth Syst. Sci. 2007, 11, 1633–1644. [Google Scholar] [CrossRef] [Green Version]
  35. AEMET-IM. Atlas Climático Ibérico. Temperatura do Ar e Precipitação (1971–2000); Departamento de Producción da Agência Estatal de Meteorologia de Espanha (Área de Climatología y Aplicaciones Operativas) e Departamento de Meteorologia e Clima (Divisão de Observação Meteorológica e Clima), do Instituto de Meteorologia: Lisboa, Portugal, 2011. [Google Scholar]
  36. ISO. ISO 5667-3:2018—Water Quality—Sampling—Part 3: Preservation and Handling of Water Samples, ISO: Geneva, Switzerland, 2018.
  37. Rice, E.W.; Baird, R.B.; Eaton, A.D. Standard Methods for the Examination of Water and Wastewater, 23rd ed.; American Public Health Association: Washington, DC, USA, 2017. [Google Scholar]
  38. Rodier, J.; Legube, B.; Merlet, N. L’analyse de l’eau, 10th ed.; Dunod: Paris, France, 2016. [Google Scholar]
  39. Borges, B.; Armindo, M.; Ferreira, I.M.P.L.V.O.; Mansilha, C. Dispersive liquid–liquid microextraction for the simultaneous determination of parent and nitrated polycyclic aromatic hydrocarbons in water samples. Acta Chromatogr. 2018, 30, 119–126. [Google Scholar] [CrossRef]
  40. Richards, L.A. Diagnosis and Improvement of Saline and Alkali Soils. Agricultural Handbook 60; United States Department of Agriculture: Washington, DC, USA, 1954. [Google Scholar]
  41. Sawyer, C.N. Chemistry of Sanitary Engineers, 2nd ed.; McGraw-Hill: New York, NY, USA, 1967. [Google Scholar]
  42. Doneen, L.D. Notes on Water Quality in Agriculture; Department of Water Science and Engineering, University of California, Davis: Oakland, CA, USA, 1964. [Google Scholar]
  43. Oni, O.; Fasakin, O. The use of water quality index method to determine the potability of surface water and groundwater in the vicinity of a municipal solid waste dumpsite in Nigeria. Am. J. Eng. Res. 2016, 5, 96–101. [Google Scholar]
  44. Tyagi, S.; Sharma, B.; Singh, P.; Dobhal, R. Water Quality Assessment in Terms of Water Quality Index. Am. J. Water Res. 2013, 1, 34–38. [Google Scholar] [CrossRef]
  45. Chatterjee, C.; Raziuddin, M. Determination of water quality index (WQI) of a degraded river in Asansol Industrial area, P.O. Raniganj, District Burdwan, West Bengal. Nat. Environ. Pollut. Technol. 2002, 1, 181–189. [Google Scholar]
  46. ATSDR. Priority List of Hazardous Substances; Agency for Toxic Substances and Disease Registry: Division of Toxicology and Human Health Sciences: Atlanta, GA, USA, 2019. [Google Scholar]
  47. Chen, D.; Feng, Q.; Liang, H.; Gao, B.; Alam, E. Distribution characteristics and ecological risk assessment of polycyclic aromatic hydrocarbons (PAHs) in underground coal mining environment of Xuzhou. Hum. Ecol. Risk Assess. 2019, 25, 1564–1578. [Google Scholar] [CrossRef]
  48. Singh, G. Impact of coal mining on mine water quality. Int. J. Mine Water 1988, 7, 49–59. [Google Scholar] [CrossRef]
  49. EPA. Guidelines for Water Reuse; U.S. Environmental Protection Agency: Washington, DC, USA, 2012. [Google Scholar]
  50. Rengasamy, P.; Marchuk, A. Cation ratio of soil structural stability (CROSS). Soil Res. 2011, 49, 280–285. [Google Scholar] [CrossRef]
  51. Oster, J.D.; Sposito, G.; Smith, C.J. Accounting for potassium and magnesium in irrigation water quality assessment. Calif. Agric. 2016, 70, 71–76. [Google Scholar] [CrossRef] [Green Version]
  52. Climatological Bulletin. Available online: https://www.ipma.pt/pt/publicacoes/boletins.jsp (accessed on 6 April 2021).
  53. Stout, S.A.; Emsbo-Mattingly, S.D. Concentration and character of PAHs and other hydrocarbons in coals of varying rank—Implications for environmental studies of soils and sediments containing particulate coal. Org. Geochem. 2008, 39, 801–819. [Google Scholar] [CrossRef]
  54. Zhang, S.; Yao, H.; Lu, Y.; Yu, X.; Wang, J.; Sun, S.; Liu, M.; Li, D.; Li, Y.-F.; Zhang, D. Uptake and translocation of polycyclic aromatic hydrocarbons (PAHs) and heavy metals by maize from soil irrigated with wastewater. Sci. Rep. 2017, 7, 12165. [Google Scholar] [CrossRef] [Green Version]
Figure 1. São Pedro da Cova location and geological setting (modified from [10]).
Figure 1. São Pedro da Cova location and geological setting (modified from [10]).
Water 13 02157 g001
Figure 2. Some aspects of the study area: (a) São Pedro da Cova mine in the mid 20th century; (b) São Pedro da Cova mine in 2021; (c) Silveirinhos stream upstream from the mining effluents discharge; (d) mining effluents discharge; (e) water monitoring in Silveirinhos stream downstream from the mining effluents discharge; (f) irrigation with polluted water from Silveirinhos stream; (g) agricultural area irrigated with polluted water from Silveirinhos stream.
Figure 2. Some aspects of the study area: (a) São Pedro da Cova mine in the mid 20th century; (b) São Pedro da Cova mine in 2021; (c) Silveirinhos stream upstream from the mining effluents discharge; (d) mining effluents discharge; (e) water monitoring in Silveirinhos stream downstream from the mining effluents discharge; (f) irrigation with polluted water from Silveirinhos stream; (g) agricultural area irrigated with polluted water from Silveirinhos stream.
Water 13 02157 g002
Figure 3. Average monthly air temperature measured between 1981 and 2010 in Serra do Pilar weather station (IPMA). Water temperature recorded at the sampling points SS-U, G1, G2 and SS-D on April, June, October and December 2017, January and November 2018, and February, May, September and December 2019.
Figure 3. Average monthly air temperature measured between 1981 and 2010 in Serra do Pilar weather station (IPMA). Water temperature recorded at the sampling points SS-U, G1, G2 and SS-D on April, June, October and December 2017, January and November 2018, and February, May, September and December 2019.
Water 13 02157 g003
Figure 4. Piper diagram of the studied waters (average values from April 2017 to December 2019, n = 9).
Figure 4. Piper diagram of the studied waters (average values from April 2017 to December 2019, n = 9).
Water 13 02157 g004
Figure 5. Stiff diagram of the studied waters (average values from April 2017 to December 2019, n = 9).
Figure 5. Stiff diagram of the studied waters (average values from April 2017 to December 2019, n = 9).
Water 13 02157 g005
Figure 6. Relation between precipitation and water quality assessed by means of WQIA. Precipitation measured at Porto meteorological station. Data from the Portuguese Institute for Sea and Atmosphere, I. P. (IPMA, IP) [52].
Figure 6. Relation between precipitation and water quality assessed by means of WQIA. Precipitation measured at Porto meteorological station. Data from the Portuguese Institute for Sea and Atmosphere, I. P. (IPMA, IP) [52].
Water 13 02157 g006
Figure 7. Composite model and average concentrations of PAHs in the studied samples.
Figure 7. Composite model and average concentrations of PAHs in the studied samples.
Water 13 02157 g007
Table 1. Water quality classification as per different water quality ratios/parameters.
Table 1. Water quality classification as per different water quality ratios/parameters.
ParameterCategoriesRanges
Sodium Adsorption Ratio (SAR) [40] SAR = Na + Ca 2 + +   Mg 2 + / 2 Excellent
Good
Doubtful
Unsuitable
<10
10–18
18–26
>26
Total Hardness (TH) [41] TH = Ca 2 + + Mg 2 + × 50 Soft
Moderately hard
Hard
Very hard
<75
75–150
150–300
>300
Residual Sodium Carbonate (RSC) [40] RSC = CO 3 2 + HCO 3 Ca 2 + + Mg 2 + Good
Medium
Not suitable
<1.25
1.25–2.5
> 2.5
Permeability index (PI) [42] PI = Na + + HCO 3 Ca 2 + +   Na + + Mg 2 + × 100 Class I
Class II
Class III
>75
25–75
<25
Table 2. Standards for irrigation water and relative weight of parameters.
Table 2. Standards for irrigation water and relative weight of parameters.
No.Parameter 1Standards 2
Si
1/SiKRelative Weight Wi
1pH6.5–8.50.11760.0070.00084
2Elect.Conduct.(EC)0.71.4290.01019
3Sodium690.0140.00010
4Chloride141.60.0070.00005
5Boron0.71.4290.01019
6Nitrate22.140.0450.00032
7Bicarbonate91.460.0110.00008
8Aluminum5.00.20.00143
9Arsenic0.1010.00.07135
10Cadmium0.011000.71350
11Chromium0.1010.00.07135
12Copper0.205.00.03568
13Fluoride1.01.00.00714
14Iron5.00.20.00143
15Manganese0.205.00.03568
16Nickel0.205.00.03568
17Lead5.00.20.00143
18Zinc2.00.50.00357
1 All values are in mg/L, except pH and EC (mS∙cm−1); 2 FAO [18]; source: own elaboration.
Table 3. Summary statistics of physical and chemical parameters analyzed in mining effluents and in Silveirinhos stream. Comparison with FAO Guidelines for irrigation.
Table 3. Summary statistics of physical and chemical parameters analyzed in mining effluents and in Silveirinhos stream. Comparison with FAO Guidelines for irrigation.
ParametersFAO Guidelines for Irrigation [18]Silveirinhos StreamMine
Wastewater
Mine
Wastewater
Silveirinhos Streamp
Potential Irrigation ProblemUnitsDegree of Restriction on Use(SS-U)(G1) (G2) (SS-D)
NoneSlight to ModerateSevere
Salinity (affects crop water availability)
ECdS/m<0.70.7–3.0>3.00.08 ± 0.02 c0.84 ± 0.14 a0.92 ± 0.18 a0.56 ± 0.21 b<0.0001 ##
Specific Ion toxicity (affects sensitive crops)
Sodium (Na)meq/L<33–9>90.32 ± 0.04 b1.03 ± 0.08 a1.03 ± 0.06 a0.78 ± 0.22 a<0.0001 ##
Chloride (Cl)meq/L<44–10>100.39 ± 0.030.87 ± 0.180.90 ± 0.120.72 ± 0.170.0003 &
Boron (B)mg/L<0.70.7–3.0>3.00.07 ± 0.02 a0.24 ± 0.07 bc0.26 ± 0.08 c0.14 ± 0.05 ab0.00127 #
Trace Elements (* Recommended Maximum Concentration)
Aluminum (Al)mg/L5.0 *0.19 ± 0.05 a0.37 ± 0.06 b0.44 ± 0.09 c0.21 ± 0.19 ab0.0359 #
Arsenic (As)mg/L0.10 *0.001 ± 0.001 d0.015 ± 0.005 a0.034 ± 0.009 b0.008 ± 0.003 c<0.0001 ##
Cadmium (Cd)mg/L0.01 *<LD<LD<LD<LD-
Chromium (Cr)mg/L0.10 *<LD<LD<LD<LD-
Copper Cu)mg/L0.20 *<LD<LD<LD<LD-
Fluoride (F)mg/L1.0 *0.04(0.00–0.07) b0.29(0.00–0.38) a0.33(0.16–0.35) a0.15(0.09–0.27) ab0.0001 &
Iron (Fe)mg/L5.0 *0.12(0.03–0.47) c52.28(45.66–62.39) ab84.66(57.65–115.08) a18.24(15.72–62.79) b<0.0001 &
Manganese (Mn)mg/L0.20 *0.06 ± 0.02 b4.96 ± 1.08 a4.82 ± 1.16 a2.93 ± 1.72 a<0.0001 ##
Nickel (Ni)mg/L0.20 *<LD b0.04(0.00–0.05) a0.02(0.00–0.03) ab0.02(0.00–0.03) ab0.0008 &
Lead (Pb)mg/L5.0 *<LD<LD<LD<LD-
Zinc (Zn)mg/L2.0 *0.01(0.00–0.90)0.00(0.00–0.08)0.00(0.00–0.03)0.00(0.00–0.03)n.s.
Miscellaneous Effects (affects susceptible crops)
Nitrate (NO3-N) 1mg/L<55–30>300.42(0.23–2.20)0.28(0.00–1.62)0.28(0.00–0.67)0.61(0.05–0.91)n.s.
Bicarbonates (HCO3)meq/L<1.51.5–8.5>8.50.12 ± 0.03 c2.99 ± 0.70 ab3.37 ± 1.04 a1.71 ± 1.10 b<0.0001 ##
Carbonates (CO3)meq/L <LD<LD<LD<LD-
pH Normal range 6.5–8.46.26(6.17–6.56) ab6.16(6.09–6.44) a6.12(6.00–6.57) a6.51(6.23–6.70) b0.0005 &
TurbidityNTU 1.1(0.3–4.9) b74.2(9.0–328.5) a111.5(1.9–342.0) a32.5(5.3–109.6) ab0.0008 &
Cyanidemg/L <LD<LD<LD<LD-
Salinity: Other cations and anionsUsual range
in irrigation water
Calcium (Ca)meq/L 0–20 0.19 ± 0.07 b4.08 ± 0.62 a4.12 ± 1.12 a2.71 ± 1.30 a<0.0001 ##
Magnesium (Mg)meq/L 0–5 0.25 ± 0.05 c5.43 ± 0.82 a5.10 ± 1.39 ab3.46 ± 1.54 b<0.0001 ##
Sulfate (SO4)meq/L 0–20 0.13(0.09–0.20) b3.56(2.96–3.98) a4.24(3.03–5.03) a1.79(1.16–4.08) ab<0.0001 &
Nutrients
Nitrite (NO2-N) 1mg/L <LD<LD<LD<LD-
Ammonium (NH4-N) 1mg/L 0–5 0.01(0.00–0.02) b0.60(0.32–1.80) a0.59(0.26–2.60) a0.52(0.13–2.06) a0.0006 &
Phosphate (PO4-P) 1mg/L 0–2 0.01(0.00–0.40)0.05(0.00–0.30)0.01(0.00–0.33)0.00(0.00–0.10)n.s.
Potassium (K)mg/L 0–2 0.51(0.00–1.17) b4.10(3.06–16.40) a5.14(3.32–19.40) a2.08(0.00–13.11) ab0.0009 &
TOC mg/L 0.80(0.23–1.55)0.56(0.00–33.16)0.40(0.00–14.47)0.63(0.60–28.56)n.s.
CODmg/L 5.65(0.10–7.42)10.20(0.00–14.70)11.60(0.00–19.10)7.96(1.20–11.00)n.s.
Suitability for irrigation water and soil water infiltration. Color Classification according to guidelines in Table 1
SARmeq/L 0.680.470.480.45-
CROSSmmol/L 1.120.820.840.76-
THmg/L 22.0 475461308-
RSCmeq/L −0.3−6.5−5.8−4.5-
PI% 87.026.028.030.0-
1 Results reported in terms of elemental nitrogen/phosphorus. EC—electrical conductivity in deciSiemens per metre at 25 °C; TOC—Total Organic Carbon; COD—Chemical Oxygen Demand; SAR—sodium adsorption ratio; CROSS—Cation Ratio of Structural Stability; TH—total hardness; RSC—residual sodium carbonate; PI—permeability index. Different letters (a, b, c) for each parameter in a row show statistically significant differences (p < 0.05) between means in normal distribution and median in non-normal distribution. Data expressed as mean ± standard deviation or as median (minimum−maximum). n.s, not significant. #—p-values from one-way ANOVA analysis. Means were compared by Tukey’s test, since homogeneity of variances was confirmed by Levene’s test (p > 0.05). ##—p-values from one-way Welch ANOVA analysis. Means were compared by Tamhane’s test, since homogeneity of variances was not confirmed by Levene’s test (p < 0.05). &—p-values from Kruskal–Wallis analysis. Medians were compared by Dunn’s post-hoc test. Analytical methods limits of detection (LD): 1.0 mg/L for COD; 0.3 for TOC; 0.8µS/cm for EC; 0.2 NTU for Turbidity; 0.02 mg/L for F; 0.005 meq/L for HCO3 and CO3; 0.03 mg/L for Cl, SO4, NO3, Na, K, Ca, Mg; 0.01 mg/L for PO4 and Cyanide; 0.001 mg/L for NO2 and NH4; 0.001 mg/L for Cd, As, Pb, Ni, Cu, Zn, Al and Cr; 0.025 mg/L for Fe; 0.05 for Al. For an easier analysis of the table the FAO‘s Guidelines are highlighted in gray.
Table 4. Calculated values of WQIA (Water Quality Status) and CI (Level of Contamination).
Table 4. Calculated values of WQIA (Water Quality Status) and CI (Level of Contamination).
Sampling
Date
WQIACI Degree of Contamination
SS-UG1G2SS-DSS-UG1G2SS-D
apr172.8 99.3 85.3 47.3−13.3 32.5 32.4 6.4
jun17 3.2 110.8 107.3 94.0−12.5 34.5 40.0 24.0
sep17-129.8140.1105.4- 41.2 55.3 30.1
dec172.0100.3101.1 47.5−13.2 28.9 40.8 3.9
nov181.7 96.1101.0 40.9−15.5 30.3 35.5 10.3
feb191.4 75.4 73.2 32.5−15.3 18.5 27.0−1.2
may191.7 81.0 81.8 35.6−15.8 26.0 30.4−1.7
sep19- 87.9 86.3 86.1- 31.6 38.8 26.5
dec19 1.5 66.2 69.416.2−16.2 20.7 22.0−7.0
Avg 2.0 94.1 93.9 56.2 −14.5 29.4 35.8 10.1
Water Quality Status (WQIA): Excellent (0–25 ); Good (26–50 ); Poor (51–75 ); Very Poor (76–100 ); Unsuitable (>100 ). Level of contamination (CI): Low (<1 ); medium (1–3 ); high (>3 ).
Table 5. Trace element toxicity index values.
Table 5. Trace element toxicity index values.
TsiSS-UG1G2SS-D
CiCi × TsiCiCi × TsiCiCi × TsiCiCi × Tsi
K (Potassium) mg/L6070.61371.145.163134.486.744089.833.081867.54
NO3 (Nitrate) mg/L6050.71428.690.46276.960.26154.610.56336.11
NO2 (Nitrite) mg/L6100.0032.090.015.690.017.590.029.83
NH4 (Ammonia) mg/L7420.013.920.76560.370.84625.880.73543.56
PO4 (Phosphate) mg/L2640.1540.430.1538.340.2772.190.1232.65
CN (Cyanide) mg/L10700.0010.790.001.120.001.080.001.11
B (Boron) mg/L4400.0730.250.24103.490.26112.550.1461.95
Al (Aluminum) mg/L6850.19129.070.37250.390.44300.780.21142.85
As (Arsenic) mg/L16760.0011.200.0225.260.0357.510.0113.65
Cd (Cadmium) mg/L13180.0000.140.000.030.000.010.000.07
Cr (Chromium) mg/L8930.0011.050.000.720.000.570.000.45
Cu (Copper) mg/L8050.0032.190.000.670.001.440.000.55
F (Fluoride) mg/L5500.0318.540.26141.110.28156.200.1584.52
Mn (Manganese) mg/L7970.0645.084.963956.934.823845.082.932338.42
Ni (Nickel) mg/L9930.015.380.0435.000.0215.460.0216.85
Pb (Lead) mg/L15310.0035.310.001.390.001.160.000.83
Zn (Zinc) mg/L9130.14126.650.0214.250.017.460.018.01
TETI 841 5412 5360 3591
Ci—mean concentration values.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mansilha, C.; Melo, A.; Flores, D.; Ribeiro, J.; Rocha, J.R.; Martins, V.; Santos, P.; Espinha Marques, J. Irrigation with Coal Mining Effluents: Sustainability and Water Quality Considerations (São Pedro da Cova, North Portugal). Water 2021, 13, 2157. https://doi.org/10.3390/w13162157

AMA Style

Mansilha C, Melo A, Flores D, Ribeiro J, Rocha JR, Martins V, Santos P, Espinha Marques J. Irrigation with Coal Mining Effluents: Sustainability and Water Quality Considerations (São Pedro da Cova, North Portugal). Water. 2021; 13(16):2157. https://doi.org/10.3390/w13162157

Chicago/Turabian Style

Mansilha, Catarina, Armindo Melo, Deolinda Flores, Joana Ribeiro, João Ramalheira Rocha, Vítor Martins, Patrícia Santos, and Jorge Espinha Marques. 2021. "Irrigation with Coal Mining Effluents: Sustainability and Water Quality Considerations (São Pedro da Cova, North Portugal)" Water 13, no. 16: 2157. https://doi.org/10.3390/w13162157

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

Article Metrics

Back to TopTop