Source apportionment of potentially toxic elements in soils of the Yellow River Delta Nature Reserve, China: The application of three receptor models and geostatistical independent simulation

The Yellow River Delta (YRD) wetland, the most important estuary wetland in eastern China, has an important ecosystem service function. Rapid and intensive development has inevitably led to the accumulation of potentially toxic elements (PTEs) in soils. Therefore, identifying quantitative sources and spatial distributions of PTEs is essential for soil environmental protection in the YRD. A total of 240 topsoil samples (0–20 cm) were collected in the Yellow River Delta Nature Reserve (YRDNR) and analyzed the PTE contents. To avoid the biases of the single receptor model, positive matrix factorization, factor analysis with nonnegative constraints, and maximum likelihood principal component analysis-multivariate curve resolution-alternating least squares were used for source apportionment of soil PTEs. To promote the efficiency of multivariate geostatistical simulation, a minimum/maximum autocorrelation factor-sequential Gaussian simulation was built to map the spatial patterns of PTEs. Three factors were derived by the three receptor models, and their contributions to the source explanation were similar. As, Cr, Cu, Mn, Ni, and Zn originated from natural sources, with contributions of 85.6%–96.4 %. A total of 61.5 % of Hg was associated with atmospheric deposition of coal combustion and wastewater from upstream. Agricultural activities and oil exploitation contributed 33.5 % and 15.9 % of the Cd and Pb concentrations. Spatial distributions of soil PTEs were controlled by sedimentary grain size. A total of 47.2 % of the total study area was identified as hazardous area for Cd, 10.3 % for As, and 5.4 % for Hg. This work is expected to provide references for soil pollution assessment and management of YRDNR.