Next Article in Journal
A Dialogical Approach to Readiness for Change towards Sustainability in Higher Education Institutions: The Case of the SDGs Seminars at the Universidad Politécnica de Madrid
Previous Article in Journal
Reading Urban Green Morphology to Enhance Urban Resilience: A Case Study of Six Southern European Cities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soil Degradation and Restoration in Southwestern Saudi Arabia through Investigation of Soil Physiochemical Characteristics and Nutrient Status as Indicators

1
College of Food and Agricultural Sciences, King Saud University, P.O. Box 2460, Riyadh 11451, Saudi Arabia
2
Department of Plant Protection, College of Food and Agricultural Science, King Saud University, P.O. Box 2460, Riyadh 11451, Saudi Arabia
3
Department of Soil Sciences, College of Food and Agricultural Sciences, King Saud University, P.O. Box 2460, Riyadh 11451, Saudi Arabia
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(16), 9169; https://doi.org/10.3390/su13169169
Submission received: 12 July 2021 / Revised: 3 August 2021 / Accepted: 9 August 2021 / Published: 16 August 2021

Abstract

:
Soil degradation is among the major global concerns resulting in declined soil fertility, loss of soil organic carbon (SOM) and nutrients, and reduced pollinator visitations. The current study was performed to explore soil quality and fertility after restoration in southwestern Saudi Arabia, as new plantation could alter soil properties. A total of 78 soil samples were collected from restored and control (without any restoration) areas from Al-Jenebeen (AJ), Queen Rearing Center (QRC), and Shekiran (SK) from depths of 0–60 cm. Results exhibited that sandy loam was the most dominant textural class in all the studied sites (~60%), followed by sandy clay loam. All the soils were slightly alkaline and nonsaline. Due to restoration practices, the pH was increased in AJ and QRC, while decreased in SK areas. With restoration, EC in AJ and QRC sites decreased, but increased in SK. Restoration practices resulted in reduced SOM in AJ and QRC, while increased in SK. Increased SOM in SK may improve floral traits after restoration. Overall, restored SK possessed higher SOM and nutrients with reduced pH, and, hence, could be a potential source for honeybee flora. Further in depth analyses of restoration activities are needed to explore the long term effects on land rehabilitation.

1. Introduction

Soil/forest degradation is one of the major problems of agricultural soils worldwide [1]. Global soil resources are at risk of degradation due to anthropogenic and natural factors, consequently resulting in lower crop yields [2]. According to some estimates, 33% of the world land surface has been degraded, resulting in a decline of about 60% of ecosystem services [3,4]. Soil/forest degradation could be elicited by some chemical, physical, biological, and/or ecological processes, which subsequently deteriorate soil quality and fertility [5]. Therefore, soil quality is of significant importance for human health and the environment, due to its significant implication for ecosystems. Consequently, the degradation of soil is disturbing the natural ecological balance, economic growth, and ecosystem sustainability. Therefore, measures must be taken to protect soil quality and restore the degraded lands to enhance ecosystem services such as food, fiber, water filtration, and elemental cycling.
Soil degradation can be evaluated through investigating the chemical, physical, biological, and hydrological characteristics of soil [6,7]. Variations to these properties could disturb the natural ecosystem balance and soil health. The deterioration of soil structure, aggregates, and pore geometry could generate crusts and compaction, consequently enhancing soil erosion, and, thus, are referred to as physical degradation. On the other hand, soil salinization, acidification, toxification, and nutrient loss could reduce soil fertility and quality and is thus referred as chemical degradation. Likewise, the loss of organic matter, declined biodiversity, and enhanced greenhouse gas emissions is referred as the biological degradation. According to an estimate, about 85% of soils have been degraded, globally, due to soil erosion [8]. Therefore, physical degradation is considered as a limiting factor in enhancing crop production and sustaining natural ecosystems. Moreover, physical degradation is the major cause of lower nutrient retention and reduced water holding capacity, especially in arid and semiarid regions of the world. However, soil degradation alone is considered to be counterproductive. Hence, a set of appropriate soil indicators should be selected to identify the problem, which is a challenge to perform in field conditions. In this context, the selection of appropriate soil quality indicators, assessing the soil quality, and sustainable land management strategies are of significant importance for soil protection and sustainable ecosystems.
Forests degradation due to extensive agricultural practices and climate has also resulted in a reduction in biodiversity [9,10,11]. It has been observed that populations of wild bees are dropping rapidly [12]. Specifically, the decline in populations of pollinators is of great concern nowadays. Therefore, efforts are being made by various researchers to conserve and enhance the populations of the wild bees in order to streamline the pollination process. One of these restoration methods involves the development of floral, woody, and shrub based landscapes in agricultural areas [13]. However, these restoration practices are costly and time consuming, hence, their adoption among growers is very low. Therefore, the restoration of degraded lands on larger scales could serve as an efficient conservation method for sustainable ecosystems and habitat rehabilitation for wild bees. However, the assessment of soil quality and restoration extent is a complex process due to the involvement of various factors, and requires detailed investigations into chemical, physical, and biological characteristics of soil, over time. Hence, keeping in mind all the above mentioned facts, a study was conducted with the objectives (i) to analyze the physical and chemical characteristics of the soil collected from restored and control (without restoration activities) areas in southwest Saudi Arabia (Al-baha Region, Baljurashi district), and (ii) to identify soil nutrient status and soil degradation to understand the major soil related problems/constraints for soil rehabilitation.

2. Materials and Methods

2.1. Study Site Description

The study was carried out in forests of southwest Saudi Arabia (Al-baha Region, Baljurashi district). The landscape constitutes the remainder of the forest stock of the country [14]. Due to the geographical diversity, this land feature has diverse climates, thus favoring the growth and flowering of greatly diversified bee plant genotypes rich in nectar and pollen, including Acacia species. As a result, most of the beekeeping activities of the country are concentrated in this region. The area is categorized as being mid to high land physiographic and climatic conditions. The humidity ranges from 52–67%, with rainfall from 229–581 mm/annum. The mean annual temperature is 22.9 °C. The major land use of the area is characterized by rangeland where free grazing and browsing is a common practice. The rangeland is also comprised of fallow land deserted decades ago.
For current study, 3 areas (Al-Jenebeen (AJ), Queen Rearing Center (QRC), and Shekiran (SK)), where restoration activities have already been practiced, were selected from the study site. The AJ and QRC are relatively flat areas, while the SK is more of hilly landscape. AJ is known for open browsing/grazing land, whereas the QRC and SK possess limited browsing/grazing land in the area despite relatively good plant cover. Similar weather prevails in QRC and SK areas, while AJ differs from these two as it receives less rainfall and is relatively warmer.

2.2. Restoration Activities

Despite its relatively good vegetation coverage, the area is severely degraded and only relict woodland patches exist. It was seen that AJ area was relatively more degraded in terms of vegetation cover and composition as compared to QRC and SK. These situations are believed to easily motivate stakeholders in the rehabilitation of the vegetation (apicultural landscapes) of the area, for both environmental and economic reasons. Therefore, some restoration activities have been conducted in the area already. These restoration activities include fencing; plantation of local plants, especially acacia species; and installation of drip irrigation systems. The fencing was developed to avoid interferences, trampling, and overgrazing of new sprouts, while plantation of new plants and drip irrigation installation were performed to enhance vegetation cover.

2.3. Site Profiling and Soil Samples Collection

The soil samples were taken from both restored and control sites located in study areas, i.e., AJ, QRC, and SK. From both groups, 3 soil profiles were dug and the samples were collected. Soil profiles were created on different localities in the study areas and samples were collected from various depths, i.e., 0–20 cm, 20–40 cm, and 40–60 cm. Additionally, some surface and composite samples were also collected to increase the possibility of analyzing the samples from different angles. Composite samples were created from each soil profile layer for both treatment groups in each site. In this way, a total of 78 samples were collected from all studied areas, comprising of different land use types, soil depths, and localities, as shown in Table 1.

2.4. Soil Samples Processing and Analyses

The collected soil samples were brought to the laboratory, air dried, ground, and passed through 2 mm screen. Thereafter, these samples were subjected to various soil analyses.

2.4.1. Soil Physiochemical Analyses

The pH of the soil was determined in a 1:1 (w/v) ratio suspension of soil to water by using a portable pH meter (pH 315i, WTW, Germany). Likewise, the filtrate of 1:1 suspension was used to analyze soil electrical conductivity (EC) by using a conductivity meter (EC meter: 4510, Jenway, UK). Cations and anions of the soil were measured by using the standard procedures, as reported earlier by [15]. Soil particle size distribution was determined by using the hydrometer method, while the texture was determined with the help of United States Department of Agriculture’s textural triangle for soil taxonomy [16]. The soils collected from the study area belong to lithosols and yermosols types. The contents of calcite (CaCO3) in collected samples were analyzed through calcimeter method [17]. Briefly, 2.0 g of soil was reacted with 0.1 N hydrochloric acid (HCl) under an airtight environment and the evolved CO2 was measured. Then, calcite percentage was determined through comparing the emission of CO2 from the sample with that of pure CaCO3. Soil organic matter was detected by following Walkley and Black (1934) method [18]. Briefly, 0.5 g soil was taken in round neck flask. Then, 10 mL of 1 N potassium dichromate solution was added, followed by 20 mL concentrated sulfuric acid. Sample was allowed to stand for 30 min. Later, 200 mL of distilled water was added to the sample followed by 3–4 drops of diphenylamine indicator. Sample was titrated with 0.5 M ferrous ammonium sulfate solution until the color changed from violet–blue to reddish brown. Additionally, a blank sample was run without soil sample.

2.4.2. Available Nutrients Analyses

Plant available phosphorus (P) and micronutrients, i.e., copper (Cu), iron (Fe), zinc (Zn), and manganese (Mn), in soil were extracted using AB-DTPA (ammonium bicarbonate-diethylenetriaminepentaacetic acid) extraction method [19]. Briefly, a mixture of 1 M ammonium bicarbonate and 0.005 M DTPA with pH of 7.6 was used to extract P and micronutrients through 1:2 soil–AB-DTPA (w/v ratio) suspension. The concentration of K in the soil samples was analyzed by using flame photometer. The concentration of P in AB-DTPA extract was measured through colorimetry method by using UV/VIS spectrophotometer (Lambda EZ 150, PerkinElmer, Norwalk, CT, USA). The contents of micronutrients in AB-DTPA extract were determined by using Inductivity coupled plasma optical emission spectrophotometer (ICP-OES, PerkinElmer Optima 4300 DV, Akron, OH, USA).

2.5. Data Analyses and Statistics

The soil analyses data were statistically analyzed by using Microsoft Excel and Statistics 8.01 program. Descriptive statistics was used to calculate minimum, maximum, mean, skewness, kurtosis, and coefficient of variation (CV) of the data. Pearson correlations were calculated to explore the relation of various soil properties and nutrients concentrations with each other.

3. Results and Discussion

3.1. Soil Physiochemical Properties as Indicators for Degradation and Restoration

3.1.1. Soil Texture

The variations in soil texture of the different soil samples collected from the control and restored areas of the studied sites are shown in Figure 1. The results exhibited that sandy loam was the most dominant textural class in all the studied sites (~60%), followed by sandy clay loam. Overall, QRC has more samples in the sandy loam textural class, followed by SK and AJ. The loamy sand class was present in QRC-C, QRC-R, and SK-C, while absent in other collected samples.
Likewise, loam can be seen only in QRC-R, AJ-C, and AJ-R, whereas clay loam was present only in QRC-C and AJ-R. It was interesting to notice that the proportion of sandy loam soils was higher in restored areas as compared to the control areas. The restored areas of QRC showed the occurrence of loam soil, which was not present in the control areas. Likewise, clay loam soils were seen in the control of QRC, and absent in its restored areas. The control areas of SK had about 16% samples of loamy sand, which disappeared in restored areas. Similarly, clay loam soils were seen in restored areas of AJ, while these were absent in AJ-C. Therefore, variations in control and restored areas evidenced the effects of restoration efforts on the physical properties of the soils in the studied area. Owing to the significant importance of soil texture, various researchers have used it as an indicator for assessing soil quality [20,21]. Various agricultural practices, especially tillage and plowing, could affect the texture of soils. Likewise, natural processes, as well as irrigation, could move the clays downward and may result in their accumulation in deep soil horizons, subsequently changing the proportions of sand, silt, and clay. Therefore, soil texture is an important parameter to investigate soil type and quality [1,22]. The predominance of sand based texture indicated lower water holding capacity, poor soil structure, and lower nutrient retention. QRC-C and AJ-R possess relatively higher contents of clay, indicating a slightly better soil structure as compared to the rest of the soil types. Therefore, these changes in the distribution of silt, clay, and sand fractions, along with changing the land-use type and locality, indicated the difference in erosion and sedimentation rates.

3.1.2. Soil pH and EC

Soil pH plays a critical role in the dissolution and availability of nutrients as well as for proper microbial activities in soil. Similarly, soil EC is an important factor for plant nutrient availability and optimum crop growth. Therefore, soil pH and EC are considered important indicators to assess the extent of soil quality. The descriptive statistical analyses for pH and EC of the collected soil samples are shown in Table 2. It can be seen that all the soils were slightly alkaline, with a mean pH range of 7.74 to 8.15. Overall, the mean pH values of AJ were higher (8.07–8.15), followed by SK (7.81–8.01) and QRC (7.74–7.99). Results showed that means increased from 8.07 in the control to 8.15 in restored areas and 7.74 in the control to 7.99 in restored areas in AJ and QRC, respectively, whereas, it decreased from 8.01 in the control to 7.81 in restored areas in SK. The increment in soil pH due to restoration processes could be due to the release of saline products by the plant’s branches and leaves, as more plants were grown during restoration [23]. Additionally, the recycling of some organic components during the restoration process might have enhanced the soil pH. Contrarily, the pH of restored areas located in SK sites decreased to 7.81 from 8.01 in the control. This decrease in pH could be caused by the addition and mineralization of organic compounds due to restoration activities [24]. This reduction in pH in SK areas after restoration could help in enhancing vegetation growth, which, in turn, could increase floral resources for honeybees. Furthermore, the decrease in pH could be owing to the uptake of exchangeable cations by plant or microorganisms in soil. The degree of variability, as indicated by CV values, was in the range of 1.99 to 3.75. All the pH data was well skewed with lower kurtosis (−1.438–−0.041), indicating a normal distribution. The differences in pH values were statistically significant, except for QRC-R and SK-C which were nonsignificantly different. The pH for AJ-R was significantly higher (8.15) as compared to all the other samples. Similarly, significant differences were also seen in EC results, where AJ-C showed the highest value (0.45 dS m−1).
The influence of soil depth on variations in soil pH can be seen in Figure 2A,C,E. Overall, the mean pH values of restored areas were higher than that of the controls for all three sites. However, this was the opposite in the topsoil of QRC and SK, where restored areas had lower pH values than the control areas. The trend for distribution of the soil pH with depth was not uniform. For instance, the mean pH values increased with increasing the depth in AJ-C, while it decreased with depth in QRC-C and SK-C. Meanwhile, this was exactly the opposite in restored areas, where mean pH values decreased with depth in AJ-R, and increased in QRC-R and SK-R. Similar contrasting results were also reported by different researchers. For example, a decrease in soil pH with an increase in soil depth was previously reported by Shukla et al. [25], whereas Zhang et al. [26] reported an increase in soil pH with an increase in the depth, suggesting that this could be due to the addition of organic amendments in the soil. Overall, the results indicated that the effects of soil depth on variations in soil pH were not significant.
Table 2 showed that the mean EC values for all the soil samples were in the range of 0.34 to 0.45 dS m−1, suggesting that EC for all the soil samples was below salinity limits, i.e., 4 dS m−1 [27]. Overall, the samples collected from AJ sites had higher EC values, followed by QRC, and SK. The results exhibited that the EC of the AJ and QRC sites reduced after restoration activities, from 0.45 to 0.36 dS m−1 and 0.43 to 0.38 dS m−1, respectively, whereas, it increased in the case of SK, from 0.34 to 0.40 dS m−1. The increase in soil EC could be the result of adding organic matter during restoration activities, while the reduction in EC clearly indicated the significant role of plantation on soil quality [28]. The values of the skewness and kurtosis showed that the distribution was not normal. Figure 2B,D,F shows variations in mean EC with soil depth. The EC values reduced with an increase in soil depth in AJ-C, while increasing with an increase in soil depth in QRC-C and SK-C. This trend was opposite in restored areas, where EC increased with soil depth in the case of AJ-R, and decreased with depth in QRC-R and SK-R. Bezborodov et al. [29] stated that EC increases with soil depth due to the presence of salts. Therefore, the distribution and extent of soil pH and EC can effectively be used as indicators for assessing soil degradation and restoration.

