Soils derived from black shale enriched in heavy metals pose risks to the environment and human health. Eighty sets of soil–rice samples were collected in northwest Zhejiang, China. The results showed significant enrichment of As, Cd, Pb, and Zn in the studied soil. Sequential extraction data indicated that Cd had relatively high bioavailability. Thirty-four percent of the rice samples exceeded the heavy metal pollution thresholds of national food safety standards. A comprehensive evaluation of soil heavy metals using the potential ecological risk index (RI) revealed low-to-moderate-grade risk. Risk assessment code (RAC) levels decreased as follows: Cd > Ni > Zn > Cr > Cu > Pb > As > Hg. The RAC of Cd was high and very high risk, whereas those of others were no or low risk. Cd was the major contributor to the environmental risk based on the total heavy metal concentration, toxicity, and chemical fractions. The health risk to the local populace was assessed using the method recommended by the USEPA. Local people may experience both non-carcinogenic and carcinogenic health risks, and children face greater health risks than adults. As, Cd, and Pb are the most significant contributors to non-carcinogenic health risks, and Cd is the main carcinogenic risk. In conclusion, special attention should be given to heavy metal pollution in black shale-derived soils.

Understanding the temporal–spatial distribution and influencing factors of heavy metals on a regional scale is crucial for assessing the anthropogenic impacts and natural variations in elemental geochemical behavior. This study evaluated the spatial distributions of the heavy metals As, Cd, Pb, and Zn as well as the driving mechanisms over the past 31 years in Guangxi, China, using three geochemical baseline projects (the Environmental Geochemical Monitoring Network Project (EGMON) project 1992–1996; the Geochemical Baseline (CGB) 1 project 2008–2012; and the CGB2 project 2015–2019). By calculating the variable importance using the random forest algorithm, it was found that natural factors are the primary drivers of the spatial distribution of heavy metals in the EGMON project, especially precipitation for As, the digital elevation model (DEM) for Cd and Pb, and temperature for Zn. Surface alluvial soils showed obvious heavy metal enrichment in the CGB1 project, with the gross domestic product (GDP) driving the spatial distribution of all heavy metals. In addition, the anomalous intensity and range of heavy metals in the CGB2 project decreased significantly compared with the CGB1 project, especially owing to the normalized difference vegetation index (NDVI) as a positive anthropogenic factor that improves the degree of rocky desertification, thus reducing the heavy metal contents of As and Pb, and the precipitation promoting the decomposition of Fe–Mn concretions and thus the migration of Cd and Zn. This research promotes an understanding of anthropogenic and natural influences on the spatiotemporal distribution of heavy metals and is of great significance for environmental monitoring and governance.