Translator Disclaimer
16 December 2020 Floristic Diversity and Distribution Patterns Along an Elevational Gradient in the Northern Part of the Ardabil Province Rangelands, Iran
Sahar Ghafari, Ardavan Ghorbani, Mehdi Moameri, Raoof Mostafazadeh, Mahmood Bidarlord, Azad Kakehmami
Author Affiliations +

This study evaluated the distribution pattern, species diversity index (richness, diversity, and evenness), importance value index (IVI), and family importance value (FIV) of all vascular plant species and various plant community life-form spectra along the Moghan Plain–Sabalan Mountain rangelands, Ardabil province, Iran. Sampling was conducted in 11 elevation classes at 300-m elevation intervals (100–3300 masl). In each elevation interval, 30 quadrats (1 × 1 m) were laid to collect vegetation data. Different vegetation attributes (density, frequency, and canopy cover) were measured for each quadrat. Regression analyses were employed to explore the interrelation of elevation with diversity, species richness, evenness, and IVI. In total, 251 species (143 genera, 38 families) were recorded across the study area. Poaceae, Fabaceae, and Asteraceae were most dominant families according to the FIV. According to the IVI, Poa bulbosa, Festuca ovina, Medicago minima, and Artemisia austriaca were the dominant species along the elevation gradient. Total diversity and species richness showed a normal distribution along the elevation gradient. Overall, hemicryptophytes and therophytes were the dominant life forms, while chamaephytes and geophytes were less frequent forms. The life-form patterns changed along elevation gradients. While the diversity and species richness values of therophytes and chamaephytes showed a decreasing trend, hemicryptophytes increased with increasing elevation. Diversity, species richness, and evenness of geophytes were not significantly different across the elevation gradient. Class I (highly palatable) and class II (mostly palatable) species declined as the elevation increased and then increased, while class III (hardly or unpalatable) species showed a reverse trend. The difference in the abundance and distribution of species in elevation classes could be related to resource availability, overlap of habitats, habitat patchiness, land area, degree of human influence, or biotic disturbances.


The Moghan Plain–Sabalan Mountain rangelands (MPSMRs) in Ardabil province, Iran, supply livestock forage and play an essential role in food provision, groundwater enrichment, clean water production, wildlife, tourism, recreation, medicinal plant provision, beekeeping, and climate change mitigation (Ghorbani et al 2014, 2017, 2018; Nazari Anbaran et al 2016; Ghafari et al 2017, 2018). According to the literature, 396 plant taxa have been identified in the area. Of these, 25 plant taxa are endemic to Iran (Ghafari et al 2017, 2018). Because of the variety of plant species and the climate, beekeeping is a good source of food in this area (330,000 kg/y honey products; Sharifi et al 2017). The MPSMRs are also rich in medicinal and commercially valuable plants, which require sustainable management to promote health, income, and employment (Ghorbani et al 2017). In recent decades, with the increasing population, low level of development, and high dependence of local communities on rangelands for their primary livelihood needs, as well as increasing livestock numbers as living standards have improved, the MPSMRs have been overused. Excessive and irregular grazing, destruction by livestock, disregard for seasonal unsuitability for grazing, conversion of rangelands, and road construction are the main causes of degradation in the MPSMRs (Ghafari et al 2017, 2018). In the early 1980s, the rangelands covered 58.2% of the land area; however, by 2015, its coverage had been reduced to 53.9% (Kakehmami et al 2017). Understanding the capabilities and limitations of rangelands is essential for appropriate management at the regional level. The rapid change in elevation, even over short distances, makes the MPSMRs an ideal place for studying the vegetation dynamics associated with biotic and abiotic factors that vary with elevation (Körner et al 2011).

Few studies conducted in the MPSMRs (Ghorbani et al 2014, 2017; Nazari Anbaran et al 2016; Ghafari et al 2017, 2018) have focused on vegetation structure and species richness. However, understanding diversity patterns is essential for developing best approaches to diversity conservation as well as the sustainable utilization and management of plant resources. Hence, the aim of this study was to analyze floristic composition, diversity, and structure along the elevation gradient of the MPSMRs. This is expected to provide baseline data along with an analysis that can be used to develop diversity conservation strategies and sustainable rangeland management plans for the study area. This study attempted to answer the following questions: (1) Do species composition, ecological importance, and structural components of vegetation vary significantly along the elevation gradient? (2) What is the relation between elevation and species richness? (3) Do the diversity index and ecological importance along the elevation gradient differ, and if so, how?

Rangeland history and biophysical specificities of the area

The Ardabil province rangelands are located at the confluence of the 3 climatic areas: the Siberian from the north, the Mediterranean from the west, and the wet from the Caspian Sea. This, together with the Sabalan and Khoroslou mountains and hills (slopes 5% to more than 100%), affects climate and plant diversity in the region (Sharifi et al 2017). Various types of rangelands, such as summer rangelands (Sabalan, Baghru, and Buzghush mountains rangelands), spring and autumn rangelands (low mountainous areas, such as Kalantar-Khoroslou, to hilly areas, such as Gharesou river landscapes, also known as the Arshagh rangelands) and winter rangelands (Moghan Plain rangelands) are part of Ardabil province.

The second largest nomad tribe, the Shahsavan, belong to Ardabil province (Hosseinzadeh et al 2017a). The nomadic livestock husbandry system is one of the most important and oldest livestock husbandry systems in Iran. It strongly depends on natural rangelands. For centuries, the Moghan Plain has served as a grazing area for nomadic tribes, with agriculture being of secondary importance. Shahsavan nomads move their livestock to make the best use of available forage and water in a seasonal cycle on their allotted spring, summer, autumn, and winter rangelands (Hosseinzadeh et al 2017b). Shahsavan nomads use the Moghan rangelands for about 6 months of the year (Arjmand Karkazloo 2016).

The climate is favorable for agriculture and animal husbandry. Traditional and semi-industrial animal husbandry are practiced in villages (paddocks) and are dependent on rangeland forage and agricultural residues. In Ardabil province, 1600 pastoralist households (11,000 people) and their livestock (about 3,400,000 animal units/sheep) contribute at least $57 million per year (milk, meat, dairy, and other by-products) to the province's economy (Sharifi et al 2017). Sheep are the dominant livestock (95%), followed by goats (3%) and cows (2%) (Nazari Anbaran et al 2016; Ghafari et al 2018).

