Kriging methods with auxiliary nighttime lights data to detect potentially toxic metals concentrations in soil

A hierarchical residual correction-based hyperspectral inversion method for soil heavy metals considering spatial heterogeneity

Promising hyperspectral remote sensing exhibits substantial potential in monitoring soil heavy metal (SHM) contamination. Nevertheless, the local spatial perturbation effects induced by environmental factors introduce considerable variability in SHM distribution. This engenders non-stationary relationship between SHM concentrations and spectral reflectance, posing challenges for accurate inversion of SHM globally. Addressing this gap, a novel Hierarchical Residual Correction-based Hyperspectral Inversion Method (HRCHIM) is proposed for SHM, considering their spatial heterogeneity. Initially, a global model is constructed using ground hyperspectral data to predict SHM concentration, capturing overarching contamination trends. Subsequently, four hierarchical levels, segmented by residual standard deviation (SD) intervals, identify critical environmental factors via Geodetector. These factors inform local residual correction models, refining global model predictions. HRCHIM aims to synergize global trends and local stochasticity to enhance prediction accuracy and interpretation of SHM spatial heterogeneity. Validated through a case study of a Cadmium(Cd)-contaminated mine area, six critical environmental factors were identified, exhibiting significant differences across hierarchical levels. By incorporating hierarchical correction models, HRCHIM demonstrated superior inversion performance compared to other conventional methods, achieving optimal prediction accuracies (Rv2 = 0.94, RMSEv = 0.21, and RPDv = 4.11). This innovative method can facilitate more precise and targeted strategies for preventing and controlling SHM contamination.

A comprehensive study of source apportionment, spatial distribution, and health risks assessment of heavy metal(loid)s in the surface soils of a semi-arid mining region in Matehuala, Mexico

BACKGROUND

This study investigates the contamination level, spatial distribution, pollution sources, potential ecological risks, and human health risks associated with heavy metal(loid)s (i.e., arsenic (As), copper (Cu), iron (Fe), manganese (Mn), lead (Pb), and zinc (Zn)) in surface soils within the mining region of Matehuala, located in central Mexico.

OBJECTIVES

The primary objectives are to estimate the contamination level of heavy metal(loid)s, identify pollution sources, assess potential ecological risks, and evaluate human health risks associated with heavy metal(loid) contamination.

METHODS

Soil samples from the study area were analysed using various indices including Igeo, Cf, PLI, mCd, EF, and PERI to evaluate contamination levels. Source apportionment of heavy metal(loid)s was conducted using the APCS-MLR and PMF receptor models. Spatial distribution patterns were determined using the most efficient interpolation technique among five different approaches. The total carcinogenic risk index (TCR) and total non-carcinogenic index (THI) were used in this study to assess the potential carcinogenic and non-carcinogenic hazards posed by heavy metal(loid)s in surface soil to human health.

RESULTS

The study reveals a high contamination level of heavy metal(loid)s in the surface soil, posing considerable ecological risks. As was identified as a priority metal for regulatory control measures. Mining and smelting activities were identified as the primary factors influencing heavy metal(loid) distributions. Based on spatial distribution mapping, concentrations were higher in the northern, western, and central regions of the study area. As and Fe were found to pose considerable and moderate ecological risks, respectively. Health risk evaluation indicated significant levels of carcinogenic risks for both adults and children, with higher risks for children.

CONCLUSION

This study highlights the urgent need for monitoring heavy metal(loid) contamination in Matehuala's soils, particularly in regions experiencing strong economic growth, to mitigate potential human health and ecological risks associated with heavy metal(loid) pollution.

Improved mapping of heavy metals in agricultural soils using machine learning augmented with spatial regionalization indices

The accurate spatial mapping of heavy metal levels in agricultural soils is crucial for environmental management and food security. However, the inherent limitations of traditional interpolation methods and emerging machine-learning techniques restrict their spatial prediction accuracy. This study aimed to refine the spatial prediction of heavy metal distributions in Guangxi, China, by integrating machine learning models and spatial regionalization indices (SRIs). The results demonstrated that random forest (RF) models incorporating SRIs outperformed artificial neural network and support vector regression models, achieving R2 values exceeding 0.96 for eight heavy metals on the test data. Hierarchical clustering for feature selection further improved the model performance. The optimized RF models accurately predicted the heavy metal distributions in agricultural soils across the province, revealing higher levels in the central-western regions and lower levels in the north and south. Notably, the models identified that 25.78 % of agricultural soils constitute hotspots with multiple co-occurring heavy metals, and over 6.41 million people are exposed to excessive soil heavy metal levels. Our findings provide valuable insights for the development of targeted strategies for soil pollution control and agricultural soil management to safeguard food security and public health.

Analysis of driving factors for potential toxic metals in major urban soils of China: a geodetetor-based quantitative study

Potential toxic metal (PTM) is hazardous to human health, but the mechanism of spatial heterogeneity of PTM at a macro-scale remains unclear. This study conducts a meta-analysis on the data of PTM concentrations in the soil of 164 major cities in China from 2006 to 2021. It utilizes spatial analysis methods and geodetector to investigate the spatial distribution characteristics of PTMs. The geographic information systems (GIS) and geodetector were used to investigate the spatial distribution characteristics of PTMs, assess the influence of natural factors (NFs) and anthropogenic factors (AFs) on the spatial heterogeneity of PTMs in urban soils, and identified the potential pollution areas of PTMs. The results indicated that the pollution levels of PTMs in urban soils varied significantly across China, with higher pollution levels in the south than in the north. Cd and Hg were the most severely contaminated elements. The geodetector analysis showed that temperature and precipitation in NFs and land use type in AFs were considered as the main influencing factors, and that both AF and NF together led to the PTM variation. All these factors showed a mutually enhancing pattern which has important implications for urban soil management. PTM high-risk areas were identified to provide early warning of pollution risk under the condition of climate change.

GIS- and Multivariate-Based Approaches for Assessing Potential Environmental Hazards in Some Areas of Southwestern Saudi Arabia

Soil contamination is a major issue that endangers the ecology in most countries. Total concentrations of As, Cd, Co, Cr, Cu, Mn, Ni, Pb, VFe, and Zn were determined by analyzing soil samples from 32 surface soil samples in southwest Saudi Arabia, including certain areas of Al-Baha. Kriging techniques were used to create maps of the distribution of metal. To assess the levels of soil contamination in the research area, principal component analysis (PCA), contamination factors (CF), and pollution load index were used. The results show the stable model gave the best fit to the As and Zn semivariograms. The circular model fits the Cd, Co, and Ni semivariograms the best, while the exponential model fits the Cr, V, and Fe semivariograms the best. For Ni and Pb, respectively, spherical and Gaussian models are fitted. The findings demonstrated two clusters containing different soil heavy metal concentrations. According to the data, there were two different pollution levels in the research region: 36.58% of it is strongly contaminated, while 63.41% of it has a moderate level of contamination (with average levels of these metals 5.28 ± 5.83, 0.81 ± 0.19, 18.65 ± 6.22, 45.15 ± 23.25, 60.55 ± 23.74, 972.30 ± 223.50, 33.45 ± 14.11, 10.05 ± 5.13, 84.15 ± 30.72, 97.40 ± 30.05, and 43,245.00 ± 8942.95 mg kg-1 for As, Cd, Co, Cr, Cu, Mn, Ni, Pb, V, Fe, and Zn, respectively). The research area's poor management practices are reflected in the current results, which raised the concentration of harmful elements in the soil's surface layers. Ultimately, the outcomes of pollution concentration and spatial distribution maps could aid in informing decision-makers when creating suitable heavy metal mitigation strategies.

Improving prediction of soil heavy metal(loid) concentration by developing a combined Co-kriging and geographically and temporally weighted regression (GTWR) model

The study of heavy metal(loid) (HM) contamination in soil using extensive data obtained from published literature is an economical and convenient method. However, the uneven distribution of these data in time and space limits their direct applicability. Therefore, based on the concentration data obtained from the published literature (2000-2020), we investigated the relationship between soil HM accumulation and various anthropogenic activities, developed a hybrid model to predict soil HM concentrations, and then evaluated their ecological risks. The results demonstrated that various anthropogenic activities were the main cause of soil HM accumulation using Geographically and temporally weighted regression (GTWR) model. The hybrid Co-kriging + GTWR model, which incorporates two of the most influential auxiliary variables, can improve the accuracy and reliability of predicting HM concentrations. The predicted concentrations of eight HMs all exceeded the background values for soil environment in China. The results of the ecological risk assessment revealed that five HMs accounted for more than 90% of the area at the "High risk" level (RQ ≥ 1), with the descending order of Ni (100%) = Cu (100%) > As (98.73%) > Zn (95.50%) > Pb (94.90%). This study provides a novel approach to environmental pollution research using the published data.

A remote sensing-based strategy for mapping potentially toxic elements of soils: Temporal-spatial-spectral covariates combined with random forest

