Spatial distribution pattern and health risk of groundwater contamination by cadmium, manganese, lead and nitrate in groundwater of an arid area.

Combining the results of base models to create a meta-model is one of the ensemble approaches known as stacking. In this study, stacking of five base learners, including eXtreme gradient boosting, random forest, feed-forward neural networks, generalized linear models with Lasso or Elastic Net regularization, and support vector machines, was used to study the spatial variation of Mn, Cd, Pb, and nitrate in Qom-Kahak Aquifers, Iran. The stacking strategy proved to be an effective substitute predictor for existing machine learning approaches due to its high accuracy and stability when compared to individual learners. Contrarily, there was not any best-performing base model for all of the involved parameters. For instance, in the case of cadmium, random forest produced the best results, with adjusted R2 and RMSE of 0.108 and 0.014, as opposed to 0.337 and 0.013 obtained by the stacking method. The Mn and Cd showed a tight link with phosphate by the redundancy analysis (RDA). This demonstrates the effect of phosphate fertilizers on agricultural operations. In order to analyze the causes of groundwater pollution, spatial methodologies can be used with multivariate analytic techniques, such as RDA, to help uncover hidden sources of contamination that would otherwise go undetected. Lead has a larger health risk than nitrate, according to the probabilistic health risk assessment, which found that 34.4% and 6.3% of the simulated values for children and adults, respectively, were higher than HQ = 1. Furthermore, cadmium exposure risk affected 84% of children and 47% of adults in the research area.

Development and applications of GIS-based spatial analysis in environmental geochemistry in the big data era.

The research of environmental geochemistry entered the big data era. Environmental big data is a kind of new method and thought, which brings both opportunities and challenges to GIS-based spatial analysis in geochemical studies. However, big data research in environmental geochemistry is still in its preliminary stage, and what practical problems can be solved still remain unclear. This short review paper briefly discusses the main problems and solutions of spatial analysis related to the big data in environmental geochemistry, with a focus on the development and applications of conventional GIS-based approaches as well as advanced spatial machine learning techniques. The topics discussed include probability distribution and data transformation, spatial structures and patterns, correlation and spatial relationships, data visualisation, spatial prediction, background and threshold, hot spots and spatial outliers as well as distinction of natural and anthropogenic factors. It is highlighted that the integration of spatial analysis on the GIS platform provides effective solutions to revealing the hidden spatial patterns and spatially varying relationships in environmental geochemistry, demonstrated by an example of cadmium concentrations in the topsoil of Northern Ireland through hot spot analysis. In the big data era, further studies should be more inclined to the integration and application of spatial machine learning techniques, as well as investigation on the temporal trends of environmental geochemical features.

Identification of high ecological risk areas with naturally high background value of soil Cd related to carbonate rocks.

The characteristics of high concentrations or high activity levels of heavy metals, especially Cd, in soils caused by the pedogenesis of rocks are attracting increased attention. Carbonate rocks and black shales often coexist during geological deposition, but the risk characteristics of heavy metals are different after their weathering into the soil. The purpose of this study was to investigate the element concentrations of a naturally high background value area, to identify patterns of different risk areas, and to make recommendations for the safe usage of farmland. The results showed that, compared with the soil in the carbonate rock area, the soil in the black shale area was more acidified and most of the heavy metal elements were leached. Based on the soil pH value and the heavy metal concentrations, an identification method for land risk areas within naturally high background values was established, and land planning was carried out using this method. The exceeding rates of Cd in rice for the preferential protection area and strict control area were 0.0 and 50.0%, respectively. Therefore, in naturally high background area, the identified lithology can apply to maximize the use of farmland resources. This method provides a basis for preliminary ecological risk screening in naturally high background value areas using the results of the soil survey. A suggestion for the prevention and control of soil pollution in areas with naturally high background values was put forward. In carbonate rock areas, the soil should be closely monitored to prevent soil acidification.

Identification of the spatial patterns and controlling factors of Se in soil and rice in Guangxi through hot spot analysis.

