Next Article in Journal
Effects of Insecticide and Herbicides on Thyroid Disturbances in Zebrafish
Next Article in Special Issue
Comparative Remediation of Arsenic and Antimony Co-Contaminated Soil by Iron- and Manganese-Modified Activated Carbon and Biochar
Previous Article in Journal
Plastic Pollution in Paradise: Analyzing Plastic Litter on Malta’s Beaches and Assessing the Release of Potentially Toxic Elements
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

GIS- and Multivariate-Based Approaches for Assessing Potential Environmental Hazards in Some Areas of Southwestern Saudi Arabia

by
Hassan Alzahrani
1,
Abdelbaset S. El-Sorogy
1,
Abdurraouf Okok
2 and
Mohamed S. Shokr
3,*
1
Geology and Geophysics Department, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
2
Earth Sciences and Engineering Department, Missouri University of Science and Technology, McNutt Hall, 1400 N. Bishop Ave, Rolla, MO 65401, USA
3
Soil and Water Department, Faculty of Agriculture, Tanta University, Tanta 31527, Egypt
*
Author to whom correspondence should be addressed.
Toxics 2024, 12(8), 569; https://doi.org/10.3390/toxics12080569
Submission received: 20 June 2024 / Revised: 30 July 2024 / Accepted: 30 July 2024 / Published: 3 August 2024

Abstract

:
Soil contamination is a major issue that endangers the ecology in most countries. Total concentrations of As, Cd, Co, Cr, Cu, Mn, Ni, Pb, VFe, and Zn were determined by analyzing soil samples from 32 surface soil samples in southwest Saudi Arabia, including certain areas of Al-Baha. Kriging techniques were used to create maps of the distribution of metal. To assess the levels of soil contamination in the research area, principal component analysis (PCA), contamination factors (CF), and pollution load index were used. The results show the stable model gave the best fit to the As and Zn semivariograms. The circular model fits the Cd, Co, and Ni semivariograms the best, while the exponential model fits the Cr, V, and Fe semivariograms the best. For Ni and Pb, respectively, spherical and Gaussian models are fitted. The findings demonstrated two clusters containing different soil heavy metal concentrations. According to the data, there were two different pollution levels in the research region: 36.58% of it is strongly contaminated, while 63.41% of it has a moderate level of contamination (with average levels of these metals 5.28 ± 5.83, 0.81 ± 0.19, 18.65 ± 6.22, 45.15 ± 23.25, 60.55 ± 23.74, 972.30 ± 223.50, 33.45 ± 14.11, 10.05 ± 5.13, 84.15 ± 30.72, 97.40 ± 30.05, and 43,245.00 ± 8942.95 mg kg−1 for As, Cd, Co, Cr, Cu, Mn, Ni, Pb, V, Fe, and Zn, respectively). The research area’s poor management practices are reflected in the current results, which raised the concentration of harmful elements in the soil’s surface layers. Ultimately, the outcomes of pollution concentration and spatial distribution maps could aid in informing decision-makers when creating suitable heavy metal mitigation strategies.

1. Introduction

Soil pollution is a major challenge that worries decision-makers worldwide [1]. There are five million locations worldwide where heavy metal or metalloid pollution of the soil exists at concentrations higher than allowed [2]. Due to their lengthy half-lives in the environment and inability to be chemically or biologically broken down once deposited into the soil, heavy metals have the potential to contaminate the environment and harm ecosystems [3,4] seriously. Providing the right nutrients for plants in the soil is largely dependent on the quality of the soil, which is why sustainable soil management is associated with sustainable agriculture [5]. Since soil is essential for producing high-quality, healthful meals, it should be improved, protected from contamination, and its fertility level raised to meet the world’s growing agricultural needs sustainably [5]. Heavy metal poisoning in soil is a global problem due to the threats to the ecosystem and biological toxicity [6,7,8]. Because of their persistence, ease of accumulation in sediments, and bioaccumulation in food chains, as well as the fact that they degrade the ecological environment and endanger human health, they are regarded as significant environmental pollutants [9,10,11]. In recent decades, there has been an increasing risk to human health and food security due to the negative impacts of pollutants on crop quality [12]. The Al-Baha region is situated in the southwest of Saudi Arabia, between the Tihama Plain and the Hijaz Mountains, with an area of 10,362 km2 [13]. Metals are elements found everywhere in the environment, whether through anthropogenic or natural processes [14]. Metals like Pb, Cr, Cd, Hg, Cu, Ni, and Zn are among those on the priority list, according to the Agency of Toxic Substances and Disease Registry [15]. Metalloids such as As have the propensity to create covalent connections with organic groups, which gives rise to their toxicological characteristics [16]. Heavy metal interpolation has been done using geostatistical methods (kriging/cokriging) in the UK [17] and Ireland [18] and empirical Bayesian kriging in California, the USA [19], and China [20]. Geostatistics is a strong technique for representing the spatial variation of soil properties [18]. Through spatial interpolation, kriging is a potent geostatistical technique that may be used to integrate data into raster maps and analyze the spatial distribution of soil properties [17,21,22]. In essence, it uses variograms and associated factors to determine the spatial structure of soil variables (nugget, sill, and range) [17,21]. Both single and combined pollution indices are seen to be useful quantitative methods for determining the concentration of heavy metals in soils and are crucial for forecasting environmental sustainability in the future, especially when it comes to agriculture [23]. Principal component analysis (PCA) has been used to identify several causes of soil pollution, including industrial and agricultural activities and the percentage of heavy metals causing soil contamination [8]. Another advantage of PCA is its capacity to handle massive volumes of data without being constrained by quantity [24,25]. Another aspect of aggregative hierarchical clustering (AHC) is the distances between samples where most similar points are concentrated in a single cluster. The process of continuously merging the two closest clusters is known as AHC, an unsupervised classification method. The current work intends to define contamination levels using PCA, assess the pollution load index (PLI) of the study region, and map the geographical distribution of several selected heavy metals in some areas of southwest Saudi Arabia (Al-Baha). The goals are to look into soil overall concentrations of As, Cd, Co, Cr, Cu, Mn, Ni, Pb, V, Fe, and Zn in certain regions of Saudi Arabia’s southwest. The statistical analysis used in this study enables the evaluation of cause-and-effect relationships and reveals surpassed levels, making it a valuable tool for locating possible sources of contamination. Our study is one of the few that evaluates heavy metal contamination in the research area, as far as we are aware. The findings offer useful details on the extent of soil contamination in the study area, which policymakers can utilize as a reference for both quantitative and qualitative management purposes.

2. Materials and Methods

2.1. Site Description

The investigated area is located in the northwest of the Al-Baha region between latitudes 19°57′0″ and 20°36′0″ N and longitudes 41°9′0″ and 41°30′0″ E (Figure 1), with an area of 1435.66 km2. The Al-Baha region is distinguished by a variety of soil types, plant covers, and topographical and climatic trends. The majority of the crystalline basement of the Al-Baha region, which is situated on the Arabian Shield, is made up of Precambrian continental crust [26]. With an average temperature between 12 and 23 °C, the summer is typically mild to hot, while the winter is warm to chilly. About 150–200 mm of precipitation falls on average each year, and between 50 and 70% of it is moist [27]. The research region is divided geologically into two major Precambrian assemblages: (1) the Jiddah and Ablah groups of metamorphosed coarse clastic and meta-andesitic assemblage, and (2) the Baish and Baha groups of metamorphosed basalt, graywacke, and chert [28], and shanti [29]. The main crops in the study area are wheat, barley, and grape trees (Figure S1).

