Environmental Filtering Process Has More Important Roles than Dispersal Limitation in Shaping Large-Scale Prokaryotic Beta Diversity Patterns of Grassland Soils.

Despite the utmost importance of microorganisms in maintaining ecosystem functioning and their ubiquitous distribution, our knowledge of the large-scale pattern of microbial diversity is limited, particularly in grassland soils. In this study, the microbial communities of 99 soil samples spanning over 3000 km across grassland ecosystems in northern China were investigated using high-throughput sequencing to analyze the beta diversity pattern and the underlying ecological processes. The microbial communities were dominated by Proteobacteria, Actinobacteria, Acidobacteria, Chloroflexi, and Planctomycetes across all the soil samples. Spearman's correlation analysis indicated that climatic factors and soil pH were significantly correlated with the dominant microbial taxa, while soil microbial richness was positively linked to annual precipitation. The environmental divergence-dissimilarity relationship was significantly positive, suggesting the importance of environmental filtering processes in shaping soil microbial communities. Structural equation modeling found that the deterministic process played a more important role than the stochastic process on the pattern of soil microbial beta diversity, which supported the predictions of niche theory. Partial mantel test analysis have showed that the contribution of independent environmental variables has a significant effect on beta diversity, while independent spatial distance has no such relationship, confirming that the deterministic process was dominant in structuring soil microbial communities. Overall, environmental filtering process has more important roles than dispersal limitation in shaping microbial beta diversity patterns in the grassland soils.

Paenibacillus tibetensis sp. nov., a psychrophilic bacterium isolated from alpine swamp meadow soil.

A novel psychrophilic strain, SSB001(T), was isolated from an alpine swamp meadow soil in Tibet, China, and identified as a representative of a novel phylogenetic subclade in the genus Paenibacillus , with Paenibacillus antarcticus (96.2%), Paenibacillus macquariensis (96.53%) and Paenibacillus glacialis (96.2%) as the most closely related species on the basis of 16S rRNA gene sequence analyses. The strain was distinguished from defined species of the genus Paenibacillus by further study of rpoB gene sequences, phenotypic characterization, cellular fatty acid composition, quinones, polar lipids and meso-diaminopimelic acid in the peptidoglycan. Based upon these results, we propose the strain as a representative of a novel species named Paenibacillus tibetensis sp. nov., with SSB001(T) ( =ACCC 19728(T) =DSM 29321(T)) as the type strain. The DNA G+C content (mol%) of strain SSB001(T) was 40.18 mol% (HPLC).

Altitudinal distribution patterns of soil bacterial and archaeal communities along mt. Shegyla on the Tibetan Plateau.

Unraveling the distribution patterns of plants and animals along the elevational gradients has been attracting growing scientific interests of ecologists, whether the microbial communities exhibit similar elevational patterns, however, remains largely less documented. Here, we investigate the biogeographic distribution of soil archaeal and bacterial communities across three vertical climate zones (3,106-4,479 m.a.s.l.) in Mt. Shegyla on the Tibetan Plateau, by combining quantitative PCR and high-throughput barcoded pyrosequencing approaches. Our results found that the ratio of bacterial to archaeal 16S rRNA gene abundance was negatively related with elevation. Acidobacteria dominated in the bacterial communities, Marine benthic group A dominated in the archaeal communities, and the relative abundance of both taxa changed significantly with elevation. At the taxonomic levels of domain, phylum, and class, more bacterial taxa than archaeal exhibited declining trend in diversity along the increasing elevational gradient, as revealed by Shannon and Faith's phylogenetic diversity indices. Unweighted UniFrac distance clustering showed that the bacterial communities from the mountainous temperate zone clustered together, whereas those from the subalpine cool temperate zone clustered together. However, the partitioning effect of elevational zones on the archaeal community was much weaker compared to that on bacteria. Redundancy analysis revealed that soil geochemical factors explained 58.3 % of the bacterial community variance and 75.4 % of the archaeal community variance. Taken together, we provide evidence that soil bacteria exhibited more apparent elevational zonation feature and decreased diversity pattern than archaea with increasing elevation, and distribution patterns of soil microbes are strongly regulated by soil properties along elevational gradient in this plateau montane ecosystem.

