Soil-Forming Factors of High-Elevation Mountains along the East African Rift Valley: The Case of the Mount Guna Volcano, Ethiopia

: The soils of the high-elevation mountains along the East African Rift Valley are poorly understood. Assessing the potential of soils for agriculture, climate change mitigation, and environmental functioning requires insight into how they relate to the factors influencing soil formation. Between 3000 and 4120 m a.s.l., 85 soil profiles of Mount Guna were described and sampled. Standard physicochemical analyses were done on all pedons. Additionally, X-ray diffraction, Al ox and Fe ox content, and P fixation were performed on six selected profiles. Soils on Mount Guna included Andosols, Phaeozems, Leptosols, Regosols, Cambisols, Luvisols, and Vertisols. With increasing elevation, clay content, bulk density, and pH decreased while the C:N ratio remained constant. In contrast, sand, silt, silt-to-clay ratio, SOC, Ntotal, and SOCS increased. With a factor analysis, the soil-forming factors’ elevation/climate could be disentangled from the factor’s parent material as these affect topsoil and subsoil differently. In the ordination based on climate/elevation and parent material, Andosols and Vertisols stood out while other Reference Soil Groups (RSG) showed indistinct patterns. Soil erosion appeared as an additional soil-forming factor not accounted for by the factor analysis. The distribution of the RSG was significantly associated with elevation belts ( p < 0.001), lithology ( p < 0.001), and landcover ( p < 0.003). On the summital ridge, the Andosols were crucial for groundwater storage due to high precipitation. Shallow and stony soils in the mid-elevation belt contributed to runoff generation. Average soil carbon stock ranged from 8.1 to 11 kg C m − 2 in the topsoil and from 29.2 to 31.9 kg C m − 2 in the upper meter, emphasizing the global importance of high-elevation areas for carbon sequestration.


Introduction
Knowledge of the soil properties of high-elevation mountains is important to understand their role in hydrological and global climate processes.As mountains provide over 50% of the world's river sources, climate changes will affect their hydrologic response with large effects on people relying on these water resources for domestic use, agricultural, energy, and industrial purposes [1].While high-elevation mountains, including the mountains along the East African Rift Valley system, play a vital role in climate, water-flow regulation, and the carbon cycle, they are also vulnerable to climate change [2][3][4].Furthermore, a recent literature review [5] highlighted that the dynamics of soil organic carbon Table 1.Number of pedons reported in peer-reviewed publications with detailed data from highelevation mountains (above approximately 3000 m a.s.l.) along the East African Rift Valley system.Mount Guna, in North-Western Ethiopia, is one of these high-elevation mountains, whose soils had not been studied before.The western slopes of this volcano are part of the Lake Tana Basin forming the upper reaches of the Blue Nile Basin.Over the last decades, Lake Tana Basin has received considerable attention from hydrologists and geomorphologists [24][25][26].Given the importance of soils in hydrological processes and their role in climate change mitigation, we aimed at characterizing the soils of Mount Guna.By studying a large number of soil profiles, we aimed at getting a better insight into the relationship between soils and soil-forming factors of high-elevation areas along the East African Rift Valley system.A related, side question, was whether Andosols occur on this volcano, as these have not been mapped for any of the high-elevation mountains in Ethiopia (Figure 1).

Study Area
Mount Guna is a complex volcano with its base at around 3000 m a.s.l., and with its summit reaching 4120 m a.s.l.The mountain has the general appearance of a shield volcano, forming a plateau between 3500 and 3750 m a.s.l. and featuring an eroded caldera crater in the north-western sector (Figure 2).The plateau is overarched by a north-southoriented summital ridge featuring several volcanic necks.The geological base of Mount Guna consists of flood basalts dating back to Mid to Late Miocene [27].The summital ridge  [22] and the Harmonized World Soil Database version 2.0 [23]).
Mount Guna, in North-Western Ethiopia, is one of these high-elevation mountains, whose soils had not been studied before.The western slopes of this volcano are part of the Lake Tana Basin forming the upper reaches of the Blue Nile Basin.Over the last decades, Lake Tana Basin has received considerable attention from hydrologists and geomorphologists [24][25][26].Given the importance of soils in hydrological processes and their role in climate change mitigation, we aimed at characterizing the soils of Mount Guna.By studying a large number of soil profiles, we aimed at getting a better insight into the relationship between soils and soil-forming factors of high-elevation areas along the East African Rift Valley system.A related, side question, was whether Andosols occur on this volcano, as these have not been mapped for any of the high-elevation mountains in Ethiopia (Figure 1).

