Diversity of Indicator and Dominant Plant Species along Elevation Gradients in Prince Mohammad Bin Salman Nature Reserve

: It is essential to protect and preserve biodiversity, especially in habitats in which natural resources are scarce. The differing ﬂora and vegetation distribution on the broad, arid landscape at the Crown Prince Mohammed Bin Salman Nature Reserve has yet to be investigated. Based on GPS coordination and the transitional zone of plant communities, 48 symmetric plots of 50 × 50 m 2 were distributed evenly over six elevations. In this study, we recorded 70 species in 33 families and elucidated ﬂoristic traits correlated with elevation. High species richness was recorded for the families Fabaceae, Poaceae, Asteraceae, and Chenopodiaceae. High numbers of chamaephyte and phanerophyte species were observed. In the chorotype, the mono-regional component contained 47% of the species, the bi-regional component 35.7%, and multi-regional and worldwide species comprised 10% and 7%, respectively. This study noted the growth habits of 23 herbs, 15 shrubs, 10 trees, and a single species of grass, vine, climber, and mistletoe. Diversity indices, indicator species, dominant plant communities, and soil proﬁles were compared for the deﬁned zones of elevation. Alpha and beta diversity were high at elevations of ≥ 1000, 800, and ≤ 100 ma.s.l., compared to elevations of 600 m, 400 m, and 200 m. The highest species richness and species turnover were recorded at elevations of ≥ 1000, 800, and ≤ 100 m, while species evenness was greater at elevations of 600, 400, and 200 m. Vegetation analyses and indicator species (based on relative abundance) showed species variation with elevation. Species domination was inﬂuenced by physical soil structure and soil chemistry. Microclimates, including temperature and relative humidity variations, were found to be a signiﬁcant driver in the ecosystem, resulting in varying plant diversity and species distribution at different elevations. Through canonical correspondence analysis (CCA), we used an autocorrelation of elevations, plant species, and soil properties to identify three phytogeographic categories that were presumed to be a proxy of microclimate change: Category I: elevations 1000 m and 800 m, including Retama raetam , Zilla Spinosa , and Vachellia gerrardii linked with sandy soil; Category II: elevations 600 m and 400 m, including species Haloxylon salicornicum, Rhazya stricta , and Leptadenia pyrotechnica linked with enriched soils containing CaCO 3 and HCO 3 and having a clay texture; and Category III: elevations 200 m and 100 m, including Zygophyllum coccineum , Tamarix nilotica , and Hyphaene thebaica , which thrived in salinity and silt soils. The spatial vegetation patterns of the xeric environment and its transition zones in Prince Mohammed Bin Salman Nature Reserve were also documented. It is recommended that microclimate effects on species nominated for vegetation restoration or afforestation be considered for the optimal management of this important nature reserve.