Responses of soil ammonia oxidizers to a short-term severe mercury stress.

The responses of soil ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) to mercury (Hg) stress were investigated through a short-term incubation experiment. Treated with four different concentrations of Hg (CK, Hg25, Hg50, and Hg100, denoting 0, 25, 50, and 100mgHg/kg dry soil, respectively), samples were harvested after 3, 7, and 28day incubation. Results showed that the soil potential nitrification rate (PNR) was significantly inhibited by Hg stress during the incubation. However, lower abundances of AOA (the highest in CK: 9.20×10(7)copies/g dry soil; the lowest in Hg50: 2.68×10(7)copies/g dry soil) and AOB (the highest in CK: 2.68×10(7)copies/g dry soil; the lowest in Hg50: 7.49×10(6)copies/g dry soil) were observed only at day 28 of incubation (P<0.05). Moreover, only the community structure of soil AOB obviously shifted under Hg stress as seen through DGGE profiles, which revealed that 2-3 distinct AOB bands emerged in the Hg treatments at day 28. In summary, soil PNR might be a very useful parameter to assess acute Hg stress on soil ecosystems, and the community structure of soil AOB might be a realistic biological indicator for the assessment of heavy metal stress on soil ecosystems in the future.

Analysis of the microbial community structure by monitoring an Hg methylation gene (hgcA) in paddy soils along an Hg gradient.

Knowledge of the diversity of mercury (Hg)-methylating microbes in the environment is limited due to a lack of available molecular biomarkers. Here, we developed novel degenerate PCR primers for a key Hg-methylating gene (hgcA) and amplified successfully the targeted genes from 48 paddy soil samples along an Hg concentration gradient in the Wanshan Hg mining area of China. A significant positive correlation was observed between hgcA gene abundance and methylmercury (MeHg) concentrations, suggesting that microbes containing the genes contribute to Hg methylation in the sampled soils. Canonical correspondence analysis (CCA) showed that the hgcA gene diversity in microbial community structures from paddy soils was high and was influenced by the contents of total Hg, SO4(2-), NH4(+), and organic matter. Phylogenetic analysis showed that hgcA microbes in the sampled soils likely were related to Deltaproteobacteria, Firmicutes, Chloroflexi, Euryarchaeota, and two unclassified groups. This is a novel report of hgcA diversity in paddy habitats, and results here suggest a link between Hg-methylating microbes and MeHg contamination in situ, which would be useful for monitoring and mediating MeHg synthesis in soils.

N-oxide reduction of quinoxaline-1,4-dioxides catalyzed by porcine aldehyde oxidase SsAOX1.

Quinoxaline-1,4-dioxides (QdNOs) are a class of quinoxaline derivatives that are widely used in humans or animals as drugs or feed additives. However, the metabolic mechanism, especially the involved enzymes, has not been reported in detail. In this study, the N-oxide reduction enzyme, porcine aldehyde oxidase SsAOX1 was identified and characterized. The SsAOX1 gene was cloned from pig liver through reverse-transcription polymerase chain reaction using degenerate primers, which encode a 147-kDa protein with typical aldehyde oxidase motifs, two [2Fe-2S] centers, a flavin adenine dinucleotide (FAD) binding domain, and a molybdenum cofactor domain. After heterologous expression in a prokaryote, purified SsAOX1 formed a functional homodimer under native conditions. Importantly, the SsAOX1 catalyzed the N-oxide reduction at the N1 position of three representative QdNOs (quinocetone, mequindox, and cyadox), which are commonly used as animal feed additives. SsAOX1 has the highest activity toward quinocetone, followed by mequindox and cyadox, with kcat/K(m) values of 1.94 ± 0.04, 1.27 ± 0.15, and 0.43 ± 0.09 minute(-1) µM(-1), respectively. However, SsAOX1 has the lowest substrate affinity for quinocetone, followed by the cyadox and mequindox, with K(m) values of 4.36 ± 0.56, 3.16 ± 0.48, and 2.96 ± 0.51 µM, respectively. In addition, using site-directed mutagenesis, we found that substitution of glycine 1019 with threonine endows SsAOX1 with N-oxide reductive activity at the N4 position. The goal of this study was to identify and characterize the N-oxide reduction enzyme for a class of veterinary drugs, QdNOs, which will aid in the elucidation of the metabolic pathways of QdNOs and will provide a theoretical basis for their administration and new veterinary drug design.