Study Area
Mount Guna is a complex volcano with its base at around 3000 m a.s.l., and with its summit reaching 4120 m a.s.l.The mountain has the general appearance of a shield volcano, forming a plateau between 3500 and 3750 m a.s.l. and featuring an eroded caldera crater in the north-western sector (Figure 2).The plateau is overarched by a north-south-oriented summital ridge featuring several volcanic necks.The geological base of Mount Guna consists of flood basalts dating back to Mid to Late Miocene [27].The summital ridge consists of the Guna Trachyte, while the slopes of the shield are dominated by the Guna Tuff [27] (Figure 3).Gridded climatic data, with a spatial resolution of about 1 km of WORLDCLIM version 2.0 was used [28].The dataset has average monthly climate variables between the years 1970 and 2000.The long-term mean annual variables were used for this study through aggregate monthly precipitation and temperature datasets (Figure 4).Gridded climatic data, with a spatial resolution of about 1 km of WORLDCLIM version 2.0 was used [28].The dataset has average monthly climate variables between the years 1970 and 2000.The long-term mean annual variables were used for this study through aggregate monthly precipitation and temperature datasets (Figure 4).The landcover units, including tree cover, shrubland, grassland, cropland, and builtup areas, were derived from the ESA WorldCover 10 m 2021 v200 (Figure 5) [29].The primary economic activities on Mount Guna and its surrounding neighbourhood involve livestock, mainly sheep, cattle, and horses, along with crop production.Potatoes and The landcover units, including tree cover, shrubland, grassland, cropland, and built-up areas, were derived from the ESA WorldCover 10 m 2021 v200 (Figure 5) [29].The primary economic activities on Mount Guna and its surrounding neighbourhood involve livestock, mainly sheep, cattle, and horses, along with crop production.Potatoes and barley are commonly cultivated up to an elevation of 3700 m above sea level.

Soil Sampling and Analysis
The field surveys consisted of three phases.First, a reconnaissance survey was conducted to acquire general insights into land use, soils, and the general environment of Mount Guna.Secondly, observations on the soils were made using a standard Edelman auger, and by studying gully and road cuts.Thirdly, based on the gained insights, sites for describing soil profiles were selected.A total of 85 soil profiles were described and sampled along seven topographic transects to cover the elevation belts, besides the lithologic and landcover units (Figure 3).The soil profiles were dug to 100-120 cm, or up to the hard rock.The soil profiles were described following the FAO guidelines for soil description [30] and georeferenced with a handheld geographical positioning system.The profiles were classified according to the 4th edition of the international soil classification system World Reference Base for Soil Resources [31].A total of 226 soil samples were collected from genetic horizons.
Soil bulk density was determined using 100 cm 3 cylindrical steel rings.From each genetic horizon, bulk samples were collected in plastic bags and analysed at the National Soil Testing Center (NSTC) in Ethiopia.The soil samples were air-dried, crushed, and passed through a 2 mm sieve.The analyses of soil characteristics were determined using standard procedures [32].The soil bulk density was determined from undisturbed soil core samples after drying to consistent weights in an oven at 105 °C.
The particle size analysis was determined by the Bouyoucos hydrometer method.Soil pH was measured in a 1:2.5 soil-to-water suspension.The Walkley and Black wet digestion method was used to estimate the organic carbon contents of the soil samples.The Kjeldahl digestion, distillation, and titration method was used to analyse total N.The Olsen method's standard technique was followed to analyse the available phosphorous (P).
Cation exchange capacity (CEC) and exchangeable bases were determined using the ammonium acetate technique.Exchangeable Ca and Mg in the extracts were analysed using an atomic absorption spectrometer (AAS), while Na and K were analysed by a flame photometer.The computation of percentage base saturation (BS) involved dividing the total of the exchangeable bases (Ca, Mg, Na, and K) by the soil's CEC and then multiplying

