3.1 Soil physical and chemical properties
The predominant soil type is Alfisol, comprising sand (60.74, 58.95, and 57.42 %), clay (28.26, 27.93, and 27.25 %), and silt (11.00, 13.11, and 15.33 %), with sand being the most abundant, followed by clay and silt. At Ìrèle and Ìpaò, the soil pH was strongly acidic, while it was moderate for Òkè-Àkò. Below the critical value of 10 mg/kg, Ìrèle exhibited the lowest available phosphorus level, whereas Òkè-Àkò and Ìpaò showed significant increases. Nitrogen and potassium (exchangeable) levels exceeded the critical levels of 1 mg/kg and 0.2 cmol/kg, respectively, although calcium and magnesium levels were lower, with the highest values recorded in Òkè-Àkò, significantly differing from the other two locations. The cation exchange capacity (CEC) was lower than the recommended range of between 10-20 cmol/kg, and exchangeable acidity (EA) was higher than the critical limit of 1.0 cmol/kg. At all locations, the base saturation level was slightly above the critical limit of 50 % and there were low levels of trace elements including copper, zinc, iron as well as manganese. Our previous work [18], further present the detailed physical and chemical properties of soil samples taken from three locations.
3.2 Precipitation and Temperature of the research area.
The presentation of precipitation and temperature patterns in the study area is outlined by [18]. The cumulative rainfall, measured in millimeters, during the period from April to August 2018 at the research site is detailed as follows: 120 mm in April, 152 mm in May, 168 mm in June, 170 mm in July, and 131 mm in August. Concurrently, the average temperatures, expressed in degrees Celsius, were recorded as 27.4 °C in April, 26.9 °C in May, 26.2 °C in June, 25.8 °C in July, and 24.9 °C in August.
3.3 The Interaction effect of sites and location on specific phosphorus and sulfur enzymes, and soil nutrient status.
Significant interaction effects between site and location were observed for the activities of phosphorus and sulfur cycling enzymes (Table 1). The recorded interaction effects were significant and varied both in size and direction of their response. Notably, phosphatase (Pho) levels were significantly higher at NPS in Òkè-Àkò and Ìpaò than sulfur enzymes. However, the peak value (2.94 mg/ml/min) was observed at NPS in Òkè-Àkò. Thiosulfate dehydrogenase (Tsd) exhibited significant site-by-location interaction effects, with the highest value recorded at NPS in Ìrèle and the lowest activity level at NPS in Òkè-Àkò, although the variation was insignificant. The recorded interaction effects were both significant and diverse, involving variations in magnitude and direction of response. Additionally, Dimethyl sulfoxide reductase (Dsr) was found to be significantly higher at CPS in Ìpaò, while low activity levels of (0.80 and 0.68 µg/ml/min) were observed at NPS in Òkè-Àkò and Ìpaò, respectively. Lastly, phytase (Phy) activity remained consistent across both sites and depths but showed significant variations among means.
Significant interaction effects between site and location were observed regarding the soil nutrient status, as indicated in Table 1. Copper and cobalt were found to be significantly highest at CPS in Òkè-Àkò, however similar trends occurred for Iron in Ìpaò. In contrast, significantly lowest values were recorded for cobalt, and zinc in Ìpaò at NPS.
3.4 Interaction effect of sites and depth on selected phosphorus and sulfur enzymes, and soil nutrient status.
With the exception of the 30-45 cm soil depth, which displayed lower values (1.78 & 0.84 µg/ml/min) at NPS for Tsd and Dsr, respectively, Table 2 highlights significant interaction effects between site and depth for the selected phosphorus and sulfur enzymes. Particularly noteworthy are the interaction effects observed, such as Tsd showing marginally higher values at the 0-15 and 30-45 cm soil depths at NPS and CPS, respectively, with the highest value recorded at the 0-15 cm soil depth at NPS. Additionally, Pho exhibited significant variation, with the highest activity at the 15-30 cm soil depth at CPS and the lowest activity at the 15-30 cm depth at NPS. Furthermore, at CPS, Cobalt, Iron, and Zinc exhibited their highest concentrations at the 15-30 cm soil depth, while the lowest amount of copper was observed at the 30-45 cm depth. Conversely, at NPS, zinc showed its significantly lowest concentration at the 30-45 cm soil depth. Interestingly, copper was found to be significantly highest at the 30-45 cm soil depth at CPS.
3.5 Interaction effects of locations and depth on selected P and S enzymes, and nutrient status
Significant interaction effects between location and soil depth were noticed in soil P and S enzymes (Table 3). Particularly, at the 0-15 cm soil depth, significantly higher activity levels (3.32 mg/ml/min & 9.41 µg/ml/min) were observed for Pho in Òkè-Àkò and Tsd in Ìrèle, respectively. A similar pattern was observed at the 30-45 cm soil depth, where Ìpaò recorded the highest activity levels for Dsr (1.78 µg/ml/min) and at Ìrèle for Tsd (12.51 µg/ml/min) at 15-30 cm depth. However, at the 30-45 cm depth, notably lower values were observed at Ìrèle for phosphatase and Òkè-Àkò for both Tsd and Dsr activity. Phytase activity remained consistent across the locations and depths but exhibited significant variations among means. In addition, location and depths showed significant interaction effect on the soil micronutrients. For instance, In Ìpaò, Iron and manganese were significantly highest at the 0-15 cm depth, whereas in Òkè-Àkò, both elements were found in lower amounts at the same depth. At the 15-30 cm depth, copper and cobalt were significantly highest, but they were located in different locations (Ìpaò & Ìrèle) respectively. Zinc exhibited the highest amount at the 0-15 cm soil depth in Ìrèle and the lowest amount at the 30-45 cm depth in both Ìrèle and Òkè-Àkò.
3.6 Interaction effects of charcoal production sites by location by soil depth on selected P and S enzymes, and soil nutrient status
At CPS, Pho, Tsd, and Dsr activity increased across all locations with increasing soil depth, except in Òkè-Àkò where it decreased (Table 4). Phy showed a consistent trend across all locations and depths at both CPS and NPS. For soil nutrients, Fe decreased across locations and depths. However, Co and Zn followed similar trends, except in Òkè-Àkò. Cu and Mg activity increased across locations with increasing soil depth, but Mn showed no significant difference. At NPS, Cu, Co, and Fe increased across two locations with increasing depth, with the exception in Ìpaò. Mn and Zn decreased down the soil profile at all locations, except in Òkè-Àkò. Pho, Tsd, and Dsr decreased with increasing soil depth, except Pho, which increased at Ìpaò across soil depth.
3.7 Interaction effects of charcoal production sites by location by soil depth on selected P and S enzymes, and soil nutrient status
Three-way analysis of variance presented in [Table 5], indicates the effect of charcoal production on some selected P and S enzymes at the CPS and NPS. Regardless of location and soil depth, there was no significant differences (P > 0.05) indicated in Phy and Tsd between CPS and NPS, but higher and significant differences (P < 0.05) were recorded in CPS with Pho and Dsr. The same table also showed the effect on some selected P and S enzymes at different locations. Irrespective of the production sites and soil depth, there were significant differences in soil P and S status amongst the three different locations. Phosphatase was significantly (P < 0.05) higher in Oke-Ako than in Irele and Ipao, higher significant differences were also recorded in Tsd and Dsr at Irele. Although no statistically detectable differences were indicated amongst the three locations. There were significant differences in soil P and S enzymes status as affected by soil depth (Table 4). Phytase and Dsr decreased in order of increasing soil depth with the lowest value recorded in 30 – 45 cm depth, however, Tsd and Pho with no significant differences recorded higher values in 0-15 and 15 – 30 cm depths respectively. There were significant interaction effects of site by location by soil depth recorded in P and S enzymes. The effect of charcoal production on some soil nutrient status at the CPS and NPS [Table 5]. Showed that Cu, Co, Fe, Mg and Zn were higher in CPS than NPS. In terms of location, Cu, Co, Mn was found to be higher in Ìpaò compared to the other two locations. However, Fe and Zn were found to be higher at Ìrèle with significant differences (P > 0.05). Furthermore, Cu, Co increased with increasing depth while Fe and Mn followed the opposite trend. Finally, Zn remained consistent across the soil profile showing significant difference (P > 0.05).
3.8 Correlation (r) among microbial biomass P, microorganisms, P and S enzyme, and micronutrients
The correlation analysis revealed significant associations among various parameters in the studied samples [Table 6]. The pairwise Pearson correlations, indicate the significance of these associations. Key observations include weak negative correlations between Mn and Zn, as well as Fe and Zn, indicating that fluctuations in Mn and Fe levels are not consistently linked to changes in Zn levels. Similarly, there is a weak positive correlation between Co and Zn, and Cu and Zn, suggesting that variations in Co and Cu levels are not significantly tied to changes in Zn levels. Furthermore, there are weak negative correlations between Nematode (103) and Zn, and Biomass P and Zn, indicating that variations in Nematode (103) and Biomass P levels do not consistently correspond to changes in Zn levels. Conversely, a strong positive correlation is detected between Dsr and Zn, and Tsd and Zn, both of which are statistically significant. This implies that as levels of Dsr and Tsd increase, Zn levels tend to increase.
Transitioning to correlations between different elements, a very weak positive correlation between Fe and Mn was observed, while weak negative correlations exist between Co and Mn, Cu and Mn, Nematode (103) and Mn, Fungi (107) and Mn, Bacteria (108) and Mn, Biomass p and Mn, Dsr and Mn, Tsd and Mn, Phy and Mn, Pho and Mn. None of these correlations are statistically significant, indicating that variations in these elements are not reliably associated with changes in Mn levels. Also, a moderate positive correlation was identified between Co and Fe, and strong negative and moderate negative correlations are noted between Cu and Fe, and Nematode (103) and Fe, respectively. These correlations are statistically significant, suggesting that as Co levels increase, Fe levels tend to increase, while an increase in Cu and a change in Nematode (103) levels are linked to a decrease in Fe levels.
The correlations between Fe and Fungi (107), Bacteria (108), and Biomass P are very weak and not statistically significant. However, a moderate positive correlation is found between Dsr and Fe, and a strong negative correlation is observed between Pho and Fe, both of which are statistically significant. These results indicate that as Dsr levels increase, Fe levels tend to increase, while an increase in Phosphatase (Pho) levels is associated with a decrease in Fe levels. The analysis extends to correlations between elements such as Co, Cu, Nematode (103), Fungi (107), Bacteria (108), and Biomass P, revealing various weak and very weak correlations, none of which are statistically significant. Shifting to Tsd and Fe, a moderate positive correlation is noted, indicating that as Tsd levels increase, Fe levels tend to increase. When exploring relationships with Cu, there is a moderate negative correlation with Co, a very weak negative correlation with Nematode (103), Fungi (107), and Biomass p, and a weak positive correlation with Bacteria (108). While the correlation with Nematode (103) is not statistically significant, correlations with Co and Fungi (107) are, suggesting that as Cu levels increase, Co levels tend to decrease, and an increase in Cu levels is associated with a decrease in Fungi (107) levels. Bacteria (108) exhibits very weak correlations with Nematode (103), Fungi (107), and Biomass P, none of which are statistically significant. Moving on to Phy, there is a weak positive correlation with Co, a very weak positive correlation with Fungi (107), and a very weak negative correlation with Nematode (103). None of these correlations are statistically significant, suggesting that variations in Phy levels are not reliably associated with changes in Co, Fungi (107), or Nematode (103) levels. Phosphatase displays strong positive correlations with Cu, a strong negative correlation with Fe, and a weak positive correlation with Fungi (107), all of which are statistically significant. These findings indicate that as Pho levels increase, Cu levels tend to increase, Fe levels tend to decrease, and there is a weak positive association with Fungi (107) levels.
The analysis concludes with correlations involving Biomass P, showing weak positive correlations with Co and Tsd, and a marginal positive correlation with Phy. None of these correlations are statistically significant, suggesting that variations in Biomass P levels are not reliably associated with changes in Co, Thiosulfate dehydrogenase, or Phy levels. Finally, the correlations between various enzymes Dsr, Tsd, Phy, and Pho are explored, revealing weak and very weak correlations, with only the correlation between Tsd and Dsr being marginally statistically significant.
3.9 Selected soil microbial activities
The biomass phosphorus content in the designated areas exhibited its highest levels at CPS in the 0-15 cm soil depth in Ìrèle, Òkè-Àkò, and Ìpaò [Fig 1]. Subsequently, it gradually declined as the soil depth increased across the profile. Conversely, at NPS, the biomass phosphorus content demonstrated higher patterns at the 0-15 cm depth and reached its peak at the 15-30 cm soil depth in Òkè-Àkò. Then it further decreased with increasing depth across the soil profile. In the study area, the bacteria activity displayed a significant increase in the 0-15 cm soil depth, followed by a gradual decline across the three locations [Fig 2]. Specifically within the NPS, Òkè-Àkò and Ìpaò exhibited significantly higher bacteria activity at a depth of 15-30 cm. At the 30-45 cm soil depth, Òkè-Àkò registered the highest bacteria activity, while Ìpaò showed the lowest. Conversely, at CPS, bacteria were relatively abundant at the 0-15 cm depth in all three locations and at the 15-30 cm depth in Ìrèle; thereafter, it progressively decreased with an increase in depth. Fungal activity was found to be abundant at CPS across various soil depths and in all three locations [Fig 2]. Furthermore, it was observed to intensify with increasing soil depth at CPS. It was significantly higher at 0-15 cm in Ìpaò, at 30-45 cm in Ìrèle, and reached its peak at 0-15 cm soil depth in Òkè-Àkò. Conversely, at NPS, Fungi exhibited a different pattern, with a significant increase at the 15-30 cm soil depth across the three locations. Additionally, at the 30-45 cm depth, fungal activity was notably lower at Ìrèle and Ìpaò. 4Nematode abundance was found to be significantly higher at natural production sites (NPS) compared to charcoal production sites (CPS) in soil layers of 0-15 cm, 15-30 cm, and 30-45 cm across all three locations [Fig 2] . However, there was less variation in nematode abundance at the 0-15 cm soil depth, unlike Ìrèle, where the highest abundance was observed at the 15-30 cm soil depth at NPS. Similar patterns were observed in Òkè-Àkò and Ìpaò.