Postnatal growth rate varies with latitude in range‐expanding geese: The role of plasticity and day length

Abstract The postnatal growth period is a crucial life stage, with potential lifelong effects on an animal's fitness. How fast animals grow depends on their life‐history strategy and rearing environment, and interspecific comparisons generally show higher growth rates at higher latitudes. However, to elucidate the mechanisms behind this gradient in growth rate, intraspecific comparisons are needed. Recently, barnacle geese expanded their Arctic breeding range from the Russian Barents Sea coast southwards, and now also breed along the Baltic and North Sea coasts. Baltic breeders shortened their migration, while barnacle geese breeding along the North Sea stopped migrating entirely. We collected cross‐sectional data on gosling tarsus length, head length and body mass, and constructed population‐specific growth curves to compare growth rates among three populations (Barents Sea, Baltic Sea and North Sea) spanning 17° in latitude. Growth rate was faster at higher latitudes, and the gradient resembled the latitudinal gradient previously observed in an interspecific comparison of precocial species. Differences in day length among the three breeding regions could largely explain the observed differences in growth rate. In the Baltic, and especially in the Arctic population, growth rate was slower later in the season, most likely because of the stronger seasonal decline in food quality. Our results suggest that differences in postnatal growth rate between the Arctic and temperate populations are mainly a plastic response to local environmental conditions. This plasticity can increase the individuals' ability to cope with annual variation in local conditions, but can also increase the potential to re‐distribute and adapt to new breeding environments.


| INTRODUC TI ON
The period of postnatal growth is a crucial stage in an animals' life cycle, with a clear link to fitness components such as reproduction and survival (Dmitriew, 2011;Haywood & Perrins, 1992;Starck & Ricklefs, 1998;. The rate of growth is a basic life-history trait, which is shaped by environmental conditions at the local breeding and rearing grounds as well as by an animal's life-history strategy (Arendt, 1997;Dmitriew, 2011). According to life-history theory, animals are expected to adapt growth rates to local conditions to maximize fitness. However, to cope with environmental variation, such as food availability and quality, growth rate also has to be flexible (Arendt, 1997;Dmitriew, 2011). Therefore, understanding which environmental factors influence growth rate and its flexibility is an important step in evaluating an animal's vulnerability and adaptability to environmental change.
In seasonal environments, the length of the growing season, which is negatively correlated with latitude, restricts the period available for growth. Correspondingly, a positive relationship between growth rate and latitude has been found in a variety of organisms, that is in fishes (Brown et al., 1998;Conover & Present, 1990), insects (Blanckenhorn et al., 2018;Kojima et al., 2020), amphibians (Lindgren & Laurila, 2005) and birds. The time constraint imposed by season length is especially pressing for Arctic-breeding migratory birds, because offspring have to be fully developed and capable of leaving the breeding area before winter sets in (Alerstam & Hedenström, 1998;Owen & Black, 1989;Tomotani et al., 2016). By migrating between breeding and wintering areas, migratory birds are able to exploit seasonally occurring food peaks and avoid local food scarcity and harsh climatic conditions (Holt & Fryxell, 2011). Thus, performing migratory journeys is expected to enable fast growth, while also imposing the need to realize it.
Birds have been shown to benefit from fast growth since it shortens the period of vulnerability to size-dependent predation (Dmitriew, 2011;Samelius & Alisauskas, 1999;Starck & Ricklefs, 1998). Furthermore, within bird populations, faster growth has been associated with larger adult size (Cooch et al., 1991;Larsson et al., 1998;Searcy et al., 2004; with fitness consequences throughout an individual's lifetime. However, growth itself might be costly. Fast growth can reduce resistance to starvation, increase cellular damage imposed by oxidative stress and reduce immune functioning, all of which may impact a bird's life span and functioning (Arendt, 1997;Dmitriew, 2011;Kim et al., 2011;Mangel & Munch, 2005). Among the environmental factors that control growth rates, food availability plays a central role. In birds, periods of food shortages have been shown to negatively affect muscle development and body mass increase (Killpack & Karasov, 2012) and, in strongly seasonal environments like the Arctic, a mismatch with the peak in food quality has been shown to result in slower growth (Brook et al., 2015;Ross et al., 2018).
The high productivity of the Arctic summer is an important prerequisite for successfully raising offspring in herbivorous and insectivorous species (Fokkema et al., 2020), and is considered a main driver of migration to the Arctic (Sedinger & Raveling, 1986). Furthermore, 24-hr daylight during Arctic summers dramatically improves the potential feeding time of animals that rely on eyesight to forage (Schekkerman et al., 2003). Combined, the high productivity and unlimited feeding time in the Arctic result in high resource availability for Arctic-breeding birds. In interspecific comparisons, higher growth rates have been reported for Arctic-breeding waders such as red knot Canutus canutus and little stint Calidris minuta as compared to temperate-breeding waders such as redshank Tringa totanus, lapwing Vanellus vanellus and black-tailed godwit Limosa limosa (Schekkerman et al., 2003;Tjørve, 2007). The same pattern is found in altricial gulls and terns (Larus and Sterna spec.), where two populations of the same species (Sterna paradisaea and Larus argentatus) show a positive relation between latitude and growth rate (Klaassen, 1994;Tjørve, 2007).
Although Schekkerman et al. (2003) mention the potential importance of day length and arthropod abundance for growing waders, the role of resource availability in explaining latitudinal differences in growth rate was not evaluated in detail in the aforementioned studies.
Interspecific comparisons, like above, suffer from the fact that species are also bound to differ in other respects than breeding environment alone (Garland & Adolph, 1994). These limitations therefore call for studies using intraspecific comparisons across environments.
Here, we make a within-species comparison of growth rates of barnacle goose goslings among three different populations (Barents Sea, Baltic Sea and North Sea) spanning 17° in latitude. These populations are genetically very similar and show substantial gene flow (Jonker et al., 2013). The Russian flyway population of barnacle geese has shown a strong increase over the past decades (over 7% annual increase since 1960; Madsen et al., 1999), and simultaneously expanded its traditionally Arctic breeding range by establishing new breeding colonies at stopover sites in the Baltic region and in the wintering area along the North Sea coast (Larsson et al., 1988;Van Der Jeugd et al., 2009). Barnacle geese breeding in the Baltic region shortened their migratory distance considerably compared to Arcticbreeding geese, while barnacle geese breeding along the North Sea coast became sedentary. Besides differences in migratory strategy, geese from these populations also experience differences in their local breeding environments such as season length, day length and feeding conditions. Outside the breeding season, geese of all three populations share common wintering grounds along the North Sea coast. The rapid range expansion of the barnacle goose can be seen as a unique natural experiment, which allows to investigate how animals cope with new or changing environments by adopting new life-history strategies. We relate the differences in growth rate to differences in environmental conditions at the breeding grounds and evaluate potential environmental constraints within the different populations.
Furthermore, we assess whether differences in gosling growth among populations can be the result of microevolution or are to be attributed to developmental phenotypic plasticity (i.e. the ability of an individual to adapt to novel circumstances through flexible expression of a trait; Dobzhansky, 1970). Finally, we compare the latitudinal gradient in growth rates observed in barnacle geese to the latitudinal gradients observed in precocial waterfowl and waders based on previously published growth rates.

| Data collection
We collected biometric data on growing goslings during longterm studies in colonies from three study populations ( Figure (Eichhorn et al., 2010;Larsson et al., 1988;Van der Jeugd et al., 2003 Table S1). Sex was determined based on cloacal inspection. Goslings were weighed in a bag using a Pesola spring scale with an accuracy of ±5 g (if <600 g) or a digital hand scale or Pesola spring scale with an accuracy of ±10 g (if >600 g). A calliper (±0.1 mm) was used to measure the outer length of the bent tarsus. Head length was measured using a ruler (±1 mm). Body mass and body size are correlated ( Figure S2), but are generally analysed separately when modelling growth (Starck & Ricklefs, 1998;Tjørve & Tjørve, 2010

| Growth models
We modelled gosling growth separately for male and female goslings using a Gompertz model (Gompertz, 1825), which is commonly used for precocial species (Schekkerman et al., 1998;Tjørve & Tjørve, 2010, with a fixed initial value (hatching size) as proposed by Tjørve and Tjørve (2017), using the following formula: In this formula, biometric size (head length in mm, tarsus length in mm or body mass in g) is modelled as a function of age (in days).
We chose to use the Gompertz equation, because it has been applied successfully to waterfowl data in the past and allows for easy comparison with other studies (Sedinger, 1986). In the equation, A represents the asymptote, which has been fixed to the average adult size for males and females (Table S2), respectively, as suggested by Austin et al. (2011). Asymptotic values were based on measurements of adult geese caught during moult in the three study areas. Because adult size and body mass were similar for all three populations (see Supporting Information), we used adult size and body mass averaged over all three populations. The size at hatch (when age = 0) is given by I, which is calculated based on measurements taken immediately upon hatch. It replaces the inflection point parameter in the original Gompertz function (Tjørve & Tjørve, 2017). We used the same value for I in models for males and females of all three populations, since we were not able to distinguish between males and females at day 0.
Differences in average hatching size were smaller than 1 g or 1 mm between the Barents Sea and North Sea populations, so we used the same averages for all three populations. Parameter k represents the growth coefficient, and is estimated by the model. In the (few) cases of multiple recaptures of an individual, only data from the first capture were used in our analyses to avoid potential bias in recapture data due to repeated handling stress. Nest (brood) ID was included as a random effect on k to account for statistical dependence due to genetic background, similar rearing environment and potential maternal effects (Sofaer et al., 2013). Since data were collected over multiple years, we also included random cohort effects to account for variation caused by annual differences in phenology and growing conditions. We nested the random effect of nest within the effect of cohort: where k i represents the random cohort effect and k ij the random nest effect. Random nest and cohort effects and their respective errors were expected to be normally distributed with a mean of zero.
Growth models were constructed using a nonlinear mixed effects model approach, using the 'nlme' package in r (Pinheiro et al., 2012; R Development Core Team, 2010).

| Comparing populations
Specific testing of differences in gosling growth rate between populations was done by adding dummy variables for the three populations to expression (2) as proposed by Sofaer et al. (2013).
Here, P Ba and P NS are the dummy variables for the Baltic and North Sea population (represented by 1 or 0) and k Ba and k NS are the . ( . population-specific differences to the k of the Barents Sea population.
In this way, we could determine parameter estimates describing the difference in growth rate among populations. We constructed separate models for male and female goslings since the sexes have different asymptotes.

| Day length
The number of daylight hours that had accumulated between hatching and capture was calculated for each gosling. Daylight was determined as the period between dawn and dusk, and was calculated based on the coordinates of the three breeding colonies using the r package 'Suncalc' (Thieurmel & Elmarhraoui, 2019). To model biometric size as a function of daylight hours, we used the same formula as expression (3), replacing 'age in days' by 'daylight hours' experienced by each individual gosling.
In addition to our analysis with fixed population effects included in a random Gompertz model (following Sofaer et al., 2013), we analysed our data using GLMMs on the residuals of non-random Gompertz models for males and females (see Supporting Information).

| Effect of hatch date on growth
The effect of hatch date on growth was analysed using the residuals of the non-random Gompertz models (expression (1)), hereafter referred to as 'residual head length', 'residual tarsus length' and 'residual body mass'. The residuals were calculated by subtracting the expected value of morphometric size of a gosling at a certain age based on the fitted growth curve from the observed size. Residuals of males and females were analysed collectively. We calculated relative hatch dates by centralizing hatch dates within each cohort, because years can differ in onset of spring and consequently in timing of breeding and hatching.
For the calculation of the relative hatch date for each gosling, we used the mean hatch date of the colonies (not only of the recaptured goslings), as established from nest monitoring (see Supporting Information for details). We constructed GLMMs with fixed effects for population, hatch date and their interaction. Sex was added as fixed effect to account for potential differences between the residuals of models for males and females. NestID and cohort were included as random effects, with nestID nested in cohort. We used a backward selection procedure using Akaike information criterion (AIC) to exclude factors that did not contribute to the fit of the model.

| Phenotypic plasticity or evolutionary response
To investigate whether any observed differences in growth rate among the study populations can be the result of microevolution or have to be (partly) attributed to phenotypic plasticity, we expressed the rate of change in haldanes (Gienapp et al., 2007). The haldane expresses the rate of change per generation in phenotypic standard deviations (SD) and is calculated with the formula given in expression (4).
Here, h represents the phenotypic change in haldanes, X 1 and X 2 are the trait mean values of two populations (synchronic comparison), S p is the pooled standard deviation from X 1 and X 2 , and g is the number of generations (Hendry & Kinnison, 1999).
We used the Gompertz growth rate of the Barents Sea and North Sea population for X 1 and X 2 , since these two populations are expected to represent the largest difference. S p is calculated using the standard deviations estimated by the growth models. The number of generations is calculated based on a generation time for barnacle geese of 7.5 years (Dillingham, 2010), and a period of change of 30 years (period from 1985 till 2015 in which the establishment of the North Sea barnacle goose colony took place).

| Population comparisons
The growth rate of gosling body mass, head length and tarsus length were found to differ among the three populations studied (Figure 1, Figure S3A,B, Table 1).
In both males and females, body mass growth in the Barents Sea The populations differences reported above are supported by the GLMM analysis on the residuals of the non-random Gompertz models. The models that included 'population' had consistently lower AICc values than models that did not include 'population' (Table S3).

| Day length
Differences in gosling growth rate among populations were largely explained by differences in day length ( Figure 2, These results were in line with the GLMM analysis on the residuals of the non-random Gompertz models. When comparing the AICc values of the models with 'daylight experienced since hatch', 'population' was not retained in the most parsimonious models for head length and tarsus length in both males and females. For body mass, 'population' was retained in the most parsimonious model for both sexes (Table S3).
The interaction effect of relative hatch date and population on residual head length was retained in the model, but was not significant (F 2, 639 = 2.76, p = 0.064; Figure 3b). The pattern was similar to the pattern found for body mass, that is, later hatched goslings tended to grow slower in the Barents Sea and Baltic Sea population, while this trend was absent in the North Sea population.
The interaction effect of relative hatch date and population was excluded from the model on residual tarsus length. However, the best model did contain relative hatch date as independent variable, showing a significant negative effect (F 1, 641 = 37.09, p < 0.001; Figure S4).

| Phenotypic plasticity versus micro-evolution
The differences in growth rate between the Barents Sea and North F I G U R E 2 Gompertz growth models for body mass and head length in relation to hours of daylight experienced by each gosling. Growth models for males are shown in panels (a) and (c); for females in panels (b)

| Intraspecific and interspecific patterns of growth rate with latitude
The intraspecific growth rate differences we found are in line with interspecific differences found among gulls, terns and waders, which also show an increase in chick growth rate with increasing breeding latitude (Klaassen, 1994;Schekkerman et al., 2003;Tjørve, 2007;Tjørve et al., 2009). When we limit our comparison of growth coefficients along a latitudinal gradient to precocial species (Table S4), because altricial species are known to grow faster (Starck & Ricklefs, 1998), and correct the growth rate for the Log(body mass) of each species, respectively, since growth rate scales with body mass (Tjørve, 2007), we find that the intraspecific linear increase in growth rate with latitude for barnacle geese is similar to the interspecific pattern in Charadriiformes and other Anseriformes ( Figure 4). Furthermore, this pattern holds regardless of foraging guild (wader chicks being insectivorous and waterfowl being herbivorous), confirming findings of Tjørve (2007) for waders and gulls.

| Resource availability and phenotypic plasticity
Our results suggest that the observed differences in growth rate are mainly the result of differences in resource availability. Continuous daylight during the Arctic summer increases potential feeding time for Arctic-breeding geese with approximately 8 hr (33%) compared to feeding time for temperate breeding geese, and with approximately 2 hr (10%) for the Baltic Sea population. Barnacle geese make use of the extended day length by adapting their circadian rhythm . The effect of daylight on growth has extensively been shown in poultry, where increasing day length led to increased food consumption resulting in higher growth rates (Kleinpeter & Mixner, 1947;Wineland, 2002). This increased food intake is expected to require a larger metabolic machinery. In order for goslings to benefit from this food peak, timing of reproduction is essential (Lameris et al., 2017;Nolet et al., 2020;Van der Graaf et al., 2006;Van der Jeugd et al., 2009), as is illustrated by the negative effect of hatch date on residual head length and tarsus length, and its interaction with population. We found the strongest negative effect of hatch date on residual body mass in the Barents Sea population, a weaker but still significant negative effect in the Baltic Sea population, and no effect of hatch date in the North Sea population. Although differences among populations were not F I G U R E 4 Relationship between latitude and Gompertz growth coefficients in precocial species in the Northern hemisphere. The growth coefficient is corrected for the LOG(body mass) to make species of different sizes comparable. Growth coefficients of Anseriformes are given by filled squares (see Table S4 for references). Red squares represent the growth coefficients of our barnacle goose study populations. Growth coefficients of Charadriiformes (open circles) have been retrieved from Tjørve (2007) and Tulp (1998 shown before in Arctic geese (Cooch, 2002;Doiron et al., 2015;Gauthier et al., 2006;Sedinger & Flint, 1991), we show here that this impact increases with latitude and is absent in the sedentary temperate population. Lindholm et al. (1994) showed experimentally that the decrease in growth rate of later-hatched goslings in the Arctic is mainly the result of a decrease in forage quality. While birds in the Barents Sea population are generally able to utilize the food quality peak, hatching in the North Sea and Baltic population occurs too late (Van der Jeugd et al., 2009). Hence, the food peak may not only be higher in the Arctic than in temperate areas, but breeding may also be better timed so goslings can profit from it. Larsson and Forslund (1991) showed that differences in food quality not only affect the growth of barnacle goose goslings but also their final adult size. Similar results were found in lesser snow geese Anser caerulescens caerulescens and black brants Branta bernicla nigrans (Cooch et al., 1991(Cooch et al., , 1996Sedinger & Flint, 1991). This developmental plasticity itself is adaptive, as it allows a growing individual important leeway when environmental conditions limit growth and the programmed size is out of reach.
Although microevolution can be fast under strong selection (Endler, 1980), especially when strong founder effects occur (Grant & Grant, 1995), the differences in growth rate observed between the Arctic and North Sea population appear too high to be caused by micro-evolution alone, in the relatively short time span between the time of establishment of the new populations and the time of our measurements. Only the smallest difference in growth rate between both populations (male head length and tarsus length) remained within the limits of plausible microevolution, with 0.077 haldanes still being higher than 75% of 2,414 evolutionary rates reported by Hendry et al. (2008).
With values over 0.10 haldanes (for differences in body mass and female head and tarsus length) being higher than 97% of the evolutionary rates reported by Hendry et al. (2008), plasticity appears to be the main mechanism behind the observed differences in growth rate.
Although our results indicate a prominent role for phenotypic plasticity in response to resource availability, it is not possible to fully disentangle the contributions of plasticity and selection. The higher growth rate found in the Barents Sea population might be the result of within season selection for faster growth, because smaller goslings may be more likely to die before recapture due to predation and adverse weather conditions. Additionally, the faster body mass growth in the Barents Sea population, even after correcting for the effect of daylight, might reflect selection for faster growth to ensure goslings fledge in time to escape harsh weather conditions with the autumn migration. Larger, faster growing goslings have higher survival prospects on their first autumn migration (Loonen et al., 1999).
Van der Jeugd et al. (2009) (Larsson & Forslund, 1991), and may obscure potential population differences. By including random effects for cohort, we corrected for variation in growing conditions among years. To account for a trend in body size over the study period as a result of an increasing mismatch with the peak in food availability due to climate change (Doiron et al., 2015;Nolet et al., 2020), we checked the residuals of the non-random Gompertz curves for the Baltic Sea population, which has the best data coverage over time. We found no trend for the residuals of head length and tarsus length, but found a negative trend for the residuals of body mass, although its effect was small (−5.6 g/year decrease), corresponding with a decrease of 87 g over the 15-year study period (around 5% of the mean asymptotic body mass).
Unfortunately, data coverage in the Barents Sea and North Sea population is not (yet) adequate to check whether this may be a general effect that can indeed be attributed to climate change but this may be confirmed in a later study.

| CON CLUS IONS
Our results show that goslings from an Arctic (migratory) population grow faster than goslings from a temperate (non-migratory) population, while goslings raised at an intermediate latitude show intermediate growth rates. Our analysis suggests that these differences are caused by a plastic response to local environmental conditions such as day length and food quality. However, it is not possible to fully disentangle these effects from micro-evolutionary adaptations of growth rate to latitude without doing experimental studies. One way to bring this further is to set up a 'common garden' and study the growth of goslings from different breeding populations under the same rearing conditions.
The differences we show in growth rates of goslings in Arctic and temperate populations of the same species help unravel the costs and benefits of a migratory lifestyle. The costs involved in completing a migratory journey should be balanced by fitness benefits. Changing conditions in both the Arctic and temperate zone can influence the cost-benefit balance of a migratory journey. In their temperate wintering and breeding sites, barnacle geese profit from managed grasslands that provide a diet of improved food quality (Eichhorn et al., 2012). Climate warming, on the other hand, pushes migratory geese to their limit to arrive at their Arctic breeding grounds in time to ensure their goslings can profit from the food peak (Lameris et al., 2018). Under these developments, the benefits of migration might not outweigh the costs any longer, whereas the costs of breeding in temperate areas may not be as high as they used to be. Plasticity in growth rates can be an important factor enabling species to be flexible enough to adapt to new or rapidly changing breeding environments.

ACK N OWLED G EM ENTS
We are grateful to a large flock of fieldworkers that helped with the catching of geese in the Netherlands and Sweden and our collaborators from the Bird Ringing Centre Russia for help with catching geese in Russia.

CO N FLI C T O F I NTE R E S T
The authors have no conflict of interest to declare.

DATA AVA I L A B I L I T Y S TAT E M E N T
Data available from the Dryad Digital Repository https://doi. org/10.5061/dryad.qjq2b vqhc .