Next Article in Journal
Diversity, Biogeography and Community Ecology of Ants: Introduction to the Special Issue
Next Article in Special Issue
Biodiversity and Conservation of Vascular Flora: A Challenge for the Survival of Humans on Earth
Previous Article in Journal
Landscape Genetics of the Yellow-Bellied Toad (Bombina variegata) in the Northern Weser Hills of Germany
Previous Article in Special Issue
Updated Checklist of Vascular Plants Endemic to Mongolia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring Effective Conservation of Charismatic Flora: Orchids in Armenia as a Case Study

1
Royal Botanic Gardens Kew, Millennium Seed Bank Building, Ardingly RH17 6TN, UK
2
A. Takhtajyan Institute of Botany of National Academy of Science of Armenia, Acharyan 1, Yerevan 0040, Armenia
3
Nature Heritage NGO of the Republic of Armenia, Yerevan 0040, Armenia
*
Author to whom correspondence should be addressed.
Submission received: 30 September 2021 / Revised: 4 November 2021 / Accepted: 23 November 2021 / Published: 28 November 2021
(This article belongs to the Special Issue Biodiversity and Conservation of Vascular Flora)

Abstract

:
Plants are the bedrock of life on Earth but are increasingly threatened with extinction. The most cost-effective way of conserving plant diversity is through Protected Areas (PAs). However, the locality, size, and management of PAs are crucial for effectively maintaining diversity and have been criticized as currently inadequate. Using Armenia as our study site and orchids as our study taxa, we sought to (1) identify spatial patterns of orchid diversity hotspots and corresponding PA network sites; (2) examine if the current PA network is effective at capturing orchid species richness and diversity and (3) explore the relationship between the range of area suitability of species and level of protection. We used data collected from herbarium, field visits and GBIF occurrence records. Using freely available mapping software, we created heatmaps of observations and species richness. We compared PA sites based on the number of species (species richness) and diversity (Shannon–Weiner Index). Species range was developed using the MaxEnt model and a correlation analysis was performed against the proportion of the range within PA. We found that 57% of PA sites have a representation of at least one species of orchid, but some threatened species are not presented within any PA site. The Tavush and Syunik province not only held the highest species richness (>10 species), but the PA network within also held high orchid diversity (2.5 diversity index value for Dilijan National Park). We did not find a significant relationship between the range of area suitability for orchids and protection; however, all our target species had less than 30% of their range under protection. Our study highlights important challenges related to the limitations of available data, and we discuss these implications towards effective conservation outcomes for orchids for the region.

1. Introduction

The urgency of protecting our world’s biodiversity is now becoming apparent, with various landmark reports attracting public attention [1,2]. Key to this is the growing understanding that biodiversity is fundamental to all life on Earth and that anthropogenic impacts are the largest contributor to its decline [1,2,3,4]. The impact of our economic growth within the last 50 years has certainly contributed towards an estimated 1 million of all described species being threatened with extinction [2]. Plant and fungal diversity are the bedrock of healthy ecosystems and are worryingly showing the same downward trend across different continents [5,6,7]; globally, 39% of all vascular plants are now threatened with extinction [8].
Arguably, the best way of curbing biodiversity decline and protecting wild species is the establishment of Protected Areas (PAs), as these provide a safe refuge for both plants and animals within their natural habitats [9,10]. International biodiversity initiatives have consistently featured various targets that relate to a percentage of land designated for protection (e.g., Aichi Biodiversity Targets [11]). Although nations have met, or, in some cases, exceeded these targets, the continued decline in biodiversity over the past few decades has cast some skepticism over their overall effectiveness [12,13,14,15]. Historically, the reasoning behind the establishment of PAs focused on areas that were remote or “wild”, their aesthetics and/or sole use of natural resources [16]. As a result, many PAs were either not adequately positioned to capture species location [12,16,17], conserving genetic diversity [18] and/or ineffectively managed [19]. In recent years, the establishment of PAs has become more strategic—for example, the Natura 2000 sites in Europe [20]. Additionally, the use of tools such as Species Distribution Models (SDMs) is starting to feature in mainstream decision-making processes, ensuring that the development of PAs is based on ecologically sound evidence [21]. Although the use of more strategic and scientifically robust methods in determining PAs is welcomed, there is a bias towards areas important for megafaunas and “umbrella species”, which can result in more cryptic species being left behind [16]. Additionally, there is increasing awareness of the phenomenon known as “plant blindness”, where plants are generally overlooked by policy-makers, conservationists and the general public [22]. As with other less well-known or well-loved species, the limited or lack of attention given to plants in conservation action can exacerbate their rate of decline, thereby risking further ecosystem collapse [22].
Some of the most charismatic taxa within the plant kingdom are those within the family Orchidaceae, the second most diverse angiosperm plant family after Asteraceae, with an estimated 28,484 species worldwide [23]. The popularity in the study of orchids, both by academics and amateurs, could be attributed to various reasons, from simply the beauty of its flowers, the curious co-evolution with pollinator systems [24,25,26], the remarkable cases of sexual deception in some genera such as Ophrys L. [27,28] and/or the complex and, at times, very specific ecology and relationship with mycorrhizal fungi [29,30]. It is their complex ecological specialization that makes orchids an ideal taxon for testing the effectiveness of various conservation strategies, as they can be used as bioindicators of ecosystem health [31].
In Europe, all orchids are of the terrestrial form, with the mountainous regions of modern Europe, Asia Minor, Middle and Central Asia and the Caucasus recognized as important areas for some taxa (e.g., Epipactis and Dactylorhiza) [32]. Despite its popularity, there are major gaps in knowledge for orchids, particularly in parts of West Asia [33], leaving species within these biodiversity-rich areas exposed to increasing threats of extinction. Although some species are found to have a broad range, stretching across Europe (e.g., Anacamptis pyramidalis (L.) Rich. [34]), others have evolved such specialist relationships with their immediate environment, restricting their range to only a small area (e.g., Goodyera macrophylla Lowe [35]). With increasing levels of habitat loss, fragmentation and the ongoing effects of a changing climate, terrestrial orchids with narrow ranges are likely to experience a greater risk of extinction in the coming decades [31,36]. Therefore, a substantial proportion of these narrow endemics will need to be accounted for within PA networks to enable the conservation of the entire species to be effective.
The Caucasus Biodiversity Hotspot holds c. 7000 species of vascular plants, of which 25% are classified as endemic to the region [37]. Within this biodiverse region is Armenia, a small landlocked country within the South Caucasus. Although the country is relatively small (territorial area of 29,743 km2), more than 3800 vascular plant species are found within, making up more than a half of the Caucasus flora. Additionally, there are elements of both horizontal and vertical habitat level zoning found here, resulting in a complex landscape structure with a wide range of different habitat types, from semi-desert zones to alpine grasslands and moist temperate forests. Therefore, we would expect that the PAs within Armenia would encapsulate a wider range of different bioclimatic and ecological types across its network. In the last decade, there has been a significant expansion of PA networks across the Caucasus, with Armenia incorporating Lake Arpi National Park [38], Arevik National Park [39], Khustup Sanctuary [40] and the Zangezur Sanctuary [40] within its network. The current PA network within Armenia consists of 34 sites, the majority of which are recognized as nationally designated Protected Areas within the IUCN Category IV (27 sites), four under IUCN Category II and three under IUCN Category 1a. There are also two internationally recognized Ramsar sites, namely Lake Arpi to the north of the country and the Khor Virap marshlands within the Ararat Province [41]. Although the continued expansion of the protected network is welcomed, there is an acknowledgement that gaps are present. Indeed, the new PAs were incorporated with the aim of enhancing the protection of the endangered Persian leopard, and there has been little attention towards key plant groups or families within the prioritization process [37].
In this paper, we aim to explore whether the current establishment of PAs is effective at conserving orchid diversity. We highlight hotspots for orchids and identify the richness and diversity within current PA networks. We hypothesize that orchid species suited to thriving within a narrow range of areas, i.e., higher extinction risk, will be underrepresented within the current PA network and therefore should be prioritized for further conservation action.