2.2. Collecting Samples and Analytical Procedures

Thirty-two soil samples were collected (Figure 1) using a plastic hand trowel; the samples were taken at a depth of less than 30 cm from the soil’s surface. To provide a representative sample, a composite sample of three subsamples was properly mixed. They were then put in plastic sample bags and kept in an ice box. Inductively coupled plasma-atomic emission spectrometry (ICP–AES) was used in ALS Arabia’s ISO/IEC 17025 (https://www.iso.org/standard/66912.html, accessed on 30 July 2024) (2017)-compliant and accredited to analyze As, Cd, Co, Cr, Cu, Mn, Ni, Pb, V, Fe, and Zn, as per USEPA 3050B [30]. Before chemical analysis, soil samples were sieved and allowed to air dry. A dry, clean Teflon beaker containing around 200 mg of samples was filled with 2 mL of HNO3, 6 mL of HCl, and 2 mL of HF after the samples were precisely weighed [31]. The samples underwent a 40-minute digestion process on a heated plate with sand at a mild temperature of 60–120 °C [13]. After filtering, the digest was then put into disposable 25 mL tubes. An empty digest was carried out identically. External calibration was used to do the ICP-AES calibration. Standard stock solutions were utilized, which included a multi-element calibration standard of 1000 mg L1. HNO3 was used to dilute both the stock and standard solutions. The ICP–AES technique was validated in terms of linearity, limits of quantification (LOQs), and limits of detection (LODs). The correlation coefficient (R2) between concentrations and the detector signal for each metal was calculated using linear regression analysis and ranged from 0.95 to 98. The LOQs LODs of each element are illustrated in Table S1.

2.3. Contamination Indices

The heavy metal (HM) contamination in soil samples was estimated using the contamination factor (CF) and pollution load index (PLI) [23,32]. According to Hökanson [33], the contamination factor (CF) can be used to express the degree of contamination. The ratio known as the CF is calculated by dividing each metal’s concentration in the sediment by the background or baseline value (Equation (1)). Table S2 shows the CF categorization that is used to assess the levels of soil pollution. Based on element abundances in sedimentary rocks (shale), the background value relates to the baseline concentrations reported by [34]. Contamination factor (CF) values are multiplied by n numbers. The integrated pollution status of the related hazardous groups at the sampling sites was ascertained using the PLI (Equation (2)). The PLI that Tomlinson [35] proposed gives the local population some knowledge about the amount of a component in the environment. The PLI of a single site is equal to the nth root of the n multiplied values of the contamination factor (CF). Table S3 displays the PLI classification used to evaluate soil pollution levels [36].
C F = C m e t a l C b a c k g r o u n d
where:
  • C m e t a l = total measured concentration of heavy metal.
  • C b a c k g r o u n d is each metal’s background value.
P L I = ( C F 1 C F 2 C F 3 .   C F n ) 1 / n
where: CF stands for contamination factor and n is the number of specific heavy metals.

2.4. Multivariate Analysis

SPSS 28.0 (IBM, Armonk, NY, USA) and Microsoft Excel 2021 (Sacramento, CA, USA) were used for the statistical analysis of the data. One-way ANOVA was utilized for the analysis of variance (ANOVA) [37]. Z-score calculations were performed to normalize the dataset. Managing environmental data, which typically exhibit aberrant distribution, requires this crucial step [22]. To investigate metal correlations in soils, Pearson’s correlation was then run on the produced z-scores. The principal component (PC) approach, along with Bartlett’s test of sphericity, and the Kaiser–Meyer–Olkin (KMO) method, the samples’ appropriateness for PCA was determined. In the event where KMO values were above 0.5, the data were suitable for PCA [38] (Table S4). A PC was only taken into consideration if its eigenvalue was greater than 1.0 [39].

2.5. Geostatistical Models

Using ordinary kriging (OK) with circular, spherical, exponential, and Gaussian models, the GIS geostatistical analyzer was utilized to produce raster maps for soil HM. To generate continuous surface layers, the OK is a sophisticated algorithm that forecasts a property’s value at an unsampled point [40]. Using Equation (3), one may estimate the predicted value Z x 0 based on measured data Z( x i ), weights of measured values ( λ i ) at a particular distance, and the number of predicted values (n) within specific neighbor samples.
Z x 0 = i = 1 n λ i × Z x i
The following is the definition of the exponential function:
Ɣ h = 0 , h = 0 C 0 + C 1 e h a ,   h > 0  
It was defined that the Gaussian function is:
Ɣ h = C 0 + C 1 exp h 2 a 2 ,   h > 0 0 , h = 0
The spherical function was defined as:
Ɣ h = C 0 + C 3 h 2 a 1 2 h a 3 ,   0 < h a C 0 + C ,    h > a 0 ,   h = 0
The definition of the circular function was given as:
Ɣ h = C 0 + C ( 1 2 π cos 1 h 2 + 1 ( h 2 / a 2 ) ) ,   0 < h a C 0 + C ,   h > a 0 ,   h = 0
The real ranges for the spherical, circular, exponential, and Gaussian functions, respectively, are represented by the letters and in these formulas. The nugget is denoted by C0, the partial sill by C, and the spatial lag by h. For these variograms, the soil sample spatial variation was isotropic. OK can provide fair approximations with the least amount of error [41]. The nugget (C0), partial sill (C), and sill (C0 + C) are the spatial variation parameters that are obtained by fitting the experimental semivariograms using multiple models. The nugget, which measures short-range variability, is the semivariogram at a lag distance of zero. The sill represents the overall sample variability and is the point at which the model flattens out [42].

2.6. Validation of Geostatistical Analysis

Heavy metal concentrations were mapped using the geostatistical models previously discussed. The mean standardized error (MSE) and the root mean square standardized error (RMSSE) were the prediction errors that were taken into account when using the cross-validation technique to test the validity and efficiency of OK models. The best-fit model was chosen based on the lowest MSE (close to zero) and RMSSE close to unity [43,44]. The mathematical expressions of these errors are as follows [45]. The models were assessed using the following two equations (Equations (8) and (9)).
Mean   standardized   error   ( MSE ) = 1 N i = 1 N Z 1 X 1 Z 2 X 2
Root   mean   square   standardized   error   ( RMSSE ) = 1 N I = 1 N Z 1 x i Z 2 x i 2
where: Z 1 ( x i ) is measured values and Z2 ( x i ) = predicted values.

3. Results and Discussion

3.1. Variation of Heavy Metals within the Study Area

Table 1 summarizes the concentrations of heavy metals in the topsoil along the study area. The As has been determined to be between 0.1 and 22 mg kg−1. The mean As concentration (7.1 ± 6.3 mg kg−1) surpasses the upper crust concentrations by Wedepohl [34], and the recommended levels of heavy metals by DEA [46] (Table 1). The element arsenic is more common in clayey soil types and is generally believed to have a geological origin—earth. However, since man-made sources release arsenic into the environment more frequently than natural ones, there is a considerable quantity of pollution from anthropogenic sources in the environment [47]. The average value of the Cd concentration was 0.8 ± 0.2 mg kg−1, with a range of 0.1 to 1.1 mg kg−1. Cadmium (Cd) is a trace element that is commonly present in the environment but is not necessary. Anthropogenic and geogenic sources can increase the levels of Cd in soils and groundwater, which is crucial for preserving a wholesome food supply and clean drinking water. High Cd levels can cause cancer in humans [48]. The cobalt concentration ranged from 9 to 35 mg kg1, with an average of 20.6 ± 6.6 mg kg−1. When cobalt is present at low concentrations, it is essential to the growth of leguminous plants. Due to its inclusion in Vitamin B12, it also offers numerous advantages for human health. However, when used excessively, it can have detrimental effects, such as restricting the nitrogen metabolism and photosynthesis of plants, and it can also have detrimental effects on the human heart and lungs [49]. The average concentration of chromium was 58.2 ± 34.4 mg kg−1, ranging from 13 to 1914 mg kg−1. One of the harmful heavy metals present in soil naturally is chromium, which is produced by the earth’s crust’s weathering of minerals or by industrial leftovers that seep into the soil. Because it is crucial for the growth of leguminous plants when it is present in low concentrations, cobalt is a crucial component of plant nutrition. Plant limitations in photosynthesis and nitrogen metabolism are just two of the several issues that arise from an increase in cobalt concentration in the soil [50]. One of the primary sources of cobalt contamination is fertilization [51]. The copper mg kg1 average is 61.1 ± 20.5. The study area’s soil samples exhibit greater amounts of Cu than both the typical upper earth crust values of Wedepohl [34] and the DEA [46] (Table 1). The primary cause of copper buildup in soils is human activity, such as mining and industrial processes. In agriculture, copper-containing chemicals are commonly utilized, especially in pesticides applied to vineyards and orchards [52]. The concentration of Mn in total varied between 684 and 1565.0 mg kg−1, with an average of 1015.2 ± 212 mg kg−1. The human body, plants, and animals all rely heavily on manganese (Mn), which is an important cofactor of many enzymes. As an alloying element, it is also an essential raw ingredient. The environment is overexposed to it when it is used extensively for industrial purposes, endangering both public health and the ecosystem [53]. The overall concentration of Ni ranged from 684 to 1565 mg kg−1, with a mean of 1015.2 ± 212 mg kg−1. The proper development of plants, animals, and soil/water bacteria are just a few of the crucial biological processes that are known to depend on nickel. However, too much nickel can poison wildlife. Researchers have discovered that nickel has an impact on higher plants’ ability to photosynthesize, significantly reducing soil fertility and contributing to a number of chronic human ailments [54]. The range of Pb content was 4–27 mg kg−1, with an average of 10.1 ± 5.5 mg kg−1. Lead pollution in the soil may be caused by pesticides and fertilizers. Lead deposition in agricultural fields can also result from the usage of lead ammunition for hunting [55]. Vanadium (V) is contained in the entire composition at an average of 89.8 ± 33 mg kg−1. According to [56,57], vanadium is a strategically significant metal that is widely used in modern civilization in the manufacturing of steel alloys and sulfuric acid. Vanadium can be found in nature in two different oxidation states: tetravalent and pentavalent [58,59]. Because pentavalent vanadium negatively affects phosphate metabolism, it is more harmful to humans, animals, and plants than tetravalent vanadium [60]. Typical soil iron concentrations range from 0.2% to 55% (20,000 to 550,000 mg/kg−1) [61]. Because of different soil types and the existence of additional sources, concentrations can vary greatly, even within localized areas. With an average of 103.08 ± 17.63 mg kg−1, the Zn values in the research region varied greatly, from 13 to 77 mg kg1. Both humans and plants need zinc, but too much of it can be dangerous [62]. It might inflict immediate damage, which could result in problems with the immune and digestive systems. Copper insufficiency symptoms can also be caused by excessive zinc levels that prevent copper absorption [63].

3.2. Heavy Metals Distribution within the Investigated Area

Table 2 provides the parameters for the semivariograms (Figure S2) that illustrate the regional distribution of the eleven metals. The stable model was the best model for the As and Zn semivariograms, with the lowest errors. The best model to fit the semivariograms of Cd, Co, and Ni was the circular model; the best model to fit the semivariograms of Cr, V, and Fe was the exponential model. Spherical and Gaussian models are fitted to Ni and Pb, respectively. There are strong correlations between the eleven metals’ measured and anticipated concentrations, according to the semivariogram models’ cross-validation (Figure S2). The prediction errors (Table 2) show that while the RMSSE values were near unity, the MSE values for each applied model were nearly zero. Table 2 illustrates that for every used OK model, an appositive nugget effect (larger than 0) and a sill value were attained. The spatial dependency (SPD) is defined by the nugget/sill ratio; values less than 0.25, 0.25–0.75, and greater than 0.75 indicate, correspondingly, a strong, moderate, and weak SPD [45]. As a result, Fe has a strong SPD, whereas the other elements have a moderate SPD. Similar spatial patterns were seen for the concentrations of Cd, Cr, and Ni, which were delineated by a northwest zone with high values. The eastern parts of the study area showed the highest levels of Cu, Cd, and Zn, while the southern part showed the highest concentrations of Mn and V, and the southern part showed the highest concentrations of Pb. The northern and middle regions of the study area showed the highest concentrations of As and Fe, respectively, Figure 2.

3.3. Correlation between Selected Heavy Metals

A linear relationship can be estimated using Pearson’s correlation coefficient. In all of the natural sciences, it is one of the most often utilized statistical values [64]. There were substantial correlations (p < 0.05) and highly significant correlations (p < 0.01) between different metals (Table 3). The concentration of As is positively correlated with Ni (r = 0.483 **) and negatively correlated with Cu (r = −0.425 *), V (r = −0.367 *), and Fe and Zn (r = −0.360 *). Cd did not significantly correlate with the other metals. Co was correlated with Cu (r= 0.642 **), V (r = 0.824 **), and Ni (r = 0.399 *). The concentration of Mn is positively correlated with Co (r = 0.380) due to their comparable chemical characteristics, cobalt (Co) and manganese (Mn) are closely related in soils [65]. There is a substantial positive association between the chemically related elements Co and Fe (0.920) in the same line of [66]. Cr is correlated positively with Co (r = 0.576 **), V (r = 0.469 **), Cu (r = 0.474 **), Ni (r = 0.716 **), and Fe (r = 0.466 **) [67]. Pb did not significantly correlate with the other metals.

3.4. Analysis of Factors Extracted from PCA

Table 4 and Table 5 display the PCA results for the metal content in the soil samples. The first eigenvalue results in Table 4 and Figure 3 indicate that factor analysis made it possible to extract four components (F1, F2, F3, and F4). CO (0.972), V (0.874), and Fe (0.945) were shown to be substantially linked with the first component (F1) in the component matrix for the data. Additionally, soil Cr and Cu (0.630 and 0.665), with moderate absolute concentrations, were linked to F1. With moderate values, As (−0.697), Ni (−0.766), and Zn (0.584) were clearly visible in the second component (F2). The levels of Cd (0.657) and Mn (0.703) had a moderate relationship with the third component (F3). (Table 5). There was a moderate correlation between Pb (0.591) and the fourth PC component (F4). F1, which accounted for 35.77% of the total variance, was dominated by CO, Cr, Cu, V, and Fe, indicating that these metals came from comparable sources (Table 4 and Figure 4) [68]. F2 mostly ascribed the variance to As, Ni, and Zn and explained 20.59% of the overall variation. As and Ni had a moderately negative association with F2 while Zn had a moderately positive association with F2, suggesting that their sources might have been distinct [68]. Cd and Mn dominated F3, which accounted for 12.77 percent of the variance in total. Pb dominated F4 and explained 10.35% of the overall variation. The principal component analysis biplot of F1, F2, F3, and F4 is displayed in Figure 4. It rescales the loadings plot and the score plot so that they overlap on a single plot. The biplot displays the variables as arrows, and the correlation between each pair of variables is determined by taking the cosine of the angle produced by the arrows. The narrower the angle between each pair of arrows, the higher the correlation between the variables [69].

3.5. Cluster Analysis of Study Area

A dendrogram representing the findings of the hierarchical cluster analysis (HCA) for the variables is displayed in Figure 5. The first cluster has 13 observations, whereas the second has 19, according to descriptive statistics presented in Table 6. The ranges, means, and standard deviations (SD) of each variable differ. Figure 5 displays the locations of every cluster observation. Using PCA, those two clusters were selected out of factors (F1, F2, F3 and F4). Cluster 1 exhibits greater amounts of Co, Cr, Cu, Mn, V, Fe, and Zn, while cluster 2 displays lower quantities of As and Cd (Table 6).

3.6. Hazard Assessment

3.6.1. Contamination Factor (CF)

For every cluster, the contamination factors (CFs) were computed to evaluate the likely ecological risk associated with a particular heavy metal. The cluster classification of each site is displayed in Figure 6. According to the results, cluster 1’s contamination factor of As indicated low contamination (38.46%), moderate contamination (38.46%), considerable contamination (15.38%), and very high contamination (7.69%); in contrast, cluster 2’s contamination level was relatively high, with 5.26% of soil samples having low contamination, 26.31% having moderate contamination, 42.10% having considerable contamination, and 26.1% having very high contamination [33], 1980) (Figure 6). In terms of Cd, CF was very high in both clusters. It is possible that human activity and fertilizer application in the research region are to blame for the rise in Cd levels in both clusters [13]. Co pollution was found in cluster 1 at moderate and considerable levels, whereas low, moderate, and considerable contamination levels were found in cluster 2. In comparison, cluster 1 recorded 7.69% low, 84.61% moderate, and 7.69% considerable contamination, whereas, the contamination levels in cluster 2 were 26.31, 68.42, and 5.26%, respectively. A decline in soil quality could have an impact on human health if the amount of chromium in the soil increases [70]. Cluster 1 has a considerable and very high level of contamination, measuring 61.53% and 38.47%, respectively, according to the Cu contamination factor. Additionally, cluster 2 had a relative percentage of 26.31 and 73.38%, respectively, indicating moderate and considerable Cu contamination. Mn pollution was found at moderate levels in both clusters. Cluster 1 has a moderate level of contamination, as indicated by the Ni contamination factor of 53.84%. In addition, cluster 2 had a relative proportion of 73.68. According to the findings, Pb contamination in cluster 1 was found to be low (100%). Cluster 2 showed low (84.21%) to moderate levels of contamination (15.78%). While cluster 1 recorded 92.3% moderate contamination and 7.69% considerable contamination by V, cluster 2 recorded 5.26% moderate contamination and 94.73% considerable contamination. Zn and Fe had moderate levels of contamination in cluster 1, whereas moderate, moderate, and considerable amounts of contamination were found in cluster 2, respectively (Figure 6).