Introduction
Elevation represents a global environmental factor linked to decreasing atmospheric pressure and temperature, while abiotic factors such as sunlight and rainfall exhibit localized patterns [1].Changes in non-living and living factors across elevation gradients can lead to different selection pressures, potentially resulting in local adaptation [2].For many Diversity 2023, 15, 1081 2 of 18 years, ecologists have studied how species richness varies along elevational gradients, considering how factors such as climate, geographical area, and habitat diversity contribute to these patterns [3].Research into large-scale ecological patterns has informed a general theory of species diversity and distributions, as well as an understanding of the various local and regional environmental factors that influence them [4,5].
Rapid climate effects coupled with anthropogenic activities diminish biodiversity values worldwide [6,7].These direct and indirect factors lead to ongoing shifts in the status of at-risk species, and other species disappearing from their native habitats [8,9].Spatial heterogeneity is an intrinsic indicator that affects species' incidence and abundance [10,11]; thus, it is essential to consider the environmental factors that impact species distribution as a precursor to assessing quantitative diversity [12].
At a wide spatial scale, elevation ranges above sea level are a critical abiotic factor, leading to variance in environmental gradients and the distribution of natural resources, such as air temperature, air humidity, water, and nutrient availability in soil, and in turn affecting the floristic diversity and spatial pattern of plant species [13].Vegetation composition and plant diversity also vary within a limited habitat range [14].The geographic site, slope, aspect orientation, and altitudinal gradients from the earth's ground level unevenly affect species distribution [15,16].
The distributional patterns of vegetation communities are also influenced by species' responses to habitat heterogeneity, such as adaptation to drought and saline stresses and outcompeting other species to colonize specific habitats [17][18][19].In xeric environments, specific species can dominate a wide range of habitats, which often occurs on a broad regional scale [20].The predominant species can also overlap with monospecific species that have distinctive microhabitats [21].This makes the evaluation of species expansion to other vegetation zones or efforts to determine the spatial scale of a given species both costly and time-consuming [22].
The elevational gradient, with its diverse topographical features, such as mountains, serves as a refuge for species richness, providing protection against the prevalent arid conditions [23,24].In the Kingdom of Saudi Arabia (KSA), biodiversity hotspots are frequently found in areas with significant elevational variations.The Sarawat Mountains in the southwestern region of KSA are a prime example, as they are home to an abundance of plant diversity.Numerous studies have highlighted this area as one of the most diverse habitats [25][26][27].Another notable elevational range is Hijazi Mountain, which is located in the western region and is known for its rich flora [28].However, the plant composition and vegetation structure across different elevational zones in the northwestern region of Saudi Arabia remain largely unexplored.
Prince Mohammad Bin Salman Nature Reserve (PMBSNR) is a new open nature reserve with public access roads comprising both agricultural and rural areas.It encompasses elevation variations involving diverse geological configurations, and environmental heterogeneity, providing a rich and biodiverse ecosystem.This study calculates the important value index to assess vegetation communities, whereas the incipient observations refer to the repetition of dominant and co-dominant species in most vegetation communities.This indicates an overlap of species distribution in the transitional zones at different elevational belts [29].To overcome the species overlap that occurs when IVI is performed, it is worth adding further ecological predictors, such as indicator species (IS), which are often undertaken for monospecific species that colonize unique habitats [30].A perfect indicator of a species is known to be exclusive to that habitat or group without incidence in other habitats [31,32].Therefore, any species in a habitat with a recorded value greater than those of another habitat would be considered a "good indicator species" for that habitat [33].Combined approaches IVI and IS are practical and robust tools that are suitable for this study.
The primary goal of this study is to investigate the floristic traits, diversity metrics, indicator species, and plant communities coupled with soil profiles along elevational gradients in PMBSNR.The assessment of elevational patterns in this study will provide insight into the current state of vegetational distribution and help evaluate the loss of plant diversity, which can drastically alter ecosystems, subsequently threatening the existence of other species.This study also offers a baseline for the development of other programs, i.e., restoration and conservation management schemes, that require urgent action to conserve keystone species, focusing on rare, endangered life forms.

Study Area and Climate
Prince Mohammad Bin Salman Nature Reserve (PMBSNR) was established by royal order in June 2018 and has an area of 16,000 sq.km.This nature reserve is located between two giga projects-NEOM and the Red Sea Project-in northwest Saudi Arabia (Figure 1).The principal aim of Prince Mohammad bin Salman Nature Reserve is to conserve the natural environment of animals and plants and restore ecological balance.PMBSNR has distinctive topography, with variable microclimates.The comparison in the available data meteorology of 5 years (July 2018-December 2022) for Almuwaylih station at a coastline area, elevation ≤ 100 m above sea level (a.s.l.), and Shigry station at a height ≥ 1000 m (a.s.l.), showed that the annual average of precipitation did not differ significantly between these two elevations.However, there were significant differences in the annual averages of relative humidity and temperatures.The general pattern of Almuwaylih station at elevation ≤ 100 m (a.s.l.) indicated a high relative humidity and temperature compared to Shigry station at an elevation ≥ 1000 m (a.s.l.) (Figure 2).gradients in PMBSNR.The assessment of elevational patterns in this study will provide insight into the current state of vegetational distribution and help evaluate the loss of plant diversity, which can drastically alter ecosystems, subsequently threatening the existence of other species.This study also offers a baseline for the development of other programs, i.e., restoration and conservation management schemes, that require urgent action to conserve keystone species, focusing on rare, endangered life forms.