Soil Sampling and Analysis
The field surveys consisted of three phases.First, a reconnaissance survey was conducted to acquire general insights into land use, soils, and the general environment of Mount Guna.Secondly, observations on the soils were made using a standard Edelman auger, and by studying gully and road cuts.Thirdly, based on the gained insights, sites for describing soil profiles were selected.A total of 85 soil profiles were described and sampled along seven topographic transects to cover the elevation belts, besides the lithologic and landcover units (Figure 3).The soil profiles were dug to 100-120 cm, or up to the hard rock.The soil profiles were described following the FAO guidelines for soil description [30] and georeferenced with a handheld geographical positioning system.The profiles were classified according to the 4th edition of the international soil classification system World Reference Base for Soil Resources [31].A total of 226 soil samples were collected from genetic horizons.
Soil bulk density was determined using 100 cm 3 cylindrical steel rings.From each genetic horizon, bulk samples were collected in plastic bags and analysed at the National Soil Testing Center (NSTC) in Ethiopia.The soil samples were air-dried, crushed, and passed through a 2 mm sieve.The analyses of soil characteristics were determined using standard procedures [32].The soil bulk density was determined from undisturbed soil core samples after drying to consistent weights in an oven at 105 • C.
The particle size analysis was determined by the Bouyoucos hydrometer method.Soil pH was measured in a 1:2.5 soil-to-water suspension.The Walkley and Black wet digestion method was used to estimate the organic carbon contents of the soil samples.The Kjeldahl digestion, distillation, and titration method was used to analyse total N.The Olsen method's standard technique was followed to analyse the available phosphorous (P).
Cation exchange capacity (CEC) and exchangeable bases were determined using the ammonium acetate technique.Exchangeable Ca and Mg in the extracts were analysed using an atomic absorption spectrometer (AAS), while Na and K were analysed by a flame photometer.The computation of percentage base saturation (BS) involved dividing the total of the exchangeable bases (Ca, Mg, Na, and K) by the soil's CEC and then multiplying the result by 100.The CEC-clay ratio was determined by dividing CEC values by the total clay content of the soil and multiplying the value by 100.The Blakemore procedure was used in determining phosphate retention capacity, whereas extractable aluminium and iron were determined by acid oxalate procedure at pH 3.0 following a four-hour shaking of acid solution in a dark room using standard methods.
Using X-ray diffraction (XRD), the nature of the colloids was determined for samples collected in the topsoil (0-25 cm) and the subsoil (25-50 cm) of six pedons (P-71, P-80, P-59, P-73, P-77, P-81) all occurring in the highest elevation belt, under Afro-alpine grassland.These analyses, conducted at Ghent University, Belgium, were aimed at verifying whether these soils could be classified as Andosols.
The samples underwent treatment with 6% NaOCl at pH 8 to remove organic matter [33].Approximately 5 g of fine earth (<2 mm) was crushed with a mortar and pestle and then micronized in a McCrone mill.After the removal of the organic matter, the sand fraction was separated by wet sieving removal.Subsequently, the clay fraction was separated from the silt fraction by successive sedimentation with 2% Na 2 CO 3 as a dispersant, followed by removal of salts by dialysis.The clay samples were analysed with X-ray diffraction (XRD) before and after ethylene-glycol solvation and after heating at 500 • C for 2 h.XRD patterns of air-dried, parallel-oriented clay samples were collected with a Bruker D8 ECO Advance system, equipped with a Cu tube anode and an energy-dispersive position-sensitive LynxEye XE detector.The incident beam was automatically collimated to an irradiated length of 15 mm.The tube operated at 40 kV and 25 mA.Patterns were collected in a θ-2θ geometry from 3 • 2θ to 30 • 2θ, at a step of 0.010 • 2θ, and a count time of 48 s per step.The obtained diffraction patterns were interpreted qualitatively using the COD database [34][35][36].

Statistical Analysis
The data from the physicochemical characteristics were summarised as mean ± standard deviation (SD).To allow us to compare data from soil profiles with different horizons each of varying depths, the calculations are based on the weighted average for three depths defined as the topsoil (0-25 cm), the subsoil (25-50 cm), and the deep subsoil (50-100 cm).The analyses were conducted using the IBM SPSS statistical software (version 26.0) package for Windows [37].
A factor analysis was performed on the same weighted averages.The values of the variables which were not normally distributed were transferred by applying either an exponential function, a square root, a logarithm function, or a combination of square root and logarithm, to ensure that the skewness was ranging between −1 and 1.The first four factor axes, all with eigenvalues > 1, were retained and subjected to varimax rotation [38].The varimax rotation was used to maximize the correlation between the factors and the soil characteristics, enhancing the interpretability of the results [39].We focused on the first four factors, on the hypothesis that these would correspond to topography (elevation), parent material, climate, and land use.The McNemar's chi-square test [40] was used to test the hypothesis that there would be a significant association between Reference Soil Groups (RSG) and elevation, lithology, and landcover units.