Initial copper stress strengthens the resistance of soil microorganisms to a subsequent copper stress.

To improve the prediction of essential ecosystem functioning under future environmental disturbances, it is of significance to identify responses of soil microorganisms to environmental stresses. In this study, we collected polluted soil samples from field plots with eight copper levels ranging from 0 to 3,200 mg Cu kg(-1) soil. Then, the soils with 0 and 3,200 mg Cu kg(-1) were selected to construct a microcosm experiment. Four treatments were set up including Cu0-C and Cu3200-C without further Cu addition, and Cu0-A and Cu3200-A with addition of 57.5 mg Cu kg(-1) soil. We measured substrate-induced respiration (SIR) and potential nitrification rate (PNR). Furthermore, the abundance of bacterial, archaeal 16S rRNA genes, ammonia-oxidizing bacteria and archaea amoA genes were determined through quantitative PCR. The soil microbial communities were investigated by terminal restriction fragment length polymorphism (T-RFLP). For the field samples, the SIR and PNR as well as the abundance of soil microorganisms varied significantly between eight copper levels. Soil microbial communities highly differed between the low and high copper stress. In the microcosm experiment, the PNR and SIR both recovered while the abundance of soil microorganisms varied irregularly during the 90-day incubation. The differences of microbial communities measured by pairwise Bray-Curtis dissimilarities between Cu0-A and Cu0-C on day 0 were significantly higher after subsequent stress than before. However, the differences of microbial communities between Cu3200-A and Cu3200-C on day 0 changed little between after subsequent stress and before. Therefore, initial copper stress could increase the resistance of soil microorganisms to subsequent copper stress.

Patterns of bacterial diversity along a long-term mercury-contaminated gradient in the paddy soils.

Mercury (Hg) pollution is usually regarded as an environmental stress in reducing microbial diversity and altering bacterial community structure. However, these results were based on relatively short-term studies, which might obscure the real response of microbial species to Hg contamination. Here, we analysed the bacterial abundance and community composition in paddy soils that have been potentially contaminated by Hg for more than 600 years. Expectedly, the soil Hg pollution significantly influenced the bacterial community structure. However, the bacterial abundance was significantly correlated with the soil organic matter content rather than the total Hg (THg) concentration. The bacterial alpha diversity increased at relatively low levels of THg and methylmercury (MeHg) and subsequently approached a plateau above 4.86 mg kg(-1) THg or 18.62 ng g(-1) MeHg, respectively. Contrasting with the general prediction of decreasing diversity along Hg stress, our results seem to be consistent with the intermediate disturbance hypotheses with the peak biological diversity under intermediate disturbance or stress. This result was partly supported by the inconsistent response of bacterial species to Hg stress. For instance, the relative abundance of Nitrospirae decreased, while that of Gemmatimonadetes increased significantly along the increasing soil THg and MeHg concentrations. In addition, the content of SO(4)(2-), THg, MeHg and soil depth were the four main factors influencing bacterial community structures based on the canonical correspondence analysis (CCA). Overall, our findings provide novel insight into the distribution patterns of bacterial community along the long-term Hg-contaminated gradient in paddy soils.

Toxicity of profenofos to the springtail, Folsomia candida, and ammonia-oxidizers in two agricultural soils.

Extensive use of organophosphorus insecticide profenofos (PFF) for agricultural and house-hold purposes has led to serious environmental pollution, with potential risk to organisms in the ecosystem. This study examined the toxicity of PFF to the soil springtail Folsomia candida and ammonia-oxidizers through a series of toxicity tests conducted on two agricultural soils. It was found that the survival, reproduction, hsp70 gene expression of F. candida and the soil potential nitrification rate (PNR) were sensitive to the PFF, whereas no apparent change was observed in the abundance of ammonia-oxidizers. The reproduction of F. candida was the most sensitive endpoint (mean 0.10 mg/kg of EC(50) value) for PFF, although the test was more time-consuming. The results of the acute toxicity tests suggested that the survival of F. candida could be considered as the most suitable bioindicator for fast screening of PFF toxicity because of its fast and easy test procedure. In addition, the hsp70 gene expression in F. candida and the PNR could be used as important parameters for assessment of PFF toxicity. The threshold concentration based on the obtained endpoints differed in the two soils, and consequently the soil property should be considered in toxicity assessments of contaminated soils.

Effects of long-term fertilization on the diversity of bacterial mercuric reductase gene in a Chinese upland soil.

Soil mercury (Hg) pollution has received considerable attention due to its neurotoxin effects and its potential risk to food safety. The microbial transformation of Hg plays a key role in reducing Hg toxicity by the mercuric reductase (MerA) conferred by genes arranged in the mer operon. This study investigated the effects of long-term fertilization on the diversity of bacterial mercuric reductase gene (merA), which specify the reduction of ionic Hg²âº to the volatile elemental form Hg°, in an agricultural soil with relatively high Hg content. The soil samples were collected from different treatments, including control without fertilizer (CK), fertilizer nitrogen (N), combined fertilizers (NPK) of N, phosphorus (P) and potassium (K), and NPK plus organic manure (NPK + OM). The merA gene diversity patterns were analyzed based on the merA clone libraries and sequencing measurements. Results showed that the merA gene diversity was influenced by soil variables depending on the fertilization practices. In particular, NH4⁺ and NO3⁻ contents had strong effect on the merA gene diversity pattern both in the N and NPK treatments, whereas the merA gene diversity pattern in NPK + OM treatment was distinctly influenced by the contents of organic matter, available P and K. These results suggested that long-term fertilization had significant influences on merA gene diversity, which could be helpful to understand the Hg reduction process and potentially serve microbial remediation of Hg contaminated soil.

Responses of activities, abundances and community structures of soil denitrifiers to short-term mercury stress.

The responses of activities, abundances and community structures of soil denitrifiers to mercury (Hg) stress were investigated through a short-term incubation experiment. Four soil treatments with different concentrations of Hg (CK, Hg25, Hg50, and Hg 100, denoted as 0, 25, 50, and 100 mg Hg/kg dry soil, respectively) were incubated for 28 days. Soil denitrification enzyme activity (DEA) was measured at day 3, 7 and 28. The abundances and community structures of two denitrification concerning genes, nirS (cd(1)-nitrite reductase gene) and nosZ (nitrous oxide reductase gene), were analyzed using real-time PCR and denaturing gradient gel electrophoresis (DGGE). Results showed that soil DEA was significantly stimulated in the treatments of Hg25 and Hg50 compared with others at day 7. Meanwhile, no difference in the abundances of soil nirS and nosZ was found between Hg spiked treatments and CK, except the lower abundance of nirS (P < 0.05) in the Hg added treatments compared with that in the CK at day 28. The community structures of denitrifiers based on nirS gene presented obvious change at day 7 along with the Hg additions, however, no variation was found in all treatments based on the nosZ gene. The results indicated that Hg (Hg25 and Hg50) had a strongly short-term stimulation on soil DEA, and nirS gene is more sensitive than nosZ gene to Hg stress.