3.6.2. Pollution Load Index (PLI)

There were two levels of pollution in the research region, as shown in Table 7 and Figure 7. Some 36.58% of the research area is significantly contaminated, whereas 63.41% of the region indicated a moderate level of pollution. Moderate concentrations of heavy metals were found in the majority of the study area, with average levels of these metals 5.28 ± 5.83, 0.81 ± 0.19, 18.65 ± 6.22, 45.15 ± 23.25, 60.55 ± 23.74, 972.30 ± 223.50, 33.45 ± 14.11, 10.05 ± 5.13, 84.15 ± 30.72, 97.40 ± 30.05, and 43,245.00 ± 8942.95 mg kg−1 for As, Cd, Co, Cr, Cu, Mn, Ni, Pb, V, Fe, and Zn, respectively. Anthropogenic sources for Cu, Zn, Cd, and Pb may have originated from sewage and agricultural practices; natural sources for As, Cr, Mn, Fe, Co, and Ni, mainly from the weathering of soil minerals and atmospheric deposition [13]. Farmers should use manure and biofertilizers, which have different levels of heavy metals and patterns, to reduce their dependency on chemical pesticides and fertilizers [13]. Perhaps they could use beneficial microbes that could drastically lower the levels of soil contamination [71].

4. Conclusions

