Translator Disclaimer
1 September 2012 Spatio-temporal relationship between calf body mass and population productivity in Fennoscandian moose Alces alces
Raisa Tiilikainen, Erling Johan Solberg, Tuire Nygrén, Jyrki Pusenius
Author Affiliations +
Abstract

Body mass is an important life history trait related to survival, mating success and fecundity in ungulates. Accordingly, we may expect that both body mass and reproductive measurements at the population level can be used as valid indices of population condition. However, several factors may modify the relationship between body mass and fecundity because of trade-offs between maturity and early body growth, and varying mortality patterns and sex/age structure among populations. To evaluate the use of such indices for population monitoring and examine the current variation in moose Alces alces population condition in Fennoscandia, we studied the spatio-temporal relationship between calving rate, twinning rate and average autumn calf body mass of moose in Norway and Finland. Calving rate and twinning rate were based on moose observations by hunters while body mass was the average carcass mass of harvested calves. We found a positive relationship between indices both within and among populations. Calves were on average heavier and the observed recruitment rates higher in Finland than in Norway, which is consistent with the higher moose density and presumably lower primary productivity (higher altitude) of moose ranges in Norway. We also found higher observed recruitment rates in populations and years with more even adult sex ratios (females per male) and low relative harvest rates of calves. This suggests that variation in recruitment rate is not only a matter of nutritional condition, but is also affected by varying hunting regulations and harvest structure. For monitoring purposes, we believe that twinning rate is best suited for ranking populations according to nutritional status as this index is closely related to fitness and is relatively insensitive to variation in perinatal and harvest mortality. However, variation in calf body mass may better reflect temporal variation in living conditions. This is because early body growth is sensitive to variation in food availability (and quality) and because body mass may respond more instantaneously than recruitment indices to adverse conditions. Accordingly, we found both calving rate and twinning rate to be best related to variation in mean calf body mass in the previous year.

Body mass is an important life history trait in ungulates as large body masses are related to high survival, mating success and fecundity (Stearns 1992, Bérubé et al. 1999, Gaillard et al. 1996, 2000b). In many species, larger females also produce larger than average offspring (Gaillard et al. 1992, Sæther & Haagenrud 1985, Sand 1996, but see Gaillard et al. 2000a), which correspondingly tend to become larger adults (Albon et al. 1987, Festa-Bianchet et al. 2000, Beckerman et al. 2002, Solberg et al. 2004, 2007). The positive relationship between early growth, and adult size and fecundity suggests that body mass variation of young individuals should influence the population dynamics (e.g. Gaillard et al. 2000b, Hewison & Gaillard 2001), and hence a positive relationship between the average size of juveniles and population productivity could be expected (Kie et al. 2003, Stewart et al. 2005).

This assumed positive relationship between life history traits is often utilised for monitoring exploited ungulates (Hanks 1981, Solberg et al. 1999, Herfindal et al. 2006, Boertje et al. 2007). For instance, by measuring the variation in average body mass (or condition) of harvested juveniles, wildlife managers can keep track of the population vigour and nutritional status (e.g. Morellet et al. 2011). However, like the relationship between life history traits within individuals, several factors may modify the relationship between body mass and fecundity at the population level. First, a trade-off between maturity and body growth may be shaped by different mortality (harvest) patterns among populations (Proaktor et al. 2007). Accordingly, females may mature at different ages and sizes among populations, leading to varying relationships between juvenile size, age at first reproduction and adult fecundity at the population level (Stearns 1992, Garel et al. 2009, Monteith et al. 2009). Second, high productivity can be associated with low offspring size in species that produce multiple young (Stearns 1992, Nygrén 2003) suggesting that populations with high productivity are not necessarily associated with larger offspring. The fact that ungulates are long-lived animals may also complicate the matter. Following improved living conditions, fecundity and offspring size may reflect past rather than current conditions, which in turn can have different effects on the size and number of offspring. Indeed, although small females tend to produce smaller offspring, they may still produce large litters when the environmental conditions are good (e.g. Monteith et al. 2009).

In Fennoscandia, the abundance of moose Alces alces has been growing during the last 50 years, leading moose to become the most important wildlife resource (Lavsund et al. 2003). Following this increase, there has been growing concern about the potential impact of density-dependent food limitation, and for the same reasons, populations are now regularly monitored for variation in population density and structure and nutritional status. Utilising data provided by hunters (i.e. harvest and observation data) has been the favoured approach for this monitoring. These data include carcass masses of harvested moose and the sex and age of moose observed during hunting from which indices of body growth, reproductive performance, density and adult sex ratio are derived (e.g. Lavsund et al. 2003, Morellet et al. 2011).

The most common indicators of population nutritional status are the mean carcass mass of harvested calves, as well as the calving rate (proportion of females with calf/calves) and twinning rate (females with twins/females with ≥ 1 calf) of observed females. To indicate population nutritional status, body mass of calves are usually preferred to the mass of moose in older age groups. This is because calves constitute the most numerous age group in the harvest (our study), can be easily aged by their size and tooth pattern (Rolandsen et al. 2008), and because the sample of harvested calves is not assumed to be affected by hunting selectivity (e.g. Moe et al. 2009). By being the youngest, smallest and arguably the weakest individuals in the population, we may also expect them to quickly respond to annual variation in living conditions (Gaillard et al. 2000b, Herfindal et al. 2006).

Large numbers of moose observations are inexpensive to collect and data can be easily converted to indices of population nutritional status, such as calving rate and twinning rate (i.e. recruitment rates). However, because these data are collected in autumn (i.e. about 4-6 months after calving) and while hunting, their ability to reflect the variation in population condition with adequate precision has been questioned. Large variation in predation pressure can for instance lead to large variation in recruitment rates among regions and over time; variations that are mostly unrelated to population nutritional status (Boertje et al. 2007). In addition, annual variation in female age structure and adult sex ratio (Solberg et al. 2002, Sæther et al. 2004, Nygrén 2009) can affect population productivity. These effects are particularly relevant for harvested species because biased harvesting of sex and age groups can create large variation in population structure (e.g. Ginsberg & Milner-Gulland 1994, Solberg et al. 1999, Mysterud et al. 2005, Nygrén 2009).

In our study, we analysed the spatio-temporal variation in moose body mass and recruitment rates at the population level, and evaluated the usefulness of these data for monitoring purposes. We did this by first testing the extent to which body masses of calves and calving/twinning rates (recruitment indices) varied among regions and countries. In general, we expected the calf body mass and recruitment indices to be higher in Finland than in Norway, given the current higher densities of moose in Norway (see below), and males to be heavier than females. Next, we tested if the recruitment indices were related to the variation in calf body mass within (temporal relationship) and among regions (spatial relationship), while simultaneously controlling for varying hunting selection and population sex ratios. To better understand the functional relationship between body mass and recruitment rates, we tested the temporal relationship at several time lags. If both calf mass and fecundity depend on resources available in the year of reproduction (as in income breeders; Jönsson 1997), we expected the recruitment rates to be positively correlated with the mean calf body mass in the same year. However, as calf body mass can have long-lasting effects on adult size and fecundity (Keech et al. 1999, Solberg et al. 2004), the temporal variation in recruitment rate may also be a time-delayed response to a change in calf body mass. The latter relationship is generally assumed among moose managers, and could favour monitoring juvenile body mass as a means to foresee changes in population productivity. However, convincing support for such a relationship has yet to be presented.