3.1.3. Soil Organic Matter, Soil Organic Carbon, and CaCO3

Soil organic matter (SOM) significantly influences the chemical, physical, and biological characteristics of soils [30]. It has been well documented that the depletion of SOM as a result of extensive agricultural practices has subsequently resulted in soil degradation [31,32]. Therefore, SOM is considered as one of the most promising and sensitive factors in assessing soil quality and the extent of its degradation. The descriptive statistical analyses for SOM for the samples collected from selected localities are shown in Table 3. Overall, the mean values of the contents of SOM were lower than 1.5% in all the locations. The samples collected from SK sites possessed higher SOM (0.13–2.56%), as compared to AJ (0.25–2.19%), and QRC (0.13–1.89%). It was noticed that the restoration practices resulted in reduced SOM in AJ (1.27% in AJ-C to 1.02% in AJ-R) and QRC (0.88% in QRC-C to 0.80% in QRC-R), while it increased with restoration practices in SK (1.16% in SK-C to 1.40% in SK-R). This data showed a higher variability, as shown with larger CV values (27.24–68.40). The skewness of the data was in the range of −0.316 to 0.244, while kurtosis was −1.326 to −0.304, indicating normality in the distribution of the data. Further, the changes in the contents of SOM with soil depth are shown in Figure 3A,C,E. The variations in SOM contents with soil depth were inconsistent. For instance, the mean SOM contents increased with an increase in soil depth in AJ-C and SK-R, while reducing with soil depth in the rest of the studied soil profiles. Interestingly, the topsoil of SK (both restored and the control) and AJ-C showed relatively lower SOM contents.
The variations in SOM of the collected soil samples suggested that restoration activities had positive impacts on soil rehabilitation in SK sites, as indicated by increased SOM contents. Contrarily, the SOM contents of AJ and QRC sites reduced as a result of restoration activities. The loss of SOM during restoration activities could be due to higher sand fractions. Generally, the SOM associated with clay fractions is more resistant to degradation. However, our study demonstrated that sandy loam was the most dominant textural class for all the collected samples, indicating a higher portion of sand [33]. Kalinina et al. [34] suggested that the loss of SOM could also be due to higher humus rich arable soils, texture, and mowing. To further understand this, the proportions of soil organic carbon (SOC), according to the soil evaluation criteria, are presented in Figure 4. It can be seen in the figure that none of the samples were in the ranges of high, very high and extremely high SOC contents. Most of the samples were in the extremely low to moderate SOC range. The proportion of samples with moderate SOC contents was higher after restoration in the SK area, while it was reduced to zero in both AJ and QRC after restoration. On the other hand, the samples with moderately low SOC contents increased in AJ and QRC, while reducing in SK after restoration activities. This clearly indicates the loss of SOC and SOM from the soil after restoration activities. It is well known that plants cannot absorb SOM/SOC from the soil, rather, SOC helps plants to uptake nutrients in soil by altering physical, chemical and biological conditions [35]. However, the mineralization of SOM may release nutrients into the soil solution for the plants to uptake; hence, it is an important factor to assess soil fertility. Moreover, SOM is of significant importance for soil structure, aggregate stability, and clay dispersibility [36]. However, the cultivation of crops may disturb larger aggregates, resulting in enhanced SOM degradation. The plantation of local plants might have resulted in the loss of organic matter contents, which can deteriorate soil health and structure. Therefore, SOM needs to be maintained for the sustainable growth and development of the planted plants, so they can be used by honeybees. A recent study demonstrated that the application of SOM enhanced the visitation of bumblebee, as compared to without SOM application [37]. In another study, it was observed that higher SOM attracted 185 pollinators, as compared to 132 pollinators attracted by the soils with lower SOM [38]. This could possibly be due to the release of micro- and macronutrients from SOM mineralization, which improve the quality of flowers via increasing sugar contents and flower longevity, consequently attracting the pollinators towards them.
Nevertheless, due to the fact that SOC is sensitive to various factors, such as agricultural practices, climate, minerals in the soil, topography, and microbial communities, there is a higher spatial and temporal variability in its contents. This makes it difficult to compare the SOC of different sites, despite its extreme importance for soil quality topography [39]. Overall, SK sites have shown the better results of land restoration practices, as compared to QRC and AJ. The retention of SOC in the SK area, as compared to AJ and QRC, could be due to the difference in various factors, such as altitude, slope gradient, and pH. It has been reported that altitude has a significant effect on SOC retention in [40]. SK was a hilly area with a relatively higher altitude, thus, the development of various microhabitats might have facilitated the growth of plant communities, subsequently accumulating the higher C contents [41]. Additionally, the chances of aboveground biomass removal are more at lower altitudes than at higher altitudes. Thus, besides other factors, the higher altitude of SK facilitated the retention of SOC in that area after restoration activities. Therefore, loss of SOC in AJ and QRC during the restoration process in the current study reflects the poor SOC accumulation ability of the soil, and lower nutrient retention due to the sand based soil texture. Hence, complete restoration of these poor textured sites would require longer periods of time under sustainable practices.
The descriptive statistical analyses for the CaCO3 of the collected soil samples are shown in Table 3. The mean CaCO3 contents were in the range of 0.35% to 0.78%, suggesting the soils were not calcareous. It was noticed that CaCO3 content reduced from 0.62% in the control to 0.35% in restored areas in the AJ sites, whereas the CaCO3 contents increased from 0.52% and 0.39% in the control to 0.78% and 0.57% in restored areas in the QRC and SK sites, respectively. The lower CV values showed that there was little variability in the data. Likewise, the data was normal, as indicated by the values of skewness (0.533–1.777) and kurtosis (−0.910–2.279). The variations in CaCO3 content with soil depth, as shown in Figure 3, demonstrated different trends. The CaCO3 content increased with an increase in soil depth in AJ-C, while decreasing in QRC-C and SK-C. On contrary, the CaCO3 content decreased with soil depth in AJ-R, whereas it increased in the QRC-R and SK-R areas. As CaCO3 has a significant role in the stabilization of soil aggregates and SOC, the higher CaCO3 contents in the QRC-R and SK-R sites showed the positive impacts of restoration activities. Our results are in agreement with those reported by Rowley et al. [42]. However, the elevated CaCO3 contents in soils may increase the pH, subsequently hindering plant growth.

3.2. Soil Nutrient Contents as Indicators for Degradation and Restoration

3.2.1. Phosphorus and Potassium