Characterization of soil heavy metal contamination and potential health risk in metropolitan region of northern China.

Soil in metropolitan region suffers great contamination risk due to the rapid urbanization especially in developing countries. Beijing and Tianjin, together with their surrounding regions, form a mega-metropolitan region in northern China. To assess the soil environmental quality, a total of 458 surface soil samples were collected from this area. Concentrations of Cr, Cu, Pb, Zn, As, Cd, and Hg were analyzed and compared to the Chinese environmental quality standards for soil. Multivariate analysis was carried out to identify the possible sources and Geographic Information Systems techniques were applied to visualize the spatial data. It was found that the primary inputs of As were due to pedogenic sources, whereas Hg was mainly of anthropogenic source. Other elements including Cr, Cu, Pb, Zn, and Cd were from both lithogenic and anthropogenic origins. Health risk assessment based on the maximum heavy metal concentration indicated that As derived from sewage irrigation area can result in carcinogenic lifetime risk due to ingestion and/or dermal contact of soil. The potential non-carcinogenic risk for children is significant for Pb and the cumulative effect of multiple metals is of concern for children in the vicinity of mining site. The results increased our knowledge for understanding natural and anthropogenic sources as well as health risk for metals in metropolitan soil.

Effects of mercury on reproduction, avoidance, and heat shock protein gene expression of the soil springtail Folsomia candida.

Based on the Cambisols of Beijing (used as agricultural soils), toxicity tests were conducted to investigate the effects of mercury (Hg) on reproduction and avoidance of Folsomia candida (Hexapoda: Collembola), as well as the transcriptional responses of the hsp70 gene, under different Hg concentrations and at different exposure times. Results showed that the hsp70 gene of the springtail was the most sensitive parameter to soil Hg stress, with a half-maximal effective concentration (EC50) of 0.42 mg/kg. The EC50 values based on reproduction and avoidance tests were 9.29 and 3.88 mg/kg, respectively. The expression level of the hsp70 gene was significantly up-regulated when soil Hg concentration was over 0.25 mg/kg (lowest-observed-effect concentration [LOEC]). In addition, responses of this gene expression were strongly induced after 48 h exposure under 1 mg/kg soil Hg, which probably was due to the fast and sensitive response of the gene transcription to Hg stress. Thus, the results suggested that the responses of the hsp70 gene and individual-level effects (reproduction and avoidance) could be integrated to provide helpful information for environmental monitoring and assessment of contaminated soils.

Characterization of PAHs contamination in soils from metropolitan region of Northern China.

The present study characterized the distribution, sources as well as carcinogenic potency of PAHs in surface soil from metropolitan region of northern China. The total PAHs in topsoil ranged from 322.6 to 23244.7 microg kg(-1). The mean and median concentrations of 16 PAHs were 1040.8 and 626.7 microg kg(-1), respectively. Source analysis revealed that pyrogenic sources played a major role at the locations and pyrogenic PAHs were mainly from incomplete combustion of coal, biomass and petroleum. The calculated mean BaP-equivalent values for individual carcinogenic PAHs were 148.4 microg kg(-1).

Microbes influence the fractionation of arsenic in paddy soils with different fertilization regimes.