Soil Morphological Characteristics
In the majority of the pedons, the particle size distribution exhibited a consistent pattern with respect to elevation and soil depth (Table 2; Figure 6).In terms of texture, the soils in the lower elevation belt (3000-3200 m a.s.l.) had more clay than at higher elevation.Within a profile, the clay content was higher in the subsurface layers.In the lower elevation belt, the moist soil colour changed from very dark grey (10YR3/1) to black (10YR2/1) in the topsoil to dark grey (10YR4/1) and very dark greyish brown (10YR3/2) in the subsoil horizons.In the middle elevation belt (3200-3500 m a.s.l.), it varied from dark brown to very dark grey (7.5YR3/3 to 10YR3/1) and brown/dark yellowish brown to dark reddish brown (5YR3/2).In the highest elevation belt (3500-4200 m a.s.l.), it varied from black (10YR2/1) to yellowish brown or even lighter colours.brown to very dark grey (7.5YR3/3 to 10YR3/1) and brown/dark yellowish brown to dark reddish brown (5YR3/2).In the highest elevation belt (3500-4200 m a.s.l.), it varied from black (10YR2/1) to yellowish brown or even lighter colours.The soil depth varied greatly depending on the topographic position.Most shallow soils occurred in the mid-elevation belt, particularly in upper shoulder and upper back slope positions.Here, soil profile development was limited to an A horizon above a C horizon, above the bedrock.In lower and midslope positions, deeper soils occurred with soil profile development typically exhibiting an A horizon, followed by a B horizon (AB, Bw, or Bt), and a C horizon above the bedrock.
The soil structure of all surface horizons was weakly developed, with a fine, granular structure.In the subsurface horizons, the structure was moderate, medium angular blocky.Likewise, the topsoil had many fine roots in the topsoil that decreased with depth.The horizon boundaries between the A and B horizons were clear and smooth, partly due to a pronounced accumulation of organic matter in the topsoil and partly due to variation in parent material.The dry consistence varied from slightly hard to very hard, whereas the moist consistence varied from friable to firm.On the other hand, the wet consistence ranged from slightly sticky/slightly plastic in the surface layers to very sticky/very plastic in the subsurface soil layers, particularly so in the Vertisols.

Variability in Physicochemical Characteristics along the Toposequence
The physicochemical characteristics of the soils of Mount Guna are summarized in Table 2. Soil texture exhibited distinct differences based on elevation: the sand content increased with higher elevation and consistently remained lower in the topsoil than in the subsoil.Conversely, the clay content rose from the topsoil to the subsoil, while the silt content decreased with soil depth.The soil's bulk density was highest in the lower and mid-elevation belts and was lower near the summit.Within a soil profile, bulk density tended to increase with soil depth.
Overall, the soil pH (H 2 O) decreased with increasing elevation, with the topsoil horizon having lower pH values than the subsoil.Soil organic carbon content (SOC), soil organic carbon stock (SOCs), and total nitrogen (Ntotal) were higher in the high-altitudinal belt compared to lower elevations.These values were consistently higher in the topsoil and decreased with soil depth (Table 2).Overall, the C:N ratio was relatively high, but hardly changed with elevation, suggesting that there was a balance between carbon and nitrogen inputs and outputs in the soil.
In contrast to the previously soil characteristics, the mean values of Ca 2+ , Mg 2+ , CEC, and percentage BS did not change significantly with an elevation and increase only marginally with soil depth.The CEC of the clay also exhibited considerable variation but showed no clear trend in relation to elevation (Table 2).The available phosphorus content was generally relatively low and did not seem to follow a particular trend related to elevation.

Factor Analysis
Table 3 shows the factor loadings of the physicochemical characteristics.For the topsoil, factor 1 (FA-1) had high positive loadings for SOC, SOCs, and Ntotal, and high negative loadings for BD.This trend aligned with the increasing elevation and climate as soil organic matter accumulates at higher elevation, where precipitation is higher and temperature lower (Figure 4).Factor 2 reflected the parent material as it had high positive loadings for sand content and high negative loadings for clay content (Table 3).On Mount Guna, soil weathered from trachyte tends to have higher sand content than soil weathered from basalt or tuff; reciprocally, soil weathered from basalt or tuff has higher clay content (Figure 6).Silt content, on the contrary, is not related to the lithology.Factor 3 had high positive loadings for Ca 2+ , Mg 2+ , CEC, and BS values.These characteristics reflect chemical soil fertility but with no obvious link to a soil-forming factor though it may be affected by parent material and land use.Factor 4 had high negative loadings for silt and silt-to-clay ratio and moderate positive loadings for soil pH, which reflects the degree of weathering.For the subsoil (25-50 cm), factor 1 was determined by high positive loadings for sand content, and reciprocally negative for clay content, which reflected the parent material (Figure 6) and corresponded to factor 1 of the topsoil (Table 4).Factor 2 was determined by chemical soil fertility (Ca 2+ , Mg 2+ , pH, BS, CEC, and av.P), and corresponded with factor 3 of the topsoil.Factor 3 reflected the accumulation of organic matter (SOC, N), which is higher at higher elevation, and corresponded to factor 1 of the topsoil.Factor 4 was determined by high negative loadings on Ntotal and moderate negative loadings on soil pH, and it had high positive loading for the C:N ratio.For the deep subsoil (50-100 cm), factor 1 corresponds again to the parent material as it has high loadings for sand and clay content.Factor 2 is also related to the chemical soil fertility (Ca 2+ , Mg 2+ , and CEC).Factor 3, with high loadings for soil organic matter (SOC, Ntotal, C:N ratio) corresponds to elevation and climate.Factor 4 is determined by high positive loading for BD and moderate positive loading on pH and Ntotal contents, with no obvious link to a soil-forming factor.

