Exploring the Biodiversity and Conservation Value of Alpine Grasslands in the Bucegi Massif, Romanian Carpathians

: The Carpathian region harbors a wide range of threatened species, making it an area of exceptional conservation value. In the Alpine belt, grasslands cover the entire region and this study aims to describe the communities in the Bucegi Massif of the Romanian Carpathians and highlight their importance for conservation. The Braun–Blanquet approach was used to record ﬂoristic data from 47 phytosociological surveys, identifying a total of 235 plant species from 40 different families, including 30 threatened species. Canonical correspondence analysis was used to analyze the data, revealing that the distribution of vegetation is mainly inﬂuenced by elevation, slope and vegetation cover. Dominant grass species in these communities include Nardus stricta , Festuca violacea , Kobresia myosuroides , Festuca amethystina , Festuca airoides , Sesleria rigida , Festuca versicolor and Festuca carpatica . The alpine and boreal siliceous grasslands of the Carpathian Mountains, identiﬁed by Natura 2000 codes 6150, 6130 and 6170, host a wide range of plant species of signiﬁcant conservation value. The higher altitude grasslands, especially, have outstanding plant species richness. We argue that although the habitats have been grazed, signiﬁcant parts of the area are still in good ecological condition, having many typical natural features.

At the EU level, the Birds and Habitats Directives [18] serve as policy instruments for safeguarding natural and semi-natural habitats, including grasslands and improving their conservation status. The Natura 2000 [19] network comprises sites of community interest, encompassing 198 habitat types listed in Annex I of the Habitats Directive, among them 26 pasture habitats and six grassland habitats threatened by the abandonment of pastoralmanagement practices. To address biodiversity threats in agricultural landscapes [17,18], legal and administrative measures have been implemented. Despite its extensive coverageover 25,000 sites and 1 million km 2 -the Natura 2000 network primarily focuses on habitats with high biological values, notably semi-natural grasslands [19][20][21]. Palearctic grasslands,

1.
Alpine grasslands have become a topic of intense discussion worldwide due to their ecological [2,3,7,10,16,22,33] and cultural significance [14,34]. The Western Carpathians and also the South-Eastern Carpathians, where our study area is located, are no exception to this trend [34,35]. The first in-depth studies of alpine grasslands in this region were carried out by Pus , cariu et al. [36] and Beldie [37], focusing on the endemic flora found in the area and the conservation value of these grasslands [36]. However, a research gap remains regarding the specific ecological impact of intensified sheep grazing on these alpine grasslands in the study area. Recent studies have shown that grasslands on steep slopes possess significant biodiversity value [38,39], making it imperative to investigate the effects of intensified sheep grazing on these grasslands [31]. Understanding the ecological consequences of this grazing practice will help identify potential conservation measures and management strategies to maintain the biodiversity and cultural significance of alpine grassland in the study area [40,41]. In addition to the negative impacts of grazing, activities such as infrastructure development, tourism [10,12,34,35], and climate change-induced effects may contribute to the fragmentation and degradation of habitats in the region [7,15,37]. Thus, the primary objectives of this study are as follows: (I) Identification of alpine grassland plant communities: The study aims to identify and characterize the alpine grassland plant communities present in the study area. (II): Assessment of remaining plant species: The research seeks to assess the diversity and distribution of plant species in the study area, particularly focusing on endemic, rare, and vulnerable species, which play a crucial role in the region's ecological and conservation significance. (III): Description of priority natural habitats: The study also aims to describe and prioritize the natural alpine grassland habitats found in the Bucegi Massif. Understanding the importance of these habitats is essential for formulating effective conservation strategies and protecting the unique biodiversity of the area. (IV) Raising Awareness: The study aims to highlight the significance of endemic, rare and vulnerable plant species; the research intends to contribute to their conservation and management.. floristic survey. The Bucegi Conglomerate, a Cretaceous formation with significant limestone content, is the predominant rock type in the mountains [45]. The location and rugged topography strongly influence the distribution and unilateral or multiple effects of general climatic or microclimatic factors [46]. As a result, the climate of the peaks is cold, with strong north-westerly winds, which drive clouds, leading to rich precipitation. The average annual temperature is about 0 °C and the average temperature of the warmest month is 10 °C [47]. The average annual rainfall at altitudes of 1400-1800 m is about 1100 mm per year [48] (Figure 2). The Bucegi Massif ( Figure 1), with an altitudinal range of 798 to 2505 m and annual average rainfall ranging from 177 to 1423 mm, is a particularly interesting area for the floristic survey. The Bucegi Conglomerate, a Cretaceous formation with significant limestone content, is the predominant rock type in the mountains [45]. The location and rugged topography strongly influence the distribution and unilateral or multiple effects of general climatic or microclimatic factors [46].
As a result, the climate of the peaks is cold, with strong north-westerly winds, which drive clouds, leading to rich precipitation. The average annual temperature is about 0 • C and the average temperature of the warmest month is 10 • C [47]. The average annual rainfall at altitudes of 1400-1800 m is about 1100 mm per year [48] (Figure 2). Given its ecological significance, the Bucegi Massif has been granted conservation status as a Natural Park and is recognized as a Natura 2000 site (RO SCI 0013 Bucegi) [49]. Given its ecological significance, the Bucegi Massif has been granted conservation status as a Natural Park and is recognized as a Natura 2000 site (RO SCI 0013 Bucegi) [49].