Material and methods

Study areas in Finland and Norway

Our study area covers populations from most of the distributional range of moose in Norway and all of Finland (Fig. 1). In Norway, our study areas are part of the national monitoring programme for cervids (Solberg et al. 2006), whereas in Finland, our study areas constitute the 15 game management districts (Nygrén & Pesonen 1993). The data for body mass of calves and moose observations were aggregated from subregion-level (municipality in Norway and game management association (mostly based on municipalities) in Finland) to these region-level mean values.

Finland and Norway are found at the same latitudes (see Fig. 1), and are covered mainly by the same vegetation types. Both countries are dominated by coniferous forests. However, Norway has a more oceanic climate than Finland, and due to the more mountainous topography, moose areas are generally found at higher altitudes. For more details on our study areas in Norway, see Solberg et al. (2002) and Garel et al. (2009), and in Finland, see Nygrén et al. (2007).

Moose population density, management and monitoring in Finland and Norway

In both Norway and Finland, moose populations began to increase in the 1960s and reached a peak in Norway in the 1990s. In Finland, the moose population has peaked twice: first in the late 1970s and early 1980s and again in the late 1990s and early 2000s (Lavsund et al. 2003, Nygrén 2009). Currently, the density of moose in Norway is almost twice the density in Finland (approximately 0.8-1.1 moose/km2 in Norway and 0.2-0.5 moose/km2 in Finland), mainly because of the official enforcement of relatively moderate density thresholds in Finland (Nygrén 1984, Nygrén & Pesonen 1993, Lavsund et al. 2003, Nygrén 2009). Otherwise, the principles of moose management are quite similar in these countries (Lavsund et al. 2003). Moose hunting occurs in autumn at much the same period (September-December), but the hunting season has varied somewhat over time and is much shorter in Norway (Nygrén et al. 2007, Garel et al. 2009). Hunting is mainly carried out by a group of hunters in a team, often with the use of dogs and/or beaters to flush out moose from cover (Koskela & Nygrén 2002, Ruusila & Pesonen 2004).

Moose harvesting is based on a quota system and the number of harvested moose is obligatorily recorded and reported according to sex and age (calf, yearling and adult (≥ 2 years old) in Norway, calf and adult (≥ 1 year old) in Finland), and moose observations (see below) are collected by hunters as a standard routine during the hunt (Lavsund et al. 2003). In both countries, body mass (carcass mass) of harvested moose has been systematically collected in many areas over several years (Solberg et al. 2006, Nygrén et al. 2007; see Fig. 1). We only used data from 1991-2006 in our study so that our study period was the same in both countries.

Observation data

Moose observation monitoring in Norway and Finland is a systematic recording of the sex and age (calf or adult) of moose observed by moose hunters during the hunting season (Lavsund et al. 2003). In both countries, hunting effort is measured as the number of hunting days, either by the hunting team (Finland) or by the individual hunters (Norway). The observations are classified according to six categories: calves (in Finland calculated from observed females with one calf or twins), adult (≥ 1 year old) males, adult females without calf, females with one calf, females with twins and individuals of unknown sex and/or age (Nygrén 1984, Nygrén & Pesonen 1993, Solberg et al. 2006). In Finland, observations are recorded on a daily basis, whereas in Norway only the sum of daily values for the entire hunting season is reported. Because the hunting season is longer in Finland (9-12 weeks) than in Norway (four weeks), we used moose observations only from the first four weeks of hunting in Finland. Currently, more than 5,000 hunting clubs with about 100,000 hunters report 200,000-400,000 moose observations in Finland. The similar figures in Norway are about 60,000 hunters reporting approximately 200,000-250,000 observations each year (Rolandsen et al. 2004, Solberg et al. 2006).

Several indices of population structure are calculated from the observation data (e.g. Nygrén & Nygrén 1976, Nygrén & Pesonen 1993, Solberg & Heim 2002, 2006, Lavsund et al. 2003, Nygrén 2009). The most important are indices of population/regional productivity; ‘calves/100 adults’ (seen calves per adult in Norway), ‘calves/100 females’ (seen calves per female in Norway), ‘seen proportion of females with calves’ and ‘seen proportion of females with twins of all females with calf/calves’, as well as an index of population sex ratio; ‘seen females per male’. In addition, population density is indexed by the number of observations per unit hunting effort (i.e. in Norway, as moose seen per hunter day, and in Finland, as moose seen per team-hunting day). In Finland, population density is also indexed by the numbers of moose that the hunters estimate to be living on their hunting grounds after the hunting season. Despite the rather crude sampling procedure and several possible confounding variables (e.g. variation in weather, hunting skills, number of hunters and hunting methods), the observation data are found to provide good information on the temporal changes in recruitment rate, population density and adult sex ratio among populations (Ericsson & Wallin 1994, 1999, Sylvén 1995, Solberg & Sæther 1999, Solberg et al. 2002, Mysterud et al. 2007, Rönnegård et al. 2008, Bjørneraas et al. 2009), provided that the number of observations is relatively high (Ericsson & Wallin 1994, Sylvén 2000).

Body mass

As a measure of variation in calf body mass, we used the carcass mass of harvested calves during 1991-2006. The carcass mass of calves constitutes about 50% of their live body mass (Wallin et al. 1996). Because carcass mass of calves tends to increase during the hunting period (e.g. Solberg et al. 2004), we adjusted all masses to 15 October by using a regression of carcass mass on kill date within country. In total, we had body mass data from 16,841 calves (8,742 males and 8,099 females) in Norway and from 152,034 calves (80,059 males and 71,975 females) in Finland. The number of calves among the combinations of year and region (N = 384) ranged between 10 and 1,834, ensuring that the mean carcass mass constituted a relatively precise reflection of the mean calf body mass in the region (in Finland, game management district, and in Norway, the monitoring area; see Fig. 1). In the forthcoming text, carcass mass is referred to as body mass.

Calf production

Calf recruitment was measured as calving rate and twinning rate. Calving rate is the seen proportion of females with calf/calves during the hunting season, whereas twinning rate is the seen proportion of females with twins of all females with calf/calves. Because moose are harvested while observations are recorded, it is likely that variation in harvest pressure and structure may affect the estimated recruitment rates (Grøtan 2003). For instance, in regions with a high harvest of calves per harvested female, both recruitment rate indices are likely to decrease during the hunting season (e.g. Grøtan 2003, Moe et al. 2009). Such effects are regularly controlled for in Finland, where harvest and observation data are recorded on a daily basis, i.e. the pre-harvest recruitment rates can be estimated (Nygrén & Pesonen 1993, Nygrén 2009). Based on the Finnish data, the effect of variation in the harvesting of calves was found to have only minor effects on the variation in twinning rate, but may have substantial effects on the variation in calving rate. Particularly pronounced effects on the calving rate are expected in regions where the number of harvested calves per harvested female deviates much from the observed calves per female, and vice versa (Grøtan 2003).