Study Area and Climate
Prince Mohammad Bin Salman Nature Reserve (PMBSNR) was established by royal order in June 2018 and has an area of 16,000 sq.km.This nature reserve is located between two giga projects-NEOM and the Red Sea Project-in northwest Saudi Arabia (Figure 1).The principal aim of Prince Mohammad bin Salman Nature Reserve is to conserve the natural environment of animals and plants and restore ecological balance.PMBSNR has distinctive topography, with variable microclimates.The comparison in the available data meteorology of 5 years (July 2018-December 2022) for Almuwaylih station at a coastline area, elevation ≤ 100 m above sea level (a.s.l.), and Shigry station at a height ≥ 1000 m (a.s.l.), showed that the annual average of precipitation did not differ significantly between these two elevations.However, there were significant differences in the annual averages of relative humidity and temperatures.The general pattern of Almuwaylih station at elevation ≤ 100 m (a.s.l.) indicated a high relative humidity and temperature compared to Shigry station at an elevation ≥ 1000 m (a.s.l.) (Figure 2).

Field and Sampling Protocol
Fieldwork was based on changes in plant species and environmental factors throughout different zones, starting from higher peaks to sea level, whereas almost all topographical configurations exhibited a west-facing aspect.Plant species and soil profiles were sampled along elevational gradients from 1100 m to 10 m a.s.l.(Figure 1).GPS coordination defined the study sites at 200 m a.s.l.intervals of elevation gradients.A total of 48 sites were distributed among six elevations, in which each elevation encompassed eight representative plots of 50 × 50 m 2 , as follows: elevation 1 (≥1000 m) a.s.l., elevation 2 (800 m a.s.l.), elevation 3 (600 m a.s.l.), elevation 4 (400 m a.s.l.), elevation 5 (200 m a.s.l.), elevation 6 (≤100 m a.s.l.).

Floristic and Vegetation Analysis
Soil sampling and analysis: Soil sampling was conducted from the same quadrats, in which vegetative sampling was performed to correlate soil functional traits and species dominance.From each quadrat (n = 8) at each elevation, five random soil samples (0-40 cm) were collected in plastic bags and pooled as one composite sample.Initially, portions from each composite soil sample were separated and collected in duly labeled moisture tins to determine their soil moisture content (MC%) using the weight-loss method.All composite samples were duly labeled and transferred to the laboratory of the Biology Department, College of Science, University of Tabuk, Saudi Arabia, for further analysis.The samples were removed from plastic bags, spread over plastic sheets, air-dried at room temperature, and filtered through a 2-mm sieve to remove any debris.After that, the samples were stored in plastic bags for further analysis.Standard techniques were used for their physical and chemical examination [34].The hydrometer method was used to examine the soil texture of the

Field and Sampling Protocol
Fieldwork was based on changes in plant species and environmental factors throughout different zones, starting from higher peaks to sea level, whereas almost all topographical configurations exhibited a west-facing aspect.Plant species and soil profiles were sampled along elevational gradients from 1100 m to 10 m a.s.l.(Figure 1).GPS coordination defined the study sites at 200 m a.s.l.intervals of elevation gradients.A total of 48 sites were distributed among six elevations, in which each elevation encompassed eight representative plots of 50 × 50 m 2 , as follows: elevation 1 (≥1000 m) a.s.l., elevation 2 (800 m a.s.l.), elevation 3 (600 m a.s.l.), elevation 4 (400 m a.s.l.), elevation 5 (200 m a.s.l.), elevation 6 (≤100 m a.s.l.).

