Impact of elevation and slope aspect on floristic composition in wadi Elkor, Sarawat Mountain, Saudi Arabia

The current research was carried out in Wadi Elkor, one of the Sarawat Mountains regions, which is a special location from an environmental standpoint and one of the only places in Saudi Arabia where a range of Palaearctic flora co-exists with Afrotropical species. The study aimed to determine the floristic composition as well as the effects of slope aspect and elevation on species, life forms, and phytogeographical elements distribution. The study area is located in Wadi Elkor, a valley in the Sarawat Mountains that cuts off the Al-Hada escarpment, 47 km southeast of Makkah City, Saudi Arabia. We conducted the research at three different locations, each with a different elevation and slope aspect. Based on floristic composition, Ward classification moreover Jaccard comparisons were performed. A total of 189 species was discovered, divided into 131 genera and 43 families. The current study identified Argyrolobium rarum Dumme as a new vascular plant in Saudi Arabia's terrestrial flora. In the current study, the Gramineae, Leguminosae, and Compositae families contributed 29% of the total plant species, whereas 14 families were represented by one species each. Solanum was the most numerous genus, with seven species, followed by Acacia and Pulicaria, each with six species, while Commicarpos and Euphorbia, each with four species. At an elevation of 1060 m above sea level, the north facing slope had the most plants, genera, and families. Therophytes had the most species, accounting for 44%, followed by Chamaephytes, which accounted for 26%. Hemicryptophytes accounted for 12% of the total, while phanerophytes accounted for 10%. In the studied area, the bioregional Saharo-Sindian and Sudano-Zambesian groups were the most well-represented (41%). The floristic composition, as well as the distribution of life forms and phytogeographical components, were found to be significantly affected by the elevation and slope aspects. The study showed that slope aspect and elevation both affected the distribution of plant species, with elevation being the most influential of the two variables.

Sample collection. An extensive survey was carried out during the period of March 2019 to March 2020.
'Three sites were studied in the valley (Fig. 1): Site I is located in the Al-Hada escarpment and is connected to the main bed of Wadi Al-Numan. At an altitude of 750 m, it is the lowest of the three sites. Site II is located on an escarpment of Al-Hada that ranges in elevation from 1060 to 1240 m above sea level with north and south facing aspects. Site III is located at the beginning of Al-Hada mountain, at an altitude of 1830 m. The near-vertical rise in elevation from 1240 to 1830 m meant it was impossible to study vegetation in this elevation belt. To study the effect of slope, aspect and elevation on growth, Site II was divided into six belts, with three different height classes assigned to both sides. The plant samples were collected following the guidelines and legislation of both the Environmental Affairs Agency and the Wildlife Authority in the Taif Governorate. Plant samples were collected from 34 stands (five stands for each belt of elevation, except at 1800 m a.s.l. with four stands), which represented the different habitats of the study area. These are usually placed perpendicular to the slope. In each stand, we www.nature.com/scientificreports/ marked five 10 m × 20 m plots, each separated by 30 m. We purposefully located the plots in areas selected with similar vegetation because of the various habitats in some elevational belts. The study area is not a nature reserve and it is allowed to collect plant samples from it for study. Species to genera ratios in the studied belts were determined and we visually estimated the dominant and codominant species in each site. The lead author, Fadl M.A., identified and named the gathered plant species according to Migahid 39 , Collenette 40 , and Chaudhary 41 .
The collected samples were deposited in the Herbarium of the Biology Department at Taif University and the ID numbers of the voucher specimens are presented in Supplementary Information Statistical analyses. Floristic similarities among elevation belts were assessed by performing a hierarchical classification analysis based on presence/absence data with Wards' (minimum variance) method and Euclidean distances as a dissimilarity measure 47 . The analysis was undertaken using the Statistica statistical software package ver. 8 (StatSoft, Inc., Tulsa, OK, USA). The Jaccard similarity index was applied to evaluate ß-diversity/similarity among stands according to the following formula: Jaccard Index = (the number in both sites)/(the number in either site) * 10), which is based on the presence/absence of species 48 .