Since the Norwegian observation data were aggregated over the entire hunting season, we were not able to adjust these observations for variation in the harvest of calves per female. We used unadjusted indices from both countries and statistically controlled varying calf harvesting by including a calf harvest index as a covariate in the analyses. The calf harvest index was calculated as the number of calves per female in the harvest divided by the number of calves per female in the observations. High values indicate that hunters shoot more calves per females than the ratio found in the region. Accordingly, the observed proportion of calves (calving/twinning rate) can be expected to be lower than in the pre-harvested population.

Statistical analyses

We first examined the variation in mean body mass of calves, as well as the calving rate and twinning rate within and among regions (for all variables, N = 384) in Finland and Norway using linear models. We expected lower body masses and recruitment rates in Norway than in Finland, and higher mean body mass of males than females. Likewise, we examined the relationship between calving rate and twinning rate within regions expecting a general positive relationship between recruitment indices over time, given the fact that both ratios are likely to be affected by female body condition (Sand 1996).

We then tested to what extent the variation in annual calving or twinning rate was related to variation in mean calf body mass. In these analyses, we used only the body mass of female calves as the body mass of females directly affects the body mass of the offspring (e.g. Keech et al. 2000). Because the calving rate and twinning rate are proportions, both variables were first arcsine transformed (arcsine (i0909-6396-18-3-304-ilm01.gif), where p is the ratio; Sokal & Rohlf 1995) to approximate a normal distribution.

To facilitate the interpretation of the spatio-temporal effects of calf body mass, we added both the annual mean body mass and the regional mean body mass in the models. These variables were calculated by first centring the annual mean body mass (Yt,r) relative to the mean body mass across all regions and years (μy):

i0909-6396-18-3-304-e01.gif
where t is year (1991-2006) and r is region. We then calculated the regional mean body mass (μr) as:
i0909-6396-18-3-304-e02.gif
where nr is number of years within region and the sum is over all years and the centred annual body mass (i0909-6396-18-3-304-ilm02.gif) as:
i0909-6396-18-3-304-e03.gif

By splitting body mass into regional mean (μr) and centred annual values (i0909-6396-18-3-304-ilm03.gif), we could specifically test to what extent the relationship between recruitment rates and body mass was due to among-region variation (spatial variation) or within-region variation (temporal variation), respectively (Singer 1998). The centred and mean body mass values were not correlated r < 0.01). Using centred values, the main effects are also biologically interpretable even when involved in interactions (Schielzeth 2010).

We analysed the variation in recruitment rates using linear mixed effect models with annual recruitment rate (twinning or calving rate) as a dependent variable and with regional mean (μr) and annual centred body mass (i0909-6396-18-3-304-ilm04.gif) of female calves as covariates. For the temporal effect, we tested the variation in recruitment rate (calving and twinning rate) against the variation in body mass in the current year (year t) as well as in the three previous years (year t-1, t-2 and t-3). We also included the observed sex ratio (females per male) and calf harvest index as covariates. We did not include country in these analyses as there was only a small overlap in annual mean body mass between countries. We expected lower recruitment rates in years and regions with a female-biased sex ratio and a high calf harvest index (more calves removed during hunting and hence decreasing calf observations during the harvest season), all else being equal.

To provide a baseline model against which we could compare the more complex models, we estimated the variance components in a model with only region included as a random factor. Two variance components were extracted from this model, representing the variances in twinning or calving rate that are attributed to variation within (residual variance) and among regions. Based on this model, we could 1) determine to what extent the annual twinning or calving rate varied among regions, and 2) estimate the fraction of explainable variation accounted for by the fixed effects in the more complex models. Following Singer (1998), we computed the proportion of explainable variation as:

i0909-6396-18-3-304-e04.gif
where VC1 and VC2 are the variance components in the baseline and the more complex model, respectively (see also e.g. Solberg et al. 2007 for the use of this method).

We considered the models with the lowest Akaike information criterion (AIC) score to be the most parsimonious and thus the ‘best’ model applied to the data (Burnham & Anderson 2002). We used AICc (corrected AIC) to account for the relatively high number of covariates compared to the number of years and regions. We started the model selection with addressing the overall model first and then all the possible nested models of the overall model. Altogether, 13 different models were tested, both for twinning rate and calving rate, from which models that differed in AICc by two or less (ΔAICc ≤ 2) were considered to be equally well supported by the data (Burnham & Anderson 2002). We also computed the Akaike weights to compare the relative performance of models, where the strength of evidence in favour of one model over another is the ratio in AIC weights between the two models (Burnham & Anderson 2002). Since we were comparing models with different fixed effect structure, we used Maximum-Likelihood (ML) for model selection (Singer 1998, Crawley 2002) and Restricted Maximum-Likelihood (REML) for parameter estimation. All the statistical analyses were performed using SPSS 19.0 for Windows.

Results

Variation on calf body mass

On average, calves were heavier in Finland than in Norway (see Figs. 1 and 2). The body mass (see Fig. 2) of male calves ranged from 25 to 120 kg in Finland (mean = 81 kg, SD = 7, N = 4,264) and from 44 to 96 kg in Norway (mean = 67 kg, SD = 7, N = 489). The corresponding figures for female calves (see Fig. 2) were 40-107 kg (mean = 77 kg, SD = 7, N = 4,245) and 28-91 kg (mean = 63 kg, SD = 7, N = 488), respectively. The variation among populations is shown in Fig. 1.

There was a strong positive relationship between annual mean body mass of male and female calves in both countries during our study period, 1991-2006 (Finland: r = 0.618, SE = 0.013, t = 48.348, P < 0.001; Norway: r = 0.677, SE = 0.039, t = 17.514, P < 0.001). Male calves were approximately 5% heavier than female calves (in Finland, the mean difference was 4.2 kg, and in Norway, it was 3.5 kg). In both countries, there was an overall negative trend in body mass (see Fig. 2; Finland: r = −0.223, SE = 0.017, t = −13.088, P < 0.001; Norway: r = −0.277, SE = 0.053, t = −5.243, P < 0.001) during our study period.

Variation in calving rate and twinning rate

Calf recruitment was generally higher in Finland than in Norway (see Figs. 1 and 2). On average, 59% of all adult females were seen with a calf/calves (i.e. calving rate) in Finland (mean = 0.59, SD = 0.07, N = 8,504) and 51% in Norway (mean = 0.51, SD = 0.82, N = 936). Similarly, 38% of the calf-rearing females in Finland recruited twin calves (i.e. twinning rate, mean = 0.38, SD = 0.09, N = 8,504) compared to only 28% in Norway (mean = 0.28, SD = 0.14, N = 936). The spatial distribution in population productivity can be seen in Fig. 1. Generally, there was a positive relationship between the annual twinning rate and calving rate within populations in both countries but with large variation among years (Finland: r = 0.821, N = 16; range among years: −0.031-0.622, N =15, P = 0.013-0.978; Norway: r = 0.725, N = 16; range among years: −0.003-0.937, N = 15, P = 0.001-0.994).