Clay Mineralogy
The X-ray diffractograms of the soils which were suspected to be Andosols are illustrated for two profiles in Figure 7.The clay samples showed similar XRD patterns after the different treatments indicating that they had a similar mineralogical composition.The XRD patterns of the untreated clay fraction showed a rational series of reflections at 1.41, 0.71, 0.472, and 0.357 nm indicating a 1.4 nm spacing of phyllosilicate minerals (Figure 7).The 1.41 nm d-spacing was not significantly affected by solvation with ethylene-glycol indicating the absence of swelling 2:1 phyllosilicates (smectites) and the persistence of the 1.41 nm peak after the heating treatment; however, with some decrease in intensity, it confirmed the presence of trioctahedral chlorite.When trioctahedral chlorite was heated to temperatures as high as 500 • C, the peak intensity of the 001 reflection (1.4 nm) normally increased.The decrease in peak intensity of the 001 reflection (1.41 nm peak) in the samples was accompanied by a clear increase in intensity of the 1.0 (0.998) nm peak indicating that some 2:1 phyllosilicate layers collapsed, suggesting the presence of vermiculite besides chlorite.0.358 nm was the resolution of the 004 chlorite and the 002 kaolinite reflections, which can be obtained for rather well-crystallized clays.The peak at 0.447 nm was a hkl-reflection of kaolinite.Kaolinite is dehydroxylated at temperatures around and above 500 °C.The peaks at 0.426 and 0.334 nm were attributed to quartz, while the small peak at 0.325 nm indicated some minor feldspar.The diffuse broad diffraction spectrum between 20 to >30° 2θ was due to the presence of amorphous silicates (opal-A and/or pyroclastic glass).

Soil Classification and Relation with Landcover and Lithology
Following the 4th edition of the international soil classification system World Reference Base for Soil Resources [31], the soils of Mount Guna belong to the Reference Soil Groups (RSG) of the Andosols, Cambisols, Leptosols, Regosols, Luvisols, Phaeozems, and Vertisols.Micas (rational series of reflections at 1.0 (0.998), 0.5 (0.498), and 0.334 nm) were also present in all the samples.The very weak intensity of the 002 reflection (0.5 nm) relative to the 001 reflection (1.0 nm) suggested a trioctahedral structure of the mica (biotite or phlogopite).The double peaks at 0.71 and 0.358 nm almost completely disappeared after heating the samples at 500 • C, indicating the presence of kaolinite.The double peak at 0.358 nm was the resolution of the 004 chlorite and the 002 kaolinite reflections, which can be obtained for rather well-crystallized clays.The peak at 0.447 nm was a hkl-reflection of kaolinite.Kaolinite is dehydroxylated at temperatures around and above 500 • C. The peaks at 0.426 and 0.334 nm were attributed to quartz, while the small peak at 0.325 nm indicated some minor feldspar.The diffuse broad diffraction spectrum between 20 to >30 • 2θ was due to the presence of amorphous silicates (opal-A and/or pyroclastic glass).