Results
Soil characters. The results of soil character showed slight differences in soil texture and electrical conductivity of all sites, including all elevation belts of Site II (  Figure 2 showed that the three sites differed in their taxa numbers, with the middle site exhibiting the highest numbers of species, genera, and families. The fewer taxa numbers were detected at the highest elevation (Site III). The effects of slope aspect and elevation on taxa distribution are presented in Fig. 3. There was no clear pattern of increase or decrease the taxa numbers by increasing the elevation on the south facing side. In contrast to the north facing side, the higher www.nature.com/scientificreports/ the elevation, the higher the taxa number. The results showed significant differences between the north face and south face directions. The numbers of species, genera, and families on the north facing slopes were higher than those in the south facing ones, which were more obvious at an elevation of 1060 m.a.s.l (Fig. 3). Aizoon canariense, Achyranthes aspera, Aerva javanica, Calotropis procera, and Abutilon bidentatum were recorded at all three sites. Within site II, Caralluma retrospiciens, Cenchrus pennisetiformis, Ocimum forsskaolii, Acacia ehrenbergiana, Acacia hamulosa, Sidda alba, Triumfetta flavescens, Coccolus pendulus and Forsskaolea tenacissima were recorded in all elevation belts on both north and south facing sides. The Ward classification of different elevation belts on both sides resulted in a dendrogram (Fig. 4). Figure 4 indicates the effect of both the elevation and the direction on floristic composition, where all belts were divided into two main groups: one represented by 1060 m elevation on the north facing side, and one comprising the remaining belts. The latter is subdivided into 1240 m north and south facing sides, and one with all other categories. The Jaccard similarities showed the one-height belts on both sides (north facing and south facing) were greater than those of different heights in the same direction. Besides, the higher the height, the greater the similarity (Table 3).    Life forms. Five life forms were recorded among the total species found in the three study sites. Therophytes exhibited the highest species with 44%, followed by chamaephytes with 26%. Hemicryptophytes and phanerophytes were represented by 12% and 10%, respectively, whereas the lowest life form, 7.8%, was geophytes with 14 species (Table 4 Chorological affinities. The chorological composition is summarized in Table 5. Generally, monoregional elements exhibited 26% of the recorded species, while biregional elements recorded 54%. Pluriregional species were represented by 7.8%. The Pantropic elements showed the highest numbers, 9.4%, within the monoregional elements; followed by Saharo-Sindian and Palaeotropical elements, each represented by 4%. The bioregional Saharo-Sindian and Sudano-Zambesian groups were the most represented in the studied area (41%). Cosmopolitan taxa accounted for 8% of the total species, whereas the endemic species, Kickxia pseudoscoparia, was recorded in the study area. The chorological affinities shown were variable in their distribution in the studied area. The Pluriregional elements in Site (III) were more than those recorded in Site (I) by 145 folds, with the  www.nature.com/scientificreports/ highest percentage recorded on the north-facing slope at 1060 m a.s.l. In contrast, the Sudano-Zambesian elements have completely disappeared from Site III, Table 6. The bioregional Saharo-Sindian/Sudano-Zambesian, which recorded the highest chorotype in the studied area, increased facing aside from the north facing side in all elevations, an increase ranging from between 13.7 and 35%. Palaeotropical elements decreased by increasing elevation (Table 6) and their proportions were higher on the north facing side than those on the south facing site. In contrast, the bioregional Mediterranean/Irano-Turanian increased at a higher elevation. A slight change was recorded for the Saharo-Sindian elements between Site I and Site III, while it was recorded with higher proportions on the second site, on all elevations for both north and south-facing sides.

Discussion
The study confirmed that both slope aspect and elevation affect the distribution of plant species, with elevation having the greater impact. In this study, 35% of the recorded families were represented by only one species per family, a prevalent characteristic of desert flora, and considered a sign of plant adaptability to xeric conditions 6 . A large number of families recorded with only one species indicated that only a few species of the many plant families adapt and survive in harsh environments. The species to genera ratio in the current study (1.46) was close to those reported in previous studies of other regions in the Kingdom of Saudi Arabia 6 . According to Magurran 49 and Pielou 50 , taxonomic diversity is greater in an environment where species are distributed across several genera than in one where the majority of species belong to the same genus. The species per genus ratio is higher than that recorded in the Al-Shafa Highland region 5 (1.3), and it is less than that reported in the entire Saudi Arabia 6 , 2.6, and lower than the estimated ratios in Khulais 4 . Consequently, the taxonomic diversity of the study area is less diverse than in the Al-Shafa heights and greater than in the Khulais region and all of Saudi Arabia. This result indicates that the taxonomic diversity is increasing with increasing height in the west of the Kingdom of Saudi Arabia. The dominance of Gramineae, Leguminosae, and Compositae and the patterns of the most different families were consistent with previous studies of other regions in Saudi Arabia, such as Khulais 4 , Al-shafa highland 5 , the Asir Mountains 51 , and Hail 52 . Moreover, the current results are like those of neighboring countries, which include Egypt 53 and Sudan, in the Marra Mountains 43 . Only one endemic species was recorded inside the studied area, Kickxia pseudoscoparia. This is in agreement with Al-Nafie 6 , who stated that endemism in Saudi Arabia is insignificant in comparison to some of the adjoining nations, such as Yemen and Oman. The highest taxa numbers in the middle part of the study area (Site II) may be attributed to the minor impact of human disturbances in this section. This site is far from the road, while the first and third sites are close to it (Fig. 1). The observed widespread species in all elevation belts in the north facing and south facing sides could be attributed to their wide ecological amplitude. The contents of total nitrogen, organic matter, and moisture in the soils of the north facing side at 1060 m a.s.l have made suitable conditions for many plants, resulting in a large increase in all taxa numbers.
It has been recorded that in areas where annual precipitation is less than 600 mm, moisture becomes an important factor in the composition, intensity, and structure of plant communities 33 . Besides, it was argued that soil moisture was the main environmental issue governing plant species composition and richness 54 .
Small differences in soil analysis were found among the studied belts, except at 1060 m a.s.l on the north facing slope. This result is not in line with Kutiel and Lavee 33 , who stated that considerable variations in soil properties between the opposing directions can be predicted in arid zones (< 400 mm of annual precipitation) due to low rainfall and high potential evaporation. Our finding is in line with Burke's 55 study of Namibia's arid Nama Karoo inselbergs.
The higher numbers of recorded taxa in all the north facing belts than those in the south facing belts can be explained because the south facing slopes receive more sunlight and become warmer and xeric, helping droughtresistant flora, and are less conducive to tree growth, while the north facing slopes preserve moisture and are humid, supporting moisture-loving plants 56 . It has been reported that insolation affects greatly the composition of plant species, particularly in arid and semiarid regions at low latitudes, through its impact on the water balance 57 . In addition, high evapotranspiration on the south facing slopes may cause drought stress during the growing season. Acacia tortilis, which was the dominant species at the lowest site, was recorded as well adapted to disturbances such as drought, fire, browsing, and pollarding 58 . Ficus palmata, the dominant species on the upper site is distributed from Egypt to Central Asia and is typical of the cold desert mountains of this area. Sometimes Ficus palmata subsp. virgata is also confused with its closely allied species, F. carica, a native of the Mediterranean region to Afghanistan 59 .
The high percentage of therophytes (44%) in the study area agrees with previous studies in other parts of Saudi Arabia 4,6,40,41,52,60-62 . A similar result was published for the northern part of the Eastern Desert of Egypt 63 . Table 6. Chorological type's percentage (actual counts between brackets) recorded in the three different sites, NF, north facing; SF, south facing.   13 , who stated that therophytes and chamaephytes increased in dominance and/or abundance with increasing elevation. In addition, the therophytes percentage in the current study mediates the proportions of the Khulais, with a lower elevation 4 , and the Alshafa highland regions, with a higher elevation 5 . Chamaephytes distribution in the current study contrasts with the observations of Dickoré and Nüsser 71 , who reported that chamaephytes decrease with increasing elevation. This discrepancy between the current study and the study of Dickoré and Nüsser 71 could be attributed to the large difference in elevation (1200-1400 m a.s.l.) to 3000-3400 m a.s.l. Phanerophytes distribution in our results agrees with the finding of Hoffmann and Hoffmann 72 , who reported that phanerophytes decrease in a temperate area and are very common in tropical regions. The percent of Phanerophytes recorded in this study agrees with White and Leonard 73 statement that "the southern and southwestern Arabian Peninsula is very poor in trees".
Numerous early studies express that the flora and vegetation of Saudi Arabia have included elements from two main phytogeographical regions that cover much of North Africa and East Asia, the Sudano-Zambezian region and the Saharo-Arabian-Sindian region [71][72][73][74][75] . The studied region belongs to the Nubo-Sindian Province, which is a part of the Sudanian region 44 . Wickens 43 recorded that the current examined area is principally impacted by the Saharo-Sindian element, whereas White and Leonard 73 stated that the Sudano-Zambezian region stretches out into southern and western Arabia. These influences clarify why the highest percentages for monoregional, bioregional, and pluriegional species have been recorded for the Saharo-Sindian elements, followed by those of the Sudano-Zambesian. In addition, the White and Leonard 73 statement explains that the bioregional Saharo-Sindian and Sudano-Zambesian groups were the most represented in the studied area. Alsherif et al. 4 recorded that Saharo-Arabian and Sudanian elements exhibited the highest chronological elements in the Khulais region, a region away from the current location by 140 km in the northwest direction. Al-sherif and Fadl 5 recorded that Saharo-Sindian had the most bioregional groups in the Alshafa highlands (1700-2300 m a.s.l). The biregional elements Saharo-Sindian and Sudano-Zambesian biregional elements were more frequent on the southern rather than northern facing slopes because the southern slopes experienced warmer and drier conditions than the northern slopes. Sudano-Zambesian elements were most frequent at lower elevations and decreased in frequency with rising elevation, possibly because these elements rely on tropical climatic conditions6. In contrast, a high percentage of Saharo-Sindian and Irano-Turanian elements were found at high elevations 6 . On the contrary, a high percent of Saharo-Sindian and Irano-Turanian elements were found at high elevations.

Conclusions
In the current study of a valley that cuts through the Sarawat Mountains, we identified a range of Palaearctic flora co-existing with Afrotropical plants. It is the first study of the influence of slope direction in the distribution of plant species in the Sarawat Mountains, and it sheds light on the study area in terms of its floristic significance. The study found that both slope aspects and elevation have an impact on plant distribution, with the influence of elevation being the most important. Since south facing slopes receive more sunlight and become warmer and xeric, the higher numbers of the reported taxa in all north facing belts than in the south facing belts can be explained. The report of the new addition to the flora of Saudi Arabia (Argyrolobium rarum) showed that the Sarawat mountains need careful botanical exploration. We recommend using LiDAR data as part of future research.