Selenium (Se) is essential to human health, anti-cancer, possessing antioxidant, and antiviral properties. In this study, the spatial patterns of rice Se and their varying relationship with soil Se on a regional scale were studied using hot spot analysis for the agricultural soils in Guangxi. According to the hot and cold spot maps, rice Se correlates positively with soil Se in Guangxi agricultural soils. High rice Se accompanies high soil Se in the central part of Guangxi (e.g., Liuzhou, Laibin), and low rice Se is in line with low soil Se in the western part (e.g., Baise). However, the hot spot analysis maps indicate that southwestern Guangxi exhibits a special characteristic of low rice Se with high soil Se (e.g., Chongzuo). This special pattern is strongly associated with the high concentrations of Fe2O3 (ferromanganese nodules) in the carbonate rock area. The hot spot analysis proves useful in revealing the spatial patterns of rice Se in Guangxi and identifying the hidden patterns.

Highly effective stabilization of Cd and Cu in two different soils and improvement of soil properties by multiple-modified biochar.

Heavy metal contamination in soil has attracted great attention worldwide. In situ stabilization has been considered an effective way to remediate soils contaminated by heavy metals. In the present research, a multiple-modified biochar (BCM) was prepared to stabilize Cd and Cu contamination in two different soils: a farmland soil (JYS) and a vegetable soil (ZZS). The results showed that BCM was a porous-like flake material and that modification increased its specific surface area and surface functional groups. The incubation experiment indicated that BCM decreased diethylenetriaminepentaacetic (DTPA)-extractable Cd and Cu by 92.02% and 100.00% for JYS and 90.27% and 100.00% for ZZS, respectively. The toxicity characteristic leaching procedure (TCLP)-extractable Cd and Cu decreased 66.46% and 100.00% for JYS and 46.33% and 100.00% for ZZS, respectively. BCM also reduced the mobility of Cd and Cu in soil and transformed them to more stable fractions. In addition, the application of BCM significantly increased the soil dehydrogenase, organic matter content and available K (p < 0.05). These results indicate that BCM has great potential in the remediation of Cd- and Cu-contaminated soil.

Long-term stabilization of Cd in agricultural soil using mercapto-functionalized nano-silica (MPTS/nano-silica): A three-year field study.

Cadmium (Cd) contamination in agricultural soil is a worldwide environmental problem. In situ stabilization has been considered an effective approach for the remediation of Cd-contaminated agricultural soil. However, information about the long-term effects of amendment on soil properties and stabilization efficiency remains limited. In the present study, mercapto-functionalized nano-silica (MPTS/nano-silica) was used to stabilize Cd in contaminated agricultural soil under field conditions for three years (with application rates of 0%, 0.2%, 0.4%, 0.6%, 0.8% and 1.0%). The application of MPTS/nano-silica reduced the soil aggregate stability (PDA0.25) (14.8%) and available K (24.9%) and significantly increased the soil dehydrogenase (DHA) (43.4%), yield of wheat grains (33.5%) and Si content in wheat tissues (55.2% in leaf, 50.4% in stem, and 37.7% in husk) (p < 0.05). More importantly, MPTS/nano-silica decreased the leachability (36.0%) and bioavailability (54.3%) of Cd in the soil and transformed Cd into a more stable fraction. The content of Cd in wheat grains decreased by 53.9%, 61.9% and 54.1% in 2017, 2018 and 2019, respectively, in comparison with the control. These results indicated that MPTS/nano-silica has long-term stabilization effects on Cd in agricultural soil and is a potential amendment for the remediation of Cd-contaminated agricultural soils.

Factors affecting simultaneous nitrification and denitrification (SND) in a moving bed sequencing batch reactor (MBSBR) system as revealed by microbial community structures.

The effects of biological factors including dissolved oxygen (DO), pH, carbon/nitrogen (C/N) and hydraulic retention times (HRT) on the performance of simultaneous nitrification and denitrification (SND) in a moving bed sequencing batch reactor (MBSBR) were investigated. A low DO was found to be advantageous to the SND in that nitrification was not inhibited, while pH and C/N ratio were shown to have positive effects on SND, and HRT needed to be controlled in a suitable range. A desirable SND efficiency was obtained at a DO of 2.5 mg L-1, pH of approximately 8.0, C/N ratio of 10 and HRT of 10 h in the MBSBR. High-throughput sequencing analysis showed that different operating conditions impacted microbial communities, resulting in different nitrogen removal mechanisms. Autotrophic and heterotrophic nitrification together contributed to the good nitrification performance, while denitrification was conducted by combined anoxic and aerobic processes. Furthermore, the results of principal component analyses (PCA) and the abundance of the predominant nitrification and denitrification genera both showed that DO and HRT might be regarded as the dominant variable factors influencing community structure analysis during SND, while the linear discriminant analysis (LDA) effect size (LEfSe) algorithm showed differences in abundance among the biofilm microbial communities with different DO. Overall, the results of this study improve our understanding of the bacterial community structure with different operating conditions in MBSBRs.

