Comparison of the common spatial interpolation methods used to analyze potentially toxic elements surrounding mining regions

Spatial distribution and risk assessment of perchlorate in raw cow milk from China

Perchlorate is a ubiquitous environmental contaminant worldwide, recognized as an emerging thyroid toxicant. This study focused on the pollution status, spatial distribution, possible sources of perchlorate in raw cow milk collected from 155 dairy farms in China, as well as the health risk of exposure to perchlorate through dairy products. The results showed that the detection rate of perchlorate in raw milk was 100% with the mean of 15.9 μg/kg, indicating the ubiquitous contamination of perchlorate in raw milk from China. The simulation of spatial distributions indicated that the levels of perchlorate in raw milk were spatially correlated, and relatively high levels of perchlorate exist in certain parts of Beijing-Tianjin-Hebei, Shanxi, Henan, and Zhejiang, suggesting potential environmental perchlorate contamination in these regions. A positive correlation was found between the perchlorate level in milk and the perchlorate level in feed, indicating the transfer process of perchlorate from feed to milk. The hazard quotient (HQ) values of exposure to perchlorate by dairy products ranged between 0 and 2.14, with the mean of 0.0188 and P95 of 0.101, indicating relatively low health risk to perchlorate through dairy products. To our knowledge, this is the first nationwide study on the spatial distribution and risk assessment of perchlorate in raw cow milk from China.

Screening and optimization of interpolation methods for mapping soil-borne polychlorinated biphenyls

There is yet no scientific consensus, and for now, on how to choose the optimal interpolation method and its parameters for mapping soil-borne organic pollutants. Take the polychlorinated biphenyls (PCBs) for instance, we present the comparison of some classic interpolation methods using a high-resolution soil monitoring database. The results showed that empirical Bayesian kriging (EBK) has the highest accuracy for predicting the total PCB concentration, while root mean squared error (RMSE) in inverse distance weighting (IDW) is among the highest in these interpolation methods. The logarithmic transformation of non-normally distributed data contributed to enhance considerably the semivariogram for modeling in kriging interpolation. The increasing of search neighborhood reduced IDW's RMSE, but slightly affected in ordinary kriging (OK), while both of them resulted in over smooth of prediction map. The existence of outliers made the difference between two points increase sharply, and thereby weakening spatial autocorrelation and decreasing the accuracy. As predicted error increased continuously, the prediction accuracy of different interpolation methods reached unanimity gradually. The attempt of the assisted interpolation algorithm did not significantly improve the prediction accuracy of the IDW method. This study constructed a standardized workflow for interpolation, which could reduce human error to reach higher interpolation accuracy for mapping soil-borne PCBs.

Exploring the accumulation capacity of dominant plants based on soil heavy metals forms and assessing heavy metals contamination characteristics near gold tailings ponds

Heavy metal contamination of soil commonly accompanies problems around gold mine tailings ponds. Fully investigating the distribution characteristics of heavy metals and the survival strategies of dominant plants in contaminated soils is crucial for effective pollution management and remediation. This study aims to investigate the contamination characteristics, sources of heavy metals (As, Cd, Pb, Hg, Cu, Zn, Cr, and Ni) in soils around gold mine tailings ponds areas (JHH and WZ) and to clarify the form distribution of heavy metals (As, Cd, Pb, Hg) in contaminated plots as well as their accumulation and translocation in native dominant plants. The results of the study showed that the concentrations of As, Pb, Cd, Cu, and Zn in soil exceeded the national limits at parts of the sampling sites in both study areas. The Nemerow pollution index showed that both study areas reached extreme high pollution levels. Spatial analysis showed that the main areas of contamination were concentrated around metallurgical plants and tailings ponds, with Cd exhibiting the most extensive area of contamination. In the JHH, As (74%), Cd (66%), Pb (77%), Zn (47%) were mainly from tailings releases, and Cu (52%) and Hg (51%) were mainly from gold ore smelting. In the WZ, As (42%), Cd (41%), Pb (73%), Cu (47%), and Zn (41%) were mainly from tailings releases. As, Cd, Pb, and Hg were mostly present in the residue state, and the proportion of water-soluble, ion-exchangeable, and carbonate-bound forms of Cd (19.93%) was significantly higher than that of other heavy metals. Artemisia L. and Amaranthus L. are the primary dominating plants, which exhibited superior accumulation of Cd compared to As, Pb, and Hg, and Artemisia L. demonstrated a robust translocation capacity for As, Pb, and Hg. Compared to the concentrations of other forms of soil heavy metals, the heavy metal content in Artemisia L correlates significantly better with the total soil heavy metal concentration. These results offer additional systematic data support and a deeper theoretical foundation to bolster pollution-control and ecological remediation efforts in mining areas.