Field Data Collection
This study aimed to conduct a survey of the Bucegi Massif, South-Eastern Carpathians, during the summers of 2018 and 2020 to identify the plant communities. The survey involved collecting floristic data from 45 (Table 1) randomly selected plots, each measuring between 4 and 100 m 2 (Table 1). Field sites were selected subjectively. Floristic and geographical data were collected for each plot in stations 1 to 8 ( Figure 1). Other parameters, such as the general topography of the area (altitude, exhibition, and inclination), were also recorded. The Braun-Blanquet approach [50] was utilized to record floristic data on each plot, with nomenclature in line with Romanian studies. Flora Illustrata of Romania [51], Flora României [52] and Flora Europaea [53] were used to identify species. The nomenclature for syntax followed the model of Mucina et al. [54].

Field Data Collection
This study aimed to conduct a survey of the Bucegi Massif, South-Eastern Carpathians, during the summers of 2018 and 2020 to identify the plant communities. The survey involved collecting floristic data from 45 (Table 1) randomly selected plots, each measuring between 4 and 100 m 2 (Table 1). Field sites were selected subjectively. Floristic and geographical data were collected for each plot in stations 1 to 8 ( Figure 1). Other parameters, such as the general topography of the area (altitude, exhibition, and inclination), were also recorded. The Braun-Blanquet approach [50] was utilized to record floristic data on each plot, with nomenclature in line with Romanian studies. Flora Illustrata of Romania [51], Flora României [52] and Flora Europaea [53] were used to identify species. The nomenclature for syntax followed the model of Mucina et al. [54].

Data Analysis
We used the Shannon diversity index (H') to analyze the distribution of species abundance across the entire dataset [55], taking into account both the number of individuals and the number of species present [56]. The evenness index was also used to assess whether the distribution of species abundance in the study area was uniform or skewed towards certain dominant species.
To reduce noise in the data, species that were present in each relevé with a count of five or less (N ≤ 5) were removed [57]. Additionally, both the species response and explanatory variables were logarithmically transformed to allow for the comparison of response variables measured on different scales.
We assessed the distribution of the data using the Shapiro-Wilk test and then performed a canonical correspondence analysis (CCA) to explain the variation in floristic composition based on environmental variables [58]. By conducting the Shapiro-Wilk test, we aimed to meet the assumptions necessary for the selected statistical methods. This allows us to make accurate inferences and draw meaningful conclusions from the data. We used the Monte Carlo permutation test with 499 permutations to determine the significant influence of environmental variables on species composition, with a significance level of p < 0.01. The statistical analysis of the dataset was conducted using the CANOCO version 4.5 software package [58]. To standardize the data for CCA, we transformed the variables to ensure that they have comparable scales and to reduce the influence of outliers using the logarithmic transformation.
Finally, we used the Kendall rank correlation coefficient to assess the correlation between environmental variables and species abundance [57].

Results
Based on our registration of 45 phytosociological relevés, we have identified eight dominant plant communities in the area. These communities are primarily dominated by N. stricta, Festuca violacea, K. myosuroides, F. amethystina, F. airoides, S. rigida, F. versicolor and F. carpatica. In total, we have identified 235 plant species from 40 different families. Among these, 30 species are endemic, rare, or vulnerable to the studied territory, highlighting the importance of conserving this region's biodiversity ( Figure 3). We also calculated the abundance of species (     Our analysis shows that around half of the identified plant species belong to four main families, which are Asteraceae (52 species), Poaceae (37 species), Caryophyllaceae (18 species) and Fabaceae (11 species) ( Figure 4). Our analysis shows that around half of the identified plant species belong to four main families, which are Asteraceae (52 species), Poaceae (37 species), Caryophyllaceae (18 species) and Fabaceae (11 species) ( Figure 4).

Species Diversity
According to Table 2 Table 1.

The Canonical Correspondence Analyses
Canonical correspondence analysis was performed to examine the relationships between floristic composition and environmental variables. The first canonical axis, with a higher eigenvalue (λ1 = 0.45), represents the sampling area as the main gradient, followed by the second and third axes with lower eigenvalues (λ2 = 0.25 and λ3 = 0.11, respectively). Vegetation coverage is represented by the second axis. The total inertia was 4.03. The highest correlations between species and environmental variables were observed for the first (r = 0.81) and second (r = 0.81) axes, whereas the third axis showed the lowest correlation (r = 0.67). The cumulative percentage variance of species data was as follows: 11

Species Diversity
According to Table 2 Table 1.

The Canonical Correspondence Analyses
Canonical correspondence analysis was performed to examine the relationships between floristic composition and environmental variables. The first canonical axis, with a higher eigenvalue (λ 1 = 0.45), represents the sampling area as the main gradient, followed by the second and third axes with lower eigenvalues (λ 2 = 0.25 and λ 3 = 0.11, respectively). Vegetation coverage is represented by the second axis. The total inertia was 4.03. The highest correlations between species and environmental variables were observed for the first (r = 0.81) and second (r = 0.81) axes, whereas the third axis showed the lowest correlation (r = 0.67). The cumulative percentage variance of species data was as follows: 11.4% for Sustainability 2023, 15, 12643 9 of 21 the first axis, 17.7% for the second axis and 20.6% for the third axis. The variance of the species-environment relationship was 46.5% for the first axis, 72.6% for the second axis and 84.2% for the third axis. The Monte Carlo test of significance of the first canonical axis was an F-ratio = 4.999 (p < 0.01) and for all canonical axes, the trace was 0.98 and F-ratio was 2.521 (p < 0.01).
Please note that the results of the canonical correspondence analysis (CCA) should be interpreted with caution due to the potential instability caused by the ratio of variables (species plus environmental variables) to relevés. It is recommended to have a minimum of three times the number of relevés compared to the number of variables for eigen analysisbased multivariate analysis. Given the specific characteristics of our dataset, where the number of variables exceeds the number of relevés, the stability of the CCA results may be compromised.
Agrostis capillaris subsp. capillaris and Veronica officinalis showed the highest positive correlation with axis 1, which is dependent on the sampling area gradient. Agrostis rupestris, F. airoides and Luzula sudetica were highly correlated with axis 2, with their abundance increasing with vegetation coverage. Dryas octopetala, on the other hand, showed a negative correlation with axis 1 ( Figure 5). Please note that the results of the canonical correspondence analysis (CCA) should be interpreted with caution due to the potential instability caused by the ratio of variables (species plus environmental variables) to relevés. It is recommended to have a minimum of three times the number of relevés compared to the number of variables for eigen analysis-based multivariate analysis. Given the specific characteristics of our dataset, where the number of variables exceeds the number of relevés, the stability of the CCA results may be compromised.
Agrostis capillaris subsp. capillaris and Veronica officinalis showed the highest positive correlation with axis 1, which is dependent on the sampling area gradient. Agrostis rupestris, F. airoides and Luzula sudetica were highly correlated with axis 2, with their abundance increasing with vegetation coverage. Dryas octopetala, on the other hand, showed a negative correlation with axis 1 ( Figure 5).      The result of the Pearson correlation coefficient (r = 0.170) indicates a very weak correlation between the plot area and the number of species in a plot (Table 3). This suggests that there is no significant linear relationship between the measured area and the number of species observed in a given area.  Significance levels for correlation coefficients are indicated as follows: * (p < 0.05), ** (p < 0.01), and *** (p < 0.001).
The results of the statistical analysis for alpine grassland variables are presented in Table S1. This legend briefly describes the contents of the table and specifies the variables analyzed. It gives an overview of the results of the statistical analysis, including test statistics such as Shapiro-Wilk W, p-values (normal and Monte Carlo), Jarque-Bera JB and Chi-square values. In addition, it highlights the result of the "OK" Chi-square test for each variable based on the criterion of having a sample size (N) greater than 20.
Based on the information provided, the canonical correspondence analysis was used to understand the relationships between floristic composition and environmental variables. The analysis revealed three canonical axes, with the first axis being the most significant, representing the sampling area as the main gradient. The second axis is represented by vegetation coverage, while the third axis has a lower correlation between species and environmental factors. The canonical correspondence analysis revealed significant relationships between floristic composition and environmental variables, with the sampling area and vegetation coverage being the main factors. The analysis highlights the species that are most strongly correlated with each axis and provides valuable insights into the factors influencing the distribution of plant species in the area. Figure 6 provides information on the elevations of each studied community. In general, grasslands above 1500 m are mainly used for grazing sheep and cattle, and are considered high-mountain pastures [5]. We also compiled a table with notes on threatened species [59] to highlight their presence once again. We also added the conservation value [60] of these grasslands. Moreover, we add a table that provides a comprehensive overview of the conservation status of various communities based on EUNIS2020 [61,62], Red List [63] and NATURA 2000 codes [64]. It includes information on threatened species within each habitat, as well as their corresponding conservative value [61,65] (Table 4). Moreover, we add a table that provides a comprehensive overview of the conservation status of various communities based on EUNIS2020 [61,62], Red List [63] and NATURA 2000 codes [64]. It includes information on threatened species within each habitat, as well as their corresponding conservative value [61,65] (Table 4).

Ecology of Studied Communities
In the Bucegi Massif of the South-Eastern Carpathians, alpine grasslands cover a wide altitudinal range, extending from the upper mountain belt around 1600 m above sea level to the alpine belt at 2400 m [43,44,49]. The dominant graminoid community in the lower alpine zone, where siliceous rocks prevail [45], is N. stricta, while F. rubra and S. rigida communities thrive on calcareous rocks. As we ascend to the upper alpine zone, K. myosuroides and F. violacea communities tend to dominate [4,31]. These alpine grasslands encompass diverse plant communities such as oligo-mesotrophic subalpine grasslands (Potentillo-Nardion) and calcareous alpine grasslands (Festuco-Seslerion bielzii) [35,38]. N. stricta. We noticed in the Western Carpathians [62], alpine grassland communities are classified into Caricetalia curvulae and Seslerietalia variae orders, depending on the bedrock composition. Caricetalia curvulae communities prevail in areas with crystalline bedrock, while Seslerietalia variae communities are more abundant in regions with carbonate bedrock [38]. However, grasslands on acidic and nutrient-poor soils associated with crystalline bedrock display lower diversity but exhibit a remarkably similar composition [62].
We firstly refer to the environmental conditions that might differ between the different communities studied. These include factors such as altitude, temperature, rainfall, soil type, and exposure to sunlight [5,7,63,66]. How do these factors vary in regions where N. stricta, F. rubra, S. rigida, K. myosuroides and F. violacea communities predominate? Are there significant correlations between certain environmental parameters and the distribution of these plant communities? We found a correlation between the abundance of different plant species in each of the studied communities and altitude, with some species increasing in abundance as altitude increases, while others decrease [67,68]. Heegaard [69] proposed that changes in species abundance are influenced by both environmental factors related to temperature as well as biological interactions. Meanwhile, decreases in values of occurrence are likely due to the increased severity of the environment with increasing altitude, which is associated with temperature decreases [70]. There is also a correlation between the abundance of certain plant species and slope, with some species increasing in abundance with slope, while others decrease. Our results are supported by other studies assessing a strong effect of slope aspect on grassland productivity and species composition, in addition to other environmental factors such as water and soil temperature [66,71]. A similar study from the Po plain, karstic and pre-Alpine mountain regions and the western part of the Pannonian plain [17] indicates that the group of soil factors had correlated with both species richness and composition, followed by climatic and topographic factors. Altitude, soil pH, geographical gradient, frequency of flooding, mean annual temperature, date of mowing, humidity, annual precipitation and soil nutrient content were found to be the most important factors in explaining the variance of plant species composition [12,61,63]. Within a community, the floristic diversity is increasing as a function of the number of species and equitability of species abundance [63][64][65]72]. The results show that there is a dependence between species number and evenness in the assessment of floristic diversity. A greater number of species and an equitable distribution of abundance of species contribute to the increasing floristic diversity. According to the coefficient of variation, there is no greater variability of the diversity and evenness, but regarding the number of species within each relevé, it a greater heterogeneity of this variable was recorded.
Vegetation coverage also appears to be a significant factor in determining the abundance of different plant species [14,73], with some species positively correlated with vegetation coverage, while others are negatively correlated. Thus, there are studies that suggest that abundance is affected by species richness and by spatial scale [64]. Anyway, environmental heterogeneity is more important than the area in some studies [74,75].
The CCA analyses reveal a significant influence of the main species gradients. They present changes in plant diversity along the environmental gradients. Great importance is attributed to the distribution pattern along the altitude gradient, for instance, Trifolium repens is more abundant at low altitudes whereas Asperula capitata is more abundant at high altitudes. The slope and vegetation coverage are also gradients that indicate changes within the structure of studied assemblages [76,77]. The result is in accordance with the autecology of the investigated plants [78]. The plants significantly correlated to high altitudes prefer low temperatures and they are growing especially in the Alpine belt, whereas those of low altitudes prefer high temperatures [72,74].
The sampling area as the main gradient described by axis 1 is highly correlated to this axis. Also, the vegetation cover is highly correlated to axis 2. The slope influences, in a positive way, the species abundance, probably due to the restricted accessibility of human activities.
F. violacea, K. myosuroides. These grasslands can either have a secondary origin, resulting from human activities such as grazing, logging or burning, or a primary origin, developing naturally over time. This information is supported by the research of Coldea and Cristea in 1998 [33], and Kyyak in 2004 [15]. We recorded a total of 235 plant species from 40 families, including 30 threatened species. This is consistent with previous studies that suggest that the highest concentration of threatened species can be found at high altitudes in the subalpine and lower alpine altitude belts [33]. Environmental factors, such as geomorphological and climatic variables, have played a significant role in fostering this richness of endemic, rare or vulnerable species in alpine regions [58].
The data suggest that the environment is moderately diverse, with 246 taxa present and a relatively even distribution of individuals among them.
Thus, relevés 1 and 3 (dominated by S. rigida) and 31, 33 and 34 (dominated by F. violacea) have the highest diversity but, at the same time, an equitable distribution of species abundance [59].
In contrast, relevés with a lower number of species (6 and 7 dominated by S. rigida, 8 dominated by F. amethystina and 28 dominated by F. airoides) have lower diversity and an uneven distribution of species abundance [59].
Analyzing the distribution of threatened species in the studied alpine grassland communities, we note an inconsistent distribution with large deviations. We try to explain this as being due to the local climate in alpine grasslands which can vary significantly over short distances, depending on factors such as exposure to wind and sun. This can create microhabitats that favor certain plant species over others [65].
Many alpine grasslands are grazed by herbivores such as cattle, sheep and goats [26,31]. From the point of view of interactions between dominant plant species and other organisms in their respective habitats, we can state, according to the literature [26,31], that there are herbivores in the studied area that preferentially feed on certain plant species, affecting their abundance and distribution [13,17]; this can affect plant distribution by altering resource availability and creating disturbances that favor certain plant species over others.
Our study of 46 alpine grassland communities has revealed a high number of threatened plant species. This is consistent with previous studies that suggest that the highest concentration of threatened plant species can be found at high altitudes in the subalpine and lower alpine altitude belt [5,10,12]. Environmental factors, such as geomorphological and climatic variables, have played a significant role in fostering this richness of endemic, rare and vulnerable species in alpine regions [65].
Comparing similar regions in southern Europe, it appears that mountain isolation has been more conducive to endemism than insularity [66]. Additionally, the degree of endemism tends to decrease as the overall number of plant species increases [67]. Rocky slopes, screes and alpine grasslands have been found to have the highest percentage of threatened plant species, consistent with the general distributional pattern of endemism in high-altitude Eurasian mountains [67].

Conservation of Studied Communities
The findings emphasize the importance of preserving the region's biodiversity. Studies in the high altitudes of the Iranian mountains have highlighted that they are an important hotspot for endemic species [68]. We also support the idea of continued botanical exploration of the mountain ranges, as well as taxonomic reviews of endemism-rich genera that have been understudied [68]. Given that our study area was contained within areas declared as protected, we consider it necessary to highlight habitat framing, even more so because it is recognized that the loss of habitat has been stopped within protected areas, but not in areas beyond their boundaries [70,71]. The EUNIS habitat classification system [72][73][74] is an important resource relied upon by professionals at various stages of nature conservation, including protected area design, resource inventories, monitoring, management planning, impact assessments and ecological restoration target setting [75,76]. In this regard, we have included the EUNIS codes for each community studied and the European Red List of Habitats [56].
We studied different grassland communities in the Alpine belt that correspond to Natura 2000 codes for habitats [55]. Specifically, we found the following.
Communities of F. violacea and of N. stricta correspond to the 6230* species-rich Nardus grasslands in siliceous substrates in the mountain area (priority habitat) code. According to Sarbu et al. [33], they are the most common grasslands in the Bucegi Mountains and they are dominated by N. stricta and F. airoides, with 80% of the species requiring low temperatures and 70% needing moist conditions. In this habitat, we found 15 protected plant species, including rare, endemic and globally protected species that depend on low temperatures and prefer humid soil [77].
Communities of F. airoides that correspond to the 6150 Siliceous alpine and boreal grassland code. These short grasslands are considered a glacial relict and are found on the high crests of the mountains at alpine level [33], and they are dominated by species like Agrostis rupestris, F. airoides and Potentilla aurea subsp. chrysocraspeda. These communities thrive in low-temperature and humid conditions. The majority of the plants in these grasslands require low temperatures and moist soil for growth and development.
Communities of K. myosuroides, F. versicolor, F. amethystina, S. rigida and F. carpatica correspond to the 6170 Alpine and subalpine calcareous grassland code. They have a high biodiversity and are scientifically important for the South-Eastern Carpathians [78,79] where they are endemic [80].
Based on Sarbu et al. [33], grassland habitats 6150 and 6230* are considered to be very sensitive to warming and somewhat sensitive to drought. Moreover, another study revealed a decrease in the biodiversity on mountaintop calcareous alpine grasslands in the central Apennines [80].
In the opinion of García-González [81], Habitat 6170 does not require active management for conservation purposes due to its high structural complexity and fragility. Therefore, the most effective management approach is to refrain from intervening in it.
According to all studies, calcareous grasslands have a high conservation value, while grasslands in siliceous substrates have a moderate conservation value [82]. Moreover, most of them are endemic in the Carpathians: K. myosuroides, F. versicolor, F. amethystina, S. rigida and F. carpatica [82]. Additional field investigations will assess the current condition of these habitats, especially the priority Carpathian habitats [83].
A study of alpine habitats in the Eastern Styrian Alps has shown that climate change will have a significant and immediate impact on their flora and habitats. Fifteen species, including Agrostis rupestris, Carex curvula, Primula minima, Dryas octopetala, Silene acaulis and Saxifraga paniculata, which are also contained in the communities we studied, were identified as particularly vulnerable. By the end of the 21st century, the habitat suitability of these species will decrease significantly [84].
The alpine calcareous and siliceous grasslands of the Carpathian Mountains (identified by Natura 2000 code 6150, 6130, and 6170) are the habitat of a wide range of plant species of significant conservation value. These plants have specific environmental requirements, such as low temperatures, snow cover and high humidity, making them particularly vulnerable to the impacts of climate change. Plant species at risk include Nigritella nigra, a rare species, and Dianthus glacialis subsp. gelidus, an endemic species [33].

Conclusions and Future Prospects
The study highlights the influence of environmental factors, such as altitude, slope and vegetation cover, on the distribution and abundance of different plant communities in the alpine grasslands of the Bucegi Massif. Understanding these relationships is necessary for the conservation and effective management of these ecosystems.
The alpine grasslands in the study area display moderate floristic diversity, hosting a wide variety of plant species, including several threatened species. These grasslands hold significant conservation value due to the presence of endemic, rare, and vulnerable plant species, making their preservation essential for maintaining biodiversity in the South-Eastern Carpathians. The observed distribution patterns of threatened species in the alpine grassland communities highlight the sensitivity of these ecosystems to even minor changes in environmental conditions [81][82][83][84]. The potential impacts of climate change on these habitats warrants further attention and necessitates long-term monitoring and conservation efforts [85][86][87][88].
The presence of high-conservation-value habitats, such as the Natura 2000-designated grassland communities [89], emphasizes the need for effective protection and management measures [84,85]. Proper management practices within protected areas [86], as well as in areas beyond their boundaries, are essential to halt habitat loss and ensure the preservation of these valuable ecosystems [89].
This study underscores the importance of continued botanical exploration and research in the mountain ranges, particularly in relation to taxonomic reviews of endemism-rich genera and threatened species. Additionally, ongoing monitoring of alpine grasslands will provide valuable insights into their responses to changing environmental conditions and help guide conservation efforts.
In conclusion, the study provides valuable insights into the distribution and composition of alpine grassland communities in the Bucegi Massif. The research highlights the importance of these ecosystems for biodiversity conservation and calls for strong conservation measures to protect them in the face of changing environmental conditions. Continued research and monitoring efforts are essential to ensure the long-term survival of these unique and ecologically significant habitats in the South-Eastern Carpathians.
The impact of herbivores on species richness is an important consideration for the conservation of protected areas. According to Fernández-Lugo et al. [2], grazing can have significant effects on vegetation patterns and community heterogeneity.
Grazing can also impact the structural diversity of plant communities. Reference [24] found that grazing had a significant effect on the structural heterogeneity of communities. Inappropriate grazing practices, such as overgrazing, can lead to soil erosion and vegetation changes [90][91][92], which in turn can negatively impact species richness and biodiversity [93].
An assessment of clues on how climate change or human activities affect the distribution and composition of these plant communities should be a priority for future research [94]. It is suspected that there has been a shift in the boundaries of some alpine areas and even an increase in the abundance of invasive species that may be influencing native plant communities [95].