Spatio-temporal covariation between calf recruitment and female calf body mass

The best linear mixed effect models explaining the variation in twinning rate and calving rate included regional mean body mass (μr), centred annual mean body mass (i0909-6396-18-3-304-ilm05.gif) and the observed adult sex ratio. In addition, region was included as a random factor (Table 1). Also the calf harvest index was included in some of the best models, but were much less supported (the best models twinning rate and calving rate were 2.9 and 1.4 times more supported than the second best models, respectively). On average, calf recruitment was higher in regions and years with high mean body mass. We also found the effect of centred annual mean body mass (i0909-6396-18-3-304-ilm06.gif) to be stronger in regions with high mean body mass (positive regional mean body mass (μr)*centred annual mean body mass interaction (i0909-6396-18-3-304-ilm07.gif). However, this effect was significant for twinning rate but not for calving rate (Table 2), probably because the high calf harvest affects the ratio of calves per females observed during harvest season. The strongest effect of centred annual mean body mass was found with a time lag of one year (year t-1), indicating that the growth conditions during the previous summer were most important for the variation in recruitment rate (see Table 2).

The negative effect of adult sex ratio (females per male) was in accordance with expectations, and indicates that higher recruitment rates occur in years and populations with a higher proportion of males. Similarly, we found that higher recruitment rates were recorded in years and regions with a low off-take of calves relative to the calves per female in the region (i.e. low calf harvest index; see Tables 1 and 2).

From a total of 13 models (both for twinning and calving rate), we found four alternative models for twinning rate and calving rate within ΔAIC ≤ 6 (see Table 1). Alternative models included or excluded the calf harvest index and regional body mass (μr) in year t-1 as a random slope (see Table 1). Hence, the general relationship between recruitment rates and mean calf body mass seemed relatively consistent among the models.

The two best models for twinning rate accounted for about 26% of the variation within regions and about 49% of the variation among regions (see Table 1). The best model for calving rate accounted for about 38% and 63% of the within- and among-region variation, respectively. Adding the calf harvest index increased the explanatory power of the model to 65% of the among-region variance while the within-region variance remained the same 38% (see Table 1). Hence, a larger proportion of the explainable variation in recruitment rates was accounted for among than within regions.

Discussion

Our results show that there is large variation in calf body mass and recruitment rates within (i.e. temporal variation) and among regions (see Figs. 1 and 2), and that the population condition, measured as recruitment rate and the mean body mass of calves, is generally higher in Finland than in Norway (see Figs. 1 and 2). Populations with larger calves are also more productive, as was expected based on the positive body size and fecundity relationship found in moose (Sand 1996, Garel et al. 2009). The annual variation in recruitment rates was best explained by a time-delayed effect of calf body mass, indicating that moose recruitment rates in Fennoscandia depends more on accumulated reserves than on resources available in the year of reproduction. We also found higher observed recruitment rates in regions and years with more even adult sex ratio and low relative harvest rates of calves. This suggests that variation in recruitment rate is not only a matter of nutritional status, but may as well be a product of varying hunting regulations and harvest structure.

The positive relationships between recruitment indices and calf body mass in moose support the notion that body condition and reproduction will be lower in populations close to carrying capacity (Kie et al. 2003, Stewart et al. 2005). In Finland, moose densities are generally lower than in Norway (Lavsund et al. 2003) and the moose ranges are found at lower altitudes. Hence, it can be assumed that the primary production on moose land is also higher in Finland than in Norway (Zheng et al. 2004). Following particularly high densities (> 2 moose/km2), many populations in southern Norway have experienced significant decline in calf body mass and recruitment rates during the last two decades, most likely due to density-dependent food limitation (Lavsund et al. 2003, Solberg et al. 2006). Further north, densities are lower and the feeding conditions better, with the result that the recruitment rates and mean calf body mass increases from south to north (Solberg et al. 2006; see Fig. 1). In Finland, the moose densities are kept relatively low by a strict harvesting regime (Lavsund et al. 2003, Nygrén 2009), providing good condition for body growth and reproduction. Probably for the same reasons, we also see much less variation in body mass and recruitment rates between Finnish than Norwegian populations. However, although the differences in life history suggest that moose densities in Finland are low compared to the carrying capacity, the effects of varying density and habitat productivity on moose population performance in Fennoscandia are still in need of closer studies.

The positive relationship between calf body mass and recruitment indices are also in accordance with previous findings showing that larger moose calves tend to become larger adults (Solberg et al. 2008), which mature at a younger age (Sæther & Haagenrud 1985, Sand 1996), start twin production earlier (Sand 1996, Solberg et al. 2008) and produce larger calves on average (Solberg et al. 2007). Large size is often correlated with high phenotypic quality (Gaillard et al. 2000a, Hewison & Gaillard 2001), and large mothers therefore can allocate more resources to offspring than smaller mothers can (e.g. Loison & Strand 2005). However, because maternal resources can be distributed on varying number of young (Stearns 1992), there is not necessarily always a positive relationship between the size of the mother and her offspring.

In moose, larger females do not seem to trade off calf size for calf number; at least not in populations living under good nutritional condition (Solberg et al. 2007). However, such trade-offs can potentially be more expressed in populations living under food limitation, as may be the case for some regions in Norway. In the populations in Oppland and Nord-Trøndelag, for instance, there are large differences in twinning rate despite hardly any differences in mean calf body mass (see Fig. 1). Possibly, this is because fewer resources available per female in Oppland (low twinning rate) are invested in a single calf, whereas more resources available to females in Nord-Trøndelag (high twinning rate) are diverted on twins. Ultimately, such differences in the trade-offs between size and fecundity may be an evolutionary product of different mortality patterns (e.g. following different environmental conditions and predation; Proactor et al. 2007). The alternative, i.e. that low calf mass and high twinning rates in Nord-Trøndelag are due to females reflecting previously poor living conditions (e.g. Monteith et al. 2009), is less likely as the population in Nord-Trøndelag has been increasing from low density during the last three decades. Likewise, we find it unlikely that the poor fit between regional mean calf body mass (μr) and calf recruitment can be explained by varying perinatal mortality among populations. Indeed, in Oppland, the ovulation rates (proportion ovulating and twin ovulating females) recorded from ovaries, were substantially lower than in Nord-Trøndelag (Solberg et al. 2006, Garel et al. 2009). This suggests that differences in recruitment rates are already settled prior to any calf loss, possibly due to a varying ovulation-to-body-mass relationship among populations (Garel et al. 2009).

Interestingly, the variation in adult sex ratio (females per male) also seemed to affect fecundity as the recruitment rates for a given calf body mass were found to be lower in regions and years with low proportion of adult males (see Table 2). The age and proportion of adult males have decreased substantially in Fennoscandian populations since the introduction of sex and age-specific harvesting in the early 1970s (Nygrén 1987, Lavsund et al. 2003, Milner et al. 2007). Like in other species (e.g. Noyes et al. 1996, Mysterud et al. 2002), such distortions of the population structure can affect reproduction if the number and/or age of males become inadequate to ensure impregnation of all females and/or affects the timing of parturition and hence the size of calves during autumn (e.g. Solberg et al. 2002, Sæther et al. 2004). However, while such a mechanism can explain why we observed a negative effect on the calving rate, it is less clear why the variation in twinning rate is also related to adult sex ratio.

The rate of population increase calculated from life table data is probably the most concise measure of nutritional status or vigour of game populations (Caughley 1977). However, for most species, age-specific fecundity and natural mortality rates are hard to obtain, and for that reason various indices of population nutritional status are often used for monitoring purposes (Hanks 1981). For adaptive management of large herbivores, temporal changes in both population and habitat features, and their interaction, is needed (Morellet et al. 2007, 2011). The twinning rate at birth is one such index that has been used for assessment of moose population nutritional status in North America (Franzmann & Schwartz 1985, Boer 1992, Boertje et al. 2007). Because twins are mainly produced by moose with the highest body mass, age and condition (e.g. Sand 1996, Solberg et al. 2008), the twinning rate at birth is probably a good reflection of the nutritional status in an area. Indeed, the population twinning rate at birth is often found to relate positively to other measurements of population nutritional status (e.g. Boer 1992, Boertje et al. 2007).

In Fennoscandia, aerial surveys of moose reproduction are seldom used for population monitoring, but instead, moose managers utilise recruitment data derived from moose observations by hunters (Morellet et al. 2011). Given that the effect of predation and harvesting is low or can be accounted for, we believe that such data can be used to generate precise indices of population nutritional status. Indeed, as such observations are provided by most hunting teams in Fennoscandia, indices based on a large number of observations can be generated at relatively high spatial resolution. Moreover, because the abundance of wolves Canis lupus and bears Ursus arctos is still quite modest in Fennoscandia (Swenson et al. 1995, Wabakken et al. 2001, Aspi et al. 2006, Kojola & Heikkinen 2006, Heikkinen et al. 2009, Nygrén 2009), predation is unlikely to have a strong impact on the calving and twinning rates from spring to autumn. An exception is the eastern Finland (Pohjois-Karjala (PK) and Kainuu (KA) in Fig. 1) where large predators are abundant enough to have regional effects on the population productivity of moose (Nygrén 2009).

Biased harvesting of calves and adults may have a stronger effect than predation on the recruitment rates, particularly if the ratio of calves to females in the harvest deviates much from the calf to female ratio in the population (Grøtan 2003). In Fennoscandia, much focus has been directed at harvesting calves and protecting productive females (Lavsund et al. 2003, Nygrén 2009), leading to decreasing calving rates during the hunting season in most populations (Moe et al. 2009). Hence, if the recruitment rates are based on the observations accumulated over the entire hunting season, the recruitment rates are likely to be underestimates. This is in accordance with our results showing that the calf harvest index had a negative effect on the variation in calving rate and twinning rate (see Table 2). This effect seemed to be particularly strong in Norway, probably because of the larger variation in the calf harvest index among Norwegian (range: 0.71-1.62) than Finnish (range: 1.34-1.65) populations.

Based on a simulation study, Grøtan (2003) showed that the variation in calving rate is much more sensitive to variation in the calf harvest rate than is the twinning rate as long as twin calves do not experience substantially higher harvest rates than singletons. This was also supported by our preliminary analyses showing that the twinning rate in Finnish populations did not change much during the hunting season despite substantially higher harvest rates of calves than adult females. The same is likely to be true if predation rates are high and twin calves are not more vulnerable than single calves. In contrast, the observed calving rate may be substantially lower at the end of the hunting season in populations where calves experience high relative harvest rates (or predation rates). Accordingly, the calving rate will probably provide a less precise reflection of the nutritional conditions of a population than do the twinning rate if the variation in calf harvest rate is large and not accounted for. Based on the experience from Finland, such effects of harvesting can best be accounted for by recording the moose observations on a daily basis, from which the pre-harvest calving and twinning rate can be estimated (Nygrén & Pesonen 1993, Nygrén 2009, Morellet et al. 2011).

While the observed twinning rate seemed to be a good index for ranking populations, it showed relatively small variation among years within regions (see Table 1). This is probably due to the resource use tactic found in moose. Although moose females partly rely on ingested food to feed their calves, they are primarily capital breeders (i.e. fecundity is a product of accumulated reserves; Stearns 1992, Jönsson 1997, Mysterud et al. 2005). Accordingly, the number and size of calves is not only a product of the feeding conditions during summer, but also the result of resources accumulated by the mother throughout her life (Keech et al. 1999, Monteith et al. 2009). Such a strategy can partly act as a buffer against short-term environmental effects on the reproductive output making the twinning rate weaker than mean calf body mass in reflecting annual variation in nutritional conditions. In this context, it is also important to note that fluctuations in female age structure (e.g. due to harvesting) may generate variation in twinning rate because older (and larger) moose are more likely to produce twins (Nygrén 1983, 2009). Thus, within a given period, the twinning rate may be high because of high nutritional conditions, high average age of females or both.

For monitoring purposes, we believe that the mean calf body mass is a better index of the short term variation in living conditions than is the twinning rate. Juveniles in ungulates are usually among the first to suffer when the conditions decline (e.g. Gaillard et al. 2000b), and the body mass of calves are found to be closely related to environmental conditions associated with food quality and quantity (e.g. Gaillard et al. 1996 on population density in roe deer, Sæther 1997 in general, Herfindal et al. 2006 for moose), as well as maternal effects (e.g. Keech et al. 1999, Monteith et al. 2009). Moreover, while the mean calf body mass reflects the body condition of one cohort, the observed recruitment rate is a composite of several cohorts that have experienced varying conditions for body growth during life. The observed positive relationship between recruitment rates and calf body mass may therefore mainly reflect that females and calves share similar conditions for body growth and fat accumulation during the summer, which in turn affect female ovulation rates and reproductive rates the next summer. This is in accordance with the result that both the calving rate and twinning rate were best related to variation in mean calf body mass in the previous year and not the current year.

From this perspective, we suggest that the twinning rate is best suited for ranking populations according to nutritional status, as this index is closely related to fitness and is relatively insensitive to variation in perinatal and harvest mortality (see Boertje et al. 2007 for similar conclusions). However, variation in calf body mass may better reflect temporal variation in living conditions. This is because early body growth is sensitive to variation in food availability (and quality) and because body mass responds more instantaneously than recruitment indices to adverse conditions. The earlier response in calf body mass than in recruitment rates also suggests that changes in mean body mass can be used by managers to foresee subsequent changes in recruitment rate. Thus, in absence of more direct censuses prior to the hunting season, the previous year's mean calf body mass can to some extent be used for fine-tuning hunting quotas and to reduce the time lag often found between changes in resource abundance and management decisions (Fryxell et al. 2010).

Acknowledgements

we thank all the Finnish and Norwegian hunters involved in collecting the data, as well as Kaarlo Nygrén, Mauri Pesonen, Riitta Tykkyläinen and Maija Wallén for providing and handling the data in Finland. We also thank Esa Koskela and Jan Lindström for their valuable comments on the manuscript. We are indebted to the University of Joensuu, Emil Aaltonen Foundation and Suomen Riistanhoito-Säätiö for funding (Raisa Tiilikainen), as well as ‘Norsk Institutt for Naturforskning’ (NINA) and the Norwegian Research Council (Miljø 2015).

References

1.

S.D Albon T.H Clutton-Brock& F.E Guinness 1987: Early development and population dynamics in red deer II. Density-independent effects and cohort variation. - Journal of Animal Ecology 56: 69–81. Google Scholar

2.

J Aspi E Roininen M Ruokonen I Kojola& C Vila 2006: Genetic diversity, population structure, effective population size, and demographic history of the Finnish wolf population. - Molecular Ecology 15: 1561–1576. Google Scholar

3.

A Beckerman T.G Benton E Ranta V Kaitala& P Lundberg 2002: Population dynamic consequences of delayed life-history effects. - Trends in Ecology & Evolution 17: 263–269. Google Scholar

4.

C.H Bérubé M Festa-Bianchet& J.T Jorgenson 1999: Individual differences, longevity, and reproductive senescence in bighorn ewes. - Ecology 80: 2555–2565. Google Scholar

5.

A.H Boer 1992: Fecundity of North American moose (Alces alces): a review. - Alces 28 (Suppl. 1): 1–10. Google Scholar

6.

R.D Boertje A.K Kalin C.T Seaton M.A Keech D.D Young B.W Dale L.G Adams& A.R Aderman 2007: Ranking Alaska moose nutrition: signals to begin liberal antlerless harvest. - Journal of Wildlife Management 71: 1494–1506. Google Scholar

7.

K Bjørneraas E.J Solberg I Herfindal& B-E Sæther 2009: Large-scale spatiotemporal variation in calf sex ratio in moose (Alces alces): an effect of density-dependent decrease in maternal condition? - Canadian Journal of Zoology 87: 346–355. Google Scholar

8.

K.P Burnham& D.R Anderson 2002: Model selection and multimodel inference: a practical information-theoretic approach. Springer Verlag, New York, New York, USA, pp. Google Scholar

9.

G Caughley 1977: Analysis of vertebrate populations. Wiley, London, UK, pp. Google Scholar

10.

M.J Crawley 2002: Statistical computing: an introduction to data analysis using S-Plus. Wiley, West Sussex, UK, pp. Google Scholar

11.

G Ericsson& K Wallin 1994: Antal älgar som ses - bara en fråga om hur många som finns? Att observera älg - en fråga om täthet, rörelser och synbarhet. Swedish University of Agricultural Sciences, Department of Animal Ecology, Umeå, Sweden, pp. (In Swedish). Google Scholar

12.

G Ericsson& K Wallin 1999: Hunter observations as an index of moose Alces alces population parameters. - Wildlife Biology 5(3): 177–185. Google Scholar

13.

M Festa-Bianchet J Jorgenson& D Réale 2000: Early development, adult mass, and reproductive success in bighorn sheep. - Behavioral Ecology 11: 633–639. Google Scholar

14.

A.W Franzmann& C.C Schwartz 1985: Moose twinning rates: a possible population condition assessment. - Journal of Wildlife Management 49: 394–396. Google Scholar

15.

J.M Fryxell C Packer K McCann E.J Solberg& B-E Sæther 2010: Resource management cycles and the sustainability of harvested wildlife populations. - Science 328: 903–906. Google Scholar

16.

J-M Gaillard D Delorme J-M Boutin G Van Laere& B Boisaubert 1996: Body mass of roe deer fawns during winter in 2 contrasting populations. - Journal of Wildlife Management 60: 29–36. Google Scholar

17.

J-M Gaillard M Festa-Bianchet D Delorme& D Jorgenson 2000a: Body mass and individual fitness in female ungulates: bigger is not always better. - Proceedings of the Royal Society, Series B 267: 471–477. Google Scholar

18.

J-M Gaillard M Festa-Bianchet D Delorme& D Jorgenson 2000b: Temporal variation in fitness components and population dynamics of large herbivores. - Annual Review of Ecology and Systematics 31: 367–393. Google Scholar

19.

J-M Gaillard A.J Sempere J-M Boutin G.V Laere& B Boisaubert 1992: Effects of age and body weight on the proportion of females breeding in a population of roe deer (Capreolus capreolus). - Canadian Journal of Zoology 70: 1541–1545. Google Scholar

20.

M Garel E.J Solberg B-E Sæther V Grøtan J Tufto& M Heim 2009: Age, size and spatio-temporal variation in ovulation patterns of a seasonal breeder, the Norwegian moose (Alces alces). - American Naturalist 173: 89–104. Google Scholar

21.

J.R Ginsberg& E.J Milner-Gulland 1994: Sex-biased harvesting and population dynamics in ungulates: implications for conservation and sustainable use. - Conservation Biology 8: 157–166. Google Scholar

22.

V Grøtan 2003: Large scale synchronization of moose (Alces alces) population dynamics trough climate and harvest. - M.Sc. thesis, Department of Biology, Norwegian University of Science and Technology, Trondheim, Norway, pp. Google Scholar

23.

J Hanks 1981: Characterization of Population Condition. In: C.W Fowler& T.D Smith (Eds.); Dynamics of Large Mammal Populations. John Wiley and Sons, New York, New York, USA, pp. 47–73. Google Scholar

24.

I Herfindal B-E Sæther E.J Solberg R Andersen& K-A Høgda 2006: Population characteristics predict responses in moose body mass to temporal variation in the environment. - Journal of Animal Ecology 75: 1110–1118. Google Scholar

25.

S Heikkinen H Hiltunen I Kojola& E Määttä 2009: Suurpetojen lukumäärä ja lisääntyminen vuonna 2008. In: M Wikman (Ed.); Riistakannat 2008: Riistaseurantojen tulokset (In Finnish with an English summary: Monitoring game abundance in Finland in 2009). Riista- ja kalatalous - Selvityksiä 18, pp. Google Scholar

26.

A.J.M Hewison& J-M Gaillard 2001: Phenotypic quality and senescence affect different components of reproductive output in roe deer. - Journal of Animal Ecology 70: 600–608. Google Scholar

27.

K.I Jönsson 1997: Capital and income breeding as alternative tactics of resource use in reproduction. - Oikos 78: 57–66. Google Scholar

28.

M.A Keech R.D Boertje R.T Bowyer& B.W Dale 1999: Effects of birth weight on growth of young moose: do low-weight neonates compensate? - Alces 35: 51–57. Google Scholar

29.

M.A Keech R.T Bowyer J.M.V Hoef R.D Boertje B.W Dale& T.R Stephenson 2000: Life-history consequences of maternal condition in Alaskan moose. - Journal of Wildlife Management 64: 450–462. Google Scholar

30.

J.G Kie R.T Bowyer& K.M Stewart 2003. Ungulates in western coniferous forests: habitat relationships, population dynamics, and ecosystem processes. In: C.J Zabel& R.G Anthony (Eds.); Mammal community dynamics. Management and conservation in the coniferous forests of Western North America. Cambridge University Press, Cambridge, UK, pp. Google Scholar

31.

I Kojola& S Heikkinen 2006: The structure of the expanded brown bear population at the edge of the Finnish range. - Annales Zoologici Fennici 43: 258–262. Google Scholar

32.

T Koskela& T Nygrén 2002: Hirvenmetsästysseurueet Suomessa vuonna 1999. (In Finnish with an English summary: Moose hunting clubs in Finland 1999). - Suomen Riista 48: 65–79. Google Scholar

33.

S Lavsund T Nygrén& E.J Solberg 2003: Status of moose populations and challenges to moose management in Scandinavia. - Alces 39: 109–130. Google Scholar

34.

A Loison& O Strand 2005: Allometry and variability of resource allocation to reproduction in a wild reindeer population. - Behavioural Ecology 16: 624–633. Google Scholar

35.

J.M Milner E.B Nilsen& H.P Andreassen 2007: Demographic side effects of selective hunting in ungulates and carnivores. - Conservation Biology 21: 36–47. Google Scholar

36.

T Moe E.J Solberg I Herfindal B-E Sæther K Bjørneraas& M Heim 2009: Sex ratio variation in harvested moose (Alces alces) calves: does it reflect population calf sex ratio or selective hunting? - European Journal of Wildlife Research 55: 217–226. Google Scholar

37.

K.L Monteith L.E Schmitz J.A Jenks J.A Delger& R.T Bowyer 2009: Growth of male white-tailed deer: consequences of maternal effects. - Journal of Mammalogy 90: 651–660. Google Scholar

38.

N Morellet J-M Gaillard A.J.M Hewison P Ballon Y Boscardin P Duncan F Klein& D Maillard 2007: Indicators of ecological change: new tools for managing populations of large herbivores. - Journal of Applied Ecology 44: 634–643. Google Scholar

39.

N Morellet F Klein E Solberg& R Andersen 2011: The census and management of populations of ungulates in Europe. In: R Putman M Apollonio& R Andersen (Eds.); Ungulate Management in Europe. Problems and Practices. Cambridge University Press, Cambridge, UK, pp. 106–143. Google Scholar

40.

A Mysterud T Coulson& N.C Stenseth 2002: The role of males in the dynamics of ungulate populations. - Journal of Animal Ecology 71: 907–915. Google Scholar

41.

A Mysterud E.L Meisingset V Veiberg R Langvatn E.J Solberg L.E Loe& N.C Stenseth 2007: Monitoring population size of red deer: an evaluation of two types of census data from Norway. - Wildlife Biology 13(3): 285–298. Google Scholar

42.

A Mysterud E.J Solberg& N.G Yoccoz 2005: Ageing and reproductive effort in male moose under variable levels of intrasexual competition. - Journal of Animal Ecology 74: 742–754. Google Scholar

43.

J.H Noyes B.K Johnson L.D Bryant S.L Findholt& J.W Thomas 1996: Effects of bull age on conception dates and pregnancy rates of cow elk. - Journal of Wildlife Management 60: 508–517. Google Scholar

44.

T Nygrén 1983: The relationship between reproduction rate and age structure, sex ratio and density in the Finnish moose population. Proceedings from XVI Congress of the International Union of Game Biologists, Vysoké Tatry, Štrbské Pleso, ČSSR, pp. 29–42. Google Scholar

45.

T Nygrén 1984: Hirvikannan arviointi ja verotuksen suunnittelu Suomessa (In Finnish with an English summary: Moose population census and planning of cropping in Finland). - Suomen Riista 31: 74–82. Google Scholar

46.

T Nygrén 1987: The history of moose in Finland. - Swedish Wildlife Research (Suppl. 1): 49-54. Google Scholar

47.

T Nygrén 2003: The potential for multiple fecundity of moose in Finland. - Alces 39: 89–107. Google Scholar

48.

T Nygrén 2009: Suomen hirvikannan säätely - biologiaa ja luonnonvarapolitiikkaa. - PhD dissertation, University of Joensuu, Joensuu, Finland, pp. (In Finnish). Google Scholar

49.

K Nygrén& T Nygrén 1976: Hirvi ja hirvenmetsästys Suomessa. - Riistantutkimusosaston Tiedonantoja 2: 1–33. (In Finnish). Google Scholar

50.

T Nygrén& M Pesonen 1993: The moose (Alces alces) population and methods of moose management in Finland, 1975-1989. - Finnish Game Research 48: 46–53. Google Scholar

51.

T Nygrén J Pusenius R Tiilikainen& J Korpelainen 2007: Moose antler type polymorphism: age and weight dependent phenotypes and phenotype frequencies in space and time. - Annales Zoologici Fennici 44: 445–461. Google Scholar

52.

G Proaktor T Coulson& E.J Milner-Gulland 2007: Evolutionary responses to harvesting in ungulates. - Journal of Animal Ecology 76: 669–678. Google Scholar

53.

C.R Rolandsen E.J Solberg& V Grøtan 2004: ‘Sett elg’-materialet i Norge 1984-2002. Hjorteviltet: 6–13. (In Norwegian). Google Scholar

54.

C Rolandsen E.J Solberg M Heim F Holmstrøm M.I Solem& B-E Sæther 2008: Accuracy and repeatability of moose (Alces alces) age as estimated from dental cement layers. - European Journal of Wildlife Research 54: 6–14. Google Scholar

55.

L Rönnegård H Sand H Andrén J Månsson& A Pehrson 2008: Evaluation of four methods used to estimate population density of moose Alces alces. - Wildlife Biology 14(3): 358–371. Google Scholar

56.

V Ruusila& M Pesonen 2004: Interspecific cooperation in human (Homo sapiens) hunting: the benefits of a barking dog (Canis familiaris). - Annales Zoologici Fennici 41: 545–549. Google Scholar

57.

H Sand 1996: Life history patterns in female moose (Alces alces): the relationship between age, body size, fecundity and environmental conditions. - Oecologia 106: 212–220. Google Scholar

58.

B-E Sæther 1997: Environmental stochasticity and population dynamics of large herbivores: a search for mechanisms. - Trends in Ecology & Evolution 12: 143–149. Google Scholar

59.

B-E Sæther& H Haagenrud 1985: Life history of the moose Alces alces: relationship between growth and reproduction. - Holarctic Ecology 8: 100–105. Google Scholar

60.

B-E Sæther E.J Solberg M Heim J.E Stacy K.S Jakobsen& R Olstad 2004: Offspring sex ratio in moose Alces alces in relation to paternal age: an experiment. - Wildlife Biology 10(1): 51–57. Google Scholar

61.

H Schielzeth 2010: Simple means to improve the interpretability of regression coefficients. - Methods in Ecology and Evolution 1: 103–113. Google Scholar

62.

J.D Singer 1998: Using SAS PROC MIXED to fit multilevel models, hierarchical models, and individual growth models. - Journal of Educational and Behavioral Statistics 24: 323–355. Google Scholar

63.

R.R Sokal& F.J Rohlf 1995: Biometry. The principles and practice of statistics in biological research. WH Freeman and Company, New York, New York, USA, pp. Google Scholar

64.

E.J Solberg M Garel M Heim V Grøtan& B-E Sæther 2008: Lack of compensatory body growth in a high performance moose Alces alces population. - Oecologia 158: 485–498. Google Scholar

65.

E.J Solberg& M Heim 2002: Monitoring moose in Norway: See them, shoot them, measure them and eat them. - Moose and Deer 2002: 16–19. Google Scholar

66.

E.J Solberg& M Heim 2006: Egenevaluering av overvåkingsprogrammet for elg. (In Norwegian with an English summary: Evaluation of the Norwegian monitoring program for moose). - Norsk institutt for naturforskning, Trondheim, Norway, pp. Google Scholar

67.

E.J Solberg M Heim V Grøtan B-E Sæther& M Garel 2007: Annual variation in maternal age and calving date generate cohort effects in moose (Alces alces) body mass. - Oecologia 154: 259–271. Google Scholar

68.

E.J Solberg A Loison J-M Gaillard& M Heim 2004: Lasting effects of conditions at birth on moose body mass. - Ecography 27: 677–687. Google Scholar

69.

E.J Solberg A Loison T.H Ringsby B-E Sæther& M Heim 2002: Biased adult sex ratio can affect fecundity in primiparous moose Alces alces. - Wildlife Biology 8(2): 117–128. Google Scholar

70.

E.J Solberg C.M Rolandsen M Heim V Grøtan B-E Sæther E.B Nilsen G Austheim& I Herfindal 2006: Elgen i Norge sett med jegerøyne: En analyse av jaktmaterialet fra overvåkningsprogrammet for elg og det samlede sett elg-materialet for perioden 1966-2004. (In Norwegian with an English summary: Moose in Norway - An analysis of material collected by moose hunters 1966-2004). Norsk institutt for naturforskning, Trondheim, Norway, pp. Google Scholar

71.

E.J Solberg& B-E Sæther 1999: Hunter observations of moose Alces alces as a management tool. - Wildlife Biology 5(2): 107–118. Google Scholar

72.

E.J Solberg B-E Sæther O Strand& A Loison 1999: Dynamics of a harvested moose population in a variable environment. - Journal of Animal Ecology 68: 186–204. Google Scholar

73.

S.C Stearns 1992: Evolution of life histories. Oxford University Press, New York, New York, USA, pp. Google Scholar

74.

K.M Stewart R.T Bowyer B.L Dick B.K Johnson& J.G Kie 2005: Density-dependent effects on physical condition and reproduction in North American elk: an experimental test. - Oecologia 143: 85–93. Google Scholar

75.

J.E Swenson P Wabakken F Sandegren A Bjarvall R Franzén& A Soderberg 1995: The near extinction and recovery of brown bears in Scandinavia in relation to the bear management policies of Norway and Sweden. - Wildlife Biology 1(1): 11–25. Google Scholar

76.

S Sylvén 1995: Moose harvest strategy to maximise yield value for multiple goal management. - A simulation study. - Agricultural Systems 49: 277–298. Google Scholar

77.

S Sylvén 2000: Effects of scale on hunter moose Alces alces observation rate. - Wildlife Biology 6(3): 157–165. Google Scholar

78.

P Wabakken H Sand O Liberg& A Bjarvall 2001: The recovery, distribution, and population dynamics of wolves on the Scandinavian peninsula, 1978-1998. - Canadian Journal of Zoology 79: 710–725. Google Scholar

79.

K Wallin G Cederlund& Å Pehrson 1996: Predicting body mass from chest circumference in moose Alces alces. - Wildlife Biology 2(1): 53–58. Google Scholar

80.

D Zhang S Price& T Hame 2004: Estimating net primary production of boreal forests in Finland and Sweden from field data and remote sensing. - Journal of Vegetation Science 15: 161–170. Google Scholar

Notes

[1] Edited by Associate Editor: Dean Anderson

Figure 1.

Our study areas in Norway (Vest-Agder, Aust-Agder, Vestfold/Telemark, Hedmark, Oppland, Nord-Trøndelag, Nordland, Beiarn and Troms) and Finland (EH = Etelä-Häme, ES = Etelä-Savo, KA = Kainuu, KS = Keski-Suomi, KY = Kymi, LA = Lappi, PH = Pohjois-Häme, PK = Pohjois-Karjala, PO = Pohjanmaa, PS = Pohjois-Savo, RP = Ruotsinkielinen Pohjanmaa, SK = Satakunta, VS = Varsinais-Suomi and UM = Uusimaa). The numerals I-V indicate the category of mean body mass of female calves and the twinning rate in the region, category I being the lowest and V the highest values, respectively.

i0909-6396-18-3-304-f01.tif

Figure 2.

Annual variation in A) mean body mass of calves (males black, females grey), B) twinning rate and C) calving rate in moose populations in Finland (black) and Norway (grey) during 1991-2006. Calf body mass is carcass mass adjusted for variation in kill date.