Exploring the Spatial Pattern of Urban Forest Ecosystem Services based on i-Tree Eco and Spatial Interpolation: A Case Study of Kyoto City, Japan

Urban forest, as an essential urban green infrastructure, is critical in providing ecosystem services to cities. To enhance the mainstreaming of ecosystem services in urban planning, it is necessary to explore the spatial pattern of urban forest ecosystem services in cities. This study provides a workflow for urban forest planning based on field investigation, i-Tree Eco, and geostatistical interpolation. Firstly, trees across an array of land use types were investigated using a sampling method. Then i-Tree Eco was applied to quantify ecosystem services and ecosystem service value in each plot. Based on the ecosystem services estimates for plots, four interpolation methods were applied and compared by cross-validation. The Empirical Bayesian Kriging was determined as the best interpolation method with higher prediction accuracy. With the results of Empirical Bayesian Kriging, this study compared urban forest ecosystem services and ecosystem service value across land use types. The spatial correlations between ecosystem service value and four types of point of interest in urban places were explored using the bivariate Moran's I statistic and the bivariate local indicators of spatial association. Our results show that the residential area in the built-up area of Kyoto city had higher species richness, tree density, ecosystem services, and total ecosystem service value. Positive spatial correlations were found between ecosystem service value and the distribution of urban space types including the tourist attraction distribution, urban park distribution, and school distribution. This study provides a specific ecosystem service-oriented reference for urban forest planning based on land use and urban space types.

Optimization of spatial prediction and sampling strategy of site contamination based on Thiessen polygon coupling interpolation

Contaminated sites pose a serious threat to the ecological environment and human health. Because of the presence of multiple peaks in the pollution data of some contaminated sites, as well as strong spatial heterogeneity and skewness in their distribution, the accuracy of spatial interpolation prediction is low. This study proposes a method for investigating highly skewed contaminated sites, which uses Thiessen polygons coupled with geostatistics and deterministic interpolation to optimize the spatial prediction and sampling strategy of sites. An industrial site in Luohe is used as an example to validate the proposed method. The results indicate that using 40 × 40 m as the minimum initial sampling unit can obtain data that is representative of the regional pollution situation. Evaluation indexes reveal that the ordinary kriging (OK) method for interpolation prediction accuracy and the radial basis function_inverse distance weighted (RBF_IMQ) method for pollution scope prediction provides the best results, which can effectively improve the spatial prediction accuracy of pollution in the study area. Each accuracy indicator is enhanced by 20-70% after supplementing 11 sampling points in the suspect region, and the identification of the pollution scope approaches 95%. This method offers a novel approach for investigating highly biased contaminated sites, which can optimize the spatial prediction accuracy of pollution and reduce economic costs.

Research on Risk Assessment and Contamination Monitoring of Potential Toxic Elements in Mining Soils

Potentially toxic element (PTE) contamination in soils has serious impacts on ecosystems. However, there is no consensus in the field of assessment and monitoring of contaminated sites in China. In this paper, a risk assessment and pollution monitoring method for PTEs was proposed and applied to a mining site containing As, Cd, Sb, Pb, Hg, Ni, Cr, V, Zn, Tl, and Cu. The comprehensive scoring method and analytical hierarchical process were used to screen the priority PTEs for monitoring. The potential ecological risk index method was used to calculate the risk index of monitoring point. The spatial distribution characteristics were determined using semi-variance analysis. The spatial distribution of PTEs was predicted using ordinary kriging (OK) and radial basis function (RBF). The results showed that the spatial distribution of As, Pd, and Cd are mainly influenced by natural factors, while Sb and RI are influenced by both natural and human factors. OK has higher spatial prediction accuracy for Sb and Pb, and RBF has higher prediction accuracy for As, Cd, and RI. The areas with high ecological risk and above are mainly distributed on both sides of the creek and road. The optimized long-term monitoring sites can achieve the monitoring of multiple PTEs.

Effects of anthropogenic and natural environmental factors on the spatial distribution of trace elements in agricultural soils

The concentrations of trace elements in agricultural soils directly affect the ecological security and quality of agricultural products. A comprehensive study aimed at quantitatively analyze the effects of anthropogenic and natural environmental factors on the spatial distribution of heavy metals (HMs) and selenium (Se) in agricultural soils in a typical grain producing area of China. Factors considered in this study were parent rock, soil physicochemical properties, topography, precipitation, mine activity, and vegetation. Results showed that the median values of Zn, Cd, Cr, and Cu of 111 topsoil samples exceeded the background values of Guangxi province but were lower than the relevant national soil quality standards, and 85% of soil samples were classified as having rich Se levels (0.40 -3.0 mg kg^-1). The potential ecological risk index of soil heavy metals as a whole was low, with Cd in 9% of the samples posing moderate ecological risk. The concentrations of heavy metals and Se were relatively high in soils from shale rock. Soil properties, mainly Fe₂O₃ and Mn played a dominant role on soil HMs and Se concentrations. Based on GeoDetector, we found that the interaction effects of two factors on the spatial differentiation of soil HMs and Se were greater than their sum effect. Among the factors, Mn enhanced the explanatory power of the model the most when interacting with other factors for soil Zn; the greatest interactive effect was between distance from mining area and Mn for Cd (q = 0.70); Fe₂O₃ significantly promoted the spatial differentiation of soil Cr, Cu and Se when interacting with other factors (q > 0.50). These findings contribute to a better understanding of the factors that drive the distribution of HMs and Se in agricultural soils.

Spatial and temporal distribution and influencing factor analysis of the malignant tumor mortality rate around the mining area

Mining activities can threaten residents' health even lives. Integrating spatial empirical Bayesian smoothing, joinpoint regression and spatiotemporal scanning methods, we analyzed aggregations and possible factors of four tumor mortality rates at township and village scales from 2012 to 2016 in Suxian district of Hunan Province, China. Results indicate: (1) Mortality rates were ranked: lung cancer > liver cancer > gastric cancer > colorectal cancer. (2) Lung cancer had a higher five-year mortality rate in the middle; relative risk (RR) of death from lung cancer from 2012 to 2015 in Xujiadong Village was 7.48. Liver cancer had a higher five-year mortality rate in the Middle West; RR in areas centered on Nanta Street with a radius of 9.87 km from 2015 to 2016, was 1.83. Gastric cancer had a higher five-year mortality rate in the east; RR in Xujiadong Village from 2012 to 2014 was 6.9. Five-year mortality rate of colorectal cancer was higher in the northwest; RR in regions centered on Huangcao Village with a radius of 12.11 km in 2016, was 2.88. (3) Pollution from ore mining and smelting, heavy metal and non-metallic, and mine transportation were the main possible factors. This research provides a method reference for studying spatiotemporal patterns of disease in China even the world.

Classification of air quality zones and fine particulate matter sensitive areas by risk assessment approach

Many studies have shown that fine particulate matter can cause health problems. Thus, effectively controlling fine particulate matter concentration is an important issue around the world. The Taiwan Environmental Protection Administration (TWEPA) divides Taiwan into seven air quality zones based on counties and cities for managing air quality and analyzing pollution transmission. However, this artificial division by administrative areas relatively poorly match natural conditions and topographical and geographic factors and hence poorly represent air quality characteristics. This study proposes an air quality sensitive map analysis framework, which uses hierarchical agglomerative clustering with empirical orthogonal function and analysis of variance methods, to provide more detailed, reasonable, and township-level air quality zones incorporating the different spatial-temporal characteristics over the region. The risk concept is introduced to evaluate PM_2.5 risk sensitivity for each administrative district, combining three aspects: hazard (PM_2.5 exceedance probability), exposure (population density of sensitive groups), and vulnerability (average wind speed). Considering air quality spatial-temporal characteristics, Taiwan can be optimally divided into 14 air quality zones. PM_2.5 risk is highest for western inland towns than western coastal towns, with eastern regions exhibiting least risk. Adopting the proposed air quality zones and clarifying high risk areas allows PM_2.5 causes to be identified for different air quality zones. This allows a targeted control strategy for high risk areas to effectively improve domestic air quality. The proposed model also provides powerful reference for environmental management and environmental impact assessment for future construction and development.

Optimizing spatial interpolation method and sampling number for predicting cadmium distribution in the largest shallow lake of North China

Cadmium (Cd) pollution has been widely recognized in lake ecosystems. Although the accurate prediction of the spatial distributions of Cd in lakes is important for controlling Cd pollution, the traditional monitoring methods of setting discrete and limited sampling points cannot actually reflect the continuous spatial distribution characteristics of Cd. In this study, we set up 93 sampling points in Baiyangdian Lake (BYDL), and collected surface water, overlying water and sediment samples from each sampling point. Cd contents were measured to predict their spatial distributions in different environmental components by three interpolation methods, inverse distance weighted (IDW), radial basis function (RBF) and ordinary kriging (OK), and the effects of different sampling numbers on the interpolation accuracy were also assessed to optimize the interpolation method and sampling number. The results showed that the interpolation accuracy of IDW decreased with increasing power values. The best basis function for RBF was IMQ, and the best semivariogram models for OK were the spherical model and stable model. The best interpolation method for the waters and sediments was RBF-IMQ compared with OK and IDW. Within the sampling number range of 50-93, the interpolation accuracy for Cd in surface water increased with the increase in sampling number. Comparatively, the interpolation accuracy was the highest for overlying water and sediments when the sampling number was 60. The findings of this work provide a combined sampling and spatial interpolation method for monitoring the spatial distribution and pollution levels of Cd in the waters and sediments of shallow lakes.

Long-term nonlinear relationship between PM2.5 and ten leading causes of death

Air pollution has become a major concern worldwide. Many epidemiological studies have proved relationships between fine particulate matter (PM_2.5) and various diseases, but most studies only use short-term and models for specific groups to derive relationships with acute diseases. This makes it difficult to understand long-term exposure, nonlinear relationships, and spatial-temporal health risks regarding chronic diseases. Therefore, this study proposed to analyze and map PM_2.5 exceedance probability from long-term spatial-temporal monitoring data using radial basis function estimation. We then constructed and compared multiple linear regression and generalized additive models to investigate linear and nonlinear relationships between long-term average PM_2.5 concentration, PM_2.5 potential probability for exceeding the standard, and standardized mortality for the top ten causes of death in all towns and villages in Taiwan nationally from 2010 to 2017. Linear models indicate that increasing PM_2.5 concentration increased malignant neoplasm, pneumonia, and chronic lower respiratory disease mortalities; chronic liver diseases; and cirrhosis; whereas heart diseases and esophagus cancer mortality decreased. For the nonlinear model results, it can be found that there were also significant nonlinear relationships between PM_2.5 concentration and malignant mortalities for neoplasm, heart disease, diabetes; and trachea, bronchus, lung, liver, intrahepatic bile duct, and esophagus cancer. Thus, long-term exposure to PM_2.5 may be a significant risk factor for multiple acute and chronic diseases. Results from this study can be directly applied worldwide to provide air quality and health management references for governments, and important information on long-term health risks for local residents in the study area.

Mapping the impact of a large municipal waste disposal area on surface water: 1993-2017, case of Laogang, Shanghai

In China, the impact of waste disposal facilities is always a cause of concern for the government and the public. Laogang Municipal Waste Disposal Area (LMDA), Shanghai, one of the largest municipal waste disposal areas in the world was selected as case in this study, and it was attempted to analyze the changes in the surface water quality, and map the impacted area by LMDA on surrounding streams from its operation period of 1993-2017. The results showed that, during the whole period, only biochemical oxygen demand (BOD₅) showed a continuous improvement with a percentage of 85.92%, however, chemical oxygen demand (COD_cr), ammonia (NH₄⁺-N) and total phosphorus (TP) significantly improved but BOD₅ slightly deteriorated began from 2013. Using spatial analysis tools and Kendall's concordance test, COD_cr and phenol at LMDA showed a significant impact on surrounding surface water; especially, the impacted area for COD_cr decreased from 106.30 km² to 22.86 km² from 1993 to 2017, which dropped from 4.3 to 0.9 times the area of LMDA. Surprisingly, NH₄⁺-N and TP at LMDA were affected by the surrounding streams, instead of having an impact on them. Interestingly, heavy metals and non-metals such as Hg, As, Zn, and Se in the surrounding streams were unlikely affected by LMDA. The driving forces for surface water quality improvement included the eco-remediation of closed unsanitary landfills, upgrade in waste shipping and terminals, operation of sanitary landfills and incineration plants for landfill diversion. Capsule: Impacted area of municipal waste disposal area is not so large.

A novel interpolation method to predict soil heavy metals based on a genetic algorithm and neural network model

To improve the prediction accuracy of soil heavy metals (HMs) by spatial interpolation, a novel interpolation method based on genetic algorithm and neural network model (GANN model), which integrates soil properties and environmental factors, was proposed to predict the soil HM content. Eleven soil HMs (Cu, Pb, Zn, Cd, Ni, Cr, Hg, As, Co, V and Mn) were predicted using the GANN model. The results showed that the model had a good prediction performance with correlation coefficients (R²) varying from 0.7901 to 0.9776. Compared with other traditional interpolation methods, including inverse distance weighting (IDW), ordinary kriging (OK), universal kriging (UK), and spline with barriers interpolation (SBI) methods, the GANN model had a relatively lower root mean square error value, ranging from 0.0497 to 77.43, suggesting that the GANN model might be a more accurate spatial interpolation method and the soil properties together with the environmental geographical factors played key roles in prediction of soil HMs.

Effects of Driving Factors on Forest Aboveground Biomass (AGB) in China’s Loess Plateau by Using Spatial Regression Models

Forests are the main body of carbon sequestration in terrestrial ecosystems and forest aboveground biomass (AGB) is an important manifestation of forest carbon sequestration. Reasonable and accurate quantification of the relationship between AGB and its driving factors is of great importance for increasing the biomass and function of forests. Remote sensing observations and field measurements can be used to estimate AGB in large areas. To explore the applicability of the panel data models in AGB and its driving factors, we compared the results of panel data models (spatial error model and spatial lag model) with those of geographically weighted regression (GWR) and ordinary least squares (OLS) to quantify the relationship between AGB and its driving factors. Furthermore, we estimated the tree height, diameter at breast height, canopy cover (CC) and species diversity index (Shannon–Wiener index) of Robinia pseudoacacia plantations in Changwu on the Loess Plateau using field data and remote sensing images by a random forest model and estimated soil organic carbon (SOC) contents using laboratory data by ordinary kriging (OK) interpolation. We estimated AGB using the already estimated tree height and diameter at breast height combined with the allometric growth equation. In this study, we estimated SOC contents by OK interpolation, and the accuracy R2 values for each soil layer were greater than 0.81. We estimated diameter at breast height (DBH), CC, SW and tree height (TH) using the random forest, and the accuracy R2 values were 0.85, 0.82, 0.76 and 0.68, respectively. We estimated AGB with random forest and the allometric growth equation and found that the average AGB was 55.80 t/ha. The OLS results showed that the residuals of the OLS regression exhibited obvious spatial correlations and rejected OLS applications. GWR, SEM and SLM were used for spatial regression analysis, and SEM was the best model for explaining the relationship between AGB and its driving factors. We also found that AGB was significantly positively correlated with CC, SW, and 0–60 cm SOC content (p < 0.05) and significantly negatively correlated with slope aspect (p < 0.01). This study provides a new idea for studying the relationship between AGB and its driving factors and provides a basis for practical forest management, increasing biomass, and giving full play to the role of carbon sequestration.

Reconstructing High-Precision Coral Reef Geomorphology from Active Remote Sensing Datasets: A Robust Spatial Variability Modified Ordinary Kriging Method

Active remote sensing technology represented by multi-beam and lidar provides an important approach for the effective acquisition of underwater coral reef geomorphological information. A spatially continuous surface model of coral reef geomorphology reconstructed from active remote sensing datasets can provide important geomorphological parameters for the research of coral reef geomorphological and ecological changes. However, the surface modeling methods commonly used in previous studies, such as ordinary kriging (OK) and natural neighborhood (NN), often represent a “smoothing effect”, which causes the strong spatial variability of coral reefs to be imprecisely reflected by the reconstructed surfaces, thus affecting the accurate calculation of subsequent geomorphological parameters. In this study, a spatial variability modified OK (OK-SVM) method is proposed to reduce the impact of the “smoothing effect” on the high-precision reconstruction of the complex geomorphology of coral reefs. The OK-SVM adopts a collaborative strategy of global parameter transformation, local residual correction, and extremum correction to modify the spatial variability of the reconstructed model, while maintaining high local accuracy. The experimental results show that the OK-SVM has strong robustness to spatial variability modification. This method was applied to the geomorphological reconstruction of the northern area of a coral atoll in the Nansha Islands, South China Sea, and the performance was compared with that of OK and NN. The results show that OK-SVM has higher numerical accuracy and attribute accuracy in detailed morphological fidelity, and is more adaptable in the geomorphological reconstruction of coral reefs with strong spatial variability. This method is relatively reliable for achieving high-precision reconstruction of complex geomorphology of coral reefs from active remote sensing datasets, and has potential to be extended to other geomorphological reconstruction applications.

An enhanced dual IDW method for high-quality geospatial interpolation

Many geoscience problems involve predicting attributes of interest at un-sampled locations. Inverse distance weighting (IDW) is a standard solution to such problems. However, IDW is generally not able to produce favorable results in the presence of clustered data, which is commonly used in the geospatial data process. To address this concern, this paper presents a novel interpolation approach (DIDW) that integrates data-to-data correlation with the conventional IDW and reformulates it within the geostatistical framework considering locally varying exponents. Traditional IDW, DIDW, and ordinary kriging are employed to evaluate the interpolation performance of the proposed method. This evaluation is based on a case study using the public Walker Lake dataset, and the associated interpolations are performed in various contexts, such as different sample data sizes and variogram parameters. The results demonstrate that DIDW with locally varying exponents stably produces more accurate and reliable estimates than the conventional IDW and DIDW. Besides, it yields more robust estimates than ordinary kriging in the face of varying variogram parameters. Thus, the proposed method can be applied as a preferred spatial interpolation method for most applications regarding its stability and accuracy.

Effectiveness of predicting spatial contaminant distributions at industrial sites using partitioned interpolation method

Soil pollution at industrial sites is an important issue in China and in most other regions of the world. The accurate prediction of the spatial distribution of pollutants at contaminated industrial sites is a requirement for the development of most soil remediation strategies, and is commonly performed using spatial interpolation methods. However, significant and abrupt variations in the spatial distribution of pollutants decrease prediction accuracy. During this study, the use of partition interpolation methods was applied to benzo fluoranthene in four soil layers at a contaminated site to determine their ability to improve prediction accuracy in comparison to unpartitioned methods. The examined methods for partitioned interpolation included inverse distance weighting (IDW), radial basis function (RBF), and ordinary kriging (OK). The prediction results of the three methods for partitioned interpolation were compared, and the applicability of partition interpolation was determined. The prediction error associated with the partitioned interpolation methods decreased by 70% compared to unpartitioned interpolation. The prediction accuracy of IDW-based partition interpolation was higher than that of RBF- and OK-based partition interpolation techniques, and it was suitable for identification of highly polluted areas. Partition interpolation is also applicable to 12 other PAHs controlled by USEPA that can be detected, and the prediction effects could also verify this interpolation choice. In addition, the results also demonstrated that the more the maximum concentration deviated from the "norm", the greater the prediction error was caused by the smoothing effects of the interpolation models. These results suggest that the partition interpolation with IDW method can be effectively used to obtain relatively accurate spatial contaminant distribution information, and to identify highly polluted areas.

Comparison of common spatial interpolation methods for analyzing pollutant spatial distributions at contaminated sites

Accurate prediction of the spatial distribution of pollutants in soils based on applicable interpolation methods is often the basis for soil remediation in contaminated sites. However, the applicable interpolation method has not been determined for contaminated sites due to the complex spatial distribution characteristics and stronger local spatial variability of pollutants. In this research, the prediction accuracies of three interpolation methods (including the different values of their parameters) for the spatial distribution of benzo[b]fluoranthene (BbF) in four soil layers were compared. These included inverse distance weighting (IDW), radial basis function (RBF), ordinary kriging (OK). The results indicated: (1) IDW1 is applicable for the first layer, RBF-IMQ is applicable to the second, third, and fourth layers. (2) For IDW, the prediction error is bigger with high weight where high values and low values intersect, while the prediction error is smaller where high (or low) values aggregated distribution. (3) For RBF, if the pollutant concentration trend at the predicted location is consistent with the known points in its neighborhood, the prediction accuracy is higher. (4) IDW is suitable for fitting more drastic curved surfaces, while RBF is more effective for relatively gentle curved surfaces and OK is reasonable for curved surfaces without local outliers. (5) The interpolation uncertainty is positively associated with the contaminant concentration and local spatial variability. Therefore, we suggest the selection of the applicable interpolation model must be based on the principle of the model and the spatial distribution characteristics of the pollutants.

Distribution Characteristics and Pollution Assessment of Soil Heavy Metals under Different Land-Use Types in Xuzhou City, China

Xuzhou, as a mining city in China, has been experiencing 130 years of coal mining and processing. To explore the spatial distribution characteristics and pollution status of soil heavy metals (Cr, Cd, As, Hg, Zn, and Pb) under different land-use types, a total of 2697 topsoil samples were collected in all of the areas (except for water) of Xuzhou in 2016. Overall, the mean concentrations of Cr (70.266 mg/kg), Cd (0.141 mg/kg), As (10.375 mg/kg), Hg (0.036 mg/kg), Zn (64.788 mg/kg), and Pb (24.84 mg/kg) in Xuzhou soils were lower than the environmental quality standard for soils (GB15618-1995). However, the mean concentrations of Cr, Hg, and Pb exceeded their corresponding background values, with the mean concentration of Hg being almost three times its background value. For different land-use types, the highest mean concentration of Cr was concentrated in grassland soils. The mean concentrations of Cd, As, Zn, and Pb in mining area soils were higher than those in the other soils. The mean concentration of Hg was the highest in the built-up area soils. Based on the potential ecological risk assessment, the forestland, garden land, grassland, and others were at low and moderate risk levels, the farmland and mining area were at low, moderate, and high risk levels, and the built-up area was at various risk levels in Xuzhou. There was a significant positive correlation between Cr, Pb, and Hg concentrations and the corresponding organic carbon contents in the farmland, built-up area, garden land, forestland, and other soils ( p < 0.01 ). A high degree of correlation was found between Cr and Hg concentrations, as well as organic carbon contents in grassland soils, with values of p < 0.05 and p < 0.01 , respectively. An obvious correlation could be seen between Hg concentrations and organic carbon contents in mining area soils ( p < 0.01 ).

Morphological Precision Assessment of Reconstructed Surface Models for a Coral Atoll Lagoon

In addition to remote-sensing monitoring, reconstructing morphologic surface models through interpolation is an effective means to reflect the geomorphological evolution, especially for the lagoons of coral atolls, which are underwater. However, which interpolation method is optimal for lagoon geomorphological reconstruction and how to assess the morphological precision have been unclear. To address the aforementioned problems, this study proposed a morphological precision index system including the root mean square error (RMSE) of the elevation, the change rate of the local slope shape (CRLSS), and the change rate of the local slope aspect (CRLSA), and introduced the spatial appraisal and valuation approach of environment and ecosystems (SAVEE). In detail, ordinary kriging (OK), inverse distance weighting (IDW), radial basis function (RBF), and local polynomial interpolation (LPI) were used to reconstruct the lagoon surface models of a typical coral atoll in South China Sea and the morphological precision of them were assessed, respectively. The results are as follows: (i) OK, IDW, and RBF exhibit the best performance in terms of RMSE (0.3584 m), CRLSS (51.43%), and CRLSA (43.29%), respectively, while with insufficiently robust when considering all three aspects; (ii) IDW, LPI, and RBF are suitable for lagoon slopes, lagoon bottoms, and patch reefs, respectively; (iii) The geomorphic decomposition scale is an important factor that affects the precision of geomorphologic reconstructions; and, (iv) This system and evaluation approach can more comprehensively consider the differences in multiple precision indices.