The current study focuses on the heavy metal contamination assessment of the soil in the southwest region of Saudi Arabia, which includes some parts of Al-Baha. This pollution is seen to be one of the major barriers to food security and sustainable development. This study proved that semivariogram models are a useful tool for forecasting the heavy metals’ spatial distribution maps in the investigated area. In addition, the combination of PCA and HCA produced unexpected outcomes by creating two zones within the study area, each with a different heavy metal content and pattern. The findings revealed two degrees of heavy metal pollution: moderate and strong. A moderate amount of pollution was found in roughly 63.41% of the research region. The majority of the pollution levels were above the threshold of their average concentration in the parent material and the crust of the planet. In order to decrease their reliance on chemical pesticides and fertilizers, farmers ought to employ manure and biofertilizers that differ in terms of heavy metal concentration and pattern. To prevent and regulate rising levels of heavy metals (HMs), especially those of Cu, Zn, Cd, and Pb, soil concentrations in the research region should be routinely monitored. It also shows that in areas with limited data availability, a combination of statistical and geostatistical approaches works well and is simple to apply. To improve the geographic interpretation of soil contamination, it is advised to increase the sampling size. The study concludes by urging the implementation of farm management legislation to curtail harmful human behaviors that worsen environmental contamination. Additionally, research in the future will concentrate on managing and lessening the consequences of soil pollution.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/toxics12080569/s1, Figure S1: Some cultivated crops in the study area; Figure S2: Semivariograms of studied elements; Table S1: Limit of detection of studied elements; Table S2: Classification of contamination factor(CF); Table S3: Classification of pollution load index (PLI); Table S4: KMO and Bartlett’s Test for studied variables.

Author Contributions

Conceptualization, A.S.E.-S. and M.S.S.; Methodology, M.S.S.; Software, M.S.S.; Validation, M.S.S.; Formal analysis, M.S.S.; Investigation, M.S.S.; Resources, H.A. and A.S.E.-S.; Writing—Original draft, A.O. and M.S.S.; Writing—Review and editing, A.O.; Supervision, H.A. and A.S.E.-S.; Project administration, H.A. and A.S.E.-S.; Funding acquisition, H.A. All authors have read and agreed to the published version of the manuscript.

Funding

Researchers Supporting Project number RSP2024R425, King Saud University, Riyadh, Saudi Arabia.

Data Availability Statement

All data are included in the manuscript and supplementary files.

Acknowledgments

The authors extend their appreciation to Researchers Supporting Project number RSP2024R425, King Saud University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Shokr, M.S.; Abdellatif, M.A.; El Behairy, R.A.; Abdelhameed, H.H.; El Baroudy, A.A.; Mohamed, E.S.; Rebouh, N.Y.; Ding, Z.; Abuzaid, A.S.J.A. Assessment of Potential Heavy Metal Contamination Hazards Based on GIS and Multivariate Analysis in Some Mediterranean Zones. Agronomy 2022, 12, 3220. [Google Scholar] [CrossRef]
  2. Li, C.; Zhou, K.; Qin, W.; Tian, C.; Qi, M.; Yan, X.; Han, W. A review on heavy metals contamination in soil: Effects, sources, and remediation techniques. Soil Sediment Contam. Int. J. 2019, 28, 380–394. [Google Scholar] [CrossRef]
  3. Rahmawati, S.; Iresha, F.M.; Wacano, D.; Fauziah, I.F.; Amrullah, M.A. Evaluation of heavy metal exposure to soil and paddy plant around the closed municipal solid waste landfill: Case study at Gunung Tugel Landfill, Banyumas-Central Java. IOP Conf. Ser. Mater. Sci. Eng. 2018, 299, 012012. [Google Scholar] [CrossRef]
  4. Srigirisetty, S.; Jayasri, T.; Netaji, C.J.T.R.O.I. Open dumping of municipal solid waste-Impact on groundwater and soil. Tech. Res. Organ. India 2017, 4, 26–33. [Google Scholar]
  5. Al-Robai, S.A.J.S. Ecological Risk Evaluation of Heavy Metals in Soils near a Water Dam in Baljurashi, KSA, and Their Accumulation in Dodonaea viscosa. Sustainability 2023, 15, 15646. [Google Scholar] [CrossRef]
  6. Gjorgieva Ackova, D. Heavy metals and their general toxicity on plants. Plant Sci. Today 2018, 5, 15–19. [Google Scholar] [CrossRef]
  7. Zhao, L.; Hu, G.; Yan, Y.; Yu, R.; Cui, J.; Wang, X.; Yan, Y. Source apportionment of heavy metals in urban road dust in a continental city of eastern China: Using Pb and Sr isotopes combined with multivariate statistical analysis. Atmospheric Environ. 2019, 201, 201–211. [Google Scholar] [CrossRef]
  8. Zhang, Y.; Hou, D.; O’Connor, D.; Shen, Z.; Shi, P.; Ok, Y.S.; Tsang, D.C.; Wen, Y.; Luo, M. Lead contamination in Chinese surface soils: Source identification, spatial-temporal distribution and associated health risks. Crit. Rev. Environ. Sci. Technol. 2019, 49, 1386–1423. [Google Scholar] [CrossRef]
  9. Zabin, S.A.; Foaad, M.; Al-Ghamdi, A.Y.J.H.; Assessment, E.R. Non-carcinogenic risk assessment of heavy metals and fluoride in some water wells in the Al-Baha Region, Saudi Arabia. Hum. Ecol. Risk Assess. 2008, 14, 1306–1317. [Google Scholar] [CrossRef]
  10. Liu, J.; Xu, Y.; Cheng, Y.; Zhao, Y.; Pan, Y.; Fu, G.; Dai, Y. Occurrence and risk assessment of heavy metals in sediments of the Xiangjiang River, China. Environ. Sci. Pollut. Res. 2017, 24, 2711–2723. [Google Scholar] [CrossRef] [PubMed]
  11. Zhao, X.; Gao, B.; Xu, D.; Gao, L.; Yin, S. Heavy metal pollution in sediments of the largest reservoir (Three Gorges Reservoir) in China: A review. Environ. Sci. Pollut. Res. 2017, 24, 20844–20858. [Google Scholar] [CrossRef] [PubMed]
  12. El Behairy, R.A.; El Baroudy, A.A.; Ibrahim, M.M.; Mohamed, E.S.; Rebouh, N.Y.; Shokr, M.S. Combination of GIS and multivariate analysis to assess the soil heavy metal contamination in some arid zones. Agronomy 2022, 12, 2871. [Google Scholar] [CrossRef]
  13. Alzahrani, Y.; Alshehri, F.; El-Sorogy, A.S.; Alzahrani, H. Environmental assessment of heavy metals in soils around Al-Janabeen Dam, southwest Saudi Arabia. J. King Saud Univ.-Sci. 2023, 35, 102503. [Google Scholar] [CrossRef]
  14. Senila, M. Metal and metalloid monitoring in water by passive sampling—A review. Rev. Anal. Chem. 2023, 42, 20230065. [Google Scholar] [CrossRef]
  15. Rolón-Cárdenas, G.A.; Hernández-Morales, A. Phytoremediation of Lead Present in Environment: A Review. In Lead Toxicity Mitigation: Sustainable Nexus Approaches; Springer: Cham, Switzerland, 2024; pp. 149–168. [Google Scholar]
  16. Briffa, J.; Sinagra, E.; Blundell, R. Heavy metal pollution in the environment and their toxicological effects on humans. Heliyon 2020, 6, e04691. [Google Scholar] [CrossRef] [PubMed]
  17. Zhen, J.; Pei, T.; Xie, S. Kriging methods with auxiliary nighttime lights data to detect potentially toxic metals concentrations in soil. Sci. Total. Environ. 2019, 659, 363–371. [Google Scholar] [CrossRef] [PubMed]
  18. Golden, N.; Zhang, C.; Potito, A.; Gibson, P.J.; Bargary, N.; Morrison, L. Use of ordinary cokriging with magnetic susceptibility for mapping lead concentrations in soils of an urban contaminated site. Sci. Total Environ. 2020, 20, 1357–1370. [Google Scholar] [CrossRef]
  19. Evans, N.; Van Ryswyk, H.; Los Huertos, M.; Srebotnjak, T. Robust spatial analysis of sequestered metals in a Southern California Bioswale. Sci. Total Environ. 2019, 650, 155–162. [Google Scholar] [CrossRef] [PubMed]
  20. Jin, Z.; Zhang, L.; Lv, J.; Sun, X. The application of geostatistical analysis and receptor model for the spatial distribution and sources of potentially toxic elements in soils. Environ. Geochem. Health 2021, 43, 407–421. [Google Scholar] [CrossRef] [PubMed]
  21. Shi, C.; Wang, Y. Non-parametric machine learning methods for interpolation of spatially varying non-stationary and non-Gaussian geotechnical properties. Geosci. Front. 2021, 12, 339–350. [Google Scholar] [CrossRef]
  22. Abuzaid, A.S.; Mazrou, Y.S.; El Baroudy, A.A.; Ding, Z.; Shokr, M.S. Multi-Indicator and geospatial based approaches for assessing variation of land quality in arid agroecosystems. Sustainability 2022, 14, 5840. [Google Scholar] [CrossRef]
  23. Kowalska, J.B.; Mazurek, R.; Gąsiorek, M.; Zaleski, T. Pollution indices as useful tools for the comprehensive evaluation of the degree of soil contamination—A review. Environ. Geochem. Health 2018, 40, 2395–2420. [Google Scholar] [CrossRef] [PubMed]
  24. Jahin, H.S.; Abuzaid, A.S.; Abdellatif, A.D. Using multivariate analysis to develop irrigation water quality index for surface water in Kafr El-Sheikh Governorate, Egypt. Environ. Technol. Innov. 2020, 17, 100532. [Google Scholar] [CrossRef]
  25. Abuzaid, A.S.; Jahin, H.S. Combinations of multivariate statistical analysis and analytical hierarchical process for indexing surface water quality under arid conditions. J. Contam. Hydrol. 2022, 248, 104005. [Google Scholar] [CrossRef] [PubMed]
  26. Al-Saud, M.; Teutsch, G.; Schüth, C.; Rausch, R. Challenges for an integrated groundwater management in the Kingdom of Saudi Arabia. Int. J. Water Resour. Arid. Environ. 2011, 1, 65–70. [Google Scholar]
  27. Al-Barakah, F.N.; Aly, A.A.; Abaakhel, E.H.; Al-Rizkid, A.M.; Alghamdi, A.G.; Al-Sewailem, M.S. Comparison and hydrochemical characterization of groundwater resources in the Arabian Peninsula: A case study of Al-Baha and Al-Qassim in Saudi Arabia. Water Resour. 2020, 47, 877–891. [Google Scholar] [CrossRef]
  28. Greenwood, W.R. Geology of the Al ‘Aqiq Quadrangle, Sheet 20/41 D, Kingdom of Saudi Arabia; Directorate General of Mineral Resources; Geological Survey: Jiddah, Saudi Arabia, 1975.
  29. Shanti, A.M. Geology of the Arabian shield of Saudi Arabia; Scientific Publishing Centre, King Abdulaziz University: Jeddah, Saudi Arabia, 2009. [Google Scholar]
  30. US EPA. Method 3050B: Acid Digestion of Sediments, Sludges, and Soils; United States Environmental Protection Agency: Washington, DC, USA, 1996.
  31. Trabzuni, D.M.; Ahmed, S.E.B.; Abu-Tarboush, H.M. Chemical composition, minerals and antioxidants of the heart of date palm from three Saudi cultivars. Food Nutr. Sci. 2014, 5, 1379. [Google Scholar] [CrossRef]
  32. El-Sorogy, A.S.; Youssef, M.; Al-Kahtany, K. Integrated assessment of the Tarut Island coast, Arabian Gulf, Saudi Arabia. Environ. Earth Sci. 2016, 75, 1336. [Google Scholar] [CrossRef]
  33. Hakanson, L. An ecological risk index for aquatic pollution control. A sedimentological approach. Water Res. 1980, 14, 975–1001. [Google Scholar] [CrossRef]
  34. Wedepohl, K.H. The composition of the continental crust. Geochim. Cosmochim. Acta 1995, 59, 1217–1232. [Google Scholar] [CrossRef]
  35. Tomlinson, D.; Wilson, J.; Harris, C.; Jeffrey, D. Problems in the assessment of heavy-metal levels in estuaries and the formation of a pollution index. Helgoländer Meeresunters. 1980, 33, 566–575. [Google Scholar] [CrossRef]
  36. Chon, H.-T.; Cho, C.-H.; Kim, K.-W.; Moon, H.-S. The occurrence and dispersion of potentially toxic elements in areas covered with black shales and slates in Korea. Appl. Geochem. 1996, 11, 69–76. [Google Scholar] [CrossRef]
  37. Abuzaid, A.S.; Jahin, H.S.; Shokr, M.S.; El Baroudy, A.A.; Mohamed, E.S.; Rebouh, N.Y.; Bassouny, M.A. A Novel Regional-Scale Assessment of Soil Metal Pollution in Arid Agroecosystems. Agronomy 2023, 13, 161. [Google Scholar] [CrossRef]
  38. Said, M.E.S.; Ali, A.M.; Borin, M.; Abd-Elmabod, S.K.; Aldosari, A.A.; Khalil, M.M.; Abdel-Fattah, M.K. On the use of multivariate analysis and land evaluation for potential agricultural development of the northwestern coast of Egypt. Agronomy 2020, 10, 1318. [Google Scholar] [CrossRef]
  39. Fan, S.; Wang, X.; Lei, J.; Ran, Q.; Ren, Y.; Zhou, J. Spatial distribution and source identification of heavy metals in a typical Pb/Zn smelter in an arid area of northwest China. Hum. Ecol. Risk Assess. Int. J. 2019, 25, 1661–1687. [Google Scholar] [CrossRef]
  40. Santra, P.; Kumar, M.; Panwar, N.; Yadav, R. Digital soil mapping: The future need of sustainable soil management. In Geospatial Technologies for Crops and Soils; Springer: Singapore, 2021; pp. 319–355. [Google Scholar]
  41. Abdellatif, M.A.; El Baroudy, A.A.; Arshad, M.; Mahmoud, E.K.; Saleh, A.M.; Moghanm, F.S.; Shaltout, K.H.; Eid, E.M.; Shokr, M.S. A GIS-based approach for the quantitative assessment of soil quality and sustainable agriculture. Sustainability 2021, 13, 13438. [Google Scholar] [CrossRef]
  42. Goenster-Jordan, S.; Jannoura, R.; Jordan, G.; Buerkert, A.; Joergensen, R.G. Spatial variability of soil properties in the floodplain of a river oasis in the Mongolian Altay Mountains. Geoderma 2018, 330, 99–106. [Google Scholar] [CrossRef]
  43. Sebei, A.; Chaabani, A.; Abdelmalek-Babbou, C.; Helali, M.A.; Dhahri, F.; Chaabani, F. Evaluation of pollution by heavy metals of an abandoned Pb-Zn mine in northern Tunisia using sequential fractionation and geostatistical mapping. Environ. Sci. Pollut. Res. 2020, 27, 43942–43957. [Google Scholar] [CrossRef] [PubMed]
  44. Mallik, S.; Mishra, U.; Paul, N. Groundwater suitability analysis for drinking using GIS based fuzzy logic. Ecol. Indic. 2021, 121, 107179. [Google Scholar] [CrossRef]
  45. Dad, J.M.; Shafiq, M.U. Spatial variability and delineation of management zones based on soil micronutrient status in apple orchard soils of Kashmir valley, India. Environ. Monit. Assess. 2021, 193, 797. [Google Scholar] [CrossRef] [PubMed]
  46. DOE International Affairs. National Norms and Standards for the Remediation of Contaminated Land and Soil Quality in the Republic of South Africa; Department of Environmental Affairs (DEA), National Environmental Management: Pretoria, South Africa, 2013.
  47. US Department of Health and Human Services; Public Health Service; Agency for Toxic Substances and Disease Registry. Toxicological Profile for Arsenic; US Department of Health and Human Services: Atlanta, GA, USA, 2000; pp. 1–14.
  48. Kubier, A.; Wilkin, R.T.; Pichler, T.J. Cadmium in soils and groundwater: A review. Appl. Geochem. 2019, 108, 104388. [Google Scholar] [CrossRef] [PubMed]
  49. ATSDR. Toxicological Profile for Cobalt; U.S. Department of Health and Human Services: Washington, DC, USA, 2004; p. 486.
  50. Mahey, S.; Kumar, R.; Sharma, M.; Kumar, V.; Bhardwaj, R. A critical review on toxicity of cobalt and its bioremediation strategies. SN Appl. Sci. 2020, 2, 1279. [Google Scholar] [CrossRef]
  51. Zaborowska, M.; Kucharski, J.; Wyszkowska, J. Biological activity of soil contaminated with cobalt, tin, and molybdenum. Environ. Monit. Assess. 2016, 188, 398. [Google Scholar] [CrossRef]
  52. Fishel, F. Pesticide toxicity profile: Copper-based pesticides. EDIS 2005, 2005. [Google Scholar] [CrossRef]
  53. Wu, R.; Yao, F.; Li, X.; Shi, C.; Zang, X.; Shu, X.; Liu, H.; Zhang, W. Manganese pollution and its remediation: A review of biological removal and promising combination strategies. Microorganisms 2022, 10, 2411. [Google Scholar] [CrossRef] [PubMed]
  54. Begum, W.; Rai, S.; Banerjee, S.; Bhattacharjee, S.; Mondal, M.H.; Bhattarai, A.; Saha, B. A comprehensive review on the sources, essentiality and toxicological profile of nickel. RSC Adv. 2022, 12, 9139–9153. [Google Scholar] [CrossRef] [PubMed]
  55. Ankush; Lamba, S.; Ritambhara; Diwedi, A.; Kumar, S.; Singh, V. Source and Distribution of Lead in Soil and Plant—A Review. In Lead Toxicity: Challenges and Solution; Springer: Cham, Switzerland, 2023; pp. 3–16. [Google Scholar]
  56. Zhang, B.; Zhao, H.; Shi, C.; Zhou, S.; Ni, J. Simultaneous removal of sulfide and organics with vanadium (V) reduction in microbial fuel cells. J. Chem. Technol. Biotechnol. 2009, 84, 1780–1786. [Google Scholar] [CrossRef]
  57. Watt, J.A.; Burke, I.T.; Edwards, R.A.; Malcolm, H.M.; Mayes, W.M.; Olszewska, J.P.; Pan, G.; Graham, M.C.; Heal, K.V.; Rose, N.L. Vanadium: A re-emerging environmental hazard. Environ. Sci. Technol. 2018, 52, 11973–11974. [Google Scholar] [CrossRef] [PubMed]
  58. Khan, S.; Kazi, T.G.; Kolachi, N.F.; Baig, J.A.; Afridi, H.I.; Shah, A.Q.; Kumar, S.; Shah, F. Hazardous impact and translocation of vanadium (V) species from soil to different vegetables and grasses grown in the vicinity of thermal power plant. J. Hazard. Mater. 2011, 190, 738–743. [Google Scholar] [CrossRef] [PubMed]
  59. Hao, L.; Zhang, B.; Feng, C.; Zhang, Z.; Lei, Z.; Shimizu, K.; Cao, X.; Liu, H.; Liu, H. Microbial vanadium (V) reduction in groundwater with different soils from vanadium ore mining areas. Chemosphere 2018, 202, 272–279. [Google Scholar] [CrossRef] [PubMed]
  60. Zhang, B.; Hao, L.; Tian, C.; Yuan, S.; Feng, C.; Ni, J.; Borthwick, A. Microbial reduction and precipitation of vanadium (V) in groundwater by immobilized mixed anaerobic culture. Bioresour. Technol. 2015, 192, 410–417. [Google Scholar] [CrossRef] [PubMed]
  61. Bodek, I.; Lyman, W.J.; Reehl, W.F.; Rosenblatt, D.H. Environmental Inorganic Chemistry: Properties, Processes, and Estimation Methods; Pergamon Press: New York, NY, USA, 1988; Volume 10. [Google Scholar]
  62. Swartjes, F.A. Dealing with Contaminated Sites: From Theory towards Practical Application; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  63. Tóth, G.; Hermann, T.; Da Silva, M.R.; Montanarella, L. Heavy metals in agricultural soils of the European Union with implications for food safety. Environ. Int. 2016, 88, 299–309. [Google Scholar] [CrossRef] [PubMed]
  64. Mudelsee, M. Correlation. In Climate Time Series Analysis; Springer: Berlin/Heidelberg, Germany, 2010; Volume 42, pp. 285–338. [Google Scholar]
  65. Uren, N.C. Cobalt and Manganese; Springer: Dordrecht, The Netherlands, 2013; pp. 335–366. [Google Scholar]
  66. Shtangeeva, I.; Alber, D.; Bukalis, G.; Stanik, B.; Zepezauer, F. Multivariate statistical analysis of nutrients and trace elements in plants and soil from northwestern Russia. Plant Soil 2009, 322, 219–228. [Google Scholar] [CrossRef]
  67. Ekpo, F.; Nzegblue, E.; Asuquo, M. A comparative study of the influence of heavy metals on soil and crops growing within quarry environment at Akamkpa, Cross River State, Nigeria. Glob. J. Agric. Sci. 2012, 11, 1–4. [Google Scholar] [CrossRef]
  68. Ağca, N. Spatial distribution of heavy metal content in soils around an industrial area in Southern Turkey. Arab. J. Geosci. 2015, 8, 1111–1123. [Google Scholar] [CrossRef]
  69. Smith, A.; Cullis, B.; Thompson, R. Exploring variety-environment data using random effects AMMI models with adjustments for spatial field trend: Part 1: Theory. In Quantitative Genetics, Genomics and Plant Breeding; CABI Publishing: Wallingford, UK, 2002; pp. 323–335. [Google Scholar]
  70. Rai, V.; Vajpayee, P.; Singh, S.N.; Mehrotra, S. Effect of chromium accumulation on photosynthetic pigments, oxidative stress defense system, nitrate reduction, proline level and eugenol content of Ocimum tenuiflorum L. Plant Sci. 2004, 167, 1159–1169. [Google Scholar] [CrossRef]
  71. Mohamed, E.S.; Jalhoum, M.E.; Hendawy, E.; El-Adly, A.M.; Nawar, S.; Rebouh, N.Y.; Saleh, A.; Shokr, M.S. Geospatial evaluation and bio-remediation of heavy metal-contaminated soils in arid zones. Front. Environ. Sci. 2024, 12, 1381409. [Google Scholar] [CrossRef]