Floristic and Vegetation Analysis
Soil sampling and analysis: Soil sampling was conducted from the same quadrats, in which vegetative sampling was performed to correlate soil functional traits and species dominance.From each quadrat (n = 8) at each elevation, five random soil samples (0-40 cm) were collected in plastic bags and pooled as one composite sample.Initially, portions from each composite soil sample were separated and collected in duly labeled moisture tins to determine their soil moisture content (MC%) using the weight-loss method.All composite samples were duly labeled and transferred to the laboratory of the Biology Department, College of Science, University of Tabuk, Saudi Arabia, for further analysis.The samples were removed from plastic bags, spread over plastic sheets, air-dried at room temperature, and filtered through a 2-mm sieve to remove any debris.After that, the samples were stored in plastic bags for further analysis.Standard techniques were used for their physical and chemical examination [34].The hydrometer method was used to examine the soil texture of the sand, silt, and clay fractions [35].The oxidization method was used to measure the soil organic matter (OM%) using K 2 Cr 2 O 7 [36].Soil water extracts (1:5) were prepared to estimate soil electrical conductivity (EC) and pH [37].A Calcimeter (Eijkelkamp, Agrisearch Equipment) was used to determine the calcium carbonate (CaCO 3 ) content.The titration method was used to examine the soluble anions (Cl and SO 4 ), while for soluble cations (Ca, Mg, Na, and K), a flame photometer was used [38].

Data Analysis
1-Floristic traits were identified based on the number of species, which were categorized by life form, growth form, chorotype, and the number of families along an elevational gradient.The nomenclature follows the KSA taxonomic database [39,40].
2-Diversity metrics were produced using (PAST-Paleontological Statistics software 4.03), which calculated alpha and beta diversity indices.The absence and presence of species was used to present the species richness Taxa-S, Shannon _H, Simpson_1-D, and Evenness_eˆH/S, and Beta-diversity Whittaker's index was used to measure species turnover [41,42].
3-This study presents a multivariate analysis-based Bray-Curtis index, hierarchy cluster analysis, and non-metric multidimensional scaling ordination (NMDS) using rank-based 48 plots via PAST software 4.03.Cluster analysis can be used to interpret the similarities between the plots.However, NMDS is more pronounced.Unlike other ordination techniques, NMDS is useful in cases of non-normality and discontinuous scales, and when some species are absent, resulting in data with zero values [43], followed by the indicator value of species (IS), which is a good tool for highlighting where the species' incidence typically occurs.The indicator value of species (IS) is calculated based on the relative abundance and frequency of the species, ranging from 0 (no indication) to 100 (a perfect indication) [44].
4-The important value index (IV) is based on the sum of relative density plus the relative covering of species.The IV determines the dominant, co-dominant, and other important species, as established by Bonham (2013) [45].
5-Canonical correspondence analysis (CCA) and positive/negative Pearson correlation (the heatmap) were applied to elucidate the relationship between the important value index (IV) and the measured physical and chemical soil variables along elevation variations.
positive correlation with MC.There was a positive correlation between measured param eters pH and salinity and most cations Na, Ca, Mg, and important species H. thebaica an T. nilotica.For other important species, S. holosericea and T. apollinea showed positive co relations with pH, Na, Ca, clay, and silt.As an important species in elevation 6 (≤100 a.s.l.), Zygophyllum coccineum was positively correlated with most soil variables, includin salinity, cations, anions, and MC.A heatmap illustrated Pearson's correlation between dominant, co-dominant, and important species and soil variables (Figure 7).Haloxylon salicornicum and R. stricta showed a positive correlation with K, HCO 3 , CACO 3 , and clay content.In contrast, they showed a negative correlation with pH and MC, while other soil variables were not significant.Vachilla tortilis and V. raddiana had negative correlations with MC.Leptadenia pyrotechnica, L. shawii, and O. baccatus were positively correlated with OM, CaCO 3 , and HCO 3 , but negatively correlated with pH, Na, Ca, and Mg.Cleome droserifolia revealed a negative correlation with OM and MC.Zilla spinosa, V. gerrardii, and R. raetam had a positive correlation with MC.There was a positive correlation between measured parameters pH and salinity and most cations Na, Ca, Mg, and important species H. thebaica and T. nilotica.For other important species, S. holosericea and T. apollinea showed positive correlations with pH, Na, Ca, clay, and silt.As an important species in elevation 6 (≤100 m a.s.l.), Zygophyllum coccineum was positively correlated with most soil variables, including salinity, cations, anions, and MC.

