Hyperspectral-based Inversion of Heavy Metal Content in the Soil of Coal Mining Areas

Remote sensing inversion on heavy metal content in salinized soil of Yellow River Delta based on Random Forest Regression-a case study of Gudao Town

To explore the potential of using the mineral alteration information extracted by remote sensing technology to indirectly estimate the heavy metal content of salinized soil, 23 sampling points were uniformly set up in the town of Gudao in the Yellow River Delta as the research area in 2022. The concentrations of seven heavy metals, Cr, Cu, Pb, Zn, As, Mn and Ni, at the sampling points were determined in laboratory tests. Spectral derivative indices, topographic factors, and mineral alteration information (iron staining, hydroxyl, and carbonate ions) were extracted and screened as modeling factors using Sentinel 2 imagery. An inverse model of heavy metal content was constructed using the random forest algorithm, and the model accuracy was evaluated using the cross-validation method. The results of the study show that: (1) Hydroxyl and carbonate ion alteration can be effectively used for the inversion of soil As and Ni content in this study area. Iron-stained alteration can be used as a modeling factor in the inversion of Cr, Cu, Pb, Zn, and Mn concentrations. (2) The inclusion of alteration information improves the accuracy of heavy metal content inversion. The Cu concentration was verified to be the best predictor, with an RMSE of 3.309, MAPE of 11.072%, and R2 of 0.904, followed by As, Ni, and Zn; the predictive value of Mn, Cr and Pb was average. (3) Based on the results of concentration inversion, the high concentration areas of As, Ni, and Mn are primarily distributed on both sides of the river and around lakes and ponds. The high-concentration areas of Zn were mainly distributed in the farmland areas on both sides of the river. Areas with high concentrations of Cu were mainly distributed in the eastern oil extraction area, both sides of the rivers, and around lakes.

Estimating the Heavy Metal Contents in Entisols from a Mining Area Based on Improved Spectral Indices and Catboost

In the study of the inversion of soil multi-species heavy metal element concentrations using hyperspectral techniques, the selection of feature bands is very important. However, interactions among soil elements can lead to redundancy and instability of spectral features. In this study, heavy metal elements (Pb, Zn, Mn, and As) in entisols around a mining area in Harbin, Heilongjiang Province, China, were studied. To optimise the combination of spectral indices and their weights, radar plots of characteristic-band Pearson coefficients (RCBP) were used to screen three-band spectral index combinations of Pb, Zn, Mn, and As elements, while the Catboost algorithm was used to invert the concentrations of each element. The correlations of Fe with the four heavy metals were analysed from both concentration and characteristic band perspectives, while the effect of spectral inversion was further evaluated via spatial analysis. It was found that the regression model for the inversion of the Zn elemental concentration based on the optimised spectral index combinations had the best fit, with R2 = 0.8786 for the test set, followed by Mn (R2 = 0.8576), As (R2 = 0.7916), and Pb (R2 = 0.6022). As far as the characteristic bands are concerned, the best correlations of Fe with the Pb, Zn, Mn and As elements were 0.837, 0.711, 0.542 and 0.303, respectively. The spatial distribution and correlation of the spectral inversion concentrations of the As and Mn elements with the measured concentrations were consistent, and there were some differences in the results for Zn and Pb. Therefore, hyperspectral techniques and analysis of Fe elements have potential applications in the inversion of entisols heavy metal concentrations and can improve the quality monitoring efficiency of these soils.

Estimation of multi-media metal(loid)s around abandoned mineral processing plants using hyperspectral technology and extreme learning machine

Hyperspectral techniques are promising alternatives to traditional methods of investigating potentially toxic metal(loid) contamination. In this study, hyperspectral technology combined with partial least squares regression (PLSR) and extreme learning machine (ELM) established estimation models to predict the contents of copper (Cu), zinc (Zn), arsenic (As), cadmium (Cd), lead (Pb) and tin (Sn) in multi-media environments (mine tailings, soils and sediments) surrounding abandoned mineral processing plants in a typical tin-polymetallic mineral agglomeration in Guangxi Autonomous Region. Four spectral preprocessing methods, Savitzky-Golay (SG) smoothing, continuum removal (CR), first derivative (FD) and continuous wavelet transform (CWT), were used to eliminate noise and highlight spectral features. The optimum combinations of spectral preprocessing and machine learning algorithms were explored, then the estimation models with best accuracy were obtained. CWT and CR were excellent spectral pretreatments for the hyperspectral data regardless of the applied algorithms. The coefficients of determination (R2) of estimation models for the best accuracy of various metals (loid) are as follows: Cu (CWT-ELM:0.85), Zn (CR-PLSR:0.93), As (CWT-ELM: 0.86), Cd (CR-PLSR: 0.89), Pb (CWT-PLSR: 0.75) and Sn (CR-ELM: 0.81). In contrast, ELM models had higher accuracy with R2 > 0.80 (except Cd and Pb). In conclusion, ELM-based spectral estimation models are able to predict metal (loid) concentrations with high accuracy and efficiency, providing a potential new combinatorial approach for estimating toxic metal contamination in multi-media environments.