Efficiency and kinetics of conventional pollutants and tetracyclines removal in integrated vertical-flow constructed wetlands enhanced by aeration.

The integrated vertical-flow constructed wetland (IVCW) is considered as a potential alternative for domestic wastewater treatment of towns and small cities. Oxygen supply is the main limitation of pollutants removal in IVCWs. In the present study, a field experiment was conducted to evaluate the capacity and kinetics of pollutants removal in IVCWs with/without artificial aeration. Two IVCWs constructed with Canna indica and Phragmites australis were running in continuous flow to remove high concentrations of conventional pollutants and low concentrations of tetracyclines (TETs), which are at similar levels of domestic wastewater. The results showed that IVCWs had a good performance on COD, phosphorus, and TETs with removal efficiencies over 80%, 64%, and 75%, respectively, with a hydraulic retention time (HRT) of 3.0 d. However, the removal of nitrogen was limited, showing as TN removal efficiency of about 30%. The IVCW with Phragmites australis had a higher removal efficiency and rate. A kinetics based on Monod Equation and solved with Matlab 2018a could describe the degradation of conventional pollutants. Artificial aeration improved the oxygen supply and remarkably raised the removal capacity for COD, N, and P in IVCWs. The q1/2 values, which was defined as the average removal loading before half of the pollutants was removed and represented the removal capacity without limitation of pollutants concentration, were increased by 5-30 times after aeration. In conclusion, IVCWs could remove conventional pollutants and TETs simultaneously showing a great potential in domestic wastewater treatment. Artificial aeration enhanced removal capacity of IVCWs on conventional pollutants while showed little influence on TETs.

Peak Chromium Pollution in Summer and Winter Caused by High Mobility of Chromium in Sediment of a Eutrophic Lake: In Situ Evidence from High Spatiotemporal Sampling.

To study the mechanisms of chromium (Cr) mobilization in sediments of eutrophic lakes, monthly sampling was performed in the Meiliang Bay of Lake Taihu, China, combined with laboratory experiments. High-resolution dialysis and diffusive gradients in thin film (DGT) sampling techniques were used. Results indicated that in July 2016 and January 2017, the concentrations of soluble Cr and DGT-labile Cr(VI) in the overlying water exceeded both drinking and fishery water quality standards, resulting from the high mobility of Cr in sediments. In July (summer), the high concentration of soluble Cr (134.04 ± 7.20 µg/L) detected in the anaerobic sediments was primarily caused by the complexation of Cr(III) with dissolved organic matter (DOM). This mechanism was supported by an observed simultaneous increase of soluble Cr and DOM under simulated anaerobic conditions. In January (winter), the high concentrations of soluble Cr (97.55 ± 9.65 µg/L) and DGT-labile Cr(VI) (25.83 ± 1.25 µg/L) in aerobic sediments were primarily caused by reoxidation of Cr(III) by Mn(III/IV) oxides as evidenced by the lowest concentrations of soluble and DGT-labile Mn(II). This study sheds light on the full-year variation and mechanisms of Cr mobilization in eutrophic lake sediments and suggests the urgent need for remediation of Cr pollution especially for winter.

Research progress and application prospect of anaerobic biological phosphorus removal.

Anaerobic biological phosphorus removal has proposed a new direction for the removal of phosphorus from wastewater, and the discovery of phosphate reduction makes people have a more comprehensive understanding of microbial phosphorus cycling. Here, from the perspective of thermodynamics, the bioreduction reaction of phosphate was analyzed and its mechanism was discussed. The research progress of phosphate reduction and the application prospects of anaerobic biological phosphorus removal from wastewater were introduced, pointing out the situation and guiding the further research in this field.

Intensification of phosphorus cycling in China since the 1600s.