Discussion
A comprehensive survey of diversity traits for all elevations involving various habitats highlighted 70 species belonging to 33 families, of which Fabaceae, Poaceae, Asteraceae, and Chenopodiaceae were the most ubiquitous families, with distinctive efficiency of prevalence, enabling their seeds to be dispersed in xeric environments [46].The studied species were also categorized according to chorological affinities, as 47% of species were deemed monoregional, whereas 11 species were identified as Saharo-Arabian (SA), and 10 as Sudanian (SUD).It was determined that 35.7% were 16 ioregional, whereas 13 species comprised the Saharo-Sindian and Sudano-Zambezian group (SS + SZ).The relatively high contribution of these phytogeographical groups reflects the general arid climate pattern of the Arabian Peninsula [47].The life form of identified species, denoted as chamaephytes and phanerophytes, and growth habits such as herbaceous and shrub forms, dominated all surveyed elevations; these patterns are very common in arid environments, especially in the north-west region of the Arabian Peninsula, in agreement with other studies [48,49].
The present work Investigated the diversity indices observed with elevational variation in PMBSNR.Alpha diversity showed that elevation ≥ 1000 m a.s.l.(El1) and elevation 800 a.s.l.(El2), and elevation ≤ 100 m a.s.l.(El6) recorded high taxa-species richness, Shannon, and Simpson values; this may be attributed to variation in climate features according to our metrology data, i.e., the humidity and temperature were significantly different at these elevational gradients compared to other elevations.Elevations 300 m and 400 m a.s.l.(El3, El4) had low alpha diversity and high evenness values.Analysis of the Bray-Curtis cluster and NMDS showed a likelihood of a resemblance in the environmental conditions of the studied sites; in these habitats, soil analysis (CCA) demonstrated the richness of HCO3, CaCO3, and K, which could be associated with similar plant species that flourished with a high percentage of these chemical elements at elevations 3 and 4. Beta diversity