The descriptive analyses for available P and K in the collected soil samples are shown in Table 4. The mean P concentration in all the samples ranged from 1.09 to 2.57 mg kg−1, suggesting poor soil fertility. The mean P concentrations in AJ and QRC are comparable, whereas, the mean P concentrations were lower in the SK sites. There was a huge variability in the data, as suggested by the larger CV values (72.55–112.36). However, the values of the skewness (−0.115–1.788) and kurtosis (−1.560–2.917) indicated normality in the data. The variations in mean P concentrations due to restoration activities were inconsistent. It was noticed that the mean P concentrations reduced from 2.27 mg kg−1 in the control to 1.21 mg kg−1 in restored areas in AJ. However, P concentrations increased from 1.09 mg kg−1 in the control to 2.57 mg kg−1 in restored areas in QRC, whereas a very slight decrease (0.02 units) was observed in the case of SK, due to the restoration process. The overall decline in available P in AJ and SK, due to restoration, could be due to the conversion of soil texture to sandy loam and variations in vegetation composition [43]. Whereas, the increased P concentrations in QRC after restoration was probably due to the mineralization of soil nutrients, variations in microbial composition and activities [44,45], as well as an uptake in nutrient uptake by plants [46]. Soil P is affected by the parent material, land-use type, and other biological as well as chemical processes, therefore, variations in these properties could result in different P retention in soil [47]. The variations in available P concentrations with soil depth are represented in Figure 5A,C,F. The highest P concentration was found in the 0–20 cm deep soil layer, as compared to the topsoil and deeper soil layers in AJ-C, AJ-R, and SK-C, while the other soil samples showed a decrease in P concentration with soil depth up to 40 cm, and then increases up to 60 cm. Similar results were also reported by Xu et al. [47], stating that >60% of the total P was detected in the upper soil layer. Moreover, the surface soil is largely influenced by climatic factors, microbial activities, and plant litter accumulation, which result in higher nutrient retention [48].
The descriptive analysis for available soil K showed a larger variability in the data (CV = 52.96–104.21), with a range of 33.23 mg kg−1 to 81.26 mg kg−1. It can be seen that the restoration activities resulted in an increment in available K concentrations. The K concertation was 56.36 mg kg−1 in the control of the AJ site, which increased to 75.58 mg kg−1 in restored areas. Likewise, it increased from 33.23 mg kg−1 in the control to 81.26 mg kg−1 in restored areas in the QRC sites, while no significant change was observed in the SK sites after restoration. The improved K in soil in all the restored areas could be due to the changes in soil structure and addition of fertilizers [39]. The variations in mean K concentrations with soil depth are shown in Figure 5B,D,F. The mean K concentrations did not follow any specific trend. For example, K concentration increased to the highest value (85.33 mg kg−1) in the 0–20 cm layer, decreased to 24.46 mg kg−1 in the 20–40 cm layer, and then increased to 54.13 mg kg−1 in the 40–60 cm layer, in AJ-C. A similar trend was seen in AJ-R and SK-C, while an opposite trend was observed in QRC-C and QRC-R. Therefore, these variations in P and K concentrations in soil with restoration activities could be due to the accumulation of organic compounds and the cultivation of plants [49]. However, in depth data is required to investigate the long term effects of the restoration activities.

3.2.2. Soil Micronutrients

The descriptive statistical analyses for the available soil micronutrients (Cu, Fe, Zn, and Mn) are shown in Table 5. Overall, the concentration of Fe was higher in all the soil types, followed by Mn, Cu, and Zn. The estimated CV values indicated a higher variability in the concentration of all the micronutrients, whereas the values of skewness and kurtosis suggested normality in the data. The effects of the restoration activities on micronutrients concentrations were not clear.
The concentration of Cu increased from 1.61 to 1.84 mg kg−1 and 1.32 to 1.49 mg kg−1 after restoration in the AJ and SK sites, respectively, while it reduced from 1.20 to 0.97 mg kg−1 after restoration in the QRC sites. In the case of Fe, a reduction from 5.04 to 4.04 mg kg−1 and 5.79 to 2.64 mg kg−1 was observed after restoration activities in the AJ and QRC sites, respectively, whereas an increment from 7.11 to 7.44 mg kg−1 was seen in the SK sites with restoration. Similarly, the concentration of Mn reduced from 2.82 to 1.80 mg kg−1 and 2.82 to 1.07 mg kg−1 in the restored areas of AJ and QRC, respectively, while, on the contrary, it increased from 2.13 to 3.31 mg kg−1 in the restored areas of SK. This trend was even different in the case of Zn, which reduced from 0.25 to 0.14 mg kg−1 after restoration in the AJ sites, while increasing from 0.08 to 0.12 mg kg−1 and 0.13 to 0.23 mg kg−1 in the QRC and SK sites, respectively, after the restoration. Therefore, these results reflect that there was no specific trend for the variations in the concentrations of micronutrients after land restoration. It is established that the availability of micronutrients is influenced by soil physiochemical and biological properties, such as speciation of the element, concentration of oxides and carbonates, SOM, and pH [50]. The availability of cationic micronutrients decreases, while the availability of anionic micronutrients increases, with a rise in pH. In the current study, the pH of the AJ and QRC sites increased with restoration, while it decreased in the SK sites. This decrease in soil pH might have resulted in the higher solubility and availability of the micronutrients (Cu, Fe, Mn, and Zn). On the contrary, the rise in pH has subsequently reduced the available portion of Fe, Mn, and Zn in the QRC and AJ sites. Moreover, the variations in the concentrations of micronutrients in soil after restoration could also be linked to an unstable soil structure, the uptake of nutrients by the plants, and leaching [51,52]. However, plant–soil interactions should be explored to understand the sustainable restoration of soil and higher availability of floral resources for honeybees.

3.3. Correlation Studies

The correlations between different physiochemical properties (pH, EC, SOC, SOM, and CaCO3) and soil nutrients (P, K, Cu, Fe, Mn, and Zn) were analyzed via Pearson’s correlation coefficients (Table 6). The values of correlation coefficient (r2) closer to 1 represent a stronger correlation, while the + or − sign shows the type of relationship at a 0.05 level of significance. The calculated correlations suggested a negative correlation of soil pH with EC (p < 0.05) and Fe contents. Moreover, soil Fe, Mn, and Zn contents were positively correlated with EC (p < 0.05). Likewise, soil Fe contents were positively correlated with SOC and SOM (p < 0.05). Nan et al. [53] reported the significant positive correlation of micronutrients with SOM. Soil Mn contents were also positively correlated with P and Fe (p < 0.05), while Zn was positively correlated with Fe and Mn (p < 0.05). Pearson correlations is a bivariant method used to investigate the relation between two parameters, thus, the results of the calculated correlation matrix showed some relations between soil available nutrients and soil properties [54]. Soil is a complex system made up of solid, liquid and gaseous phases, as well as various organisms. Therefore, the interpretation of different soil properties and of different types of soils is difficult due to the diversity of characteristics [55]. Therefore, due to varying soil locations, soil characteristics, and nutrient status, the correlations between the various parameters of soil collected from different localities were inconsistent.

4. Conclusions

Soil degradation extent in forests of southwestern Saudi Arabia was evaluated in this study. About 78 samples from Al-Jenebeen (AJ), Queen Rearing Center (QRC), and Shekiran (SK) areas were collected, both from control (no treatment) and restored (restoration activities) areas. Soil physiochemical properties, such as phosphorus (P), potassium (K), pH, electrical conductivity (EC), soil organic matter (SOM), soil organic carbon (SOC), calcite (CaCO3), texture, and micronutrients (Cu, Zn, Fe, and Mn), were determined. The data of soil physical and chemical properties showed inconsistent trends and larger variations. Overall, all the soil had lower SOM and nutrient contents, however, the soils were slightly alkaline and nonsaline. The results exhibited that sandy loam was the most dominant textural class in all the studied sites (~60%), followed by sandy clay loam. The variations in soil properties with restoration activities were inconstant. SK soil exhibited relatively prominent positive effects of restoration practices, as compared to QRC and AJ soils, as predicted by the highest SOC retention, reduced pH, and retention of nutrients. Therefore, currently, restored areas in SK could be beneficial for honeybee production due to the higher availability of floral resources. However, identifying the existing plant resources, soil physiochemical characteristics, and restoration potential of the study areas must be explored more for sustainable honeybee flourishing. Therefore, it was concluded that further temporal studies are needed to explore the long term effects of restoration activities on rehabilitation potential in the Al-baha region of Saudi Arabia.

Author Contributions

Conceptualization, project administration, resources and methodology, A.A.A.-G.; investigation and review and editing, Y.T.; software and visualization, N.A.; validation, formal analysis, writing and review and editing, A.G.A. All authors have read and agreed to the published version of the manuscript.

Funding

