Location of studies and evidence of effects of herbivory on Arctic vegetation: a systematic map

Background: Herbivores modify the structure and function of tundra ecosystems. Understanding their impacts is necessary to assess the responses of these ecosystems to ongoing environmental changes. However, the effects of herbivores on plants and ecosystem structure and function vary across the Arctic. Strong spatial variation in herbivore effects implies that the results of individual studies on herbivory depend on local conditions, i


Background
Herbivores are a central component of northern socioecological systems, both as subsistence resources and due to their influence on ecosystem structure and function [1,2]. For example, many northern communities rely on vertebrate herbivores to support local economies through hunting (e.g., geese, ptarmigan, caribou) or herding (e.g., reindeer) [3,4]. Similarly, outbreaks of invertebrate herbivores can influence the livelihood of northern communities [5] through their dramatic impacts on vegetation and the subsequent cascading effects on local wildlife [6,7]. Herbivores also play an important ecological role in Arctic tundra ecosystems by influencing soils [8] and plant communities [9,10], with consequences for ecosystem processes [11]. Moreover, herbivores have the potential to offset some of the effects of ongoing warming in tundra ecosystems by counteracting increases in deciduous shrub cover [12,13] or advances of the treeline [14]. Thus, understanding how herbivory influences the responses of Arctic ecosystems to environmental change can help guide appropriate adaptive strategies to sustain ecosystem functioning and the provision of ecosystem services [1].
However, the effects of herbivores on plants and ecosystem functioning are far from homogeneous across the Arctic. Rather, herbivore effects depend on the local conditions, i.e., ecological context [15,16]. We here define ecological context as the range of physical and biological conditions that can affect the interplay between plants and herbivores, and include descriptions of geography, climate, food web composition and structure and human use. For example, plant responses to herbivory have been shown to depend on local productivity [10,17] and temperature [12]. The effects of herbivores on plants are thus expected to differ along climatic gradients within the Arctic [18]. Distance to the coast can matter due to marine subsidies to the food web [19] or salinity [20,21] and distance to the treeline due to spillover of boreal herbivores [22][23][24]. The magnitude of recent warming and the composition of herbivore communities are also important [25,26]. The extent to which the literature on Arctic herbivory covers the range of possible ecological contexts will determine which conclusions can be drawn about the role of herbivores in shaping Arctic vegetation. While uneven geographical cover of the Arctic has been identified for environmental studies more generally [27], we lack a systematic overview of which ecological contexts are covered by the existing literature on herbivore impacts.
To better understand the present "state of affairs" of Arctic herbivory research, we systematically mapped the coverage of ecological context in studies considering the effects of herbivores on vegetation. Individual studies that address the effects of herbivores on tundra ecosystems are often conducted at a specific location and within a specific ecological context. Although these studies are designed in ways that are adapted to their local context, they can also be considered as sample points within the larger, discipline-wide study of the effects of herbivores within the circumpolar Arctic. A first attempt to assess whether generalizable conclusions can be produced from these studies includes gathering information about the number of studies and assessing how well they cover the underlying ecological gradients across the Arctic. Uneven distribution of studies across a particular ecological gradient indicates that enough evidence may be available for some parts of the gradient but not for others. In such case conclusions on the effects of herbivores cannot be generalized across the Arctic. This systematic map will thus assess our ability to make generalizable and robust conclusions on Arctic herbivory. Furthermore, identifying areas that are underrepresented in the current evidence base will facilitate targeted research efforts to fill these gaps.
The topic for this systematic map was originally identified at a specialist workshop aimed at identifying data gaps in herbivory research, where the main stakeholder group of scientists working on Arctic herbivory was invited to participate. A core team developed the protocol that is used in this systematic map [28]. Active participation in the systematic map and co-authorship was encouraged through an open call for collaboration through the Herbivory Network (https:// herbi vory. lbhi. is/), which currently features more than 200 members and includes relevant stakeholders for the topic of this systematic map.

Objective of the review
We used a published protocol to map the coverage of ecological contexts in which herbivory has been studied in the Arctic [28]. Hence, the primary question of the systematic map was "What evidence exists on the effects of herbivores on Arctic vegetation?" and it covers terrestrial Arctic plants and plant communities (population) responses (outcome) to herbivory (exposure). Note that we did not assess the strength of the herbivore impacts on vegetation but focused on the coverage of the existing studies across the underlying environmental conditions across the Arctic. In addition to vertebrate and invertebrate herbivores' direct effects on plants through consumption, exposure also included their disturbance and fertilization effects. Studies also had to compare a certain level of herbivory to no herbivory or alternative levels of herbivory (comparator), such as experimental exclusion or herbivore abundance changes across space or time. Thus, the primary question of the systematic map contains four key elements, needed for the question to be answerable: the population of interest (terrestrial Arctic plants and plant communities), the exposure (herbivory); the comparator (different levels of herbivory) and the outcome (any response from the population to the exposure). This PECO-type question defines the eligibility criteria for individual studies, as outlined below.

Methods
This systematic map follows methods described in a protocol by Soininen et al. [28] to map the existing evidence on the effects of herbivory on Arctic vegetation. Throughout this map, we use the following terminology: an article refers to any publication/document matching our search criteria (below), and a study is a report of herbivory effects using a certain study design in a given ecological context. An article could thus contain from one to several studies. The systematic map conforms to ROSES reporting standards (see Additional file 1).

Deviations from the protocol
Some minor deviations from the original protocol were necessary when developing the systematic map and are described in the sections below. These deviations included (i) the use of other maps for delimiting the study region as the original one was not available and (ii) limiting the check of redundant studies to those that were conducted in the same location. The latter was deemed sufficient as only studies from a shared location could be redundant. In addition, (iii) assessing repeatability of the data coding strategy at the full-text stage by replicate extraction of data from at least 10% of the studies was deemed unnecessary as all coded data were checked for consistency by two reviewers.

Search for articles
We searched for relevant literature first from the global databases of Web of Science Core Collection (all years search within Topic) and Scopus (article title, abstract and keyword search with no further limitations applied). Thereafter, we searched literature that these databases might have missed, and grey literature, from Google Scholar (title search standardized to disregard search history, including only the first 300 search results) and local/ regional and specialist databases or sources (listed in Additional file 2).
The search string for the systematic map was optimized during scoping exercises as described by Soininen et al. [28] and comprised two substrings: one targeted at delimiting the study region and system, and the other targeted at the exposure component of our primary question (i.e., herbivory). The full search string (formatted for Web of Science) was: (arctic OR subarctic OR tundra) AND (herbivor* OR graz* OR browser OR browsing OR grubb* OR trampl* OR defolia* OR ((invertebrate OR insect) AND (gall* OR mining OR miner))) The specificity of the search string, that is the proportion of relevant studies returned by the search, was relatively low, with 46% of the documents excluded at title screening stage in the current search (Fig. 1). Still, the specificity of the search string was deemed adequate, as it was not possible to include additional search terms to narrow down the population of interest [28]. The sensitivity of the search string, that is its ability to find all relevant studies was assessed by comparing the results obtained with the full search and the two substrings separately, to identify potentially relevant documents that could have been missed. The sensitivity of the search string was very high, with only one article out of a subset of 1000 articles deemed as potentially relevant [28].
We conducted our search in English in global search sources, and English together with relevant local languages (Russian, French, Finnish, Swedish, Norwegian, Icelandic and Danish) in searches from local/regional sources. In some cases, the full search strings had to be simplified because some search engines did not allow Boolean operators or the use of "wildcards" (*) or could handle only a reduced number of terms. Full details on the search sources and strings used in each search engine and how they were implemented in other languages, as well as access information on institutional access to paywalled databases are provided in Additional file 2.

Screening process
Screening of the articles for relevancy was structured in four stages: (i) title screening, (ii) abstract screening, (iii) locality screening, and (iv) screening during data coding from the full text. Each article was assessed by two reviewers at the title and abstract screening stages, except for those found through Google Scholar and local/regional databases that were assessed by one reviewer. Locality and full text screening were assessed by one reviewer, with additional reviewers in cases that were not clear. These numbers refer to reviewers per RepOrting standards for Systematic Evidence Syntheses (ROSES) diagram indicating the inclusion and exclusion process, including the number of studies retained at each stage of the process. At the full text stage, studies were excluded if they were located outside the study region (locality), if they did not target Arctic terrestrial vegetation (population), did not assess effects of herbivores on plants (exposure), did not compare different levels of herbivory (comparator), focused on ecosystem components other than plants (outcome), their design did not address current ecological contexts in the Arctic (study design), represented unsuitable text types or presented the same data as other studies (redundant). The ROSES diagram is based on [97] article; several reviewers were involved in each screening phase.
During screening, we erred on the side of inclusion. Exclusion of an article at title and abstract screening stages occurred upon the agreement of two independent reviewers. The reviewers first noted whether an article should be excluded based on the title. If not, they screened the abstract and noted either "excluded based on abstract" or "included". In the case of reviewer disagreement about inclusion, the article was kept for further stages of screening. No reviewers assessed articles they had authored or co-authored. When one reviewer had scored an article to be excluded based on title and the other reviewer based on abstract, the article was assigned the score "exclude based on abstract". Additional articles that were found through Google Scholar or from local/ regional sources were screened for both title and abstract by one reviewer and were only included in our database if they passed this screening stage.
Articles that passed the title and abstract screening stage were then assessed for their geographical locality by one reviewer. Locality screening was based on geographical coordinates that were given in the text, extracted from maps provided in the publication, or extracted from Google Maps based on study location names mentioned in the text. To delimit the study region, we used the southern limit of the Subarctic as defined by the Conservation of Arctic Flora and Fauna (CAFF, [29]). This region contains tundra, the tundra-forest ecotone and boreal forests. To limit the scope to tundra and tundra-forest ecotone, we excluded areas south of the Subarctic zone, based on the definition used in the Arctic Biodiversity Assessment [30] and excluded studies in boreal forests based on full text screening.
When extracting locality coordinates, the reviewer also did a first screening of the full text to identify any obvious reasons for exclusion, such as being a review paper with no primary data. Articles that passed the locality screening were thereafter read in full and their content judged by at least one reviewer. The full text screening stage was in practice combined with the data coding. Reviewers carefully checked the full text against the eligibility criteria and reasons for exclusion were recorded. When in doubt, the article was reviewed by an additional reviewer. The complete database of retrieved articles, with information on exclusion stage and reasons (for full text stage only) is provided in Additional file 3 (see also Fig. 1). Locality screening diverges from the description in the published protocol [28], where the CAFF vegetation zone map ("CAFF Map Nr. 10", [31]) was planned to be used to delimit the Arctic region. That map was however no longer available, leading us to use the approach described here.
We assessed the repeatability of title and abstract screening by measuring the consistency among reviewers to either include or exclude an article. For this, we used percent agreement and Cohen's kappa (κ) statistic [32]. Percent agreement is calculated as the number of agreements divided by the total number of screened articles, and kappa is a ratio of the relative observed agreement to the probability of random allocation. Consistency among two reviewers to exclude an article based on title screening alone was low (49.3%, κ = 0.54 indicating 'weak agreement' according to McHugh [33], N = 2704). After abstract screening, consistency of scoring was considerably higher: two independent reviewers agreed to exclude an article in 74.1% of cases, include an article in 13.4% of cases and disagreed only in 12.5% of cases (κ = 0.61, 'moderate agreement' [33], N = 2704). The majority (79.6%) of articles where reviewers disagreed after the title and abstract screening stage were excluded based on their locality or full-text evaluation. One of the main reasons for exclusion at the full-text stage was the lack of an eligible outcome.
Assessing the repeatability of the location screening stage was deemed not necessary, as coordinates were extracted from the study (either from the text, figures, maps or inferred from site names) and checked against the study area using maps. Full text screening was combined with data coding and conducted by a single reviewer, except in unclear cases where the opinion of additional reviewers was sought. Reasons for exclusion at full text screening are listed in Additional file 3. The consistency of this stage was ensured by checking all coded data by two reviewers (see below). In addition, 45 articles excluded at full text stage (corresponding to 10% of articles excluded at this stage for other reasons than location) were independently assessed by a second reviewer. Consistency among the two reviewers to exclude an article based on full text was high; all but one of the articles were unanimously excluded. For the last article, the second reviewer sought for the opinion of an additional reviewer, whereafter also this article was deemed to be excluded.

Eligibility criteria
Eligible studies were defined as those as reporting primary data on herbivore impacts on Arctic plants and plant communities. To be included, studies had to include an eligible population, exposure, comparator and outcome (as defined in the objective of the review above). In addition, studies had to fulfill some requirements regarding their study design as described in this section.
Eligible population Studies had to target an eligible population, i.e., studies were excluded if they did not study Arctic terrestrial vegetation. This filtering included additional locality information, such as describing the study sites as boreal forests or temperate grasslands and, in some cases, inspecting Google Earth images of the study location to assess the extent of canopy cover and proximity to open tundra. Studies were also excluded on the grounds of not targeting an eligible population if the ecological context did not correspond to any current climatic context of the Arctic, such as palaeoecological studies or greenhouse studies not mimicking Arctic growing conditions.
Eligible exposure Studies had to target an eligible exposure to herbivores, that is, they had to assess the effects of herbivores on plants (consumptive effects, fertilizing, trampling and other types of disturbance). To be included, studies mimicking effects that could be related to herbivore activities needed to explicitly state their link to herbivory.
Eligible comparator Furthermore, studies had to assess the effect of herbivores by comparing a given level of herbivory or other herbivore-related effects to either no herbivory or another level(s) of herbivory; that is, studies needed to have an eligible comparator. The type of comparison was not restricted, and thus different types (e.g., factor levels of experimental treatments, continuous variable changing across spatial or temporal gradients) and intensities of the comparison (e.g., number of herbivory levels, magnitude of the difference in herbivory along a gradient) were eligible. We included studies that measured herbivore presence in direct and indirect ways such as biting marks or galls but excluded those that did not relate herbivore presence to any measurable effect of herbivores on plants.
Eligible outcome We excluded studies lacking an eligible outcome if they assessed the effects of herbivory on ecosystem components other than plants or plant communities. We included lichens but excluded soil properties and plant symbionts such as microbial endophytes or mycorrhizal fungi. Studies were also excluded on the grounds of lacking an eligible outcome when they presented no quantitatively analyzable primary data, as in the case of many reviews, book chapters, synthesis papers, or perspectives.
Eligible types of study design Studies had to address herbivore effects on plants in a given ecological context within the Arctic; thus, some modeling studies were excluded based on their lack of an eligible study design.
Additional criteria We excluded studies for which we could not access electronic or paper copies of the full text, as this prevented assigning studies to an ecological context. Similarly, we excluded studies that were not in a suitable format (Fig. 1), such as corrections of published articles, maps and raw datasets (cf. "Unsuitable text types" in ref. [34]). We further excluded studies that presented data already reported in another article. For chapters of academic theses that were published as separate studies, we only included the published peer-reviewed versions, and for annual reports that presented incremental data, we only included the longest time-series. For all other studies, we checked for potential redundancy based on the study site. We first edited the spelling of study sites to be consistent across studies and then cross-checked all studies conducted at a given study site. We coded a study as redundant if it presented a temporal or spatial subset of data that was also included in another article. In cases of data duplication, we included the most recent study.

Study validity assessment
Because the primary aim of the systematic map was to identify the ecological contexts where herbivory has been studied, critical appraisal of the validity of the studies based on aspects such as replication, confounding factors, or issues with interpretation of measured variables is not of paramount importance (but see examples of systematic reviews on herbivores where such appraisal has been done [18,35]). Thus, we did not exclude any study based on such criteria but included an appraisal of (i) the spatio-temporal extent and resolution of the study designs, and (ii) lack of information on study design. This information was extracted from the studies as part of the data coding (Table 1), and allowed us to explore whether local (i.e., study area size less than 10 × 10 km) or shortterm (i.e., 1-year) studies were more common for some of the ecological contexts or herbivore groups than others. The definitions of these spatial and temporal groupings were, as described in the protocol [28], based on exploring the ranges of these variables in the studies retrieved by our search.

Data coding strategy
For each study that passed all the screening stages, we extracted information on variables reporting the study design and study method, the target population (plant), exposure (type of herbivore) and comparator (how the contrasting levels of herbivory were obtained) and variables describing the ecological context (Table 1). Additional information for these variables, such as potential values, examples and specifications, and other variables extracted from the studies but not presented in detail in this paper are given in Additional file 4. Table 1 Ecological context variables extracted from the studies included in the systematic map database during data coding Values in column "source" describe where the data is extracted from: P for publication (i.e., data explicitly available in the published document), D for digital spatial data layers (i.e., data extracted from a data layer based on the spatial coordinates given in the publication), and C for classified by the reviewers based on information available in the publication. For variables denoted with P and D, we used the D-version in our analysis to reduce the heterogeneity in the accuracy reported across studies. See Additional file 4 for a full list of variables extracted from the studies a Variable added after publication of the protocol b Variable data source differs from the planned source (i.e., from the reference given in protocol) c Note that all 19 WorldClim variables were extracted and explored, but here we list only the variables for which results are presented in the text. The variables differ from the protocol, where use of three main PCA axis of variation across these variables was outlined The values for the ecological context variables were recorded as they were reported in the study, either with data explicitly available in the publication (variables indicated with P in Table 1), classified by the reviewers based on information available in the publication (C in Table 1) or extracted from a data layer based on the spatial coordinates given in the publication (D in Table 1). We did not attempt to contact the authors of the studies to obtain missing information or confirm unclear points, as this could have biased the quality of information extracted from studies where corresponding authors were no longer available.
One article could contain information about several studies when it included more than one study question, methodological approach or ecological context. For instance, when an article presented different datasets of plant-herbivore pairs, such as one experimental and one observational dataset, these were considered as separate studies. Similarly, when an article presented data from several locations for which ecological contexts were described separately, datasets from the different locations were considered as separate studies. In cases where an article contained more than one study, information about the variables was recorded for each study.
After all the data were coded, two reviewers (EMS, ICB) edited possible dissimilarities in spelling, checked accuracy and corrected obvious errors (such as deviations from pre-defined factor levels), and identified additional uncertainties in the inclusion of studies. At this stage, lengthy descriptions for some variable values were summarized into more general categories. For example, when detailed descriptions of the habitats were provided in the text, they were replaced by a shorter, more general description. In the original protocol [28] it was stated that data from at least 10% of the studies would be extracted by at least two reviewers to reduce discrepancies and inconsistencies between reviewers. Since all coding data were checked for consistency by two reviewers, this step was considered unnecessary.

Data mapping method
To code data from the full text of each study, we used the Excel spreadsheet published as an appendix of the systematic map protocol [28]. The data were collated using R-scripts (available on GitHub https:// github. com/ James DMSpe ed/ Arcti cHerb ivory Syste matic Map). We then extracted geospatial ecological context variables (Table 1) and added them to the dataset.
We grouped the ecological context variables according to whether they described the geography, climate, climate change, food web, or human use (see Table 1 for variables, their definitions, and how they were extracted).
We used univariate and pairwise plots to visually explore the available ecological space within the Arctic (i.e., the range of values of the variables that are represented in the study region) and to assess how much of that space is covered by the available studies. We delimited the Arctic study region by coastline but did not remove other areas such as glaciers or lakes. In our narrative synthesis of the results, we highlight biases in the available studies across the ecological contexts, as they point to areas for which abundant information exists and where more studies are needed.
The raw data are available for download (Additional file 5), as well as through a Shiny App: https:// shiny. vm. ntnu. no/ users/ speed/ Arcti cHerb ivory Syste matic Map/, created with the R-package "Shiny" [36]. The Shiny App helps the reader navigate the database and visualize the contribution of individual studies. The app allows plotting studies on a geographic map, filtering studies based on different criteria (e.g., selecting country, species, study language, herbivore type, and different aspects of the study design) and interactive exploration of the variables using bar plots, histograms and pairwise plots. This allows the reader to explore the data freely, including variables that are not the focus of the narrative synthesis of results in this study. The pairwise plots allow simultaneous visualization of study design variables together with two ecological context variables to identify which contexts have been studied using different approaches.

Review descriptive statistics
The literature search underpinning this review was conducted on 13th February 2019 (Web of Science and Scopus, refer to Additional file 2 for dates in local searches). We identified 3200 articles from searches in Web of Science and Scopus, an additional 311 from local/regional and specialist databases/sources as well as an additional 75 from Google Scholar ( Fig. 1; Additional file 3). After evaluation of eligibility across all screening stages, the final evidence base included 309 articles, split into 662 studies. The full list of included articles can be found in Additional file 3, and the full list of studies is presented in the Shiny App.

Mapping the quantity of studies relevant to the questions Approaches, study designs, and methods
Experimental and observational field studies were equally common (n = 299 and 306 studies, respectively, for definitions see Additional file 4). In addition to field studies, the evidence base included remote sensing studies (n = 28), modeling studies (n = 19) and greenhouse studies (n = 12). Most of the experimental field studies had no additional exposures other than herbivory, but among those studies that had additional experimental treatments, warming (n = 26 studies) and nutrient manipulation (n = 23 studies) were the most common.
Most studies were published after 1980, with a stark rise in studies occurring around 2000 and peaking around 2010 ( Fig. 2A, B). The majority of studies reported on data collected only once (n = 363) and/or during only one year (n = 262), although data sampling of 112 studies spanned more than a decade (Fig. 2C-D). The results were most often reported at finegrained spatial resolution (n = 253 up to 1 × 1 m plots and n = 123 from 1 × 1 m to 10 × 10 m plots). The study area size was most often small (n = 170 at 1 × 1 km or less, n = 121 between 1 × 1 km and 10 × 10 km, and n = 65 between 10 × 10 km and 100 × 100 km). Note that we defined the study area size as a square based on the study area description and coordinates, but many studies provided so little information that we were not able to assign them to a study area category (n = 156).

Arctic plants and plant communities (population)
To gain an overview of the studied plant types, we grouped the plant species and groups reported by the studies into plant functional groups (for definitions, see Additional file 6). Graminoids were the most commonly studied plants, with 343 studies including data on this functional group. They were followed by deciduous dwarf shrubs (n = 283), forbs (n = 235), evergreen dwarf shrubs (n = 224), bryophytes (n = 214), lichens (n = 174), deciduous trees (n = 114), the vascular plant community as a whole (n = 94), deciduous tall shrubs (n = 88), ferns and allies (n = 61), evergreen trees (n = 9) and evergreen tall shrubs (n = 4) (Fig. 3A). Note that deciduous and evergreen trees are included here, as the study region also included the tundra-forest ecotone. Fig. 2 The distribution of studies across time and temporal study design. All y-axes show a count of studies. A shows study publication year and B the first year when both plant and herbivore elements were measured. C shows the temporal extent of the studies, defined as the number of years when both plant and herbivore elements were measured. Note that the studies with temporal extent longer than 1000 years come from one modeling study. D shows the temporal resolution of the studies, defined as the interval between measurements of the outcome (plant), if regularly measured or not and for how long. Sample size refers to studies for which the variable was feasible to assign Often, plant response to herbivory was measured at one level of biological organization but results were reported at a more aggregated level. For example, some studies conducted measurements on individual plants but reported their results for the plant population/species (for instance recording survival of marked plant individuals and reporting survival rate per species per treatment). Thus, measurements most often targeted plant data at the population/species level (n = 320 studies), followed by individual (n = 198), plant functional or taxonomic group (n = 93) and plant community (n = 84) levels. In contrast, results were reported at the population/species (n = 315), plant functional or taxonomic group (n = 192), plant community (n = 178) and individual (n = 63) levels (Fig. 3B).

Herbivores (exposure)
The vast majority of studies focused on vertebrate herbivores (n = 582) compared to invertebrates (n = 75) (Fig. 3). We divided vertebrate herbivore species into functional groups according to Speed et al. 2019 [37]. This classification is based on a suite of functional traits reflecting the ecology of the herbivore species and recognizes three functional groups: (i) waterfowl (limnic-habitat associated species migrating outside the Arctic for winter exemplified by snow goose Anser caerulescens), (ii) small rodents and pikas (hindgut-fermenter, burrowing species exemplified by Norwegian lemming Lemmus lemmus), and (iii) other vertebrates (larger-bodied, facultative-generalist species exemplified by reindeer Rangifer tarandus). The functional group "other vertebrates" was most studied (385 studies), and studies of this group were by far dominated by Rangifer tarandus, which alone was the focus of 333 studies (Fig. 3). The group "small rodents and pikas" was studied in 153 studies. These studies were dominated by genera Lemmus and Myodes (81 and 73 studies, respectively). The group "waterfowl" was studied in 128 studies, of which almost half targeted the genus Branta (n = 60) (Fig. 3).
Among invertebrate herbivores, we separated functional groups based on feeding mode: defoliating invertebrates were addressed in 67 studies (mainly Operophtera brumata, n = 32, and Epirrita autumnata, n = 18), while other types of invertebrate herbivores were very rarely studied (n = 3 studies for galling invertebrates, one for phloem feeders, 2 for root-feeding invertebrates and 2 for invertebrates feeding on reproductive structures).
While studies targeting the effects of invertebrate herbivores were mostly observational, there was a balanced mix of experimental and observational approaches for studies on vertebrate herbivores (Fig. 3C). Studies targeting multiple herbivores (n = 90) used experimental approaches more often than observational ones (n = 71 experimental and 19 observational).
Most studies assessing the effects of invertebrate herbivores reported their outcomes at the plant individual or population/species levels, while those assessing the effects of vertebrates reported their outcomes mainly at higher biological organization levels, including plant population/species, groups of species and communities (Fig. 3B).

Comparison between levels of herbivore impact (comparator)
Field studies included a variety of approaches to assess the impacts of herbivory. The most frequent approach in experimental studies was the use of exclosures that eliminated a certain type of herbivory (n = 169 studies), followed by simulated herbivory (n = 83) and to a lesser extent the use of fences with different herbivory levels on either side (n = 34) or enclosures (n = 22). In observational field studies, spatial contrasts in herbivory (n = 176 studies including spatial gradients), fences (n = 46), and herbivore outbreaks (n = 41) were the most common approaches. Temporal contrast in herbivory was also used, but only 18 studies had increasing or decreasing herbivore populations as their main means to compare herbivore intensities. Control-impact or beforeafter-control-impact designs were used in 420 studies, representing 99% of field experiments and 36% of observational field studies.

Response types to herbivory (outcome)
Almost half of the studies (n = 317) measured more than one type of outcome following exposure to herbivory. Overall, plant biomass was the most common response type (n = 247), followed by morphological measures (n = 188, including for example canopy height or leaf demography), vegetation cover (n = 157, including the cover of plant species or groups), measures of plant community diversity (n = 133) and physiological responses (n = 115, including measures of plant chemistry).

Ecological contexts covered by the evidence base
Most studies were clustered in northern Fennoscandia (Fig. 4). Together, Norway (excluding Svalbard), Sweden and Finland represented 40% of all studies (263 out of 662). Given that this region is mainly Subarctic, this contributed to the Subarctic bioclimatic zone being represented by 318 studies, as compared to the 291 studies from actual Arctic bioclimatic zones A-E (see Table 1 for definitions and references for ecological context variables). Other large clusters of studies were found in Svalbard and coastal Alaska (Fig. 4A).  Fig. 3 Herbivore groups, plant groups, and study types. A combines the studies per herbivore group with the plant groups these studies investigated. B shows the count of studies across herbivore group and the level at which results on plants are reported. C shows the count of studies according to herbivore type and study approach. Herbivore group "defoliating inv. " refers to defoliating invertebrate herbivores and "other inv. " to other types of invertebrates. The vertebrate herbivores were first divided into three functional groups according to Speed et al. 2019 [37]; the group "small rodents and pikas" is denoted as "small rodents" in B and C. The functional group of larger vertebrates that overwinter above snow was further sub-divided to genera. Sample size refers to unique combinations; for instance, if a study that included several herbivores it is plotted in each herbivore category Only 9% (n = 59) of studies came from areas that were further than 100 km inland from the coast (Fig. 5A); thus, herbivore impacts are less well covered by the evidence base in inland Arctic areas. The elevational range of the Arctic was covered relatively well (Fig. 5A). Over onethird (39%) of the studies were located below 100 m a.s.l. while the uppermost ones were as high as 1870 m a.s.l., which reflects well the distribution of elevation across the Arctic. Subarctic studies were conducted at higher elevations, on average, than Arctic studies, as the median elevation for Subarctic studies was 330 m a.s.l., and 54 m a.s.l. for studies in the Arctic bioclimatic zones A-E. We used the limit between Subarctic and Arctic to delineate the treeline ( Table 1). Frequency of studies was highest in the vicinity of the treeline, illustrated by a median distance of 163 km.  [92]) and the Subarctic (CAFF/ABA [30]). The thick grey line shows the Arctic boundary defined by CAFF [29]. The study region was delimited within this boundary and areas outside Subarctic areas filtered removed; areas outside study region are shown as white in both panels Sites with continuous permafrost (n = 295) were most common in the evidence base, with the remaining sites being divided into sporadic (n = 124), discontinuous (n = 98), and isolated patches (n = 12) of permafrost. However, the permafrost data layer (Table 1, [38]) yielded"NA" values for sites south of the permafrostcovered Arctic (n = 38 for Subarctic) and for sites close to the coast in many places of the High Arctic (n = 38 for bioclimatic zones A,B, and C combined). The most common soil types in the evidence base were lithosols  ) show data for studies, while blue points (and histograms) describe the range of all ecological contexts across the Arctic study region at 100 × 100 km resolution. Sample size refers to the number of studies with data (i.e., non NA-values). When sample size is given next to the figure subtitle it refers to both variables. When it is given along x and y axis titles, it refers to the specific variable. Note that in these cases the histograms are showing all available data for each of the variables, even though the points are not. For example, histogram for "change in temperature" includes all 662 data points even though only 455 can be plotted against change in NDVI, due to missing data for NDVI (n = 158, distributed rather evenly between Alaska, Canada, Finland, Norway, Russia and Sweden), podzols (n = 148, mainly in Norway and Finland), gleysols (n = 99, mainly in Alaska and Russia) and regosols (n = 63, only in Canada and Greenland). However, the soil data layer (Table 1, [39]) also resulted in all of the 74 studies on Svalbard having"NA" values for soil data. Thus, the permafrost data and soil type data are not entirely informative at the circumpolar scale.
Studies were overrepresented in warm and moderately wet areas (Fig. 5B, for definitions of the climate variables, see Table 1). For instance, half of the studies were in areas with a mean annual temperature above − 3.7 °C, while the mean annual temperature across the entire study region's land surface was − 11.4 °C. Mean annual precipitation has a bimodal distribution (Fig. 5B), with the upper precipitation mode having a higher proportion of studies than the lower mode, which in turn is more characteristic of the study region. Furthermore, studies were more common in parts of the Arctic with less intra-annual temperature variation. For example, median annual temperature range for the evidence base was 34 °C, but 47 °C for the study region. Finally, although the studies did cover the range of growing season lengths present in the Arctic, they were scarce in areas with very long or short growing seasons (Fig. 5C). Areas with low vegetation greenness (based on Normalized Difference Vegetation Index, NDVI, Table 1) were also less well covered (Fig. 5C).
The magnitude of recent climate change was represented by the change in average temperature, vegetation greenness and growing season length over the last decades (Table 1). While the studies overrepresented areas of moderate temperature increase, they covered the range of the observed warming (Figs. 4B, 5D). Furthermore, there was good coverage of the range of changes in vegetation greenness and growing season length change (for patterns of the latter, see the Shiny App).
The evidence base covered the range of Arctic vertebrate herbivore diversity in terms of species richness, phylogenetic diversity and functional diversity (Fig. 5E, for the figure of phylogenetic diversity, see Shiny App). Most studies included no description of the food web of the study site beyond vegetation; 490 studies did not mention non-target herbivores, and 631 did not mention predators.
The ecological contexts created by humans were represented by human population density and the human footprint index, which combines human population density with infrastructure and other land uses (Table 1, [40][41][42]). For these variables, the studies covered the range that is present across the Arctic but were concentrated in areas with a higher human footprint and population density (Fig. 5F). In addition, 148 studies reported current and/or historical management of the area, while 22 were reportedly done in protected areas and 35 mentioned a conservation focus of the work.
When looking at the ecological contexts covered by studies targeting different herbivores, we found important geographical biases. For example, in the case of studies on R. tarandus, which accounted for approximately half of the entire evidence base of herbivores (n = 349), 75% of these data came from just three Arctic regions-Fennoscandia, Svalbard and the Yamal Peninsula (combined n = 264) and are thereby biased towards warm and wet areas close to the forest and coast. Similar geographical biases were also found for other groups of herbivores, like waterfowl: among the 128 studies on waterfowl, only six were on the Eurasian continent. These studies were highly clustered at a few intensively studied locations such as Svalbard, La Pérouse Bay, Bylot Island, Akimiski Island, and the Yukon-Kuskokwim Delta. Studies for waterfowl came mainly from coastal lowlands (119 studies at < 150 m a.s.l., 123 studies within 50 km from the coast), from the Arctic bioclimatic zones (A = 1, B = 4, C = 51, D = 5, E = 30, Subarctic = 21), and in permafrost covered sites (n = 96 studies with continuous permafrost). Studies addressing the functional group of "small rodents and pikas" were most common on the Eurasian continent (n = 64 out of the 79 for this group, note that no pikas inhabit this continent). These studies were mainly close to the coast (65 studies within 50 km from the coast). The 385 studies on vertebrate herbivores other than waterfowl and small rodents (i.e., the functional group "other vertebrates") were spread throughout the Arctic, with only one evidence point further east in Siberia than the Yamal Peninsula (work in the Lena Delta on marmots). These studies covered mainly Subarctic and Low Arctic areas (n = 272 studies in bioclimatic zone E and Subarctic) and areas close to the coast (229 studies within 50 km from the coast). The 67 studies on defoliating invertebrates were predominantly one-time studies (n = 49) relating to herbivore outbreaks (n = 38), at lowland sites with trees (n = 46 below 300 m a.s.l.), and clustered in the warm (n = 55 with higher than − 5 °C annual mean temperature) parts of the Arctic with high precipitation (n = 60 with more than 400 mm annual precipitation), mainly in northern Fennoscandia (n = 47) and eastern Canada (n = 13).

Mapping the validity of studies relevant to the question
Mapping the validity of the studies was not the aim of our systematic map, but we recorded studies with incomplete information regarding their study design. Several studies omitted information describing data sampling to the extent that we were unable to assign the first year of study (n = 45), temporal resolution (i.e., how often the data were sampled, n = 27), how large the study area was (n = 156), or the scale of sampling plots/units (n = 63 and 46 for reporting scale and recording scale). We explored visually and with linear regressions whether the proportion of studies with missing information for these variables changed over time, but found support for a decreasing trend only for the extent of spatial scale (effect of year − 0.01 [− 0.02, − 0.004 95%CI]).
We also explored visually whether local (i.e., study area size less than 10 × 10 km) or short-term (i.e., one-year) studies were more common for some of the ecological contexts (Additional file 7: Figures S1 and S2). This comparison revealed no striking patterns. However, Subarctic regions had an evidence base with a balanced mix of small-and large-scale studies, whereas small-scale studies dominated in the Arctic zones C and E. This difference was further reflected in some of the other variables (e.g. distance to treeline, vegetation greenness and herbivore species richness), presumably because these variables correlate with the bioclimatic zonation.

Limitations of the map
Accessibility of the grey literature is likely to vary greatly among the Arctic countries. For example, MSc-theses from Russia are usually not available online, whereas MSc-theses from the Nordic countries are routinely deposited in institutional and national online portals (e.g., Norway: bora.uib.no, nora.openaccess.no; Iceland: skemman.is). We assessed this source of bias by tabulating the number of MSc and PhD theses per country, and the number of other grey literature studies per country. A total of 39 student theses were retrieved in our search, 4 BSc theses, 15 MSc theses and 20 PhD theses. Most theses were affiliated with universities in Canada (13) and Sweden (12), followed by Iceland (4) and Norway (3). Other countries included the United States (3), United Kingdom (2), Denmark (1) and Germany (1). No theses were retrieved from Russia. This pattern differs from the total number of studies per country (e.g., Norway having more studies than Canada). However, it is also possible that countries differ in terms of how often theses are published as peer-reviewed articles, and therefore excluded as duplicates.
The data reported by the authors was highly heterogeneous for some context variables. For example, information on temperature and precipitation varied in terms of which variable was reported (such as annual mean, seasonal means, means of warmest and coldest months) and from where data were acquired (the study site, the nearest weather station or from gridded meteorological data). To overcome such heterogeneity, and to acquire data on additional variables that are generally not reported in publications, we used several geospatial data layers to extract variables at the evidence point coordinates (see Table 1 for a list of these variables). While this allows for uniformly defined data across all studies, it also comes with some disadvantages. The precision of geographic coordinates varies among studies, and some studies may thus have been assigned ecological context variable values less accurately than others. Further, geospatial data layers can also have systematic errors in data quality across the Arctic [43], and some geospatial layers gave no data for a substantial proportion of the studies due to data layer resolution or coastline placement. The number of studies with NA values ranged from 0 (extent of recent warming) to 209 (vegetation greenness and growing season length).
Additional potential biases identified during the study coding process were related to the definition of tundra and forest and to inconsistency of defining studies. First, the authors of the reviewed studies described forest in various ways, likely leading to various interpretations by the reviewers coding information from the studies. When in doubt, reviewers inspected Google Earth to assess whether the habitat could be considered an open ecotone forest versus boreal forest, potentially creating some subjectivity in the exclusion/inclusion of particular studies. However, we do not consider this a major concern since only 1% of studies were excluded due to the study being carried out in boreal forests. Second, our definition of studies left some space for interpretation, especially from complex multi-site multi-habitat studies. While this was ameliorated by having two authors to check and edit all coded data, we acknowledge that it might affect the repeatability of our study.
Our approach does not use within-study replication to weigh the amount of evidence provided within a study. For instance, larger spatial and temporal extents did not correspond to more studies per article, unless the article included several clearly stated different ecological contexts. Thus, several studies with replicates spread across large spatial extents [12] or long time-frames [44] were coded as one study, similarly to local, short-term studies [45,46].

Conclusions
We present the first synthesis of the circumpolar coverage of the available evidence about the effects of herbivores on Arctic vegetation. We found that there is a growing number of studies in the field of Arctic herbivory, but our analysis revealed clear gaps in existing coverage. Similar to previous studies reporting geographic biases in Arctic research [27], we found that the evidence base is geographically biased towards Fennoscandia and Subarctic regions. This translates to an evidence base that is, in terms of ecological context, biased towards areas that are warmer, wetter, relatively close to the coast and treeline and have moderately changing climatic conditions. This bias has direct consequences for the extent to which our current understanding of herbivore impacts on vegetation can be used to guide policy and management decisions across the Arctic and how to focus future research.

Implication for policy/management
The bias in ecological contexts covered in the current evidence base may limit our ability to generalize the conclusions on herbivore impacts on Arctic vegetation, and more broadly, the management of plant-herbivore interactions. Herbivory by ungulates has repeatedly been suggested to counteract some of the impacts of climate change on tundra vegetation [1,13,[47][48][49], and grazing management has been proposed as a potential management tool [7,12,50]. Impacts of large herbivores on vegetation are context-dependent [18]. For instance, Rangifer tarandus herbivory is related to different vegetation state changes in Subarctic tall-shrub-tundra and High Arctic graminoid-tundra [12,51], while precipitation appears to be an important determinant of the impacts of Rangifer trampling on lichens [52,53]. Thus, generalizations based on an evidence base that is biased towards Subarctic vegetation types and high precipitation, as shown in the current systematic map, are unlikely to be valid throughout the Arctic. Accordingly, for successful implementation of R. tarandus grazing as a vegetation management tool in colder and drier parts of the Arctic, as well as further away from the coast and/or forests, such as continental Siberia and North America, a better understanding of R. tarandus impact on vegetation in such areas is critical.
However, not all geographical biases in the evidence base represent an ecologically meaningful lack of studies. For example, for actively managed Arctic geese [54][55][56] studies were scarce in inland areas, at higher elevations, in areas with no permafrost, areas with high temperature seasonality and on the Eurasian continent. However, in this case, the evidence base generally reflects actual goose habitat use, which is typically concentrated on coastal lowlands in the Arctic [57,58]. Thus, the existing spatial clusters of studies can provide a good base for local and regional goose management. However, the knowledge gained from these clusters of studies should be applied carefully to management of goose grazing in sites with less represented environmental contexts, such as in permafrost-free Iceland or northern Fennoscandia.

Implication for research
A rather worrying result in terms of understanding the combined effects of herbivores and climate change is that the evidence base has a notable bias towards areas with moderate warming. Consequently, even though studies do exist across the range of climate change contexts, general conclusions of herbivore impacts on vegetation are likely to be biased to the knowledge gained from moderate climate changes. Arctic herbivores [59][60][61][62] and vegetation [63,64] are affected by climate change directly and indirectly, e.g., via changes in phenology and food webs. Thus, there are numerous potential trajectories of plantherbivore interactions in a changing climate, highlighting the need for good coverage across the entire gradient of climatic changes. Our systematic map calls for future studies to prioritize the areas experiencing the most pronounced climate changes but also in those that are more resilient (i.e., those changing the least) to increase our broader understanding and management of the effects of herbivores (see Fig. 4B).
By far, most studies on the effects of herbivores on tundra vegetation focused on vertebrate herbivores, while the evidence base for the effects of invertebrate herbivores on tundra plants remains scant. This knowledge gap on the role of Arctic invertebrate herbivores has been emphasized by recent studies [65,66] and has been traditionally attributed to the lower diversity and abundance of invertebrate herbivores at higher latitudes compared to lower latitudes [67,68]. As well, several studies assessed patterns of invertebrate herbivory but not their impacts on plants (e.g., [65,66,[69][70][71][72]), and were therefore excluded from our systematic map. Others, like [48] reported invertebrate herbivory as an additional "disturbance", rather than as the focal herbivore of their study. Most studies on invertebrates focused on outbreaks of defoliating herbivores, mainly in northern Fennoscandia and Canada, east of the Hudson Bay, with only one article from west Greenland [73]. Outbreaks of invertebrate herbivores have been well documented in the forest-tundra ecotone [61,74] and less frequently in tundra [73]. Background invertebrate herbivory (i.e., herbivory at non-outbreak densities) is widespread across the tundra [65,66] but its effects on tundra plants have yet to be thoroughly investigated [70]. Our systematic map also underscores the need for future studies on invertebrate herbivores in the Arctic, in support of other studies that have emphasized the need for long-term monitoring of insect herbivory across the Arctic [75,76].
According to our systematic map, field experiments commonly use exclosures to study the effects of herbivores on tundra vegetation. Consequently, our understanding of the effects of herbivores is largely biased to understanding what happens after herbivores are removed from the system. However, the effects on vegetation of reducing or completely excluding herbivores may not parallel those of increasing herbivore densities [77]. Consequently, exclosure studies to make management recommendations related to herbivores or vegetation should be used with caution. An additional consideration is that, even if most experimental field studies excluded herbivores from experimental plots, they used very different exclosure methods. Methods ranged from fencing off areas to prevent access by larger herbivores (e.g., [78]), to using size-selective exclosures that prevent access by some herbivores but not others (e.g., [47]), or the use of pesticides to exclude invertebrate herbivores (e.g., [79]). Such differences in methodology make syntheses challenging and prevent the exploration of broad-scale patterns and context dependencies in the processes under study, even if the evidence is available (e.g., [80]). Coordinated experiments that use standardized methodologies across broad geographical ranges can partly help overcome this problem [81]. Harmonizing research efforts through the development and use of common protocols will help make results of different studies more comparable, allowing for more robust generalizations across sites and plant-herbivore systems [82,83]. However, the implementation of such protocols at existing study sites alone will not solve the issue of some ecological contexts being underrepresented, given the current spatial distribution of Arctic research documented by this study and others [27]. Targeting inclusion of study sites in the underrepresented areas (e.g., High Arctic, inland areas) would provide much-needed data and should be a priority in the future planning of coordinated large-scale studies.
Surprisingly, we found that several studies had major gaps in the description of study design, such as omitting information on which year the fieldwork was done, coordinates of the field site, or which parts of the data were already published. Furthermore, for most variables there appeared to be no improvement in the reporting standards over time, a pattern also observed by Haddaway and Verhoeven [84]. The lack of basic, well-structured study information needed for comparisons, such as covariates or metadata [84,85], makes individual studies less useful for general knowledge synthesis. This insufficient reporting of methods and results is a common problem in ecology and evolution preventing transparency and repeatability [84,86]. We strongly urge authors of future studies on herbivory (and any other ecological investigation) to give clear descriptions of study designs (maps showing hierarchical sites, etc.) and to clearly state whether some parts of the data are spatial or temporal subsets of (or expansions to) already published data. We further suggest that authors consider adding an appendix providing information similar to Additional file 4, describing the study design, measured variables, and the ecological context of the study.
We mapped the ecological contexts in which research on Arctic herbivory has been conducted and demonstrate how it can be used to identify gaps in our knowledge and direct future studies on herbivory in the Arctic. Yet, the analyses conducted offer only broad patterns. The evidence base can be used as a basis for systematic reviews addressing more detailed questions and zooming in to specific ecological contexts, herbivores, plants or herbivore effect types. Furthermore, our list of context variables is not exhaustive. For instance, climate change driven changes in snow conditions affect both herbivores and plants [87][88][89], and soil temperatures experienced by plants can differ substantially from gridded air temperature data based on measurements 2 m above ground [90]. Data related to these variables were not available in a format that could be used in a circumpolar synthesis, but once they become available their addition would provide additional insights. As we have made the results available through the interactive Shiny App, they are easily available for exploration and to generate further hypotheses.
Finally, a remaining challenge will be keeping the results up to date. A systematic map represents a substantial investment in time spent writing the protocol, searching and summarizing, and communicating results. Yet, with the ever-increasing number of new studies, systematic maps tend to become outdated fairly quickly. The protocol for this systematic map outlines an idea of a system by which new studies could be submitted to the dataset [28], to increase the longevity and relevance of the systematic map. However, given (i) the amount of manual work related to the preparation of the coded data (minor corrections for typos, grouping of variables, verifications for redundancy), and (ii) the increasing availability of new, relevant geospatial layers, we now believe that such an automatic approach is unlikely to provide robust results with new insights. Rather, regular revisits and updated versions of the original systematic map will be needed (Additional file 1).