Small birds, big effects: the little auk (Alle alle) transforms high Arctic ecosystems

Ivan González-Bergonzoni, Kasper L. Johansen, Anders Mosbech, Frank Landkildehus, Erik Jeppesen, Thomas A. Davidson
Published 15 February 2017.DOI: 10.1098/rspb.2016.2572

Abstract

In some arctic areas, marine-derived nutrients (MDN) resulting from fish migrations fuel freshwater and terrestrial ecosystems, increasing primary production and biodiversity. Less is known, however, about the role of seabird-MDN in shaping ecosystems. Here, we examine how the most abundant seabird in the North Atlantic, the little auk (Alle alle), alters freshwater and terrestrial ecosystems around the North Water Polynya (NOW) in Greenland. We compare stable isotope ratios (δ15N and δ13C) of freshwater and terrestrial biota, terrestrial vegetation indices and physical–chemical properties, productivity and community structure of fresh waters in catchments with and without little auk colonies. The presence of colonies profoundly alters freshwater and terrestrial ecosystems by providing nutrients and massively enhancing primary production. Based on elevated δ15N in MDN, we estimate that MDN fuels more than 85% of terrestrial and aquatic biomass in bird influenced systems. Furthermore, by using different proxies of bird impact (colony distance, algal δ15N) it is possible to identify a gradient in ecosystem response to increasing bird impact. Little auk impact acidifies the freshwater systems, reducing taxonomic richness of macroinvertebrates and truncating food webs. These results demonstrate that the little auk acts as an ecosystem engineer, transforming ecosystems across a vast region of Northwest Greenland.

1. Introduction

Migratory animals translocate energy and nutrients between ecosystems and may support productivity and biomass in otherwise unproductive systems [1]. Species responsible for such translocation of nutrients are often termed ‘ecosystem engineers’ as they may profoundly change the recipient ecosystem (e.g. [2–4]). For example, Pacific salmon species are responsible for large-scale transport of marine-derived nutrients (MDN) to freshwater and terrestrial ecosystems in temperate and arctic regions of North America and Asia [5–8]. Productivity and biodiversity increase in systems with Pacific salmon as MDN are assimilated into stream biofilms and terrestrial vegetation [7,9].

In many locations around the globe, seabirds feeding at sea and breeding in colonies in terrestrial systems are known to bring nutrients to land [9–13]. However, evidence of the extent to which the seabird-MDN subsidy alters ecosystems is sparse and what exists is limited to detecting the presence of MDN in terrestrial soil, vegetation and a few soil invertebrates [14–17]. These studies rely on nitrogen stable isotope ratios (δ15N), exploiting the fact that marine δ15N is higher than terrestrial δ15N and thus easily traceable in freshwater and terrestrial biota (e.g. [17]). For example, δ15N was found to be higher in soil and vegetation at bird-affected sites compared with bird-free sites for a range of colonial seabird species in Antarctica [16], Florida, USA [18], New Zealand [15,19] and Svalbard [20].

The few studies considering change in ecosystem structure induced by seabirds (e.g. [10,11,13]) suggest that seabird colonies alter terrestrial vegetation structure, increase primary productivity and fuel terrestrial food webs both above and below ground [13,14]. Seabird guano may also alter the chemical properties of fresh waters by changing nutrient concentrations and pH [21]. However, these findings apply to relatively small catchments and are mostly valid only for ecosystems located immediately adjacent to specific bird colonies on islands [13,18,19].

Here, we sought to elucidate the nature and extent of the marine nutrient subsidy from the extensive seabird colonies along the shores of the North Water Polynya (NOW) in Northwest Greenland (figure 1). This area is the main breeding ground of the little auk (Alle alle), a small (approx. 160 g), zooplanktivorous alcid, which is the most abundant seabird in the North Atlantic [22,23]. Within a range of approx. 400 km, approximately 80% of the global little auk population, or 33 million breeding pairs, have their nesting sites in dense colonies on talus slopes up to 10 km inland from the coast [24–26]. Colonies are attended from early-May to mid-August when parents, performing round-trips to at-sea foraging areas up to 100 km from the breeding site, raise a single chick on a diet of lipid-rich Arctic copepods [27]. In the NOW, little auks are estimated to be capable of consuming up to 24% of the copepod standing stock [28], bringing vast quantities of MDN to land.

Figure 1.
Figure 1.

The investigation area in Northwest Greenland encompassing the entire range of little auk breeding colonies associated with the North Water Polynya (NOW). Red areas: little auk breeding colonies after (24). Black dots: sites sampled in 2014–2015. White dots: sites sampled in 2001 and used for comparison. Black squares: permanently inhabited settlements. On the overview map, the solid blue area represents the approximate late winter/early spring extent of the NOW.

The overall aim of the study was to determine the contribution of seabird-MDN to the biomass of primary producers and consumers in both terrestrial and aquatic habitats. Specifically, we sought to identify the effect of bird colonies on terrestrial and freshwater primary productivity, freshwater physical–chemical characteristics and biological community composition and also to investigate the potential mechanisms behind any differences between affected and unaffected systems. We hypothesize that: (i) a very large proportion of the biomass of aquatic and terrestrial primary producers and consumers in bird colony areas is fuelled by MDN; (ii) little auk colonies increase terrestrial and aquatic primary production, and alter physical–chemical properties and community structure in recipient freshwater ecosystems, resulting in increased nutrient concentrations, algal biomass and taxa richness.

2. Material and methods

(a) Study area

The NOW polynya in Smith Sound, Northern Baffin Bay lies between Greenland and Canada. It covers around 85 000 km2 and is the largest polynya in the Arctic [29] (figure 1). The combination of year-round nutrient rich, open waters and constant light in the summer makes the NOW one of the most productive marine areas in the Arctic [30].

(b) Sampling campaigns