This project was funded by the National Plan for Science, Technology and Innovation (MAARIFAH), King Abdulaziz City for Science and Technology, Kingdom of Saudi Arabia, Award Number (13-AGR2118-02).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors extend their appreciation to the National Plan for Science, Technology and Innovation (MAARIFAH), King Abdulaziz City for Science and Technology, Kingdom of Saudi Arabia, Award Number (13-AGR2118-02).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mesfin, S.; Taye, G.; Desta, Y.; Sibhatu, B.; Muruts, H.; Mohammedbrhan, M. Short-Term Effects of Bench Terraces on Selected Soil Physical and Chemical Properties: Landscape Improvement for Hillside Farming in Semi-Arid Areas of Northern Ethiopia. Environ. Earth. Sci. 2018, 77, 399. [Google Scholar] [CrossRef]
  2. Alves, L.F.; Vieira, S.A.; Scaranello, M.A.; Camargo, P.B.; Santos, F.A.M.; Joly, C.A.; Martinelli, L.A. Forest Structure and Live Aboveground Biomass Variation along an Elevational Gradient of Tropical Atlantic Moist Forest (Brazil). For. Ecol. Manag. 2010, 260, 679–691. [Google Scholar] [CrossRef]
  3. Lamb, D.; Erskine, P.D.; Parrotta, J.A. Restoration of Degraded Tropical Forest Landscapes. Science 2005, 310, 1628–1632. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. León, J.D.; Osorio, N.W. Role of Litter Turnover in Soil Quality in Tropical Degraded Lands of Colombia. Sci. World J. 2014, 2014, e693981. [Google Scholar] [CrossRef]
  5. Lal, R. Restoring Soil Quality to Mitigate Soil Degradation. Sustainability 2015, 7, 5875–5895. [Google Scholar] [CrossRef] [Green Version]
  6. Martínez, F.J.; Sierra, M.; Sierra, C.; Roca, A. Assessing Sustainable Use of Land under Olive Cultivation in Alcala La Real (Jaen, Spain) Using GIS. Adv. Geoecol. 2005, 36, 75–84. [Google Scholar]
  7. Khresat, S.A.; Qudah, E.A. Formation and Properties of Aridic Soils of Azraq Basin in Northeastern Jordan. J. Arid Environ. 2006, 64, 116–136. [Google Scholar] [CrossRef]
  8. Mbagwu, J.S.C. Aggregate Stability and Soil Degradation in the Tropics; ETDEWEB: Trieste, Italy, 2004; pp. 246–252. Available online: http://users.ictp.it/~pub_off/lectures/lns018/22Mbagwu1.pdf (accessed on 10 December 2020).
  9. Kremen, C.; Williams, N.M.; Thorp, R.W. Crop Pollination from Native Bees at Risk from Agricultural Intensification. Proc. Natl. Acad. Sci. USA 2002, 99, 16812–16816. [Google Scholar] [CrossRef] [Green Version]
  10. Isaacs, R.; Williams, N.; Ellis, J.; Pitts-Singer, T.L.; Bommarco, R.; Vaughan, M. Integrated Crop Pollination: Combining Strategies to Ensure Stable and Sustainable Yields of Pollination-Dependent Crops. Basic Appl. Ecol. 2017, 22, 44–60. [Google Scholar] [CrossRef]
  11. Koh, I.; Lonsdorf, E.V.; Williams, N.M.; Brittain, C.; Isaacs, R.; Gibbs, J.; Ricketts, T.H. Modeling the Status, Trends, and Impacts of Wild Bee Abundance in the United States. Proc. Natl. Acad. Sci. USA 2016, 113, 140–145. [Google Scholar] [CrossRef] [Green Version]
  12. Biesmeijer, J.C.; Roberts, S.P.M.; Reemer, M.; Ohlemüller, R.; Edwards, M.; Peeters, T.; Schaffers, A.P.; Potts, S.G.; Kleukers, R.; Thomas, C.D.; et al. Parallel Declines in Pollinators and Insect-Pollinated Plants in Britain and the Netherlands. Science 2006, 313, 351–354. [Google Scholar] [CrossRef]
  13. Kovács-Hostyánszki, A.; Espíndola, A.; Vanbergen, A.J.; Settele, J.; Kremen, C.; Dicks, L.V. Ecological Intensification to Mitigate Impacts of Conventional Intensive Land Use on Pollinators and Pollination. Ecol. Lett. 2017, 20, 673–689. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Al-Ghamdi, A.A.; Tadesse, Y.; Adgaba, N. Evaluation of Major Acacia Species in the Nursery towards Apicultural Landscape Restoration around Southwestern Saudi Arabia. Saudi J. Biol. Sci. 2020, 27, 3385–3389. [Google Scholar] [CrossRef]
  15. Ryan, J.; Estefan, G.; Rashid, A. Soil and Plant Analysis Laboratory Manual; ICARDA—Jointly published by International Center for Agricultural Research in Dry areas, Syria and The National Agriculture Research Center: Islamabad, Pakistan, 2001; ISBN 978-92-9127-118-4.
  16. Koehler, F.E.; Moodie, C.D.; McNeal, B.L. Laboratory Manual for Soil Fertility; Washington State University: Pulman, WA, USA, 1984. [Google Scholar]
  17. Horváth, B.; Opara-Nadi, O.; Beese, F. A Simple Method for Measuring the Carbonate Content of Soils. Soil Sci. Soc. Am. J. 2005, 69, 1066–1068. [Google Scholar] [CrossRef]
  18. Walkley, A.; Black, I.A. An examination of the degtjareff method for determining soil organic matter, and a proposed modification of the chromic acid titration method. Soil Sci. 1934, 37, 29–38. [Google Scholar] [CrossRef]
  19. Soltanpour, P.N.; Schwab, A.P. A New Soil Test for Simultaneous Extraction of Macro- and Micro-nutrients in Alkaline Soils. Commun. Soil Sci. Plant Anal. 1977, 8, 195–207. [Google Scholar] [CrossRef]
  20. Pattison, A.B.; Moody, P.W.; Badcock, K.A.; Smith, L.J.; Armour, J.A.; Rasiah, V.; Cobon, J.A.; Gulino, L.-M.; Mayer, R. Development of Key Soil Health Indicators for the Australian Banana Industry. Appl. Soil Ecol. 2008, 40, 155–164. [Google Scholar] [CrossRef]
  21. Khadka, D.; Lamichhane, S.; Shrestha, K.; Joshi, S.; Karna, M.; Pant, B.B.; Yadav, S. Soil Fertility Assessment and Mapping of Agricultural Research Station, Jaubari, Illam, Nepal. Int. J. Environ. 2017, 6, 46–70. [Google Scholar] [CrossRef] [Green Version]
  22. Tadesse, B.; Mesfin, S.; Tesfay, G.; Abay, F. Effect of Integrated Soil Bunds on Key Soil Properties and Soil Carbon Stock in Semi-Arid Areas of Northern Ethiopia. S. Afr. J. Plant Soil 2016, 33, 297–302. [Google Scholar] [CrossRef]
  23. Wang, X.P.; Li, X.R.; Xiao, H.L.; Pan, Y.X. Evolutionary Characteristics of the Artificially Revegetated Shrub Ecosystem in the Tengger Desert, Northern China. Ecol. Res. 2006, 21, 415–424. [Google Scholar] [CrossRef]
  24. Hartemink, A.E. Soil Chemical and Physical Properties as Indicators of Sustainable Land Management under Sugar Cane in Papua New Guinea. Geoderma 1998, 85, 283–306. [Google Scholar] [CrossRef]
  25. Shukla, A.; Vyas, D.; Anuradha, J. Soil Depth: An Overriding Factor for Distribution of Arbuscular Mycorrhizal Fungi. J. Soil Sci. Plant Nutr. 2013, 13, 23–33. [Google Scholar] [CrossRef]
  26. Zhang, Y.L.; Chen, L.J.; Chen, X.H.; Tan, M.L.; Duan, Z.H.; Wu, Z.J.; Li, X.J.; Fan, X.H. Response of Soil Enzyme Activity to Long-Term Restoration of Desertified Land. Catena 2015, 133, 64–70. [Google Scholar] [CrossRef]
  27. Scianna, J. Salt-Affected Soils: Their Causes, Measure, and Classification. Res. Method Hort. Note 2002, 5, 1–3. [Google Scholar]
  28. He, B.; Cai, Y.; Ran, W.; Jiang, H. Spatial and Seasonal Variations of Soil Salinity Following Vegetation Restoration in Coastal Saline Land in Eastern China. Catena 2014, 118, 147–153. [Google Scholar] [CrossRef]
  29. Bezborodov, G.A.; Shadmanov, D.K.; Mirhashimov, R.T.; Yuldashev, T.; Qureshi, A.S.; Noble, A.D.; Qadir, M. Mulching and Water Quality Effects on Soil Salinity and Sodicity Dynamics and Cotton Productivity in Central Asia. Agric. Ecosyst. Environ. 2010, 138, 95–102. [Google Scholar] [CrossRef]
  30. Krull, E.S.; Skjemstad, J.O.; Baldock, J.A. Functions of Soil Organic Matter and the Effect on Soil Properties; Cooperative Research Centre for Greenhouse Accounting: Canberra, Australia, 2004. [Google Scholar]
  31. Obalum, S.E.; Okpara, I.M.; Obi, M.E.; Wakatsuki, T. Short Term Effects of Tillage-Mulch Practices under Sorghum and Soybean on Organic Carbon and Eutrophic Status of a Degraded Ultisol in Southeastern Nigeria. Trop. Subtrop. Agroecosyst. 2011, 14, 393–403. [Google Scholar]
  32. Spaccini, R.; Piccolo, A.; Mbagwu, J.S.C.; Teshale, A.Z.; Igwe, C.A. Influence of the Addition of Organic Residues on Carbohydrate Content and Structural Stability of Some Highland Soils in Ethiopia. Soil Use Manag. 2002, 18, 404–411. [Google Scholar] [CrossRef]
  33. du Preez, C.C.; van Huyssteen, C.W.; Mnkeni, P.N.S. Land Use and Soil Organic Matter in South Africa 2: A Review on the Influence of Arable Crop Production. S. Afr. J. Sci. 2011, 107, 35–42. [Google Scholar] [CrossRef] [Green Version]
  34. Kalinina, O.; Cherkinsky, A.; Chertov, O.; Goryachkin, S.; Kurganova, I.; Lopes de Gerenyu, V.; Lyuri, D.; Kuzyakov, Y.; Giani, L. Post-Agricultural Restoration: Implications for Dynamics of Soil Organic Matter Pools. Catena 2019, 181, 104096. [Google Scholar] [CrossRef]
  35. Hoyle, F.C.; Baldock, J.A.; Murphy, D.V. Soil Organic Carbon—Role in Rainfed Farming Systems. In Rainfed Farming Systems; Tow, P., Cooper, I., Partridge, I., Birch, C., Eds.; Springer: Dordrecht, The Netherlands, 2011; pp. 339–361. ISBN 978-1-4020-9132-2. [Google Scholar]
  36. Dexter, A.R.; Pagliai, M.; Jones, R. Soil Structure: The Key to Soil Function. Adv. GeoEcol. 2002, 35, 57–69. [Google Scholar]
  37. Banaszak-Cibicka, W.; Takacs, V.; Kesy, M.; Langowska, A.; Blecharczyk, A.; Sawinska, Z.; Sparks, T.H.; Tryjanowski, P. Manure application improves both bumblebee flower visitation and crop yield in intensive farmland. Basic Appl. Ecol. 2019, 36, 26–33. [Google Scholar] [CrossRef]
  38. Sáez, A.; Morales, C.L.; Ramos, L.Y.; Aizen, M.A. Extremely frequent bee visits increase pollen deposition but reduce drupelet set in raspberry. J. Appl. Ecol. 2014, 51, 1603–1612. [Google Scholar] [CrossRef]
  39. Baldock, J.; Skjemstad, J.O. Soil Organic Carbon/Soil Organic Matter; CSIRO Publishing: Collingwood, VIC, Australia, 1999; ISBN 978-0-643-06376-1. [Google Scholar]
  40. Tashi, S.; Singh, B.; Keitel, C.; Adams, M. Soil carbon and nitrogen stocks in forests along an altitudinal gradient in the eastern Himalayas and a meta-analysis of global data. Glob. Chang. Biol. 2016, 22, 2255–2268. [Google Scholar] [CrossRef] [PubMed]
  41. do Carmo, F.F.; Jacobi, C.M. Diversity and plant trait-soil relationships among rock outcrops in the Brazilian Atlantic rainforest. Plant Soil 2016, 403, 7–20. [Google Scholar] [CrossRef]
  42. Rowley, M.C.; Grand, S.; Verrecchia, É.P. Calcium-Mediated Stabilisation of Soil Organic Carbon. Biogeochemistry 2018, 137, 27–49. [Google Scholar] [CrossRef] [Green Version]
  43. Dong, S.K.; Wen, L.; Li, Y.Y.; Wang, X.X.; Zhu, L.; Li, X.Y. Soil-Quality Effects of Grassland Degradation and Restoration on the Qinghai-Tibetan Plateau. Soil Sci. Soc. Am. J. 2012, 76, 2256–2264. [Google Scholar] [CrossRef]
  44. Chen, D.D.; Zhang, S.H.; Dong, S.K.; Wang, X.T.; Du, G.Z. Effect of Land-Use on Soil Nutrients and Microbial Biomass of an Alpine Region on the Northeastern Tibetan Plateau, China. Land Degrad. Dev. 2010, 21, 446–452. [Google Scholar] [CrossRef]
  45. Feng, R.; Long, R.; Shang, Z.; Ma, Y.; Dong, S.; Wang, Y. Establishment of Elymus Natans Improves Soil Quality of a Heavily Degraded Alpine Meadow in Qinghai-Tibetan Plateau, China. Plant Soil 2010, 327, 403–411. [Google Scholar] [CrossRef]
  46. Wu, G.L.; Liu, Z.H.; Zhang, L.; Hu, T.M.; Chen, J.M. Effects of Artificial Grassland Establishment on Soil Nutrients and Carbon Properties in a Black-Soil-Type Degraded Grassland. Plant Soil 2010, 333, 469–479. [Google Scholar] [CrossRef]
  47. Xu, H.; Qu, Q.; Li, P.; Guo, Z.; Wulan, E.; Xue, S. Stocks and Stoichiometry of Soil Organic Carbon, Total Nitrogen, and Total Phosphorus after Vegetation Restoration in the Loess Hilly Region, China. Forests 2019, 10, 27. [Google Scholar] [CrossRef] [Green Version]
  48. Barreto, P.A.B.; Gama-Rodrigues, E.F.; Gama-Rodrigues, A.C.; Fontes, A.G.; Polidoro, J.C.; Moço, M.K.S.; Machado, R.C.R.; Baligar, V.C. Distribution of Oxidizable Organic C Fractions in Soils under Cacao Agroforestry Systems in Southern Bahia, Brazil. Agrofor. Syst. 2011, 81, 213–220. [Google Scholar] [CrossRef]
  49. Ruiz-Sinoga, J.D.; Diaz, A.R. Soil Degradation Factors along a Mediterranean Pluviometric Gradient in Southern Spain. Geomorphology 2010, 118, 359–368. [Google Scholar] [CrossRef]
  50. Seenivasan, R.; Prasath, V.; Mohanraj, R. Restoration of Sodic Soils Involving Chemical and Biological Amendments and Phytoremediation by Eucalyptus Camaldulensis in a Semiarid Region. Environ. Geochem. Health 2015, 37, 575–586. [Google Scholar] [CrossRef]
  51. Nicholson, F.A.; Smith, S.R.; Alloway, B.J.; Carlton-Smith, C.; Chambers, B.J. Quantifying Heavy Metal Inputs to Agricultural Soils in England and Wales. Water Environ. J. 2006, 20, 87–95. [Google Scholar] [CrossRef]
  52. Qadir, M.; Tubeileh, A.; Akhtar, J.; Larbi, A.; Minhas, P.S.; Khan, M.A. Productivity Enhancement of Salt-Affected Environments through Crop Diversification. Land Degrad. Dev. 2008, 19, 429–453. [Google Scholar] [CrossRef]
  53. Nan, Z.; Zhao, C.; Li, J.; Chen, F.; Sun, W. Relations Between Soil Properties and Selected Heavy Metal Concentrations in Spring Wheat (Triticum Aestivum L.) Grown in Contaminated Soils. Water Air Soil Pollut. 2002, 133, 205–213. [Google Scholar] [CrossRef]
  54. Sharma, M.R.; Raju, N.S. Correlation of Heavy Metal Contamination with Soil Properties of Industrial Areas of Mysore, Karnataka, India by Cluster Analysis. Int. J. Environ. Sci. 2013, 2, 22–27. [Google Scholar]
  55. Ludovisi, A.; Minozzo, M.; Pandolfi, P.; Taticchi, M.I. Modelling the Horizontal Spatial Structure of Planktonic Community in Lake Trasimeno (Umbria, Italy) Using Multivariate Geostatistical Methods. Ecol. Modell. 2005, 181, 247–262. [Google Scholar] [CrossRef]
