The role of indirect interspecific effects in the stochastic dynamics of a simple trophic system.

Indirect interspecific effects (IIEs) occur when one species affects another through a third intermediary species. Understanding the role of IIEs in population dynamics is key for predicting community-level impacts of environmental change. Yet, empirically teasing apart IIEs from other interactions and population drivers has proven challenging and data-demanding, particularly in species-rich communities. We used stochastic population models parameterized with long-term time series of individual data to simulate population trajectories and examine IIEs in a simple high-arctic vertebrate food chain consisting of the wild Svalbard reindeer, its scavenger (the Arctic fox) and the barnacle goose, a migratory prey of the fox. We used the simulated population trajectories to explore co-fluctuations between the species within the food chain. Additionally, we adjusted the model in two ways: first, to isolate the impact of fluctuations in the abundance of a species by keeping its abundance constant; and second, to isolate the impact of a trophic interaction on the dynamics of other species by setting the abundance of the influencing species to zero. We found that fluctuations in reindeer carcasses shaped fox abundance fluctuations, which subsequently affected goose population dynamics. Reindeer and goose population growth rates were nevertheless only weakly correlated, probably in part due to demographic and environmental stochasticity, density dependence and lagged dynamics in the geese. However, removing the fluctuations in reindeer abundance or setting reindeer abundance to zero indeed demonstrated strong underlying IIEs on goose population dynamics and extinction probability. This study thus highlights the importance of species interactions, including IIEs, on species coexistence and communities in the long-term, that is beyond immediate effects and covariation in short-term fluctuations.

Mercury-Modulated Immune Responses in Arctic Barnacle Goslings (Branta leucopsis) upon a Viral-Like Immune Challenge.

Historical mining activities in Svalbard (79°N/12°E) have caused local mercury (Hg) contamination. To address the potential immunomodulatory effects of environmental Hg on Arctic organisms, we collected newborn barnacle goslings (Branta leucopsis) and herded them in either a control or mining site, differing in Hg levels. An additional group at the mining site was exposed to extra inorganic Hg(II) via supplementary feed. Hepatic total Hg concentrations differed significantly between the control (0.011 ± 0.002 mg/kg dw), mine (0.043 ± 0.011 mg/kg dw), and supplementary feed (0.713 ± 0.137 mg/kg dw) gosling groups (average ± standard deviation). Upon immune challenge with double-stranded RNA (dsRNA) injection, endpoints for immune responses and oxidative stress were measured after 24 h. Our results indicated that Hg exposure modulated the immune responses in Arctic barnacle goslings upon a viral-like immune challenge. Increased exposure to both environmental as well as supplemental Hg reduced the level of natural antibodies, suggesting impaired humoral immunity. Hg exposure upregulated the expression of proinflammatory genes in the spleen, including inducible nitric oxide synthase (iNOS) and interleukin 18 (IL18), suggesting Hg-induced inflammatory effects. Exposure to Hg also oxidized glutathione (GSH) to glutathione disulfide (GSSG); however, goslings were capable of maintaining the redox balance by de novo synthesis of GSH. These adverse effects on the immune responses indicated that even exposure to low, environmentally relevant levels of Hg might affect immune competence at the individual level and might even increase the susceptibility of the population to infections.

Environmental change reduces body condition, but not population growth, in a high-arctic herbivore.

Environmental change influences fitness-related traits and demographic rates, which in herbivores are often linked to resource-driven variation in body condition. Coupled body condition-demographic responses may therefore be important for herbivore population dynamics in fluctuating environments, such as the Arctic. We applied a transient Life-Table Response Experiment ('transient-LTRE') to demographic data from Svalbard barnacle geese (Branta leucopsis), to quantify their population-dynamic responses to changes in body mass. We partitioned contributions from direct and delayed demographic and body condition-mediated processes to variation in population growth. Declines in body condition (1980-2017), which positively affected reproduction and fledgling survival, had negligible consequences for population growth. Instead, population growth rates were largely reproduction-driven, in part through positive responses to rapidly advancing spring phenology. The virtual lack of body condition-mediated effects indicates that herbivore population dynamics may be more resilient to changing body condition than previously expected, with implications for their persistence under environmental change.

State dependence explains individual variation in nest defence behaviour in a long-lived bird.

Parental care, such as nest or offspring defence, is crucial for offspring survival in many species. Yet, despite its obvious fitness benefits, the level of defence can consistently vary between individuals of the same species. One prominent adaptive explanation for consistent individual differences in behaviours involves state dependency: relatively stable differences in individual state should lead to the emergence of repeatable behavioural variation whereas changes in state should lead to a readjustment of behaviour. Therefore, empirical testing of adaptive state dependence requires longitudinal data where behaviour and state of individuals of the same population are repeatedly measured. Here, we test if variation in states predicts nest defence behaviour (a 'risky' behaviour) in a long-lived species, the barnacle goose Branta leucopsis. Adaptive models have predicted that an individual's residual reproductive value or 'asset' is an important state variable underlying variation in risk-taking behaviour. Hence, we investigate how nest defence varies as a function of time of the season and individual age, two state variables that can vary between and within individuals and determine asset. Repeated measures of nest defence towards a human intruder (flight initiation distance or FID) of females of known age were collected during 15 breeding seasons. Increasing values of FID represent increasing shyness. We found that females strongly and consistently differed in FID within- and between-years. As predicted by theory, females adjusted their behaviour to state by decreasing their FID with season and age. Decomposing these population patterns into within- and between-individual effects showed that the state-dependent change in FID was driven by individual plasticity in FID and that bolder females were more plastic than shyer females. This study shows that nest defence behaviour differs consistently among individuals and is adjusted to individual state in a direction predicted by adaptive personality theory.

Contrasting consequences of climate change for migratory geese: Predation, density dependence and carryover effects offset benefits of high-arctic warming.

Climate change is most rapid in the Arctic, posing both benefits and challenges for migratory herbivores. However, population-dynamic responses to climate change are generally difficult to predict, due to concurrent changes in other trophic levels. Migratory species are also exposed to contrasting climate trends and density regimes over the annual cycle. Thus, determining how climate change impacts their population dynamics requires an understanding of how weather directly or indirectly (through trophic interactions and carryover effects) affects reproduction and survival across migratory stages, while accounting for density dependence. Here, we analyse the overall implications of climate change for a local non-hunted population of high-arctic Svalbard barnacle geese, Branta leucopsis, using 28 years of individual-based data. By identifying the main drivers of reproductive stages (egg production, hatching and fledging) and age-specific survival rates, we quantify their impact on population growth. Recent climate change in Svalbard enhanced egg production and hatching success through positive effects of advanced spring onset (snow melt) and warmer summers (i.e. earlier vegetation green-up) respectively. Contrastingly, there was a strong temporal decline in fledging probability due to increased local abundance of the Arctic fox, the main predator. While weather during the non-breeding season influenced geese through a positive effect of temperature (UK wintering grounds) on adult survival and a positive carryover effect of rainfall (spring stopover site in Norway) on egg production, these covariates showed no temporal trends. However, density-dependent effects occurred throughout the annual cycle, and the steadily increasing total flyway population size caused negative trends in overwinter survival and carryover effects on egg production. The combination of density-dependent processes and direct and indirect climate change effects across life history stages appeared to stabilize local population size. Our study emphasizes the need for holistic approaches when studying population-dynamic responses to global change in migratory species.

High-Arctic family planning: earlier spring onset advances age at first reproduction in barnacle geese.

Quantifying how key life-history traits respond to climatic change is fundamental in understanding and predicting long-term population prospects. Age at first reproduction (AFR), which affects fitness and population dynamics, may be influenced by environmental stochasticity but has rarely been directly linked to climate change. Here, we use a case study from the highly seasonal and stochastic environment in High-Arctic Svalbard, with strong temporal trends in breeding conditions, to test whether rapid climate warming may induce changes in AFR in barnacle geese, Branta leucopsis. Using long-term mark-recapture and reproductive data (1991-2017), we developed a multi-event model to estimate individual AFR (i.e. when goslings are produced). The annual probability of reproducing for the first time was negatively affected by population density but only for 2 year olds, the earliest age of maturity. Furthermore, advanced spring onset (SO) positively influenced the probability of reproducing and even more strongly the probability of reproducing for the first time. Thus, because climate warming has advanced SO by two weeks, this likely led to an earlier AFR by more than doubling the probability of reproducing at 2 years of age. This may, in turn, impact important life-history trade-offs and long-term population trajectories.

Flower-visitor communities of an arcto-alpine plant-Global patterns in species richness, phylogenetic diversity and ecological functioning.

Pollination is an ecosystem function of global importance. Yet, who visits the flower of specific plants, how the composition of these visitors varies in space and time and how such variation translates into pollination services are hard to establish. The use of DNA barcodes allows us to address ecological patterns involving thousands of taxa that are difficult to identify. To clarify the regional variation in the visitor community of a widespread flower resource, we compared the composition of the arthropod community visiting species in the genus Dryas (mountain avens, family Rosaceae), throughout Arctic and high-alpine areas. At each of 15 sites, we sampled Dryas visitors with 100 sticky flower mimics and identified specimens to Barcode Index Numbers (BINs) using a partial sequence of the mitochondrial COI gene. As a measure of ecosystem functioning, we quantified variation in the seed set of Dryas. To test for an association between phylogenetic and functional diversity, we characterized the structure of local visitor communities with both taxonomic and phylogenetic descriptors. In total, we detected 1,360 different BINs, dominated by Diptera and Hymenoptera. The richness of visitors at each site appeared to be driven by local temperature and precipitation. Phylogeographic structure seemed reflective of geological history and mirrored trans-Arctic patterns detected in plants. Seed set success varied widely among sites, with little variation attributable to pollinator species richness. This pattern suggests idiosyncratic associations, with function dominated by few and potentially different taxa at each site. Taken together, our findings illustrate the role of post-glacial history in the assembly of flower-visitor communities in the Arctic and offer insights for understanding how diversity translates into ecosystem functioning.

Density-dependent population dynamics of a high Arctic capital breeder, the barnacle goose.

Density regulation of the population growth rate occurs through negative feedbacks on underlying vital rates, in response to increasing population size. Here, we examine in a capital breeder how vital rates of different life-history stages, their elasticities and population growth rates are affected by changes in population size. We developed an integrated population model for a local population of Svalbard barnacle geese, Branta leucopsis, using counts, reproductive data and individual-based mark-recapture data (1990-2017) to model age class-specific survival, reproduction and number of individuals. Based on these estimates, we quantified the changes in demographic structure and the effect of population size on age class-specific vital rates and elasticities, as well as the population growth rate. Local density regulation at the breeding grounds acted to reduce population growth through negative effects on reproduction; however, population size could not explain substantial variation in survival rates, although there was some support for density-dependent first-year survival. With the use of prospective perturbation analysis of the density-dependent projection matrix, we show that the elasticities to different vital rates changed as population size increased. As population size approached carrying capacity, the influence of reproductive rates and early-life survival on the population growth rate was reduced, whereas the influence of adult survival increased. A retrospective perturbation analysis revealed that density dependence resulted in a positive contribution of reproductive rates, and a negative contribution of the numbers of individuals in the adult age class, to the realised population growth rate. The patterns of density dependence in this population of barnacle geese were different from those recorded in income breeding birds, where density regulation mainly occurs through an effect on early-life survival. This indicates that the population dynamics of capital breeders, such as the barnacle goose, are likely to be more reproduction-driven than is the case for income breeders.

Potential Effect of Migration Strategy on Pollutant Occurrence in Eggs of Arctic Breeding Barnacle Geese ( Branta leucopsis).

Arctic-breeding geese acquire resources for egg production from overwintering grounds, spring stopover sites and breeding grounds, where pollutant exposure may differ. We investigated the effect of migration strategy on pollutant occurrence of lipophilic polychlorinated biphenyls (PCBs) and protein-associated poly- and perfluoroalkyl substances (PFASs) and mercury (Hg) in eggs of herbivorous barnacle geese ( Branta leucopsis) from an island colony on Svalbard. Stable isotopes (δ13C and δ15N) in eggs and vegetation collected along the migration route were similar. Pollutant concentrations in eggs were low, reflecting their terrestrial diet (∑PCB = 1.23 ± 0.80 ng/g ww; ∑PFAS = 1.21 ± 2.97 ng/g ww; Hg = 20.17 ± 7.52 ng/g dw). PCB concentrations in eggs increased with later hatch date, independent of lipid content which also increased over time. Some females may remobilize and transfer more PCBs to their eggs, by delaying migration several weeks, relying on more polluted and stored resources, or being in poor body condition when arriving at the breeding grounds. PFAS and Hg occurrence in eggs did not change throughout the breeding season, suggesting migration has a greater effect on lipophilic pollutants. Pollutant exposure during offspring production in arctic-breeding migrants may result in different profiles, with effects becoming more apparent with increasing trophic levels.

Climate warming may affect the optimal timing of reproduction for migratory geese differently in the low and high Arctic.

Rapid climate warming is driving organisms to advance timing of reproduction with earlier springs, but the rate of advancement shows large variation, even among populations of the same species. In this study, we investigated how the rate of advancement in timing of reproduction with a warming climate varies for barnacle goose (Branta leucopsis) populations breeding at different latitudes in the Arctic. We hypothesized that populations breeding further North are generally more time constrained and, therefore, produce clutches earlier relative to the onset of spring than southern populations. Therefore, with increasing temperatures and a progressive relief of time constraint, we expected latitudinal differences to decrease. For the years 2000-2016, we determined the onset of spring from snow cover data derived from satellite images, and compiled data on egg laying date and reproductive performance in one low-Arctic and two high-Arctic sites. As expected, high-Arctic geese laid their eggs earlier relative to snowmelt than low-Arctic geese. Contrary to expectations, advancement in laying dates was similar in high- and low-Arctic colonies, at a rate of 27% of the advance in date of snowmelt. Although advancement of egg laying did not fully compensate for the advancement of snowmelt, geese laying eggs at intermediate dates in the low Arctic were the most successful breeders. In the high Arctic, however, early nesting geese were the most successful breeders, suggesting that high-Arctic geese have not advanced their laying dates sufficiently to earlier springs. This indicates that high-Arctic geese especially are vulnerable to negative effects of climate warming.

Stress behaviour and physiology of developing Arctic barnacle goslings ( Branta leucopsis) is affected by legacy trace contaminants.

Natural populations are persistently exposed to environmental pollution, which may adversely impact animal physiology and behaviour and even compromise survival. Responding appropriately to any stressor ultimately might tip the scales for survival, as mistimed behaviour and inadequate physiological responses may be detrimental. Yet effects of legacy contamination on immediate physiological and behavioural stress coping abilities during acute stress are virtually unknown. Here, we assessed these effects in barnacle goslings ( Branta leucopsis) at a historical coal mine site in the Arctic. For three weeks we led human-imprinted goslings, collected from nests in unpolluted areas, to feed in an abandoned coal mining area, where they were exposed to trace metals. As control we led their siblings to feed on clean grounds. After submitting both groups to three well-established stress tests (group isolation, individual isolation, on-back restraint), control goslings behaved calmer and excreted lower levels of corticosterone metabolites. Thus, legacy contamination may decisively change stress physiology and behaviour in long-lived vertebrates exposed at a young age.

Greenland Norse walrus exploitation deep into the Arctic.

Walrus ivory was a prized commodity in medieval Europe and was supplied by Norse intermediaries who expanded across the North Atlantic, establishing settlements in Iceland and Greenland. However, the precise sources of the traded ivory have long remained unclear, raising important questions about the sustainability of commercial walrus harvesting, the extent to which Greenland Norse were able to continue mounting their own long-range hunting expeditions, and the degree to which they relied on trading ivory with the various Arctic Indigenous peoples that they were starting to encounter. We use high-resolution genomic sourcing methods to track walrus artifacts back to specific hunting grounds, demonstrating that Greenland Norse obtained ivory from High Arctic waters, especially the North Water Polynya, and possibly from the interior Canadian Arctic. These results substantially expand the assumed range of Greenland Norse ivory harvesting activities and support intriguing archaeological evidence for substantive interactions with Thule Inuit, plus possible encounters with Tuniit (Late Dorset Pre-Inuit).

Extreme events, trophic chain reactions, and shifts in phenotypic selection.

Demographic consequences of rapid environmental change and extreme climatic events (ECEs) can cascade across trophic levels with evolutionary implications that have rarely been explored. Here, we show how an ECE in high Arctic Svalbard triggered a trophic chain reaction, directly or indirectly affecting the demography of both overwintering and migratory vertebrates, ultimately inducing a shift in density-dependent phenotypic selection in migratory geese. A record-breaking rain-on-snow event and ice-locked pastures led to reindeer mass starvation and a population crash, followed by a period of low mortality and population recovery. This caused lagged, long-lasting reductions in reindeer carrion numbers and resultant low abundances of Arctic foxes, a scavenger on reindeer and predator of migratory birds. The associated decrease in Arctic fox predation of goose offspring allowed for a rapid increase in barnacle goose densities. As expected according to r- and K-selection theory, the goose body condition (affecting reproduction and post-fledging survival) maximising Malthusian fitness increased with this shift in population density. Thus, the winter ECE acting on reindeer and their scavenger, the Arctic fox, indirectly selected for higher body condition in migratory geese. This high Arctic study provides rare empirical evidence of links between ECEs, community dynamics and evolution, with implications for our understanding of indirect eco-evolutionary impacts of global change.

No divergence in Cassiope tetragona: persistence of growth response along a latitudinal temperature gradient and under multi-year experimental warming.

BACKGROUND AND AIMS: The dwarf shrub Cassiope tetragona (Arctic bell-heather) is increasingly used for arctic climate reconstructions, the reliability of which depends on the existence of a linear climate-growth relationship. This relationship was examined over a high-arctic to sub-arctic temperature gradient and under multi-year artificial warming at a high-arctic site. METHODS: Growth chronologies of annual shoot length, as well as total leaf length, number of leaves and average leaf length per year, were constructed for three sites. Cassiope tetragona was sampled near its cold tolerance limit at Ny-Ålesund, Svalbard, at its assumed climatic optimum in Endalen, Svalbard, and near its European southern limit at Abisko, Sweden. Together these sites represent the entire temperature gradient of this species. Leaf life span was also determined. Each growing season from 2004 to 2010, 17 open top chambers (OTCs) were placed near Ny-Ålesund, thus increasing the daily mean temperatures by 1·23°C. At the end of the 2010 growing season, shoots were harvested from OTCs and control plots, and growth parameters were measured. KEY RESULTS: All growth parameters, except average leaf length, exhibited a linear positive response (R(2) between 0·63 and 0·91) to mean July temperature over the temperature gradient. Average leaf life span was 1·4 years shorter in sub-arctic Sweden compared with arctic Svalbard. All growth parameters increased in response to the experimental warming; the leaf life span was, however, not significantly affected by OTC warming. CONCLUSIONS: The linear July temperature-growth relationships, as well as the 7 year effect of experimental warming, confirm that the growth parameters annual shoot length, total leaf length and number of leaves per year can reliably be used for monitoring and reconstructing temperature changes. Furthermore, reconstructing July temperature from these parameters is not hampered by divergence.

Keeping the clock set under the midnight sun: diurnal periodicity and synchrony of avian Isospora parasites cycle in the High Arctic.

For Isospora (Protozoa: Eimeriidae) parasites of passerine birds, diurnal periodicity of oocyst output is a well-described phenomenon. From the temporal zone to the tropics, oocyst production is correlated with the light-dark cycle, peaking in the afternoon hours. However, nothing is known about the existence of diurnal periodicity of these parasites in the birds of High Arctic environments, under permanent light during summer. We sampled free-ranging Snow Bunting (Aves: Passeriformes), on Svalbard in summer and tested oocysts output of Isospora plectrophenaxia. Here we show that under the permanent light conditions of Arctic summer in the wild, Isospora plectrophenaxia, a parasite of the Snow Bunting, still keeps the 24-h rhythm of oocyst output with the peak in the post-meridiem hours, despite the absence of diurnal periodicity in host's activity. Our findings prove the ability of avian Isospora to invoke alternative cues for synchronizing the circadian rhythms. Possible cues and adaptive significance of diurnal periodicity of parasite output in High Arctic are discussed. The maintenance of synchronization and timing of the parasite life-cycle stages is under positive selection pressure even in permanent daylight in the Arctic.

Judging a reindeer by its teeth: A user-friendly tooth wear and eruption pattern recording scheme to estimate age-at-death in reindeer (Rangifer tarandus).

Reindeer (Rangifer tarandus) have shaped the cultures and provided livelihood to peoples of the Northern Hemisphere for thousands of years. They are still the socio-economic cornerstone of many northern cultures. Insight into reindeer mortality patterns is important for understanding past human-reindeer interactions and reindeer population fluctuations in relation to climatic and environmental change. Beyond archaeology, assessing the age structures of modern reindeer populations is important for developing wildlife management strategies. This paper presents a quick, non-destructive and cheap method to estimate age in reindeer in both modern and ancient populations based on tooth wear and eruption patterns of mandibular teeth. We devised the method using a large sample of Svalbard reindeer (Rangifer tarandus platyrhynchus) of known age. We blind-tested the method and tested its applicability on another known-age Svalbard reindeer mandible assemblage. The tests demonstrate our methods' user-friendliness and reliability to generate reproducible, reusable datasets and accuracy in estimating reindeer age-at-death.

Seabird-affected taluses are denitrification hotspots and potential N2O emitters in the High Arctic.

In High Arctic tundra ecosystems, seabird colonies create nitrogen cycling hotspots because of bird-derived labile organic matter. However, knowledge about the nitrogen cycle in such ornithocoprophilous tundra is limited. Here, we determined denitrification potentials and in-situ nitrous oxide (N2O) emissions of surface soils on plant-covered taluses under piscivorous seabird cliffs at two sites (BL and ST) near Ny-Ålesund, Svalbard, in the European High Arctic. Talus soils at both locations had very high denitrification potentials at 10 °C (2.62-4.88 mg N kg-1 dry soil h-1), near the mean daily maximum air temperature in July in Ny-Ålesund, with positive temperature responses at 20 °C (Q10 values, 1.6-2.3). The talus soils contained abundant denitrification genes, suggesting that they are denitrification hotspots. However, high in-situ N2O emissions, indicating the presence of both active aerobic nitrification and anaerobic denitrification, were observed only at BL (max. 16.6 µg N m-2 h-1). Rapid nitrogen turnover at BL was supported by lower carbon-to-nitrogen ratios, higher nitrate content, and higher δ15N values in the soils at BL compared with those at ST. These are attributed to the 30-fold larger seabird density at BL than at ST, providing the larger organic matter input.

Serosurvey for Toxoplasma gondii in arctic foxes and possible sources of infection in the high Arctic of Svalbard.

Samples (blood or tissue fluid) from 594 arctic foxes (Alopex lagopus), 390 Svalbard reindeer (Rangifer tarandus platyrhynchus), 361 sibling voles (Microtus rossiaemeridionalis), 17 walruses (Odobenus rosmarus), 149 barnacle geese (Branta leucopsis), 58 kittiwakes (Rissa tridactyla), and 27 glaucous gulls (Larus hyperboreus) from Svalbard and nearby waters were assayed for antibodies against Toxoplasma gondii using a direct agglutination test. The proportion of seropositive animals was 43% in arctic foxes, 7% in barnacle geese, and 6% (1 of 17) in walruses. There were no seropositive Svalbard reindeer, sibling voles, glaucous gulls, or kittiwakes. The prevalence in the arctic fox was relatively high compared to previous reports from canid populations. There are no wild felids in Svalbard and domestic cats are prohibited, and the absence of antibodies against T. gondii among the herbivorous Svalbard reindeer and voles indicates that transmission of the parasite by oocysts is not likely to be an important mechanism in the Svalbard ecosystem. Our results suggest that migratory birds, such as the barnacle goose, may be the most important vectors bringing the parasite to Svalbard. In addition to transmission through infected prey and carrion, the age-seroprevalence profile in the fox population suggests that their infection levels are enhanced by vertical transmission.

Balancing ecosystem function, services and disservices resulting from expanding goose populations.

As goose populations increase in abundance, their influence on ecological processes is increasing. We review the evidence for key ecological functions of wild goose populations in Eurasia and North America, including aquatic invertebrate and plant propagule transport, nutrient deposition in terrestrial and aquatic ecosystems, the influence of goose populations on vegetation biomass, carbon storage and methane emission, species diversity and disease transmission. To estimate the implications of their growing abundance for humans, we explore how these functions contribute to the provision of ecosystem services and disservices. We assess the weight, extent and trends among such impacts, as well as the balance of their value to society. We examine key unresolved issues to enable a more balanced assessment of the economic costs or benefits of migratory geese along their flyways, including the spatial and temporal variation in services and their contrasting value to different user groups. Many ecological functions of geese are concluded to provide neither services nor disservices and, ecosystem disservices currently appear to outweigh services, although this varies between regions. We consider an improved quantification of ecosystem services and disservices, and how these vary along population flyways with respect to variation in valuing certain cultural services, and under different management scenarios aimed at reducing their disservices, essential for a more balanced management of goose populations.