Influence of Climate Change and Trophic Coupling across Four Trophic Levels in the Celtic Sea

Climate change has had profound effects upon marine ecosystems, impacting across all trophic levels from plankton to apex predators. Determining the impacts of climate change on marine ecosystems requires understanding the direct effects on all trophic levels as well as indirect effects mediated by trophic coupling. The aim of this study was to investigate the effects of climate change on the pelagic food web in the Celtic Sea, a productive shelf region in the Northeast Atlantic. Using long-term data, we examined possible direct and indirect ‘bottom-up’ climate effects across four trophic levels: phytoplankton, zooplankton, mid-trophic level fish and seabirds. During the period 1986–2007, although there was no temporal trend in the North Atlantic Oscillation index (NAO), the decadal mean Sea Surface Temperature (SST) in the Celtic Sea increased by 0.66±0.02°C. Despite this, there was only a weak signal of climate change in the Celtic Sea food web. Changes in plankton community structure were found, however this was not related to SST or NAO. A negative relationship occurred between herring abundance (0- and 1-group) and spring SST (0-group: p = 0.02, slope = −0.305±0.125; 1-group: p = 0.04, slope = −0.410±0.193). Seabird demographics showed complex species–specific responses. There was evidence of direct effects of spring NAO (on black-legged kittiwake population growth rate: p = 0.03, slope = 0.0314±0.014) as well as indirect bottom-up effects of lagged spring SST (on razorbill breeding success: p = 0.01, slope = −0.144±0.05). Negative relationships between breeding success and population growth rate of razorbills and common guillemots may be explained by interactions between mid-trophic level fish. Our findings show that the impacts of climate change on the Celtic Sea ecosystem is not as marked as in nearby regions (e.g. the North Sea), emphasizing the need for more research at regional scales.


Introduction
Human-induced climate change has profoundly impacted marine ecosystems across the globe. These impacts have had wide-ranging effects upon the physiology, distribution, phenology and abundance of species, resulting in long-term threats to biodiversity [1]. A key feature of these climate-induced impacts is a high degree of spatial heterogeneity [2], [3]. Understanding the nature of this variation is a key goal for assessing and mitigating the impacts of global climate change.
Environmental change may impact different trophic levels in varying ways [4], [5], such that marine food webs may be impacted both directly and indirectly [6]. Direct effects of climate change include the influence of temperature change, particularly for ectothermic organisms (i.e. fish and invertebrates), or extreme weather events, which can impact endothermic organisms [7]. Taken together these direct effects can influence physiology, morphology and behaviour, leading to a suite of emergent ecological responses [8]. Indirect effects are typically mediated via trophic coupling. This is normally manifested via bottom-up control where climate-mediated changes in the availability of lower trophic levels have knock on consequences for higher trophic levels [4], [9]. Marine higher trophic level predators can also govern the abundance of lower trophic levels by top-down control [10], or mid-trophic level species may exert both top-down and bottom-up effects in a process known as wasp-waist control [11]. It is still unclear, however, how the nature of these effects, as well as ecosystem responses, varies across regions [12], [13].
The aim of this study was to look for climate-related influences across four trophic levels in the Celtic Sea, a productive shelf region in the northeast Atlantic. This is an extremely important area in terms of fish and invertebrate biodiversity [14] and it supports a large community of apex predators in the form of seabirds [15] and marine mammals, as well as several important European fisheries [16]. Profound climate-mediated changes to the nearby North Sea have led to concerns about the long-term viability of certain populations of apex predators [4], [17] so a key question is to determine how climate change might be impacting the Celtic Sea food web. This region is quite different from the North Sea in terms of physical characteristics and oceanography (i.e. general circulation pattern, depth and sea temperature) [18], [19] and this may influence trophic responses to climate change.
We tested for direct and indirect effects of climate across four trophic levels. Direct climate effects at each level of the food web were expected via a significant correlation between abundance/ biomass/demography and environmental variables: the North Atlantic Oscillation index (NAO) and Sea Surface Temperature (SST). In the North Atlantic variations in the NAO index have induced changes in the zooplankton abundance, distribution and community assemblage [2], [20], [21]. In addition, the NAO strongly influences the frequency of extreme weather events, which may directly impact some seabirds [7]. SST may have direct effects through changes to the biology and distribution of ectothermic fish and invertebrates [6].
Indirect effects were expected via bottom-up processes, characterised by a positive correlation between a measure of predator abundance/biomass/demography and prey. We also tested for potential top-down effects typified by a negative correlation between predator and prey. We modelled the Celtic Sea pelagic food web simplified into four trophic levels: four species of piscivorous seabird (black-legged kittiwake Rissa tridactyla, hereafter kittiwake; common guillemot Uria aalge, hereafter guillemot; razorbill Alca torda and Atlantic puffin Fratercula artica, hereafter puffin), pelagic fish (Atlantic herring Clupea harengus, hereafter herring), zooplankton and phytoplankton. For each trophic level, long-term data from 1986 to 2007 were collated and analysed along with measures of environmental conditions (SST and NAO).

Study area
The Celtic Sea is an area of the northeast Atlantic continental shelf, southwest of the United Kingdom ( Figure 1). It represents a transition zone between the Atlantic Ocean and coastal waters of the Bristol Channel and Irish Sea. There is a persistent northwestwards current running from Brittany to the Bristol Channel, as well as oceanographic fronts (the Irish Shelf, the Celtic Sea and Ushant fronts), which tend to inhibit lateral dispersal of phytoplankton [19].

Phytoplankton
The annual phytoplankton bloom in the Celtic Sea typically occurs from April to October [22]. Continuous Plankton Recorder (CPR) data suggest a steady increase in phytoplankton over at least the last 20 years in the region [22]. Phytoplankton productivity and taxonomic composition in the Celtic Sea depend on water column structure: diatoms dominate areas with high nutrient content and display high rates of productivity, while dinoflagellates and microflagellates are found in stratified waters exhibiting lower rates of productivity. Productivity is reasonably high on the shelf with a rapid decrease west of the shelf break [22].

Zooplankton
In the Celtic Sea the large copepod Calanus helgolandicus is an important component of the zooplankton community [23]; CPR data suggest an overall decline in the abundance of zooplankton with Calanus abundance falling below the long term mean in the region [22]. Long-term studies reveal that spatial patterns of zooplankton have changed significantly over the past 40 years, possibly as a result of climate-related reorganization of the zooplankton [24]. The ecological mechanism responsible for these changes remains unclear, however, and further analysis of CPR data in relation to environmental change is needed to clarify the situation.

Pelagic fish
This region is particularly important for pelagic fish such as herring, sardine Sardina pilchardus and sprat Sprattus sprattus, and it is an important spawning ground for key migratory species notably mackerel Scomber scombrus, horse mackerel Trachurus trachurus, and blue whiting Micromesistius potassou [14], [25]. Like many coastal seas, the size-structure of the fish community has changed significantly over recent decades: there has been a decrease in the relative abundance of larger fish with a concomitant increase in the numbers of smaller fish [16]. Henderson (2007) [26] reports two main events in the 1980s and 1990s representing changes in the fish community composition, coinciding with climate-induced changes of plankton community in some regions of the North Atlantic [2].

Seabirds
The Celtic Sea is an extremely important area for seabirds, supporting ,300,000 breeding pairs of 15 species [15], including internationally important populations of northern gannet (Morus bassanus) and Manx shearwater (Puffinus puffinus), as well as nationally or regionally important populations of guillemot, lesser black-backed gull (Larus fuscus), herring gull (Larus argentatus), kittiwake, puffin and razorbill.

Environmental variables
Two environmental predictors were used to test for direct and indirect effects of climate change: the North Atlantic Oscillation (NAO) and Sea Surface Temperature (SST). The description of the ecological mechanism associated with each climate predictor is summarised in Table 1. NAO influences water circulation and sea temperature, which can result in changes to plankton communities [21]. These changes are likely to have effects upon higher trophic levels such as fish and seabirds by influencing food availability [5] or affecting wind, rainfall and air temperature which may consequently influence seabird populations through survival, e.g. by increasing extreme events occurrence [27], [7]; or reducing prey availability within range of breeding colonies which, in turn, can affect the length of seabird foraging trips at the expense of chick survival [28].
Two different versions of the NAO index were used in this study: spring NAO during the seabird breeding season (SNAO, March-June) and winter NAO (WNAO, December-March) in order to test respectively for direct and indirect effects of climate change respectively. Monthly data from 1986 to 2007 were downloaded from the University Corporation from theAtmospheric Research (UCAR) website, Climate Analysis section http://www.cgd.ucar. edu/cas/jhurrell/indices.data.html#naostatmon) and an annual value for each index calculated (Fig. 2a, b).
Variation in SST can affect the marine ecosystem from plankton communities [23], to mid-trophic level fish [29], [9], up to apex predators such as seabirds [30] via match-mismatch events between predators and prey [12], [13]. Both winter and spring SST (WSST, December-March; SSST, March-June) were used as climate predictors. WSST was used to test for direct effects on plankton, and lagged SSST (1-2 years) to test for direct and indirect effects on pelagic fish and indirect effects on seabirds ( Table 2-3). SST data were derived from satellite images and collated from the POET database available at http://poet.jpl.
nasa.gov with a spatial resolution of 0.04u longitude 60.04u latitude.

Plankton data
The Sir Alister Hardy Foundation for Ocean Science (SAHFOS) hosts the Continuous Plankton Recorder (CPR), the world's largest plankton dataset, which SAHFOS has been collecting since 1931 (for more details see [31]). This dataset represents a consistent semi-quantitative index of phytoplankton and zooplankton abundance and has captured seasonal and annual changes in plankton communities (e.g. [4], [2]).
In the present study, a total of 2299 CPR samples, taken between 1986 and 2007, were analysed to investigate possible changes in the plankton community in the Celtic Sea. The abundance of diatoms and copepods was determined in each sample and spring means calculated by averaging across samples taken in the period March-June for the Celtic Sea (49u-53uN, 4u-10uW). Copepod biomass was calculated by multiplying the abundance of each copepod taxon (mainly calanoid; 28 in total between species and taxonomic group see Table S1) by its average mass estimated from an allometric relationship based on size [31].
Diatom abundance was used as a proxy of copepod food availability [4], and copepod biomass was used as a proxy for pelagic fish prey. Copepods were divided into two groups: small copepods (,2 mm) and large copepods (.2 mm) because of changes in pelagic fish diet such as herring during different life stages. Herring larvae (8-10 mm at hatching) feed mainly on copepods and other small planktonic organisms (i.e. fish eggs) [32], [33]. Herrings are selective predators and Calanoid copepods (i.e. Acartia spp., Pseudocalanus) are the predominant prey items during the early juvenile (,3 cm) stage of life with larger herring also consuming predominantly larger copepods (i.e. Calanus spp.) with small fish as gape width increases [34], [35]. Smaller copepods biomass (,2 mm) was used as food proxy for young herring (1 year old, #10 cm Total Length (TL); hereafter 0-group), while . The Celtic Sea region considered in this study (49-53u N, 4-10uW). Including oceanographic features and the seabird colony investigated in this study (Skomer Island: 51u40N, 05u15W). The Irish Shelf Front occurs to the south and west of Ireland and exists all year-round. This front marks the boundary between waters of the shelf (often mixed vertically by the tide) and offshore North Atlantic waters. In addition there are two seasonal fronts systems which tend to develop during spring: the Celtic Sea front (dividing the Celtic Sea from the Irish Sea) and the Ushant Front, which develops from the coast of Brittany and extends to the western English Channel (dividing the Celtic Sea from the English Channel). doi:10.1371/journal.pone.0047408.g001 larger size copepods biomass (.2 mm) for older herring (.1 year old, #15 cm TL; hereafter 1-group).

Fish data
Mid-trophic level fish are vectors for the transfer of energy from low trophic levels to apex predators [10]. In this study long-term data of herring 0-and 1-group abundance were used as a proxy for seabird food availability. Data spanned the period 1986-2007 and were extracted from the Herring Assessment Working Group (HAWG-ICES) acoustic survey designed to evaluate the state of pelagic fish species (sprat and herring) around the UK coast and, in particular, herring stock in the Celtic Sea ( [36]; Fig. 3c).
Long-term data on other small pelagic fish such as sprat or lesser-sandeels (Ammodytes spp.) (hereafter, sandeel) were not available for the region. Nevertheless, data from the Centre for Environment Fisheries and Aquaculture Science (CEFAS) ground fish survey (trawl survey designed to study the distribution, composition and abundance of all fish, commercial shellfish and cephalopod species in the Celtic Sea) has shown that sprat and herring are often caught together; therefore these pelagic species can be considered ecologically equivalent (with similar habitats and diet composition) [37]. In order to test if herring was a good proxy of small pelagic fish species in the Celtic Sea, we analysed the occurrence of sprat and herring in the environment by using the only dataset available of landings from CEFAS ground fish surveys [38] covering the period 1986-2002. This showed that these two species seem to have similar fluctuations in the Celtic Sea region ( Figure S1). Herring is one of the most abundant planktivorous fish in the Celtic Sea, and juvenile stages (0-and 1group) along with other small schooling pelagic fish, such as sandeel or sprat, are an important prey of many seabird species [39]. In the Celtic Sea herring juveniles tend to remain in shallow coastal areas (nursery) for the first two years of their lives [36]. For these reasons we believe that herring (0-and 1-group) represents a good proxy for seabird food availability.

Seabird data
Data on breeding success and population estimates for kittiwake, guillemot, razorbill and puffin breeding on Skomer Island, Wales (51u40N, 05u15W;  average of 3 eggs per breeding attempt, while guillemots, razorbills and puffins are all pursuit divers and lay just a single egg [15]. These four species also differ somewhat in their foraging range: kittiwake, razorbill and guillemot forage mainly inshore [40], [41], [42], whilst puffins tend to forage further offshore [43]. The species also differ in their prey loading: kittiwakes, puffins and razorbills are multiple prey-loaders whereas guillemots are single prey-loaders. Given these ecological differences and their possible diverse vulnerability to food shortage and hence climate change, we modelled their responses separately.

Model construction and statistical analyses
Colinearity among explanatory variables may increase the probability of type I errors, therefore we tested for possible temporal trends in our data using linear regression. Multiple regression models were used to identify the main predictors for each trophic level (diatoms, copepods, herring and seabirds). When the dependent variable was not normally distributed, it was log e transformed. For the four seabird species two demographic rates were used: (1) breeding success, expressed as the number of fledged chicks per breeding pair, weighted by the number of pairs sampled and (2) population growth rate (R t ), which was calculated from the time series as log 10 (N t+1 )-log 10 (N t ) [44].
The following explanatory environmental variables were examined: (1) spring NAO and winter NAO to test for direct and indirect climate effects (spring NAO for seabirds, winter NAO for plankton and fish); (2) winter SST for plankton, lagged spring SST (1,2 years lag) for pelagic fish to test for direct climate effects on each age class, lagged spring SST (1 year lag) to test for indirect climate effects on seabirds. When there was a temporal trend in the dependent variable, year was included as a continuous covariate. The variables used for the model construction are shown in Table 2.
Starting from the full model, the most parsimonious model for each trophic level was selected on the basis of the lowest Akaike Information Criterion (AIC), corrected for small sample size (AICc). This approach selects the model with the best balance between bias and precision and avoids problems of, for example, multiple testing among explanatory variables ( [45]. A set of candidate models was compared using difference in AIC c between the top-ranked and current model (delta AIC c ), and by calculating the AIC c weight (the scaled likelihood that each model is the best description of the data; Burnham and Anderson 2002). Competing models were selected when having their AIC c within 2 of the minimum [45] and are presented in Tables S3, S4, S5. Model goodness of fit was compared using the deviance and coefficient of determination (R 2 ). Covariates were considered statistically significant when the p value was ,0.05. Model residuals were evaluated to check for non-normality, heteroscedasticity and autocorrelation of errors. All modelling was carried out using R version 2.8 [46].

Correlations between covariates
Preliminary explanatory analyses showed that weak correlations occurred in some cases (Table S2); however no significant correlation was found between the environmental variables (both winter and spring measures of NAO and SST) and the other covariates.

Environmental variables
Temporal trends in environmental variables (winter NAO, spring NAO; winter and spring SST) are shown in figure 2. There was no linear trend for either spring or winter NAO over time, but there was considerable inter-annual variability. SST increased significantly over time (winter SST: p,0.001, slope = 0.00660.01; spring SST: p = 0.004, slope = 0.0560.01) with the minimum value in 1986 and the maximum in 2007 (8.5uC and 11.5uC for winter SST, 9.7uC and 12.1uC for spring SST respectively).

Phytoplankton
Diatom abundance (cell count) increased over the period 1986-2007 (p = 0.03, slope = 3.4761.53; Fig. 3a). None of the environ-  mental covariates (winter NAO or winter SST) explained a significant amount of the variation in diatom abundance (Table 3).

Pelagic fish
Herring abundance fluctuated over time, with no clear linear trend (Fig. 3c). Both herring 0-and 1-group did not appear to be regulated by our proxies of food availability (i.e. the biomass of small and large copepods; Table 3), but there was a weak negative climate effect for both groups, with 0-group herring showing a significant negative effect of 1 year lagged spring SST (p = 0.02, slope = 20.30560.125) and 1-group herring a significant effect of 2 year lagged spring SST (p = 0.04, slope = 20.41060.193; Table 3).

Seabirds
Kittiwakes showed highly variable reproductive success, ranging from 0.21 chicks pair 21 in 2001 to 1.01 chicks pair 21 in 2005 (Fig. 4a), but was not significantly correlated with any covariates ( Table 3). The number of breeding pairs has declined significantly on Skomer over the last 22 years (Fig. 4b; p = 0.002, slope = 222.7166.62, R 2 = 0.34) and population growth rate (R t ) was significantly positively correlated with spring NAO (p = 0.03, slope = 0.031460.014).
Puffin productivity did not show a temporal trend (p = 0.37, R 2 = 20.009) with an annual average of 0.75 chicks pair 21 (Fig. 4g), although the number of breeding pairs increased significantly on Skomer over time (p = 0.008, slope = 122.2641.8, R 2 = 0.29; Fig. 4h). None of the explanatory covariates was related to productivity or population growth rate (Table 3).

Discussion
Our results showed both direct and indirect effects of climate change on the Celtic Sea food web, suggesting a weak climate impact from mid-trophic levels to seabirds. In particular, we found that despite changing environmental conditions in the Celtic Sea (i.e. SST is increasing), the response of organisms differed across trophic levels. Increasing SST, for example, had negative impacts on pelagic fish (herring), but did not show any effect on copepods during the study period (1986-2007) and among apex predators, only affected razorbill productivity. Possible mechanisms and explanations for these findings, as well as comparisons with climate related patterns in other regions are discussed below.

Direct Climate change effects on the Celtic Sea pelagic food web
Our results show that diatom abundance in the Celtic Sea also increased during the period 1986 to 2007 (Fig. 3a), although in contrast to other regions in the northeast Atlantic (i.e. Irish Sea [47] this could not be directly linked with changes in climate predictors (winter NAO and winter SST). In the North Sea a species-specific response of diatoms mean annual abundance was shown to be influenced by increasing winds and SST during the last fifty years [48], therefore further research exploring factors influencing changes in phytoplankton in the Celtic Sea could focus on additional climatic indices such as wind and precipitation.
As with the nearby Irish Sea [47], our study suggests that copepod biomass in the Celtic Sea were not significantly related to changes in SST and NAO ( Table 3). As with phytoplankton abundance, this is in contrast to the North Sea where mean annual abundance of calanoid copepods abundance is positively correlated with winter NAO [20], although previous work has revealed strong regional variation in this relationship within the Northeast Atlantic [3]. Previously, Pitois and Fox (2006) [24] argued that climate change had led to a structural reorganisation of zooplankton communities in the Celtic Sea region during the period 1958-2003. Our lack of a strong climate signal in copepod biomass (Table 3) over the period 1986-2007 indicates that either these changes occurred prior to the period of our study , or it is only possible to detect a climate change signal over four decades [2], [49]. Previous studies have found a negative correlation between copepods biomass and rainfall [47], further Previous work has shown that both NAO and SST can strongly impact upon fish growth and abundance [29], [50]. Our results did not indicate an effect of winter NAO on herring 0 and 1group, although both age classes were negatively correlated with spring SST (Table 3). Increasing SST has been found to have both positive and negative effects on small pelagic fish in the Northeast Atlantic [51], [52], [53]. The response of herring to climate change is likely to be latitude-dependent with positive responses at high latitudes and negative at lower latitudes, such as in the Celtic Sea latitudes [54]. This negative effect is likely to be explained by the direct influence of sea temperature on herring spawning and recruitment.
We found a weak positive effect of spring NAO on black-legged kittiwake population growth rate and a weak negative effect on razorbill breeding success ( Table 3). The main driver for this effect is unclear but may be related to the direct effects of wind-speed or storm-frequency [7], [55], both of which are correlated with NAO. Under this scenario strong winds associated with positive NAO phases may differentially affect species such as kittiwake and razorbill because of variations in wing shape and foraging strategies.

Trophic coupling in the Celtic Sea
Multi-trophic level studies from the nearby North Sea [12], [4] as well as in other high latitude shelf regions [5], [30], [56], [57] have shown strong evidence for bottom-up forcing (Table 4). Given the strong increase in SST over this period in the Celtic Sea (Fig. 2c), we had anticipated bottom-up control, however, our results suggest that during 1986-2007 both the plankton community and herring in the Celtic Sea were not strongly regulated via bottom-up forcing (Table 3). Our study supports previous findings, which highlighted that regional variability in the strength of bottom-up control is common, however, and there is evidence that strong variation exists within the North Sea [39], [3]. In general there was no direct evidence that herring was a limiting factor for seabirds in the Celtic Sea, instead a negative correlation between seabird demographics (guillemot R t and razorbill breeding success) and herring abundance (Table 3) suggested top-down control. While we cannot exclude this possibility, an alternative explanation is that this age-class of herring may exert strong top-down effects on other pelagic fish such as sandeels that form an important part of the diet of these two Alcids. A similar trophic mechanism was previously proposed in the North Sea, where herring abundance was negatively correlated with sandeel stocks [58].
A significant negative relationship between razorbill productivity and spring SST lagged by one year (Table 3), suggests indirect bottom-up forcing, since SST is unlikely to directly impact razorbills, but may instead influence the availability of mid trophic level forage fish. Under this situation, SST might be a reliable proxy for overall abundance of forage fish, rather than herring abundance alone. It is unclear why this effect was not shown by the other seabird species, but differences in foraging range and behaviour may explain this.

Conclusions
Previous studies have suggested a strong negative impact of climate change on seabirds elsewhere in the North Atlantic (Table 4; [4] [5]); however the situation in the Celtic Sea appears to be much less clear. Although previous works have demonstrated links between climatic conditions and seabird demographics in this region [55], [59], [60], these have mostly been connected with changes in NAO indices and are consistent with direct weather effects mediated by changes in storm frequency or wind conditions. Nevertheless, a recent study revealed a link between warmer waters in the Celtic Sea and offspring condition in Manx shearwater (Puffinus puffinus) [61]. Therefore, the role of climate change on the Celtic Sea remains unclear but it certainly does not appear to share the same very strong signal exhibited elsewhere in the North Atlantic. Moreover, although kittiwake numbers have decreased significantly at Skomer in the past two decades, the numbers of the three Alcids has increased (Fig. 4). However, declines in the breeding success of these Alcids are perhaps reason for concern and this could be linked with density-dependent effects.
We reiterate the call for future research to focus upon multitrophic level, region-wide research in order to understand the ecological processes regulating marine food webs in response to climate change. However, data availability is a common limitation in this approach and there are still only a small number of studies that have used these combined long-term datasets (representing all trophic levels from plankton to apex predators) in the North Atlantic (e.g. [12], [13], [4]). However, we also urge that such ecosystem level approaches should also investigate the potential for synergistic effects of fisheries on climatic impacts. Marine ecosystems are not equally sensitive to climate change, with some regions more vulnerable than others [62]. This study has important implications for our understanding of climate change impacts on marine ecosystems and in particular on apex predators, highlighting the importance of regional variability even within a relatively small geographic area (i.e. North Sea and Celtic Sea).  Significance is indicate as follow: pvalue,0.001 ***, pvalue,,0.01**, pvalue,0.05* SNAO: spring North Atlantic Oscillation index; WNAO: winter North Atlantic Oscillation index; SSST: spring Sea Surface Temperature; WSST: winter Sea Surface Temperature; Small cop: small copepods (,2 mm); Large cop: large copepods (.2 mm); KittBS: blacklegged kittiwake productivity; GuiBS: guillemot productivity; RazBS: razorbill productivity; PufBS: puffin productivity; Her 0g: herring 0-group; Her-1g: herring 1-group; Large cop: large copepod; Small cop: small copepods; KittRt: black-legged kittiwake population growth rate; GuiRt: guillemot population growth rate; RazRt: razorbill population growth rate; PufRt: puffin population growth rate. (DOCX)