Figure 1. Site description and sample distribution.
Figure 1. Site description and sample distribution.
Toxics 12 00569 g001
Figure 2. Kriging interpolation maps of studied heavy metals. As (a), Cd (b), Co (c), Cr (d), Cu (e), Mn (f), Ni (g), Pb (h), V (i), Fe (j) and Zn (k).
Figure 2. Kriging interpolation maps of studied heavy metals. As (a), Cd (b), Co (c), Cr (d), Cu (e), Mn (f), Ni (g), Pb (h), V (i), Fe (j) and Zn (k).
Toxics 12 00569 g002aToxics 12 00569 g002bToxics 12 00569 g002c
Figure 3. Scree plot of different principal components.
Figure 3. Scree plot of different principal components.
Toxics 12 00569 g003
Figure 4. Biplot of principal components.
Figure 4. Biplot of principal components.
Toxics 12 00569 g004
Figure 5. Hierarchical cluster analysis (HCA) extracted from PCA; blue expresses C1, and red represents C2.
Figure 5. Hierarchical cluster analysis (HCA) extracted from PCA; blue expresses C1, and red represents C2.
Toxics 12 00569 g005
Figure 6. CF for each studied heavy metal in cluster 1 (a) and cluster 2 (b) (low contamination(LC), moderate contamination (MC), considerable contamination(CC), and VH (very high contamination).
Figure 6. CF for each studied heavy metal in cluster 1 (a) and cluster 2 (b) (low contamination(LC), moderate contamination (MC), considerable contamination(CC), and VH (very high contamination).
Toxics 12 00569 g006
Figure 7. Map of PLI in the study area.
Figure 7. Map of PLI in the study area.
Toxics 12 00569 g007
Table 1. Descriptive statistics of the studied heavy metals.
Table 1. Descriptive statistics of the studied heavy metals.
Recommended Concentrations
(mg kg−1)
ElementsConcentrationsNMin.Max.MeanSkewnessKurtosisUpper Crust Concentrations [34]Recommended Concentrations [46]
Asmg kg−1320.122.07.1 ± 6.30.8−0.125.8
Cd0.11.10.8 ± 0.2−3.316.30.17.5
Co935.020.6 ± 6.60.4−0.611.6300
Cr13191.058.2 ± 34.42.06.1356.5
Cu30110.061.1 ± 20.50.90.014.316
Mn6841565.01015.2 ± 2120.60.0527740
Ni1272.042.3 ± 17.20.0−0.96291
Pb427.010.1 ± 5.51.83.51720
V39163.089.8 ± 330.7−0.553150
Fe31,70070,70045,846.9 ± 10,041.70.6−0.559,100n/a
Zn60191.096.7 ± 24.71.96.252240
Table 2. The spatial data modelling semivariogram.
Table 2. The spatial data modelling semivariogram.
ElementsConcentrationsModelNuggetPartial
Sill
SillNugget/
Sill
SPDMSERMSSE
Asmg kg−1Stable0.50.40.900.55Moderate−0.0080.98
CdCircular0.510.410.920.55Moderate0.0341.22
CoCircular0.540.4610.54Moderate0.0050.99
CrExponential0.410.641.050.39Moderate0.0111.07
CuSpherical0.70.270.970.72Moderate−0.0081.008
MnGaussian0.500.450.950.52Moderate0.0061.004
NiCircular0.510.480.990.51Moderate0.0051.01
PbGaussian0.230.530.760.30Moderate−0.0030.99
VExponential0.380.631.010.37Moderate0.0170.98
FeExponential0.010.950.960.01Strong0.0190.91
ZnStable0.550.290.840.65Moderate−0.031.05
Table 3. Pearson correlation between different studied metals.
Table 3. Pearson correlation between different studied metals.
VariablesAsCdCoCrCuMnNiPbVZnFe
As1−0.007−0.3090.099−0.425 *−0.0140.483 **−0.060−0.367 *−0.360 *−0.360 *
Cd−0.00710.134−0.2540.2580.320−0.2050.0470.0180.0680.086
Co−0.3090.13410.576 **0.642 **0.380 *0.399 *−0.2400.824 **0.0950.920 **
Cr0.099−0.2540.576 **10.2350.0040.716 **−0.1370.469 **−0.0480.466 **
Cu−0.425 *0.2580.642 **0.23510.238−0.0500.0330.474 **0.391 *0.539 **
Mn−0.0140.3200.380 *0.0040.23810.029−0.1300.046−0.0240.211
Ni0.483 **−0.2050.399 *0.716 **−0.0500.0291−0.2090.251−0.1560.402 *
Pb−0.0600.047−0.240−0.1370.033−0.130−0.2091−0.3110.320−0.277
V−0.367 *0.0180.824 **0.469 **0.474 **0.0460.251−0.3111−0.0180.871 **
Zn−0.360 *0.0680.095−0.0480.391 *−0.024−0.1560.320−0.01810.189
Fe−0.360 *0.0860.920 **0.466 **0.539 **0.2110.402 *−0.2770.871 **0.1891
Note: * p < 0.05; ** p < 0.01. Bold = significance.
Table 4. Eigenvalue of factors extracted by PCA.
Table 4. Eigenvalue of factors extracted by PCA.
Factor Loading
F1F2F3F4F5F6F7F8F9F10F11
Eigenvalue3.9472.2651.4051.1390.6850.5890.4250.2990.1680.0570.022
Variability (%)35.87720.59412.77110.3586.2235.3573.8632.7151.5240.5180.199
Cumulative %35.87756.47269.24379.60185.82491.18095.04497.75999.28299.801100.000
Bold = significance.
Table 5. Factors extracted by PCA for studied heavy metals.
Table 5. Factors extracted by PCA for studied heavy metals.
ElementsConcentrationsComponent
1234
As mg kg −1−0.344−0.6970.2760.402
Cd0.0820.4200.6570.256
Co0.9720.0180.0950.045
Cr0.630−0.525−0.2740.256
Cu0.6650.490−0.0140.187
Mn0.2880.1490.7030.250
Ni0.431−0.766−0.0650.367
Pb−0.3000.369−0.3700.591
V0.874−0.013−0.077−0.302
Zn0.1510.584−0.4120.409
Fe0.9450.035−0.027−0.066
Note: Bold values match for each variable to the factor for which the squared cosine is the largest.
Table 6. Quantitative data of the studied elements for two clusters (C1 and C2).
Table 6. Quantitative data of the studied elements for two clusters (C1 and C2).
ClustersStatisticsAsCdCoCrCuMnNiPbVFeZn
mg kg−1
C1n13
Min0.10.12213527691458549,90077
Max151.1351911101565701616370,700144
Mean4.19 ± 4.54 a0.86 ± 0.25 a27.08 ± 4.21 a75.08 ± 42.51 a74.77 ± 20.21 a1109.31 ± 178.73 a48.31 ± 18.08 a8.62 ± 3.50 a116.23 ± 27.2 a56,276 ± 6075.24 a103.08 ± 17.63 a
Skewness9.16 ± 6.57 b0.82 ± 0.05 a16.16 ± 3.51 b46.58 ± 21.91 b51.74 ± 14.93 b950.84 ± 212.88 b38.11 ± 15.80 a11.11 ± 6.36 a71.79 ± 22.92 b38,710.53 ± 4090.35 b92.32 ± 28.11 a
Kurtosis1.138−2.6160.3871.5880.6020.861−0.5470.9860.3841.0980.936
C2n19
Min 0.1 0.7 9 17 30 684 12 4 39 31,700 60
Max 22 0.9 24 107 91 1425 72 27 139 47,300 191
Mean 9.16 ± 6.57 b 0.82 ± 0.05 a 16.16 ± 3.51 b 46.58 ± 21.91 b 51.74 ± 14.93 b 950.84 ± 212.88 b 38.11 ± 15.80 a 11.11 ± 6.36 a 71.79 ± 22.92 b 38,710.53 ± 4090.35 b 92.32 ± 28.11 a
Skewness 0.57 0.22 0.24 1.34 1.04 1.01 0.28 1.57 1.7 0.38 2.46
Kurtosis −0.72 0.31 0.79 2.33 1.24 0.01 −0.18 2.09 3.42 −0.27 8.54
Note: There is a considerable discrepancy between the means of variables with different letters.
Table 7. The mean concentration of heavy metals in every level of pollution in the research area.
Table 7. The mean concentration of heavy metals in every level of pollution in the research area.
Pollution LevelAsCdCoCrCuMnNiPbVFeZnArea km2, %
mg kg−1
M5.28 ± 5.830.81 ± 0.1918.65 ± 6.2245.15 ± 23.2560.55 ± 23.74972.30 ± 223.5033.45 ± 14.1110.05 ± 5.1384.15 ± 30.7243,245.00 ± 8942.9597.40 ± 30.05 910.41, (63.41%)
S10.25 ± 5.920.883 ± 0.0823.83 ± 6.1779.83 ± 39.6762.00 ± 14.421086.75 ± 177.2356.92 ± 10.9010.17 ± 6.2099.33 ± 35.6850,183.33 ± 10,639.96 95.50 ± 12.34 525.19, (36.58%)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alzahrani, H.; El-Sorogy, A.S.; Okok, A.; Shokr, M.S. GIS- and Multivariate-Based Approaches for Assessing Potential Environmental Hazards in Some Areas of Southwestern Saudi Arabia. Toxics 2024, 12, 569. https://doi.org/10.3390/toxics12080569

AMA Style

Alzahrani H, El-Sorogy AS, Okok A, Shokr MS. GIS- and Multivariate-Based Approaches for Assessing Potential Environmental Hazards in Some Areas of Southwestern Saudi Arabia. Toxics. 2024; 12(8):569. https://doi.org/10.3390/toxics12080569

Chicago/Turabian Style

Alzahrani, Hassan, Abdelbaset S. El-Sorogy, Abdurraouf Okok, and Mohamed S. Shokr. 2024. "GIS- and Multivariate-Based Approaches for Assessing Potential Environmental Hazards in Some Areas of Southwestern Saudi Arabia" Toxics 12, no. 8: 569. https://doi.org/10.3390/toxics12080569

APA Style

Alzahrani, H., El-Sorogy, A. S., Okok, A., & Shokr, M. S. (2024). GIS- and Multivariate-Based Approaches for Assessing Potential Environmental Hazards in Some Areas of Southwestern Saudi Arabia. Toxics, 12(8), 569. https://doi.org/10.3390/toxics12080569

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop