Vascular plant communities in the polar desert of Alert (Ellesmere Island, Canada): Establishment of a baseline reference for the 21st century

ABSTRACT

The startling warming of the Arctic is driving important environmental changes, but vegetation responses have been spatially heterogeneous and difficult to predict. In this context, establishing new ecological baselines and initiating monitoring schemes are essential. We conducted systematic plot-based surveys in the polar desert surrounding Alert (Nunavut, Canada). We aimed at (1) identifying distinct plant communities, (2) characterizing community attributes, including diversity and abundance, as well as environmental variables associated with each community, and (3) establishing a georeferenced baseline with permanent field markers allowing robust re-surveying. We used hierarchical clustering to categorize cover values of vascular plant species, cryptogams, and ground substrates from 1,320 quadrats (1 m² each) surveyed in 264 vegetation plots. Five plant communities were identified, with one community associated with each of the barren and mesic habitats, and three communities associated with wetlands. The mean biotic covers were generally higher at Alert (13–98%) compared to other polar deserts in the Canadian Arctic Archipelago. A total of 250 quadrats from 50 vegetation plots were permanently marked, and a database describing all plots is available online. This study improves our understanding of High-Arctic plant communities and establishes an important vegetation monitoring reference at the northernmost permanently inhabited settlement on Earth.

Résumé

Le réchauffement sans précédent de l’Arctique entraîne des changements environnementaux importants, mais les réponses de la végétation sont spatialement hétérogènes et difficiles à prédire. Dans ce contexte, l’établissement de nouvelles bases de référence écologiques et la mise en place de programmes de suivi sont essentiels. Nous avons effectué des relevés systématiques par quadrats dans le désert polaire entourant Alert (Nunavut, Canada). Nous avions pour objectifs (1) d’identifier les communautés végétales distinctes, (2) de caractériser les attributs de chaque communauté, y compris la diversité et l’abondance, ainsi que les variables environnementales associées à chacune et (3) d’établir une base de données géoréférencée avec des quadrats marqués de façon permanente permettant un rééchantillonnage rigoureux. Nous avons utilisé une analyse de groupement hiérarchique pour catégoriser les valeurs de recouvrement des espèces végétales vasculaires, des cryptogames et des types de substrats de 1320 quadrats (1 m² chacun) étudiés dans 264 parcelles de végétation. Cinq communautés végétales ont été identifiées, une associée aux habitats arides, une associée aux habitats mésiques, et trois associées aux habitats humides. Les recouvrements biotiques moyens étaient généralement plus élevés dans les communautés à Alert (13 à 98%) que dans celles d’autres déserts polaires de l’archipel Arctique canadien. Au total, 250 quadrats répartis dans 50 parcelles de végétation ont été marqués avec des clous, et une base de données détaillée décrivant toutes les parcelles est disponible en ligne. Cette étude contribue à une meilleure compréhension des communautés végétales du Haut-Arctique et établit une référence importante en matière de surveillance de la végétation à l’endroit habité en permanence le plus nordique de la planète.

KEYWORDS:

• Arctic
• vegetation
• monitoring
• point-intercept method
• plant cover
• rare species

Mots clés:

• Arctique
• végétation
• suivi
• méthode de point d’interception
• recouvrement végétal
• espèces rares

Introduction

Vegetation is essential to the functioning of Arctic terrestrial ecosystems and, through complex feedbacks and cascades, changes in plant communities greatly affect animal consumers and ecological processes such as water exchange, litter decomposition, nutrient cycling and soil fertility (Wookey et al. 2009; Callaghan et al. 2011; Bjorkman et al. 2020). Over the past two decades, average Arctic air surface temperatures have increased by more than twice the global average (IPCC 2019), with already clear responses from Arctic plant communities. Climate-induced shifts in vegetation diversity, biomass and abundance have been reported (Elmendorf et al. 2012; Epstein et al. 2012; Bjorkman et al. 2020), but are not universal across the Arctic (Tape et al. 2012; Pattison et al. 2015). Spatial heterogeneity likely stems from the diversity of drivers determining plant community composition, including macroclimate, topography, soil conditions, and biotic interactions (Nabe-Nielsen et al. 2017; Berner et al. 2020; Bjorkman et al. 2020), thus making predictions difficult. Efforts to address heterogeneity in tundra responses to warming have been hindered by lack of data. In addition, experiments were short term and conducted at few sites (Elmendorf et al. 2012). In this context, establishing baselines and implementing monitoring schemes to track the responses of Arctic plant communities to climate change are critical (Mihoub et al. 2017; Christensen et al. 2020).

A coordination of circumpolar monitoring effort has arisen from the Arctic Council working group Conservation of Arctic Flora and Fauna (CAFF), which established the Circumpolar Biodiversity Monitoring Programme to help implementing adaptive and standardized monitoring of Arctic environments (Petersen et al. 2004; CAFF 2018). This initiative included a monitoring plan specific to Arctic terrestrial biodiversity (Christensen et al. 2013), where vegetation is considered one of four key biotic groups. Within each biotic group, entities designated as Focal Ecosystem Components (FECs) were identified, with their associated core attributes, as targets for monitoring. Vascular and non-vascular plants (trees, shrubs, forbs, graminoids, mosses, and lichens) were all considered to be FECs (Christensen et al. 2013; Appendix A), whereas FEC essential attributes for vegetation included diversity, abundance, composition, phenology, temporal cycles, and ecosystem functions and processes. These attributes were chosen based on ecological significance, relevance to ecosystem services, value to local people, and importance for management and legislation needs (Christensen et al. 2013). The most recent report on status and trends of Arctic terrestrial biodiversity confirmed the validity of the above vegetation FECs and attributes (Aronsson et al. 2021).

More long-term monitoring programs have been implemented in the Low Arctic than in the High Arctic (Bjorkman et al. 2020), and additional vegetation monitoring in the High Arctic is thus greatly promoted (Christensen et al. 2020; Ravolainen et al. 2020). The High Arctic tundra biome corresponds to subzones A, B, and C of the Circumpolar Arctic Vegetation Map (CAVM Team 2003) and includes the Canadian Arctic Archipelago, the northern margins of Russia, the northern half of coastal Greenland, and Svalbard (Halliday 2002). The High Arctic tundra covers a variety of landforms and soil deposits, and is exposed to a range of climatic conditions promoting the development of three broad vegetation types: polar oases (relatively high plant biomass and density; Yurtsev 1994), polar semi-deserts (5–20% vascular plant cover; Bliss et al. 1973; Bliss and Svoboda 1984) and polar deserts (<5% vascular plant cover; Bliss et al. 1973; Bliss and Svoboda 1984). Polar deserts encompass 850,000 km² worldwide and 75% of them are found in the Canadian High Arctic (Bliss and Matveyeva 1992). Despite a relatively sparse vegetation cover, a low species diversity and the ubiquitous distribution of some species, polar deserts support various plant assemblages (Bay 1997; Lévesque 1997). They, however, remain less studied than polar oases and semi-deserts (Bliss and Svoboda 1984; Muc et al. 1989).

The vegetation of polar deserts is constrained by continuous permafrost, a short-growing season, low annual and summer mean temperatures, strong winds, low precipitations, and severe nutrient limitations (Billings 1987; Peterson 2014). Arctic tundra ecosystems are also characterized by a large variety of natural disturbances, either directly or indirectly climatically driven, that occur over a large range of spatio-temporal scales (Walker and Walker 1991). At a local scale, natural disturbance regimes include annual and daily freeze-thaw cycles, growth and erosion of ice-wedges, annual snow and associated water run-off, and snowbank formation and melting (Walker and Walker 1991). Intense frost action processes (congeliturbation) create patterned ground and landforms typical of Arctic environments and heavily influence the soil and vegetation of Arctic regions (Bliss 1962; Anderson and Bliss 1998). In such extreme environments, only a few species of vascular plants, bryophytes and lichens can establish (Bliss 1962; Billings 1987). Biological soil crusts consisting of soil surface and subsurface associations of cyanobacteria, algae, lichen, bryophyte, microfungus, and bacteria are often the first organisms to colonize bare soil, thereby facilitating the growth of vascular plants by promoting high moisture and nutrient levels (Belnap and Lange 2003; Breen and Lévesque 2006). Early vascular colonizers are often herbaceous plants such as graminoids and forbs that are very tolerant to adverse conditions and can produce a relatively large number of seeds (Lévesque 1997; Jones and Henry 2003; Breen and Lévesque 2006; O’Kane 2018). Long-lived, slow-growing woody species, such as Salix arctica, Dryas integrifolia, and especially Cassiope tetragona, often arrive in later successional stages (Lévesque 1997; Jones and Henry 2003; Breen and Lévesque 2006; O’Kane 2018). Plant succession in the High Arctic can occur over hundreds or even thousands of years (Billings and Peterson 1980; Svoboda and Henry 1987) and can be of different types. Along glacier forelands near polar oases, directional succession with species replacement generally occurs through recognizable stages. In more extreme polar desert environments, directional succession without species replacement arises when invading species become established and slowly expand over time (but without changes in community composition), whereas non-directional succession without species replacement occurs when few species succeed at establishing during favorable periods, although their populations are reduced or eradicated when conditions deteriorate (Svoboda and Henry 1987; Jones and Henry 2003). However, plant succession in the High Arctic is highly variable in space and time, and may take several trajectories depending on individual species traits and site-specific environmental conditions, such as temperature, moisture, micro-topography and substrate, leading to the emergence of different plant communities (Jones and Henry 2003; O’Kane 2018).

In order to fill geographic gaps in data coverage in Canada (Bjorkman et al. 2020) while improving harmonization of data collection and monitoring (Christensen et al. 2020), a systematic, plot-based survey of vascular plants was conducted at Alert (Ellesmere Island, Nunavut, Canada) in 2018–2019 by Desjardins et al. (2021a). This survey estimated the total number of vascular species at 77 in the Alert region and recorded no introduced species. Following the detailed inventory of vascular plants at Alert, the objectives of the present study were to (1) classify vegetation into distinct plant communities, (2) characterize plant community attributes, including diversity (alpha diversity, rare species, community composition) and abundance (percent cover), as well as environmental variables associated with each community, and (3) establish a georeferenced baseline with permanent field markers allowing future replication of our work.

Material and methods

Study area

The 170-km² study area surrounds the Canadian Forces Station (CFS) Alert (82°30ʹN, 62°20ʹW), the northernmost permanently inhabited settlement on Earth, which is located on the north-eastern tip of Ellesmere Island, Nunavut, Canada. It is roughly delimited by the Lincoln Sea to the North and the boundaries of CFS Alert property in other directions. The local topography consists of rugged and undulating terrain with mountains (maximum of 525 m a.s.l.), hills, valleys and creeks. The overburden, ranging from 2.4 to 4 m thick, is composed of till and shattered rock filled with ice (Taylor et al. 1982), while the underlying bedrock is highly calcareous, composed of argillite with greywacke in some places (Smith et al. 2012). The area is also characterized by a thick (>600 m) layer of permafrost (Smith et al. 2012). Isostatic rebound following deglaciation raised the marine limit to 135 m above sea level (England 1976) and uncovered many beds of marine silts and clays mixed with masses of recent marine shells (Bruggemann and Calder 1953). A more detailed description of the study area can be found in Desjardins et al. (2021a).

The sun remains under the horizon from mid-October to late February and a 24-h sunlight period occurs from early April to early September. Summer mean air temperature average 3.4°C in July while annual precipitations average 158 mm (Government of Canada 2010). The growing season for plants is short (2.5 months), with a snow-free season limited to mid-June to early September.

The study area is part of subzone A, which is the coldest bioclimate subzone in the Canadian Arctic, with an average July temperature of 3°C (Christensen et al. 2013). This subzone is generally characterized by a polar desert landscape that is mostly barren with some lichen and moss cover, and a vascular plant cover <5% (Christensen et al. 2013). Alert, nevertheless, supports several herbivore populations, including the collared lemming (Dicrostonyx groenlandicus), Arctic hare (Lepus arcticus), muskox (Ovibos moschatus), peary caribou (Rangifer tarandus pearyi), and rock ptarmigan (Lagopus muta), as well as their associated predators such as the Arctic wolf (Canis lupus arctos), Arctic fox (Vulpes lagopus), ermine (Mustela erminea), snowy owl (Bubo scandiacus), and long-tailed jaeger (Stercorarius longicaudus).

Data collection

The systematic vegetation survey was carried out from 28 July to 24 August 2018, and from 3 July to 2 September 2019. A random stratified design was used, based on a habitat map described in Desjardins et al. (2021a). The habitat map was constructed from a multispectral satellite picture of the study area (40-cm resolution pan-sharpened multispectral picture from WorldView-2/3; taken on 2 and 13 August 2015) and consisted of six classes, including snow/water and five broad vegetation gradients classified into habitat types: barren ground, xeric, xeric-mesic, mesic, and wetland. A 2 × 2 km grid was overlaid on the study area and, for each grid cell, vegetation surveys were performed at ≥1 random location for each of the five habitat types present in the grid (Figure 1(a)). ArcGIS 10.6.1 (ESRI 2018) was used to complete these steps.

Figure 1. a. Location of 264 vegetation plots (squares) randomly distributed in five habitat types in the study area: barren ground, xeric habitat, xeric-mesic habitat, mesic habitat and wetland. Closed squares indicate locations of permanently marked plots. b. Set up of the vegetation plot with the five 1 × 1 m quadrats. c. Top view of a quadrat. d. Vegetation sampling by a scribe (left), and an observer (right). Barren ground, xeric habitat, xeric-mesic habitat, mesic habitat, and wetland correspond, respectively, to 12%, 33%, 32%, 15%, and 7% of the total land cover

Each vegetation survey corresponded to a vegetation plot consisting of five 1 m × 1 m quadrats, each located 5 m from a central point and at equal distance from one another (Bay 1998; Figure 1(b)). GPS coordinates were recorded over each quadrat. A quadrat consisted of a frame with a distended rope grid in two layers, forming 100 intersections per layer (Figure 1(c)). This followed the intersection quadrat protocol of the International Tundra Experiment (Walker 1996). Following a point-intercept method, the frame was placed horizontally above the vegetation using metal poles at each corner. At each rope intersection, a metal pin was lowered along the two rope layers and the first plant encountered was recorded. Vascular plants (live or dead) were identified to the species, whereas cryptogams (i.e., spore-bearing organisms also known as non-vascular plants; Belnap et Lange 2003) were identified as biological soil crust (black or white), blue-green algae (macroscopic sheet colonies dominated by Nostoc sp.), lichen, macrofungus (defined as fungus with visible sporocarp; Dahlberg and Bültmann 2013), or moss. Vascular plants that could not be identified in the field were sampled just outside the quadrat for further identification using a microscope, identification guides, and expert opinion. Species names are based on VASCAN (Brouillet et al. 2010; Desmet and Brouillet 2013). For vascular species and cryptogams, we also recorded if touched plants were alive or dead (Desjardins et al. 2021b). If no plant was touched by the lowered pin, the ground substrate was noted (bare soil or rock). Following Walker (1996), an index of absolute cover for each vascular species, cryptogam, and ground substrate was then calculated for each quadrat as the total number of times the category was touched by the pin at a rope intersection, divided by 100. All vascular species and cryptogams found inside the quadrats but not touched by the pin were identified and assigned a cover value of 0.5%. The measure of percent cover (absolute abundance) obtained using the ITEX point-intercept method is highly correlated with the relative above-ground biomass obtained by harvest methods (Jonasson 1988; Shaver et al. 2001; Mamet et al. 2016). Using a Panasonic FZ70 (resolution: 180 dots per inch), Samsung ST150F (resolution: 72 dots per inch), or an iPhone SE (resolution: 72 dots per inch), oblique color photographs of quadrat were taken from hand and archived as Joint Photographic Experts Group (JPEG)-file. For 250 quadrats from 50 vegetation plots, holes left in the ground after removing the metal poles from the quadrat corners were permanently marked using two 20-cm metal nails hammered into the ground at opposite corners of the quadrats. One nail was tagged using a small numbered metal plate. The 50 marked vegetation plots were spread across the five habitat types (5 marked plots in barren ground, 17 marked plots in xeric habitat, 8 marked plots in mesic habitats and 20 marked plots in wetlands). This permanent marking and the associated georeferenced archive Nordicana D (Desjardins et al. 2021b) allow finding of plots with a metal detector and replication of our quantitative vegetation survey with the intersection quadrat protocol.

Data analysis

Plant community classification was performed through hierarchical cluster analysis using the complete linkage method and Bray–Curtis dissimilarity. The Bray–Curtis dissimilarity (Motyka et al. 1950) was chosen over other dissimilarity indices because it is not biased by species richness and yields index values that are higher among sites without species in common (Legendre et al. 2013). The index is sensitive to differences in abundance between species, and abundant species weight more than rare species (Chao et al. 2006; Ricotta and Podani 2017). Vegetation plots were the sampling unit for the plant community classification; thus, the input data for the analysis was the absolute cover of vascular plant species and cryptogams (dead and live plants combined), and one combined ground category (rock and soil) averaged over the five quadrats in each vegetation plot and transformed with Hellinger. For cryptogams, black biological crust, often embedded in moss, was merged with the moss category, while white biological crust, known to be composed of lichen (Belnap et al. 2001), was merged with the lichen category. Therefore, cryptogams were grouped into four categories instead of six: blue-green algae, lichen (including white biological crust), macrofungus, and moss (including black biological crust). In the Hellinger transformation recommended for clustering or ordination of species abundance data (Rao 1995; Legendre et al. 2001), the relative abundance values are square-rooted, which reduces the differences between the highest and lowest abundance values (Borcard et al. 2018).

We used the indicator value index (IndVal) to determine the appropriate number of groups in the cluster analysis (Dufrêne and Legendre 1997). IndVal assesses the strength of the association between a single species and several groups of sites, each group representing here a plant community. It is based on two components, specificity and fidelity, where specificity is the mean abundance of the species in the target site group divided by the sum of the mean abundance values over all groups, and fidelity is the relative frequency of occurrence of the species inside the target site group (Dufrêne and Legendre 1997; De Cáceres and Legendre 2009). Specificity is highest when the species is only present in the group and fidelity is highest when the species is present in all sites of the group (Dufrêne and Legendre 1997). The untransformed absolute cover values (with uncombined cryptogams) were used to estimate IndVal and significance was determined using Monte Carlo tests with 999 permutations. The highest number of groups with at least one vascular species or cryptogam with a significant IndVal in each group was chosen. The cluster analysis and IndVal calculations were performed with R Software (R Core Team 2020) using the vegan and indicspecies packages. For each group, we calculated the vascular plant species richness as well as the mean cover percentage of each vascular species, cryptogam, and ground substrate. The frequency of occurrence (number of vegetation plots with occurrence divided by the total number of plots multiplied by 100) of each vascular species, cryptogam, and ground substrate was also calculated.

Canonical Correspondence Analysis (CCA) was used to identify primary relationships between cover data and environmental variables (Legendre and Legendre 2012). Environmental variables included: (1) Normalized difference vegetation index (NDVI) computed from the pan-sharpened multispectral satellite picture of the study area, (2) elevation from the Canadian Digital Elevation Model (DEM) (Natural Resources Canada 2015), (3) slope derived from the DEM using ArcGIS, (4) distance to the nearest perennial snowbank derived from the habitat map (Desjardins et al. 2021a), (5) location with respect to a perennial snowbank (binary state with 1 indicating location below a snowbank and 0 location above or beside a snowbank); this binary variable was added since the distance to the nearest perennial snowbank may not fully capture the effect of snowbanks, considering that sites spatially close to a snowbank but located upslope or besides the snow patch do not actually benefit from melt water that flows downslope, (6) distance to the nearest lake or pond derived from the habitat map and from the Topographic Data of Canada – CanVec Series (Natural Resources Canada 2017), (7) distance to the nearest river derived from Topographic Data of Canada – CanVec Series (Natural Resources Canada 2017), and (8) distance to the nearest ocean shore derived from Topographic Data of Canada – CanVec Series (Natural Resources Canada 2017). These eight environmental variables were not collinear, as shown by the variance inflation factor (VIF; regclass package in R software, R Core Team 2020) being <10.

Multivariate analysis of variance based on Welch MANOVA, test (Hamidi et al. 2019), with 999 permutations, followed by pairwise comparisons, $T_w^2$ test (Alekseyenko 2016), were performed on species cover data and environmental variables separately to test for differences between communities using R software (R Core Team 2020). The multivariate analysis W^*_d has the advantage, over permutational multivariate analysis of variance (PERMANOVA), of being robust to unbalanced designs and to heteroscedasticity in the data (i.e., differences in multivariate group dispersion; Hamidi et al. 2019), which was our case.

For all statistical analyses, we considered the significance level at p < 0.05.

Results

Overview of floristic composition

A total of 264 vegetation plots were sampled. Of the 77 vascular plant species ever recorded in the study area (Desjardins et al. 2021a), 58 were found during this study, as well as five cryptogams (biological soil crust, blue-green algae, lichen, macrofungus, and moss). Among the 58 vascular species, 54 were found within the vegetation plots while four rare species were found only through opportunistic sampling (i.e., collected along marked trails and while traveling within the study area, outside of vegetation plots), namely Carex aquatilis Wahlenberg var. minor Boott, Eriophorum scheuchzeri Hoppe subsp. arcticum Novoselova, Puccinellia phryganodes (Trinius) Scribner & Merrill subsp. neoarctica (Á. Löve & D. Löve) Elven, and Saxifraga tricuspidata Rottbøll. Desjardins et al. (2021a) detail the floristic composition of the study area.

The five most commonly found vascular species, each encountered in >70% of the vegetation plots, were Saxifraga oppositifolia Linnaeus subsp. oppositifolia, Draba corymbosa R. Brown ex de Candolle, Stellaria longipes Goldie subsp. longipes, Cerastium arcticum Lange, and Papaver dahlianum Nordhagen (Appendix B). Fourteen species were found in ≤10 vegetation plots, namely Carex fuliginosa Schkuhr, Draba lactea Adams, Draba pauciflora R. Brown, Equisetum arvense Linnaeus, Equisetum variegatum Schleicher ex F. Weber & D. Mohr subsp. variegatum, Festuca baffinensis Polunin, Festuca hyperborea Holmen ex Frederiksen, Festuca viviparoidea Krajina ex Pavlick subsp. viviparoidea, Pleuropogon sabinei R. Brown, Poa arctica R. Brown subsp. arctica, Poa pratensis Linnaeus subsp. colpodea (Th. Fries) Tzvelev, Ranunculus hyperboreus Rottbøll, Ranunculus sabinei R. Brown, and Ranunculus sulphureus Solander (Appendix B). When added to the four species found only through opportunistic sampling, this gives a total of 18 rare vascular plant species found in the study area.

Description of plant communities

Five plant communities (Table 1) were identified using cluster analysis (Figure 2) and indicator species values (Table 2). The multivariate analysis of variance $W_d^{\ast }$ and pairwise comparisons $T_w^2$ confirmed a significant difference in the cover data among each of the five communities ($W_d^{\ast }$ statistic = 99.03; p = 0.001).

Table 1. Dominant vegetation, main habitat class, vascular plant species richness and mean percentage cover of ground, cryptogams, and vascular species groups for the five plant communities identified at Alert (Nunavut, Canada). Habitat classes correspond to the dominant vegetation classes identified by Desjardins et al. 2021a. Vascular species groups correspond to the Focal Ecosystem Components identified by Christensen et al. (2013). Community name was based on Bridgland (1986), where biotic cover <40% is associated with barrens, while biotic cover ranging between 40% and 70% is associated to tundras. We arbitrary associated biotic cover ≥70% to wetlands. No distinction was made among the wetland communities on whether the substrate was saturated with standing water during the entire growing season (marsh) or if it was drained and remained moist (meadow)

Plant community		Dominant vegetation	Main habitat class	Vascular species richness	Cover (%)
					Soil/Rock	Biological sol crust	Lichen	Moss	Algae/ Macrofungus	Graminoid	Forb	Shrub
I	Forb-dominated tundra	Saxifraga oppositifolia subsp. oppositifolia Moss Stellaria longipes subsp. longipes	Xeric-Mesic	46	57.2	1.3	0.7	8.7	0.1	6.5	25.0	2.1
II	Forb-dominated barren	Saxifraga oppositifolia subsp. oppositifolia Salix arctica Moss	Barren	35	88.0	0.2	1.5	1.6	0	2.1	5.6	1.8
III	Sedge-dominated wetland	Eriophorum triste Moss Salix arctica	Wetland	40	4.0	0.2	< 0.1	20.5	0.1	58.2	7.6	10.3
IV	Moss-dominated wetland	Moss Saxifraga cernua Luzula nivalis	Wetland	32	3.7	4.1	0.5	53.0	0.3	15.7	24.2	0.4
V	Grass-dominated wetland	Moss Alopecurus magellanicus Juncus biglumis	Wetland	50	21.0	3.6	0.2	22.3	0.5	35.8	13.9	4.3

Table 2. Significant indicator values (IndVal) and their associated specificity (Spec.) and fidelity (Fid.) components, for five plant communities and their indicator species identified at Alert (Nunavut, Canada) using an indicator species analysis. Indicator species (including vascular species and cryptogams) are ordered by decreasing IndVal within each community. * p < 0.05, ** p < 0.01, *** p < 0.001

Plant community	Indicator species	Family	IndVal	Spec.	Fid.	p
I	Stellaria longipes subsp. longipes	Caryophyllaceae	0.75	0.57	0.99	0.001 ***
	Saxifraga oppositifolia subsp. oppositifolia	Saxifragaceae	0.72	0.53	0.99	0.001 ***
	Sabulina rubella	Caryophyllaceae	0.71	0.54	0.94	0.001 ***
	Puccinellia bruggemannii	Poaceae	0.69	0.64	0.74	0.004 **
	Cerastium arcticum	Caryophyllaceae	0.68	0.48	0.95	0.003 **
	Draba corymbosa	Brassicaceae	0.58	0.36	0.96	0.026 *
II	Braya thorild-wulffii subsp. thorild-wulffii	Brassicaceae	0.42	0.65	0.29	0.048 *
III	Eriophorum triste	Cyperaceae	0.99	0.99	1.00	0.001 ***
	Arctagrostis latifolia subsp. latifolia	Poaceae	0.95	0.97	0.92	0.001 ***
	Bistorta vivipara	Polygonaceae	0.84	0.86	0.81	0.001 ***
	Salix arctica	Salicaceae	0.73	0.55	0.96	0.002 **
	Equisetum variegatum subsp. variegatum	Equisetaceae	0.43	0.97	0.19	0.007 **
	Pedicularis hirsuta	Orobranchaceae	0.42	0.45	0.38	0.044 *
	Carex fuliginosa	Cyperaceae	0.37	0.91	0.15	0.048 *
IV	Micranthes tenuis	Saxifragaceae	0.94	0.88	1.00	0.001 ***
	Saxifraga cernua	Saxifragaceae	0.93	0.86	1.00	0.001 ***
	Cardamine bellidifolia	Brassicaceae	0.87	0.93	0.80	0.001 ***
	Phippsia algida	Poaceae	0.79	0.78	0.80	0.001 ***
	Saxifraga cespitosa	Saxifragaceae	0.76	0.58	1.00	0.001 ***
	Papaver dahlianum	Papaveraceae	0.73	0.53	1.00	0.001 ***
	Moss		0.71	0.50	1.00	0.001 ***
	Black biological soil crust		0.68	0.46	1.00	0.006 **
	Ranunculus sulphureus	Ranunculaceae	0.62	0.97	0.40	0.001 ***
	Luzula nivalis	Juncaceae	0.59	0.87	0.40	0.008 **
	Draba micropetala	Brassicaceae	0.51	0.43	0.60	0.040 *
	Poa arctica R. Brown subsp. arctica	Poaceae	0.45	1.00	0.20	0.022 *
	Festuca hyperborea	Poaceae	0.42	0.89	0.20	0.016 *
V	Alopecurus magellanicus	Poaceae	0.85	0.76	0.96	0.001 ***
	Juncus biglumis	Juncaceae	0.66	0.45	0.97	0.004 **
	Puccinellia vahliana	Poaceae	0.64	0.69	0.59	0.007 **
	Cerastium regelii	Caryophyllaceae	0.60	0.44	0.81	0.017 *
	×Pucciphippsia vacillans	Poaceae	0.40	0.72	0.22	0.035 *

Figure 2. Five community types (I to V) revealed by a hierarchical cluster analysis of 264 vegetation plots, based on the complete linkage method and Bray–Curtis dissimilarity index with Hellinger transformed-data. Colors under each plot refer to the five vegetation classes described in Desjardins et al. (2021a)

When examining how the five habitat map classes were distributed among the plant communities (Figure 2), we observed that 100% of barren and 64% of xeric sampled habitats were associated with Community II (Figure 3(a)), while 79% of xeric-mesic, 44% of mesic and 36% of xeric sampled habitats were associated to Community I (Figure 3(b)). Wetlands were almost exclusively distributed among the three remaining communities, with 35% associated with Community III (Figure 3(c)), 7% with Community IV (Figure 3(d)) and 56% with Community V (Figure 3(e)). Community V also comprised 54% of mesic sampled habitats.

Figure 3. General appearance of five community types identified at Alert (Nunavut, Canada). a. Forb-dominated barren (Community II). b. Forb-dominated tundra (Community I). c. Sedge-dominated wetland (Community III). d. Moss-dominated wetland (Community IV). e. Grass-dominated wetland (Community V)

Community I had the second highest species richness (46) and an average mean biotic cover (44%, with a 34% vascular plant cover) when compared with other communities (Table 1; Figure 3(b)). Cover was mostly represented by Saxifraga oppositifolia subsp. oppositifolia (11%), moss (9%), Stellaria longipes subsp. longipes (5%), Cerastium arcticum (3%), Puccinellia bruggemannii T.J. Sørensen (2%), Salix arctica Pallas (2%), Draba corymbosa (1%), Papaver dahlianum (1%), Sabulina rubella (Wahlenberg) Dillenberger & Kadereit (1%), Juncus biglumis Linnaeus (1%), black and white soil biological crusts (1% each), lichen (1%), Puccinellia vahliana (Liebmann) Scribner & Merrill (1%), Poa abbreviata R. Brown subsp. abbreviata (1%), and Phippsia algida (Solander) R. Brown (1%; Appendix B). There were six indicator species, but the most significant ones included two chickweeds, Sabulina rubella and Stellaria longipes subsp. longipes, and one saxifrage, Saxifraga oppositifolia subsp. oppositifolia (Table 2). These three indicator species had a moderately high IndVal (0.71–0.75), mostly due to a very high fidelity (0.94–0.99; Table 2).

Community II had the second lowest species richness (35) and the lowest mean biotic cover (13%, with a 10% vascular plant cover; Table 1; Figure 3(a)). Cover was mostly represented by Saxifraga oppositifolia subsp. oppositifolia (3%), Salix arctica (2%), moss (2%), lichen (1%), Poa abbreviata subsp. abbreviata (1%), and Puccinellia bruggemannii (1%; Appendix B). The only indicator species significantly associated with this community is the crucifer, Braya thorild-wulffii Ostenfeld subsp. thorild-wulffii, which had moderate IndVal (0.42; Table 2). This indicator species had an absolute cover that was relatively low (0.3%), but still higher than in the other four communities (Appendix B).

Community III had an intermediate species richness (40) and the highest mean biotic cover (97%, with a 76% vascular plant cover; Table 1; Figure 3(c)). Cover was mostly dominated by Eriophorum triste (Th. Fries) Hadac & Á. Löve (38%), moss (20%), Salix arctica (10%), Arctagrostis latifolia (R. Brown) Grisebach subsp. latifolia (9%), Juncus biglumis (5%), Alopecurus magellanicus Lamarck (5%), Bistorta vivipara (Linnaeus) Delarbre (5%), Saxifraga oppositifolia subsp. oppositifolia (1%), and Equisetum variegatum subsp. variegatum (1%; Appendix B). There were seven indicator species, the most significant ones including one sedge, Eriophorum triste, one grass, Arctagrostis latifolia subsp. latifolia, and a knotweed, Bistorta vivipara (Table 2). Eriophorum triste, Arctagrostis latifolia subsp. latifolia, and Bistorta vivipara had high IndVal (0.84–0.99), due to both high specificity (0.86–0.99) and high fidelity (0.81–1.00; Table 2).

Community IV had the lowest species richness (32), but a very high mean biotic cover (98%, with a 40% vascular plant cover; Table 1; Figure 3(d)). Cover was mostly represented by moss (53%), Saxifraga cernua Linnaeus (6%), Luzula nivalis (Laestadius) Sprengel (6%), black biological soil crust (4%), Phippsia algida (3%), Juncus biglumis (3%), Papaver dahlianum (3%), Cardamine bellidifolia Linnaeus (2%), Alopecurus magellanicus (2%), Cerastium arcticum (2%), Saxifraga cespitosa Linnaeus (2%), Saxifraga oppositifolia subsp. oppositifolia (2%), Micranthes tenuis (Wahlenberg) Small (2%), Stellaria longipes subsp. longipes (1%), Cerastium regelii (1%), Draba corymbosa (1%), Ranunculus sulphureus (1%), and Poa arctica subsp. arctica (1%; Appendix B). There were 13 indicator species, the most significant ones being the moss, three saxifrages (Saxifraga cernua, Saxifraga cespitosa, and Micranthes tenuis), one crucifer (Cardamine bellidifolia), one grass (Phippsia algida), one poppy (Papaver dahlianum), and one buttercup (Ranunculus sulphureus; Table 2). Two of the saxifrages and the crucifer had the highest IndVal (0.87–0.94), mainly due to a high fidelity and specificity (Table 2). Phippsia algida, Saxifraga cespitosa, and Papaver dahlianum had relatively high IndVal (0.73–0.79) due to being specific to this community (specificity = 0.80–1.00; Table 2). The buttercup had an intermediate IndVal (0.62) with a high specificity (0.97; Table 2).

Community V had the highest species richness (50) and a high mean biotic cover (81%, with a 54% vascular plant cover; Table 1; Figure 3(e)). Cover was mostly represented by moss (22%), Alopecurus magellanicus (22%), Juncus biglumis (7%), Salix arctica (4%), Saxifraga oppositifolia subsp. oppositifolia (4%), black biological soil crust (3%), Puccinellia vahliana (2%), Stellaria longipes subsp. longipes (2%), Cerastium regelii (1%), Draba corymbosa (1%), Saxifraga cespitosa (1%), Papaver dahlianum (1%), Saxifraga cernua (1%), Cerastium arcticum (1%), ×Pucciphippsia vacillans (T. Fries) Tzvelev (1%), Sabulina rubella (1%), Puccinellia bruggemannii (1%), and blue-green algae (1%; Appendix B). There were five indicator species, the most significant one being one grass, Alopecurus magellanicus, which presented a relatively high IndVal (0.85) and a high specificity (0.96; Table 2).

Several rare species (fidelity <0.25) were highly restricted to one community (III or IV). This included Equisetum variegatum subsp. variegatum (specificity to III = 0.97), Carex fuliginosa (specificity to III = 0.91), Poa arctica subsp. arctica (specificity to IV = 1.00), and Festuca hyperborea (specificity to IV = 0.89; Table 2).

Environmental characteristics of plant communities

Based on the multivariate analysis of variance $W_d^{\ast }$, there was a significant difference in the measured environmental variables among the five communities ($W_d^{\ast }$ statistic = 99.03; p = 0.001). The pairwise comparisons $T_w^2$ revealed that all communities were different in terms of environmental characteristics, except for communities I and V that were not different from one another ($T_w^2$ statistic = 1.40; p = 0.19). The overall model including the eight environmental variables was significant (F_8,255 = 9.38, p = 0.001). However, only NDVI, elevation, distance to a snowbank, location with respect to a snowbank, and distance to the ocean shore were significantly related to the cover values of vascular species, cryptogams and ground substrates (Appendix D). Permutation tests with 999 runs in the CCA indicated that the first three canonical axes were significant (p = 0.001 for CCA1, CCA2, and CCA3). The eigenvalues for the first and second axes of the species-environment plot were 0.56 and 0.11, respectively, representing 16.18% and 3.13% of the total variation in the dataset. NDVI was most highly correlated with the first axis (0.99), elevation with the second axis (0.79), distance to the ocean shore with the second axis (0.68), and distance to a snowbank as well as location with respect to a snowbank were correlated with the third axis (−0.66 and 0.61, respectively; CCA3 not shown).

The CCA ordination biplot showed that plant communities were rather well separated (Figure 4). Communities were arranged from highest to lowest NDVI values on the first axis, with III, IV, and V associated with higher NDVI values (right side of the plot) and I and II associated with lower NDVI values (left side of the plot). Sites located at higher elevations were associated with Community IV. Community III sites were the most associated with location with respect to a snowbank, being therefore mainly located downslope of a snowbank. Community III sites were also negatively correlated with elevation, distance to snowbanks, distance to lakes and ponds, and distance to ocean, as corroborated by their low average elevation (55 m), relatively short distance to snowbanks (351 m), relatively short distance to lakes and ponds (855 m), and relatively short distance to ocean shores (1219 m; Appendix C).

Figure 4. Relations between cover of vascular plant species, cryptogams and ground substrates, and eight environmental variables, as revealed by a Canonical correspondence analysis. Dots represent sampling sites, and colors identify the five plant communities revealed by an indicator species analysis. Arrows represent the environmental variables, with bold variables being significant

Marking of vegetation plots

Despite the a posteriori classification of plant communities, the permanent marking of 250 quadrats (50 vegetation plots) spread across five habitats allowed permanent identification of ≥10 quadrats (two plots) in each plant community (Table 3). As a result, each plant community is represented in the pool of permanently marked plots by ≥1% of the surveyed plots. Not surprisingly, communities I and II, which appeared the most in surveyed plots, are also the most represented in permanently marked plots. Conversely, Community IV appeared the least in surveyed plots and is also the least represented in permanently marked plots.

Table 3. Number of quadrats (and plots) permanently marked in each of the five plant communities identified at Alert, Nunavut (Canada). The number of quadrats (and plots) surveyed is also indicated for each community, as well as the proportion of marked quadrats with respect to all quadrats surveyed in each community and in the study area

Plant community	I	II	III	IV	V	Total
Marked quadrats (plots)	70 (14)	70 (14)	40 (8)	10 (2)	60 (12)	250 (50)
Surveyed quadrats (plots)	390 (78)	435 (87)	130 (26)	25 (5)	340 (68)	1320 (264)
Marked quadrats in proportion to all quadrats surveyed in the community	17.9%	16.1%	30.8%	40.0%	17.6%
Marked quadrats in proportion to all quadrats surveyed in the study area	5.3%	5.3%	3.0%	0.8%	4.5%	18.9%

Discussion

Following an important field effort to survey vegetation in five habitats identified from a multispectral satellite picture of the Alert study area, we successfully identified five plant communities and assessed their correlations with eight environmental variables. We first compare the rare and frequent species found at Alert with those found at other sites, we then compare the identified plant communities with each other and in the context of previously published work (mostly in the Canadian Arctic), we discuss their environmental characteristics, and we finally present the benefits of establishing a baseline reference for the evaluation of future changes in High-Arctic vegetation.

Overview of floristic composition

The five most common vascular plant species at Alert (Saxifraga oppositifolia subsp. oppositifolia, Draba corymbosa, Stellaria longipes subsp. longipes, Cerastium arcticum, and Papaver dahlianum) are among the 25 most frequent species at the circumpolar scale within subzone A (Daniëls et al. 2016). Cerastium arcticum is also a diagnostic species associated with patterned grounds in subzone A (Walker et al. 2011). Except Stellaria longipes subsp. longipes, the remaining common species at Alert were also among the 10 most commonly found by Lévesque (1997) in polar desert sites in central Ellesmere Island. These widespread Arctic species are highly tolerant and effective colonizers of isolated and poor sites, due to their flexible metabolic rate, growth allocation, and type of reproduction (Lévesque 1997; Peterson 2014).

Among the 18 rare species that we recorded, Poa arctica subsp. arctica, Saxifraga tricuspidata, Draba lactea, Carex fuliginosa (previously called Carex misandra), and Festuca hyperborea were also found by Lévesque (1997) in central Ellesmere Island, but only the three latter species were considered as rare. Two of the rare species, Carex aquatilis var. minor and Eriophorum scheuchzeri, are dominant species in more southern marshes and meadows such as Cape Herschel on Ellesmere Island (Bridgland 1986) and Polar Bear Pass on Bathurst Island (Sheard and Geale 1983). The two Equisetum Linnaeus species are more abundant in moist to wet, coarse mineral soils at lower latitudes on Ellesmere Island, such as at Lake Hazen and Okse Bay (Hodson et al. 2008). In addition, Alert is currently the northernmost reported location on Ellesmere Island for six of the rare species that we found (Festuca baffinensis, Festuca hyperborea, Poa arctica subsp. arctica, Puccinellia phryganodes, Ranunculus sabinei, and Saxifraga tricuspidata), as well as the only reported location on Ellesmere Island of Festuca viviparoidea subsp. vivparoidea (Aiken et al. 2007; GBIF 2020), a species which is listed as imperiled in Nunavut by NatureServe (2020). Interestingly, we found Puccinellia phryganodes strictly on the ocean shore. Stolons of this salt marsh grass embedded in sea ice have reportedly been recovered and grown, providing a potential mechanism for the distribution of this species across the Arctic (Peterson 2014). Rare plants are often ecological specialists that occupy microhabitats where competition with zonal vegetation is absent or limited (Talbot et al. 1999).

Plant communities

The five communities showed clear differences with each other in terms of dominant species, biotic cover, species richness, and associated habitats. Despite the relatively sparse vegetation characterizing the Alert region, mosaics of diverse communities within the landscape have also been reported in other polar deserts (Crawford 2008; Peterson 2014; Daniëls et al. 2016; Zwolicki et al. 2020).

Barren habitats – Community II, here called forb-dominated barren and dominated by Saxifraga oppositifolia subsp. oppositifolia, Salix arctica, and moss, had the lowest mean biotic cover (13%), but surprisingly not the lowest species richness. The sites associated with Community II were found in five types of barren habitats: small/medium non-sorted polygons made of till with or without frost-shattered rocks on the surface, boulder/rock fields, riverbeds with gravel, saline areas near the ocean shore, and very few desert pavements. The non-sorted polygon landform was the most widespread and vegetation was concentrated in the cracks delimiting the polygons, which allow seedling establishment (Bell and Bliss 1980; Sohlberg and Bliss 1984; Elberling 2000) through protection from desiccation, exposure to wind and soil disturbance (Lévesque 2001; Walker et al. 2011).

Community II was similar in terms of vascular plant species to the plateau above Truelove lowland on Devon Island (Bliss et al. 1994), to the ridge community found in Polar Bear Pass (Sheard and Geale 1983), and to the cushion plant community in polar deserts on six islands in the central Canadian Arctic Archipelago (Bliss et al. 1984), but their respective biotic cover compared to Alert is either higher (52% in Polar Bear Pass, a polar oasis) or lower (5% on Devon Island plateau, a polar desert, and 7% in polar deserts on six islands). Community II had many species in common with these three communities, such as Saxifraga oppositifolia, Salix arctica, Puccinellia angustata (R. Brown) E.L. Rand & Redfield, Draba corymbosa, Cerastium sp. Linnaeus, Papaver sp. Linnaeus, Sabulina rubella (previously called Minuartia rubella) (Sheard and Geale 1983; Bliss et al. 1984, 1994). Puccinellia angustata was present in Community II, but it was rather Puccinellia bruggemannii, one of its close relatives, that was encountered more often. As for Papaver, the previous studies reported P. radicatum Rottbøll, but this name is strictly associated with European species whereas three Papaver species have been recognized recently on Ellesmere Island, including P. dahlianum (Solstad 2008) that we found at Alert. As for Cerastium, we found the arcticum subspecies instead of C. alpinum Linnaeus (a very close relative). C. alpinum and C. arcticum have traditionally been lumped into the alpinum subspecies due to overlapping characters (Hultén 1956; Böcher 1977; Brysting and Elven 2000).

When we subdivided Community II into two sub-communities by selecting a lower threshold in the dendrogram (dissimilarity = 0.55), we observed that the first sub-community (5% biotic cover) was similar to cushion plant barrens found by Bridgland (1986) at Cape Herschel (a polar semi-desert) by having Saxifraga oppositifolia and Poa abbreviata as dominant species, with a high abundance of lichen. The second sub-community (20% biotic cover) was similar to the dwarf shrub barrens also found at Cape Herschel by having Saxifraga oppositifolia and Salix arctica as dominant species.

Community II had therefore many dominant species in common with the other barren communities of Ellesmere Island, Prince of Wales, Somerset, Devon, Cornwallis, Bathurst, and King Christian islands (Sheard and Geale 1983; Bliss et al. 1984, 1994; Bridgland 1986), but the biotic cover was mostly higher at Alert than in other polar deserts. Moreover, a recurring difference between Alert and the other polar desert sites was the higher cover of Phippsia algida and Luzula nivalis in the barren communities, whereas these species were only abundant in a wetland community in Alert.

Mesic habitats – Community I, here called forb-dominated tundra, was mostly found in xeric-mesic and mesic habitats, and had the highest species richness and a moderate biotic cover (44%) and vascular plant cover (34%). As for Community II, Community I was growing mostly on non-sorted patterned grounds made of till, but showed a higher biotic cover spreading from the cracks. Sometimes, vegetation was covering the entire cracks giving a flat appearance to the ground with a bare till circle in the centre (landforms also known as non-sorted circles). Some sites associated with this community also had hummocks as landforms. As a result of an interaction between physical processes and biotic properties of vegetation, patterned-ground changes occur during vegetation primary succession (Shur et al. 2008; Walker et al. 2011). When vegetation starts colonizing small non-sorted polygon cracks, the soil temperature differential between the centre and margins of the polygons increases the thaw-layer thickness in the centre of the polygons (Shur et al. 2008). This increases the subsurface flow and pooling of water and organic matter basins beneath the polygons (Shur et al. 2008). Then, small cryptogamic and forb species followed by thicker moss mats and more diverse species start colonizing the polygon margins and spread over the active centres. At this point, the non-sorted polygons become non-sorted circles. As the vegetation and soil organic matter build in the centres of the non-sorted circles, the active layer is reduced in thickness, and aggradational ice (new ice added to the top of the permafrost table) develops beneath the polygons. This causes the centres of the circles to become permanently elevated relative to the cracks, forming the earth mounds, called hummocks (Walker et al. 2011). Community I may therefore represent a more advanced succession stage than Community II, with well-drained soils. Indeed, Community I presented the same ubiquitous species as Community II, namely Saxifraga oppositifolia subsp. oppositifolia, Stellaria longipes subsp. longipes, Cerastium arcticum, Puccinellia bruggemannii, Salix arctica, Draba corymbosa, Papaver dahlianum, and Sabulina rubella, although with a higher cover. In addition, three other dominant species mostly associated with moister soils were also present in the Community I: Juncus biglumis, Phippsia algida, and Puccinellia vahliana (Aiken et al. 2007).

In terms of species, very few communities similar to Community I have been documented in the High Arctic. The most similar ones are a hummock community reported in Polar Bear Pass (Sheard and Geale 1983) and the dwarf shrub tundra in Cape Herschel (Bridgland 1986). The hummock community (24% mean vascular plant cover and 72% mean biotic cover) has 10 out of 12 dominant species in common with Community I (Saxifraga oppositifolia, Juncus biglumis, Stellaria longipes, Papaver sp. (previously recorded as P. radicatum), Salix arctica, Puccinellia bruggemannii, Draba corymbosa (previously called Draba bellii Holm), Poa abbreviata, and Cerastium alpinum; Sheard and Geale 1983), while the dwarf shrub tundra (63% mean biotic cover) has four dominant species (Salix arctica, Saxifraga oppositifolia, Draba sp. Linnaeus, and Stellaria longipes; Bridgland 1986). As observed by Sheard and Geale (1983), this mesic community with hummock landforms had the highest vascular plant cover of the upland area, with the exception of meadows and marshes. The mean biotic cover of Community I was, however, lower compared to the hummock and the dwarf shrub tundra communities found, respectively, in a polar oasis and a polar semi-desert.

Wet habitats – The three remaining communities, mostly wetlands, had very high mean biotic covers (81–98%), and differed among each other in their dominant vegetation. Community III (Sedge-dominated wetland) was mostly dominated by Eriophorum triste (38%), while Community IV (Moss-dominated wetland) was dominated by moss (53%), and Community V (Grass-dominated wetland) by Alopecurus magellanicus (22%). The indicator species in Community III and IV were very specific, which also help to distinguish the three communities in the field. These communities resembled polar oases commonly found in more southern latitudes and dominated by graminoids (Bliss 1975, 1981). The wetland communities presented here were similar, in terms of dominant species, to the snowflush communities described by Bliss et al. (1984). They discerned different patterns within the snowflush communities, including those dominated by moss, Eriophorum triste, or Alopecurus magellanicus (previously called Alopecurus alpinum Smith), which correspond to the three wetland communities found in this study. Like Bliss et al. (1984), we observed that mosses were abundant but lichens were uncommon, and Salix arctica was present in some sites. The snowflush communities are associated with late-lying snow, while in Alert the wetland communities were associated with perennial snowbanks, ponds, lakes, streams, and probably with underground water. The biotic cover reported by Bliss et al. (1984) was, however, lower (36%), suggesting favourable hydrological regimes at Alert, since water availability is necessary for plant communities to develop relatively extensive and lush vegetation (Muc et al. 1989; Le Roux et al. 2013; Nabe-Nielsen et al. 2017).

Community III is very similar to the cotton grass meadows of Cape Herschel. The mean biotic cover there was very high (79%), reaching 100% at some sites, and the community was dominated by two similar species: Eriophorum triste and Salix arctica (Bridgland 1986).

Community IV is similar to two communities: the marsh community found at Cape Herschel (Bridgland 1986), and the mossy snowbed in Amdrup Land, a semi-polar desert in Greenland (Bay 1997). The vascular plant species richness was lower at Cape Herschel, but 9 out of 15 vascular species present were found in Community IV, including few rare grasses and buttercup: Ranunculus sulphureus, Poa arctica, and Pleuropogon sabinei (Bridgland 1986). The biotic cover (93%) and the number of vascular plant species (25 spp.) reported by Bay (1997) was similar to the one reported here (96% and 32 spp.), and we share half of the dominant species found in the mossy snowbed community (7 out of 13 vascular species). This type of community, also known as wet moss tundra, is commonly described in the literature as a peatland dominated by a thick and dense moss layer (Vanderpuye et al. 2002). Usually, graminoids and forbs, such as Carex aquatilis var. minor, Pedicularis hirsuta, and Pleuropogon sabinei which are later successional species, are emerging from the moss and their roots are attached directly to the moss layer. Wet moss tundra is usually formed in moist slopes and concave terrains and relies on nutrient inputs from seabird and large herbivore feces (Nakatsubo et al. 2015). It can also be restricted to the borders of permanent ponds and streams (Bridgland 1986). In the Alert region, the moss-dominated wetland occurred in different areas, including gentle slopes at the foot of a mountain, along a river stream, and on the margins of a pond or a lake. Peary caribou feces were also very frequently observed in this community (Desjardins E., unpublished data).

Community V is similar to the moss graminoid meadows, dominated by Alopecurus magellanicus (previously called Alopecurus alpinus) found in valleys on King Christian Island (Bliss and Svoboda 1984). At Alert, 18 sites associated with this community had a lower cover (mean biotic cover 47–70%) with drier soils and some seemed to be in transition toward luxuriant vegetation or others in a retrogression toward the Community I (due to high number of dead stems from previous years).

We report higher biotic cover in the communities at Alert compared to other polar deserts in the Canadian Arctic Archipelago (Bliss et al. 1984, 1994; Lévesque 1997). This might be partly explained by differences in methodology; we used the point-intercept method, while previous studies were based on visual estimations of cover. We reported higher vascular species richness than previous studies in similar barren communities, which could have led us to report higher cover as we attributed a 0.5% cover to every species present in the quadrat but not hit at one of the 100 rope intersections. Heterogeneity in regional plant communities can also result from local variations in soil attributes (moisture, pH, organic matter content, soil chemistry, soil texture, active-layer dynamics), topography (elevation, water movement, as well as sun, wind and snow exposure), and time since the site became free of permanent snow or ice (Walker 1995; Lévesque 1997; Christensen et al. 2013; Nabe-Nielsen et al. 2017; Zwolicki et al. 2020). Indeed, at Alert, barren grounds made of till have a high proportion of clay and many coastal sites are clay-rich soils; the small particle size of clay favours a better water retention (Liptzin 2006) and consequently supports plant colonization and growth (Nabe-Nielsen et al. 2017). Alert was also among the first regions to be free of the Innuitian Ice Sheet on Ellesmere Island (Dalton et al. 2020), leaving more time for plant development.

Successional pathways

As highlighted by Bliss et al. (1984), Yurtsev (1994) and Daniëls et al. (2016), polar deserts are not always dry and sparsely vegetated, but some areas with limited extent can be very moist and rather continuously vegetated. We may be observing different successional stages at Alert. The early one corresponds to the first sub-community of the ford-dominated barren (Community II), which was dominated by few purple saxifrages, grasses, and lichen. The second sub-community, where expansion of the first-subcommunity species and colonization by new tolerant species (e.g., crucifers) led to a biotic cover increase, may represent the next successional stages. The Arctic willow, Salix arctica, was also present at some sites in Community II, but became more abundant in the forb-dominated tundra (Community I), indicating a stability of 250 years (Cray and Pollard 2015). Additional species, mostly grasses, chickweeds, and saxifrages, also colonized this mesic community. These species were also present in the more advanced successional stages that could be the wetland communities, but were supplanted by highly dominant, moisture-associated species (sedge, grass, rush, and moss). The moss-dominated wetland may be the community in the study area that reached the higher successional stage owing to peatland below the thick moss layer serving as a substrate for forb and graminoid species that were not encountered anywhere else in the study area (e.g., buttercups and Pleuropogon sabinei). Based on our observations, we suggest that primary succession at Alert corresponds to directional succession without species replacement in drier areas, and to directional succession with species replacement where favourable conditions prevail (e.g., constant water sources and fertilization such as animal feces). Radiocarbon dating of macrofossils and paleosol analysis would be necessary to determine the actual chronosequences (Bliss and Gold 1994) and confirm these hypotheses.

Environmental characteristics

The ordination biplot showed that communities differed regarding NDVI, elevation, distance to the ocean shore, and location with respect to the nearest snowbank (distance and downslope positioning). NDVI has been extensively used in Arctic studies as an index of vegetation ‘greenness’ (Stow et al. 2007; Walker et al. 2012; Berner et al. 2020), as a proxy for net primary productivity (Goward and Dye 1987), or to reflect biophysical properties such as biomass, leaf area index, and CO₂ flux (Hope et al. 1993; Gower et al. 1999; McMichael 1999). Communities III, IV and V were characterized by higher NDVI values than other communities, reflecting their higher proportion of vascular plant and moss cover, and probably a higher productivity. Hope et al. (1993) reported that NDVI could be community specific, with variance in the NDVI being mainly related to photosynthetic biomass and vegetation cover. Distinct spectral signatures can be used to detect vegetation communities over larger areas based on remotely sensed data (Hope et al. 1993; Stow et al. 2004). Further analyses of spectral signatures validated by ground truthing are, however, still needed at Alert. Wetland habitats represented only 7% of the landscape in the study area (Figure 1), but their rich, productive plant communities may be important for herbivores (Bliss et al. 1984).

Elevation encompasses multiple factors that influence vegetation, including microclimate, exposure, timing of snowmelt, length of the growing season and seed availability (Lévesque 1997). Upland areas in the High Arctic are generally the most barren (Bliss et al. 1984, 1994; Lévesque 1997). Although the elevation gradient was limited from 0 to 365 m a.s.l. among the surveyed plots, plant communities in wet habitats were segregated according to elevation. All sites of Community IV, and some sites of Community V, found at higher elevations probably benefit from relatively favorable temperature and moisture conditions due to local topography (mountain slopes, valleys) or a source of water (snowmelt, underground water), as also seen in other upland polar desert areas in the Canadian Arctic Archipelago (Bliss et al. 1984, 1994). Sites of Community III were mainly located in lowlands as they were associated with lower altitude and closer to ocean shores. Sites of Community III were also the ones most associated with snowbanks; the rolling lowlands may favor the formation of perennial snowbanks.

Surprisingly, if we except the variables related to snowbanks and ocean, the other variables related to a source of water (distance to a river and lake/pond) were not significant in the CCA. Distance to a waterbody may not correlate with actual soil moisture; thus, accurately measuring soil water content with a probe will be important in future vegetation studies in the High Arctic. It is also possible that the wetland communities are mostly ground water-fed wetlands. These wetlands are found in topographical depressions and slope concavities favoring drainage input and water retention, and are known to be the most predominant wetlands in High Arctic (Woo and Young 2006). However, our results highlight the importance of perennial snowbanks for sustaining some vegetation communities at Alert. As stated previously, soil moisture plays a key role in tundra community composition (Le Roux et al. 2013; Nabe-Nielsen et al. 2017). Snowbanks gradually melting during the summer create a gradient of increasing water availability closer to the remaining snow patch (Billings and Bliss 1959; Woo and Young 2014). The perpetual inflow of melt water to areas downslope of snowbanks influences the richness and composition of tundra plant communities, as well as plant growth and productivity (Billings and Bliss 1959; Canaday and Fonda 1974). High soil moisture present until the end of the growing season favors species that have a low tolerance to drought (Billings and Bliss 1959). Community III may therefore be the most vulnerable to a potential loss of snowbanks due to climate warming. Shrinkage or disappearance of snowbanks in the High Arctic, and subsequent modifications of thermal and hydrological conditions causing shifts in vegetation, have already been reported in Canada (Young and Woo 2003; Woo and Young 2014). Increased summer precipitations linked to climate change (Collins et al. 2013) may compensate for the loss of snowbank melt water to some extent, but as mentioned by Nabe-Nielsen et al. (2017), it is unlikely that precipitation alone will recreate the continually wet conditions found near snowbanks.

Monitoring

Two FEC attributes identified as essential by the Circumpolar Biodiversity Monitoring Programme (Christensen et al. 2020) were quantified in this study, namely diversity (alpha diversity, rare species, community composition) and abundance (percent cover). Along with voucher specimens representative of plot vegetation species and deposited at the Louis-Marie Herbarium (Quebec City, Quebec, Canada) and at the Canadian Museum of Nature Herbarium (Ottawa, Ontario, Canada) (Desjardins et al. 2021a), we archived in Nordicana D (Desjardins et al. 2021b) the complete floristic database of 264 georeferenced and photographed vegetation plots, including 50 plots with permanently marked quadrats. Communities currently under-represented will have additional vegetation plots marked and added to the database in the coming years. Permanently marked quadrats surveyed with the point-intercept method establish a strong baseline for the continuous monitoring of vegetation or for future comparisons at specific points in time (Elmendorf et al. 2012; Bjorkman et al. 2020). Indeed, whereas re-surveying provides the best source of information on the status and trends of taxa, the ability to detect trends strongly depends on sampling standardization. Our baseline reference is also important given the evident geographical gaps in vegetation monitoring that exist in Siberia and in large parts of the Canadian Arctic (Aronsson et al. 2021), and given the large variation among sites and species in the direction and magnitude of change in plant abundance (Bjorkman et al. 2020).

Based solely on vascular plant species and cryptogams at the highest taxonomic level, we observed distinct plant communities among the surveyed plots. However, we recognize that including cryptogams identified to species to the classification would increase the discriminating power of the analysis (Sheard and Geale 1983; Bliss and Svoboda 1984; Maycock and Fahselt 1992), as cryptogam diversity is known to significantly differ between locations and habitat types (Wietrzyk-Pełka et al. 2021). This is especially true considering that cryptogams are a major component of Arctic plant communities, with vascular plants, bryophytes, and lichens contributing, respectively, 1%, 6%, and 10% of overall Arctic plant diversity (Lewis et al. 2017; Zwolicki et al. 2020). The small size and difficult taxonomic identification of cryptogams contribute to their underrepresentation in ecological studies (Belnap and Lange 2003; Wietrzyk-Pełka et al. 2021), although they are very important for ecosystem functioning (Shively et al. 2001; Garcia-Pichel and Wojciechowski 2009; Belnap and Lange 2003; Pushkareva et al. 2017). To fully understand the impact of climate change on vascular plant and cryptogam communities and on ecosystem functioning, cryptogam species should be included in future vegetation studies at Alert.

Long-term monitoring of vegetation is especially important where anthropogenic pressures (frequent use of heavy vehicles, recreational trampling, infrastructure expansion) may negatively impact the landscape, which could be the case in the vicinity of CFS Alert. We highlight that 18 vascular species, found only opportunistically or in less than 10 vegetation plots, were rare and had a very local distribution. Several rare species were also highly restricted to a specific wetland community, which increases their vulnerability. These rare species and species of concern will be important targets for future monitoring (Christensen et al. 2013).

To accommodate changes in vegetation in response to outside influences, monitoring programmes must include conceptual models on expected responses and their drivers (Lindenmayer and Likens 2010). These models should establish which vegetation parameters should change, what main drivers are likely to be at play, and how variables can be monitored to inform trends and causal relationships (Lindenmayer and Likens 2010). Vegetation monitoring should thus be integrated with monitoring of other ecosystem components and environmental conditions, such as climate. It is our ambition to develop an ecosystem-based monitoring programme at Alert, similar to programmes maintained elsewhere in Canada (Gauthier et al. 2013), as well as in Norway and Greenland (Ims and Ehrich 2013).

Conclusion

This study adds knowledge to the composition and structure of Arctic plant communities in an understudied but highly vulnerable ecosystem. Along with typical polar desert plant communities found on drier grounds, we found that wetland habitats at Alert, despite representing only a small fraction of the total landscape, were occupied by three distinct communities. A detailed knowledge of local plant communities and species is a key prerequisite for ecological research and environmental assessments.

Clear baseline knowledge of species and communities is also essential in providing a template for recovery from disturbances (e.g., human activities and climate change), including identifying native species that could be used for revegetation efforts (Cray and Pollard 2015). Plot-scale floristic surveys combined with high-resolution imagery allow mapping of vegetation at local and regional scales and contribute to delineating essential wildlife foraging and breeding habitats, thereby helping to manage human activities in these areas. In addition, data on the occurrence of species that are rare, endemic, or of conservation concern will be a useful international contribution to the comprehensive ‘Red list of the Arctic vascular plants’, in development through the CAFF Group of the Arctic Council (Aronsson et al. 2021), which evaluates the status of rare and endangered Arctic plants based on the criteria used by the International Union for Conservation of Nature. Finally, the current status and distribution of plant species and communities could also help interpreting past climatic variations and will be crucial for predicting the ecological impact of the ongoing climate change in the Arctic.

Acknowledgments

We thank Marie-Pier Poulin, Gabrielle Roy, Michaël Fortier, Marie-Jeanne Rioux, Kevin Young, Jacob Caron Carrier, Francis Robitaille, Justine Drolet, Roxanne Gauvreau, and Marc Landry for helping collect the data. We thank Esther Lévesque for helping with biological soil crust classification and for commenting on the manuscript. We thank Martin Dubé, Paul Sokoloff, Annie St-Louis, Luc Sirois and Jade Legros for their help with species identification. We thank Alain Caron and Alexis Grenier-Potvin for helping with the satellite picture analysis. We thank Nathan Koutroulides, Station Warrant Officers Patrick Marceau and Dwayne Fox, as well as all the CFS station personnel for the support during the field work.

Disclosure statement

No potential conflict of interest was reported by the authors

Data availability statement

The data that support the findings of this study will be openly available in ‘Nordicana D’ at http://www.cen.ulaval.ca/nordicanad/dpage.aspx?doi=45711CE-0B46454CA9594B53.

Additional information

Funding

This work was financially supported by the Canada Research Chairs Program, Department of National Defence of Canada, Natural Sciences and Engineering Research Council of Canada (NSERC) under Grant numbers [RGPIN-2019-05292 and RGPNS-2019-305531], Fonds de recherche du Québec – Nature et technologies (FRQNT), Network of Centers of Excellence of Canada ArcticNet, Northern Scientific Training Program (Polar Knowledge Canada), BIOS² NSERC CREATE program, and Orville Erickson memorial scholarship fund.

Appendix A.

Definitions of Focal Ecosystem Components (FEC) as modified from the Arctic terrestrial biodiversity monitoring plan of the Arctic Council working group Conservation of Arctic Flora and Fauna (CAFF; Christensen et al. 2013).

Table

FEC	Description
Tree	Single woody stem, tall (greater than snow depth).
Deciduous shrub	Multiple-stemmed woody plant, generally > 60 cm, with leaves not persistent in winter (e.g., Salix Linnaeus).
Evergreen shrub	Multiple-stemmed woody plant with persistent green leaves (e.g., Dryas Linnaeus).
Forb	Broad-leaved herbaceous plant, mostly dicots, non-woody, annual or perennial; also includes horsetails.
Graminoid	Grass-like herbaceous plant with leaves mostly very narrow or linear in outline; includes grasses, rushes and sedges.
Moss	Small, soft plant growing in mats or clumps; includes peat moss (e.g., Sphagnum Linnaeus).
Lichen	Low-growing fungus-like organism, not forming leafy stems like moss, often in continuous mats or crusts.

Appendix B.

Percentage cover of the ground substrates, cryptogams, and vascular species for five plant communities identified at Alert (Nunavut, Canada). Significant (p < 0.05) and highly significant (p < 0.001) indicator species are identified with light grey and dark grey shading, respectively (detailed results of the indicator species analysis appear in Table 2). Frequency of occurrence represents the percentage of vegetation plots in which a given vascular species, cryptogam or ground substrate occurred, out of 264 vegetation plots each comprised of five 1-m² quadrats. *Arctic endemic species according to Daniëls et al. (2020).

Table

Plant community	I	II	III	IV	V
	Cover (%)					Frequency of occurrence (%)
Ground substrates
Soil	30.37	33.28	3.28	0.36	13.16	92.42
Rock	26.80	54.73	0.76	3.36	7.80	95.45
Cryptogams
Black biological soil crust	0.64	0.00	0.12	3.54	3.21	19.70
White biological soil crust	0.64	0.20	0.05	0.52	0.39	56.82
Blue-green algae	0.01	0.00	0.09	0.24	0.51	12.12
Lichen	0.68	1.47	0.00	0.48	0.23	60.98
Macrofungus	0.01	0.00	0.03	0.06	0.04	9.09
Moss	8.70	1.64	20.46	53.00	22.29	91.29
Vascular plants
Asteraceae
Taraxacum phymatocarpum J. Vahl	0.03	0.10	0.00	0.00	0.01	6.06
Brassicaceae
Braya purpurascens (R. Brown) Bunge ex Ledebour	0.14	0.22	0.01	0.12	0.07	24.24
Braya thorild-wulffii Ostenfeld subsp. thorild-wulffii*	0.15	0.30	0.01	0.00	< 0.01	16.67
Cardamine bellidifolia Linnaeus	< 0.01	0.00	< 0.01	2.12	0.14	8.33
Cochlearia groenlandica Linnaeus	0.02	0.06	0.02	0.32	0.14	12.12
Draba corymbosa R. Brown ex de Candolle	1.50	0.36	0.12	1.24	1.00	77.27
Draba lactea Adams	0.01	0.00	0.13	0.14	0.01	3.79
Draba micropetala Hooker*	0.27	0.03	0.03	0.50	0.33	32.20
Draba pauciflora R. Brown*	< 0.01	0.00	0.00	0.00	0.00	0.38
Draba subcapitata Simmons	0.13	0.02	< 0.01	0.14	0.16	14.77
Caryophyllaceae
Cerastium arcticum Lange*	2.55	0.31	0.09	1.70	0.64	70.83
Cerastium regelii Ostenfeld	0.20	0.03	0.12	1.24	1.24	39.39
Sabulina rossii (R. Brown ex Richardson) Dillenberger & Kadereit*	0.11	< 0.01	0.16	0.00	0.26	22.72
Sabulina rubella (Wahlenberg) Dillenberger & Kadereit	1.16	0.20	0.05	0.16	0.56	67.05
Silene uralensis (Ruprecht) Bocquet subsp. uralensis*	0.03	0.03	0.12	0.00	0.05	12.88
Stellaria longipes Goldie subsp. longipes	5.39	0.48	0.16	1.44	1.93	77.27
Cyperaceae
Carex fuliginosa Schkuhr	0.00	0.00	0.04	0.00	< 0.01	2.27
Eriophorum triste (Th. Fries) Hadac & Á. Löve	< 0.01	0.00	38.20	0.00	0.46	10.98
Equisetaceae
Equisetum arvense Linnaeus	0.00	0.00	0.15	0.00	0.00	0.38
Equisetum variegatum Schleicher ex F. Weber & D. Mohr subsp. variegatum	0.00	0.00	0.70	0.00	0.02	2.65
Juncaceae
Juncus biglumis Linnaeus	1.00	0.02	5.14	2.72	7.38	53.41
Luzula nivalis (Laestadius) Sprengel	0.03	0.00	0.45	6.26	0.42	10.23
Orobanchaceae
Pedicularis hirsuta Linnaeus	0.01	0.03	0.21	0.00	0.22	9.47
Papaveraceae
Papaver dahlianum Nordhagen	1.47	0.21	0.01	2.68	0.73	70.83
Poaceae
Alopecurus magellanicus Lamarck	0.44	0.02	4.67	2.02	22.13	39.77
Arctagrostis latifolia (R. Brown) Grisebach subsp. latifolia	< 0.01	0.00	9.30	0.00	0.31	11.74
Deschampsia cespitosa subsp. septentrionalis Chiapella	0.02	0.12	0.17	0.00	0.35	8.33
Festuca baffinensis Polunin	0.07	0.00	0.00	0.00	< 0.01	1.89
Festuca brachyphylla Schultes & Schultes f. subsp. brachyphylla	0.11	0.00	0.00	0.08	0.08	4.92
Festuca edlundiae S.G. Aiken, Consaul, & Lefkovitch*	0.17	0.02	0.00	0.00	0.16	12.50
Festuca hyperborea Holmen ex Frederiksen*	0.01	0.00	0.00	0.12	< 0.01	1.89
Festuca viviparoidea Krajina ex Pavlick subsp. viviparoidea	0.08	< 0.01	0.00	0.00	0.16	3.03
Phippsia algida (Solander) R. Brown	0.61	< 0.01	0.00	3.00	0.21	10.98
Pleuropogon sabinei R. Brown*	0.00	0.00	0.00	0.12	0.36	0.76
Poa abbreviata R. Brown subsp. abbreviata	0.64	0.71	0.05	0.12	0.07	44.70
Poa arctica R. Brown subsp. arctica	0.00	0.00	0.00	1.00	0.00	0.38
Poa pratensis subsp. colpodea (Th. Fries) Tzvelev*	0.01	0.01	0.00	0.00	0.00	1.52
Puccinellia angustata (R. Brown) E.L. Rand & Redfield*	0.26	0.42	< 0.01	0.00	0.24	5.30
Puccinellia bruggemannii T.J. Sørensen*	2.23	0.55	0.07	0.12	0.53	40.91
Puccinellia vahliana (Liebmann) Scribner & Merrill*	0.65	0.20	0.07	0.12	2.34	31.06
×Pucciphippsia vacillans (T. Fries) Tzvelev*	0.17	< 0.01	0.05	0.00	0.57	7.58
Polygonaceae
Bistorta vivipara (Linnaeus) Delarbre	0.21	0.09	4.62	0.00	0.44	17.80
Oxyria digyna (Linnaeus) Hill	0.35	0.11	0.00	0.10	0.03	17.42
Ranunculaceae
Ranunculus hyperboreus Rottbøll	0.00	0.00	0.00	0.00	0.04	0.38
Ranunculus sabinei R. Brown*	0.00	0.00	0.00	0.00	0.03	1.52
Ranunculus sulphureus Solander	0.00	0.00	< 0.01	1.04	0.03	2.27
Rosaceae
Dryas integrifolia Vahl subsp. integrifolia	< 0.01	0.06	0.16	0.00	< 0.01	4.17
Potentilla pulchella R. Brown	0.08	0.10	< 0.01	0.08	0.02	13.26
Salicaceae
Salix arctica Pallas	2.06	1.71	10.13	0.04	4.33	45.83
Saxifragaceae
Micranthes tenuis (Wahlenberg) Small	0.02	0.00	0.02	1.52	0.17	12.12
Saxifraga cernua Linnaeus	0.21	< 0.01	0.05	6.30	0.72	41.67
Saxifraga cespitosa Linnaeus	0.42	< 0.01	0.03	1.66	0.77	40.91
Saxifraga flagellaris subsp. platysepala (Trautvetter) A.E. Porsild	0.01	0.00	0.02	0.16	0.12	13.64
Saxifraga oppositifolia Linnaeus subsp. oppositifolia	10.58	3.00	0.72	1.58	4.03	89.39

Appendix C.

Summary of eight environmental variables for five plant communities identified at Alert (Nunavut, Canada)

Table

Plant community	I	II	III	IV	V
Environmental variables	Mean ± SD					Min-Max
NDVI	0.01 ± 0.05	−0.08 ± 0.04	0.22 ± 0.06	0.28 ± 0.07	0.13 ± 0.06	−0.19–0.38
Elevation (m)	157.79 ± 86.57	104.13 ± 89.90	55.31 ± 40.71	277.60 ± 86.85	128.40 ± 82.87	0–365
Slope (%)	4.50 ± 3.65	4.13 ± 3.90	6.41 ± 4.61	8.59 ± 7.63	4.72 ± 3.49	0–21.31
Distance to a snowbank (m)	695.24 ± 526.48	700.13 ± 633.72	351.35 ± 268.44	494.60 ± 153.27	776.10 ± 731.40	5.87–3101.17
Below a snowbank	0/78 (0%)	3/87 (3%)	12/26 (46%)	1/5 (20%)	8/68 (12%)
Distance to a lake/pond (m)	1622.93 ± 1184.09	1195.72 ± 1075.71	855.04 ± 703.24	2894.47 ± 2098.12	1292.76 ± 992.20	12.66–5235.68
Distance to a river (m)	305.20 ± 296.10	288.52 ± 292.95	315.88 ± 249.15	101.22 ± 130.80	257.13 ± 219.50	0.5–1602.63
Distance to an ocean shore (m)	3203.83 ± 2261.34	2179.78 ± 2333.66	1219.18 ± 1006.12	6890.71 ± 1139.00	2734.32 ± 2227.54	58.87–8226.52

Appendix D.

Relations between eight environmental variables and cover of vascular plant species, cryptogams, and ground substrates, as assessed from a Canonical Correspondence Analysis. *** p < 0.001, * p < 0.05

Table

Environmental variables	CCA1	CCA2	Chi Square	F statistic	p
NDVI	0.99	0.15	0.55	52.32	0.001 ***
Elevation	−0.18	0.79	0.10	9.17	0.001 ***
Distance to nearest snowbank	−0.08	0.21	0.04	3.80	0.001 ***
Location with respect to a snowbank	0.47	−0.43	0.05	4.32	0.001 ***
Distance to nearest ocean shore	−0.13	0.68	0.02	1.71	0.050 *
Slope	0.16	−0.14	0.01	1.36	0.172
Distance to nearest lake or pond	−0.14	0.48	0.01	1.23	0.228
Distance to nearest river	−0.01	−0.14	0.01	1.15	0.268