Two-dimensional empirical mode decomposition of heavy metal spatial variation in agricultural soils, Southeast China.

The distribution of heavy metals in agricultural soils is affected by various anthropogenic activities and environmental factors occurring at different spatial scales. This paper introduced the two-dimensional empirical mode decomposition (2D-EMD) to separate the spatial variability in soil heavy metals into different scales. Geostatistics and multivariate analysis were also utilized to quantify their spatial structure and identify their potential influencing factors. The study was conducted in an arable land in southeastern China where 260 surface soil samples were collected and measured for total contents of cadmium (Cdtotal), mercury (Hgtotal), and sulfur (TS); pH; and soil organic carbon content (SOC). The results showed that both Cdtotal and Hgtotal had high coefficients of variation. The overall variation in all five soil variables was separated into three intrinsic mode functions (IMFs) and spatial residues. All three IMFs had short-range spatial correlations (1-8 km), while the spatial residues had moderate-large spatial ranges (13-39 km). IMF1 of Cdtotal was strongly correlated with IMF1 of SOC and TS, which was consistent with the principal component analysis. This indicated that IMF1 of Cdtotal represented local variations which were influenced by agricultural activities. IMFs of Hgtotal showed clustered distributions in the study area, with IMF1 and IMF2 of Hgtotal correlated in one principal component, and IMF3 of Hgtotal and IMF3 of soil pH in another component. This indicated that all three IMFs of Hgtotal might be influenced by different industrial activities or different pathways of the same industrial activities. The residues of Cdtotal and Hgtotal, representing the regional trends, only accounted for 26% of the total variance, whereas IMF1 contributed about half of the total variance. It can be concluded that agricultural activities and industrial activities were the dominant contributors of the overall variations in Cdtotal and Hgtotal in the study area, respectively.