Sequential extraction procedures were used to investigate the influence of the microbes on the distribution of arsenic in a Chinese paddy soil under different long-term fertilization treatments. The paddy soil with four long-term fertilization treatments (CK, M, NPK and NPK+M) and three levels of arsenate addition (0, 50, 100 mg As kg(-1) dry soil), were selected to construct microcosms for laboratory incubation. After the incubation, soil samples were sequentially extracted to determine As in various fractions, i.e. water soluble (F0), exchangeable (F1), bound to carbonates (F2), bound to Fe and Mn oxides (F3), bound to organic matter and sulfides (F4), and residual (F5, mineral matrix). Results showed that most of the As was fixed by mineral matrix (F5, ratios ranging from 46.22% to 96.37%), followed by As bound to Fe and Mn oxides (F3, ratios ranging from 3.14% to 28.18%), and the ratios of the other four fractions (F0, F1, F2 and F4) were mostly less than 10%. The microbes in the paddy soil could make As transform from inactive fraction (F5) to relatively active fractions (F0, F1, F2 and F3) and thus increase its environmental risk. With the increase of the As addition levels and with the application of manure or chemical NPK fertilizers, As was distributed more in the relatively active fractions (F0, F1, F2, F3 and F4) in the paddy soil mediated by the microbes. In addition, Fe and Mn oxides could play an important role in decreasing the As leaching potential from the mineral matrix to soil solution and thus abate the As risk to human health.

Assessing the health risk of heavy metals in vegetables to the general population in Beijing, China.

A systematic survey of As, Cd, Cr, Cu, Ni, Pb and Zn concentrations in vegetables from 416 samples (involving 100 varieties) in Beijing was carried out for assessing the potential health risk to local inhabitants. The results indicated that the metal concentrations in vegetables ranged from < 0.001 to 0.479 microg/g fresh weight (fw) (As), < 0.001 to 0.101 microg/g fw (Cd), < 0.001 to 1.04 microg/g fw (Cr), 0.024 to 8.25 microg/g fw (Cu), 0.001 to 1.689 microg/g fw (Ni), < 0.001 to 0.655 microg/g fw (Pb) and 0.01 to 25.6 microg/g fw (Zn), with average concentrations of 0.013, 0.010, 0.023, 0.51, 0.053, 0.046 and 2.55 microg/g fw, respectively. The results showed that the concentrations of As, Cr, Cu, Cd, Pb and Ni in vegetables from open-fields were all significantly higher than those grown in greenhouses. In addition, in local-produced vegetables, all HMs except Zn were significantly higher than those in provincial vegetables. The estimated daily intake (EDI) of As, Cd, Cr, Cu, Ni, Pb and Zn from vegetables was 0.080, 0.062, 0.142, 3.14, 0.327, 0.283 and 15.7 microg/(kg body weight (bw) x d) for adults, respectively. Arsenic was the major risk contributor for inhabitants since the target hazard quotient based on the weighted average concentration (THQw) of arsenic amounted to 44.3% of the total THQ (TTHQ) value according to average vegetable consumption. The TTHQ was lower than 1 for all age groups, indicating that it was still safe for the general population of Beijing to consume vegetables.

Effects of environmental governance in mining areas: The trend of arsenic concentration in the environmental media of a typical mining area in 25 years.

China has considered different environmental management measures (EMMs) in mining areas. However, their effects remain unclarified. In this study, the achievements and limitations of different EMMs of a typical mining area-Huangchang realgar mine-located in Hunan province were explored. The variations in the arsenic concentrations in the soil, agricultural products, drinking water, and atmosphere in 25 years of EMM implementation were investigated. Source control measures, such as ceasing mining and smelting activities, disposal of waste residues, and purifying wastewater, significantly reduced the arsenic concentrations in the atmosphere and surface water by more than 99%-from 68 µg m-3 and 0.42 mg L-1 to 3.63 ng m-3 and 4.31 µg L-1, respectively. The arsenic concentrations in agricultural products decreased by more than 78.8%-from 1.32 mg kg-1 in wheat to 0.28 mg kg-1 in vegetable and 0.13 mg kg-1 in maize-after the planting structure adjustment (PSA). However, the chronic daily intake of arsenic via product ingestion was 1.5 times higher than the benchmark dose lower confidence limit. Natural attenuation measures exerted limited effects on soil remediation; the arsenic concentration in the soil decreased insignificantly from 291.9 mg kg-1 to 213.3 mg kg-1. With the current attenuation rate, decreasing the soil arsenic concentration to under 30 mg kg-1 would require 47,900 years. The exceeding contaminant concentration in the resuspended dust, surface runoff, and agricultural products from the contaminated soil must be considered. China's EMMs in mining areas have achieved significant results, but the contaminated soil requires more attention and the PSA should accommodate the dietary habits and economic limits.

Mercury in soils of three agricultural experimental stations with long-term fertilization in China.

Mercury (Hg) in the agricultural ecosystem is a global concern because of its high potential toxicity. The objectives of this study were to determine the concentration and distribution of Hg in soils from three long-term experimental stations, i.e., Taoyuan (TY) and Qiyang (QY) in Hunan Province and Fengqiu (FQ) in Henan Province of China, and thus to assess the possible food and health risk of long-term applications of fertilizers. Soil samples at each site were collected from different fertilization plots and also from soil profiles with depths 0-100 cm. There were significant differences in soil Hg concentrations in 0-20 cm (A) or 20-40 cm (B) horizon among the three experimental stations. QY station showed significantly higher Hg concentrations than TY and FQ stations. However, there were no significant differences in soil Hg concentrations between A and B horizons at each station. It was concluded that the soil Hg concentrations at the three sites were mainly controlled by the parent materials. Moreover, chemical fertilizer, especially phosphorous fertilizers, could influence the soil Hg concentrations to some extent at the station with lower soil Hg concentrations, for example, at TY station. There were minimal amounts of Hg resulting from applications of the other chemical fertilizers and organic manure, and thus the fertilization had very low risk to the food security of the agro-ecosystems in the terms of Hg inputs and contamination.

Method for identifying outliers of soil heavy metal data.

Artificial errors in the experimental process may lead to some outliers, which reduce data quality and cause erroneous judgment in soil pollution assessment. Based on this, a method for detecting outliers of soil heavy metal data was proposed in this study. The As, Cd, and Pb concentrations of the soil in Beijing, China, were taken as samples to verify the validity of the method. Results showed that there were 8, 34, and 38 outliers for the As, Cd, and Pb concentrations in the Beijing soil, respectively. The result of re-analyzed revealed that 75.0, 76.5, and 92.1% of the As, Cd, and Pb outliers, respectively, were caused by artificial errors. After correcting, the interpolation accuracy for data was improved significantly. The mean relative error (MRE) of the As, Cd, and Pb outliers decreased by 48.0, 44.6, and 54.7%, while the mean square error of these outliers decreased by 34.2, 33.3, and 46.4%, respectively. The MRE values of the nearest neighboring points which were influenced by the outliers decreased by 5.2, 20.6, and 27.6%, while the mean square error of these points decreased by 5.3, 17.3, and 33.2%, respectively. To our knowledge, this is the first study on detecting outliers of soil heavy metal data. The method considers both spatial and numerical outliers, which avoids the limitation of single method, and can effectively improve the data quality of soil heavy metal concentrations with a finite sample size and analysis time.

Landscape ecology of the Guanting Reservoir, Beijing, China: multivariate and geostatistical analyses of metals in soils.

Surface soil samples were collected from 52 sites around Guanting Reservoir in Beijing, China, and contents of 'total' metals (Cu, Zn, Pb, Cd, As, Ni and Cr) were determined. The results indicate that the degree of heavy metal pollution in the soils declines in the order of Cd>Cr>Zn>As>Cu>Ni>Pb. Based on the results of a combination of multivariate statistics and geostatiscal analyses, it was concluded that land application of phosphate fertilizer, wastewater and sludge were the primary sources of Cd and Zn in soils. Whereas As, Cu, Cr and Ni in some soils were due to natural rock weathering. The sources of Pb in soils only partially originated from land application of phosphate fertilizer, but mainly from vehicle exhaust. The greatest concentrations of all metals, except for Pb, were found in Huailai County and the towns of Yanghedaqiao and Guanting.