Phosphorus (P) is an essential nutrient for living systems with emerging sustainability challenges related to supply uncertainty and aquatic eutrophication. However, its long-term temporal dynamics and subsequent effects on freshwater ecosystems are still unclear. Here, we quantify the P pathways across China over the past four centuries with a life cycle process-balanced model and evaluate the concomitant potential for eutrophication with a spatial resolution of 5 arc-minutes in 2012. We find that P cycling in China has been artificially intensified during this period to sustain the increasing population and its demand for animal protein-based diets, with continuous accumulations in inland waters and lands. In the past decade, China's international trade of P involves net exports of P chemicals and net imports of downstream crops, specifically soybeans from the United States, Brazil, and Argentina. The contribution of crop products to per capita food P demand, namely, the P directly consumed by humans, declined from over 98% before the 1950s to 76% in 2012, even though there was little change in per capita food P demand. Anthropogenic P losses to freshwater and their eutrophication potential clustered in wealthy coastal regions with dense populations. We estimate that Chinese P reserve depletion could be postponed for over 20 y by more efficient life cycle P management. Our results highlight the importance of closing the P cycle to achieve the cobenefits of P resource conservation and eutrophication mitigation in the world's most rapidly developing economy.

Medical geological study of disease-causing elements in Wassa area of Southwest Ghana.

Insufficient information on the link between health data and geology in developing countries is a major barrier to identify sources of some emerging public health problems. A total of 2868 soil samples were collected from field sheet 0503B in Ghana to evaluate the concentrations and distributions of trace elements and their effects on human health. The samples were sieved to < 106 µm fraction and analysed for elements, As, Ba, K, Zn, Co, Cr, Cu, Mn, Ni, Pb, Mg and Fe by XRF technique and Au by fire assay method. The study identified disparities in averages of As, Cr, Fe and Mg, which resulted in enrichment and deficiencies when compared with the worldwide background average. The measured averages for As and Cr were 17.27 mg/kg and 89.25 mg/kg, respectively, for the entire area. Both averages exceeded the worldwide background values of 10 mg/kg and 8 mg/kg of As and Cr. The four traditional towns with varied activities recorded As concentrations ranging from 6.11 mg/kg at Samreboi, 16.29 mg/kg at Asankragwa, 17.42 mg/kg at Akropong and 25.99 mg/kg at Bogoso. Principal component analysis revealed a good association among Ba, Cr, Cu, Fe, K, Ni, Pb and Zn in Group 1, and their main source was interpreted as the underlying geology. Arsenic, Cr and Mg in Group 2 show a relatively weak correlation, and their sources were ascribed to a combination of geologic and anthropogenic sources. Gold had a good correlation with As, which was associated with the hydrothermal veins in the underlying rocks. The spatial plots generated from transformed soil data by Getis Ord Gi* treatments were visual methods to clearly identify geographically the hotspots and coldspots of elements that cause diseases.

Simultaneous measurements of arsenic and sulfide using diffusive gradients in thin films technique (DGT).

Diffusive gradients in thin films technique (DGT) is a dynamically passive sampling technique which has been applied increasingly to the environmental monitoring field. In the preliminary period, the DGT with zirconium hydroxide-silver iodide as the binding phase (ZrO-AgI DGT) has been developed for the determination of sulfide (S(II)). On this basis, this paper developed its determination method for inorganic arsenite (As(III)) to further realize the simultaneous and high-resolution measurements of labile inorganic As and S(II) in sediments. ZrO-AgI binding gel had a strong ability in adsorbing and fixing As(III), showing a linear increase in the initial 12.5 min. After saturation of S(II) on ZrO-AgI binding gel, the adsorption rate and adsorption capacity of As(III) reduced by 8 and 14%, respectively. A stable elution rate (89.1 ± 2.2%) was obtained by extraction of As(III) on the binding gel using a mixture solution of 1.0 M NaOH and 1.0 M H2O2 (1:1). The DGT capacity of As(III) determined by the ZrO-AgI DGT was 23.6 µg cm-2. DGT uptakes of As(III) were independent of pH (4.0-9.0) and ionic strength (0.01-100 mM), and they did not interfere with each other during the uptake process. Simultaneous measurements of labile As and S(II) in four sediment cores of Taihu Lake (China) with ZrO-AgI DGT showed that they had similarly vertical distributions in the top 16-mm layer in one core and in the whole profile up to the 35 mm depth in two cores. It likely reflected a simultaneous release of As and S(II) in sediments by synchronous reduction of As-hosted oxidized iron and sulfate, respectively. The simultaneous decreases of labile As and S(II) from their co-precipitation (e.g., As2S3) were not obvious in deeper sediment layer through the measurement with ZrO-AgI DGT.