i0909-6396-18-3-304-f02.tif

Table 1.

The best candidate models explaining the spatio-temporal variation in twinning rate and calving rate (both arcsine transformed) in Fennoscandian moose using linear mixed effect models and with region as random factor (intercept). X indicates variables included in the model. BM is the baseline model with only region included as random factor (intercept). Models where centred annual mean body mass in year t-1 was included as random slope are marked with an underlined X. ΔAICc refers to the difference in AICc values between the best model and the candidate model. All candidate models with ΔAICc ≤ 6 are presented in the table. The best twinning rate model (model 1) explained 26% of the within-region (i.e. temporal) and 49% of the among-region variance compared to the baseline model (BM). The corresponding numbers for calving rate were 38% and 63% (model 1), respectively. Regional mean body mass (μr) = rbm and centred annual mean body mass (i0909-6396-18-3-304-ilm08.gif) in year t-1 = cbmt-1.

i0909-6396-18-3-304-t01.eps

Table 2.

Parameter estimates and test statistics for the best models in Table 1 (i.e. models with the lowest AICc).

i0909-6396-18-3-304-t02.eps
Raisa Tiilikainen, Erling Johan Solberg, Tuire Nygrén, and Jyrki Pusenius "Spatio-temporal relationship between calf body mass and population productivity in Fennoscandian moose Alces alces," Wildlife Biology 18(3), 304-317, (1 September 2012). https://doi.org/10.2981/10-116
Received: 3 November 2010; Accepted: 1 February 2012; Published: 1 September 2012
JOURNAL ARTICLE
14 PAGES


SHARE
ARTICLE IMPACT
Back to Top