Mapping soil arsenic pollution at a brownfield site using satellite hyperspectral imagery and machine learning

Heavy metal contamination is ubiquitous in brownfields. Traditional site investigation employs geostatistical interpolation methods (GIMs) to predict the distribution of soil pollutants after soil sampling and chemical analysis. However, the heterogeneity of soil pollution in brownfields makes the assumptions of GIMs no longer valid and further undermines the accuracy of soil investigation. In the present study, a satellite hyperspectral image processing and machine learning method was developed to map arsenic pollution at a brownfield site. To eliminate the noise caused by atmospheric factors and increase the efficiency of spectral data, 1.3 million spectral indexes (SIs) were constructed and 1171 of them were selected due to their high correlations with soil arsenic. Five machine learning methods, i.e., Random forest (RF), ExtraTrees, Adaptive Boosting, Extreme Gradient Trees, and Gradient Descent Boosting Trees (GDB) were built to predict soil arsenic. The RF method was found to render the best performance (r = 0.78), reducing 30 % of prediction errors compared with traditional GIMs. RF also maintained a relatively higher level of accuracy (r = 0.56) when the sampling grids increased to 100 m, which was higher than that of GIMs under a 50 m sampling grid (r = 0.42), revealing that the proposed method can provide more accurate results with fewer sampling points, namely less investigation cost. It was indicated that the second derivate was the most efficient preprocessing method to remove spectral noise and normalized difference (ND) was the most reliable spectral index construction strategy. Based on uncertainty analysis, the heterogeneity of soil arsenic distribution was considered the most influential factor causing prediction errors. This study demonstrates that machine learning based on satellite visible and near-infrared reflectance spectroscopy (VNIR) is a promising approach to map soil arsenic contamination at brownfield sites with high accuracy and low cost.

Field hyperspectral data and OLI8 multispectral imagery for heavy metal content prediction and mapping around an abandoned Pb-Zn mining site in northern Tunisia

Mining and smelting releases toxic contaminants such as zinc (Zn), lead (Pb) or cadmium (Cd) into the soil thereby poisoning it and rendering it unproductive. Remotely alternatives have been widely employed in the attempt of estimating heavy metal content within soils. The present study provides a methodological approach based on VNIR-SWIR field hyperspectral data and multispectral Landsat OLI 8 imageries for the prediction and mapping of Pb, Zn and Cd heavy metal contents around the abandoned Jebel Ressas mine site in Northern Tunisia. Thus, eighty-seven soil and tailing samples were collected from the study site and VNIR-SWIR field reflectances were measured on the same collection points, as well. All samples were analysed by atomic absorption for the estimation of heavy metal concentrations. The partial least squares regression PLSR was conducted considering the measured heavy metal concentrations and using multi-scale data: VNIR-SWIR field hyperspectral data and multispectral Landsat OLI 8 imagery. Standard normal variable (SNV) and multiple scatter correction (MSC) preprocessing methods were applied for further mapping improvement. Thus, this work aims to automate the estimation of the heavy metal contents in contaminated soils, by carrying out: a modeling approach based on the PLSR using VNIR-SWIR field hyperspectral data, ii) the mapping of Pb and Zn contents thanks to the exploitation of Landsat OLI8 multispectral imagery and iii) the application of both MSC and SNV preprocessing methods to optimize the performance of the developed models, when using such spectrally and spatially degraded data.

Retrieving soil heavy metals concentrations based on GaoFen-5 hyperspectral satellite image at an opencast coal mine, Inner Mongolia, China

Soil heavy metals pollution has been becoming one of the severely environmental issues globally. Previous studies reported laboratory-measured spectra could be used to infer soil heavy metals concentrations to some extent. However, using field-obtained spectra to estimate soil heavy metals concentrations is still a great challenge due to the low precision and weak efficiency at large scales. The present study collected 110 topsoil samples from an Opencast Coal Mine of Ordos, Inner Mongolia, China. Then, the spectra and soil heavy metals concentrations of samples were measured under laboratory conditions. The direct standardization (DS) algorithm was introduced to calibrate the Gaofen-5 (GF-5) hyperspectral image based on the measured spectra of samples. The spectral reflectance of the GF-5 hyperspectral image was reconstructed using continuous wavelet transform (CWT) at different scales. The characteristic bands of GF-5 for estimating heavy metals concentrations were selected by the Boruta algorithm. Finally, the random forest (RF), the extreme learning machine (ELM), the support vector machine (SVM), and the back-propagation neural network (BPNN) algorithms were used to predict the heavy metals concentrations. Some findings were achieved. First, CWT can effectively eliminate the noise of satellite hyperspectral data. The characteristic bands of Zn (480-677, 827-1029, 1241-1334, 1435-1797, and 1949-2500 nm), Ni (514-630, 835-985, 1258-1325, 1460-1578, and 1949-2319 nm), and Cu (822-831; 1029-1300, 1486-1595, and 1730-2294 nm) can be effectively retrieved via the Boruta algorithm. Second, the estimation accuracy was significantly improved by using the DS algorithm. For zinc (Zn), nickel (Ni), and copper (Cu), the determination coefficients of the validation dataset (R_v²) were 0.77 (RF), 0.62 (RF), and 0.56 (ELM), respectively. Third, the distribution trends of heavy metals were almost consistent with the results of actual ground measurements. This paper revealed that the GF-5 can be one of the reliable satellite hyperspectral imagery for mapping soil heavy metals.

Diagnosis of cadmium contamination in urban and suburban soils using visible-to-near-infrared spectroscopy

Previous studies have mostly focused on using visible-to-near-infrared spectral technique to quantitatively estimate soil cadmium (Cd) content, whereas little attention has been paid to identifying soil Cd contamination from a perspective of spectral classification. Here, we developed a framework to compare the potential of two spectral transformations (i.e., raw reflectance and continuum removal [CR]), three optimization strategies (i.e., full-spectrum, Boruta feature selection, and synthetic minority over-sampling technique [SMOTE]), and three classification algorithms (i.e., partial least squares discriminant analysis, random forest [RF], and support vector machine) for diagnosing soil Cd contamination. A total of 536 soil samples were collected from urban and suburban areas located in Wuhan City, China. Specifically, Boruta and SMOTE strategies were aimed at selecting the most informative predictors and obtaining balanced training datasets, respectively. Results indicated that soils contaminated by Cd induced decrease in spectral reflectance magnitude. Classification models developed after Boruta and SMOTE strategies out-performed to those from full-spectrum. A diagnose model combining CR preprocessing, SMOTE strategy, and RF algorithm achieved the highest validation accuracy for soil Cd (Kappa = 0.74). This study provides a theoretical reference for rapid identification of and monitoring of soil Cd contamination in urban and suburban areas.

Estimation of Heavy Metals in Tailings and Soils Using Hyperspectral Technology: A Case Study in a Tin-Polymetallic Mining Area

Rapid assessment of heavy metal (HM) pollution in mining areas is urgently required for further remediation. Here, hyperspectral technology was used to predict HM contents of multi-media environments (tailings, surrounding soils and agricultural soils) in a mining area. The correlation between hyperspectral data and HMs was explored, then the prediction models were established by partial least squares regression (PLSR) and back propagation neural networks (BPNN). The determination coefficients (R²), root mean squared error and ratios of performance to interquartile range (RPIQ) were used to evaluate the performance of the models. Results show that: (1) both PLSR and BPNN had good prediction ability, and (2) BPNN had better generalization ability (Cu (R² = 0.89, RPIQ = 3.05), Sn (R² = 0.86, RPIQ = 4.91), Zn (R² = 0.74, RPIQ = 1.44) and Pb (R² = 0.70, RPIQ = 2.10)). In summary, this study indicates that hyperspectral technology has potential application in HM estimation and soil pollution investigation in polymetallic mining areas.

The Spatial Distribution and Prediction of Soil Heavy Metals Based on Measured Samples and Multi-Spectral Images in Tai Lake of China

Soil is an important natural resource. The excessive amount of heavy metals in soil can harm and threaten human health. Therefore, monitoring of soil heavy metal content is urgent. Monitoring soil heavy metals by traditional methods requires many human and material resources. Remote sensing has shown advantages in the field of monitoring heavy metals. Based on 971 heavy metal samples and Sentinel-2 multi-spectral images in Tai Lake, China, we analyzed the correlation between six heavy metals (Cd, Hg, As, Pb, Cu, Zn) and spectral factors, and selected As and Hg as the input factors of inversion model. The correlation coefficient of the best model of As was 0.53 (p < 0.01), and of Hg was 0.318 (p < 0.01). We used the methods of partial least squares regression (PLSR) and back propagation neural network (BPNN) to establish inversion models with different combinations of spectral factors by using 649 measured samples. In addition, 322 measured samples were used for accuracy evaluation. Compared with the PLSR model, the BP neural network builds the model with higher accuracy, and B1-B4 combined with LnB1-LnB4 builds the model with the highest accuracy. The accuracy of the best model was verified, with an average error of 19% for As and 45% for Hg. Analyzing the spatial distribution of heavy metals by using the interpolation method of Kriging and IDW. The overall distribution trend of the two interpolations is similar. The concentration of As elements tends to increase from north to south, and the relatively high value of Hg elements is distributed in the east and west of the study area. The factories in the study area are distributed along rivers and lakes, which is consistent with the spatial distribution of heavy metal enrichment areas. The relatively high-value areas of heavy metal elements are related to the distribution of metal products factories, refractory porcelain factories, tile factories, factories and mining enterprises, etc., indicating that factory pollution is the main reason for the enrichment of heavy metals.

Retrieving zinc concentrations in topsoil with reflectance spectroscopy at Opencast Coal Mine sites

Heavy metals contaminations in mining areas aroused wide concerns globally. Efficient evaluation of its pollution status is a basis for further soil reclamation. Visible and near-infrared reflectance (Vis-NIR) spectroscopy has been diffusely used for retrieving heavy metals concentrations. However, the reliability and feasibility of calibrated models were still doubtful. The present study estimated zinc (Zn) concentrations via the random forest (RF) and partial least squares regression (PLSR) using ground in-situ Zn concentrations as well as soil spectral reflectance at an Opencast Coal Mine of Ordos, China in February 2020. The coefficient of determination (R2), root mean square error (RMSE), mean absolute error (MAE), and the ratio of performance to deviation (RPD) were selected to assess the robustness of the methods in estimating Zn contents. Moreover, the characteristic bands were chosen by Pearson correlation analysis and Boruta Algorithm. Finally, the comparison between RF and PLSR combined with eight spectral reflectance transformation methods was conducted for four concentration groups to determine the optimal model. The results indicated that: (1) Zn contents represented a skewed distribution (coefficient of variation (CV) = 33%); (2) the spectral reflectance tended to decrease with the increase of Zn contents during 580-1850 nm based on Savitzky-Golay smoothing (SG); (3) the continuous wavelet transform (CWT) demonstrated higher effectiveness than other spectral reflectance transformation methods in enhancing spectral responses, the R2 between Zn contents and the soil spectral reflectance achieved the highest (R2 = 0.71) by using CWT; (4) the RF combined with CWT exhibited the best performance than other methods in the current study (R2 = 0.97, RPD = 3.39, RMSE = 1.05 mg kg-1, MAE = 0.79 mg kg-1). The current study supplied a scientific scheme and theoretical support for predicting heavy metals concentrations via the Vis-NIR spectral method in possible contaminated areas such as coal mines and metallic mineral deposit areas.

VIRS based detection in combination with machine learning for mapping soil pollution

Widespread soil contamination threatens living standards and weakens global efforts towards the Sustainable Development Goals (SDGs). Detailed soil mapping is needed to guide effective countermeasures and sustainable remediation operations. Here, we review visible and infrared reflectance spectroscopy (VIRS) based detection methods in combination with machine learning. To date, proximal, airborne and spaceborne carrier devices have been employed for soil contamination detection, allowing large areas to be covered at low cost and with minimal secondary environmental impact. In this way, soil contaminants can be monitored remotely, either directly or through correlation with soil components (e.g. Fe-oxides, soil organic matter, clay minerals). Observed vegetation reflectance spectra has also been proven an effective indicator for mapping soil pollution. Calibration models based on machine learning are used to interpret spectral data and predict soil contamination levels. The algorithms used for this include partial least squares regression, neural networks, and random forest. The processes underlying each of these approaches are outlined in this review. Finally, current challenges and future research directions are explored and discussed.