Discussion
A comprehensive survey of diversity traits for all elevations involving various habitats highlighted 70 species belonging to 33 families, of which Fabaceae, Poaceae, Asteraceae, and Chenopodiaceae were the most ubiquitous families, with distinctive efficiency of prevalence, enabling their seeds to be dispersed in xeric environments [46].The studied species were also categorized according to chorological affinities, as 47% of species were deemed monoregional, whereas 11 species were identified as Saharo-Arabian (SA), and 10 as Sudanian (SUD).It was determined that 35.7% were 16 ioregional, whereas 13 species comprised the Saharo-Sindian and Sudano-Zambezian group (SS + SZ).The relatively high contribution of these phytogeographical groups reflects the general arid climate pattern of the Arabian Peninsula [47].The life form of identified species, denoted as chamaephytes and phanerophytes, and growth habits such as herbaceous and shrub forms, dominated all surveyed elevations; these patterns are very common in arid environments, especially in the north-west region of the Arabian Peninsula, in agreement with other studies [48,49].
The present work Investigated the diversity indices observed with elevational variation in PMBSNR.Alpha diversity showed that elevation ≥ 1000 m a.s.l.(El1) and elevation 800 a.s.l.(El2), and elevation ≤ 100 m a.s.l.(El6) recorded high taxa-species richness, Shannon, and Simpson values; this may be attributed to variation in climate features according to our metrology data, i.e., the humidity and temperature were significantly different at these elevational gradients compared to other elevations.Elevations 300 m and 400 m a.s.l.(El3, El4) had low alpha diversity and high evenness values.Analysis of the Bray-Curtis cluster and NMDS showed a likelihood of a resemblance in the environmental conditions of the studied sites; in these habitats, soil analysis (CCA) demonstrated the richness of HCO 3 , CaCO 3 , and K, which could be associated with similar plant species that flourished with a high percentage of these chemical elements at elevations 3 and 4. Beta diversity (Whittaker index) showed that elevation ≥ 1000 m a.s.l.(El1) and elevation 800 a.s.l.(El2) had the highest species turnover; respective species colonized at a significant distance from the coastline.This suggests that the unique climate variables (lower relative humidity and temperatures) at these elevations have notably influenced the plant species that can thrive in such conditions.These elevations are far from the coastline, so plant species adapted accordingly, exhibiting low salt tolerance and more xerophytes [18,50].
The current study provides insight into indicator species, important values of plant communities, and relative soil in terms of elevational variation.The highest contribution of indicator species (six plant species) recorded at elevation ≤ 100 m (a.s.l.) (El6) may be the result of increased environmental heterogeneity of micro-habitats represented in this range compared to other elevations, which had only one or two species [51].The elevational gradient is often a proxy of climate variation on a broad spatial scale [52].It was found in the present study that elevation ≤ 100 m (a.s.l.) (El6) had the highest relative humidity, incentivizing the development of several indicator species and supporting the co-occurrence of species typical of contrasting habitat types, even within a couple of meters [51].Therefore, changes in habitat patches and species interactions across PMBSNR can affect the distribution and abundance of indicator species [53].
Four species were widespread and overlapping for plant communities at all elevational gradients: Haloxylon salicornicum, Vachellia tortilis subsp.Raddiana, Rhazya stricta, and Vachellia tortilis.These species are unique to desert lands in the Arabian Peninsula [47,54]; it is worth explaining that environmental conditions, particularly soil characteristics, can significantly affect plant communities' colonization efficiency.CCA appears to be more pronounced; accordingly, both H. salicorncum and R. stricta in Elevation 600 (EL3) may be tolerant of soils with high CO 3 , HCO 3 , and CaCO 3 .Moreover, the competitive ability, a merit of these two species, often plays a role in their prevalence, promoting domination in most studied elevations [55].
One key factor is the dominant and co-dominant species' ability to tolerate a broad range of environmental conditions.For example, some plants and animals can thrive in both wet and dry environments or in hot and cold climates.This flexibility may be attributed to their phenotypic plasticity, which is linked to their functional diversity in terms of nutritional strategies [56].Therefore, super-dominant species among elevations, such as Haloxylon salicornicum, can occupy diverse habitats and exploit different resources, increasing their chances of survival at the most elevated gradients [57].
Notably, elevation 3, at which alpha and beta diversity were very low and no species indicator was found, may reflect the increment in environmental stress compared to other elevations.In other words, only species presented by vegetation communities as dominant and co-dominant species may survive in harsh, low-nutrient soil.Other important species colonized at either elevation ≤ 100 m or elevation 200 m (a.s.l.), where the vegetation zone meets the Red Sea shore and is influenced by the coastal atmosphere, exhibited positive correlations with EC, Na, SO 4 , Cl, Ca, and Mg.Most of these species are classified as salttolerant plants (halophytes), which grow in saline environments, in line with the results obtained by [58,59].Other important species, such as Retama retam, Vachellia gerrardii, Zilla spinosa, Lycium shwaii, and Leptadenia pyrotechnica, were correlated with sandy soil as prominent soil properties in elevation ≥ 1000 m and elevation 800 m (a.s.l.).According to field observations at elevations above 1000 m and 800 m a.s.l., the topography contains a high sand percentage, appearing as patches among rugged rocks; thus, these plant species' deep root systems have adapted to grow and develop in these elevational conditions [50].
This study found that PMBSNR had two distinguished microclimates, according to the data obtained from two weather stations: one in the uplands (Shigry station, 1000 m a.s.l.).and one on the coastline (Almuwaylih station, 100 m a.s.l.).Phytogeography of sites at 1000 m a.s.l.belonging to the Saharo-Arabian region as a general pattern of the Arabian Peninsula, which is characterized by high aridity as a part of interior climate zonation, identified species at elevation (≤100 m a.s.l.) would represent the Sudanese region's climate, with high temperatures and high humidity throughout the year [27,47].These variations in microclimates at different elevations lead to disparities in plant diversity metrics (alpha diversity indices, indicator species, and soil characteristics).To clarify, the highest species richness was the obvious diversity index for sites at these comparable elevations, as well as the indicator species (Table 1) Retama raetam and Zilla spinosa, which had distinctive microclimates as they grew on sand soil at elevations above 1000 m; correspondingly, eight indicators were identified as belonging to saline-tolerant species that flourish in high-salinity soil as a distinctive microclimate, and these existed at elevations beginning at 100 m from the seashore.The meteorology data were limited for other elevations (800 m, 600 m, 400 m, 200 m), at which there was no mini-weather station that could help interpret the distribution of the remaining species.However, the soil profile of these elevations can be an indicator of the effect of the local climate on vegetation distribution.

