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. This study provides new ideas for accurate soil HM spatial concentration prediction and ecological risk based on Co-kriging and GTWR models, which assign spatial and temporal heterogeneity of socio-economic activities to the process of exploring the influencing factors of HM accumulation. Our study introduces a fresh perspective for conducting large-scale pollutant studies using published data, and also provides scientific support for decision makers to develop pollution management strategies and reduce management costs.