2. Materials and Methods

2.1. Data Collection and Cleaning

Data on species name, locality (GPS, locality notes and/or descriptions), collection dates, altitude and, where possible, collector and verifier were gathered from herbarium specimens of orchids deposited in the Herbarium of the A. Takhtajyan Institute of Botany of NAS RA (ERE) and the Komarov Botanical Institute of the Russian Academy of Sciences, St. Petersburg Herbarium (LE) (herbaria acronyms according to [42]). We used the QGIS GBIF plugin [43] to download data on orchids from the Global Biodiversity Information Facility (GBIF) occurrence database [44]. We used data (taxonomy, locality, altitude and data source) derived from preserved specimens (i.e., herbarium vouchers) and human observations. Species that were found only within the GBIF human observation dataset were removed to reduce the possibility of misidentification. Additional fieldwork between the months of May and June 2020 was conducted to confirm orchid locality and identify new populations. Ecology and flowering times were obtained through herbarium studies and the published literature [45,46]. The nomenclature of the resulting species list of orchids was validated through name matching using the World Checklist of Selected Plant Families (WCSP) [23] and the World Checklist of Vascular Plants (WCVP) [47].
For herbarium collections, labels of 931 herbarium specimens, representing 39 species and 16 genera, were initially transcribed, and by using the published literature on the habitat preferences of each species, and location description, geographic coordinates for each specimen were obtained through Google Earth (version 2020). A total of 58 vouchers that were collected between 1929 and 1980 were omitted from the subsequent analysis as their locality could not be verified. For field observations, we omitted records where the field observers could not readily identify the plant to species due to lack of flowers during the time of the visit. Finally, for the GBIF records, we limited our investigation to those identified to species.
The Protected Areas (PAs) shapefile was downloaded from the World Database of Protected Areas (WDPA) through the Protected Planet website [41] and cropped for the Armenia territory shapefile obtained from the freely available resource, the DIVA-GIS website [48].

2.2. Sampling Density and Species Richness Analysis

We used a freely available mapping software, QGIS (version 3.10.12-A Coruña), to project the 1101 geocoordinates from the resulting orchid dataset. A heatmap of the collection data was produced using the Heatmap Tool to identify the pattern of orchid distribution according to presence-only data. The tool uses Kernel Density Estimation with a radial distance buffer of 5 km. Where 10 or more points and their radial buffer intersect, these would be identified as a cluster. We overlaid the PA network shapefile with the point coordinates of orchid species and, using the “point sampling” tool in QGIS, extracted orchid occurrence points within each PA site to identify which PA was better at capturing orchid species diversity using the Shannon–Weiner diversity index found within the “vegan” package [49] in R software 4.0.2 [50].
We identified areas of high species richness across the study area by overlaying 5 km × 5 km grids across the study area and using a readily available algorithm within the QGIS toolbox (“Count points in polygon”), and we counted the number of unique species within each grid cell. We specified the accepted taxonomic names for each species as the unique class field.

2.3. Species Distribution Models (SDMs) of Conservation Priority Species

We used Species Distribution Models (SDMs) to develop the potential areas suitable for orchid species found in Armenia. We used a maximum entropy (MaxEnt) approach to model orchid species found within our study system using the MAXENT software (version 3.4.3) [51], which is the preferred method for presence-only data (e.g., herbarium specimens, etc.), and handles small datasets relatively well [52]. SDMs were done for species that had 10 occurrence data or over.
Previous studies have found various factors that can influence orchid distribution across a landscape, including pollination specificity and mycorrhizal association [31]. As both factors are also found to correlate with altitude [53,54,55,56], we used this as proxy variables within the SDMs. Alongside altitudinal data, we also retrieved climatic variables relating to precipitation and temperature from the World Clim dataset (worldclim.org, 2020). These were clipped to the Armenia territory outline and bioclimatic variables from this monthly climate data were extracted using the biovars function within the “dismo” package [57] within R statistical software [50]. Environmental variables of orchid localities were extracted using the “Point Sampling Tool” on QGIS. A Principal Component Analysis (PCA) was used to identify correlations across variables using R software. To reduce the risk of overfitting our models, we took the variables strongly associated with the first four Principal Component axes as our environmental variables for the subsequent models.
The MAXENT software generated 10,000 random points within our study area as “absence” points and developed predicted distribution using the Jacknifing method. For model validation, we used 25% of the presence-only data and cloglog as the output format, as recommended by [58]. We used linear, quadratic and hinge features to run the models using the auto feature turned on within the MAXENT software. To identify the model with the most appropriate regularization parameter, we ran models with regularization multiplier values from 1 to 10 for each species and compared the models using the Akaike’s information criterion [59] corrected for small sample size (AICc), following the methodology based on Warren and Seifert [60]. The model with the lowest AICc value was taken as the one to move forward with in subsequent analysis as this represented the model that was the most parsimonious.
We calculated area suitability ranges with high predicted conditions by converting the output raster with prediction values above 0.7 into a vector and calculating this area using the “Add vector geometry” tool in QGIS, which gives the area in km2 per pixel. The mean range size for each species within Armenia was calculated. An overlap analysis was done with the PA shapefile layer to extract the percentage of potential species distribution within PA networks. A correlation analysis between mean potential species range and mean percentage of range protected was done using R software, and Pearson’s correlation coefficient was calculated to test the null hypothesis outlined in the Introduction.

3. Results

3.1. Sampling Density and Overall Orchid Diversity

The oldest specimen found within the ERE herbarium was collected in 1916 by B. Shishkin and N. Abzianidze and identified as Orchis palustris Jacq., later recognized as a synonym of Anacamptis palustris (Jacq.) R. M. Bateman, Pridgeon & M. W. Chase [23]. The most recent herbarium specimen used for the study was Anacamptis morio (L.) R. M. Bateman, Pridgeon & M. W. Chase (syn Orchis morio (L.)), collected in 2017 by A. Nersesyan, A. Papikyan, S. Galstyan and N. Hayrapetyan near Goris, a town in the southeast of Armenia. Based on herbarium vouchers and the published literature, the earliest flowering orchid was A. pyramidalis (April), while the latest flowering orchids were Epipactis persica (Soó) Hausskn. ex Nannf. and Epipogium aphyllum Sw. (August). Generally, orchids found in Armenia tended to bloom between May and July [45].
After data cleaning and accounting for synonyms, an initial list of 48 species was reduced to a total of 43 accepted species and subspecies from 14 genera of orchids found in Armenia, with 18 identified as synonyms according to the WCSP and WCVP (Supplementary Materials Table S1 Orchid species occurrence data). Dactylorhiza maculata subsp. fuchsii (Druce) Hyl., D. maculata (L.) Soó. and Orchis x calliantha Renz & Taubenheim records were omitted from subsequent analysis as these were the only species found in the GBIF human observation dataset and not verified elsewhere. There were 1101 occurrence data captured from the resulting list of species from the various sources. We found 17 clusters, where 10 or more points and their 5 km radial buffer overlapped. The majority were found to the south of Armenia within Syunik province, which also had the largest cluster, followed by Tavush province, towards the north-east of the country (Figure 1a). We also found a small cluster within Lori province, close to the Gyulagarak PA site, and in the northern part of Sevan Lake, within Gegharkunik province. A similar pattern was found in relation to species richness, with Syunik having the highest number of species (16) within a 5 km-by-5-km grid, followed by Tavush province with 12 species and the site in Lori with eight species within a 5 km-by-5-km grid square (Figure 1b). In the Ararat and northern Kotayk provinces, we found small observational clusters, but with relatively low species numbers (2 to 5 species) (Figure 1a,b).

3.2. Protected Area Network Comparisons

Out of the current PA network in Armenia, 57% of the sites held at least one species of orchid (20 out of 35) (Table S2 Species presence absence within PA network in Armenia). The Dilijan National Park, classified as Category II under the IUCN Protected Area Categories System, had the highest orchid diversity (number and evenness) at 2.498 (Table 1). PA sites that were categories under strict protection (Category Ia) had variable richness and diversity, with Shikahogh State Reserve holding 13 species and a diversity index of 2.26, compared to Khosrov Forest State Reserve, holding only six species and a diversity index of 1.50, despite having a larger area (Table 1). Moreover, 14 out of the 43 orchid species (33%) were not found to be represented within any of the protected sites (Table S2 of species presence within PA network in Armenia), including Steveniella satyrioides (Spreng.) Schltr., a species listed as EN under the Armenian Red Data Book of Plants [61]. Eight species were only found within one PA site across the entire network, including two species also included within the Armenian Red Data Book of Plants [61], Cephalanthera kurdica Bornm. ex Kraenzl. (VU*), from herbarium records and Corallorhiza trifida Chatel. (CR) from GBIF (human observation) records. In the case of C. kurdica, all known observations for this species were only within the territory of the Arevik National Park.

3.3. Potential Areas of Suitability and Protection

Out of the 43 species of orchids found, only 21 had sufficient data points for modeling. PCA analysis of the bioclimatic variables selected showed that 93% of the variation could be explained by the first four Principal Components, with the first (PC1) and second (PC2) Principal Components explaining 77.4% of the overall variation (Figure S2 PCA plot). The maximum temperature of the warmest month (Bio5) and precipitation of the warmest quarter (Bio18) corresponded to PC1 and PC2, respectively, while precipitation of the driest quarter (Bio17) and altitude (Alt) corresponded to PC3 and PC4, respectively (Table S3 Principal Component Analysis results).
Test AUC scores (i.e., the real test of the model’s predictive power) for the models varied from relatively low (i.e., model performs close to a random model), in the case of Epipactis persica (Test AUC = 0.572), Platanthera chlorantha (Custer) Rchb. (0.590) and Gymnadenia conopsea (L.) R. Br. (0.597), to those with high predictive power, in the case of Limodorum abortivum (L.) Sw. (0.908), Cephalanthera damasonium Mill. (0.866) and Anacamptis pyramidalis (0.844) (Table 2). Species-specific predicted distributions and response curves can be found in the Supplementary Materials (Figure S3 SDM models per species). Maximum temperature of the warmest month (Bio5) was the best predictor for six of the 21 species, namely for A. coriophora (L.) R. M. Bateman, Pridgeon & M. W. Chase and Orchis punctulata Steven ex Lindl. Precipitation of the warmest quarter (Bio18) best predicted five species, namely P. chlorantha and Neottia nidus-avis (L.) Rich., and with regard to Traunsteinera sphaerica (M.Bieb.) Schltr., both Bio5 and Bio18 contributed roughly 50/50 to its distribution (Table 2). Precipitation of the driest quarter (Bio17) did not rank highly for any of the species, although for Dactylorhiza urvilleana (Steud.) H.Baumann & Künkele, together with Bio5, the variables contributed over 80% of the model’s variation, with similar sigmoidal response curves (Figure S3 Species-specific distribution models from MaxEnt). Altitude appeared to be the best predictor for 10 out of the 21 species, although this variable contributed the least within the PCA (Table 2).
According to the MaxEnt models generated, species that had the highest average range of environmentally suitable area were Epipactis helleborine (L.) Crantz (1070 km2), followed by Neottia nidus-avis (933 km2) and Dactylorhiza urvilleana (793 km2). Those with the narrowest area of suitability area were Neotinea tridentata (Scop.) R. M. Bateman, Pridgeon & M. W. Chase (192 km2), followed by Traunsteinera sphaerica (221 km2) and Gymnadenia conopsea (238 km2) (Figure 2). N. tridentata is currently an accepted name for two synonyms of nationally threatened species, Orchis simia Vill. (NT) and O. tridentata Scop. (EN) (Table S1 Orchid species occurrence data).
We found no significant correlation between the average potential range of suitable area versus proportion of said area under protection (Pearson’s correlation: r = −0.211, p-value = 0.358). None of the orchid species found had more than 30% of its potential area within the PA network (Figure 2). The two species with the largest potential distribution varied in their levels of protection, with E. helleborine having, on average, 1.65% of its range under protection, but D. urvilleana showing a higher proportion of its range currently under protection (21.9%) (Figure 2). We found several species with narrow suitability range showing similar protection values to D. urvilleana, namely N. tridentata (19.2%), Anacamptis coriophora (22.8%) and D. viridis (L.) R. M. Bateman, Pridgeon & M. W. Chase (21.3%). The study also highlighted species with relatively small suitability ranges that had less than 5% of this under protection, including Cephalanthera rubra (L.) Rich. (mean range = 271 km2; mean range protected = 0.9%) and Ophrys scolopax subsp. cornuta (Steven) E. G. Camus, P. Bergon & A. A. Camus (mean range = 280 km2; mean range protected = 2.48%) (Figure 2), recognized in Armenia as Ophrys oestrifera M. Bieb. and a Near Threatened species according to the National Red List (Table S1 Orchid species occurrence data).

4. Discussion

Compared to its neighbors, Armenia holds a relatively high number of orchid species for its territorial size. Iran, for example, holds approximately 46 species and subspecies [62], but is more than 14 times larger in land mass than Armenia. The diversity per area of land is comparable to neighboring Georgia, with ~57 species of orchids within a territory that is approximately twice the size [63]. This relatively high diversity of orchids within Armenia, therefore, warrants further protection, both through in situ and ex situ means.
Overall, current PA networks performed well in conserving the orchid species we found, with approximately half of the PA sites capturing orchid species. The distribution pattern of high orchid richness largely focuses on the central and south of the country (Tavush and Syunik provinces), which also correspond to the typical focus for the enhancement of PA networks within the country. After Lake Sevan, protected based on its importance as a key freshwater ecosystem resource, the two largest PAs are within the Syunik and Tavush provinces, respectively [37]. The relatively large size, coupled with PA locality within these provinces, can potentially explain the high orchid richness and diversity captured within. However, over a third of orchids were not represented within any PA site, with a further eight with only one occurrence record within PA. Of particular concern are those recognized as locally threatened, such as Steveniella satyrioides. Further uploads onto the WDPA database of newly endorsed PA sites for Armenia are pending, [41]; therefore, it is possible that these new sites will capture the orchid taxa left behind. A clear example is the Khustsup Nature Reserve, which is established, but not yet uploaded into the WDPA system. In addition to enhancing the PA network to include these species, additional conservation measures should be explored, such as ex situ strategies [64], although we recognize that orchids are known as “exceptional species” and may require the use of unconventional seed banking techniques (i.e., cryopreservation and/or in vitro) to ensure long-term conservation of orchid species diversity. Similarly, species only found within one PA and associated with known threatened habitats should also be prioritized for additional studies to fully determine their threatened status within the national context. For example, Corallorhiza trifida, recorded in the 1970s and found within forests, could now be lost due to much of the forest being cleared during the economic crises in the 1990s. Therefore, we would prioritize the need to revisit the locality of older records for species with low representation within PAs.
One major issue that we faced was around the taxonomic resolution of orchid species, as several species were recognized as synonyms, thereby reducing the overall species richness and diversity and/or creating uncertainty around conservation status. In our case, two morphologically distinct species, Orchis simia Vill., classified as Near Threatened, and O. tridentata Scop., classified as Endangered nationally, are both known as synonyms of Neotinea tridentata. These taxonomic resolution issues have potential implications within species modeling as the data of two species recognized as synonyms are combined under one globally accepted name, thus potentially overinflating the suitability range for those known locally as narrow endemics, and for N. tridentata, the already narrow suitability range found within our study would be of greater concern if this was the case. Difficulties in species-level taxonomy come at a time where taxonomic certainly is increasingly sought after by policy-makers and conservationists. For example, the Global Strategy for Plant Conservation has a target to assess 25% of all plant species, which is largely not reached due to incomplete taxonomic information [65,66]. Complex families, such as Orchidaceae, particularly in plant-rich regions, would need to be prioritized for resolving taxonomically ambiguous taxa. One way that this can be achieved is through integrated taxonomy, which considers both traditional morphological studies and DNA sequence data. Since its inception almost two decades ago, the use of DNA taxonomy has increased substantially with regard to reconstructing phylogenies, particularly for complex families [67,68]. In the Caucasus, there is increasing focus towards the use of DNA within plant-related studies [68]. The region-wide barcoding project currently underway in the Caucasus (CaBOL project: https://ggbc.eu/ (accessed on 28 September 2021)) is considered a first step in capturing DNA data for use in this way, but continued investment in integrated taxonomic studies has the potential to enhance conservation outcomes within biodiverse regions [69,70,71].
We need to acknowledge the risk of sampling bias with our study, as the occurrence data used were derived predominantly from herbarium records and a publicly accessed georeferenced database (GBIF). As these kinds of data typically do not utilize systematic methodologies, the resulting dataset can be subject to various aspects of sampling bias, such as temporal [72] and geographic bias (e.g., accessibility or known biodiverse areas such as outskirts of PAs) [73,74]. Using species richness estimates does control for some of these biases; however, we still have the risk of false negatives across our study site (i.e., zero occurrences due to lack of visitation rather than true absence), which is an issue across the study of conservation [73]. Indeed, research continues to discover new populations/localities of orchid species across Armenia [75,76]. Nonetheless, our study did find a potential orchid richness hotspot within the Lori province using our observational data. The Lori province lakes/wetlands have also been highlighted as a key area for plant diversity overall [77]; however, the current PA network within the province is much reduced in comparison with Tavush and Syunik. Therefore, we recommend further investigation within this province to estimate the potential of extending current conservation measures. A more comprehensive recommendation is to enhance PA networks using Key Biodiversity Areas (KBAs), which was described for the Caucasus region in 2020 [37]. There are currently 22 KBA sites within Armenia, with only 3% covered under strict protection (i.e., IUCN Category I) [37]. The areas that we proposed within our study coincide with the current KBAs planned. There is also an opportunity to enhance species protection through incorporating orchid-specific management plans for PAs that are already holding populations of orchids within, as highlighted by our study.
The use of SDMs gives an indication of a species’ potential ecological niche and can be used to predict current distribution [78,79,80]. We found that, overall, the models that were used fitted well, with over 60% showing an AUC value of above 0.7. However, estimating model suitability using AUC has been shown to be problematic, particularly with evidence of sampling bias, and, therefore, model interpretations should be done with caution [81]. Models of three orchid species did show low predictive power, namely Epipactis persica, Platanthera chlorantha and Gymnadenia conopsea. For E. persica, model improvement could be achieved through increasing the input of presence data; however, for the other two species, further analysis to incorporate additional environmental variables could be warranted. Nonetheless, our study makes an initial attempt at identifying potential suitability areas for orchids in Armenia, and, in general, our models reflect the known ecology of the target species. For example, Traunsteinera sphaerica was shown to be related to areas where precipitation rates remained high, with low maximum temperatures during the warmest months, typical of wetlands/marshes [45]. This close association could also explain their relatively restricted range of suitable areas across Armenia, with highly predicted sites found towards the northern part of the country, where it is relatively humid, with no distinct dry season. Therefore, conservation measures could include increasing the protection and effective management of existing wetland areas and/or restoring degraded marshland areas within the north of the country. The limitation for our model is the exclusion of variables related to habitat, which is commonly used in other modeling studies [82,83]. However, due to the wide temporal range of occurrence data used, some from the 1990s, and the rapidly changing landscape in Armenia due to economic crises, illegal logging and the rise of mining activities, the use of habitat variables could reduce the overall accuracy of predictions. We recognize that further work is needed to refine our models, including increasing field observations for model validation, the incorporation of more refined ecological variables and gathering recent landscape-level changes. Additionally, only a subset of the orchids found in Armenia had enough data points to be suitable for modeling. With increasing threats relating to habitat loss within the region [37], there is an urgent need to accelerate the gathering of locality, ecological and population-related data for these species.
We did not find a significant relationship between species range and proportion under protection, as species with the potential of inhabiting a wider climatic niche had varied protection levels, and vice versa. It is currently uncertain whether additional data for species missing from current SDM analysis will shed further light on a significant relationship between range and protection, which is seen elsewhere [12]. Additionally, the great advantage of SDMs is the ability to predict future distributions and threats regarding changes in species’ potential range, which can assist planning for a sustainable conservation strategy in the face of a changing climate [13,84].

5. Conclusions

In conclusion, the establishment and management of PA networks can be a costly process for many nations. Targeting areas where there is high diversity can be an effective strategy to ensure the maximum amount of protection. For Armenia, the current in situ conservation strategies seem to be effective at providing conservation impact for orchids; however, improvements can be made through increasing the proportion of strictly protected sites within already established KBA networks, coupled with enhancing the capacity for ex situ conservation for orchids, particularly those that are found outside of protected sites.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/d13120624/s1, Figure S1: Map of genera distribution; Figure S2: PCA plot; Figure S3: Species-specific distribution models from MaxEnt; Table S1: Orchid species occurrence data; Table S2: Species presence/absence within PA network in Armenia; Table S3: Principal Component Analysis results; Table S4: AICc values for model selection using regularization.

Author Contributions

Conceptualization, A.F. and A.N.; methodology, A.F.; resources, A.P. and A.N.; data curation, A.F., A.N. and A.P.; writing—original draft preparation, A.F.; writing—review and editing, A.N. and A.P.; visualization, A.F.; supervision, A.N.; project administration, A.F.; funding acquisition, A.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the RBG Kew Pilot Fund.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data will be available upon request to the corresponding author due to the sensitivity of threatened species’ locality.

Acknowledgments

We would like to express our sincere gratitude to the staff of the herbarium of the Komarov Botanical Institute of the Russian Academy of Sciences, St. Petersburg and personally to L. Averyanov, the Head of LE, and V. Shvanova, a curator of the Caucasian section of LE. We are thankful to our colleague, zoologist Vasil Ananian, for providing us with important photo evidence of orchid species, and to Jenny Williams from RBG Kew for her advice on species distribution and mapping. We also thank volunteer Esther Man for assisting in the curation of the tables for this work. Finally, we thank the five reviewers who took their time to review, comment and make suggestions regarding the manuscript, which allowed us to greatly improve it.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. WWF. Living Planet Report 2020—Bending the Curve of Biodiversity Loss; Almond, R.E.A., Grooten, M., Petersen, T., Eds.; WWF: Gland, Switzerland, 2020; ISBN 978-2-940529-99-5. [Google Scholar]
  2. IPBES. Global Assessment Report on Biodiversity and Ecosystem Services of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services; IPBES Secretariat: Bonn, Germany, 2019. [Google Scholar]
  3. Khapugin, A.A.; Kuzmin, I.V.; Silaeva, T.B. Anthropogenic drivers leading to regional extinction of threatened plants: Insights from regional Red Data Books of Russia. Biodivers. Conserv. 2020, 29, 2765–2777. [Google Scholar] [CrossRef]
  4. Le Roux, J.J.; Hui, C.; Castillo, M.L.; Iriondo, J.M.; Keet, J.H.; Khapugin, A.A.; Médail, F.; Rejmánek, M.; Theron, G.; Yannelli, F.A.; et al. Recent Anthropogenic Plant Extinctions Differ in Biodiversity Hotspots and Coldspots. Curr. Biol. 2019, 29, 2912–2918.e2. [Google Scholar] [CrossRef] [PubMed]
  5. Štípková, Z.; Kindlmann, P. Orchid Extinction over the Last 150 Years in the Czech Republic. Diversity 2021, 13, 78. [Google Scholar] [CrossRef]
  6. Kougioumoutzis, K.; Kokkoris, I.P.; Panitsa, M.; Strid, A.; Dimopoulos, P. Extinction Risk Assessment of the Greek Endemic Flora. Biology 2021, 10, 195. [Google Scholar] [CrossRef]
  7. Knapp, W.M.; Frances, A.; Noss, R.; Naczi, R.F.C.; Weakley, A.; Gann, G.D.; Baldwin, B.G.; Miller, J.; McIntyre, P.; Mishler, B.D.; et al. Vascular plant extinction in the continental United States and Canada. Conserv. Biol. 2021, 35, 360–368. [Google Scholar] [CrossRef]
  8. Nic Lughadha, E.; Bachman, S.P.; Leão, T.C.C.; Forest, F.; Halley, J.M.; Moat, J.; Acedo, C.; Bacon, K.L.; Brewer, R.F.A.; Gâteblé, G.; et al. Extinction risk and threats to plants and fungi. Plants People Planet 2020, 2, 389–408. [Google Scholar] [CrossRef]
  9. Gray, C.L.; Hill, S.L.; Newbold, T.; Hudson, L.N.; Börger, L.; Contu, S.; Hoskins, A.J.; Ferrier, S.; Purvis, A.; Scharlemann, J.P. Local biodiversity is higher inside than outside terrestrial protected areas worldwide. Nat. Commun. 2016, 7, 12306. [Google Scholar] [CrossRef] [Green Version]
  10. Aubele, D.T.; Wang, D.; Post, K.W.; Hahn, S.; Wincott, F. Effectiveness of Parks in Protecting Tropical Biodiversity. Nucleic Acids Res. 1997, 17, 5191. [Google Scholar]
  11. CBD. Zero Draft of Post—2020 Biodiversity Framework. 2020. Available online: https://www.cbd.int/conferences/post2020/wg2020-02/documents (accessed on 29 September 2021).
  12. Rodrigues, A.S.L.; Andelman, S.J.; Bakarr, M.I.; Boitani, L.; Brooks, T.M.; Cowling, R.M.; Fishpool, L.D.C.; da Fonseca, G.A.B.; Gaston, K.J.; Pilgrim, J.D.; et al. Effectiveness of the global protected area network in representing species diversity. Nature 2004, 428, 640–643. [Google Scholar] [CrossRef]
  13. Dobrowski, S.Z.; Littlefield, C.E.; Lyons, D.S.; Hollenberg, C.; Carroll, C.; Parks, S.A.; Abatzoglou, J.T.; Hegewisch, K.; Gage, J. Protected-area targets could be undermined by climate change-driven shifts in ecoregions and biomes. Commun. Earth Environ. 2021, 2, 198. [Google Scholar] [CrossRef]
  14. Beyer, H.L.; Venter, O.; Grantham, H.S.; Watson, J.E. Substantial losses in ecoregion intactness highlight urgency of globally coordinated action. Conserv. Lett. 2020, 13, e12692. [Google Scholar] [CrossRef]
  15. Mancheno, C.S.M.; Zazanashvili, N.; Beruchashvili, G. Effectiveness of the network of protected areas of the South Caucasus at representing terrestrial ecosystems after the dissolution of the Soviet Union. Environ. Conserv. 2017, 44, 158–165. [Google Scholar] [CrossRef]
  16. Delso, Á.; Fajardo, J.; Muñoz, J. Protected area networks do not represent unseen biodiversity. Sci. Rep. 2021, 11, 12275. [Google Scholar] [CrossRef]
  17. Milchev, B.; Georgiev, V. The effect of N atura 2000 network on the Eurasian Eagle-owl (Bubo bubo) population in Southeast Bulgaria: Implications for conservation. Ornis Hung. 2021, 29, 170–178. [Google Scholar] [CrossRef]
  18. Koskela, J.; Lefèvre, F.; Schueler, S.; Kraigher, H.; Olrik, D.C.; Hubert, J.; Longauer, R.; Bozzano, M.; Yrjänä, L.; Alizoti, P.; et al. Translating conservation genetics into management: Pan-European minimum requirements for dynamic conservation units of forest tree genetic diversity. Biol. Conserv. 2013, 157, 39–49. [Google Scholar] [CrossRef] [Green Version]
  19. Oberosler, V.; Tenan, S.; Zipkin, E.F.; Rovero, F. Poor management in protected areas is associated with lowered tropical mammal diversity. Anim. Conserv. 2020, 23, 171–181. [Google Scholar] [CrossRef]
  20. Moreno-Saiz, J.C.; Albertos, B.; Ruiz-Molero, E.; Mateo, R.G. The European Union can afford greater ambition in the conservation of its threatened plants. Biol. Conserv. 2021, 261, 109231. [Google Scholar] [CrossRef]
  21. Guisan, A.; Tingley, R.; Baumgartner, J.B.; Naujokaitis-Lewis, I.; Sutcliffe, P.R.; Tulloch, A.I.T.; Regan, T.J.; Brotons, L.; Mcdonald-Madden, E.; Mantyka-Pringle, C.; et al. Predicting species distributions for conservation decisions. Ecol. Lett. 2013, 16, 1424–1435. [Google Scholar] [CrossRef] [PubMed]
  22. Balding, M.; Williams, K.J.H. Plant blindness and the implications for plant conservation. Conserv. Biol. 2016, 30, 1192–1199. [Google Scholar] [CrossRef] [PubMed]
  23. Govearts, R.; Bernet, P.; Kratochvil, K.; Gerlach, G.; Carr, G.; Alrich, P.; Pridgeon, A.M.; Pfahl, J.; Campacci, M.A.; Holland Baptista, D.; et al. World Checklist of Orchidaceae; The Board of Trustees of the Royal Botanic Gardens: Richmond, UK, 2008. [Google Scholar]
  24. Esposito, F.; Merckx, T.; Tyteca, D.; Esposito, F.; Merckx, T.; Tyteca, D. Noctuid moths as potential hybridization agents for Platanthera orchids. Lankesteriana 2017, 17, 383–393. [Google Scholar] [CrossRef] [Green Version]
  25. Ayasse, M.; Schiestl, F.P.; Paulus, H.F.; Ibarra, F.; Francke, W. Pollinator attraction in a sexually deceptive orchid by means of unconventional chemicals. Proc. R. Soc. Lond. B 2003, 270, 517–522. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Paulus, H.F. Speciation, pattern recognition and the maximization of pollination: General questions and answers given by the reproductive biology of the orchid genus Ophrys. J. Comp. Physiol. A 2019, 205, 285–300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Turco, A.; Medagli, P.; Wagensommer, R.P.; D’Emerico, S.; Gennaio, R.; Albano, A. A morphometric study on Ophrys sect. Pseudophrys in Apulia (Italy) and discovery of Ophrys japigiae sp. nov. (Orchidaceae). Plant Biosyst. Int. J. Deal. All Asp. Plant Biol. 2021. [Google Scholar] [CrossRef]
  28. Schiestl, F.P. On the success of a swindle: Pollination by deception in orchids. Naturwissenschaften 2005, 92, 255–264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Calevo, J.; Voyron, S.; Adamo, M.; Alibrandi, P.; Perotto, S.; Girlanda, M. Can orchid mycorrhizal fungi be persistently harbored by the plant host? Fungal Ecol. 2021, 53, 101071. [Google Scholar] [CrossRef]
  30. Esposito, F.; Jacquemyn, H.; Waud, M.; Tyteca, D. Mycorrhizal fungal diversity and community composition in two closely related Platanthera (orchidaceae) species. PLoS ONE 2016, 11, e0164108. [Google Scholar] [CrossRef] [PubMed]
  31. Swarts, N.D.; Dixon, K.W. Terrestrial orchid conservation in the age of extinction. Ann. Bot. 2009, 104, 543–556. [Google Scholar] [CrossRef] [Green Version]
  32. Averyanov, L.V. Averyanov, L.V. A Review of the Genus Dactylorhiza. In Orchid Biology: Reviews and Perspectives; Arditti, J., Ed.; Timber Press, Inc.: Portland, OR, USA, 1990; Volume V, pp. 161–204. ISBN 0-88192-170-X. [Google Scholar]
  33. Khapugin, A.A. A global systematic review on orchid data in protected areas. Nat. Conserv. Res. 2020, 5, 19–33. [Google Scholar] [CrossRef]
  34. POWO Plants of the World Online. Facilitated by the Royal Botanic Gardens, Kew. Available online: https://powo.science.kew.org/ (accessed on 1 September 2021).
  35. Rankou, H. Goodyera Macrophylla. Available online: https://www.iucnredlist.org/species/162070/5527443 (accessed on 29 September 2021). [CrossRef]
  36. Wagensommer, R.P.; Medagli, P.; Turco, A.; Perrino, E.V. IUCN red list evaluation of the orchidaceae endemic to apulia (Italy) and considerations on the application of the IUCN protocol to rare species. Nat. Conserv. Res. 2020, 5, 90–101. [Google Scholar] [CrossRef]
  37. Zazanashvili, N.; Sanadiradze, G.; Garforth, M.; Bitsadze, M.; Manvelyan, K.; Askerov, E.; Mousavi, M.; Krever, V.; Shmuk, V.; Kalem, S.; et al. Ecoregional Conservation Plan (ECP) for The Caucasus 2020 Edition; WWF, KfW: Tbilisi, Georgia, 2020. [Google Scholar]
  38. Schuerholz, G. Lake Arpi National Park and Sanctuaries: Akhuryan Gorge, Ardenis and Alvar; WWF, KfW: Gyumri, Armenia, 2009. [Google Scholar]
  39. Zazanashvili, N.; Mallon, D. Status and Protection of Globally Threatened Species in the Caucasus; Zazanashvili, N., Mallon, D., Eds.; Contour Ltd.: Tbilisi, Georgia, 2009. [Google Scholar]
  40. Mousavi, M.; Moqanaki, E.; Hamidi, A.K.; Breitenmoser, U.; Breitenmoser-Würsten, C.; Zazanashvili, N.; Heidelberg, A. Conservation of the Leopard in the Caucasus: Iran Chapter International Experts Workshop “Conservation of the Leopard in the Caucasus” Workshop Report; Tbilisi, 2014. Available online: https://www.researchgate.net/publication/274256377_Conservation_of_the_leopard_in_the_Caucasus_Iran_chapter (accessed on 29 September 2021).
  41. IUCN; UNEP-WCMC. The World Database on Protected Areas (WDPA). Available online: https://www.protectedplanet.net/en/thematic-areas/wdpa?tab=WDPA (accessed on 15 October 2021).
  42. Thiers, B. Index Herbariorium. Available online: http://sweetgum.nybg.org/science/ih/herbarium-list/?NamOrganisationAcronym=MW (accessed on 15 October 2021).
  43. Noé, N. GBIF Occurrences Plugin for QGIS 3. 2019. Available online: https://plugins.qgis.org/plugins/qgisgbifapi/ (accessed on 10 September 2021).
  44. GBIF.org Occurrence Download. Available online: https://www.gbif.org/occurrence/download/0014012-210819072339941 (accessed on 10 September 2021).
  45. Averianov, L.V.; Nersesian, A.A. Orchidaceae Juss. In Flora of Armenia V. 10; Takhtajan, A.L., Ed.; A. R. G. Gantner Verlag KG: Ruggell, Leichtenstein, 2001; pp. 165–225. [Google Scholar]
  46. Averianov, L.V. Orchidaceae Juss. In Conspectus Florae Caucasi; Takhtajan, A.L., Ed.; Saint-Petersburg University Press: Saint-Petersburg, Russia, 2006; pp. 84–101. [Google Scholar]
  47. WCVP World Checklist of Vascular Plants, Version 2.0. Facilitated by the Royal Botanic Gardens, Kew. Available online: http://wcvp.science.kew.org/ (accessed on 9 September 2021).
  48. Hijmans, R.J.; Guarino, L.; Jarvis, A.; O’Brian, R.; Mathur, P.; Busink, C.; Cruz, M.; Barrantes, I.; Rojas, E. DIVA-GIS: Freely Available Resource. Available online: www.diva-gis.org (accessed on 15 September 2021).
  49. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Roeland, K.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, B.R.; Simpson, G.L.; Solymos, P.; et al. Vegan: Community Ecology Package. R Package Version 2.5-7. 2020. Available online: https://CRAN.R-project.org/package=vegan (accessed on 20 September 2021).
  50. R Version 4.0.2 Copyright © The R Foundation for Statistical Computing Platform 2020. Available online: https://www.r-project.org/about.html (accessed on 10 October 2021).
  51. Phillips, S.J.; Dudik, M.; Schapire, R. Maxent Software for Modeling Species Niches and Distributions (Version 3.4.1). 2021. Available online: https://biodiversityinformatics.amnh.org/open_source/maxent/ (accessed on 15 September 2021).
  52. Deb, C.R.; Jamir, N.S.; Kikon, Z.P. Distribution Prediction Model of a Rare Orchid Species (Vanda bicolor Griff.) Using Small Sample Size. Am. J. Plant Sci. 2017, 8, 1388–1398. [Google Scholar] [CrossRef]
  53. Ma, X.; Chen, T.; Zhang, G.; Wang, R. Microbial community structure along an altitude gradient in three different localities. Folia Microbiol. 2004, 49, 105–111. [Google Scholar] [CrossRef] [PubMed]
  54. Siles, J.A.; Margesin, R. Seasonal soil microbial responses are limited to changes in functionality at two Alpine forest sites differing in altitude and vegetation. Sci. Rep. 2017, 7, 2204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Baumann, K.; Keune, J.; Wolters, V.; Jauker, F. Distribution and pollination services of wild bees and hoverflies along an altitudinal gradient in mountain hay meadows. Ecol. Evol. 2021, 11, 11345–11351. [Google Scholar] [CrossRef] [PubMed]
  56. Viana, T.A.; Martins, F.M.; Lourenço, A.P. The Orchid Bee Fauna (Hymenoptera: Apidae: Euglossini) of a Neotropical Savanna: An Efficient Protocol to Assess Bee Community and Diversity Along Elevational and Habitat Complexity Gradients. Neotrop. Entomol. 2021, 50, 748–758. [Google Scholar] [CrossRef]
  57. Hijmans, R.J.; Phillips, S.J.; Leathwick, J.R.; Elith, J. Dismo: Species Distribution Modeling, R Package Version 1.3-3; 2020. Available online: https://CRAN.R-project.org/package=dismo (accessed on 15 September 2021).
  58. Phillips, S.J.; Anderson, R.P.; Dudík, M.; Schapire, R.E.; Blair, M.E. Opening the black box: An open-source release of Maxent. Ecography 2017, 40, 887–893. [Google Scholar] [CrossRef]
  59. Akaike, H. A New Look at the Statistical Model Identification. IEEE Trans. Automat. Contr. 1974, 19, 716–723. [Google Scholar] [CrossRef]
  60. Warren, D.L.; Seifert, S.N. Ecological niche modeling in Maxent: The importance of model complexity and the performance of model selection criteria. Ecol. Appl. 2011, 21, 335–342. [Google Scholar] [CrossRef] [Green Version]
  61. Tamanyan, K.; Fayvuh, G.; Nanagulyan, S. Higher Plants and Fungi. In The Red Data Book of Plants of the Republic of Armenia; Danielyan, T., Ed.; Printing-House of “Zangak-97”, LLC: Yerevan, Armenia, 2010; p. 592. [Google Scholar]
  62. Shahsavari, A. Flora of Iran. Part 57: Orchidaceae; Research Institute of Forests and Rangelands: Tehran, Iran, 2008. [Google Scholar]
  63. Akhalkatsi, M.; Arabuli, G.; Lorenz, R. Orchids as indicator species of forest disturbances on limestone quarry in Georgia (South Caucasus). J. Eur. Orchid. 2014, 46, 123–160. [Google Scholar]
  64. IUCN/SSC Orchid Specialist Group. Status Survey and Conservation Action Plan Orchids; IUCN: Gland, Switzerland; Cambridge, UK, 1996. [Google Scholar]
  65. Feeley, K.J.; Silman, M.R. Keep collecting: Accurate species distribution modelling requires more collections than previously thought. Divers. Distrib. 2011, 17, 1132–1140. [Google Scholar] [CrossRef]
  66. Meyer, C.; Weigelt, P.; Kreft, H. Multidimensional biases, gaps and uncertainties in global plant occurrence information. Ecol. Lett. 2016, 19, 992–1006. [Google Scholar] [CrossRef]
  67. Muñoz-Rodríguez, P.; Carruthers, T.; Wood, J.R.I.; Williams, B.R.M.; Weitemier, K.; Kronmiller, B.; Goodwin, Z.; Sumadijaya, A.; Anglin, N.L.; Filer, D.; et al. A taxonomic monograph of Ipomoea integrated across phylogenetic scales. Nat. Plants 2019, 5, 1136–1144. [Google Scholar] [CrossRef] [PubMed]
  68. Korotkova, N.; Parolly, G.; Khachatryan, A.; Ghulikyan, L.; Sargsyan, H.; Akopian, J.; Borsch, T.; Gruenstaeudl, M. Towards resolving the evolutionary history of Caucasian pears (Pyrus, Rosaceae)—Phylogenetic relationships, divergence times and leaf trait evolution. J. Syst. Evol. 2018, 56, 35–47. [Google Scholar] [CrossRef]
  69. Cristescu, M.E. From barcoding single individuals to metabarcoding biological communities: Towards an integrative approach to the study of global biodiversity. Trends Ecol. Evol. 2014, 29, 566–571. [Google Scholar] [CrossRef] [PubMed]
  70. Vogel Ely, C.; Bordignon, S.A.D.L.; Trevisan, R.; Boldrini, I.I. Implications of poor taxonomy in conservation. J. Nat. Conserv. 2017, 36, 10–13. [Google Scholar] [CrossRef]
  71. Sheth, B.P.; Thaker, V.S. DNA barcoding and traditional taxonomy: An integrated approach for biodiversity conservation. Genome 2017, 60, 618–628. [Google Scholar] [CrossRef] [Green Version]
  72. Haque, M.D.M.; Nipperess, D.A.; Baumgartner, J.B.; Beaumont, L.J. A journey through time: Exploring temporal patterns amongst digitized plant specimens from Australia. Syst. Biodivers. 2018, 16, 604–613. [Google Scholar] [CrossRef]
  73. Hughes, A.C.; Orr, M.C.; Ma, K.; Costello, M.J.; Waller, J.; Provoost, P.; Yang, Q.; Zhu, C.; Qiao, H. Sampling biases shape our view of the natural world. Ecography 2021, 44, 1259–1269. [Google Scholar] [CrossRef]
  74. Fisher-Phelps, M.; Cao, G.; Wilson, R.M.; Kingston, T. Protecting bias: Across time and ecology, open-source bat locality data are heavily biased by distance to protected area. Ecol. Inform. 2017, 40, 22–34. [Google Scholar] [CrossRef]
  75. Asatryan, A. New data on distribution of some rare plant species (Pyrus gergerana Gladkova, P. daralagezi Mulk., P. voronovii Rubtzov, Orchis punctulata Steven ex Lindl.) in Armenia. Takhtajania 2018, 4, 51–56. [Google Scholar]
  76. Nersesyan, A. New data on families Apiaceae, Caryophyllaceae, Orchidaceae, Poaceae of the Armenian flora. Takhtajania 2018, 4, 53–56. [Google Scholar]
  77. Fayvush, G.; Tamanyan, K.; Kalahsyan, M.; Vitek, E. “Biodiversity Hotspots” in Armenia. Ann. Naturhist. Mus. Wien B 2013, 115, 11–20. [Google Scholar]
  78. Fois, M.; Cuena-Lombraña, A.; Fenu, G.; Cogoni, D.; Bacchetta, G. Does a correlation exist between environmental suitability models and plant population parameters? An experimental approach to measure the influence of disturbances and environmental changes. Ecol. Indic. 2018, 86, 1–8. [Google Scholar] [CrossRef]
  79. Hosseinzadeh, M.S.; Fois, M.; Zangi, B.; Kazemi, S.M. Predicting past, current and future habitat suitability and geographic distribution of the Iranian endemic species Microgecko latifi (Sauria: Gekkonidae). J. Arid Environ. 2020, 183. [Google Scholar] [CrossRef]
  80. Abdelaal, M.; Fois, M.; Fenu, G.; Bacchetta, G. Using MaxEnt modeling to predict the potential distribution of the endemic plant Rosa arabica Crép. in Egypt. Ecol. Inform. 2019, 50, 68–75. [Google Scholar] [CrossRef]
  81. Fourcade, Y.; Engler, J.O.; Rödder, D.; Secondi, J. Mapping species distributions with MAXENT using a geographically biased sample of presence data: A performance assessment of methods for correcting sampling bias. PLoS ONE 2014, 9, e97122. [Google Scholar] [CrossRef] [Green Version]
  82. Wan, J.; Wang, C.; Han, S.; Yu, J. Planning the priority protected areas of endangered orchid species in northeastern China. Biodivers. Conserv. 2014, 23, 1395–1409. [Google Scholar] [CrossRef]
  83. Štípková, Z.; Kindlmann, P. Factors determining the distribution of orchids—A review with examples from the Czech Republic. Eur. J. Environ. Sci. 2021, 11, 21–30. [Google Scholar] [CrossRef]
  84. Schwager, P.; Berg, C. Global warming threatens conservation status of alpine EU habitat types in the European Eastern Alps. Reg. Environ. Chang. 2019, 19, 2411–2421. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Map showing (a) a heatmap of orchid observations (genera distribution can be found in Figure S1: Map of genera distribution) and (b) species richness within the territory of Armenia, derived from herbarium vouchers, field visits and GBIF records (1101 total data points). Heat map developed using the Kernel Density Estimation with a radial distance buffer of 5 km is graded from yellow to red showing increasing observational findings. Location of Protected Areas derived from the World Database of Protected Areas is also mapped across the study site for comparison.
Figure 1. Map showing (a) a heatmap of orchid observations (genera distribution can be found in Figure S1: Map of genera distribution) and (b) species richness within the territory of Armenia, derived from herbarium vouchers, field visits and GBIF records (1101 total data points). Heat map developed using the Kernel Density Estimation with a radial distance buffer of 5 km is graded from yellow to red showing increasing observational findings. Location of Protected Areas derived from the World Database of Protected Areas is also mapped across the study site for comparison.
Diversity 13 00624 g001aDiversity 13 00624 g001b
Figure 2. Scatterplot showing species potential climatic suitability range (km2) against the mean proportion of species range within a PA represented as a percentage. Colors represent genera according to the World Checklist of Vascular Plants.
Figure 2. Scatterplot showing species potential climatic suitability range (km2) against the mean proportion of species range within a PA represented as a percentage. Colors represent genera according to the World Checklist of Vascular Plants.
Diversity 13 00624 g002
Table 1. Orchid species richness, diversity (Shannon–Weiner Index), area and designation of Protected Area (PA) site found within the territory of Armenia listed by province from north to south. Number of species per PA was calculated using the georeferenced data from herbarium vouchers, field visits and GBIF records. Shannon–Weiner Index was calculated using the “vegan” package [49] in R statistical software [50]. PA names, their reported area and designation are from the World Database of Protected Areas, downloaded from the Protected Planet website [41].
Table 1. Orchid species richness, diversity (Shannon–Weiner Index), area and designation of Protected Area (PA) site found within the territory of Armenia listed by province from north to south. Number of species per PA was calculated using the georeferenced data from herbarium vouchers, field visits and GBIF records. Shannon–Weiner Index was calculated using the “vegan” package [49] in R statistical software [50]. PA names, their reported area and designation are from the World Database of Protected Areas, downloaded from the Protected Planet website [41].
Protected AreaOrchid Species RichnessOrchid Diversity IndexReported Area (km2)DesignationIUCN Category
Shirak
Lake Arpi30.85212National ParkII
Lori
Gyulagarak42.4125.8SanctuaryIV
Zikatar10.001.50SanctuaryIV
Margahovit41.3633.7SanctuaryIV
Tavush
Ijevan142.4159.1SanctuaryIV
Gandzakar10.0068.1SanctuaryIV
Dilijan172.50377National ParkII
Akhnabat Yew Grove10.000.25SanctuaryIV
Gegharkunik
Sevan61.551474National ParkII
Juniper Open Woodland20.6933.1SanctuaryIV
Kotayk
Arzakan-Meghradzor41.21135SanctuaryIV
Ararat
Khosrov Forest61.50232State ReserveIa
Vayots Dzor
Jermuk Forest10.0038.7SanctuaryIV
Herher Open Woodland10.0061.4SanctuaryIV
Syunik
Arevik172.30344National ParkII
Boghaqar10.0027.3SanctuaryIV
Shikahogh132.26121State ReserveIa
Zangezur20.69259SanctuaryIV
Goris10.0018.5SanctuaryIV
Plane Grove10.000.64SanctuaryIV
Table 2. Results from MaxEnt species distribution models for orchids in Armenia. The presence records are the number of points used to build each model, while 25% of the data were set aside to be used to test the model (Test records). The features (l = linear, q = quadratic, h = hinge) and the regularization multiplier used for the best fit model per species are shown. Best-fit model was chosen based on comparing the AICc values of models with regularization multiplier from 1 to 10 for each species (see Table S4: AICc values for model selection using regularization for results). Training and test AUC, with corresponding standard deviation (SD), are also presented, along with the estimate of contributions of the environmental variables (normalized to percentage). Environmental variables were chosen based on a Principal Component Analysis, with max temperature of warmest month (Bio5) corresponding to Principal Component 1 (PC1), precipitation of warmest quarter (Bio18) corresponding to PC2, precipitation of driest quarter (Bio17) corresponding to PC3 and altitude (Alt) corresponding to PC4.
Table 2. Results from MaxEnt species distribution models for orchids in Armenia. The presence records are the number of points used to build each model, while 25% of the data were set aside to be used to test the model (Test records). The features (l = linear, q = quadratic, h = hinge) and the regularization multiplier used for the best fit model per species are shown. Best-fit model was chosen based on comparing the AICc values of models with regularization multiplier from 1 to 10 for each species (see Table S4: AICc values for model selection using regularization for results). Training and test AUC, with corresponding standard deviation (SD), are also presented, along with the estimate of contributions of the environmental variables (normalized to percentage). Environmental variables were chosen based on a Principal Component Analysis, with max temperature of warmest month (Bio5) corresponding to Principal Component 1 (PC1), precipitation of warmest quarter (Bio18) corresponding to PC2, precipitation of driest quarter (Bio17) corresponding to PC3 and altitude (Alt) corresponding to PC4.
SpeciesPresence RecordsTest
Records
FeaturesRegularization MultiplierAICc ValueTraining AUCTest AUCSDVariable Contribution (%)
Bio5Bio18Bio17Alt
Anacamptis coriophora246hlq14770.7200.6120.08498.50.11.30.1
Anacamptis pyramidalis287hlq54820.7160.8440.0540.00.00.0100
Cephalanthera damasonium164hlq22850.860 0.8660.0500.035.30.064.6
Cephalanthera longifolia113lq11970.8190.7830.0240.018.50.081.5
Cephalanthera rubra338hlq36010.8340.7790.0603.513.10.083.5
Dactylorhiza euxina154hlq22990.6860.6980.1020.00.40.099.6
Dactylorhiza incarnata
subsp. cilicica
308hlq35740.7950.7100.08420.550.20.029.4
Dactylorhiza urvilleana4511hlq58590.7380.6460.06456.213.630.10.0
Dactylorhiza viridis113lq12070.8350.7400.17385.38.85.90.0
Epipactis
helleborine
144lq12560.8100.7510.0880.026.10.073.9
Epipactis
persica
92l31620.7810.5720.0300.00.00.0100
Gymnadenia conopsea246hlq34590.6630.5970.0830.00.00.0100
Limodorum abortivum113lq11920.857 0.9080.0170.00.012.287.8
Neotinea
tridentata
154hlq22800.8940.8040.05979.02.518.40.0
Neottia nidus-avis164hlq22700.7820.7990.0830.083.51.614.9
Ophrys
scolopax subsp. cornuta
92l11690.8340.8390.01168.70.00.131.2
Orchis
mascula
4211hlq28130.7400.6070.08520.940.310.528.3
Orchis
punctulata
92l11700.8530.6600.08790.80.09.20.0
Orchis
purpurea
92l21600.7740.8400.0350.00.00.0100
Platanthera chlorantha236hlq34470.7020.5900.1220.083.80.016.2
Traunsteinera sphaerica72l11350.7700.7660.10849.850.20.00.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Faruk, A.; Papikyan, A.; Nersesyan, A. Exploring Effective Conservation of Charismatic Flora: Orchids in Armenia as a Case Study. Diversity 2021, 13, 624. https://0-doi-org.brum.beds.ac.uk/10.3390/d13120624

AMA Style

Faruk A, Papikyan A, Nersesyan A. Exploring Effective Conservation of Charismatic Flora: Orchids in Armenia as a Case Study. Diversity. 2021; 13(12):624. https://0-doi-org.brum.beds.ac.uk/10.3390/d13120624

Chicago/Turabian Style

Faruk, Aisyah, Astghik Papikyan, and Anush Nersesyan. 2021. "Exploring Effective Conservation of Charismatic Flora: Orchids in Armenia as a Case Study" Diversity 13, no. 12: 624. https://0-doi-org.brum.beds.ac.uk/10.3390/d13120624

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop