<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Mar. Sci.</journal-id>
<journal-title>Frontiers in Marine Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Mar. Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-7745</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmars.2024.1343966</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Marine Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Seabird-vectored pelagic nutrients integrated into temperate intertidal rocky shores</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Healing</surname><given-names>Samuel</given-names>
</name>
<xref ref-type="author-notes" rid="fn001"><sup>*</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2287229"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Benkwitt</surname><given-names>Cassandra E.</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/2362620"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Dunn</surname><given-names>Ruth E.</given-names>
</name>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Graham</surname><given-names>Nicholas A. J.</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/1261884"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
</contrib>
</contrib-group>
<aff id="aff1"><institution>Lancaster Environment Centre, Lancaster University</institution>, <addr-line>Lancaster</addr-line>, <country>United Kingdom</country></aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Vitor H. Paiva, University of Coimbra, Portugal</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Gustavo Martins, University of the Azores, Portugal</p>
<p>Ryan Andrades, Federal University of Espirito Santo, Brazil</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Samuel Healing, <email xlink:href="mailto:samuelthealing@gmail.com">samuelthealing@gmail.com</email>
</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>01</day>
<month>05</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>11</volume>
<elocation-id>1343966</elocation-id>
<history>
<date date-type="received">
<day>24</day>
<month>11</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>08</day>
<month>04</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Healing, Benkwitt, Dunn and Graham</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Healing, Benkwitt, Dunn and Graham</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Seabirds provide inter-habitat connectivity by transporting vast quantities of nutrients from their pelagic feeding grounds to terrestrial and marine ecosystems via their nitrogen and phosphorous rich guano. However, it remains unclear whether seabird nutrients are cycling through many nearshore, temperate systems, such as United Kingdom (UK) rocky shores. Furthermore, it also remains unknown how seabird nutrient subsidies impact biodiversity, productivity, and recolonisation rates in UK rocky intertidal systems. Here, the impact of seabird nutrient subsidies to the Farne Islands, Northumberland, areas of high seabird densities during their breeding seasons, was compared to nearby low seabird density areas on the Northumberland coast. Estimated seabird derived nitrogen inputs and seabird densities were 420-1,025 and 949-2,279 times higher, respectively, on seabird islands than on the mainland. Seabird derived nutrient inputs led to substantially enriched nitrogen stable isotope (&#x3b4;<sup>15</sup>N) values and total nitrogen content in terrestrial soil, plants, and lichen, and intertidal barnacles, and limpets, and higher &#x3b4;<sup>15</sup>N values in macroalgae and turf algae. The seabird islands contained marginally greater algal species richness and higher densities of intertidal predators. However, other drivers could not be ruled out, such as nutrient inputs from mainland sources and marginally higher wave exposure at the mainland sites. No difference in limpet size, biomass, or abundance existed between the seabird islands and the mainland. On plots artificially cleared to simulate a disturbance event, barnacles recolonised faster on mainland sites. These findings provide evidence that seabird vectored pelagic nutrients are incorporated into terrestrial and intertidal ecosystems of UK islands, and therefore play an important role in facilitating ecosystem connectivity.</p>
</abstract>
<kwd-group>
<kwd>cross-ecosystem nutrients</kwd>
<kwd>intertidal environment</kwd>
<kwd>nutrient subsidy</kwd>
<kwd>seabird guano</kwd>
<kwd>marine ecology</kwd>
<kwd>stable isotope analysis</kwd>
</kwd-group>
<counts>
<fig-count count="7"/>
<table-count count="1"/>
<equation-count count="2"/>
<ref-count count="98"/>
<page-count count="16"/>
<word-count count="8869"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Marine Ecosystem Ecology</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>The vectoring of nutrients and organic matter via the movement of motile consumers has long been recognised as a mechanism that links ecosystems (<xref ref-type="bibr" rid="B41">Hutchinson, 1950</xref>). Allochthonous nutrient subsidies support many ecosystems, influencing productivity, functioning, diversity, biomass, and community dynamics, boosting food availability and quality, and enhancing ecosystem resilience (<xref ref-type="bibr" rid="B76">Polis et&#xa0;al., 1997</xref>; <xref ref-type="bibr" rid="B58">Lundberg and Moberg, 2003</xref>; <xref ref-type="bibr" rid="B56">Loreau and Holt, 2004</xref>; <xref ref-type="bibr" rid="B12">Benkwitt et&#xa0;al., 2019</xref>). Seabirds are important nutrient vectors, transporting nutrients from their pelagic ocean feeding grounds to coastlines and islands via marine prey items that drop on the colony, lost feathers, abandoned or destroyed eggs, and individuals that die on land. However, the greatest nutrient subsidy is in the form of nitrogen- and phosphorous-rich guano; nutrients that are often limiting resources across many ecosystems. At local scales, these levels of nutrient input dictate the structure and dynamics of communities to the extent that seabirds have been described as &#x201c;keystone species&#x201d; (<xref ref-type="bibr" rid="B66">Mulder et&#xa0;al., 2011b</xref>). At global scales, seabirds may link marine and terrestrial environments, contributing to inter-habitat connectivity and providing multiple ecosystem services (<xref ref-type="bibr" rid="B71">Otero et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B85">Signa et&#xa0;al., 2021</xref>).</p>
<p>Whilst the functional role of seabird nutrient inputs on terrestrial habitats has received a great deal of attention, empirical studies investigating the effects of bottom-up nutrients supplied by seabirds on coastal marine ecosystems remain comparatively limited (<xref ref-type="bibr" rid="B49">Kazama, 2019</xref>; <xref ref-type="bibr" rid="B85">Signa et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B33">Grant et&#xa0;al., 2022</xref>). Evidence does exist, though, of nutrient inputs from seabirds moving through systems, such as consistently higher nitrogen signals in primary producers and invertebrates that are in close proximity to seabird colonies (<xref ref-type="bibr" rid="B47">Jones, 2010</xref>; <xref ref-type="bibr" rid="B50">Kolb et&#xa0;al., 2010</xref>; <xref ref-type="bibr" rid="B31">Gagnon et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B78">Rankin and Jones, 2021</xref>). However, whilst it is important to understand the extent that seabird nutrient inputs propagate through ecosystems, it is also critical to establish whether these allochthonous nutrients impact biodiversity and productivity. These influences have received recent attention in tropical systems: seabird nutrients led to increases in the growth rates of coral, damselfish, and parrotfish, as well as the biomass of fish communities (<xref ref-type="bibr" rid="B32">Graham et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B81">Savage, 2019</xref>; <xref ref-type="bibr" rid="B10">Benkwitt et&#xa0;al., 2021a</xref>, <xref ref-type="bibr" rid="B11">Benkwitt et&#xa0;al., 2021b</xref>, <xref ref-type="bibr" rid="B9">Benkwitt et&#xa0;al., 2023</xref>). Some research on the bottom-up effects of seabirds on temperate environments has been performed, examining intertidal community structure, and seabird nutrients impacts on macroalgae. For example, on the Mercury Islands, New Zealand, macroalgal diversity was highest on islands with the longest history of seabird presence (<xref ref-type="bibr" rid="B78">Rankin and Jones, 2021</xref>), whilst in the Stockholm archipelago, there was a shift towards increased epiphytic algae and invertebrate biomass near seabird islands (<xref ref-type="bibr" rid="B50">Kolb et&#xa0;al., 2010</xref>). Further, in South Africa, higher levels of nutrients due to seabirds resulted in higher algal production (<xref ref-type="bibr" rid="B18">Bosman and Hockey, 1986</xref>) and growth rates (<xref ref-type="bibr" rid="B17">Bosman et&#xa0;al., 1986</xref>). Whilst research in both tropical and temperate environments exists, more is needed to understand the influence of seabird nutrients to recipient systems, with the United Kingdom (UK) amongst the countries where research is deficient.</p>
<p>Despite their importance, seabirds are one of the most threatened bird groups, with almost half of all seabird species threatened with extinction (<xref ref-type="bibr" rid="B43">IUCN, 2022</xref>), and over half of Britain and Ireland&#x2019;s seabirds are in decline (<xref ref-type="bibr" rid="B46">JNCC, 2023</xref>). The threats to seabirds do not only impact their own populations, but could also cascade through ecosystems, threatening communities and food webs that are dependent upon seabirds&#x2019; nutrients (<xref ref-type="bibr" rid="B61">McCann et&#xa0;al., 2005</xref>; <xref ref-type="bibr" rid="B80">Rooney et&#xa0;al., 2006</xref>; <xref ref-type="bibr" rid="B89">Thoresen et&#xa0;al., 2017</xref>). Current anthropogenic threats to seabirds include commercial fishing, pollution, habitat degradation, human disturbance, climate change, severe weather events, and human-introduced invasive predators (<xref ref-type="bibr" rid="B27">Croxall et&#xa0;al., 2012</xref>). Indeed, invasive predator eradication programmes are becoming an increasingly important conservation technique, with approximately half of all successful European invasive alien vertebrate eradications on inhabited islands having occurred within the UK, including on Staple island of the Farne Islands, Northumberland (<xref ref-type="bibr" rid="B86">Stanbury et&#xa0;al., 2017</xref>).</p>
<p>Cold-water systems, such as those surrounding the UK, make up a large proportion of global oceans, yet the impact of seabird nutrient input in these environments remains understudied. Strong seasonal differences exist in these environments, and they are prone to stormy weather, which is predicted to increase under future climate predictions (<xref ref-type="bibr" rid="B42">IPCC, 2022</xref>). The UK is of international importance for breeding assemblages of seabirds and shorebirds, with many of its islands designated as Special Protection Areas under EU legislation (<xref ref-type="bibr" rid="B86">Stanbury et&#xa0;al., 2017</xref>). These breeding seabirds excrete vast quantities of nutrients along the UK&#x2019;s coastlines (<xref ref-type="bibr" rid="B79">Riddick et&#xa0;al., 2012</xref>). Given the global threats to seabird populations, that cold-water systems make up a large proportion of global oceans, and the international importance of the UK&#x2019;s seabirds, it is important to determine the influence of seabird nutrient input on UK nearshore marine systems for both regional and global understanding.</p>
<p>To test the impact of seabirds upon terrestrial and intertidal systems in the UK, islands with high densities of breeding seabirds (the Farne Islands) were compared to mainland sites on the adjacent Northumberland coast, where no seabirds breed. It was hypothesised that seabird nutrient input would cycle through the terrestrial and intertidal ecosystems of the Farne Islands, boosting intertidal biodiversity. We also predicted that seabird nutrient input would increase intertidal productivity through increased limpet size and biomass. Rocky shore intertidal systems are subject to regular disturbance and recolonisation by intertidal organisms, which are expected to increase under future climate change. Furthermore, experimentally adding nutrients to rocky intertidal habitats has been shown to increase recolonisation (<xref ref-type="bibr" rid="B52">Kraufvelin et&#xa0;al., 2006</xref>; <xref ref-type="bibr" rid="B51">Kraufvelin, 2007</xref>; <xref ref-type="bibr" rid="B60">Masterson et&#xa0;al., 2008</xref>). Therefore, following experimental clearances it was predicted seabird nutrient input would support a more rapid recolonisation of algae and barnacles recolonisation.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>Methods</title>
<sec id="s2_1">
<label>2.1</label>
<title>Study sites</title>
<p>The Farne Islands are an archipelago situated off the Northumberland coast, located in the North Sea (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1</bold></xref>; 55.63&#xb0; N, 1.63&#xb0; W). The Farne Islands have National Nature Reserve, Special Protection Area, and Site of Special Scientific Interest (SSSI) designations and have been under the care of the National Trust since 1925. The three study islands were Inner Farne, Staple, and Brownsman, all of which have large seasonal seabird populations. Staple island underwent successful brown rat (<italic>Rattus norvegicus</italic>) eradications in 2006 and 2009. To the best of our knowledge, Inner Farne and Brownsman have no history of rat invasion. All data were collected at one site on each island, with partial observational data also collected at a fourth site (a second location on Inner Farne; <xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1</bold></xref>). The four comparison mainland sites were situated within two regions on the Northumberland coast, which, compared to the Farne Islands, have negligible seabird populations. Two sites were located at Blackrocks Point and Harkness Rocks, around Bamburgh lighthouse (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1</bold></xref>; 55.62&#xb0; N, 1.72&#xb0; W). This region has both Northumberland Coast SSSI and Bamburgh Coast and Hills SSSI designations. A further two sites were located at Greenhill Rocks and Monks House Rocks (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1</bold></xref>; 55.59&#xb0; N, 1.68&#xb0; W), which has Northumberland Coast SSSI designation. As on the Farne Islands, all data were collected at three sites, and partial observational data were also collected at a fourth site. All sites had similar rocky shore headland habitats; the seabird islands and mainland sites were as geographically close to one another as possible, had a low level of incline, and similar geology (<xref ref-type="bibr" rid="B13">BGS, 1896</xref>). Despite attempts to reduce the impact of confounding variables, the mainland sites may have been subject to nutrient enrichment from agricultural run-off and sewage treatment facilities, which have been reported at Lindisfarne National Nature reserve, approximately eight and ten kilometres north of the two mainland study sites, respectively. Furthermore, external influences from the North Sea and the River Tweed have been reported as sources of nitrogen within the area (<xref ref-type="bibr" rid="B69">NCC, 2020</xref>). The River Tweed is approximately 25&#xa0;km north of the mainland study sites. The Ballantine exposure scales (<xref ref-type="bibr" rid="B5">Ballantine, 1961</xref>) was used to describe the wave exposure levels for the study sites, where one on the scale describes extremely exposed sites and eight describes extremely sheltered sites. On the Farne Islands, two sites scored seven on the Ballantine exposure scale, and one site scored six; all mainland sites were five on the scale.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Map of the study sites. Location of the seabird islands (blue) and mainland (red) sites within the North Sea, with a close-up of the Farne Islands and Northumberland mainland.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1343966-g001.tif"/>
</fig>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Seabird surveys</title>
<p>Breeding seabird densities on the Farne Islands were counted annually from 2017-2021 by the <xref ref-type="bibr" rid="B67">National Trust (2022a)</xref>. Species-specific survey methods were used, comprising surveys for apparently occupied nests (AON), apparently occupied burrows (AOB), apparently occupied sites (AOS), and whole-colony census counts (<xref ref-type="bibr" rid="B91">Walsh et&#xa0;al., 1995</xref>; <xref ref-type="bibr" rid="B14">Bibby, 2000</xref>). Repeated surveys were performed throughout single breeding seasons for all species to verify counts.</p>
<p>Non-breeding seabird densities on the mainland were counted annually from 2015-2019 by the British Trust Ornithology (BTO) as part of the Wetland Bird Survey (<xref ref-type="bibr" rid="B22">BTO, 2022a</xref>; note: data from seabird surveys during 2020 and 2021 were excluded due to disruptions caused by COVID-19). A mean seabird density for the entire survey area was estimated as all four mainland sites were located within a single BTO survey sector.</p>
<p>Specific species included in the Farne Island and mainland surveys were: Atlantic puffins (<italic>Fratercula arctica</italic>); black-legged kittiwakes (<italic>Rissa tridactyla</italic>); common guillemots (<italic>Uria aalge</italic>); European shags (<italic>Gulosus aristotelis</italic>); great cormorants (<italic>Phalacrocorax carbo</italic>); northern fulmars (<italic>Fulmarus glacialis</italic>); razorbills (<italic>Alca torda</italic>); black<italic>-</italic>headed gulls (<italic>Chroicocephalus ridibundus</italic>); common gulls (<italic>Larus canus</italic>); European herring gulls (<italic>Larus argentatus</italic>); great black-backed gulls (<italic>Larus marinus</italic>); lesser black<italic>-</italic>backed gulls (<italic>Larus fuscus</italic>); little gulls (<italic>Hydrocoloeus minutus</italic>), Arctic terns (<italic>Sterna paradisaea</italic>); common terns (<italic>Sterna hirundo</italic>); little terns (<italic>Sternula albifrons</italic>); roseate terns (<italic>Sterna dougallii</italic>); sandwich terns (<italic>Thalasseus sandvicensis</italic>). Non-seabird species were excluded from counts as they likely do not contribute to allochthonous nutrient inputs. Importantly, these seabird species feed in the open ocean, and not in the intertidal, and thus are not expected to exert top-down control on rocky intertidal systems. In addition, shorebird predation was presumed to exert a negligible influence on limpet populations on both the Farne Islands and the mainland due population sizes being comparable between the Farne Islands and the mainland.</p>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Total estimated nitrogen input from guano</title>
<p>Total nitrogen input from guano per hectare per year was estimated for each site using previously published methods (<xref ref-type="bibr" rid="B96">Young et&#xa0;al., 2010</xref>; <xref ref-type="bibr" rid="B32">Graham et&#xa0;al., 2018</xref>):</p>
<disp-formula>
<mml:math display="block" id="M1">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mtext>NI</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>ij</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mtext>=&#xa0;</mml:mtext>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mtext>N</mml:mtext>
<mml:mtext>g</mml:mtext>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:mtext>&#xa0;x&#xa0;Dr</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>i&#xa0;</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:mtext>x&#xa0;Bd</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>ij</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:mtext>&#xa0;x&#xa0;P</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>ij</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mtext>Area</mml:mtext>
</mml:mrow>
<mml:mtext>j</mml:mtext>
</mml:msub>
<mml:mtext>&#xa0;</mml:mtext>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where nitrogen input per hectare per year (NI) is estimated from the nitrogen content of guano (N<sub>g</sub>), the defecation rate in g per species of bird (<italic>i</italic>) per day (Dr), the number of that species of bird (Bd) on the site (<italic>j</italic>), the number of days of the year that the species is present on the site (P), and the area of the survey site (Area). Nitrogen content of guano was held at 15.46%, as the mean of other seabird species (<xref ref-type="bibr" rid="B55">Lindeboom, 1984</xref>; <xref ref-type="bibr" rid="B87">Staunton Smith and Johnson, 1995</xref>; <xref ref-type="bibr" rid="B96">Young et&#xa0;al., 2010</xref>). Defecation rate was based on the European herring gull and scaled for all species by species&#x2019; biomass, assuming allometric relationships with body size (<xref ref-type="bibr" rid="B77">Portnoy and Soukup, 1990</xref>; <xref ref-type="bibr" rid="B96">Young et&#xa0;al., 2010</xref>). The number of days a year that a species was present was the sum of each species&#x2019; incubation period, fledging period, and colony-based post-fledging care period (in days; <xref ref-type="bibr" rid="B82">Schreiber and Burger, 2001</xref>). In instances of data deficiency, the number of post-fledging care days of the most closely related bird species was used (or a mean of closely related species, where multiple species were present).</p>
<p>The area of mainland survey sites and each seabird island was calculated in ArcGIS Pro (<xref ref-type="bibr" rid="B30">Esri, 2022</xref>; version 2.5.2). The mainland&#x2019;s transect area was provided by the <xref ref-type="bibr" rid="B23">BTO (2022b)</xref>, with all four sites falling within a single BTO sector, estimated to be 722&#xa0;ha. The Farne Islands were estimated using data from Operational Land Imager (OLI) Collection 2 Level-1, Landsat 8 satellite imagery (30&#xa0;m spatial resolution; <xref ref-type="bibr" rid="B90">USGS, 2022</xref>). The seabird islands were estimated to be 7.3&#xa0;ha, 5&#xa0;ha, and 6.9&#xa0;ha for Inner Farne, Staple, and Brownsman, respectively. Note, these areas refer to the size of the areas in which bird survey areas were undertaken, not the size of the site in which the data were collected.</p>
</sec>
<sec id="s2_4">
<label>2.4</label>
<title>Isotope sampling</title>
<p>Nitrogen is commonly used in stable isotope analysis, where the ratio of <sup>15</sup>N to <sup>14</sup>N is compared, expressed as &#x3b4;<sup>15</sup>N. Stable isotope analysis provides trophic level estimates, with greater &#x3b4;<sup>15</sup>N present in predator diets, such as seabirds, due to their diets of fish, squids, and other relatively high trophic level species (<xref ref-type="bibr" rid="B38">Hobson and Welch, 1992</xref>). &#x3b4;<sup>15</sup>N and total nitrogen content (%N) were used to infer the uptake of nutrients through terrestrial and intertidal food chains, using samples from three high seabird density island sites (seabird islands) and three low seabird density mainland sites (the mainland) to determine spatial variation. These samples were collected during the early seabird breeding season (4-8 May 2022).</p>
<p>Terrestrial samples included soil, terrestrial plants, and lichen. To collect topsoil, loose leaf litter and vegetation was cleared to expose the soil. New-growth leaves were taken from non-legume plants that use the C3 photosynthetic pathway to fix carbon dioxide during photosynthesis: sea campion (<italic>Silene maritima</italic>), dock (<italic>Rumex</italic> spp.), nettles (<italic>Urtica dioica</italic>), fiddlenecks (<italic>Amsinckia</italic> spp.), small bugloss (<italic>Anchusa arvensis</italic>), common orache (<italic>Atriplex patula</italic>), or red goosefoot (<italic>Chenopodium rubrum</italic>). The lichen collected was common orange lichen (<italic>Xanthoria parietina</italic>). Intertidal samples included macroalgae, epiphytic algae, turf algae, barnacles, and limpets. For the intertidal primary producers, macroalgae species were egg wrack (<italic>Ascophyllum nodosum</italic>), bladder wrack (<italic>Fucus vesiculosus</italic>), or serrated wrack (<italic>Fucus serratus</italic>); epiphytic algae species were red algae (<italic>Polysiphonia</italic> spp.) or tiny wrack bush (<italic>Elachista fucicola</italic>); turf algae species were pepper dulse (<italic>Osmundea</italic> spp.), tufted coral weed (<italic>Corallina</italic> spp.), green branched weeds (<italic>Cladophera</italic> spp.), sea lettuce (<italic>Ulva</italic> spp.), Irish moss (<italic>Chrondus crispus</italic>), or false Irish moss (<italic>Mastocarpus stellatus)</italic>. For the intertidal consumers, we collected a mollusc, the common limpet (<italic>Patella vulgata</italic>) and barnacles, the acorn barnacle (<italic>Semibalanus balanoides</italic>). Five samples of each terrestrial and intertidal sample type were taken, except for epiphytic algae on Brownsman, due to its absence. Samples were taken a minimum of 5&#xa0;m apart from one another to ensure each sample was unique.</p>
<p>All samples were dried in a food dehydrator at 60&#xb0;C for 36&#xa0;h, or until fully dry, and stored in plastic sample vials. Soil, lichen, macroalgae, epiphytic algae, turf algae, and barnacles were acid washed with hydrochloric acid to dissolve any calcareous matter or sediments that may have contaminated the samples, whereas limpets were extracted from their shells. Stable isotope analysis of nitrogen for all samples was carried out at Lancaster Environment Centre, Lancaster University, UK, using an Elementar Vario MICRO cube Elementar Analyser and analysed in an Isoprime 100 Isotope Ratio Mass Spectrometer, with international standards IAEA 600 and USGS 41. Isotope ratios and %N were calculated from 2 to 91 mg of each sample added to tin capsules. Accuracy based on international standards was 0.2 per mil standard deviation, with selected samples run in duplicate or triplicate to ensure precision of readings.</p>
</sec>
<sec id="s2_5">
<label>2.5</label>
<title>Community composition</title>
<p>Intertidal community composition was surveyed at four seabird island and four mainland study sites. Five replicate 0.5 m<sup>2</sup> quadrats were positioned a minimum of 5&#xa0;m apart at each site in the mid-shore intertidal zone. To determine algae and sessile filter feeder community structure, all algae and sessile filter feeder species were identified and their percent coverage estimated by two observers, where the final value was the mean of the two estimates; the percent coverage of all vacant rocky surface was estimated using the same method. All invertebrate fauna were identified and counted to calculate their abundance. To determine the relative similarity of intertidal community structure between the rocky shores at the mainland sites and seabird islands, intertidal species were classed into four functional groups based on their ecology: macroalgae, turf algae, intertidal invertebrate herbivores (e.g., periwinkles, but excluding limpets) and intertidal invertebrate predators (e.g., dogwhelks). Limpets were excluded from this analysis as they were analysed separately (see section 2.6). To compare species diversity and taxonomic relatedness, diversity measures were calculated for both algal and faunal taxa to determine species richness (total taxa), Shannon&#x2019;s diversity index (H&#x2019;; <xref ref-type="bibr" rid="B83">Shannon, 1948</xref>), and Pielou&#x2019;s evenness index (J&#x2019;; <xref ref-type="bibr" rid="B74">Pielou, 1966</xref>).</p>
</sec>
<sec id="s2_6">
<label>2.6</label>
<title>Limpet size, abundance, and biomass</title>
<p>The size distribution of 100 common limpet individuals at each site was compared between four mainland and four seabird island sites during the early breeding season visit. Vernier callipers were used to measure limpet maximum shell length to an accuracy of 1&#xa0;mm. Limpet abundance (count per 0.5 m<sup>2</sup>) was counted in the same five quadrats per site used for the community analysis data (see section 2.5). As the quantity of algae grazed is more closely related to the biomass of limpets than the number present, limpet biomass was calculated using the following mass/length regression equation (<xref ref-type="bibr" rid="B20">Branch, 1971</xref>; <xref ref-type="bibr" rid="B19">Bosman and Hockey, 1988</xref>):</p>
<disp-formula>
<mml:math display="block" id="M2">
<mml:mrow>
<mml:mtext>M</mml:mtext>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>=</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mn>1.9</mml:mn>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xd7;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mn>10</mml:mn>
</mml:mrow>
<mml:mn>6</mml:mn>
</mml:msup>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xd7;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msup>
<mml:mtext>L</mml:mtext>
<mml:mrow>
<mml:mn>3.5</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>;</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where dry somatic and gonadal mass (M) is a factor of the length in mm (L). Limpet biomass was site specifically calculated using the size distribution of the 100 limpets measured at each respective site, with limpets within each quadrat assumed to show the same pattern of size as the 100 individuals measured. Biomass was then calculated as a factor of limpet length and abundance.</p>
</sec>
<sec id="s2_7">
<label>2.7</label>
<title>Recolonisation patterns</title>
<p>To determine re-colonisation rates following a recent experimental disturbance, five 15 cm<sup>2</sup> plots were artificially cleared in the mid-shore intertidal zone on three mainland sites and three seabird islands. Prior to clearance, percentage cover of algal species and sessile filter feeders was determined, and all invertebrate intertidal fauna were identified and counted. Plots were then cleared of all algae, sessile filter feeders, and common limpets. Plot clearances took place during 4-8 May 2022, around the period of acorn barnacles settlement (<xref ref-type="bibr" rid="B35">Hawkins and Hartnoll, 1982</xref>), using paint scrapers and wire brushes, leaving a substratum of bare rock. The plots were checked at two eight-week intervals (totalling 16 weeks) to determine algae and sessile filter feeder recolonisation levels (percent cover), and the number of all invertebrates present (counts).</p>
</sec>
<sec id="s2_8">
<label>2.8</label>
<title>Statistical analysis</title>
<p>To test for differences in the responses of terrestrial and intertidal communities around high seabird density islands versus low seabird density mainland sites, a combination of univariate and multivariate statistics were used. Linear mixed models (LMMs) were used to analyse univariate responses to seabird presence for &#x3b4;<sup>15</sup>N and %N, with site included as a random effect to account for spatial non-independence among samples.</p>
<p>A multivariate non-metric multidimensional scaling (NMDS) was performed on Bray-Curtis similarity matrices to analyse algae and sessile filter feeder percentage coverage data, as well as the intertidal fauna count data (<xref ref-type="bibr" rid="B53">Kruskal, 1964</xref>). Multivariate analyses were performed on untransformed data, to determine the contributions of the dominant species, and on presence-absence data, to determine the contribution of intermediate and rarer species (<xref ref-type="bibr" rid="B26">Clarke et&#xa0;al., 2014</xref>).</p>
<p>PERMANOVA tests were conducted to test for community differences in intertidal species (<xref ref-type="bibr" rid="B2">Anderson and Walsh, 2013</xref>). Seabird islands were compared with the mainland (considered as a fixed factor) using data averaged by site to account for any non-independence of multiple samples at the same site. Indicator species analysis (IndVal; <xref ref-type="bibr" rid="B29">Dufr&#xea;ne and Legendre, 1997</xref>) can be used to help determine when one or several species characterise a habitat. Here, it was used to help determine which species drove similarities between communities and was performed on both percentage coverage and count data. LMM univariate tests were then performed on IndVal selected algae and intertidal fauna species on both the untransformed and presence-absence data, with site included as a random effect. LMMs were also used to analyse univariate responses to seabird presence in limpet lengths, abundances, and biomasses, again with site included as a random effect to account for spatial non-independence among samples.</p>
<p>Univariate analyses were not performed when there were extreme unequal variances between mainland and seabird island sites (e.g., abundance of intertidal dogwhelks, <italic>Nucella lapillus</italic>). All univariate analyses were assessed for normality and homogeneity, with square root transformations used when those assumptions were not met. This was the case for individual species densities of serrated wrack, pepper dulse, mussels (<italic>Mytilus edulis</italic>), and common winkles (<italic>Littorina littorea</italic>), as well as the &#x2018;turf algae&#x2019; and &#x2018;other herbivores&#x2019; functional groups, and limpet biomass. Likelihood ratio tests were used to determine the effect of seabird density (low versus high) on each response (<xref ref-type="bibr" rid="B98">Zuur et&#xa0;al., 2009</xref>).</p>
<p>To determine whether the presence of seabirds influenced recolonisation responses of barnacles (percent coverage) and limpets (counts), the effect of seabirds (high versus low seabird abundances), time (prior to clearance of plots, zero weeks, eight weeks, and 16 weeks), and their interaction were analysed. As above, site was included as a random effect to account for spatial non-independence among samples (<xref ref-type="bibr" rid="B98">Zuur et&#xa0;al., 2009</xref>), and post-hoc tests with a Bonferroni correction performed to determine the drivers of differences (<xref ref-type="bibr" rid="B39">Holm, 1979</xref>).</p>
<p>All statistical analyses were conducted in R version 4.0.2 with associated packages ggplot2, dplyr, vegan, ape, indicspecies, rstatix, and lme4 (<xref ref-type="bibr" rid="B28">De Caceres and Legendre, 2009</xref>; <xref ref-type="bibr" rid="B7">Bates et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B93">Wickham, 2016</xref>; <xref ref-type="bibr" rid="B94">Wickham et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B73">Paradis and Schliep, 2019</xref>; <xref ref-type="bibr" rid="B70">Oksanen et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B48">Kassambara, 2021</xref>).</p>
<p>All results were interpreted using the language suggested by <xref ref-type="bibr" rid="B65">Muff et&#xa0;al. (2022)</xref>.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<sec id="s3_1">
<label>3.1</label>
<title>Seabird density and estimated nitrogen input</title>
<p>Mean seabird density, averaged across a five-year period, was 1563, 2279, and 949 times greater for seabird islands Inner Farne, Staple, and Brownsman, respectively, than on the mainland (<xref ref-type="fig" rid="f2"><bold>Figure&#xa0;2A</bold></xref>; mean &#xb1; standard error (SE); 3778.0 &#xb1; 125.1, 5509.3 &#xb1; 232.2, and 2334.3 &#xb1; 93.1 birds ha<sup>-1</sup> for Inner Farne, Staple, and Brownsman, respectively; 2 &#xb1; 0.4 birds ha<sup>-1</sup> for the mainland). Mean estimated seabird nitrogen input, averaged across a five-year period, was 669, 1025, and 420 times greater for Inner Farne, Staple, and Brownsman, respectively, than on the mainland (<xref ref-type="fig" rid="f2"><bold>Figure&#xa0;2B</bold></xref>; mean &#xb1; SE; 923.68 &#xb1; 49.8, 1414.6 &#xb1; 76.6, and 579.4 &#xb1; 21.2 kg<sup>-1</sup> ha<sup>-1</sup> for Inner Farne, Staple, and Brownsman, respectively; 1.38 &#xb1; 0.5 kg<sup>-1</sup> ha<sup>-1</sup> for the mainland).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Seabird densities and nitrogen input to seabird islands (Inner Farne, Staple, and Brownsman, the Farne Islands) and low seabird density mainland sites (Northumberland coast). <bold>(A)</bold> Seabird density on Inner Farne, Staple, and Brownsman and the mainland. <bold>(B)</bold> Nitrogen input by seabirds for Inner Farne, Staple, Brownsman, and the mainland. Values are mean ( &#xb1; SE, n = 5). Box limits represent the first and third quartiles (25% and 75% percentiles), middle line represents the median (50% percentile), and whiskers represent the smallest and largest observations less than or equal to 1.5 &#xd7; inter-quartile range. Note differences in y-axes.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1343966-g002.tif"/>
</fig>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Nitrogen stable isotope and nutrient analysis</title>
<p>Soils, terrestrial plants, lichen, macroalgae, turf algae, barnacles, and limpets on seabird islands had enriched &#x3b4;<sup>15</sup>N in comparison to those from mainland sites (<xref ref-type="fig" rid="f3"><bold>Figures&#xa0;3A&#x2013;H</bold></xref>, <xref ref-type="table" rid="T1"><bold>Table&#xa0;1</bold></xref>). &#x3b4;<sup>15</sup>N was much higher for terrestrial samples soil, plants, and lichen than intertidal samples, reflecting the direct guano signature. Soil, terrestrial plants, lichen, turf algae, barnacles, and limpets also had significantly higher %N on high seabird density islands (<xref ref-type="fig" rid="f3"><bold>Figures&#xa0;3I&#x2013;P</bold></xref>, <xref ref-type="table" rid="T1"><bold>Table&#xa0;1</bold></xref>).</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Early-breeding season nitrogen isotope signals and total nitrogen content in terrestrial and intertidal areas on the low seabird density mainland sites and seabird islands. A-H, &#x3b4;<sup>15</sup>N values for <bold>(A)</bold> soil, <bold>(B)</bold> new growth leaves, <bold>(C)</bold> lichen, <bold>(D)</bold> macroalgae, <bold>(E)</bold> epiphytic algae, <bold>(F)</bold> turf algae, <bold>(G)</bold> barnacles, and <bold>(H)</bold> limpets. I-P, total nitrogen content for <bold>(I)</bold> soil, <bold>(J)</bold> new growth leaves, <bold>(K)</bold> lichen, <bold>(L)</bold> macroalgae, <bold>(M)</bold> epiphytic algae, <bold>(N)</bold> turf algae, <bold>(O)</bold> barnacles, and <bold>(P)</bold> limpets. Values are mean (&#xb1; SE, n = 15 for all sample types except for seabird island epiphytic algae, where n = 10). Box limits represent the first and third quartiles (25% and 75% percentiles), middle line represents the median (50% percentile), and whiskers represent the smallest and largest observations less than or equal to 1.5 &#xd7; inter-quartile range.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1343966-g003.tif"/>
</fig>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Early breeding season nitrogen isotope and total nitrogen statistical results in terrestrial and intertidal areas on low seabird density mainland sites and seabird islands.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" rowspan="2" align="left">Sample type</th>
<th valign="top" colspan="2" align="center">&#x3b4;<sup>15</sup>N isotopes</th>
<th valign="top" colspan="2" align="center">Total nitrogen content</th>
</tr>
<tr>
<th valign="middle" align="center">&#x3c7;<sup>2</sup>
</th>
<th valign="middle" align="center"><italic>p</italic>
</th>
<th valign="middle" align="center">&#x3c7;<sup>2</sup>
</th>
<th valign="middle" align="center"><italic>p</italic>
</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Soil</td>
<td valign="top" align="center">22.56</td>
<td valign="top" align="center"><bold>&lt; 0.001</bold>
</td>
<td valign="top" align="center">8.51</td>
<td valign="top" align="center"><bold>0.004</bold>
</td>
</tr>
<tr>
<td valign="top" align="left">Terrestrial plant</td>
<td valign="top" align="center">16.58</td>
<td valign="top" align="center"><bold>&lt; 0.001</bold>
</td>
<td valign="top" align="center">4.27</td>
<td valign="top" align="center"><bold>0.039</bold>
</td>
</tr>
<tr>
<td valign="top" align="left">Lichen</td>
<td valign="top" align="center">8.99</td>
<td valign="top" align="center"><bold>0.003</bold>
</td>
<td valign="top" align="center">11.77</td>
<td valign="top" align="center"><bold>&lt; 0.001</bold>
</td>
</tr>
<tr>
<td valign="top" align="left">Macroalgae</td>
<td valign="top" align="center">9.67</td>
<td valign="top" align="center"><bold>0.002</bold>
</td>
<td valign="top" align="center">1.66</td>
<td valign="top" align="center">0.197</td>
</tr>
<tr>
<td valign="top" align="left">Epiphytic algae</td>
<td valign="top" align="center">1.02</td>
<td valign="top" align="center">0.312</td>
<td valign="top" align="center">0.14</td>
<td valign="top" align="center">0.711</td>
</tr>
<tr>
<td valign="top" align="left">Turf algae</td>
<td valign="top" align="center">15.88</td>
<td valign="top" align="center"><bold>&lt; 0.001</bold>
</td>
<td valign="top" align="center">3.14</td>
<td valign="top" align="center">0.076</td>
</tr>
<tr>
<td valign="top" align="left">Barnacles</td>
<td valign="top" align="center">8.60</td>
<td valign="top" align="center"><bold>0.003</bold>
</td>
<td valign="top" align="center">15.35</td>
<td valign="top" align="center"><bold>&lt; 0.001</bold>
</td>
</tr>
<tr>
<td valign="top" align="left">Limpets</td>
<td valign="top" align="center">23.41</td>
<td valign="top" align="center"><bold>&lt; 0.001</bold>
</td>
<td valign="top" align="center">13.50</td>
<td valign="top" align="center"><bold>&lt; 0.001</bold>
</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>n = 15 for all sample types except for mainland epiphytic algae, where n = 10 (see Methods section 2.4). Bold p is significant (&lt; 0.05).</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="s3_3">
<label>3.3</label>
<title>Community composition analysis</title>
<sec id="s3_3_1">
<label>3.3.1</label>
<title>Algae and sessile filter feeders</title>
<p>There was no difference in intertidal community structure [all algal types, sessile filter feeders, and vacant surfaces (i.e., rock and sand)] between mainland sites and the seabird island sites when using the untransformed data (so as to determine the contributions of the dominant species; <xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4A</bold></xref>; PERMANOVA, <italic>F</italic> = 1.94, <italic>p</italic> = 0.122). There was weak evidence of a difference when analysing the presence-absence data to determine contributions of intermediate and rarer species (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4B</bold></xref>; PERMANOVA, <italic>F</italic> = 1.50, <italic>p</italic> = 0.058).</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Non-metric dimensional analysis (NMDS) biplots and boxplots comparing indicator species (those driving community dissimilarities) on low seabird density mainland sites and seabird islands. <bold>(A, B)</bold>, NMDS for <bold>(A)</bold> raw, untransformed algae and sessile filter feeder percent cover data (highlights dominant species importance), and <bold>(B)</bold> presence-absence transformed data (highlights intermediate and rarer species importance). Mainland sites are displayed as M1-4, seabird sites are displayed as S1-4. <bold>(C&#x2013;F)</bold>, algae and sessile filter feeder percent cover (untransformed data) for <bold>(C)</bold> serrated wrack, <bold>(D)</bold> Irish moss, <bold>(E)</bold> dulse, and <bold>(F)</bold> pink encrusting algae. Values are mean (&#xb1; SE, n = 20). Box limits represent the first and third quartiles (25% and 75% percentiles), middle line represents the median (50% percentile), and whiskers represent the smallest and largest observations less than or equal to 1.5 &#xd7; inter-quartile range.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1343966-g004.tif"/>
</fig>
<p>When analysing the contributions of the dominant species, indicator species analysis (IndVal) determined which species characterised each habitat. IndVal determined that whilst some species were driving differences between the seabird islands and the mainland, none of the differences were significant. These species were the macroalga <italic>Fucus serratus</italic> (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4C</bold></xref>; mean &#xb1; SE; Chi-squared test; mainland = 12.7 &#xb1; 5.9% per 0.5 m<sup>2</sup>, seabird islands = 2.6 &#xb1; 1.3% per 0.5 m<sup>2</sup>; <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>0.95, <italic>p</italic> = 0.329), and turf algae <italic>Chrondus crispus</italic> (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4D</bold></xref>; mean &#xb1; SE; not statistically analysed; mainland = 0.0 &#xb1; 0.0% per 0.5 m<sup>2</sup>, seabird islands = 0.5 &#xb1; 0.2% per 0.5 m<sup>2</sup>), <italic>Palmaria palmata</italic> (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4E</bold></xref>; mean &#xb1; SE; Chi-squared test; mainland = 0.4 &#xb1; 0.2% per 0.5 m<sup>2</sup>, seabird islands = 1.9 &#xb1; 0.8% per 0.5 m<sup>2</sup>; <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>0.73, <italic>p</italic> = 0.394), and <italic>Corallinaceae</italic> (<xref ref-type="fig" rid="f4"><bold>Figure&#xa0;4F</bold></xref>; mean &#xb1; SE; Chi-squared test; mainland = 3.5 &#xb1; 1.4% per 0.5 m<sup>2</sup>, seabird islands = 6.5 &#xb1; 1.7% per 0.5 m<sup>2</sup>; <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>1.14, <italic>p</italic> = 0.286).</p>
<p>For intermediate and rarer species, IndVal again determined the species that characterised the seabird island mainland habitats. Again, non-significant community differences were driven by the green turf algae <italic>Cladophora</italic> (proportion of samples where species was present &#xb1; SE; Chi-squared test; mainland = 0.3 &#xb1; 0.1%, seabird islands = 0.6 &#xb1; 0.1; <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>0.95, <italic>p</italic> = 0.329), and <italic>Mastocarpus stellatus</italic> (proportion of samples where species was present &#xb1; SE; Chi-squared test; mainland = 0.2 &#xb1; 0.1, seabird islands = 0.3 &#xb1; 0.1; <italic>df</italic> = 1, <italic>N</italic> = 15, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>0.43, <italic>p</italic> = 0.511).</p>
</sec>
<sec id="s3_3_2">
<label>3.3.2</label>
<title>Intertidal invertebrate fauna</title>
<p>Moderate evidence was found of a difference between mainland and seabird island intertidal invertebrate faunal communities when analysing presence-absence data to determine the contributions of intermediate and rarer species (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5B</bold></xref>; PERMANOVA, <italic>F</italic> = 2.64, <italic>p</italic> = 0.026). Dogwhelks drove dissimilarities between the mainland and the seabird island communities (proportion of samples where species was present &#xb1; SE; not statistically analysed; mainland = 0.0 &#xb1; 0.0, seabird islands = 0.4 &#xb1; 0.1).</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Non-metric dimensional analysis (NMDS) biplots and boxplots comparing indicator species (those driving community dissimilarities) on the low seabird density mainland sites and seabird islands. A-B, NMDS for <bold>(A)</bold> untransformed data (highlights dominant species importance), and <bold>(B)</bold> presence-absence transformed data (highlights intermediate and rarer species importance). Mainland sites are displayed as M1-4, seabird sites are displayed as S1-4. C, count (untransformed data) for <bold>(C)</bold> dogwhelks, <bold>(D)</bold> common winkles, and <bold>(E)</bold> beadlet anemones. Values are mean (&#xb1; SE, n = 20). Box limits represent the first and third quartiles (25% and 75% percentiles), middle line represents the median (50% percentile), and whiskers represent the smallest and largest observations less than or equal to 1.5 &#xd7; inter-quartile range.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1343966-g005.tif"/>
</fig>
<p>There was no difference in intertidal invertebrate fauna when analysing the untransformed data to determine the contributions of the dominant species (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5A</bold></xref>; PERMANOVA, <italic>F</italic> = 2.39, <italic>p</italic> = 0.113). The species that drove dissimilarities between the mainland and seabird island communities were dogwhelks (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5C</bold></xref>; mean &#xb1; SE; not statistically analysed; mainland = 0.0 &#xb1; 0.0 no. per 0.5 m<sup>2</sup>, seabird islands = 1.1 &#xb1; 0.4 no. per 0.5 m<sup>2</sup>) and common winkle (<xref ref-type="fig" rid="f5"><bold>Figure&#xa0;5D</bold></xref>; mean &#xb1; SE; Chi-squared test; mainland = 1.8 &#xb1; 1.0 no. per 0.5 m<sup>2</sup>, seabird islands = 0.4 &#xb1; 0.2 no. per 0.5 m<sup>2</sup>; <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>0.77, <italic>p</italic> = 0.381).</p>
</sec>
<sec id="s3_3_3">
<label>3.3.3</label>
<title>Ecological community functional groups</title>
<p>Strong evidence was found of higher intertidal predator densities on seabird islands than on the mainland (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure S1E</bold></xref>; Chi-squared test; mean &#xb1; SE; mainland = 0.1 &#xb1; 0.1 no. per 0.5 m<sup>2</sup>, seabird islands = 1.3 &#xb1; 0.4 no. per 0.5 m<sup>2</sup>; <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>7.22, <italic>p</italic> = 0.007). However, there were no significant differences in densities in macroalgae (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure S1A</bold></xref>; Chi-squared test; mean &#xb1; SE; mainland = 48.9 &#xb1; 9.1% per 0.5 m<sup>2</sup>, seabird islands = 63.0 &#xb1; 6.7% per 0.5 m<sup>2</sup>; <italic>df</italic> = 1, <italic>N</italic> = 15, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>1.38, <italic>p</italic> = 0.240), turf algae (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure S1B</bold></xref>; Chi-squared test (square root transformed); mean &#xb1; SE; mainland = 11.1 &#xb1; 3.2% per 0.5 m<sup>2</sup>, seabird islands = 14.3 &#xb1; 3.3% per 0.5 m<sup>2</sup>; <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>0.28, <italic>p</italic> = 0.597), sessile filter feeders (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure S1C</bold></xref>; Chi-squared test; mean &#xb1; SE; mainland = 28.1 &#xb1; 6.6% per 0.5 m<sup>2</sup>, seabird islands = 31.6 &#xb1; 4.6% per 0.5 m<sup>2</sup>; <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>0.06, <italic>p</italic> = 0.812), and other herbivores (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure S1D</bold></xref>; Chi-squared test; mean &#xb1; SE; mainland = 3.4 &#xb1; 1.3 no. per 0.5 m<sup>2</sup>, seabird islands = 1.1 &#xb1; 0.3 no. per 0.5 m<sup>2</sup>, <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>1.66, <italic>p</italic> = 0.198).</p>
</sec>
<sec id="s3_3_4">
<label>3.3.4</label>
<title>Diversity measures</title>
<p>There was weak evidence of a trend towards higher algal species richness on seabird islands compared to the mainland (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure S2A</bold></xref>; Chi-squared test; mean &#xb1; SE; mainland = 4.8 &#xb1; 0.5 no. per 0.5 m<sup>2</sup>, seabirds = 6.6 &#xb1; 0.5 no. per 0.5 m<sup>2</sup>; <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>3.23, <italic>p</italic> = 0.072). However, there was no difference in the total number of intertidal fauna species (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure S2B</bold></xref>; mean &#xb1; SE; Chi-squared test; mainland = 1.7 &#xb1; 0.2 no. per 0.5 m<sup>2</sup>, seabird islands = 2.0 &#xb1; 0.3 no. per 0.5 m<sup>2</sup>; <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>0.55, <italic>p</italic> = 0.458). No difference in Shannon&#x2019;s Diversity Index (H&#x2019;) between seabird islands and the mainland was detected for algal species (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure S2C</bold></xref>; mean &#xb1; SE; Chi-squared test; mainland = 0.6 &#xb1; 0.1 H&#x2019;, seabird islands = 0.8 &#xb1; 0.1 H&#x2019;; <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>0.41, <italic>p</italic> = 0.521) or fauna species (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure S2D</bold></xref>; mean &#xb1; SE; Chi-squared test; mainland = 0.6 &#xb1; 0.1 H&#x2019;, seabird islands = 0.6 &#xb1; 0.1 H&#x2019;; <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>0.19, <italic>p</italic> = 0.661). No distinct variation in Pielou&#x2019;s Evenness Index (J&#x2019;) between the seabird islands and the mainland was detected for algal species (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure S2E</bold></xref>; mean &#xb1; SE; Chi-squared test; mainland = 0.2 &#xb1; 0.0 J&#x2019;, seabird islands = 0.2 &#xb1; 0.0 J&#x2019;; <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>0.41, <italic>p</italic> = 0.521) or fauna species (<xref ref-type="supplementary-material" rid="SM1"><bold>Supplementary Figure S2F</bold></xref>; mean &#xb1; SE; Chi-squared test; mainland = 0.2 &#xb1; 0.0 J&#x2019;, seabird islands = 0.2 &#xb1; 0.0 J&#x2019;; <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>0.19, <italic>p</italic> = 0.662).</p>
</sec>
</sec>
<sec id="s3_4">
<label>3.4</label>
<title>Limpet size, abundance, and biomass</title>
<p>Due to high variation, there was no difference in limpet length despite a higher effect size on the seabird islands compared to the mainland (<xref ref-type="fig" rid="f6"><bold>Figure&#xa0;6A</bold></xref>; Chi-squared test; mean &#xb1; SE, mainland = 29.7 &#xb1; 0.6&#xa0;mm, seabird islands = 35.8 &#xb1; 0.4&#xa0;mm; <italic>df</italic> = 1, <italic>N</italic> = 400, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>2.17, <italic>p</italic> = 0.141). Similarly, there was no difference in limpet abundance despite a high effect size (<xref ref-type="fig" rid="f6"><bold>Figure&#xa0;6C</bold></xref>; Chi-squared test; mean &#xb1; SE; mainland = 12.7 &#xb1; 3.6 no. per 0.5 m<sup>2</sup>; seabird islands = 5.0 &#xb1; 1.1 no. per 0.5 m<sup>2</sup>; <italic>df</italic> = 1, <italic>N</italic> = 20, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>2.83, <italic>p</italic> = 0.160). No difference in limpet biomass was detected between the mainland and seabird islands (<xref ref-type="fig" rid="f6"><bold>Figure&#xa0;6B</bold></xref>; Chi-squared test; mean &#xb1; SE; mainland = 4.8 &#xb1; 1.4&#xa0;g per 0.5 m<sup>2</sup>, seabird islands = 3.1 &#xb1; 0.7&#xa0;g per 0.5 m<sup>2</sup>; <italic>df</italic> = 1, <italic>N</italic> = 400, <italic>&#x3c7;<sup>2</sup>&#xa0;=&#xa0;</italic>1.19, <italic>p</italic> = 0.275).</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Limpet distribution boxplots comparing low seabird density mainland sites and seabird islands. Distributions for <bold>(A)</bold> limpet length, <bold>(B)</bold> limpet abundance, and <bold>(C)</bold> limpet biomass. Values are mean (&#xb1; SE, n = 400 for A; n = 20 for <bold>(B, C)</bold>. Box limits represent the first and third quartiles (25% and 75% percentiles), middle line represents the median (50% percentile), and whiskers represent the smallest and largest observations less than or equal to 1.5 &#xd7; inter-quartile range. Note differences in y-axis.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1343966-g006.tif"/>
</fig>
</sec>
<sec id="s3_5">
<label>3.5</label>
<title>Recolonisation patterns</title>
<p>Recolonisation rates of barnacles following experimental clearances were higher on the mainland than on seabird islands, and this effect did not vary with time (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7A</bold>
</xref>; seabirds*time; Chi-squared test; df = 2, N = 15, c2 = 1.31, p = 0.520; seabirds; Chi- squared test; df = 2, N = 15, c2 = 4.38, p = 0.036). Barnacle abundance was significantly higher at 16 weeks post-clearance (Bonferroni; F = 10.240, p = 0.009), but there was only weak evidence of a difference in abundance in pre-clearance plots (Bonferroni; F = 5.97, p = 0.063) and little evidence at eight weeks post-clearance (Bonferroni; F = 4.86, p = 0.108).</p>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>Recolonisation rates of clearance for low seabird density mainland sites and seabird islands. Recolonisation rates from pre-clearance, zero weeks, eight weeks, and 16 weeks for mainland sites and seabird islands for <bold>(A)</bold> barnacles, and <bold>(B)</bold> limpets. Values are mean (&#xb1; SE; n = 5). Box limits represent the first and third quartiles (25% and 75% percentiles), middle line represents the median (50% percentile), and whiskers represent the smallest and largest observations less than or equal to 1.5 &#xd7; inter-quartile range.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-11-1343966-g007.tif"/>
</fig>
<p>Limpet recolonisation differed between mainland/seabirds and varied by time (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7B</bold>
</xref>; seabirds*time; Chi-squared test; df = 2, N = 15, c2 = 5.98, p = 0.050). Strong evidence was found of higher limpet abundance at seabird sites at eight weeks post-clearance (Bonferroni; F = 11.10, p = 0.006), but there was no difference in abundance in pre-clearance plots (Bonferroni; F = 2.38, p = 0.402) or at 16 weeks post-clearance (Bonferroni; F = 0.35, p = 1.000).</p>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<p>The findings of our study provide evidence that that seabirds are inputting pelagic nutrients to the terrestrial and intertidal ecosystems of temperate UK islands. Thus, seabirds provide a key cross-ecosystem link. However, there was limited evidence for seabird nutrients affecting intertidal community structure. There was some evidence that intertidal diversity was marginally greater on seabird islands (i.e., differences in intertidal fauna communities and functional group densities), while barnacle recolonisation rates were faster at low seabird density mainland sites. Disentangling the effects of seabird nutrients from other potential differences between the seabird islands and mainland locations will enable a better understanding of the role of seabird nutrient subsidies in temperate rocky intertidal systems.</p>
<p>Higher densities of seabirds were found on the Farne Islands, Northumberland, compared to the Northumberland mainland by four orders of magnitude. By foraging in pelagic feeding areas, seabirds were estimated to have vectored between three and four orders of magnitude more nitrogen inputs to the seabird islands compared to the mainland. The seabird density and nitrogen inputs revealed by our study are comparable with those observed in a tropical environment, where seabird density and nitrogen input were both three orders of magnitude greater on rat-free islands than those that had rats present on them (<xref ref-type="bibr" rid="B32">Graham et&#xa0;al., 2018</xref>). The higher density of seabirds on the Farne Islands is most likely due to these islands being predator free, in comparison to the mainland sites where predators including rats, cats and foxes would decimate seabird populations in accessible habitats.</p>
<p>The nutrient inputs from seabirds to the Farne Islands are used by terrestrial and intertidal organisms, evidenced by higher &#x3b4;<sup>15</sup>N signals being found in soil, new-growth leaves, and lichen, as well as macroalgae, turf algae, barnacles, and limpets on seabird islands compared to low density seabird sites. The third highest &#x3b4;<sup>15</sup>N values were found in limpets, which is likely explained by &#x3b4;<sup>15</sup>N values increasing with trophic levels, and microphagous limpets grazing on <sup>15</sup>N enriched early macroalgal stages and microalgal film (<xref ref-type="bibr" rid="B45">Jenkins and Hartnoll, 2001</xref>). This stepwise enrichment is consistent with other research in marine systems (<xref ref-type="bibr" rid="B38">Hobson and Welch, 1992</xref>; <xref ref-type="bibr" rid="B97">Zmudczy&#x144;ska-Skarbek et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B32">Graham et&#xa0;al., 2018</xref>), although other studies have found mixed responses (<xref ref-type="bibr" rid="B31">Gagnon et&#xa0;al., 2013</xref>), or instances of primary consumers having higher <sup>15</sup>N values than secondary consumers (<xref ref-type="bibr" rid="B40">Hong et&#xa0;al., 2019</xref>; <xref ref-type="bibr" rid="B3">Andrades et&#xa0;al., 2024</xref>). Significantly elevated nitrogen content was also found in soil, terrestrial plants, lichen, barnacles, and limpets, as well as marginally elevated turf algae nitrogen content. These findings support the hypothesis that seabird nutrient subsidies are cycling through the terrestrial and intertidal habitats. Previous studies in saltwater systems show that the proportion of seabird-derived nutrients decreases with increasing distance from seabird colonies (<xref ref-type="bibr" rid="B47">Jones, 2010</xref>; <xref ref-type="bibr" rid="B62">McCauley et&#xa0;al., 2012</xref>; <xref ref-type="bibr" rid="B57">Lorrain et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B81">Savage, 2019</xref>; <xref ref-type="bibr" rid="B10">Benkwitt et&#xa0;al., 2021a</xref>). Indeed, a similar pattern was generally observed in this study, with the samples closest to the seabird colonies (terrestrial plants and soil) having the greatest &#x3b4;<sup>15</sup>N values on seabird islands, and macroalgae, epiphytic algae, and turf algae that inhabit intertidal areas further from the colony having lower values. By contrast, total nitrogen content was generally higher within intertidal samples than terrestrial ones, with all the intertidal samples except macroalgae having a higher %N value. This may indicate that nitrogen is also being input directly from the marine environment, such as from the North Sea. Whilst the difference in the patterns of &#x3b4;<sup>15</sup>N and %N could be argued to contradict each other to an extent, the high &#x3b4;<sup>15</sup>N values show that, despite other sources of nutrient input to the system, pelagic-derived seabird nutrients make up a dominant component of the Farne Island terrestrial and intertidal nutrient budgets.</p>
<p>Compared with tropical or polar regions, cold temperate regions are typically very productive and may provide high levels of background nutrients to the region (<xref ref-type="bibr" rid="B84">Sigman and Hain, 2012</xref>). Furthermore, two possible sources of nitrogen input exist in two areas close to the mainland study sites. First, there are reported marine nitrate toxicity issues in the River Tweed, around eight to ten kilometres north of the mainland study sites (<xref ref-type="bibr" rid="B44">Jarvie et&#xa0;al., 2002</xref>). Second, there is eutrophication at Lindisfarne National Nature Reserve, around 25&#xa0;km north of the mainland study sites, where sources of nitrogen include diffuse agricultural runoff, misconnections in public and private sewage treatment works, and seasonal tourism pressures on dated sewage treatment facilities (<xref ref-type="bibr" rid="B69">NCC, 2020</xref>). Sewage &#x3b4;<sup>15</sup>N signatures are comparable to those of animal waste and are almost always over 10&#x2030;, reaching up to around 20&#x2030; (<xref ref-type="bibr" rid="B36">Heaton, 1986</xref>). In our study sites, mainland &#x3b4;<sup>15</sup>N signals were below 10&#x2030; in each sample type, providing evidence that sewage is not a major source of nitrogen input to these areas. Fertiliser &#x3b4;<sup>15</sup>N signatures typically range from -2 to 2&#x2030;, though can reach values of 4&#x2030; (<xref ref-type="bibr" rid="B6">Bateman and Kelly, 2007</xref>). A global study of soil &#x3b4;<sup>15</sup>N trends estimated that UK &#x3b4;<sup>15</sup>N values range between 3.5 and 6.2&#x2030; (<xref ref-type="bibr" rid="B1">Amundson et&#xa0;al., 2003</xref>), and indeed, on the mainland study sites, &#x3b4;<sup>15</sup>N signatures were 5.4 &#xb1; 0.3&#x2030;. Given that fertiliser &#x3b4;<sup>15</sup>N values are lower than those of both typical soil and those observed on the mainland, &#x3b4;<sup>15</sup>N cannot reliably be used to determine whether fertiliser is a source of nitrogen input to the mainland sites.</p>
<p>As well as being seabird breeding colonies, seals also inhabit the Farne Islands, predominantly grey seals (<italic>Halichoerus grypus</italic>), who use the islands as rookeries during their breeding season from October to January. Foraging distance in grey seals from the Farne Islands varies from short trips to known haul-out sites (mean distance 39.8&#xa0;km) to long and distant journeys (2100&#xa0;km; <xref ref-type="bibr" rid="B63">McConnell et&#xa0;al., 1999</xref>). Therefore, it is likely that seals are also providing an allochthonous source of nutrients to the Farne Islands. Furthermore, seal biomass exceeds that of seabird biomass on two of the three study islands (<xref ref-type="bibr" rid="B68">National Trust, 2022b</xref>). However, seabird abundance is approximately fifty times higher than seal abundance across the three islands (<xref ref-type="bibr" rid="B92">Webber et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B67">National Trust, 2022a</xref>, <xref ref-type="bibr" rid="B68">National Trust, 2022b</xref>). In terms of seal proximity to the sites on the seabird islands, it was confirmed by National Trust rangers working on the Farne Islands that none of the study sites were seal rookeries, and that the seabird colonies were closer to each study site than seal rookeries were. Some seabird species, like the gulls, also travel across the islands frequently, whereas the seals tend to haul out in one spot. Whilst the particularly high &#x3b4;<sup>15</sup>N in the terrestrial samples suggest that it is the seabirds are the dominant cause of increased nitrogen loads to the terrestrial environments, it is clear that the seals are also providing a key source of nutrient input to the Farne Islands.</p>
<p>Nutrient input from seabirds can increase primary producer biomass, resulting in heightened algal herbivory and herbivore biomass, leading to increased resources for intertidal predators (<xref ref-type="bibr" rid="B64">Methratta, 2004</xref>). The similar densities of macroalgae, turf algae, sessile filter feeders, and other herbivores observed in this study indicate that the intertidal communities on seabird islands and the mainland are functionally similar. However, dogwhelks comprised 81% of the studied predators on seabird islands, yet were absent from mainland sites. Moreover, all sites were relatively sheltered and fairly similar in terms of wave exposure (on the Ballantine Index all mainland sites were considered slightly more exposed than those on the seabird islands: mainland sites = 5; Farne Islands = 6, 7 and 7). Given that typically dogwhelks are more common on exposed rocky shores than sheltered coasts (<xref ref-type="bibr" rid="B5">Ballantine, 1961</xref>), dogwhelks would have been expected to be more common on the mainland sites, but here, we found the opposite to be true. Thus, another cause of higher densities of intertidal predators could be that seabird nutrient subsidies result in bottom-up interactions, possibly cascading productivity up the food web to create greater predator densities (<xref ref-type="bibr" rid="B95">Young et&#xa0;al., 2011</xref>). While it is possible that the increased seabird nutrients on the Farne Islands are supporting this additional trophic level, on its own this evidence falls short of substantiating this effect.</p>
<p>When comparing species diversity and taxonomic relatedness, we found that the intertidal communities surrounding seabird islands and mainland sites were similar and that there was no observed difference in fauna species richness, untransformed community species diversity, or evenness for either algal or fauna species. A non-significant trend towards higher algal species richness was observed, though, similar to higher macroalgal diversity on islands with the longest history of seabird presence in New Zealand (<xref ref-type="bibr" rid="B78">Rankin and Jones, 2021</xref>). Overall, though, the diversity effects of seabird nutrients in the present study were small. Drivers of ecological community diversity are highly complex and it is possible confounding variables exist between the seabird islands and the mainland in this study, as well as high levels of nutrients already being present in the cold-water system. In addition, whilst the exposure levels between the mainland and the seabird islands were relatively similar (Ballantine Index of mainland sites = 5; Farne Islands = 6 and 7), those differences are large enough for some different community patterns to exist (<xref ref-type="bibr" rid="B5">Ballantine, 1961</xref>). Island biogeography may also have influenced differences in community patterns (<xref ref-type="bibr" rid="B59">MacArthur and Wilson, 1967</xref>). However, island sites were less than five kilometres from the mainland and were close to one another as well as the other surrounding islands of the Farne Islands archipelago, thereby reducing the likelihood of this impact (<xref ref-type="fig" rid="f1"><bold>Figure&#xa0;1</bold></xref>). The effects of island biogeography are also more relevant to terrestrial species than marine or intertidal species, due to differences in their dispersal mechanisms, which formed the focus of this study (<xref ref-type="bibr" rid="B25">Carr et&#xa0;al., 2003</xref>; <xref ref-type="bibr" rid="B34">Hachich et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B88">Strathmann, 2015</xref>; <xref ref-type="bibr" rid="B75">Pinheiro et&#xa0;al., 2017</xref>).</p>
<p>Although no statistical difference was detected in limpet size, abundance or biomass between seabird islands and the mainland, the effect size between the mainland and seabird islands was high, with mean limpet size 20% larger on seabird islands. It is possible that greater sampling is needed to reduce variability and detect statistically significant differences. The relationship between body size and population density within communities is central to explaining community resource use (<xref ref-type="bibr" rid="B72">Pagel et&#xa0;al., 1991</xref>), and indeed, the quantity of algae grazed is more closely related to the biomass of limpets than their abundance (<xref ref-type="bibr" rid="B20">Branch, 1971</xref>; <xref ref-type="bibr" rid="B19">Bosman and Hockey, 1988</xref>). Shell growth in common limpets is rapid in their first two years, with growth rates between 0.4 and 2.0 mm/month. This rate then decreases over time, with a negative linear relationship existing between initial length and increment each year. By five years of age, they exhibit no measurable growth rate. Indeed, common limpet individuals that are five years and older are few, and form only negligible proportions of limpet populations (<xref ref-type="bibr" rid="B15">Blackmore, 1969</xref>). Therefore, limpets attaining a larger size may suggest that their growth rates are more rapid (<xref ref-type="bibr" rid="B21">Branch, 1974</xref>; <xref ref-type="bibr" rid="B4">Balaparameswara Rao, 1976</xref>). If indeed seabird nutrient input is increasing limpet growth rate on the Farne Islands, then this is consistent with previous findings of faster limpet growth in the presence of seabirds (<xref ref-type="bibr" rid="B19">Bosman and Hockey, 1988</xref>). Enhanced growth rates have been observed in other marine organisms such as corals and fishes due to the bottom-up effects of allochthonous nutrients (<xref ref-type="bibr" rid="B32">Graham et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B81">Savage, 2019</xref>; <xref ref-type="bibr" rid="B11">Benkwitt et&#xa0;al., 2021b</xref>, <xref ref-type="bibr" rid="B9">Benkwitt et&#xa0;al., 2023</xref>). Another possibility is that limpet growth rates is related to population density in both populations. The lower population densities and higher body size on the seabird islands, and vice versa on the mainland, could be explained by density-dependent growth rates in response to the limited resources of space and food. This pattern has been observed in the closely related species, <italic>Patella depressa</italic>, and is further supported by the similar biomass between the seabird islands and the mainland (<xref ref-type="bibr" rid="B16">Boaventura et&#xa0;al., 2003</xref>).</p>
<p>Contrary to predictions, recolonisation rates were faster on the mainland than seabird islands for both barnacles and limpets, which was inconsistent with other studies (<xref ref-type="bibr" rid="B52">Kraufvelin et&#xa0;al., 2006</xref>; <xref ref-type="bibr" rid="B51">Kraufvelin, 2007</xref>; <xref ref-type="bibr" rid="B60">Masterson et&#xa0;al., 2008</xref>). Whilst we sought to minimise the impact of confounding variables in the study, the mainland and seabird islands are inherently physically and biologically different, and it is likely that seabird nutrient input is not the dominant driver of demographic differences between the communities. For example, lower recolonisation rates are expected on islands with increasing isolation (<xref ref-type="bibr" rid="B59">MacArthur and Wilson, 1967</xref>), but all of our study islands were close to the mainland (&lt;5&#xa0;km) so it is unlikely that the Theory of Island Biogeography alone explains the higher rates of recolonisation observed on the mainland. Barnacle abundances increase with heightened wave exposure levels, with the most sheltered shores being largely absent of barnacles (<xref ref-type="bibr" rid="B5">Ballantine, 1961</xref>). Indeed, the three sites on the seabird islands were more sheltered from the prevailing UK west to south-westerly winds compared to the mainland. Furthermore, biological barnacle settlement cues will also impact community differences, these including the presence of algae species and diatoms, the abundance/fineness of detrital matter, microheterogeneity, food supply, and chlorophyll <italic>a</italic> concentration (<xref ref-type="bibr" rid="B54">Le Tourneux and Bourget, 1988</xref>; <xref ref-type="bibr" rid="B37">Hills and Thomason, 1998</xref>; <xref ref-type="bibr" rid="B24">Burrows et&#xa0;al., 2010</xref>). Faster limpet recolonisation was also observed on mainland sites than on seabird islands, with limpet abundance at week eight driving the difference. It is likely that limpets were more abundant on the recently disturbed plots as this disturbance drove increases in the abundance of microphagous limpets&#x2019; preferred food type within the plots: a microalgal film largely consisting of diatoms and cyanobacteria (<xref ref-type="bibr" rid="B45">Jenkins and Hartnoll, 2001</xref>). It is also possible that the observed increase in barnacles recolonisation was in part driven by the increased limpet recolonisation, as has been observed in other empirical research (<xref ref-type="bibr" rid="B8">Benedetti-Cecchi, 2000</xref>).</p>
<p>To conclude, strong evidence was found of seabird nutrients moving through UK terrestrial and intertidal ecosystems, through elevated &#x3b4;<sup>15</sup>N signals and total nitrogen content. This study suggests a clear link by seabirds from pelagic to terrestrial and intertidal habitats. Intertidal diversity differences between the Farne Islands and the low seabird abundance mainland were small, though higher numbers of intertidal dogwhelks and weak evidence of greater algal species richness did offer some support of community differences on the seabird islands. However, cold-water systems typically have high levels of background nutrients, and anthropogenic nutrient inputs from the mainland could not be quantified, so community differences between the islands and the mainland cannot definitively be attributed to seabird nutrients. This research shows the importance of seabirds in connecting ecosystems, but further work is necessary to disentangle the influence of seabird nutrient inputs compared to other factors that structure temperate rocky intertidal ecosystems.</p>
</sec>
<sec id="s5" sec-type="data-availability">
<title>Data availability statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s6" sec-type="ethics-statement">
<title>Ethics statement</title>
<p>The animal study was approved by Lancaster University Animal Welfare and Ethics Review Board. The study was conducted in accordance with the local legislation and institutional requirements.</p>
</sec>
<sec id="s7" sec-type="author-contributions">
<title>Author contributions</title>
<p>SH: Writing &#x2013; review &amp; editing, Writing &#x2013; original draft, Visualization, Software, Resources, Project administration, Methodology, Investigation, Formal analysis, Data curation, Conceptualization. CB: Writing &#x2013; review &amp; editing, Supervision, Resources, Methodology, Investigation, Funding acquisition. RD: Writing &#x2013; review &amp; editing, Visualization, Supervision, Resources, Methodology, Investigation, Funding acquisition. NG: Writing &#x2013; review &amp; editing, Validation, Supervision, Resources, Project administration, Methodology, Funding acquisition, Conceptualization.</p>
</sec>
</body>
<back>
<sec id="s8" sec-type="funding-information">
<title>Funding</title>
<p>The author(s) declare financial support was received for the research, authorship, and/or publication of this article. CB and RD were supported by the Bertarelli Foundation as part of the Bertarelli Programme in Marine Science (grant number: 2021-819047).</p>
</sec>
<ack>
<title>Acknowledgments</title>
<p>We thank Natural England, the National Trust, and Northumberland County Council for granting us permission to undertake the research. Note, initial agreement was made with Natural England for three sites on the Farne Islands for all methods, though partial observational data were also collected at a fourth site and communicated to Natural England for inclusion in the analyses. We also thank Dave Hughes for performing the stable isotope analysis.</p>
</ack>
<sec id="s9" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s10" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s11" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fmars.2024.1343966/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fmars.2024.1343966/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Presentation_1.pdf" id="SM1" mimetype="application/pdf"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Amundson</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Austin</surname> <given-names>A. T.</given-names>
</name>
<name>
<surname>Schuur</surname> <given-names>E. A. G.</given-names>
</name>
<name>
<surname>Yoo</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Matzek</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Kendall</surname> <given-names>C.</given-names>
</name>
<etal/>
</person-group>. (<year>2003</year>). <article-title>Global patterns of the isotopic composition of soil and plant nitrogen</article-title>. <source>Global Biogeochemical Cycles</source> <volume>17</volume>, <fpage>(1)</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1029/2002GB001903</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Anderson</surname> <given-names>M. J.</given-names>
</name>
<name>
<surname>Walsh</surname> <given-names>D. C. I.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>PERMANOVA, ANOSIM, and the Mantel test in the face of heterogeneous dispersions: What null hypothesis are you testing</article-title>? <source>Ecol. Monogr.</source> <volume>83</volume>, <fpage>557</fpage>&#x2013;<lpage>574</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1890/12-2010.1</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Andrades</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Cardozo-Ferreira</surname> <given-names>G. C.</given-names>
</name>
<name>
<surname>Benevides</surname> <given-names>L. J.</given-names>
</name>
<name>
<surname>Pimentel</surname> <given-names>C. R.</given-names>
</name>
<name>
<surname>Mancini</surname> <given-names>P. L.</given-names>
</name>
<name>
<surname>Ferreira</surname> <given-names>C. E.</given-names>
</name>
<etal/>
</person-group>. (<year>2024</year>). <article-title>Seabird guano reshapes intertidal reef food web in an isolated oceanic islet</article-title>. <source>Coral Reefs</source>. <volume>43</volume>, <fpage>1</fpage>&#x2013;<lpage>9</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00338-024-02469-w</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Balaparameswara Rao</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>1976</year>). <article-title>Studies on the growth of the limpet Cellana radiata (Born) (Gastropoda: Prosobranchia)</article-title>. <source>J. Molluscan Stud.</source> <volume>42</volume>, <fpage>136</fpage>&#x2013;<lpage>144</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/oxfordjournals.mollus.a065314</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ballantine</surname> <given-names>W. J.</given-names>
</name>
</person-group> (<year>1961</year>). <article-title>A biologically-defined exposure scale for the comparative description of rocky shores</article-title>. <source>Field Stud. J.</source> <volume>1</volume>, <fpage>1</fpage>&#x2013;<lpage>19</lpage>.</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bateman</surname> <given-names>A. S.</given-names>
</name>
<name>
<surname>Kelly</surname> <given-names>S. D.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Fertilizer nitrogen isotope signatures</article-title>. <source>Isotopes Environ. Health Stud.</source> <volume>43</volume>, <fpage>237</fpage>&#x2013;<lpage>247</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/10256010701550732</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bates</surname> <given-names>D.</given-names>
</name>
<name>
<surname>M&#xe4;chler</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Bolker</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Walker</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Fitting linear mixed-effects models using lme4</article-title>. <source>J. Stat. Software</source> <volume>67</volume>, <fpage>1</fpage>&#x2013;<lpage>48</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.18637/jss.v067.i01</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Benedetti-Cecchi</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Predicting direct and indirect interactions during succession in a mid-littoral rocky shore assemblage</article-title>. <source>Ecol. Monogr.</source> <volume>70</volume>, <fpage>45</fpage>&#x2013;<lpage>72</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1890/0012-9615(2000)070[0045:PDAIID]2.0.CO;2</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Benkwitt</surname> <given-names>C. E.</given-names>
</name>
<name>
<surname>D&#x2019;Angelo</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Gunn</surname> <given-names>R. L.</given-names>
</name>
<name>
<surname>Healing</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Mardones</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Wiedenmann</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2023</year>). <article-title>Seabirds boost coral reef resilience</article-title>. <source>Sci. Adv.</source> <volume>9</volume>, <elocation-id>eadj0390</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1126/sciadv.adj0390</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Benkwitt</surname> <given-names>C. E.</given-names>
</name>
<name>
<surname>Gunn</surname> <given-names>R. L.</given-names>
</name>
<name>
<surname>Le Corre</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Carr</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Graham</surname> <given-names>N. A. J.</given-names>
</name>
</person-group> (<year>2021</year>a). <article-title>Rat eradication restores nutrient subsidies from seabirds across terrestrial and marine ecosystems</article-title>. <source>Curr. Biol.</source> <volume>31</volume>, <fpage>2704</fpage>&#x2013;<lpage>2711.e2704</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.cub.2021.03.104</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Benkwitt</surname> <given-names>C. E.</given-names>
</name>
<name>
<surname>Taylor</surname> <given-names>B. M.</given-names>
</name>
<name>
<surname>Meekan</surname> <given-names>M. G.</given-names>
</name>
<name>
<surname>Graham</surname> <given-names>N. A. J.</given-names>
</name>
</person-group> (<year>2021</year>b). <article-title>Natural nutrient subsidies alter demographic rates in a functionally important coral-reef fish</article-title>. <source>Sci. Rep.</source> <volume>11</volume>, <fpage>12575</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41598-021-91884-y</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Benkwitt</surname> <given-names>C. E.</given-names>
</name>
<name>
<surname>Wilson</surname> <given-names>S. K.</given-names>
</name>
<name>
<surname>Graham</surname> <given-names>N. A. J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Seabird nutrient subsidies alter patterns of algal abundance and fish biomass on coral reefs following a bleaching event</article-title>. <source>Glob Chang Biol.</source> <volume>25</volume>, <fpage>2619</fpage>&#x2013;<lpage>2632</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/gcb.14643</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>BGS</collab>
</person-group>. (<year>1896</year>). <source>Quarter-inch geological map of England and Wales [Old Series]. Northumberland, and parts of Cumberland and Durham</source>. Available at: <uri xlink:href="https://webapps.bgs.ac.uk/data/maps/maps.cfc?method=viewRecord&amp;mapId=11707">https://webapps.bgs.ac.uk/data/maps/maps.cfc?method=viewRecord&amp;mapId=11707</uri> (Accessed <access-date>February 2023</access-date>).</citation>
</ref>
<ref id="B14">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Bibby</surname> <given-names>C. J.</given-names>
</name>
</person-group> (<year>2000</year>). <source>Bird census techniques</source> (<publisher-loc>London, San Diego</publisher-loc>: <publisher-name>Academic</publisher-name>).</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Blackmore</surname> <given-names>D. T.</given-names>
</name>
</person-group> (<year>1969</year>). <article-title>Studies of Patella vulgata L. I. Growth, reproduction and zonal distribution</article-title>. <source>J. Exp. Mar. Biol. Ecol.</source> <volume>3</volume>, <fpage>200</fpage>&#x2013;<lpage>213</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/0022-0981(69)90018-5</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Boaventura</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Da Fonseca</surname> <given-names>L. C.</given-names>
</name>
<name>
<surname>Hawkins</surname> <given-names>S. J.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Size matters: competition within populations of the limpet Patella depressa</article-title>. <source>J. Anim. Ecol.</source> <volume>72</volume>, <fpage>435</fpage>&#x2013;<lpage>446</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1046/j.1365-2656.2003.00713.x</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bosman</surname> <given-names>A. L.</given-names>
</name>
<name>
<surname>Du Toit</surname> <given-names>J. T.</given-names>
</name>
<name>
<surname>Hockey</surname> <given-names>P. A. R.</given-names>
</name>
<name>
<surname>Branch</surname> <given-names>G. M.</given-names>
</name>
</person-group> (<year>1986</year>). <article-title>A field experiment demonstrating the influence of seabird guano on intertidal primary production</article-title>. <source>Estuarine Coasta Shelf Sci.</source> <volume>23</volume>, <fpage>283</fpage>&#x2013;<lpage>294</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/0272-7714(86)90028-4</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bosman</surname> <given-names>A. L.</given-names>
</name>
<name>
<surname>Hockey</surname> <given-names>P. A. R.</given-names>
</name>
</person-group> (<year>1986</year>). <article-title>Seabird guano as a determinant of rocky intertidal community structure</article-title>. <source>Mar. Ecol. Prog. Ser.</source> <volume>32</volume>, <fpage>247</fpage>&#x2013;<lpage>257</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3354/meps032247</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bosman</surname> <given-names>A. L.</given-names>
</name>
<name>
<surname>Hockey</surname> <given-names>P. A. R.</given-names>
</name>
</person-group> (<year>1988</year>). <article-title>The influence of seabird guano on the biological structure of rocky intertidal communities on islands off the west coast of southern Africa</article-title>. <source>South Afr. J. Mar. Sci.</source> <volume>7</volume>, <fpage>61</fpage>&#x2013;<lpage>68</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2989/025776188784379134</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Branch</surname> <given-names>G. M.</given-names>
</name>
</person-group> (<year>1971</year>). <article-title>The ecology of patella linnaeus from the Cape Peninsula, South Africa I. Zonation, movements and feeding</article-title>. <source>Zoologica Africana</source> <volume>6</volume>, <fpage>1</fpage>&#x2013;<lpage>38</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/00445096.1971.11447402</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Branch</surname> <given-names>G. M.</given-names>
</name>
</person-group> (<year>1974</year>). <article-title>The ecology of <italic>patella</italic> linneaus from the Cape Peninsula, South Africa, 3. Growth rates</article-title>. <source>Trans. R. Soc. South Afr.</source> <volume>41</volume>, <fpage>161</fpage>&#x2013;<lpage>193</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/00359197409520069</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>BTO</collab>
</person-group>. (<year>2022</year>a). <source>Wetland Bird Survey (WeBS)</source>. Available at: <uri xlink:href="https://www.bto.org/our-science/projects/wetland-bird-survey">https://www.bto.org/our-science/projects/wetland-bird-survey</uri> (Accessed <access-date>July 2021</access-date>).</citation>
</ref>
<ref id="B23">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>BTO</collab>
</person-group>. (<year>2022</year>b). <source>Wetland Bird Survey (WeBS) Sites</source>. Available at: <uri xlink:href="https://app.bto.org/websonline/sites/data/sites-data.jsp">https://app.bto.org/websonline/sites/data/sites-data.jsp</uri> (Accessed <access-date>September 2022</access-date>).</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Burrows</surname> <given-names>M. T.</given-names>
</name>
<name>
<surname>Jenkins</surname> <given-names>S. R.</given-names>
</name>
<name>
<surname>Robb</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Harvey</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Spatial variation in size and density of adult and post-settlement Semibalanus balanoides: effects of oceanographic and local conditions</article-title>. <source>Mar. Ecol. Prog. Ser.</source> <volume>398</volume>, <fpage>207</fpage>&#x2013;<lpage>219</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3354/meps08340</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Carr</surname> <given-names>M. H.</given-names>
</name>
<name>
<surname>Neigel</surname> <given-names>J. E.</given-names>
</name>
<name>
<surname>Estes</surname> <given-names>J. A.</given-names>
</name>
<name>
<surname>Andelman</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Warner</surname> <given-names>R. R.</given-names>
</name>
<name>
<surname>Largier</surname> <given-names>J. L.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Comparing marine and terrestrial ecosystems: implications for the design of coastal marine reserves</article-title>. <source>Ecol. Appl.</source> <volume>13</volume>, <fpage>90</fpage>&#x2013;<lpage>107</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1890/1051-0761(2003)013[0090:CMATEI]2.0.CO;2</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Clarke</surname> <given-names>K. R.</given-names>
</name>
<name>
<surname>Gorley</surname> <given-names>R. N.</given-names>
</name>
<name>
<surname>Somerfield</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Warwick</surname> <given-names>R. M.</given-names>
</name>
</person-group> (<year>2014</year>). <source>Change in marine communities: an approach to statistical analysis and interpretation</source>. <publisher-loc>Plymouth, UK</publisher-loc>: <publisher-name>Primer-E Ltd</publisher-name>.</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Croxall</surname> <given-names>J. P.</given-names>
</name>
<name>
<surname>Butchart</surname> <given-names>S. H. M.</given-names>
</name>
<name>
<surname>Lascelles</surname> <given-names>B. E. N.</given-names>
</name>
<name>
<surname>Stattersfield</surname> <given-names>A. J.</given-names>
</name>
<name>
<surname>Sullivan</surname> <given-names>B. E. N.</given-names>
</name>
<name>
<surname>Symes</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group>. (<year>2012</year>). <article-title>Seabird conservation status, threats and priority actions: a global assessment</article-title>. <source>Bird Conserv. Int.</source> <volume>22</volume>, <fpage>1</fpage>&#x2013;<lpage>34</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1017/S0959270912000020</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>De Caceres</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Legendre</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Associations between species and groups of sites: indices and statistical inference</article-title>. <source>Ecology</source> <volume>90</volume>, <fpage>3566</fpage>&#x2013;<lpage>3574</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1890/08-1823.1</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dufr&#xea;ne</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Legendre</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>1997</year>). <article-title>Species assemblages and indicator species: the need for a flexible asymmetrical approach</article-title>. <source>Ecol. Monogr.</source> <volume>67</volume>, <fpage>345</fpage>&#x2013;<lpage>366</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2307/2963459</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>Esri</collab>
</person-group>. (<year>2022</year>). <source>ArcGIS Pro</source>. Available at: <uri xlink:href="https://www.esri.com/en-us/home">https://www.esri.com/en-us/home</uri> (Accessed <access-date>August 2021</access-date>).</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gagnon</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Rothausler</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Syrjanen</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Yli-Renko</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Jormalainen</surname> <given-names>V.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Seabird guano fertilizes Baltic Sea littoral food webs</article-title>. <source>PLoS One</source> <volume>8</volume>, <elocation-id>e61284</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal.pone.0061284</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Graham</surname> <given-names>N. A. J.</given-names>
</name>
<name>
<surname>Wilson</surname> <given-names>S. K.</given-names>
</name>
<name>
<surname>Carr</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Hoey</surname> <given-names>A. S.</given-names>
</name>
<name>
<surname>Jennings</surname> <given-names>S.</given-names>
</name>
<name>
<surname>MacNeil</surname> <given-names>M. A.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Seabirds enhance coral reef productivity and functioning in the absence of invasive rats</article-title>. <source>Nature</source> <volume>559</volume>, <fpage>250</fpage>&#x2013;<lpage>253</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41586-018-0202-3</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Grant</surname> <given-names>M. L.</given-names>
</name>
<name>
<surname>Bond</surname> <given-names>A. L.</given-names>
</name>
<name>
<surname>Lavers</surname> <given-names>J. L.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>The influence of seabirds on their breeding, roosting and nesting grounds: A systematic review and meta-analysis</article-title>. <source>J. Anim. Ecol.</source> <volume>91</volume>, <fpage>1266</fpage>&#x2013;<lpage>1289</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/1365-2656.13699</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hachich</surname> <given-names>N. F.</given-names>
</name>
<name>
<surname>Bonsall</surname> <given-names>M. B.</given-names>
</name>
<name>
<surname>Arraut</surname> <given-names>E. M.</given-names>
</name>
<name>
<surname>Barneche</surname> <given-names>D. R.</given-names>
</name>
<name>
<surname>Lewinsohn</surname> <given-names>T. M.</given-names>
</name>
<name>
<surname>Floeter</surname> <given-names>S. R.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Island biogeography: patterns of marine shallow-water organisms in the Atlantic Ocean</article-title>. <source>J. Biogeography</source> <volume>42</volume>, <fpage>1871</fpage>&#x2013;<lpage>1882</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/jbi.12560</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hawkins</surname> <given-names>S. J.</given-names>
</name>
<name>
<surname>Hartnoll</surname> <given-names>R. G.</given-names>
</name>
</person-group> (<year>1982</year>). <article-title>Settlement patterns of Semibalanus balanoides (L.) in the Isle of Man, (1977&#x2013;1981)</article-title>. <source>J. Exp. Mar. Biol. Ecol.</source> <volume>62</volume>, <fpage>271</fpage>&#x2013;<lpage>283</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/0022-0981(82)90207-6</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Heaton</surname> <given-names>T. H. E.</given-names>
</name>
</person-group> (<year>1986</year>). <article-title>Isotopic studies of nitrogen pollution in the hydrosphere and atmosphere: A review</article-title>. <source>Chem. Geology: Isotope Geosci. section</source> <volume>59</volume>, <fpage>87</fpage>&#x2013;<lpage>102</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/0168-9622(86)90059-X</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hills</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Thomason</surname> <given-names>J. C.</given-names>
</name>
</person-group> (<year>1998</year>). <article-title>The effect of scales of surface roughness on the settlement of barnacle (<italic>Semibalanus balanoides</italic>) cyprids</article-title>. <source>Biofouling</source> <volume>12</volume>, <fpage>57</fpage>&#x2013;<lpage>69</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/08927019809378346</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hobson</surname> <given-names>K. A.</given-names>
</name>
<name>
<surname>Welch</surname> <given-names>H. E.</given-names>
</name>
</person-group> (<year>1992</year>). <article-title>Determination of trophic relationships within a high Arctic marine food web using &#x3b4;13C and &#x3b4;15N analysis</article-title>. <source>Mar. Ecol. Prog. Ser.</source> <volume>84</volume>, <fpage>9</fpage>&#x2013;<lpage>18</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3354/meps084009</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Holm</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>1979</year>). <article-title>A simple sequentially rejective multiple test procedure</article-title>. <source>Scandinavian J. Stat</source>, <fpage>65</fpage>&#x2013;<lpage>70</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2307/4615733</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hong</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Gu</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Lee</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>The fates of nitrogen in an experimental wetland food web: a stable isotope study</article-title>. <source>Wetlands</source> <volume>39</volume>, <fpage>303</fpage>&#x2013;<lpage>310</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s13157-018-1085-7</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Hutchinson</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>1950</year>). <source>The biogeochemistry of vetebrate excretion</source> (<publisher-loc>New York, USA</publisher-loc>: <publisher-name>American Museum of Natual History, New York</publisher-name>).</citation>
</ref>
<ref id="B42">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>IPCC</collab>
</person-group>. (<year>2022</year>). <article-title>Climate Change 2022: Impacts, Adaptation and Vulnerability</article-title>. <source>Contribution of Working Group II to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change</source> [<person-group person-group-type="editor">
<name>
<surname>P&#xf6;rtner</surname> <given-names>H.-O.</given-names>
</name>
<name>
<surname>Roberts</surname> <given-names>D.C.</given-names>
</name>
<name>
<surname>Tignor</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Poloczanska</surname> <given-names>E. S.</given-names>
</name>
<name>
<surname>Mintenbeck</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Alegr&#xed;a</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group> (eds.)]. <publisher-loc>Cambridge, UK and New York, NY, USA</publisher-loc>: <publisher-name>Cambridge University Press</publisher-name>, <fpage>3056</fpage> pp. doi:&#xa0;<pub-id pub-id-type="doi">10.1017/9781009325844</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>IUCN</collab>
</person-group>. (<year>2022</year>). <source>2021 IUCN Red List of Threatened Species</source>. Available at: <uri xlink:href="https://www.iucnredlist.org/">https://www.iucnredlist.org/</uri> (Accessed <access-date>30 November 2021</access-date>).</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jarvie</surname> <given-names>H. P.</given-names>
</name>
<name>
<surname>Wade</surname> <given-names>A. J.</given-names>
</name>
<name>
<surname>Butterfield</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Whitehead</surname> <given-names>P. G.</given-names>
</name>
<name>
<surname>Tindall</surname> <given-names>C. I.</given-names>
</name>
<name>
<surname>Virtue</surname> <given-names>W. A.</given-names>
</name>
<etal/>
</person-group>. (<year>2002</year>). <article-title>Modelling nitrogen dynamics and distributions in the River Tweed, Scotland: an application of the INCA model</article-title>. <source>Hydrol. Earth Syst. Sci.</source> <volume>6</volume>, <fpage>433</fpage>&#x2013;<lpage>454</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.5194/hess-6-433-2002</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jenkins</surname> <given-names>S. R.</given-names>
</name>
<name>
<surname>Hartnoll</surname> <given-names>R. G.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Food supply, grazing activity and growth rate in the limpet Patella vulgata L.: a comparison between exposed and sheltered shores</article-title>. <source>J. Exp. Mar. Biol. Ecol.</source> <volume>258</volume>, <fpage>123</fpage>&#x2013;<lpage>139</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/S0022-0981(01)00211-8</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>JNCC</collab>
</person-group> (<year>2023</year>). <source>Seabirds Count</source> (<publisher-loc>Peterborough</publisher-loc>: <publisher-name>JNCC</publisher-name>). Available at: <uri xlink:href="https://jncc.gov.uk/our-work/seabirds-count/">https://jncc.gov.uk/our-work/seabirds-count/</uri> (Accessed <access-date>February 2024</access-date>).</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jones</surname> <given-names>H. P.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Prognosis for ecosystem recovery following rodent eradication and seabird restoration in an island archipelago</article-title>. <source>Ecol. Appl.</source> <volume>20</volume>, <fpage>1204</fpage>&#x2013;<lpage>1216</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1890/09-1172.1</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="web">
<person-group person-group-type="author">
<name>
<surname>Kassambara</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2021</year>). <source>rstatix: Pipe-Friendly Framework for Basic Statistical Tests.</source> (<publisher-name>R package version 0.7.0</publisher-name>). Available at: <uri xlink:href="https://CRAN.R-project.org/package=rstatix">https://CRAN.R-project.org/package=rstatix</uri> (Accessed <access-date>June 2021</access-date>).</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kazama</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Bottom-up effects on coastal marine ecosystems due to nitrogen input from seabird feces</article-title>. <source>Ornithol Sci.</source> <volume>18</volume>, <fpage>117</fpage>&#x2013;<lpage>126</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2326/osj.18.126</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kolb</surname> <given-names>G. S.</given-names>
</name>
<name>
<surname>Ekholm</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Hamb&#xe4;ck</surname> <given-names>P. A.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Effects of seabird nesting colonies on algae and aquatic invertebrates in coastal waters</article-title>. <source>Mar. Ecol. Prog. Ser.</source> <volume>417</volume>, <fpage>287</fpage>&#x2013;<lpage>300</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3354/meps08791</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kraufvelin</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Responses to nutrient enrichment, wave action and disturbance in rocky shore communities</article-title>. <source>Aquat. Bot.</source> <volume>87</volume>, <fpage>262</fpage>&#x2013;<lpage>274</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.aquabot.2007.06.011</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kraufvelin</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Moy</surname> <given-names>F. E.</given-names>
</name>
<name>
<surname>Christie</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Bokn</surname> <given-names>T. L.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Nutrient addition to experimental rocky shore communities revisited: Delayed responses, rapid recovery</article-title>. <source>Ecosystems</source> <volume>9</volume>, <fpage>1076</fpage>&#x2013;<lpage>1093</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10021-005-0188-1</pub-id>
</citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kruskal</surname> <given-names>J. B.</given-names>
</name>
</person-group> (<year>1964</year>). <article-title>Multidimensional scaling by optimizing goodness of fit to a nonmetric hypothesis</article-title>. <source>Psychometrika</source> <volume>29</volume>, <fpage>1</fpage>&#x2013;<lpage>27</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/BF02289565</pub-id>
</citation>
</ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Le Tourneux</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Bourget</surname> <given-names>E.</given-names>
</name>
</person-group> (<year>1988</year>). <article-title>Importance of physical and biological settlement cues used at different spatial scales by the larvae of Semibalanus balanoides</article-title>. <source>Mar. Biol.</source> <volume>97</volume>, <fpage>57</fpage>&#x2013;<lpage>66</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/BF00391245</pub-id>
</citation>
</ref>
<ref id="B55">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lindeboom</surname> <given-names>H. J.</given-names>
</name>
</person-group> (<year>1984</year>). <article-title>Nitrogen pathway in a penguin rookery</article-title>. <source>Ecology</source> <volume>65</volume>, <fpage>269</fpage>&#x2013;<lpage>277</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2307/1939479</pub-id>
</citation>
</ref>
<ref id="B56">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Loreau</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Holt</surname> <given-names>R. D.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Spatial flows and the regulation of ecosystems</article-title>. <source>Am. Nat.</source> <volume>163</volume>, <fpage>606</fpage>&#x2013;<lpage>615</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1086/382600</pub-id>
</citation>
</ref>
<ref id="B57">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lorrain</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Houlbreque</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Benzoni</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Barjon</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Tremblay-Boyer</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Menkes</surname> <given-names>C.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>). <article-title>Seabirds supply nitrogen to reef-building corals on remote Pacific islets</article-title>. <source>Sci. Rep.</source> <volume>7</volume>, <fpage>3721</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41598-017-03781-y</pub-id>
</citation>
</ref>
<ref id="B58">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lundberg</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Moberg</surname> <given-names>F.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Mobile link organisms and ecosystem functioning: Implications for ecosystem resilience and management</article-title>. <source>Ecosystems</source> <volume>6</volume>, <fpage>87</fpage>&#x2013;<lpage>98</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10021-002-0150-4</pub-id>
</citation>
</ref>
<ref id="B59">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>MacArthur</surname> <given-names>R. H.</given-names>
</name>
<name>
<surname>Wilson</surname> <given-names>E. O.</given-names>
</name>
</person-group> (<year>1967</year>). <source>The theory of island biogeography</source> (<publisher-loc>Princeton, New Jersey, United States</publisher-loc>: <publisher-name>Princeton university press</publisher-name>).</citation>
</ref>
<ref id="B60">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Masterson</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Arenas</surname> <given-names>F. A.</given-names>
</name>
<name>
<surname>Thompson</surname> <given-names>R. C.</given-names>
</name>
<name>
<surname>Jenkins</surname> <given-names>S. R.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Interaction of top down and bottom up factors in intertidal rockpools: Effects on early successional macroalgal community composition, abundance and productivity</article-title>. <source>J. Exp. Mar. Biol. Ecol.</source> <volume>363</volume>, <fpage>12</fpage>&#x2013;<lpage>20</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jembe.2008.06.001</pub-id>
</citation>
</ref>
<ref id="B61">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McCann</surname> <given-names>K. S.</given-names>
</name>
<name>
<surname>Rasmussen</surname> <given-names>J. B.</given-names>
</name>
<name>
<surname>Umbanhowar</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>The dynamics of spatially coupled food webs</article-title>. <source>Ecol. Lett.</source> <volume>8</volume>, <fpage>513</fpage>&#x2013;<lpage>523</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1461-0248.2005.00742.x</pub-id>
</citation>
</ref>
<ref id="B62">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McCauley</surname> <given-names>D. J.</given-names>
</name>
<name>
<surname>Desalles</surname> <given-names>P. A.</given-names>
</name>
<name>
<surname>Young</surname> <given-names>H. S.</given-names>
</name>
<name>
<surname>Dunbar</surname> <given-names>R. B.</given-names>
</name>
<name>
<surname>Dirzo</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Mills</surname> <given-names>M. M.</given-names>
</name>
<etal/>
</person-group>. (<year>2012</year>). <article-title>From wing to wing: the persistence of long ecological interaction chains in less-disturbed ecosystems</article-title>. <source>Sci. Rep.</source> <volume>2</volume>, <elocation-id>409</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/srep00409</pub-id>
</citation>
</ref>
<ref id="B63">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McConnell</surname> <given-names>B. J.</given-names>
</name>
<name>
<surname>Fedak</surname> <given-names>M. A.</given-names>
</name>
<name>
<surname>Lovell</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Hammond</surname> <given-names>P. S.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Movements and foraging areas of grey seals in the North Sea</article-title>. <source>J. Appl. Ecol.</source> <volume>36</volume>, <fpage>573</fpage>&#x2013;<lpage>590</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1046/j.1365-2664.1999.00429.x</pub-id>
</citation>
</ref>
<ref id="B64">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Methratta</surname> <given-names>E. T.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Top-down and bottom-up factors in tidepool communities</article-title>. <source>J. Exp. Mar. Biol. Ecol.</source> <volume>299</volume>, <fpage>77</fpage>&#x2013;<lpage>96</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jembe.2003.09.004</pub-id>
</citation>
</ref>
<ref id="B65">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Muff</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Nilsen</surname> <given-names>E. B.</given-names>
</name>
<name>
<surname>O&#x2019;Hara</surname> <given-names>R. B.</given-names>
</name>
<name>
<surname>Nater</surname> <given-names>C. R.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Rewriting results sections in the language of evidence</article-title>. <source>Trends Ecol. Evol.</source> <volume>37</volume>, <fpage>203</fpage>&#x2013;<lpage>210</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.tree.2021.10.009</pub-id>
</citation>
</ref>
<ref id="B66">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Mulder</surname> <given-names>C. P. H.</given-names>
</name>
<name>
<surname>Anderson</surname> <given-names>W. B.</given-names>
</name>
<name>
<surname>Towns</surname> <given-names>D. R.</given-names>
</name>
<name>
<surname>Bellingham</surname> <given-names>P. J.</given-names>
</name>
</person-group> (<year>2011</year>). &#x201c;<article-title>An introduction to seabird islands</article-title>,&#x201d; in <source>Seabird Islands: Ecology, Invasion, and Restoration.</source> (<publisher-name>Oxford University Press</publisher-name>, <publisher-loc>New York, USA</publisher-loc>), <fpage>3</fpage>&#x2013;<lpage>26</lpage>.</citation>
</ref>
<ref id="B67">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>National Trust</collab>
</person-group>. (<year>2022</year>a). <source>Farne Island Seabird Data.</source>, Unpublished data.</citation>
</ref>
<ref id="B68">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>National Trust</collab>
</person-group>. (<year>2022</year>b). <source>Farne Island Seal Data.</source>, Unpublished data.</citation>
</ref>
<ref id="B69">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>NCC</collab>
</person-group>. (<year>2020</year>). <source>Northumberland Coast AONB Management Plan</source>. Available at: <uri xlink:href="https://www.northumberlandcoastaonb.org/management-plan/">https://www.northumberlandcoastaonb.org/management-plan/</uri> (Accessed <access-date>February 2023</access-date>).</citation>
</ref>
<ref id="B70">
<citation citation-type="web">
<person-group person-group-type="author">
<name>
<surname>Oksanen</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Guillaume Blanchet.</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Friendly</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Kindt</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Legendre</surname> <given-names>P.</given-names>
</name>
<name>
<surname>McGlinn</surname> <given-names>D.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <source>vegan: Community Ecology Package.</source> (<publisher-name>R package version 2.5-7</publisher-name>). Available at: <uri xlink:href="https://CRAN.R-project.org/package=vegan">https://CRAN.R-project.org/package=vegan</uri> (Accessed <access-date>June 2021</access-date>).</citation>
</ref>
<ref id="B71">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Otero</surname> <given-names>X. L.</given-names>
</name>
<name>
<surname>La Pe&#xf1;a-Lastra</surname> <given-names>D.</given-names>
</name>
<name>
<surname>P&#xe9;rez-Alberti</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Ferreira</surname> <given-names>T. O.</given-names>
</name>
<name>
<surname>Huerta-Diaz</surname> <given-names>M. A.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Seabird colonies as important global drivers in the nitrogen and phosphorus cycles</article-title>. <source>Nat. Commun.</source> <volume>9</volume>, <fpage>1</fpage>&#x2013;<lpage>8</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41467-017-02446-8</pub-id>
</citation>
</ref>
<ref id="B72">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pagel</surname> <given-names>M. D.</given-names>
</name>
<name>
<surname>Harvey</surname> <given-names>P. H.</given-names>
</name>
<name>
<surname>Godfray</surname> <given-names>H. C. J.</given-names>
</name>
</person-group> (<year>1991</year>). <article-title>Species-abundance, biomass, and resource-use distributions</article-title>. <source>Am. Nat.</source> <volume>138</volume>, <fpage>836</fpage>&#x2013;<lpage>850</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1086/285255</pub-id>
</citation>
</ref>
<ref id="B73">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Paradis</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Schliep</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R</article-title>. <source>Bioinformatics</source> <volume>35</volume>, <fpage>526</fpage>&#x2013;<lpage>528</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/bioinformatics/bty633</pub-id>
</citation>
</ref>
<ref id="B74">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pielou</surname> <given-names>E. C.</given-names>
</name>
</person-group> (<year>1966</year>). <article-title>The measurement of diversity in different types of biological collections</article-title>. <source>J. Theor. Biol.</source> <volume>13</volume>, <fpage>131</fpage>&#x2013;<lpage>144</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/0022-5193(66)90013-0</pub-id>
</citation>
</ref>
<ref id="B75">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pinheiro</surname> <given-names>H. T.</given-names>
</name>
<name>
<surname>Bernardi</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Simon</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Joyeux</surname> <given-names>J.-C.</given-names>
</name>
<name>
<surname>Macieira</surname> <given-names>R. M.</given-names>
</name>
<name>
<surname>Gasparini</surname> <given-names>J. L.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>). <article-title>Island biogeography of marine organisms</article-title>. <source>Nature</source> <volume>549</volume>, <fpage>82</fpage>&#x2013;<lpage>85</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nature23680</pub-id>
</citation>
</ref>
<ref id="B76">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Polis</surname> <given-names>G. A.</given-names>
</name>
<name>
<surname>Anderson</surname> <given-names>W. B.</given-names>
</name>
<name>
<surname>Holt</surname> <given-names>R. D.</given-names>
</name>
</person-group> (<year>1997</year>). <article-title>Toward an integration of landscape and food web ecology: The dynamics of spatially subsidized food webs</article-title>. <source>Annu. Rev. Ecol. Systematics</source> <volume>28</volume>, <fpage>289</fpage>&#x2013;<lpage>316</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1146/annurev.ecolsys.28.1.289</pub-id>
</citation>
</ref>
<ref id="B77">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Portnoy</surname> <given-names>J. W.</given-names>
</name>
<name>
<surname>Soukup</surname> <given-names>M. A.</given-names>
</name>
</person-group> (<year>1990</year>). <article-title>Gull contributions of phosphorus and nitrogen to a Cape Cod kettle pond</article-title>. <source>Hydrobiologia</source> <volume>202</volume>, <fpage>61</fpage>&#x2013;<lpage>69</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/BF02208127</pub-id>
</citation>
</ref>
<ref id="B78">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rankin</surname> <given-names>L. L.</given-names>
</name>
<name>
<surname>Jones</surname> <given-names>H. P.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Nearshore ecosystems on seabird islands are potentially influenced by invasive predator eradications and environmental conditions: a case study at the Mercury Islands, New Zealand</article-title>. <source>Mar. Ecol. Prog. Ser.</source> <volume>661</volume>, <fpage>83</fpage>&#x2013;<lpage>96</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3354/meps13590</pub-id>
</citation>
</ref>
<ref id="B79">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Riddick</surname> <given-names>S. N.</given-names>
</name>
<name>
<surname>Dragosits</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Blackall</surname> <given-names>T. D.</given-names>
</name>
<name>
<surname>Daunt</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Wanless</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Sutton</surname> <given-names>M. A.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>The global distribution of ammonia emissions from seabird colonies</article-title>. <source>Atmospheric Environ.</source> <volume>55</volume>, <fpage>319</fpage>&#x2013;<lpage>327</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.atmosenv.2012.02.052</pub-id>
</citation>
</ref>
<ref id="B80">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rooney</surname> <given-names>N.</given-names>
</name>
<name>
<surname>McCann</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Gellner</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Moore</surname> <given-names>J. C.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Structural asymmetry and the stability of diverse food webs</article-title>. <source>Nature</source> <volume>442</volume>, <fpage>265</fpage>&#x2013;<lpage>269</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nature04887</pub-id>
</citation>
</ref>
<ref id="B81">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Savage</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Seabird nutrients are assimilated by corals and enhance coral growth rates</article-title>. <source>Sci. Rep.</source> <volume>9</volume>, <fpage>4284</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41598-019-41030-6</pub-id>
</citation>
</ref>
<ref id="B82">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Schreiber</surname> <given-names>E. A.</given-names>
</name>
<name>
<surname>Burger</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2001</year>). <source>Biology of Marine Birds</source> (<publisher-loc>Baton Rouge</publisher-loc>: <publisher-name>CRC Press</publisher-name>). doi:&#xa0;<pub-id pub-id-type="doi">10.1201/9781420036305</pub-id>
</citation>
</ref>
<ref id="B83">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shannon</surname> <given-names>C. E.</given-names>
</name>
</person-group> (<year>1948</year>). <article-title>A mathematical theory of communication</article-title>. <source>Bell system Tech. J.</source> <volume>27</volume>, <fpage>379</fpage>&#x2013;<lpage>423</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/bltj.1948.27.issue-3</pub-id>
</citation>
</ref>
<ref id="B84">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sigman</surname> <given-names>D. M.</given-names>
</name>
<name>
<surname>Hain</surname> <given-names>M. P.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>The biological productivity of the ocean</article-title>. <source>Nat. Educ. Knowledge</source> <volume>3</volume>, <fpage>21</fpage>.</citation>
</ref>
<ref id="B85">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Signa</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Mazzola</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Vizzini</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Seabird influence on ecological processes in coastal marine ecosystems: An overlooked role? A critical review</article-title>. <source>Estuarine Coast. Shelf Sci.</source> <volume>250</volume>, <elocation-id>107164</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ecss.2020.107164</pub-id>
</citation>
</ref>
<ref id="B86">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Stanbury</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Thomas</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Aegerter</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Brown</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Bullock</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Eaton</surname> <given-names>M.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>). <article-title>Prioritising islands in the United Kingdom and crown dependencies for the eradication of invasive alien vertebrates and rodent biosecurity</article-title>. <source>Eur. J. Wildlife Res.</source> <volume>63</volume>, <fpage>(1)</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10344-017-1084-7</pub-id>
</citation>
</ref>
<ref id="B87">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Staunton Smith</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Johnson</surname> <given-names>C. R.</given-names>
</name>
</person-group> (<year>1995</year>). <article-title>Nutrient inputs from seabirds and humans on a populated coral cay</article-title>. <source>Mar. Ecol. Prog. Ser.</source> <volume>124</volume>, <fpage>189</fpage>&#x2013;<lpage>200</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3354/meps124189</pub-id>
</citation>
</ref>
<ref id="B88">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Strathmann</surname> <given-names>R. R.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Why life histories evolve differently in the sea1</article-title>. <source>Am. Zoologist</source> <volume>30</volume>, <fpage>197</fpage>&#x2013;<lpage>207</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/icb/30.1.197</pub-id>
</citation>
</ref>
<ref id="B89">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thoresen</surname> <given-names>J. J.</given-names>
</name>
<name>
<surname>Towns</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Leuzinger</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Durrett</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Mulder</surname> <given-names>C. P. H.</given-names>
</name>
<name>
<surname>Wardle</surname> <given-names>D. A.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Invasive rodents have multiple indirect effects on seabird island invertebrate food web structure</article-title>. <source>Ecol. Appl.</source> <volume>27</volume>, <fpage>1190</fpage>&#x2013;<lpage>1198</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/eap.1513</pub-id>
</citation>
</ref>
<ref id="B90">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>USGS</collab>
</person-group>. (<year>2022</year>). <source>EarthExplorer</source>. Available at: <uri xlink:href="https://earthexplorer.usgs.gov/">https://earthexplorer.usgs.gov/</uri> (Accessed <access-date>August 2021</access-date>).</citation>
</ref>
<ref id="B91">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Walsh</surname> <given-names>P. M.</given-names>
</name>
<name>
<surname>Halley</surname> <given-names>D. J.</given-names>
</name>
<name>
<surname>Harris</surname> <given-names>M. P.</given-names>
</name>
<name>
<surname>del Nevo</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Sim</surname> <given-names>I. M. W.</given-names>
</name>
<name>
<surname>Tasker</surname> <given-names>M. L.</given-names>
</name>
</person-group> (<year>1995</year>). <source>Seabird monitoring handbook for Britain and Ireland</source> (<publisher-loc>Peterborough, UK</publisher-loc>: <publisher-name>JNCC / RSPB / ITE / Seabird Group</publisher-name>).</citation>
</ref>
<ref id="B92">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Webber</surname> <given-names>M. A.</given-names>
</name>
<name>
<surname>Jefferson</surname> <given-names>T. A.</given-names>
</name>
<name>
<surname>Pitman</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Jefferson</surname> <given-names>T. A.</given-names>
</name>
</person-group> (<year>2015</year>). <source>Marine Mammals of the World: A Comprehensive Guide to Their Identification</source> (<publisher-loc>San Diego</publisher-loc>: <publisher-name>San Diego: Elsevier Science &amp; Technology</publisher-name>).</citation>
</ref>
<ref id="B93">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Wickham</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2016</year>). <source>ggplot2: Elegant Graphics for Data Analysis</source> (<publisher-loc>New York, USA</publisher-loc>: <publisher-name>Springer-Verlag</publisher-name>). doi:&#xa0;<pub-id pub-id-type="doi">10.1007/978-3-319-24277-4</pub-id>
</citation>
</ref>
<ref id="B94">
<citation citation-type="web">
<person-group person-group-type="author">
<name>
<surname>Wickham</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Fran&#xe7;ois</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Henry</surname> <given-names>L.</given-names>
</name>
<name>
<surname>M&#xfc;ller</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2018</year>). <source>dplyr: A Grammar of Data Manipulation.</source> (<publisher-name>R package version 1.0.8</publisher-name>). Available at: <uri xlink:href="https://CRAN.R-project.org/package=dplyr">https://CRAN.R-project.org/package=dplyr</uri> (Accessed <access-date>June 2021</access-date>).</citation>
</ref>
<ref id="B95">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Young</surname> <given-names>H. S.</given-names>
</name>
<name>
<surname>Hurrey</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Kolb</surname> <given-names>G. S.</given-names>
</name>
</person-group> (<year>2011</year>). &#x201c;<article-title>Effects of seabird-derived nutrients on aquatic systems</article-title>,&#x201d; in <source>Seabird Islands: Ecology, Invasion, and Restoration</source>. Eds. <person-group person-group-type="editor">
<name>
<surname>Mulder</surname> <given-names>C. P. H.</given-names>
</name>
<name>
<surname>Anderson</surname> <given-names>W. B.</given-names>
</name>
<name>
<surname>Towns</surname> <given-names>D. R.</given-names>
</name>
<name>
<surname>Bellingham</surname> <given-names>P. J.</given-names>
</name>
</person-group> (<publisher-name>Oxford University Press</publisher-name>, <publisher-loc>New York, USA</publisher-loc>).</citation>
</ref>
<ref id="B96">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Young</surname> <given-names>H. S.</given-names>
</name>
<name>
<surname>McCauley</surname> <given-names>D. J.</given-names>
</name>
<name>
<surname>Dunbar</surname> <given-names>R. B.</given-names>
</name>
<name>
<surname>Dirzo</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Plants cause ecosystem nutrient depletion via the interruption of bird-derived spatial subsidies</article-title>. <source>Proc. Natl. Acad. Sci.</source> <volume>107</volume>, <fpage>2072</fpage>&#x2013;<lpage>2077</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1073/pnas.0914169107</pub-id>
</citation>
</ref>
<ref id="B97">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zmudczy&#x144;ska-Skarbek</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Balazy</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Kuklinski</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>An assessment of seabird influence on Arctic coastal benthic communities</article-title>. <source>J. Mar. Syst.</source> <volume>144</volume>, <fpage>48</fpage>&#x2013;<lpage>56</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jmarsys.2014.11.013</pub-id>
</citation>
</ref>
<ref id="B98">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Zuur</surname> <given-names>A. F.</given-names>
</name>
<name>
<surname>Ieno</surname> <given-names>E. N.</given-names>
</name>
<name>
<surname>Walker</surname> <given-names>N. J.</given-names>
</name>
<name>
<surname>Saveliev</surname> <given-names>A. A.</given-names>
</name>
<name>
<surname>Smith</surname> <given-names>G. M.</given-names>
</name>
</person-group> (<year>2009</year>). <source>Mixed effects models and extensions in ecology with R</source> (<publisher-loc>New York</publisher-loc>: <publisher-name>Springer</publisher-name>). doi:&#xa0;<pub-id pub-id-type="doi">10.1007/978-0-387-87458-6</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>