Impact of grass cover on the magnetic susceptibility measurements for assessing metal contamination in urban topsoil.

In recent decades, magnetic susceptibility monitoring has developed as a useful technique in environmental pollution studies, particularly metal contamination of soil. This study provides the first ever examination of the effects of grass cover on magnetic susceptibility (MS) measurements of underlying urban soils. Magnetic measurements were taken in situ to determine the effects on κ (volume magnetic susceptibility) when the grass layer was present (κgrass) and after the grass layer was trimmed down to the root (κno grass). Height of grass was recorded in situ at each grid point. Soil samples (n=185) were collected and measurements of mass specific magnetic susceptibility (χ) were performed in the laboratory and frequency dependence (χfd%) calculated. Metal concentrations (Pb, Cu, Zn and Fe) in the soil samples were determined and a gradiometry survey carried out in situ on a section of the study area. Significant correlations were found between each of the MS measurements and the metal content of the soil at the p<0.01 level. Spatial distribution maps were created using Inverse Distance Weighting (IDW) and Local Indicators of Spatial Association (LISA) to identify common patterns. κgrass (ranged from 1.67 to 301.00×10-5 SI) and κno grass (ranged from 2.08 to 530.67×10-5 SI) measured in situ are highly correlated [r=0.966, n=194, p<0.01]. The volume susceptibility datasets in the presence and absence of grass coverage share a similar spatial distribution pattern. This study re-evaluates in situ κ monitoring techniques and the results suggest that the removal of grass coverage prior to obtaining in situ κ measurements of urban soil is unnecessary. This layer does not impede the MS sensor from accurately measuring elevated κ in soils, and therefore κ measurements recorded with grass coverage present can be reliably used to identify areas of urban soil metal contamination.

Simultaneous Measurements of Eight Oxyanions Using High-Capacity Diffusive Gradients in Thin Films (Zr-Oxide DGT) with a High-Efficiency Elution Procedure.

A zirconium oxide binding gel-based diffusive gradients in thin films (Zr-oxide DGT) was developed for simultaneous measurements of P(V), As(V), Cr(VI), Mo(VI), Sb(V), Se(VI), V(V), and W(VI). All of the oxyanions were rapidly bound to Zr-oxide gel with differences in binding affinity. The eight bound oxyanions were successfully recovered by one-step elution using a mild reagent of 0.2 M NaOH-0.5 M H2O2 by overcoming the problems in analyses of the oxyanions. The optimized elution time was reduced to 3-5 h from 24-48 h required by other DGTs. DGT uptakes of all the oxyanions were independent of pH (4.42-8.45) and ionic strength (0.1-500 mM). The DGT capacities for six oxyanions detected in multioxyanion solution were only 0.19 to 0.35 times of those detected in single-oxyanion solution, reflecting a strong competition among the oxyanions during DGT uptake. Except for Se(VI) in seawater, Zr-oxide DGT accurately measured all of the oxyanions in synthetic freshwater and seawater, with the capacities ∼29 to >2397 times and ∼7.5 to 232 times those of two commonly used DGTs (Metsorb and precipitated ferrihydrite (PF) DGTs) in freshwater and seawater, respectively. Measurements by Zr-oxide DGT in contaminated sediments were in agreement with only two oxyanions with the two commonly used DGTs; the two DGTs accumulated less or no mass of other oxyanions. This study demonstrates significant advantage of Zr-oxide DGT over the other DGTs in simultaneous measurements of the eight oxyanions due to the former's high capacity and a wide tolerance to environmental interferences, together with a high efficiency in elution.

New trends in removing heavy metals from wastewater.

With the development of researches, the treatments of wastewater have reached a certain level. Whereas, heavy metals in wastewater cause special concern in recent times due to their recalcitrance and persistence in the environment. Therefore, it is important to get rid of the heavy metals in wastewater. The previous studies have provided many alternative processes in removing heavy metals from wastewater. This paper reviews the recent developments and various methods for the removal of heavy metals from wastewater. It also evaluates the advantages and limitations in application of these techniques. A particular focus is given to innovative removal processes including adsorption on abiological adsorbents, biosorption, and photocatalysis. Because these processes have leaded the new trends and attracted more and more researches in removing heavy metals from wastewater due to their high efficency, pluripotency and availability in a copious amount. In general, the applicability, characteristic of wastewater, cost-effectiveness, and plant simplicity are the key factors in selecting the most suitable method for the contaminated wastewater.

A combined biological removal of Cd(2+) from aqueous solutions using Phanerochaete chrysosporium and rice straw.

The removal of Cd(2+) from aqueous solutions by agricultural residues rice straw combined with white rot fungus Phanerochaete chrysosporium (P. chrysosporium) was investigated. The results showed that over 99% of the total Cd(2+) (initial concentration of 150mgL(-1)) was removed at the optimal operating conditions (pH 5.0 at 35°C). We also found that P. chrysosporium could survive under Cd(2+) stress even with an initial Cd(2+) concentration of 250mgL(-1). But when Cd(2+) concentration increased to 250mgL(-1), fungus growth and reproduction were remarkably restrained, and as a result, Cd(2+) removal dropped to 59.2%. It was observed that the fungus biomass and activities of ligninolytic enzymes decreased at some degree under high concentration of Cd(2+) (above 100mgL(-1)). Also, we found that a moderate Cd(2+) stress (below 150mgL(-1)) could stimulate P. chrysosporium's production of the heavy metals chelator - oxalate. This study will provide useful information for the application of biological removal of heavy metal irons from wastewater.

Effects of constant pH and unsteady pH at different free ammonia concentrations on shortcut nitrification for landfill leachate treatment.

On the basis of achieving shortcut nitrification in a lab-scale SBR, the effects of constant pH and unsteady pH at different free ammonia concentrations on shortcut nitrification for landfill leachate treatment was investigated. The results indicate that under the condition of DO of 0.5 ± 0.2 mg/L and temperature of 30 ± 2 °C, the absolute value of nitrite accumulation increased significantly with the increase in free ammonia (FA) concentration from 5.30 to 48.67 mg/L; however, the nitrite accumulation rate remained almost constant at a constant pH of 8.0 ± 0.1. Ammonia oxidation and the nitrite accumulation become slow with the pH decreased from 8.0 ± 0.1 to 7.5 ± 0.2, and the activities of ammonia-oxidizing bacteria (AOB) and nitrite-oxidizing bacteria (NOB) were severely inhibited when the pH further decreased to 6.5. More importantly, this study confirmed that the pH decrease from 8.0 to 6.5 within a short time exhibited significant negative effect on the ammonia oxidation rather than the FA concentration.

Pollution profiles and risk assessment of PBDEs and phenolic brominated flame retardants in water environments within a typical electronic waste dismantling region.

The aim of this study was to assess the pollution profiles of various typical brominated flame retardants in water and surface sediment near a typical electronic waste dismantling region in southern China. We found that polybrominated diphenyl ethers (PBDEs), 2,4,6-tribromophenol (TBP), pentabromophenol (PeBP), tetrabromobisphenol A (TBBPA), and bisphenol A (BPA) were ubiquitous in the water and sediment samples collected in the study region. In water, Σ19PBDEs (sum of all 20 PBDE congeners studied except BDE-209, which was below the detection limit) levels ranged from 0.31 to 8.9 × 10(2) ng L(-1). TBP, PeBP, TBBPA, and BPA concentrations in the water samples ranged from not being detectable (nd-under the detection limit) to 3.2 × 10(2) (TBP), from nd to 37 (PeBP), from nd to 9.2 × 10(2) (TBBPA) and from nd-8.6 × 10(2) ng L(-1) (BPA). In sediment, Σ19PBDEs ranged from nd to 5.6 × 10(3) ng g(-1), while BDE-209 was the predominant congener, with a range of nd to 3.5 × 10(3) ng g(-1). Tri- to hepta-BDE concentrations were significantly (p < 0.01) correlated with each other, except for BDE-71 and BDE-183, and octa- to nona-BDEs concentrations were significantly (p < 0.05) correlated with each other, except for BDE-208. BDE-209 was not significantly correlated with tri- to nona-BDEs. Risk assessments indicated that the water and sediment across the sampling sites posed no estrogenic risk. However, different eco-toxicity risk degrees at three trophic levels did exist at most sampling sites.