The selection of predictor variables is a crucial issue in building a digital mapping model of potentially toxic elements (PTEs) in soil. Traditionally, the predictor variables for mapping models of soil PTEs have been chosen from sets of spatial parameters or spectral parameters derived from geographical environmental data. However, the enrichment of soil PTEs exhibits significant variations in both spatial and temporal dimensions, with the temporal dimension often being overlooked in the selection of predictor variables for digital mapping models. This limitation hampers the robustness and generalizability of the models. Therefore, multi-source geographical data were used in this study to determine three temporal indices for characterizing the enrichment process of soil PTEs in temporal dimensions, and additionally to construct the temporal-spatial-spectral (TSS) covariate combinations. The random forest (RF) algorithm was used to map soil PTEs at a regional scale. Results showed that: (1) When using spatial parameters or spectral parameters as predictor variables and measured Pb content as the dependent variable, the values of the model performance indicator RPIQ (ratio of performance to inter-quartile range) were 2.66 and 2.27, respectively. After incorporating the temporal parameters into the model input, values of RPIQ for the RF model reached 3.55 (using spatial-temporal covariates) and 3.21 (using spectral-temporal covariates), representing performance improvements of 33.46% and 41.41%, respectively. (2) The RF model constructed with the temporal-spatial-spectral covariates achieved satisfactory mapping accuracy (R2 = 0.85; RMSE = 0.80 mg kg-1; RPIQ = 4.09). (3) The soil Pb content in the western and northeastern regions was relatively high, while the remaining areas exhibited lower Pb levels, mainly due to industrial activities. (4) The mapping results of Pb obtained in this study were superior to other mapping methods, such as ordinary kriging, artificial neural networks, and multivariate linear regression methods. The soil PTE mapping technique employed in this study that combined TSS covariates with the RF provided an effective methodological approach for preventing soil pollution, controlling environmental risk, and improving soil management.

An interpolation method incorporating the pollution diffusion characteristics for soil heavy metals - taking a coke plant as an example

The existing spatial interpolation methods in the prediction of soil heavy metal distribution are generally based on spatial auto correlation theory, rarely considering the pollution patterns. By contrast, in polluted sites, heavy metals have a strong heterogeneity even within a very small area, which is not exactly in line with auto correlation theory. This contradiction may lead to inaccuracy in spatial prediction. Atmospheric diffusion and deposition are one of the main sources of soil heavy metal pollution caused by coal-related production activities. To improve the prediction accuracy, the diffusion patterns of pollutants were considered in this paper by integrating Geodetector, Co-Kriging (COK), and partition interpolation. Geodetector was used to identify the main driving factors of soil pollution, based on which, the main driving factors were used as covariates introduced into the interpolation method (COK). Specifically, the amount of particulate matter deposition obtained by a pollutant diffusion model (AERMOD) was used as a covariate. For comparison, the distances to quenching, coke oven, and ammonium sulfate section were also used as covariates. Compared with the Ordinary Kriging method, the method COK-AERMOD established here decreased the root mean square error values of As (2.05 reduced to 1.89), Cd (0.18 reduced to 0.16), Cr (19.07 reduced to 12.97), Cu (6.92 reduced to 4.72), Hg (0.32 reduced to 0.28), Ni (16.92 reduced to 16.10), Pb (18.29 reduced to 16.62), and Zn (159.68 reduced to 153.66). This method in this paper is informative for the interpolation of soil elements in contaminated areas with known pollution source and diffusion patterns.

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.

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.

Effectiveness of predicting the spatial distributions of target contaminants of a coking plant based on their related pollutants

The prediction accuracy of the spatial distribution of soil pollutants at a site is relatively low. Related pollutants can be used as auxiliary variables to improve the prediction accuracy. However, little relevant research has been conducted on site soil pollution. To analyze the prediction accuracy of target pollutants combined with auxiliary pollutants, Cu, toluene, and phenanthrene were selected as the target pollutants for this study. Based on geostatistical analysis and spatial analysis, the following results were obtained. (1) The reduction in the root mean square errors (RMSEs) for Cu, toluene, and phenanthrene with multivariable cokriging was 68.4%, 81.6%, and 81.2%, respectively, which are proportional to the correlation coefficient of the relationship between the auxiliary pollutants and the target pollutants. (2) The RMSEs calculated for the multivariable cokriging were lower than those obtained by only combining one related pollutants, and two co-variables should be better. (3) The predicted results for Cu, phenanthrene, and toluene and their corresponding related pollutants are more accurate than the results obtained not using the related pollutants. (4) In the interpolation process, the RMSEs for Cu, toluene, and phenanthrene with multivariable cokriging basically increase as the neighborhood sample data increases, and then they become stable. (5) When 84, 61, and 34 sample points were removed, the RMSEs for Cu, toluene, and phenanthrene, respectively, with multivariable cokriging were close to the RMSEs of the target pollutants based on the total samples. The results are of great significance to improving the prediction accuracy of the spatial distribution of soil pollutants at coking plant sites.

Spatial methods to analyze the relationship between Spanish soil properties and cadmium content

In this study, concentrations of cadmium using 3778 samples encompassing the total size of Spain (about 505 km²) were investgated. Two novel spatial methods namely Moran eigenvector spatially varying coefficient (MESVC) and spatially filtered unconditional quantile regression (SF-UQR) were employed with the aim of avoiding the problem of local collinearity which is prevalent in regression models. Additionally, the spatially varying coefficients methods were applied to assess the influence of soil properties together with soil texture on the spatial variations of cadmium. It was indicated that the overall level of cadmium is low compared to the concentrations found around the world. In particular, the values of Cd varied between 0.01 and 2.00 mgkg^-1, with the median of 0.23 mgkg^-1. The residual standard error and adjusted R² produced by MESVC were 0.16 and 0.69, respectively which are better than 0.21 and 0.39 yielded by the SF-UQR model. Both of these models outperformed compared to the geographically weighted regression (GWR) and the performance of MESVC was also better than the traditional method of kriging. For instance, in terms of willmott index (d) and root mean squared relative error (RMSRE), the MESVC had superior performance with values equal to 0.612 and 0.275 compared to 0.399 and 0.379 obtained for the ordinary kriging. The MESVC and GWR demonstrated that CaCO₃, sand, silt and clay had a negligible influence on spatial variations of cadmium whereas, EC had the largest contribution followed by SOM and pH.

Determinants of Residents’ Willingness to Accept and Their Levels for Ecological Conservation in Ganjiang River Basin, China: An Empirical Analysis of Survey Data for 677 Households

Using the contingent valuation method and the Heckman two-stage model, we explore residents’ willingness to accept (WTA) compensation and their WTA level for ecological conservation compensation in the upstream of the Ganjiang River Basin in China. The findings reveal that 86.26% of the respondents are willing to accept compensation, and the average compensation level is ¥789.60/household per year. The residents’ gender, annual disposable income, residential location, decision on whether or not the watershed environment is important, and their satisfaction with water quality and quantity are significantly related to their WTA. The influencing factors that significantly affect compensation level are residents’ occupation, educational background, annual disposable income, family size, residential location, decision on whether or not the watershed environment is important, and their satisfaction with water quality and quantity. The results of this empirical research have important policy implications: the government should strengthen advocacy and education of watershed ecological environment protection, intensify farming and other agricultural activities, establish a differentiated and diversified compensation strategy, so as to protect and improve the ecological environment of the Ganjiang River Basin.

Prediction of Soil Organic Carbon based on Landsat 8 Monthly NDVI Data for the Jianghan Plain in Hubei Province, China

High-precision maps of soil organic carbon (SOC) are beneficial for managing soil fertility and understanding the global carbon cycle. Digital soil mapping plays an important role in efficiently obtaining the spatial distribution of SOC, which contributes to precision agriculture. However, traditional soil-forming factors (i.e., terrain or climatic factors) have weak variability in low-relief areas, such as plains, and cannot reflect the spatial variation of soil attributes. Meanwhile, vegetation cover hinders the acquisition of the direct information of farmland soil. Thus, useful environmental variables should be utilized for SOC prediction and the digital mapping of such areas. SOC has an important effect on crop growth status, and remote sensing data can record the apparent spectral characteristics of crops. The normalized difference vegetation index (NDVI) is an important index reflecting crop growth and biomass. This study used NDVI time series data rather than traditional soil-forming factors to map SOC. Honghu City, located in the middle of the Jianghan Plain, was selected as the study region, and the NDVI time series data extracted from Landsat 8 were used as the auxiliary variables. SOC maps were estimated through stepwise linear regression (SLR), partial least squares regression (PLSR), support vector machine (SVM), and artificial neural network (ANN). Ordinary kriging (OK) was used as the reference model, while root mean square error of prediction (RMSEP) and coefficient of determination of prediction (R2P) were used to evaluate the model performance. Results showed that SOC had a significant positive correlation in July and August (0.17, 0.29) and a significant negative correlation in January, April, and December (−0.23, −0.27, and −0.23) with NDVI time series data. The best model for SOC prediction was generated by ANN, with the lowest RMSEP of 3.718 and highest R2P of 0.391, followed by SVM (RMSEP = 3.753, R2P = 0.361) and PLSR (RMSEP = 4.087, R2P = 0.283). The SLR model was the worst model, with the lowest R2P of 0.281 and highest RMSEP of 3.930. ANN and SVM were better than OK (RMSEP = 3.727, R2P = 0.372), whereas PLSR and SLR were worse than OK. Moreover, the prediction results using single-data NDVI or short time series NDVI showed low accuracy. The effect of the terrain factor on SOC prediction represented unsatisfactory results. All these results indicated that the NDVI time series data can be used for SOC mapping in plain areas and that the ANN model can maximally extract additional associated information between NDVI time series data and SOC. This study presented an effective method to overcome the selection of auxiliary variables for digital soil mapping in plain areas when the soil was covered with vegetation. This finding indicated that the time series characteristics of NDVI were conducive for predicting SOC in plains.