Recently, more extensive irrigated farmlands and agroindustrial schemes have been introduced in the winter rangelands of the Moghan Plain, which, together with the expansion of farming, has further reduced the rangelands and increased the settlement rate among nomads, causing drastic social and economic changes in their way of life (Arjmand Karkazloo 2016). In spring (April to late June) and autumn (mid-September to end of November) the herds feed on the heavily grazed transition routes between winter and summer rangelands. In recent years, spring and autumn rangelands have been cultivated in residential areas, irrigated, dry-farmed, and used for horticulture. The small remaining area is rangeland (Ghafari et al 2018). While overall the rangelands of Ardabil province are now becoming better understood, little is known about its spring and autumn rangelands. At the end of April, as the snow recedes and rangeland becomes available, the nomads move to their summer Sabalan Mountain rangelands at higher elevations. Summer rangelands have prolonged snow cover in winter and comparatively high annual precipitation (400–700 mm). They mark the upper limit of agriculture and permanent human settlement. The harsh climate (extremely low temperature and strong winds), rugged terrain, and geological and hydrological conditions make the Sabalan Mountain rangelands particularly vulnerable to unsound practices of natural resource management and environmental degradation. Erosion rates are higher and the loss of fertility due to nutrient leaching is more pronounced than elsewhere. Because of the prevailing low temperatures, plant growth and soil formation are slower, and vegetation cover is significantly less abundant than in lowland areas (Sharifi et al 2017).

Study sites

The study was carried out in the area stretching from the Moghan Plains to Sabalan Mountain in Ardabil province, Iran, located at 47°45′ to 48°23′E and 38°18′ to 39°27′N (Figure 1). Data collection sites have the following characteristics: (1) stratified at approximately 300 m elevation intervals; (2) located in the natural rangelands (depending on road accessibility corridors selected); and (3) a minimum distance of 300 m from the road, agricultural lands, and residential and other institutional compounds to avoid edge effects. According to the criteria used in selecting sites, the study sites were in the elevational range of 100–3300 masl (anthropogenic ecosystems, such as agricultural lands and residential and industrial areas, dominate elevations below 100 masl, and there are no access roads to sites above 3300 masl). For this study, we categorized the elevation zone between 100 and 3300 masl into 3 elevation categories: low elevation (100–900 masl), mid-elevation (900–2100 masl), and high elevation (2100–3300 masl). The climate at higher elevations is cold semiarid, with cold winters and temperate summers; however, lower elevations have a temperate semiarid climate and most of the precipitation is during winter. Higher elevations are covered with snow and ice for 3–4 months of the year (Tavosi and Delara 2011). More details about the study sites are presented in Table 1.


Geographical location of the area under study (Moghan Plain–Sabalan Mountain rangelands).



Environmental characteristics of 11 elevation classes along elevation gradient of Moghan Plain–Sabalan Mountain rangelands.



The selected rangelands were sampled from April 2015 to July 2016. Within each of the elevation classes, three 100-m transects were set at 50-m intervals parallel to each other and perpendicular to the main slope. Along each transect, 10 quadrats (1 × 1 m2) (at 10-m intervals) were set for systematic random sampling. Since diversity and species richness are quadrat size and abundance dependent (Magurran 2004), the same quadrat size was used for both. Vegetation attributes (density, canopy cover, and frequency), along with spatial coordinates (latitude, longitude, and elevation) using a global positioning system, were recorded for each quadrat.

Species were identified using literature such as Assadi et al (1988–2012), Davis (1965–1988), Ghahreman (1979–2003), Komarov (1934–1954), and Rechinger (1963–2010). The name of taxon authors was verified using IPNI (International Plant Names Index; 2018) and The Plant List (2018). The life forms were identified using Raunkiær classification (Raunkiaer 1934). Species were sorted into palatability classes: highly palatable (class I)—forage plants, most preferred by domestic livestock; mostly palatable (class II)—plants averagely preferred by livestock; hardly or unpalatable (class III)—plants least preferred or not grazed by livestock at any stage; and harmful and possibly toxic (Hussain and Durrani 2009). The palatability classes for each species were determined using literature such as FRWO (1982) and Ghafari et al (2017).

A complete picture of the ecological importance of each species in the community was obtained by calculating the importance value index (IVI) and the family importance value (FIV). To determine the IVI, the percentage values of relative density, relative dominance, and relative frequency were added together (Curtis and McIntosh 1950). The IVI was calculated for each species in each of the 11 elevation classes to determine the vegetation structure. The plant species with the highest importance value at each site were considered to be dominant species, and the plant species with the lowest importance value at each site were considered rare (Curtis and McIntosh 1950).

The FIV for botanical families was calculated by summing the IVI for different species of the same family (Mori et al 1983). IVI and FIV were calculated using the equations in Table 2.


Equations for importance value index (IVI), family importance value (FIV), species diversity, species richness, and evenness.


The vegetation was analyzed for species diversity, richness, and evenness to reduce the complexity of data (Whittaker 1953). Species diversity (Shannon and Weaver 1963), richness (Menhinick 1964), and evenness (Pielou 1969) were determined using the equations in Table 2.

Species diversity indices (species diversity, richness, and evenness) for the 4 life forms (hemicryptophyte, therophyte, chamaephyte, and geophyte) were obtained using Past 3.1 software (Øyvind Hammer, Natural History Museum, University of Oslo, Oslo, Norway, 2013). Regression models were fitted to assess the interrelations between elevation and the 3 diversity indices (diversity, richness, and evenness) and IVI of total plant and variety life form. For each model, the usual statistics such as R2 and partial coefficients were obtained. Regression analyses were performed using Minitab 17 (Minitab Inc. State College, PA, USA, 2013).


Floristic composition, family dominance, and their distribution

Overall, 251 species of 143 genera belonging to 38 families were recorded along the elevation gradient on the MPSMRs (Supplemental material,  Table S1 (mred-40-01-09_s01.pdf); The highest species richness was observed for the families Asteraceae, Poaceae, and Fabaceae. The dominant genera in the sampling sites were Astragalus (11 species) and Allium, Trigonella, and Veronica (6 species each). The number of species varied significantly along the elevation gradient, with the highest value, 74, being recorded at the mid-elevation sites (1800–2100 masl), and the lowest value, 28, being recorded at the low-elevation sites (100–300 masl) (Table 3). Five families (Asteraceae, Caryophyllaceae, Fabaceae, Poaceae, and Scrophulariaceae) showed great potential to adapt across a wide elevational range (Supplemental material,  Table S1 (mred-40-01-09_s01.pdf);


Number of species, genera, families, life forms, and palatability classes a) along elevation gradient of the Moghan Plain–Sabalan Mountain rangelands.


Among families, Poaceae, followed by Fabaceae, Asteraceae, Caryophyllaceae, and Brassicaceae, had the highest FIV. Malvaceae, Apocynaceae, Violaceae, Berberidaceae, and Dipsacaceae were recorded with the lowest FIV in the study area (Table 4).


Family importance values along the elevation gradient.a)


Trifolium subterraneum had the highest IVI value in a single elevation class, reaching 124.95 at 600–900 masl (Supplemental material,  Table S1 (mred-40-01-09_s01.pdf); Poa bulbosa was distributed across the widest range of elevation classes and reached the highest total IVI value (184.66), followed by Festuca ovina, Medicago minima, and Artemisia austriaca. The lowest IVI values across the entire study area were for Erodium hoefftianum, Cirsium sp, Malva nicaeensis, Astragalus stevenianus, and Linaria simplex.

Of the 251 species that were recorded, 44 species were class I, 55 species were class II, and 152 species were class III (Table 3). In class III, Artemisia austriaca, Colpodium parviflorum, and Bromus tectorum had the highest overall IVI values (108.35, 72.06, and 68.47, respectively), making them the most dominant unpalatable species. Agrostis canina, Sanguisorba minor, Geranium pusillum, Astragalus stevenianus, and Erodium hoefftianum were recorded as preferred plants (low IVI; Supplemental material,  Table S1 (mred-40-01-09_s01.pdf); The IVI of species class I declined with increasing elevation and then increased, while species class III showed the reverse trend (Figure 2).


Variation in IVI by palatability class along the elevation gradient of the Moghan Plain–Sabalan Mountain rangelands. Palatability classes: highly palatable (class I)—forage plants, most preferred by domestic livestock; mostly palatable (class II)—plants averagely preferred by livestock; hardly or unpalatable (class III)—plants least preferred or not grazed by livestock at any stage; and harmful and possibly toxic. ns indicates no significant difference.*Significant difference at P < 0.05. **Significant difference at P < 0.01.


Diversity index

Analyses of the diversity of species showed that diversity indices were significantly different across the elevation gradient (Figure 3). Diversity and species richness showed similar trends. They increased significantly with increasing elevation and then declined (Figure 4). There was no significant difference in species evenness along the elevation gradient (Figure 5).


Variation in species diversity by life form along the elevation gradient of the Moghan Plain–Sabalan Mountain rangelands. Life forms: He indicates hemicryptophyte; Th, therophyte; Ge, geophyte; Ch, chamaephyte. ns indicates no significant difference. *Significant difference at P < 0.05. **Significant difference at P < 0.01.



Variation in species richness by life form along the elevation gradient of the Moghan Plain–Sabalan Mountain rangelands. Life forms: He indicates hemicryptophyte; Th, therophyte; Ge, geophyte; Ch, chamaephyte. ns indicates no significant difference. *Significant difference at P < 0.05. **Significant difference at P < 0.01.



Variation in species evenness by life form along the elevation gradient of the Moghan Plain–Sabalan Mountain rangelands. Life forms: He indicates hemicryptophyte; Th, therophyte; Ge, geophyte; Ch, chamaephyte. ns indicates no significant difference. *Significant difference at P < 0.05. **Significant difference at P < 0.01.


The dominant species of hemicryptophytes with high IVI values were Festuca ovina and Colpodium parviflorum. Silene sp and Astragalus stevenianus were less abundant hemicryptophytes, with the lowest IVI in the study sites (Figure 6; Supplemental material,  Table S1 (mred-40-01-09_s01.pdf); A statistically significant second-order polynomial interrelation between species richness, evenness, and IVI with elevation was found, but the interrelation between diversity and elevation was not significant. Species richness and evenness were normally distributed along the elevation gradient, while IVI increased with increasing elevation and showed J-shaped interrelations with elevation. The diversity of hemicryptophytes increased with elevation but was not significantly different across the elevation gradient (Figures 35).


Variation in IVI by life form along the elevation gradient of the Moghan Plain–Sabalan Mountain rangelands. Life forms: He indicates hemicryptophyte; Th, therophyte; Ge, geophyte; Ch, chamaephyte. ns indicates no significant difference. *Significant difference at P < 0.05. **Significant difference at P < 0.01.


The most abundant therophytes of the study area were Trifolium subterraneum and Medicago minima, with high IVI. The rarest therophytes recorded with minimum IVI were Linaria simplex and Malva nicaeensis (Figure 6; Supplemental material,  Table S1 (mred-40-01-09_s01.pdf); The overall pattern of therophyte IVI shows an inverse interrelation with increasing elevation. A statistically significant second-order polynomial interrelation between diversity, species richness, and evenness with elevation was found for therophytes. Diversity and species richness were normally distributed along the elevation gradient (Figures 35).

The dominant species of chamaephytes with the highest IVI values were Artemisia austriaca and Thymus kotschyanus, whereas Helichrysum oligocephalum and Noaea mucronata were considered rare (Figure 6; Supplemental material,  Table S1 (mred-40-01-09_s01.pdf); The diversity and species richness of chamaephytes were normally distributed along the elevation gradient, while evenness and IVI showed U-shaped interrelations with elevation (Figures 35).

The dominant species of geophytes were Poa bulbosa and Colchicum trigynum, while the rare geophytes were Allium affine and Bongardia chrysogonum, which had the lowest IVI in the region (Figure 6; Supplemental material,  Table S1 (mred-40-01-09_s01.pdf); Species richness and evenness decreased with elevation, while diversity and IVI increased with increasing elevation; however, neither of these indicators was significantly different across the elevation gradient (Figures 35).


Interrelation between species composition and elevation gradient

The most species-rich families were Asteraceae, Poaceae, and Fabaceae. The dominance of these families is due to their high compatibility with the semiarid climate conditions (Ghafari et al 2018). Several studies (eg Teimoorzadeh et al 2015; Nazari Anbaran et al 2016; Ghafari et al 2017; Taheri Niari et al 2018) have also reported similar results in other habitats of Ardabil province. The high FIV values of Poaceae, Fabaceae, Asteraceae, and Caryophyllaceae might be due to their widespread dispersal mechanism (Muhammad et al 2016). Asteraceae, Caryophyllaceae, Fabaceae, Poaceae, and Scrophulariaceae have been able to adapt to different environmental conditions and resource limitations (Kargar Chigani et al 2017; Noroozi et al 2018), which results in their dominance in all selected sites. However, members of other families were restricted to a particular elevation, suggesting that specific physiographic and elevation conditions support different species. Asteraceae was the most species-rich family, showing high compatibility with the difficult conditions (shallow soil depth, steep slopes, and high solar radiation) in cold arid and semiarid climates (Ghafari et al 2017; Taheri Niari et al 2018).

Some species were mainly distributed at lower elevations, such as Pimpinella kotschyana, Lepidium draba, Phlomis herbaventi, and Senecio glaucus, whereas Ficaria kochii, Agrostis canina, and Lolium rigidum were recorded at the higher elevations. Plant species composition gradually changed as a function of variation in environmental factors; therefore, species distribution and composition differ across elevation classes along an ecological gradient (Rahman et al 2016). The species distribution showed that no single species has a wide distribution along the whole elevational range. Poa bulbosa (300–3000 masl) and Festuca ovina (900–3300 masl) were the only species that had a wide range of compatibility over the study area. The occurrence of these species across the elevation gradient suggests their wider ecological amplitude and tolerance to biotic pressures, in addition to their good reproductive capacity (Gairola et al 2011).

The IVI of species differed along the elevation, which indicates the relative importance of each species in the spectrum of elevation classes (Pandey 2015). Among the species with a high IVI, Poa bulbosa, Festuca ovina, Medicago minima, Artemisia austriaca, and Bromus tectorum spanned the widest elevational range. This indicates that they are the most compatible with available resources. Some species showed no distinctive distribution patterns along the gradient; some were distributed over a restricted elevational range, whereas others showed wider elevational ranges. The number of rare species in this study confirms the commonly cited concept that most of the species in plant communities are rare, rather than dominant (Magurran 2004). The rarity of species may be due to (1) strong density dependence (population growth rates are regulated by the density of a population); (2) uneven distribution of a resource, which leads to occupation by species at various locations resulting in variation in abundance and distribution; (3) natural or anthropogenic disturbances; and (4) competition effects (Scherrer and Körner 2011; Kacholi 2014; Körner et al 2017).

Interrelation of diversity, species richness, and evenness with elevation gradient

Species richness and diversity significantly increased with increasing elevation to 1200–1500 masl but then start to decline to 3300 masl. The increase in diversity value and species richness at mid-elevations could be because of (1) appropriate environmental conditions, such as low thermal and water stress and high humidity; (2) overlapping habitats; (3) the wide range (larger area) and proportionally more habitat area availability; and (4) potential high productivity (Ren et al 2012).

Due to the optimal combination of resources, more species coexist in the mid-elevation ranges (Ahmad Dar and Sundarapandian 2016); thus, more species were found in this elevation range. A positive impact of increased precipitation combined with the negative impact of decreased temperature produces a complex double gradient resulting in an optimum climate at the mid-elevation that favors diversity and species richness. While the limitations of low precipitation at low elevations or low temperature at high elevations would lead to different low-diversity ecosystems, dominated respectively by species adapted to limited precipitation or low temperatures (Sanchez et al 2013; Moradi et al 2017; Körner and Paulsen 2017).

Overlapping habitats and resources in mid-elevation areas could be partially responsible for the high diversity and species richness at this elevation range, which may indicate the existence of a transitional region for plants located at the mid-elevation (Senbeta et al 2014). Thus, in the mid-elevation a mix of high-elevation and low-elevation species give greater species richness than at higher or lower elevations (Ghafari et al 2018). The maximum diversity at mid-elevations might be attributed to the intermediate disturbance hypothesis. which suggests that the highest diversity will occur at levels of intermediate disturbance (Connell 1978; Mohamed Neji et al 2018). Most residential areas and agricultural lands are at low elevations, where there is a heavy disturbance. At higher elevations, degradation involves seasonal grazing and tourism activities; thus, disturbance is less than at low and mid-elevations. The reasoning is that diversity and species richness should be increased at median levels of disturbance (mid-elevation) because at high-disturbance levels (low elevation) most species will go locally extinct, and only the most resistant species will survive, whereas, at low levels of disturbance (high elevation), competitive exclusion increases and excludes other species (Zhang et al 2013). At the median disturbance, diversity increases because of disturbance-tolerant species and because highly competitive species coexist (Lee and Chun 2016).

Usually, with an increase of elevation, the total habitat decreases (McCain and Grytnes 2010; Noroozi and Körner 2018). Hence, the wider area in the lower elevations would provide a larger area for species distribution (Lomolino 2001). Though the area was larger at the lower elevation gradient, because of high temperature, low humidity, high human population density, changes in land use, road construction, overgrazing, and conversion of rangelands for agricultural purposes, only a few plant species are favored, with correspondingly decreased diversity and species richness (Nogués-Bravo et al 2008; Noroozi et al 2018).

The effect of fewer surfaces with increasing isolation results in a decreasing number of species at higher elevations (Lomolino 2001). In addition, the reduction in diversity and species richness at higher elevations could be related to ecophysiological constraints and harsh climatic conditions (extremely low temperature and strong winds; increase in snow cover and solar radiation), nutrient availability, water holding capacity, short growing seasons, geographical obstacles, and low ecosystem productivity (Mahdavi et al 2013; Sanchez et al 2013; Körner and Paulsen 2017). These restrict the physiological development of species to adapt and survive (Namgail et al 2012). The degree of disturbance decreased as elevation increased because tourism density and agriculture activities decreased with increased elevations (Sharifi et al 2017). However, heavy snowfall in winter, high grazing intensity in early spring and summer, and uncontrolled collection of medicinal and fodder plants has decreased species diversity (Sharifi et al 2017). With increasing elevation, people are increasingly dependent on rangelands because access to goods and services decreases (Thapa et al 2016). Sources of income and occupations other than livestock husbandry in the MPSMRs are rare, and, because of this, the number of livestock per unit area held by each household has increased (Sharifi et al 2017).

The evenness index showed no significant difference along the elevation gradient. Although the total evenness of all plants showed some fluctuation, there was no clear peak. According to Huang et al (1997), variations in the evenness of plant communities could be independent of the environmental gradient.

Hemicryptophytes were dominant (the highest IVI) at higher elevations. The change in dominance of life forms with elevation is probably due to the greater frequency of some species along the elevation gradient. For instance, the abundance of Festuca ovina, Colpodium parviflorum, Alopecurus vaginatus, Trifolium medium, and Poa pratensis resulted in a large number of hemicryptophytes at 2100–3300 masl elevation. Mota et al (2017) and Moradi and Attar (2019) found that hemicryptophytes are the most dominant life forms at high elevations. The lower disturbance at higher elevations favors hemicryptophytes. Thus, elevation and human activities play an important role in determining the diversity, life-form, and biological distribution spectrum of a habitat (Ren et al 2012).

The high proportion of therophytes in the investigated area reflects the cold winters and dry summers (Aghaei et al 2013). The spread of therophytes in the study area, especially unpalatable weeds, is likely a response to overgrazing (Balpinar et al 2018; Gamoun et al 2018). The preponderance of therophytes at mid-elevations might be the result of agricultural activities, overgrazing, and human activities (Ahmadi et al 2013). The dominance of therophytes at lower elevation can be partly explained by low soil moisture. At the lower elevation, although temperatures and growing season length are suitable, soil moisture is low (Barone et al 2008). Regardless of the type of ecosystem, overgrazing, especially at low elevations, led to the predominance of invasive and weed therophytes, indicating degradation (Carey et al 2017; Ghafari et al 2017; Wronski et al 2017; Lake and Minteer 2018; Shapero et al 2018).

The higher percentage of chamaephytes in the mid-elevation may be related to the impact of human activities (Balpinar et al 2018). Chamaephytes affect other species associated with stress tolerance and a strong competitive-response ability (Arila and Gupta 2016). Additionally, most chamaephytes are relatively resistant to livestock grazing (Terradas 1991) because of protective or defensive organs (such as spines, hairs, thorns, bitter taste, and unpleasant odor) and unpalatability. The species that provide the greatest protection to the buds are compact cushion plants, such as Acantholimon, Astragalus, and Onobrychis (Balpinar et al 2018; Moradi and Attar 2019). The higher elevation is characterized by higher solar radiation, which is unfavorable for the development of chamaephytes (Traquillini 1979).

The low proportion of geophytes in the region indicates that they are not adapted to the climatic conditions. The dominance of geophytes at higher elevations results from the high average annual rainfall, lengthy wet period during the growing season, and long duration of snow cover (Naqinezhad and Esmailpoor 2017). According to Meço et al (2018), geophytes are mostly associated with cold winters.

Plant palatability variation with elevation gradient

IVI species class I declined as the elevation increased (1800–2100 masl) and then increased, while species class III showed a reverse trend. At lower elevations, the rangeland conversion to agricultural land decreases rural livestock grazing intensity. Most of the rangeland is grazed by nomads' livestock. In addition, at low elevations, fodder crops, such as Medicago, Trifolium, and Onobrychis, are planted as livestock forage. Moreover, post-harvest crop residues are given to ranchers during the dry season for hay feeding, which provides part of the livestock needs in these areas. Therefore, in these regions, the presence of plants with good and medium palatability is higher than class III plants. In other words, although the diversity and species richness are low, rangeland degradation is medium at low elevation.

Optimal humidity conditions at mid-elevations and the high productivity in these habitats resulted in an optimal combination of resource availability (Hua 2004). Because of human population density, excessive grazing (rural and nomads' livestock), and conversion of rangelands for agricultural purposes, there is more degradation at mid-elevations. Although diversity and species richness increased, the presence of class III unpalatable plants was greater than the presence of plants with good and medium palatability at the mid-elevation. This is consistent with other studies that found that heavier grazing intensity was associated with higher species richness (Espinoza et al 2017). Moreover, many plants were found to be highly preferred, whereas some plants were avoided by the grazing livestock, leading these species to be widespread. This leads to the dominance of many unpalatable species, such as Astragalus microcephalus, Bromus arvensis, Bromus tectorum, Taeniatherum caput-medusae, and Artemisia austriaca as well as the disappearance of many other species that are less well adapted to unfavorable conditions and intense biotic influences. According to Renne and Tracy (2007), unpalatable plants increased in sites under grazing pressure. The major palatable species in the mid-elevation are Astragalus stevenianus, Astragalus pinetorum, Sanguisorba minor, Trifolium repens, and Trigonella arcuata, which are scarcely present at mid-elevations. Because they are highly palatable, these species have been greatly affected by grazing pressure. Most of the rare species present at mid-elevations were palatable and faced the excessive pressure of grazing. The grazing caused substitution of palatable plants by unpalatable plants, affecting abundance, dominance, plant frequency, and thus overall plant species composition (Espinoza et al 2017).

At higher elevations, the dominance of palatable species might be due to low human population density, lack of cultivation and agricultural activities, and limited seasonal grazing due to the shortened growing period (mid-May to mid-September). This rangeland is also only grazed by nomads' livestock. Nogués-Bravo et al (2008) explained that the higher elevations probably harbor large numbers of species because anthropogenic activities and destruction decrease with increased elevation.


Results showed that species diversity in plant communities of the region was considerable. Species richness and diversity peaked at mid-elevations along the elevational gradient. However, relative richness and diversity of species within different life-form groups varied differently along the gradient. The IVI of species differed along the elevation, indicating the relative importance of each species in a spectrum of elevation classes. The results showed that the analysis of plants using diversity indicators did not provide sufficient information on the actual condition of the rangeland plants. We believe the present approach should be encouraged in future studies for better understanding the structural status of rangelands using both species diversity and ecological importance indices of plants. By considering the results of this study, a proper conservation strategy, multipurpose use, and comprehensive regional planning could be developed to promote the sustainable use of MPSMRs.


This research was funded by the Grant and Research Program of the University of Mohaghegh Ardabili. Samaneh Mohammadi Moghadam, Zahra Abdolalizade, and Ruhollah Dabiri assisted with the field work.



Aghaei R, Alvaninejad S, Zolfaghari R, Mirzaei Gharehlar MR . 2013. Flora, life-form and geographical distribution of plants in west south forests of Iran (case study: Vezg, Yasouj). Journal of Farming and Allied Sciences 23(2):1046–1053. Google Scholar


Ahmad Dar J, Sundarapandian S . 2016. Patterns of plant diversity in seven temperate forest types of Western Himalaya, India. Journal of Asia-Pacific Biodiversity 9:280–292. Scholar


Ahmadi F, Mansory F, Maroofi H, Karimi K . 2013. Study of flora, life-form and chorotypes of the forest area of West Kurdistan (Iran). Bulletin of Environment, Pharmacology and Life Sciences 10(2):33–44. Google Scholar


Arila KE, Gupta A . 2016. Life-forms and biological spectrum along the elevation gradient in Montane forests of Senapati district of Manipur in Northeast India. Pleione 10(1):80–89. Google Scholar


Arjmand Karkazloo KA. 2016. The Effect of Different Grazing Gradient on Density and Production of Sage-Brush in Different Distances from Critical Spots in Moghan Plain [MSc dissertation] [in Persian]. Ardabili, Iran: University of Mohaghegh Ardabili. Google Scholar


Assadi M, Maassoumi AA, Khatamsaz M, Mozaffarian V . 1988–2012. Flora of Iran [in Persian]. Tehran, Iran: Research Institute of Forests and Rangelands Publication. Google Scholar


Balpinar N, Kavgaci A, Bingol Mu, Ketenoglu O . 2018. Diversity and gradients of vegetation of Sivrihisar Mountains (Eskişehir-Turkey). Acta Botanica Croatica 77(1):18–27.–0016Google Scholar


Barone JA, Thomlinson J, Cordero PA, Zimmermann JK . 2008. Metacommunity structure of tropical forest along an elevation gradient in Puerto Rico. Journal of Tropical Ecology 24(5):1–10. Scholar


Carey CJ, Blankinship JC, Eviner VT, Malmstrom CM, Hart SC . 2017. Invasive plants decrease microbial capacity to nitrify and denitrify compared to native California grassland communities. Biological Invasions 19(10):2941–2957. Scholar


Connell JH. 1978. Diversity in the tropical rainforest and coral reefs. Science 199:1301–1310. Google Scholar


Curtis JT, McIntosh RP . 1950. The interrelations of certain analytic and synthetic phytosociological characters. Ecology 31:434–455. Scholar


Davis PH. 1965–1988. Flora of Turkey and East Aegean Islands. Vols 1–10. Edinburgh, United Kingdom: Edinburgh University Press. Google Scholar


Espinoza JJO, Ayala CC, Castillón EE, Saldivar FG, Sauceda JU, Jurado E, Chapavargas L, Jaramillo EM, Hernández EO . 2017. Livestock effect on floristic composition and vegetation structure of two desert scrublands in northwest Coahuila, Mexico. Southwestern Naturalist 62(2):138–145. Scholar


FRWO [Forests, Rangeland and Watershed Organization]. 1982. Code of Rangeland Plant Species in Iran. Tehran, Iran: Research Institute of Forests and Rangelands Press. Google Scholar


Gairola S, Sharma CM, Ghildiyal SK, Suyal S . 2011. Tree species composition and diversity along an altitudinal gradient in moist tropical montane valley slopes of the Garhwal Himalaya, India. Forest Science and Technology 7(3):91–102. Scholar


Gamoun M, Belgacem O, Louhaichi M . 2018. Diversity of desert rangelands of Tunisia. Plant Diversity 40(5):217–225. Scholar


Ghafari S, Ghorbani A, Arjmand K, Teymorzadeh A, Hashemi Majd K, Jafari S, Dabiri R . 2017. Effect of grazing intensity on vegetation and soil physiochemical properties (case study: in rangelands of Kolash village, Parsabad, Ardabil province) [in Persian]. Journal of Plant Ecosystem Conservation 10:183–204. Google Scholar


Ghafari S, Ghorbani A, Moameri M, Mostafazadeh R, Bidarlord M . 2018. Composition and structure of species along elevation gradient in Moghan-Sabalan rangelands, Iran. Journal of Mountain Science 15(6):1209–1228. Scholar


Ghahreman A. 1979–2003. Color Flora of Iran. Vols 1–26 [in Persian]. Tehran, Iran: Research Institute of Forests and Rangelands Press. Google Scholar


Ghorbani A, Ahmadalei V, Asghari A . 2014. Study the effect of distance from village on plant diversity and composition in rangeland of Southeastern Sabalan [in Persian]. Journal of Rangeland 8(2):178–191. Google Scholar


Ghorbani A, Ghafari S, Sattarian A, Akbarlou M, Bidar Lord M . 2017. Medicinal plants of Sabalan rangeland ecosystem in Ardabil province [in Persian]. Journal of Plant Ecosystem Conservation 9(4):77–96. Google Scholar


Ghorbani A, Mohammadi Moghaddam S, Hashemi Majd K, Dadgar N . 2018. Spatial variation analysis of soil properties using spatial statistics: A case study in the region of Sabalan mountain, Iran. eco.mont 10(1):60–70. Scholar


Hosseinzadeh A, Heidari GA, Barani H, Zali H . 2017a. Effects of beneficiaries' social issues on rangeland ecological sustainability (Case study: Shahsavan nomads of Meshginshahr city) [in Persian]. Journal of Rangeland 10(4):465–473. Google Scholar


Hosseinzadeh A, Heidari GA, Barani H, Zali H, Zandi Esfahan E . 2017b. Assessing the effects of socio-economic issues on social cohesion of nomads (case study: Shahsavan nomads, Sabalan Mountain Range) [in Persian]. Journal of Range and Watershed Management 69(4):851–861. Google Scholar


Hua Y. 2004. Distribution of Plant Species Richness Along Elevation Gradient in Hubei Province, China. Nanjing, China: International Institute for Earth System Science, Nanjing University. Google Scholar


Huang JH, Gao XM, Ma KP, Chen LZ . 1997. A comparative study on species diversity in zonal forest communities. Acta Ecologica Sinica 17(6):611–618. Google Scholar


Hussain F, Durrani MJ . 2009. Seasonal availability, palatability and animal preferences of forage plants in Harboi arid ange land, Kalat, Pakistan. Pakistan Journal of Botany 41(2):539–554. Google Scholar


IPNI [The International Plant Names Index]. 2018. The International Plant Names Index. London, United Kingdom: Royal Botanic Gardens, Kew, Harvard University Herbaria, and Australian National Herbarium.; accessed on 2 May 2018. Google Scholar


Kacholi DS. 2014. Analysis of structure and diversity of the Kilengwe Forest in the Morogoro Region, Tanzania. International Journal of Biodiversity 48:1–8. Scholar


Kakehmami A, Ghorbani A, Kayvan BF, Mirzaei MA . 2017. Comparison of visual and digital interpretation methods of land use/cover mapping in Ardabil province. Journal of RS and GIS for Natural Resources 8(3):121–134. Google Scholar


Kargar Chigani H, Akbarjavadi S, Zahedi Amiri G, Jafari M, Khajeddin JS . 2017. The floristic composition and biological spectrum of vegetation in the Meymeh region of Northern Isfahan province, Iran. Applied Ecology and Environmental Research 15(1):415–428. Google Scholar


Komarov VL. 1934–1954. Flora of USSR. Vols 1–21. Botanicheskii institut (Akademiia nauk SSSR). Mifal tirgume ha-mada ha-Yisreeli. Jerusalem, Israel: Israel Program for Scientific Translations. Google Scholar


Körner C, Jetz W, Paulsen J, Payne D, Rudmann-Maurer KM, Spehn E . 2017. A global inventory of mountains for bio-geographical applications. Alpine Botany 127(1):1–15. Scholar


Körner C, Paulsen JA . 2017. Geostatistical and bioclimatological comparison of the Central Great Caucasus and the Central Alps. In : Nakhutsrisvili G, Abdaladze O, Batsatsashvili K, Spehn EM, Körner C, editprs. Plant Diversity in the Central Great Caucasus: A Quantitative Assessment. Cham, Switzerland: Springer, pp 1–9. Scholar


Körner C, Paulsen J, Spehn EM . 2011. A definition of mountains and their bioclimatic belts for global comparisons of biodiversity data. Alpine Botany 121(2):73–78. Scholar


Lake EC, Minteer CR . 2018. A review of the integration of classical biological control with other techniques to manage invasive weeds in natural areas and rangelands. BioControl 63(1):71–86. Scholar


Lee CB, Chun JH . 2016. Environmental drivers of patterns of plant diversity along a wide environmental gradient in Korean temperate forests. Forests 7:3–19. Scholar


Lomolino MV. 2001. Elevation gradients of species-density: Historical and prospective views. Global Ecology and Biogeography 10:3–13. Scholar


Magurran AE. 2004. Measuring Biological Diversity. Oxford, United Kingdom: Wiley. Google Scholar


Mahdavi P, Akhani H, Van der Maarel E . 2013. Species diversity and life-form patterns in steppe vegetation along a 3000 m altitudinal gradient in the Alborz Mountains, Iran. Folia Geobotanica 48:7–22. Scholar


McCain CM, Grytnes JA . 2010. Elevational gradients in species richness. In : Encyclopedia of Life Sciences. Chichester, United Kingdom: John Wiley & Sons. Scholar


Meço M, Mullaj A, Barina Z . 2018. The vascular flora of the Valamara mountain range (SE Albania), with three new records for the Albanian flora. Flora Mediterranea 28:5–20. Scholar


Menhinick EF. 1964. A comparison of some species–individuals diversity indices applied to samples of field insects. Ecology 45(4):859–861. Scholar


Misra R. 1968. Ecology Work Book. Calcutta, India: Oxford and IBH. Google Scholar


Mohamed Neji CD, Serbajib MM, Hardyc O, Chaieba M . 2018. Floristic diversity and vegetation patterns along disturbance gradient in arid coasts in southern Mediterranean: Case of the Gulf of Gabès, southern Tunisia. Arid land Research and Management 32(3):291–315. Scholar


Moradi H, Attar F . 2019. Comparative study of floristic diversity along altitude in the northern slope of the central Alborz Mountains, Iran. Biodiversitas 20(1):305–312. Scholar


Moradi H. Attar F. Oldeland J. 2017. Plant functional type approach for a functional interpretation of altitudinal vegetation zones in the Alborz Mts., Iran. Journal of Mountain Science 14:2257–2269. Scholar


Mori SA, Boom BM, de Carvalho AM, dos Santos TS . 1983. Southern Bahian moist forests. Botanical Review 49(2):155–232. Scholar


Mota GS, Luz GR, Mota NM, Coutinho E, Veloso MDM, Fernandes GW, Nunes YRF . 2017. Changes in species composition, vegetation structure, and life-forms along an altitudinal gradient of rupestrian grasslands in southeastern Brazil. Flora 238 (1):32–42. Scholar


Muhammad Z, Khan N, Ali Sh, Ullah A, Khan ShM. 2016. Density and taxonomic diversity of understory vegetation in relation to site conditions in natural stands of Acacia modesta in Malakand Division, Khyber Pakhtunkhwa, Pakistan. Science, Technology and Development 35(1):26–34. Google Scholar


Namgail T, Rawat GS, Mishra C, Evan Wieren S, Prins HHT . 2012. Biomass and diversity of dry alpine plant communities along altitudinal gradient in the Himalayas. Journal of Plant Research 125:93–101. Scholar


Naqinezhad A, Esmailpoor A . 2017. Flora and vegetation of rocky outcrops/cliffs near the Hyrcanian forest timberline in the Mazandaran mountains, northern Iran. Nordic Journal of Botany 35(4):449–466. Scholar


Nazari Anbaran F, Ghorbani A, Azimi Motem F, Teymorzadeh A, Asghari A, Hashemimajd K . 2016. Floristic and species diversity in altitudinal gradient of Lahrod-Shabil (North Sabalan). [in Persian]. Journal of Plant Ecosystem Conservation 7(3):1–18. Google Scholar


Nogués-Bravo D, Araújo MB, Romdal T, Rahbek C . 2008. Scale effects and human impact on the elevational species richness gradients. Nature 453:216–220. Google Scholar


Noroozi J, Körner C . 2018. A bioclimatic characterization of high elevation habitats in the Alborz mountains of Iran. Alpine Botany 128(1):1–11. Scholar


Noroozi J. Talebi A. Doostmohammadi M. Rumpf SB. Linder HP. Schneeweiss GM. 2018. Hotspots within a global biodiversity hotspot: Areas of endemism are associated with high mountain ranges. Scientific Reports 8:10345. Scholar


Pandey KP. 2015. Assessment of Structure, Composition and Carbon Stocks of Forests Along the Elevational Gradient in the Langtang National Park, Nepal [MSc dissertation]. Munich, Germany: Technical University of Munich. Google Scholar


Pielou EC. 1969. An Introduction to Mathematical Ecology. New York, NY: Wiley. Google Scholar


Plant List. 2018. The Plant List.; accessed on 2 May 2018. Google Scholar


Rahman AU, Khan ShM, Khan S, Hussain A, Rahman IU, Iqbal Z, Ijaz F . 2016. Ecological assessment of plant communities and associated edaphic and topographic variables in the Peochar Valley of the Hindu Kush Mountains. Mountain Research and Development 36(3):332–341. Scholar


Raunkiaer C. 1934. The Life-Forms of Plants and Statistical Plant Geography. Oxford, United Kingdom: Clarendon Press. Google Scholar


Rechinger KH. 1963–2010. Flora Iranica. Vols 1–178. Graz, Austria: Akademische Druck-u Verlagsanstalt. Google Scholar


Ren X, Yang G, Zhu F, Qin X, Wang D, Liu Z, Feng Y . 2012. Plant communities, species richness and life-forms along elevational gradients in Taibai Mountain, China. African Journal of Agricultural Research 7(12):1834–1848. Scholar


Renne IJ, Tracy BF . 2007. Disturbance persistence in managed grasslands: Shifts in aboveground community structure and the weed seed bank. Plant Ecology 190:71–80. Scholar


Sanchez M, Pedronia F, Eisenlohrb PV, Oliveira-Filhob AT . 2013. Changes in tree community composition and structure of Atlantic rain forest on a slope of the Serra do Mar range, southeastern Brazil, from near sea level to 1000 m of elevation. Flora 208:184–196. Scholar


Scherrer D, Körner C . 2011. Topographically controlled thermal-habitat differentiation buffers alpine plant diversity against climate warming. Journal of Biogeography 38(2):406–416. Scholar


Senbeta F, Schmitt Ch, Woldemariam T, Boehmer HJ, Denich M . 2014. Plant diversity, vegetation structure and relationship between plant communities and environmental variables in the Afromontane forests of Ethiopia. Ethiopian Journal of Science 37(2):113–130. Google Scholar


Shannon CE, Weaver W . 1963. The Mathematical Theory of Communication. Urbana, IL: University of Illinois Press. Google Scholar


Shapero MW, Huntsinger L, Becchetti TA, Mashiri FE, James JJ . 2018. Land manager perceptions of opportunities and constraints of using livestock to manage invasive plants. Rangeland Ecology and Management 71(5):603–611. Scholar


Sharifi J, Rostami Kia Y, Azimi F, Ashoori P . 2017. Ecological Regions of Iran Vegetation Types of Ardabil Province [in Persian]. Tehran, Iran: Research Institute of Forests and Rangelands Press. Google Scholar


Taheri Niari MM, Ghorbani A, Bidar Lord M, Ghaffari S . 2018. Study the structural of flora composition in elevation gradient of QezelOzan - Kosars' altitudes rangelands, Ardabil province [in Persian]. Iranian Journal of Taxonomy and Biosystematics 34:53–72. Google Scholar


Tavosi T, Delara Gh . 2011. Climatic zoning of Ardabil Provinces [in Persian]. Journal of Nivar 70–71:47–52. Google Scholar


Teimoorzadeh A, Ghorbani A, Kavianpoor AH . 2015. Study on the flora, life forms and chorology of the south eastern of Namin forests (Asi-Gheran, Fandoghloo, Hasani and Bobini), Ardabil province [in Persian]. Journal of Plant Research (Iranian Journal of Biology) 28(2):264–275. Google Scholar


Terradas J. 1991. Mediterranean woody plant growth-forms, biomass and production in the eastern part of the Iberian Peninsula. Oecologia Aquatica 10:337–349. Google Scholar


Thapa S, All J, Yadav RKP . 2016. Effects of livestock grazing in pastures in the Manaslu Conservation Area, Nepalese Himalaya. Mountain Research and Development 36(3):311–319.  Google Scholar


Traquillini W. 1979. Growth of trees at timberline. In : Traquillini W. Physiological Ecology of the Alpine Timberline. Ecological Studies Vol 31. Berlin, Germany: Springer, pp 19–38. Scholar


Whittaker RH. 1953. A consideration of climax theory: The climax as a population and pattern. Ecological Monographs 23(1):41–78. Scholar


Wronski T, Bariyanga J, Sun P, Plath M, Apio A . 2017. Pastoralism versus agriculturalism-how do altered land-use forms affect the spread of invasive plants in the degraded Mutara rangelands of north-eastern Rwanda?. Plants 6(2):2–19. Scholar


Zhang JT, Xu B, Li M . 2013. Vegetation patterns and species diversity along elevational and disturbance gradients in the Baihua Mountain Reserve, Beijing, China. Mountain Research and Development 33(2):170–178. Scholar


Supplemental material

 TABLE S1 (mred-40-01-09_s01.pdf) Importance value index (IVI) along elevation gradient in the Moghan Plain–Sabalan Mountain rangelands.

Found at:

© 2020 Ghafari et al. This open access article is licensed under a Creative Commons Attribution 4.0 International License ( Please credit the authors and the full source.
Sahar Ghafari, Ardavan Ghorbani, Mehdi Moameri, Raoof Mostafazadeh, Mahmood Bidarlord, and Azad Kakehmami "Floristic Diversity and Distribution Patterns Along an Elevational Gradient in the Northern Part of the Ardabil Province Rangelands, Iran," Mountain Research and Development 40(1), R37-R47, (16 December 2020).
Received: 1 September 2019; Accepted: 1 May 2020; Published: 16 December 2020

Back to Top