Figure 1. Proportion of various textural classes in the soil samples collected from different localities in southwestern Saudi Arabia (QRC-C: Queen Rearing Center Control, QRC-R: Queen Rearing Center Restored, SK-C: Shekiran control, SK-R: Shekiran Restored, AJ-C: Al-Jenebeen Control, and AJ-R: Al-Jenebeen Restored).
Figure 1. Proportion of various textural classes in the soil samples collected from different localities in southwestern Saudi Arabia (QRC-C: Queen Rearing Center Control, QRC-R: Queen Rearing Center Restored, SK-C: Shekiran control, SK-R: Shekiran Restored, AJ-C: Al-Jenebeen Control, and AJ-R: Al-Jenebeen Restored).
Sustainability 13 09169 g001
Figure 2. Variations in mean pH and electrical conductivity (EC) with soil depth for the samples collected from different localities in southwestern Saudi Arabia ((A): pH and (B): EC in Al–Jenebeen, (C): pH and (D): EC in Queen Rearing Center, and (E): pH and (F): EC in Shekiran). The error bars represent the standard deviations of the mean values.
Figure 2. Variations in mean pH and electrical conductivity (EC) with soil depth for the samples collected from different localities in southwestern Saudi Arabia ((A): pH and (B): EC in Al–Jenebeen, (C): pH and (D): EC in Queen Rearing Center, and (E): pH and (F): EC in Shekiran). The error bars represent the standard deviations of the mean values.
Sustainability 13 09169 g002
Figure 3. Variations in mean soil organic matter (SOM) contents and CaCO3 with soil depth for the samples collected from different localities in southwestern Saudi Arabia ((A): SOM and (B): CaCO3 in Al–Jenebeen, (C): SOM and (D): CaCO3 in Queen Rearing Center, and (E): SOM and (F): CaCO3 in Shekiran). The error bars represent the standard deviations of the mean values.
Figure 3. Variations in mean soil organic matter (SOM) contents and CaCO3 with soil depth for the samples collected from different localities in southwestern Saudi Arabia ((A): SOM and (B): CaCO3 in Al–Jenebeen, (C): SOM and (D): CaCO3 in Queen Rearing Center, and (E): SOM and (F): CaCO3 in Shekiran). The error bars represent the standard deviations of the mean values.
Sustainability 13 09169 g003
Figure 4. Proportions of soil organic carbon according to the soil quality evaluation criteria for the samples collected from different localities in southwestern Saudi Arabia (QRC-C: Queen Rearing Center Control, QRC-R: Queen Rearing Center Restored, SK-C: Shekiran control, SK-R: Shekiran Restored, AJ-C: Al-Jenebeen Control, and AJ-R: Al-Jenebeen Restored). The used criteria are: extremely low = <0.4%, moderately low = 0.4–0.59%, low = 0.6–0.99%, moderate = 1–1.59%, high = 1.6–1.99%, very high = 2–2.99%, and extremely high = >3% (Hazelton and Murphy, 2007).
Figure 4. Proportions of soil organic carbon according to the soil quality evaluation criteria for the samples collected from different localities in southwestern Saudi Arabia (QRC-C: Queen Rearing Center Control, QRC-R: Queen Rearing Center Restored, SK-C: Shekiran control, SK-R: Shekiran Restored, AJ-C: Al-Jenebeen Control, and AJ-R: Al-Jenebeen Restored). The used criteria are: extremely low = <0.4%, moderately low = 0.4–0.59%, low = 0.6–0.99%, moderate = 1–1.59%, high = 1.6–1.99%, very high = 2–2.99%, and extremely high = >3% (Hazelton and Murphy, 2007).
Sustainability 13 09169 g004
Figure 5. Variations in mean phosphorus and potassium concentrations with soil depth for the samples collected from different localities in southwestern Saudi Arabia ((A): concentration of phosphorus and (B): concentration of potassium in Al–Jenebeen, (C): concentration of phosphorus and (D): concentration of potassium in Queen Rearing Center, and (E): concentration of phosphorus and (F): concentration of potassium in Shekiran). The error bars represent the standard deviations of the mean values.
Figure 5. Variations in mean phosphorus and potassium concentrations with soil depth for the samples collected from different localities in southwestern Saudi Arabia ((A): concentration of phosphorus and (B): concentration of potassium in Al–Jenebeen, (C): concentration of phosphorus and (D): concentration of potassium in Queen Rearing Center, and (E): concentration of phosphorus and (F): concentration of potassium in Shekiran). The error bars represent the standard deviations of the mean values.
Sustainability 13 09169 g005
Table 1. Details of the collected samples along with the coordinates of the sampling locations.
Table 1. Details of the collected samples along with the coordinates of the sampling locations.
No.AreaCodeLand UseTypeNo. of SamplesLatitudeLongitude
1Al-JenebeenAJ-RRestoredProfile-1319.8604541.69706
2Al-JenebeenAJ-RRestoredProfile-2319.8600941.696915
3Al-JenebeenAJ-RRestoredProfile-3319.8600241.696861
4Al-JenebeenAJ-RRestoredComposite and surface4----
5Al-JenebeenAJ-CControlProfile-1319.8599941.69712
6Al-JenebeenAJ-CControlProfile-2319.8603741.69768
7Al-JenebeenAJ-CControlProfile-3319.8583041.68455
8Al-JenebeenAJ-CControlComposite and surface4----
9Queen Rearing CenterQRC-RRestoredProfile-1319.8486441.58691
10Queen Rearing CenterQRC-RRestoredProfile-2319.8480441.58657
11Queen Rearing CenterQRC-RRestoredProfile-3319.8482341.58625
12Queen Rearing CenterQRC-RRestoredComposite and surface4----
13Queen Rearing CenterQRC-CControlProfile-1319.8475541.58692
14Queen Rearing CenterQRC-CControlProfile-2319.8480841.58704
15Queen Rearing CenterQRC-CControlProfile-3319.8485441.58651
16Queen Rearing CenterQRC-CControlComposite and surface 4----
17ShekiranSK-RRestoredProfile-1319.8503141.58733
18ShekiranSK-RRestoredProfile-2319.8503741.58776
19ShekiranSK-RRestoredProfile-3319.8503241.58773
20ShekiranSK-RRestoredComposite and surface4----
21ShekiranSK-CControlProfile-1319.8488941.58808
22ShekiranSK-CControlProfile-2319.8491641.58823
23ShekiranSK-CControlProfile-3319.8495541.58846
24ShekiranSK-CControlComposite and surface 4----
Table 2. Descriptive statistics of the pH and EC data of soil samples collected from Al-Jenebeen (AJ), Queen Rearing Center (QRC), and Shekiran (SK) in southwestern Saudi Arabia.
Table 2. Descriptive statistics of the pH and EC data of soil samples collected from Al-Jenebeen (AJ), Queen Rearing Center (QRC), and Shekiran (SK) in southwestern Saudi Arabia.
ParameterLocationLand UseMeanMedianMin.Max.SDCVSkewKurtosis
pHAJControl (n = 13)8.078.037.828.271.161.99−0.092−1.438
Restored (n = 13)8.158.247.878.310.672.07−0.993−0.718
QRCControl (n = 13)7.747.707.288.250.293.750.406−0.622
Restored (n = 13)7.998.027.608.290.212.63−0.619−0.646
SKControl (n = 13)8.017.967.818.360.182.230.914−0.303
Restored (n = 13)7.817.827.398.070.22.51−0.678−0.041
EC (dS m−1)AJControl (n = 13)0.450.390.241.180.2556.572.0153.602
Restored (n = 13)0.360.300.201.090.2464.642.4505.181
QRCControl (n = 13)0.430.480.200.600.1431.25−0.503−1.103
Restored (n = 13)0.380.370.190.550.1128.34−0.013−0.971
SKControl (n = 13)0.340.310.240.560.0926.971.0260.666
Restored (n = 13)0.400.360.240.600.1332.460.305−1.423
Table 3. Descriptive statistics for soil organic matter (SOM) and CaCO3 of the soil samples collected from Al-Jenebeen (AJ), Queen Rearing Center (QRC), and Shekiran (SK) in southwestern Saudi Arabia.
Table 3. Descriptive statistics for soil organic matter (SOM) and CaCO3 of the soil samples collected from Al-Jenebeen (AJ), Queen Rearing Center (QRC), and Shekiran (SK) in southwestern Saudi Arabia.
ParameterLocationLand UseMeanMedianMin.Max.SDCVSkewKurtosis
SOM (%)AJControl (n = 13)1.271.130.252.190.5543.470.183−0.304
Restored (n = 13)1.020.960.561.460.2827.240.242−1.007
QRCControl (n = 13)0.880.630.131.890.6068.400.244−1.326
Restored (n = 13)0.800.830.261.320.3037.25−0.214−0.556
SKControl (n = 13)1.161.210.531.780.3933.99−0.059−1.212
Restored (n = 13)1.401.530.132.560.6646.83−0.361−0.375
CaCO3 (%)AJControl (n = 13)0.620.480.521.520.0083.890.533−0.910
Restored (n = 13)0.350.300.250.960.0473.070.8860.460
QRCControl (n = 13)0.520.260.622.220.04120.551.7772.279
Restored (n = 13)0.780.350.792.390.00101.470.898−0.699
SKControl (n = 13)0.390.300.411.390.00105.511.1230.570
Restored (n = 13)0.570.260.701.960.00122.701.039−0.521
Table 4. Descriptive statistics for phosphorus (P) and potassium (K) concentrations of the soil samples collected from Al-Jenebeen (AJ), Queen Rearing Center (QRC), and Shekiran (SK) in southwestern Saudi Arabia.
Table 4. Descriptive statistics for phosphorus (P) and potassium (K) concentrations of the soil samples collected from Al-Jenebeen (AJ), Queen Rearing Center (QRC), and Shekiran (SK) in southwestern Saudi Arabia.
ParameterLocationLand UseMeanMedianMin.Max.SDCVSkewKurtosis
P (mg kg−1) AJControl (n = 13)2.271.650.049.462.55112.361.7882.917
Restored (n = 13)1.211.560.042.570.9376.50−0.115−1.560
QRCControl (n = 13)1.090.280.113.491.21111.120.907−0.726
Restored (n = 13)2.572.570.046.241.9776.780.172−0.916
SKControl (n = 13)1.821.650.043.491.3272.55−0.014−1.448
Restored (n = 13)1.801.190.046.241.7295.771.3031.347
K (mg kg−1) AJControl (n = 13)56.3631.0011.2021056.72100.631.6452.065
Restored (n = 13)75.5843.0012.0020867.4589.241.106−0.332
QRCControl (n = 13)33.2327.609.4070.2017.5952.960.623−0.410
Restored (n = 13)81.2632.408.20216.0084.68104.210.797−1.255
SKControl (n = 13)41.2132.4018.2098.8025.1961.131.3680.503
Restored (n = 13)41.1025.007.60137.6036.9989.991.6431.601
Table 5. Descriptive statistics for heavy metals (Cu, Fe, Mn, and Zn) concentrations of the soil samples collected from Al-Jenebeen (AJ), Queen Rearing Center (QRC), and Shekiran (SK) in southwestern Saudi Arabia.
Table 5. Descriptive statistics for heavy metals (Cu, Fe, Mn, and Zn) concentrations of the soil samples collected from Al-Jenebeen (AJ), Queen Rearing Center (QRC), and Shekiran (SK) in southwestern Saudi Arabia.
ParameterLocationLand useMeanMedianMin.Max.SDCVSkewKurtosis
Cu (mg kg−1)AJControl (n = 13)1.611.210.674.121.0565.251.4660.854
Restored (n = 13)1.841.550.514.211.1662.990.805−0.420
QRCControl (n = 13)1.201.160.312.640.6654.500.5830.058
Restored (n = 13)0.970.820.312.020.5860.150.732−0.665
SKControl (n = 13)1.321.480.362.460.7153.990.088−1.270
Restored (n = 13)1.491.640.382.290.5335.73−0.657−0.330
Fe (mg kg−1)AJControl(n = 13)5.044.530.5212.454.0179.610.797−0.646
Restored (n = 13)4.042.850.028.632.5964.150.185−1.185
QRCControl (n = 13)5.796.820.009.282.6145.12−0.804−0.211
Restored (n = 13)2.642.860.005.432.0176.060.0081.497
SKControl (n = 13)7.117.882.6713.182.9942.060.312−0.600
Restored (n = 13)7.447.153.7411.532.5634.460.265−1.085
Mn (mg kg−1)AJControl (n = 13)2.821.380.5610.213.27115.711.6031.000
Restored (n = 13)1.801.570.154.431.2770.470.631−0.502
QRCControl (n = 13)2.822.610.097.672.4085.130.618−0.776
Restored (n = 13)1.070.890.103.470.9286.011.4251.630
SKControl (n = 13)2.131.440.695.541.6577.480.993−0.406
Restored (n = 13)3.311.440.4815.644.28129.112.0843.403
Zn (mg kg−1)AJControl(n = 13)0.250.030.001.390.45176.191.8291.802
Restored (n = 13)0.140.000.000.640.20147.051.4160.947
QRCControl (n = 13)0.080.040.000.360.10135.601.5811.983
Restored (n = 13)0.120.000.001.090.32262.622.4314.495
SKControl (n = 13)0.130.030.000.580.18135.611.2790.787
Restored (n = 13)0.230.000.001.480.48206.451.9202.018
Table 6. Pearson correlation coefficients (r2) and probability (p) as estimated for the relation of different parameters (pH, electrical conductivity (EC), heavy metals (Cu, Fe, Mn, and Zn), soil organic carbon (SOC), soil organic matter (SOM), CaCO3, phosphorus (P), and potassium (K)) of soil samples collected from Al-Jenebeen (AJ), Queen Rearing Center (QRC), and Shekiran (SK) in southwestern Saudi Arabia.
Table 6. Pearson correlation coefficients (r2) and probability (p) as estimated for the relation of different parameters (pH, electrical conductivity (EC), heavy metals (Cu, Fe, Mn, and Zn), soil organic carbon (SOC), soil organic matter (SOM), CaCO3, phosphorus (P), and potassium (K)) of soil samples collected from Al-Jenebeen (AJ), Queen Rearing Center (QRC), and Shekiran (SK) in southwestern Saudi Arabia.
p DepthpHEC SOC SOM CaCO3PKCuFeMn
pHr2−0.128
p-value0.265
ECr2−0.119−0.295
p-value0.3010.009 *
SOCr2−0.068−0.148−0.005
p-value0.5530.1970.967
SOMr2−0.069−0.147−0.0031.000
p-value0.5470.1990.9790.000 *
CaCO3r20.1080.139−0.0060.1130.113
p-value0.3480.2260.9580.3250.324
Pr20.1860.047−0.0140.1020.101−0.017
p-value0.1020.6850.9020.3720.3790.881
Kr20.0700.209−0.130−0.032−0.033−0.1060.637
p-value0.5450.0670.2570.7830.7730.3560.000 *
Cur2−0.0230.037−0.041−0.001−0.0010.0190.1450.105
p-value0.8410.7460.7230.9910.9920.8690.2040.360
Fer20.073−0.4460.3920.2790.280−0.0200.003−0.1800.297
p-value0.5250.000 *0.000 *0.013 *0.013 *0.8650.9800.1150.008 *
Mnr20.137−0.1910.2930.1070.107−0.1160.3530.1660.0810.496
p-value0.2330.0950.009 *0.3510.3530.3110.002 *0.1470.4830.000 *
Znr2−0.069−0.0340.4140.0010.001−0.1550.1840.0610.0290.3920.638
p-value0.5510.7710.000 *0.9920.9910.1760.1070.5940.8000.000 *0.000 *
* Correlation is significant at the 0.05 level (two-tailed).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Al-Ghamdi, A.A.; Tadesse, Y.; Adgaba, N.; Alghamdi, A.G. Soil Degradation and Restoration in Southwestern Saudi Arabia through Investigation of Soil Physiochemical Characteristics and Nutrient Status as Indicators. Sustainability 2021, 13, 9169. https://doi.org/10.3390/su13169169

AMA Style

Al-Ghamdi AA, Tadesse Y, Adgaba N, Alghamdi AG. Soil Degradation and Restoration in Southwestern Saudi Arabia through Investigation of Soil Physiochemical Characteristics and Nutrient Status as Indicators. Sustainability. 2021; 13(16):9169. https://doi.org/10.3390/su13169169

Chicago/Turabian Style

Al-Ghamdi, Ahmad A., Yilma Tadesse, Nuru Adgaba, and Abdulaziz G. Alghamdi. 2021. "Soil Degradation and Restoration in Southwestern Saudi Arabia through Investigation of Soil Physiochemical Characteristics and Nutrient Status as Indicators" Sustainability 13, no. 16: 9169. https://doi.org/10.3390/su13169169

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