In late-July and early-August 2014 and 2015, terrestrial and freshwater ecosystems were sampled along the Greenlandic coastline of NOW from Savissivik in the south to Siorapaluk in the north (figure 1). In the field, sampling sites were classified as either ‘colony’ or ‘control’ sites. Sites located in the drainage catchment of a little auk colony or under a flight corridor of birds commuting between a colony and at-sea foraging areas (areas receiving bird droppings) were classified as colony sites. Sites located in catchments without colonies or overflying little auks were classified as control sites. Additionally, we used data on stable isotopes and nutrient concentrations in lakes sampled in an area without little auks near Pituffik in 2001 (electronic supplementary material, Methods S1; figure 1).

(c) Stable isotope sampling

At each site, samples were collected for analysis of C and N stable isotopes (δ15N and δ13C). Samples were obtained from soil, terrestrial mosses, pooled terrestrial plant leaves, excrement (little auk, geese, arctic hare (Lepus arcticus) and musk ox (Ovibos moschatus)), hair from arctic hare, and a skull from an arctic fox (Vulpes lagopus). In freshwater habitats, filamentous algae, aquatic mosses, debris (conditioned leaf litter), benthic biofilm, macroinvertebrates, seston, zooplankton, profundal lake sediments and fish were sampled. Samples were analysed at UC Davis Stable Isotope Facilities, California, USA (http://stableisotopefacility.ucdavis.edu). Isotopic data from animals were lipid-corrected based on their C : N ratio following eqn 3 of Post et al. [31]. Details of the stable isotope sampling and corrections are given in the electronic supplementary material, Methods S1.

(d) Freshwater consumer taxa richness and physical–chemical properties

The richness of consumer taxa was measured as the number of aquatic consumer taxa obtained during the sampling for stable isotopes, taxa being defined at family level. Family level richness was used as species-level identification of macroinvertebrates was not possible in the field. Family-level richness generally correlates well with species richness, thus constituting a valid measure of richness (e.g. [32]). Sampling effort was similar at the different sites, allowing comparison between sites (see electronic supplementary material, Methods S1 for details). In each aquatic system, we also took water samples for measuring nutrient concentrations and algal biomass samples (chlorophyll-a). Key physical–chemical variables were recorded using a multi-parameter probe (for details, see electronic supplementary material, Methods S1). Information about the physical–chemical characteristics of sites is available in electronic supplementary material, table S2 and details of taxa collected are given in electronic supplementary material, table S3.

(e) Terrestrial productivity

As a measure of terrestrial productivity, we used an enhanced vegetation index (EVI) image from MODIS Terra, 28/7–12/8 2015 [33,34]. Owing to the coarse resolution of the image (250 × 250 m), and the heterogeneous nature of the landscape in the vicinity of the sampling sites (e.g. patches of vegetation in places where soil formation is possible, interspaced with areas of bare rock and water), the maximum EVI value within a 500 m radius of each sampling site was used in the statistical analyses (for details, see electronic supplementary material, Methods S1).

(f) Comparison of colony and control sites

To avoid making type II errors, univariate tests of differences between colony and control sites were preceded by two PERMANOVAs [35], one including δ15N and δ13C of terrestrial and aquatic primary producer and consumer groups, and another including grouped freshwater physical–chemical parameters. In the univariate tests, variance was often heterogeneous and generalized least squares models (GLS, α = 0.05, [36]) were used with the appropriate error structures to account for this. Residual plots were checked for remaining heterogeneity and for spatial autocorrelation [36]. Where residual spatial autocorrelation was detected, a spatial weights matrix was integrated into the model and residuals were re-checked. Full details of the models can be found in electronic supplementary material, Methods S1.

In relation to aquatic habitats, comparisons of parameters between colony and control sites were made for all system types (lakes, ponds and streams) pooled. However, in the case of algal biomass, separate comparisons were made for lotic (streams) and lentic (lakes + ponds) systems due to different units being used: in lotic systems benthic algal biomass was in micrograms per square centimetre whereas in lentic systems phytoplankton biomass was expressed in micrograms per litre.

(g) Estimation of the contribution of seabird-derived nutrient to biomass

Enhanced δ15N at colony sites relative to control sites provides an unequivocal marker of the presence of seabird-MDN. Following the procedure of Harding et al. [11], the contribution of seabird-MDN to the biomass of different terrestrial and aquatic primary producer and consumer groups at colony sites was estimated by means of mass balance models for N. Details of methodology are provided in electronic supplementary material, Methods S1.

(h) Changes along a gradient of bird impact

To study how terrestrial and freshwater ecosystems were affected along a gradient of bird impact, we investigated the relationships between distance to nearest little auk colony and EVI, aquatic algal biomass (lotic and lentic Chl-a) and δ15N of freshwater benthic algae. Further, in a detailed case study of Savissivik Island, where GPS tracking of breeding little auks was conducted, we examined how EVI and freshwater benthic algal δ15N varied in relation to overflight intensity of little auks (proxy of bird dropping intensity) at five sample sites at varying distances from the tracking colony. We also modelled the drainage pattern from the Savissivik colony to evaluate its effect on the spread of nutrients in the landscape. All details are provided in electronic supplementary material, Methods S1.

The combined results of these analyses strongly indicated that benthic algal δ15N is a good proxy of the relative magnitude of bird nutrient input in fresh waters, reflecting true impact much better than distance to nearest little auk colony (see Results and Discussion). We therefore used benthic algal δ15N to detect changes in freshwater physical–chemical and community characteristics (pH, total nitrogen, total phosphorus, algal biomass and consumer taxa richness) along a gradient of bird impact. The use of δ15N as an indicator of MDN input has been supported in diverse studies (e.g. [7,12,21,37]). Details of the models used to test the relationships are provided in electronic supplementary material, Methods S1.

(i) Potential drivers of changes in freshwater ecosystems along a gradient of bird impact

Finally, in order to identify potential mechanisms behind changes in freshwater community structure along a gradient of bird impact, we tested for relationships between environmental variables changing with bird impact, i.e. nutrient concentrations and pH, versus algal biomass (phytoplankton Chl-a in lentic and benthic algal Chl-a in lotic systems) and consumer taxa richness. Statistical procedures are described in electronic supplementary material, Methods S1.

3. Results and discussion

(a) Comparison of colony and control sites

Unequivocal evidence of fertilization by seabird-derived nutrients was reflected in the different isotopic fingerprints of C and N in terrestrial and aquatic primary producers and consumers at little auk colony sites compared with control sites (PERMANOVA F = 9.9; d.f.res = 11, p < 0.01), in particular by the approximate 10-fold difference in their δ15N (GLS: p < 0.05, table 1; electronic supplementary material, figure S1). The greater statistical significance of the difference in δ15N compared with δ13C between colony and control sites reflects the fact that while marine-derived nitrogen is incorporated in both terrestrial and freshwater ecosystems, carbon of marine origin is only incorporated in freshwater systems (e.g. [15]). Specifically aquatic mosses and chironomids were enriched in δ13C at colony sites. The pattern of elevated δ15N at colony sites was significant in terrestrial systems and across all fresh waters—lakes, ponds and streams (table 1; electronic supplementary material, figure S1). Lake sediment, seston, zooplankton and hair from arctic hare (Lepus arcticus) also had higher δ15N values at colony sites, although this could not be tested statistically due to lack of replicates (table 1; electronic supplementary material, figure S1). The observed 10-fold increase in δ15N is larger than that recorded in systems where migratory fish transfer MDN to terrestrial and aquatic ecosystems (i.e. a threefold to fourfold δ15N increase at fish impacted versus control sites [7,37]).

View this table:
Table 1.

Nitrogen and carbon isotopic signatures (δ15N and δ13C, ‰, mean ± s.d.) at colony and control sites, and estimated contribution of marine-derived nitrogen (MDN) to the biomass of primary producers and consumers at colony sites. Calculations of the contribution of MDN to biomass were performed with mass balance models using as marine nitrogen source δ15N of peat at colony sites, adding 2‰ enrichment in the case of aquatic resources and consumers. This 2‰ enrichment represents the hydrolysable proportion of nitrogen reaching aquatic ecosystems. This estimate is tentative as the variability in δ15N reaching freshwaters is high due to diverse microbial processes affecting guano and soil δ15N, implying that values more than 100% may occur. Significant differences in isotopic fingerprints are marked in bold and marginal p-values are given in italics.

In agreement with our first hypothesis, the modelling of the MDN subsidy indicated that an overwhelming majority of both terrestrial and freshwater primary producer and consumer biomass was fuelled by MDN at little auk colony sites (table 1). The mass balance models always yielded values larger than 85% at colony sites. While directly comparable with other studies (e.g. [5,11]), there are uncertainties associated with these estimates due to possible variation in the fractionation of marine nitrogen from guano to its final uptake product, which is dependent on various microbial processes [11]. For example, the process of conversion of uric acid to ammonia involves the volatilization of ammonia, a powerful fractionation process leaving the remaining substrate enriched by approximately 40‰ in δ15N [38]. By contrast, the fractionation of nitrogen during nitrification depletes δ15N of the substrate with about −25‰ [39]. For seston and zooplankton, the uncertainty is higher due to low sample size from colony sites (table 1). Notwithstanding these uncertainties, the data provide strong evidence of a very large MDN subsidy of terrestrial and aquatic ecosystem production at colony sites. This is highlighted by the fact that the proportions of MDN assimilated into freshwater organisms described here (always more than 85%) are much higher than those reported for biota in New Zealand streams related to seabird colonies using directly comparable methods (28 to 38% of biomass generated from MDN) [11]. Our proportions are also much higher than those estimated in MDN subsidy studies of Pacific salmon (e.g. 23 and 25% of marine-derived biomass in aquatic organisms and terrestrial vegetation, respectively; [7]). The higher values compared with other studies probably reflect both the large quantity of the MDN input in little auk colonies and the paucity of other nutrient sources at these high latitudes (76–78 deg. N).

In agreement with our second hypothesis, we found the physical–chemical characteristics of freshwater systems differed significantly between colony and control sites (PERMANOVA F = 8.8, d.f.res = 26; p < 0.01; table 2). Nutrient concentrations were significantly higher at colony sites, the only exception being the marginal significance of phosphate concentrations (GLS: p = 0.07; table 2). Algal biomass was also significantly higher at colony sites—approximately 20-fold for phytoplankton biomass in lentic systems and approximately 10-fold for benthic algal biomass in lotic systems (GLS: p < 0.05; table 2). The nutrient levels and algal biomass observed in the aquatic systems at colony sites are the highest reported for Greenland, where most systems are characterized by nutrient limitation (e.g. [40,41]). Correspondingly, in terrestrial systems, there were significantly higher EVI values (approximately twofold higher) at colony sites (GLS: p < 0.01; table 2). However, contrary to expectations, freshwater consumer taxa richness was lowest at the nutrient enriched, bird-impacted sites (GLS: p < 0.05; table 2; electronic supplementary material, table 3 in supplementary appendix S2).

View this table:
Table 2.

Terrestrial vegetation, physical–chemical and biotic characteristics of freshwater systems in catchments with and without little auk colonies (colony versus control sites). Values are given as mean ± s.d. PERMANOVA and pairwise GLS test parameters are given to allow comparison of each parameter between colony and control sites. Significant differences are marked in bold and marginal p-values are given in italics. nt = not tested due to lack of replicates. Generally, the GLS tests were conducted for all system types pooled (lakes, ponds streams). However, differences in algal biomasses were tested separately for lentic and lotic systems due to different measurement methods (see Material and methods).

Among the chemical properties of the fresh waters, pH differed significantly with colony sites being more acidic than control sites (GLS: p < 0.05; table 2). This effect appears to be a particular characteristic of little auk colonies, contrasting with findings from Devon Island, Canada, where pools under a northern fulmar (Fulmarus glacialis) colony exhibited increased pH values relative to control sites without birds [21]. In Svalbard, it has been observed that zooplanktivorous seabirds, such as the little auk, promote soil acidification, whereas piscivorous seabirds do not [42].

(b) Changes along a gradient of bird impact

At a broad scale, a decrease in the distance to the nearest little auk colony was associated with an increase in EVI (GLS: t = −2.6 p < 0.05, pseudo r2 = 0.06, n = 29), benthic algae biomass (lotic Chl-a) (GLS: t = −4.7 p < 0.0001, pseudo r2 = 0.07, n = 17) and δ15N of freshwater benthic algae (GLS: t = −5.6 p < 0.00001, pseudo r2 = 0.29, n = 27), whereas the relationship with phytoplankton biomass (lentic Chl-a) was not significant (electronic supplementary material, figure S2). When considering only sites closer than 2500 m from colonies the relationships became more clear: EVI increased strongly with proximity to colony (GLS: t = −5.1 p < 0.001, pseudo r2 = 0.52, n = 27) as did benthic algal biomass (GLS: t = −4.5 p < 0.0001, pseudo r2 = 0.12, n = 15) and benthic algal δ15N (GLS: t = −4.6 p < 0.0001, pseudo r2 = 0.17, n = 21) (electronic supplementary material, figure S2). However, from the scatterplots it is evident that distance to colony is a relatively poor predictor as a site can be proximal to a colony and remain relatively unaffected by MDN.

An explanation for this is provided by a case study of Savissivik Island, where, due to GPS tracking, we were in the unique position of being able to relate EVI and δ15N of benthic algae with an estimate of overflight intensity of birds (proxy of bird dropping intensity) (figure 2). Acknowledging that we only have five sample sites on Savissivik Island, these indicate strong, positive correlations between overflight intensity of little auks and both EVI (LM: p < 0.01; r2 = 0.88; n = 5) and δ15N of freshwater benthic algae (LM: p < 0.01; r2 = 0.88; n = 5). Further, benthic algal δ15N and EVI were strongly, positively correlated at the sample sites (LM: p < 0.001; r2 = 0.94; n = 5). The drainage pattern from the colony seemingly also corresponds with the EVI and benthic algal δ15N values at the sample sites (figure 2). Thus, while it is clear that the little auk colony is the main source of the nutrients on Savissivik Island, the relative importance of overflight versus drainage in the spread of nutrients could not be discerned.

Figure 2.
Figure 2.

Relationships between flight intensity of little auks, enhanced vegetation index (EVI) and δ15N of freshwater benthic algae at Savissivik Island where GPS tracking of little auks was conducted. (a) Relative flight intensity of little auks (kilometre of track line per square kilometre) on a colour scale from blue (low) to red (high). Coastline depicted as black line, little auk colony as hatched area, and deployment site used during GPS tracking as a white dot. Freshwater sampling sites from Savissivik Island (NOW 9–13) are labelled on the map and the symbol size is scaled according to the benthic algal δ15N value of the site. (b) EVI from MODIS Terra, 28/7–12/8 2015 (33), on a colour scale from blue (low) to red (high). The black, hatched region extending from the little auk colony indicates drainage from the colony, modelled on the basis of a digital elevation model of Savissivik Island (see electronic supplementary material, Methods S1, for details).

The Savissivik case clearly demonstrates that little auks use distinct flight paths when commuting between their breeding colony and offshore foraging areas, and that local drainage systems transport the nutrients in particular directions. Thus, it is possible to have a low impact site close to a large colony, if the site is located outside the flight path and upstream of the drainage from the colony/flight path. This appears to be the main reason for the inadequacy of distance to colony as a predictor of aquatic and terrestrial primary producer biomass. The EVI evidence suggests that relatively fertile terrestrial areas do exist outside the influence of little auk colonies, primarily in conjunction with meadows. However, EVI values more than 0.25 are almost exclusively found in association with little auk colonies, and terrestrial productivity is clearly elevated over large areas of the Greenlandic coastal forelands of the NOW due to the extensive little auk colonies (figure 3).

Figure 3.
Figure 3.

Terrestrial productivity is elevated in little auk colonies. Little auk colonies (black polygons; after (24)) and enhanced vegetation index (EVI) values below 0.25 (light green) and above 0.25 (dark green) in the northern part of the investigation area. The EVI image is from MODIS Terra, 28/7–12/8 2015 (33). Despite many EVI values missing along the coastline due to mixed pixels, the map clearly shows that EVI values above 0.25 are almost exclusively found in association with little auk colonies. In inserts (a–d), landscape photos illustrate the contrasting terrestrial productivity at colony and control sites. (a) Productive landscape in the catchment of the little auk colony on Savissivik Island (site 9); (c) productive landscape in the catchment of the little auk colony at Annikitsoq on the south coast of Cape York (sites 25–32); (b) barren landscape at the control site on Savissivik Island (site 10); (d) barren landscape at the control site close to Booth Sound (site 15–16). See figure 1 and electronic supplementary material, table S2 for exact positions.

The fact that freshwater benthic algal δ15N values were approximately 15-fold higher at colony sites compared with control sites, decrease with distance to colony, and, in the case of Savissivik island, are tightly coupled to overflight intensity and drainage input from the colony, indicated that it is a robust indicator of the intensity of MDN impact. This is in line with findings in Canadian seabird colonies [12,21] and relationships observed between δ15N of terrestrial vegetation and distance to a salmon stream [7,37] or relative salmon carcass density [37]. Thus, benthic algal δ15N was employed as a proxy of little auk impact to investigate changes of freshwater physical–chemical properties and biotic community structure along a gradient of impact. As δ15N of benthic algae rose there was a significant increase in concentrations of total nitrogen (GAM: F = 2.5, p < 0.0001, r2 = 0.63, n = 34) and phosphorus (GAM: F = 3.5, p < 0.0001, r2 = 0.36, n = 34) and system acidity (LM: p < 0.0001, r2 = 0.53, n = 32) (figure 4). With increasing benthic algal δ15N, algal biomass rose in both lentic (GAM: F = 2.1, p < 0.05, r2 = 0.39) and lotic systems (GAM: F = 2.8, p < 0.001, r2 = 0.55), whereas freshwater consumer taxa richness decreased (GAM: F = 10.8, p < 0.01, r2 = 0.21) (figure 4).

Figure 4.
Figure 4.

Changes in environmental and biotic characteristics induced by increasing little auk impact (using δ15N of benthic algae as a proxy of impact). Sampling sites in catchments with little auk colonies are marked with black dots and sites in control areas with open circles. (a) Decrease in pH with increasing bird influence. (b) Increase in total nitrogen concentrations with increasing bird influence. (c) Increase in total phosphorus concentrations with increasing bird influence. (d) Increase in phytoplankton biomass in lentic systems with increasing bird influence. (e) Increase in stream benthic algal biomass with increasing bird influence. (f) Decrease in consumer taxa richness with increasing bird influence. Model details are shown in each panel figure.

(c) Potential drivers of changes in freshwater ecosystems along a gradient of bird impact

We attribute the increase in nutrients to be the main driver of increased algal biomass in lakes and streams (electronic supplementary material, figure S3). In both lotic and lentic systems, the increase in total nitrogen concentrations promoted increased algal biomass (LMs: p < 0.0001, r2 = 0.77, n = 16 and p < 0.0001, r2 = 0.89, n = 11 for lotic and lentic systems, respectively) as did the total phosphorus concentration, which was also strongly related to algal biomass in both the lentic (LM: p < 0.0001, r2 = 0.99, n = 16) and lotic systems (LM: p < 0.0001, r2 = 0.75, n = 11) (electronic supplementary material, figure S3). In nutrient-poor systems, increased nutrient concentrations enhance primary producer biomass, resulting in bottom-up effects [13,43] that increase the abundance and biomass of primary and secondary consumers [4], often creating higher taxonomic richness [44]. This matches our field observations of terrestrial systems where sites located below bird colonies were the most productive and greenest with most observations of foxes, hares, geese and muskoxen, whereas control sites were largely barren (figure 3). However, in the freshwater systems, the enhanced algal biomass with increasing bird impact was associated with decreasing taxonomic richness. The acidification associated with little auk impact is a potential driver as pH was found to be negatively correlated with consumer taxa richness in both lotic (GAM: p < 0.001, r2 = 0.44, n = 17) and lentic systems (LM: p = 0.06, r2 = 0.49, n = 12) (electronic supplementary material, figure S3). Some bird-impacted sites had extremely low pH, down to 3.4, and such low levels cause biodiversity loss in lakes and streams [45,46]. At two colony sites, the harsh environment with low pH values meant that neither fish nor macroinvertebrates were found during the sampling (electronic supplementary material, table S3).

4. Conclusion

In the study region, the approximately 60–70 million breeding little auks [26] alter both terrestrial and freshwater ecosystems by promoting primary and secondary production. While these findings are similar to studies of MDN subsidy produced by migrating Pacific salmon species, the magnitude of the MDN subsidy reported here is much higher and the consequences for freshwater consumers differ significantly. In association with little auk impact on freshwater systems, we found a decrease in species richness of higher consumers and truncated food webs without fish. This highlights the key relevance of the identity of the vector of nutrient subsidies in order to understand and predict ecosystem-wide consequences of engineer species that translocate nutrients between ecosystems.

As the total horizontal extent of breeding colonies is approximately 400 km [24], a significant proportion the coastal forelands around the NOW has been transformed by this single species: the little auk. Similar significant changes in the terrestrial environment related to the presence of seabird colonies have also been reported for islands in Alaska where in total more than 10 million birds nest [10]. Here, the introduction of arctic fox (Vulpes lagopus) on some islands in the late nineteenth century resulted in decreased bird abundance, reducing the nutrient subsidies by seabirds to terrestrial productivity, and consequently the landscape shifted from grasslands to tundra [10]. During the breeding season, little auks depend on lipid-rich copepods species associated with cold water, and consequently it has been suggested that little auk populations will decline in response to the current warming of the Arctic [47,48]. If so, a landscape shift comparable to the one observed in [10] may be expected for the NOW.

Ethics

The procedures used conform to the legal requirements of the country and institutional guidelines.

Data accessibility

The datasets supporting this article have been uploaded as part of the electronic supplementary material.

Authors' contributions

I.G.-B., A.M., K.L.J., E.J. and T.A.D. participated in the conception and design of the study. All authors carried out the fieldwork. I.G.-B., T.A.D. and K.L.J. analysed the data. I.G.-B. carried out laboratory work and drafted the manuscript. A.M., K.L.J., F.L., E.J. and T.A.D. helped draft the manuscript. T.A.D. coordinated the study. All authors gave final approval for publication.

Competing interests

We have no competing interests

Funding

This study is part of The North Water Project (NOW.KU.DK) funded by the Velux Foundation, the Villum Foundation and the Carlsberg Foundation of Denmark. E.J. was further supported by the MARS project (Managing Aquatic ecosystems and water Resources under multiple Stress) funded under the 7th EU Framework Programme, Theme 6 (Environment including Climate Change), contract no.: 603378 (http://www.mars-project.eu). I.G.-B. was supported by SNI (Agencia Nacional de Investigación e Innovación, ANII, Uruguay).

Acknowledgements

We thank Anne Mette Poulsen for manuscript editing, and we wish to extend our warmest thanks for help during fieldwork to the local communities Siorapaluk, Savissivik and Qaanaaq, and, at Thule Air Base, to liaison officer Kim Mikkelsen, and Tony Rønne Pedersen and Erland Søndergaard of Greenland Contractors. We would also like to thank the crew of the ships Minna Martek, Blue Jay and Hot Totty.

Footnotes

  • Received November 20, 2016.
  • Accepted January 25, 2017.
http://royalsocietypublishing.org/licence

Published by the Royal Society. All rights reserved.

References

  1. ↵
    1. Polis GA,
    2. Anderson WB,
    3. Holt RD
    . 1997 Toward an integration of landscape and food web ecology: the dynamics of spatially subsidized food webs. Ann. Rev. Ecol. Sys. 28, 289–316. (doi:10.2307/2952495)
    OpenUrl
  2. ↵
    1. KB Gido,
    2. D Jackson
    1. Flecker AS,
    2. McIntyre PB,
    3. Moore JW,
    4. Anderson JT,
    5. Taylor BW,
    6. Hall Jr RO
    . 2010 Migratory fishes as material and process subsidies in riverine ecosystems. In Community ecology of stream fishes: concepts, approaches, and techniques (eds KB Gido, D Jackson), pp. 559–592. Bethesda, MD: American Fisheries Society, Symposium.
    1. Jones CG,
    2. Lawton JH,
    3. Shachak M
    . 1994 Organisms as ecosystem engineers. Oikos 69, 373–386. (doi:10.2307/3545850)
    OpenUrlCrossRefWeb of Science
  4. ↵
    1. Polis GA,
    2. Power MA,
    3. Huxel GR
    . 2004 Food webs at the landscape level. Chicago, IL: University of Chicago Press.
  5. ↵
    1. Chaloner DT,
    2. Wipfli MS
    . 2002 Influence of decomposing Pacific salmon carcasses on macroinvertebrate growth and standing stock in southeastern Alaska streams. J. N. Am. Bent. Soc. 21, 430–442. (doi:10.2307/1468480)
    OpenUrl
    1. Gende SM,
    2. Edwards RT,
    3. Willson MF,
    4. Wipfli MS
    . 2002 Pacific salmon in aquatic and terrestrial ecosystems: Pacific salmon subsidize freshwater and terrestrial ecosystems through several pathways, which generates unique management and conservation issues but also provides valuable research opportunities. BioScience 52, 917–928. (doi:10.1641/0006-3568(2002)052[0917:psiaat]2.0.co;2)
    OpenUrlFREE Full Text
  7. ↵
    1. Koshino Y,
    2. Kudo H,
    3. Kaeriyama M
    . 2013 Stable isotope evidence indicates the incorporation into Japanese catchments of marine-derived nutrients transported by spawning Pacific salmon. Freshw. Biol. 58, 1864–1877. (doi:10.1111/fwb.12175)
    OpenUrl
  8. ↵
    1. Zhang Y,
    2. Negishi JN,
    3. Richardson JS,
    4. Kolodziejczyk R
    . 2003 Impacts of marine-derived nutrients on stream ecosystem functioning. Proc. R. Soc. Lond. B 270, 2117–2123. (doi:10.1098/rspb.2003.2478)
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Field RD,
    2. Reynolds JD
    . 2011 Sea to sky: impacts of residual salmon-derived nutrients on estuarine breeding bird communities. Proc. R. Soc. B 278, 3081–3088. (doi:10.1098/rspb.2010.2731)
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Croll DA,
    2. Maron JL,
    3. Estes JA,
    4. Danner EM,
    5. Byrd GV
    . 2005 Introduced predators transform subarctic islands from grassland to tundra. Science 307, 1959–1961. (doi:10.1126/science.1108485)
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. Harding JS,
    2. Hawke DJ,
    3. Holdaway RN,
    4. Winterbourn MJ
    . 2004 Incorporation of marine-derived nutrients from petrel breeding colonies into stream food webs. Freshw. Biol. 49, 576–586. (doi:10.1111/j.1365-2427.2004.01210.x)
    OpenUrl
  12. ↵
    1. Michelutti N,
    2. Keatley BE,
    3. Brimble S,
    4. Blais JM,
    5. Liu H,
    6. Douglas MSV,
    7. Mallory ML,
    8. Macdonald RW,
    9. Smol JP
    . 2009 Seabird-driven shifts in Arctic pond ecosystems. Proc. R. Soc. B 276, 591–596. (doi:10.1098/rspb.2008.1103)
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Sánchez-Piñero F,
    2. Polis GA
    . 2000 Bottom-up dynamics of allochthonous input: direct and indirect effects of seabirds on islands. Ecology 81, 3117–3132. (doi:10.1890/0012-9658(2000)081[3117:BUDOAI]2.0.CO;2)
    OpenUrlWeb of Science
  14. ↵
    1. Callaham Jr MA,
    2. Butt KR,
    3. Lowe CN
    . 2012 Stable isotope evidence for marine-derived avian inputs of nitrogen into soil, vegetation, and earthworms on the isle of Rum, Scotland, UK. Eur. J. Soil Biol. 52, 78–83. (doi:10.1016/j.ejsobi.2012.07.004)
    OpenUrl
  15. ↵
    1. Hawke DJ,
    2. Newman J
    . 2007 Carbon-13 and nitrogen-15 enrichment in coastal forest foliage from nutrient-poor and seabird enriched sites in southern New Zealand, NZ. New Zealand J. Botany 45, 309–315. (doi:10.1080/00288250709509719)
    OpenUrl
  16. ↵
    1. Huang T,
    2. Sun L,
    3. Wang Y,
    4. Chu Z,
    5. Qin X,
    6. Yang L
    . 2014 Transport of nutrients and contaminants from ocean to island by emperor penguins from Amanda Bay, East Antarctic. Sci. Total Environ. 468–469, 578–583. (doi:10.1016/j.scitotenv.2013.08.082)
    OpenUrl
  17. ↵
    1. Wainright SC,
    2. Haney JC,
    3. Kerr C,
    4. Golovkin AN,
    5. Flint MV
    . 1998 Utilization of nitrogen derived from seabird guano by terrestrial and marine plants at St Paul, Pribilof Islands, Bering Sea, Alaska. Mar. Biol. 131, 63–71. (doi:10.1007/s002270050297)
    OpenUrl
  18. ↵
    1. Irick DL,
    2. Gu B,
    3. Li YC,
    4. Inglett PW,
    5. Frederick PC,
    6. Ross MS,
    7. Wright AL,
    8. Ewe SML
    . 2015 Wading bird guano enrichment of soil nutrients in tree islands of the Florida Everglades. Sci. Total Environ. 532, 40–47. (doi:10.1016/j.scitotenv.2015.05.097)
    OpenUrl
  19. ↵
    1. Markwell TJ,
    2. Daugherty CH
    . 2002 Invertebrate and lizard abundance is greater on seabird-inhabited islands than on seabird-free islands in the Marlborough Sounds, New Zealand. Ecoscience 9, 293–299. (doi:10.1080/11956860.2002.11682715)
    OpenUrl
  20. ↵
    1. Zmudczyńska-Skarbek K,
    2. Balazy P,
    3. Kuklinski P
    . 2015 An assessment of seabird influence on Arctic coastal benthic communities. J. Mar. Syst. 144, 48–56. (doi:10.1016/j.jmarsys.2014.11.013)
    OpenUrl
  21. ↵
    1. Keatley BE,
    2. Douglas MSV,
    3. Blais JM,
    4. Mallory ML,
    5. Smol JP
    . 2008 Impacts of seabird-derived nutrients on water quality and diatom assemblages from Cape Vera, Devon Island, Canadian High Arctic. Hydrobiologia 621, 191–205. (doi:10.1007/s10750-008-9670-z)
    OpenUrl
  22. ↵
    1. Barrett RB,
    2. Chapdelaine G,
    3. Anker-Nilssen T,
    4. Mosbech A,
    5. Montevecch IWA,
    6. Reid J,
    7. Veit R
    . 2006 Seabird numbers and prey consumption in the North Atlantic. ICES J. Mar. Sci. 63, 1145–1158. (doi:10.1016/j.icesjms.2006.04.004)
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Stempniewicz L
    . 2001 Alle alle Little auk. BWP Update. Birds of the Western Palearctic 3, 175–201.
    OpenUrl
  24. ↵
    1. Boertmann D,
    2. Mosbech A
    . 1998 Distribution of little auk (Alle alle) breeding colonies in Thule District, Northwest Greenland. Polar Biol. 19, 206–210. (doi:10.1007/s003000050236)
    OpenUrl
    1. Kampp K,
    2. Falk K,
    3. Pedersen CE
    . 2000 Breeding density and population of little auks (Alle alle) in a Northwest Greenland colony. Polar Biol. 23, 517–521. (doi:10.1007/s003000000115)
    OpenUrlCrossRef
  26. ↵
    1. Egevang C,
    2. Boertmann D,
    3. Mosbech A,
    4. Tamstorf MP
    . 2003 Estimating colony area and population size of little auks (Alle alle) at Northumberland Island using aerial images. Polar Biol. 26, 8–13. (doi:10.1525/pol.2003.26.2.8)
    OpenUrl
  27. ↵
    1. Frandsen M,
    2. Fort J,
    3. Rigét FF,
    4. Galatius A,
    5. Mosbech A
    . 2014 Composition of chick meals from one of the main little auk (Alle alle) breeding colonies in Northwest Greenland. Polar Biol. 37, 1055–1060. (doi:10.1007/s00300-014-1491-0)
    OpenUrl
  28. ↵
    1. Karnovsky N,
    2. Hunt GL
    . 2002 Estimation of carbon flux to dovekies (Alle alle) in the North Water. Deep-Sea Res. PT II 49, 5117–5130. (doi:10.1016/S0967-0645(02)00181-9)
    OpenUrl
  29. ↵
    1. Stirling I
    . 1980 The biological importance of polynyas in the Canadian Arctic. Arctic 33, 303–315. (doi:10.14430/arctic2563)
    OpenUrl
  30. ↵
    1. Melling H,
    2. Gratton Y,
    3. Ingram G
    . 2001 Ocean circulation within the North Water Polynya of Baffin Bay. Atmos. Ocean 39, 301–325. (doi:10.1080/07055900.2001.9649683)
    OpenUrl
  31. ↵
    1. Post DM,
    2. Layman CA,
    3. Arrington DA,
    4. Takimoto G,
    5. Quattrochi J,
    6. Montaña CG
    . 2007 Getting to the fat of the matter: models, methods and assumptions for dealing with lipids in stable isotope analyses. Oecologia 152, 179–189. (doi:10.1007/s00442-006-0630-x)
    OpenUrlCrossRefPubMedWeb of Science
  32. ↵
    1. Balmford A,
    2. Green MJB,
    3. Murray MG
    . 1996 Using higher-taxon richness as a surrogate for species richness: I. Regional tests. Proc. R. Soc. Lond. B 263, 1267–1274. (doi:10.1098/rspb.1996.0186)
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Didan K
    . 2015 MOD13Q1 MODIS/Terra Vegetation Indices 16-Day L3 Global 250 m SIN Grid V006. (NASA EOSDIS Land Processes DAAC, https://lpdaac.usgs.gov/dataset_discovery/modis/modis_products_table/mod13q1
  34. ↵
    1. Huete A,
    2. Didan K,
    3. Miura T,
    4. Rodriguez EP,
    5. Gao X,
    6. Ferreira LG
    . 2002 Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sensing of Environment 83, 195–213. (doi:10.1016/S0034-4257(02)00096-2)
    OpenUrlCrossRef
  35. ↵
    1. Anderson MJ
    . 2001 A new method for non-parametric multivariate analysis of variance. Austr. Ecol. 26, 32–46.
    OpenUrlCrossRef
  36. ↵
    1. Zuur AF,
    2. Ieno EN,
    3. Walker N,
    4. Saveliev AA,
    5. Smith GM
    . 2009 Mixed effects models and extensions in ecology with R. New York, NY: Springer.
  37. ↵
    1. J Stockner
    1. Reimchen TE,
    2. Mathewson D,
    3. Hocking MD,
    4. Moran J,
    5. Harris D
    . 2002 Isotopic evidence for nrichment of salmon- derived nutrients in vegetation, soil and insects in riparian zones in coastal British Columbia. In Nutrients in salmonid ecosystems: sustaining production and biodiversity (ed. J Stockner), pp. 59–69. Bethesda, MD: American Fisheries Society Symposium.
  38. ↵
    1. Mizutani H,
    2. Hasegawa H,
    3. Wada E
    . 1986 High nitrogen isotope ratio for soils of seabird rookeries. Biogeochemistry 2, 221–247. (doi:10.1007/BF02180160)
    OpenUrlGeoRef
  39. ↵
    1. Robinson D
    . 2001 δ15N as an integrator of the nitrogen cycle. TREE 16, 153–162. (doi:10.1016/S0169-5347(00)02098-X)
    OpenUrl
  40. ↵
    1. Anderson NJ,
    2. Bennike O,
    3. Christoffersen K,
    4. Jeppesen E,
    5. Markager S,
    6. Miller G,
    7. Renberg I
    . 1999 Limnological and palaeolimnological studies of lakes in South-western Greenland. Geol. Greenland Survey Bull. 183, 68–73.
    OpenUrl
  41. ↵
    1. González-Bergonzoni I,
    2. Landkildehus F,
    3. Meerhoff M,
    4. Lauridsen TL,
    5. Özkan K,
    6. Davidson TA,
    7. Mazzeo N,
    8. Jeppesen E
    . 2014 Fish determine macroinvertebrate food webs and assemblage structure in Greenland subarctic streams. Freshw. Biol. 59, 1830–1842. (doi:10.1111/fwb.12386)
    OpenUrl
  42. ↵
    1. Zwolicki A,
    2. Małgorzata K,
    3. Zmudczyńska-Skarbek K,
    4. Iliszko L,
    5. Stempniewicz L
    . 2013 Guano deposition and nutrient enrichment in the vicinity of planktivorous and piscivorous seabird colonies in Spitsbergen. Polar Biol. 36, 363–372. (doi:10.1007/s00300-012-1265-5)
    OpenUrl
  43. ↵
    1. Huryn AD
    . 1998 Ecosystem-level evidence for top-down and bottom-up control of production in a grassland stream system. Oecologia 115, 173–183. (doi:10.1007/s004420050505)
    OpenUrlCrossRefWeb of Science
  44. ↵
    1. VanderMeulen MA,
    2. Hudson AJ,
    3. Scheiner SM
    . 2001 Three evolutionary hypotheses for the hump-shaped productivity–diversity curve. Evol. Ecol. Res. 3, 379–392.
    OpenUrl
  45. ↵
    1. Layer K,
    2. Hildrew AG,
    3. Jenkins GB,
    4. Riede J,
    5. Rossiter SJ,
    6. Townsend CR,
    7. Woodward G
    . 2011 Long-term dynamics of a well-characterised food web: four decades of acidification and recovery in the Broadstone Stream model system. Adv. Ecol. Res. 44, 69–117. (doi:10.1016/B978-0-12-374794-5.00002-X)
    OpenUrlCrossRef
  46. ↵
    1. Schindler DW
    . 1990 Experimental perturbations of whole lakes as tests of hypotheses concerning ecosystem structure and function. Oikos 57, 25–41. (doi:10.2307/3565733)
    OpenUrlCrossRefWeb of Science
  47. ↵
    1. Grémillet D,
    2. Fort J,
    3. Amélineau F,
    4. Zakharova E,
    5. Le Bot T,
    6. Sala E,
    7. Gavrilo M
    . 2015 Arctic warming: nonlinear impacts of sea-ice and glacier melt on seabird foraging. Glob. Change Biol. 21, 1116–1123. (doi:10.1111/gcb.12811)
    OpenUrl
  48. ↵
    1. Jakubas D,
    2. Trudnowska E,
    3. Wojczulanis-Jakubas K,
    4. Iliszko L,
    5. Kidawa D,
    6. Darecki M,
    7. Błachowiak-Samołyk K,
    8. Stempniewicz L
    . 2013 Foraging closer to the colony leads to faster growth in little auks. MEPS 489, 263–278. (doi:10.3354/meps10414)
    OpenUrl
View Abstract
Back to top
22 February 2017
Volume 284, issue 1849
Proceedings of the Royal Society B: Biological Sciences: 284 (1849)

Keywords

marine-derived nutrients
nutrient subsidies
stable isotopes
arctic food webs
ecosystem engineer
seabird colonies
Email
Print
Manage alerts
Citation tools
Download

Article reuse

Get Permission

See related subject areas:

Related articles

Cited by