Soil Classification and Relation with Landcover and Lithology
Following the 4th edition of the international soil classification system World Reference Base for Soil Resources [31], the soils of Mount Guna belong to the Reference Soil Groups (RSG) of the Andosols, Cambisols, Leptosols, Regosols, Luvisols, Phaeozems, and Vertisols.
The distribution of the RSG per elevation belt, lithology, and their occurrence per landcover unit is presented in Table 4.The most common RSG are the Luvisols (34 pedons), followed by Leptosols (16 pedons) and Vertisols (11 pedons).To ensure a sufficient number of observations in each category for performing the chi-square test, Cambisols, Leptosols, and Regosols were regrouped as one set of shallow or weakly developed soils, while Luvisols and Phaeozems were grouped as another set with well-developed soils.
Andosols are predominantly found in the highest elevation belt (3500-4200 m a.s.l.) formed on trachyte.These soils are mostly under grassland and partly under cropland.
The only two pedons classified as Cambisols are found in the lower elevation belt (3000-3200 m a.s.l.) and are under cropland.The Leptosols and Regosols are mainly found in the lower and partly in the mid-elevation belt (3200-3500 m a.s.l.).The Leptosols formed on tuff can still be used for cropland, whereas the Regosols are more common under tree cover.The Regosols are predominantly formed on tuff, but some are formed on basalt and trachyte.
The Luvisols and Phaeozems occur predominantly in the lower elevation belt and are mostly used as cropland, although some are under tree cover and grassland.Vertisols are mainly in the lower elevation belt, though some still occur in the mid-elevation belt and are mainly found under cropland, with some occurring in shrubland and tree cover.
Both Silandic Andosols, characterized by a low bulk density (≤0.9 kg dm −3 ) and high P retention, and Vitric Andosols, distinguished by the occurrence of volcanic glass, occur.Silandic Andosols require to have oxalate extractable Al + ½ Fe value of at least 2% and a P retention capacity of at least 85%, which is the case for pedons P-71 and P-80 (Table 5).The Vitric Andosols require, besides the presence of volcanic glass in the fine earth fraction, an Al + ½ Fe oxalate extractable value of at least 0.4% and a phosphate retention of at least 25%.This is the case for the other four soil profiles suspected to be Andosols (Table 5).

Discussion
The SCORPAN paradigm of digital soil mapping is that spatial prediction of soil characteristics and soil classes can be made based on covariates relating to the soil-forming factors climate (C), organism (O), relief (R), parent material (P), age (A), besides prior information on soils (S), and the soil's spatial correlation (N) [11].Digital soil mapping at a global scale yet encounters challenges in data availability, scalability, model selection, expertise, and costs.Limited soil data in remote areas lead to incomplete or inaccurate maps [41].By focussing on a remote and under-researched area such as Mount Guna, this study sheds light on the nature of data that are missing from such areas.The results also provide insights into the importance of these areas to the regional hydrology and their role in climate change through carbon sequestration.

Soil-Forming Factors
The soil characteristics of Mount Guna vary greatly depending on the climate, relief, parent material, and landcover.When multiple factors change, it becomes difficult to partition the importance of different factors of soil formation [42].The factor analysis allows us to partly disentangle the intricate effect of the various soil-forming factors.Soil characteristics of the topsoil are predominantly determined by the combined effect of climate and elevation, whereas the characteristics of the subsoils are predominantly determined by the parent material (Table 3).As the climate is directly influenced by the elevation, with higher precipitation and lower temperatures at higher elevation (Figure 4), it is not possible to disentangle the combined effect of climate and topography, whereby we particularly focused on elevation.However, although the lithology of Mount Guna also varies with elevation (Figure 3), the factor analysis allows us to separate the factor parent material from elevation/climate, as it affects topsoil and subsoil characteristics differently.It has been shown that lithology can provide a powerful and easily used covariate to complement other parent material-related covariates and improve the statistical performance of digital soil models and maps [43].As our results show that the topsoil characteristics are more influenced by the combined effect of elevation and climate, this will particularly be the case for subsoil characteristics.
The ordination of the Reference Soil Groups (RSG), based on one hand on the factor climate/elevation and on the other hand the parent material (Figure 8), reveals that Andosols are distinct from the other RSG.Although Andosols are predominantly found on the trachytes of the summital ridge of Mount Guna, it is evident that their distinctiveness is primarily linked to elevation/climate-related soil characteristics, setting them apart from the other soil groups.In terms of these two soil-forming factors, the set of Cambisols, Leptosols, and Regosols, which are soils with limited profile development that are often shallow and stony, and the set of Luvisols and Phaeozems, which are soils with pronounced soil development, seem to be principally distinct from each other in terms of topsoil characteristics.The soils of both sets of soils are predominantly formed in tuff of the lower or middle elevation belt and are extensively used for cropland (Table 4).The shallow and stony soils represent strongly eroded soils, hinting to soil erosion as a significant soil-forming factor.Following the SCORPAN paradigm, conceptually, past soil erosion can be accounted for by the time factor age [11].Practically, in digital soil mapping, geomorphometric covariates representing, e.g., slope steepness, slope curvature, besides landcover, can be used to reflect past soil erosion [44,45].Lastly, the Vertisols also appear as a distinct group.On Mount Guna, these soils are primarily found in the lower elevation belt formed from tuff, and where climate is warm and dry enough for developing the typical shrink and swell features.The topographic association of Andosols at higher elevation and Vertisols at lower elevation is commonly found in volcanic regions [46].

Implication for the Hydrology
Taking hydropedological characteristics into account has been shown to improve the modelling accuracy in Afromontane catchments [47,48].The soil characteristics typically taken into account in hydrological modelling include soil texture, soil structure, soil organic matter, soil water retention, soil depth, and hydraulic conductivity.Although a hydropedological characterisation of the soils of Mount Guna was beyond the scope of the current research, our findings allow us to reflect on the implication of the soils' characteristics on the hydrology.
In the highest elevation belt, on the summital ridge, soils are relatively deep.Andosols, which are the dominant Reference Soil Group (RSG) here, are known for their high porosity, high water-holding capacity, and high hydraulic conductivity [46,49].Given the high precipitation in this area, these soils should be expected to play a key role in groundwater storage and recharge.The soil hydraulic properties of Andosols are however prone to deterioration when used for cultivation or livestock grazing [50,51], and therefore should be protected.

Implication for the Hydrology
Taking hydropedological characteristics into account has been shown to improve the modelling accuracy in Afromontane catchments [47,48].The soil characteristics typically taken into account in hydrological modelling include soil texture, soil structure, soil organic matter, soil water retention, soil depth, and hydraulic conductivity.Although a hydropedological characterisation of the soils of Mount Guna was beyond the scope of the current research, our findings allow us to reflect on the implication of the soils' characteristics on the hydrology.
In the highest elevation belt, on the summital ridge, soils are relatively deep.Andosols, which are the dominant Reference Soil Group (RSG) here, are known for their high porosity, high water-holding capacity, and high hydraulic conductivity [46,49].Given the high precipitation in this area, these soils should be expected to play a key role in groundwater storage and recharge.The soil hydraulic properties of Andosols are however prone to deterioration when used for cultivation or livestock grazing [50,51], and therefore should be protected.
The shallow and stony soils, predominantly found in the mid-elevation belt, should be expected to have a low water-holding capacity, and where slopes are steeper, substantial sources for runoff generation.As these soils are already strongly affected by erosion, soil conservation measures should be promoted there.However, there are concerns about the potential adverse effects resulting from efforts to increase infiltration rates on steep slopes without simultaneously addressing the consequences of saturation excess overland flow development in lower-lying foot-slope positions [52].Inappropriate land management interventions, especially in Vertisol areas, have been shown to dramatically increase gully erosion [53,54].This research also supports previous claims that more emphasis should be placed on understanding the variability of soil properties with soil depth, and the influence this has on catchment hydrology [55].

Potential for Soil Carbon Sequestration
Soil possesses significant carbon sequestration potential, particularly in permafrost areas and soils under grassland, such as Chernozems and Kastanozems [56].Recent studies have emphasized past depletions from the soil organic carbon pool and have indicated a heightened risk of accelerated future losses in response to a warming climate [57].A recent literature review of the Andes region has shed light on the high potential of tropical highelevation areas as both a carbon sink and a carbon source [58].The paucity of soil data for the high-elevation areas of Africa, does not permit a comparable assessment.Soils of Mount Guna harbour in the topsoil 8.1 to 11 kg C m −2 , and 29.2 to 31.9 kg C m −2 in the upper meter (Table 2).For the Afroalpine vegetation belt of the Abune Yosef Mountains, a highelevation mountain to the north-east of Mount Guna (Figure 1), 8.0 kg C m −2 was reported for the topsoil [59].For Mount Kilimanjaro in Tanzania, maximum soil organic carbon stocks in the upper meter have been reported in the range of 94.9 to 142.6 kg m −2 [18].
The values both from Mount Guna and from the Abune Yosef Mountains are lower than those reported for Mount Kilimanjaro but are comparable to SOC stocks reported for Andosols in Hawaii ranging between 30.9 and 62.5 kg C m −2 [60].These findings emphasize the global significance of high-elevation tropical areas for soil organic carbon sequestration.However, they also indicate the need for further research in these regions.

Conclusions
This study, focusing on a remote area, addressed the challenge of data scarcity in global-scale digital soil mapping.By confirming the presence of Andosols in high-elevation areas of Ethiopia, the study filled a gap in existing soil maps for the region.It identified climate/elevation and parent material as key factors influencing the distribution of Reference Soil Groups.Although not explicitly expressed by the factor analysis, a significant association was also found between the Reference Soil Groups and landcover units, reflecting the soil-forming factor of organisms.Soil erosion also emerged as an important soil-forming process, conceptually only accounted for in the SCORPAN paradigm by age.
The results provide insights into the importance of high-elevation areas for regional hydrology, particularly emphasizing the role of Andosols in water storage and groundwater recharge.Furthermore, the data of Mount Guna hint at the importance of high-elevation areas along the East African Rift valley for carbon sequestration.

Figure 1 .
Figure 1.Thirteen high-elevation mountains, or mountain ranges, with summits above 4000 m a.s.l., occur along the East African Rift Valley system.Except for the Rwenzori mountain range, all are volcanoes.The extent of the Andosols is taken from the 2nd version of the Harmonized World Soil Database (authors' cartography based on the Global Multi-resolution Terrain Elevation Data 2010 [22] and the Harmonized World Soil Database version 2.0 [23]).

Figure 1 .
Figure 1.Thirteen high-elevation mountains, or mountain ranges, with summits above 4000 m a.s.l., occur along the East African Rift Valley system.Except for the Rwenzori mountain range, all are volcanoes.The extent of the Andosols is taken from the 2nd version of the Harmonized World Soil Database (authors' cartography based on the Global Multi-resolution Terrain Elevation Data 2010 [22] and the Harmonized World Soil Database version 2.0 [23]).
Soil Syst.2024, 8, x FOR PEER REVIEW 5 of 22 consists of the Guna Trachyte, while the slopes of the shield are dominated by the Guna Tuff[27] (Figure3).

Figure 3 .
Figure 3. Lithological units and location of the 85 pedons on Mount Guna (authors' cartography based on the Geology, Geochemistry and Gravity Survey of the Debre Tabor Area [27]).

Figure 3 .
Figure 3. Lithological units and location of the 85 pedons on Mount Guna (authors' cartography based on the Geology, Geochemistry and Gravity Survey of the Debre Tabor Area [27]).

Figure 4 .
Figure 4. Average annual and monthly rainfall and temperature of Mount Guna for the years 1970-2000 (authors' cartography based on WorldClim 2 data [28]).

Figure 4 .
Figure 4. Average annual and monthly rainfall and temperature of Mount Guna for the years 1970-2000 (authors' cartography based on WorldClim 2 data [28]).

Figure 6 .
Figure 6.Box plots showing the variation in the soil texture at different depths in relation to the lithology of 85 pedons of Mount Guna.

Figure 6 .
Figure 6.Box plots showing the variation in the soil texture at different depths in relation to the lithology of 85 pedons of Mount Guna.

Figure 7 .
Figure 7. XRD patterns of the oriented clay fraction in the topsoil (a,b) and subsoil (c,d) of two suspected Andosols profiles (P71 and P73) reveal the presence of weatherable minerals, including trioctahedral chlorite, vermiculite, micas (biotite or phlogopite), and feldspar, alongside quartz and kaolinite.Additionally, amorphous silicates such as opal-A and/or pyroclastic glass are observed.

Figure 7 .
Figure 7. XRD patterns of the oriented clay fraction in the topsoil (a,b) and subsoil (c,d) of two suspected Andosols profiles (P71 and P73) reveal the presence of weatherable minerals, including trioctahedral chlorite, vermiculite, micas (biotite or phlogopite), and feldspar, alongside quartz and kaolinite.Additionally, amorphous silicates such as opal-A and/or pyroclastic glass are observed.

Figure 8 .
Figure 8. Ordination of the Reference Soil Groups of Mount Guna according to the soil-forming factors elevation/climate and parent material, based on a factor analysis of the physicochemical characteristics of the topsoil, subsoil, and deep subsoil of 85 pedons.The confidence ellipsoids (p = 0.683) are also indicated.Legend: Andosols (AN), soils with noncrystalline minerals of volcanic origin; Cambisols (CM), Leptosols, and Regosols (RG), soils with limited depth or profile development.Luvisols (LV) and Phaeozems (PH), soils with pronounced soil development.Vertisols (VR), soils with shrink-swell clay minerals.

Figure 8 .
Figure 8. Ordination of the Reference Soil Groups of Mount Guna according to the soil-forming factors elevation/climate and parent material, based on a factor analysis of the physicochemical characteristics of the topsoil, subsoil, and deep subsoil of 85 pedons.The confidence ellipsoids (p = 0.683) are also indicated.Legend: Andosols (AN), soils with noncrystalline minerals of volcanic origin; Cambisols (CM), Leptosols, and Regosols (RG), soils with limited depth or profile development.Luvisols (LV) and Phaeozems (PH), soils with pronounced soil development.Vertisols (VR), soils with shrink-swell clay minerals.

Table 2 .
Summary statistics of the physicochemical characteristics (mean ± SD) of 85 soil profiles along seven topographic transects of Mount Guna, Ethiopia.

Table 4 .
Association between Reference Soil Groups (RSG), elevation belts, lithology, and landcover units based on 85 pedons of Mount Guna.
1Set of RSG with shallow or weakly developed soils; 2 set of RSG with well-developed soils.

Table 5 .
Acid oxalate extractable Al and Fe values, and P fixation capacity of topsoil horizons six Andosols of Mount Guna.