Accounting for the Compositional Nature of Geochemical Data to Improve the Interpretation of Their Univariate and Multivariate Spatial Patterns: A Case Study from the Campania Region (Italy)

This study investigates the application of Compositional Data Analysis (CoDA) and multivariate statistical techniques to geochemical data from the soils of the Campania region. The dataset examined includes 3571 soil samples analyzed for 37 chemical elements. Principal Component Analysis (PCA) was employed to reduce the dataset’s dimensionality and identify key relationships between elements. The first PCA identified groups of highly correlated variables, which were then reduced to 20 representative elements for a second PCA. The three most significant principal components (PC1, PC2, and PC3) explained approximately 65% of the total variability. PC1 (accounting for 29.97% of variability) revealed an anticorrelation between Ti, La, and Sc with Au, Hg, and Ag, with positive scores primarily located in the inland Apennine areas. PC2 (21.8%) was dominated by Na, K, and Cu, with positive scores corresponding to volcanic deposits, aligning with the dispersion patterns of historical Vesuvian eruption products. PC3 (11%) was associated with Ca and S, with higher scores found in the alluvial plains and inland areas. These results demonstrate the efficacy of CoDA in minimizing spurious correlations and uncovering latent relationships between elements, thereby enhancing the interpretation of natural and anthropogenic processes influencing soil variability in the region.