Conclusions
PMBSNR embraces distinguished climate ecosystems, particularly elevations ≥ 1000 m and ≤100 m (a.s.l.), which can create heterogeneity in microhabitats, exhibiting floristic diversity, distinct indicator species, and species domination in plant communities.Ongoing environmental changes, such as variability in soil characteristics, can manipulate plant functional traits, affecting species assemblage across interzonal elevations.Understanding ambient environmental factors is essential when determining ways to manage the xeric ecosystem.This study focused on a field survey and established microsite references from high elevations to coastline zones through phytosociology and soil profiles.Future work could investigate the long-term effects of microclimates at varying elevations.Replanting native species next to parent plants in similar populations and microhabitats using a habitat restoration protocol, such as the soil-seed bank technique, and manipulating seeds and seedlings through a symmetric plot design at each elevation is a possible strategy for creating a sustainable ecosystem.Monitoring microsites is necessary to explore vegetation dynamics and species distribution patterns and to determine how they respond to and survive the heterogeneity of the local climate.Such research could help us identify specific strategies and patterns of restoration that relate to successful management and conservation in PMBSNR.

Figure 2 .
Figure 2. Multiple years of monthly average values of precipitation, relative humidity, and temperature from July 2018 to December 2022.Shed bar denotes data taken from Almuwaylih station at elevation ≤ 100 m a.s.l., and the open bar denotes data taken from Shigry station at elevation ≥ 1000 m a.s.l.(ns) denotes a non-significant p-value; an asterisk (*) denotes a significant p-value ≤ 0.05.Multiple t-test analysis, GraphPad Prism software version 9.5.1.(Data source, National Center Meteorology of KSA).

Figure 2 .
Figure 2. Multiple years of monthly average values of precipitation, relative humidity, and temperature from July 2018 to December 2022.Shed bar denotes data taken from Almuwaylih station at elevation ≤ 100 m a.s.l., and the open bar denotes data taken from Shigry station at elevation ≥ 1000 m a.s.l.(ns) denotes a non-significant p-value; an asterisk (*) denotes a significant p-value ≤ 0.05.Multiple t-test analysis, GraphPad Prism software version 9.5.1.(Data source, National Center Meteorology of KSA).

Table 2 .
The important value of species domination (IVs) at the studied elevations, the most dominant, co-dominance, and important species.Values are the average of the important values of species based on the relative cover and density.