Statistical population reconstruction of moose (Alces alces) in northeastern Minnesota using integrated population models

Abstract

Given recent and abrupt declines in the abundance of moose (Alces alces) throughout parts of Minnesota and elsewhere in North America, accurately estimating statewide population trends and demographic parameters is a high priority for their continued management and conservation. Statistical population reconstruction using integrated population models provides a flexible framework for combining information from multiple studies to produce robust estimates of population abundance, recruitment, and survival. We used this framework to combine aerial survey data and survival data from telemetry studies to recreate trends and demographics of moose in northeastern Minnesota, USA, from 2005 to 2020. Statistical population reconstruction confirmed the sharp decline in abundance from an estimated 7,841 (90% CI = 6,702–8,933) in 2009 to 3,386 (90% CI = 2,681–4,243) animals in 2013, but also indicated that abundance has remained relatively stable since then, except for a slight decline to 3,163 (90% CI = 2,403–3,718) in 2020. Subsequent stochastic projection of the population from 2021 to 2030 suggests that this modest decline will continue for the next 10 years. Both annual adult survival and per-capita recruitment (number of calves that survived to 1 year per adult female alive during the previous year) decreased substantially in years 2005 and 2019, from 0.902 (SE = 0.043) to 0.689 (SE = 0.061) and from 0.386 (SE = 0.030) to 0.303 (SE = 0.051), respectively. Sensitivity analysis revealed that moose abundance was more sensitive to fluctuations in adult survival than recruitment; thus, we conclude that the steep decline in 2013 was driven primarily by decreasing adult survival. Our analysis demonstrates the potential utility of using statistical population reconstruction to monitor moose population trends and to identify population declines more quickly. Future studies should focus on providing better estimates of per-capita recruitment, using pregnancy rates and calf survival, which can then be incorporated into reconstruction models to help improve estimates of population change through time.

Introduction

Effective management and conservation of wildlife species requires an accurate understanding of population abundance, recruitment, survival, and age- and sex-ratios, and how these parameters change over time and in response to various extrinsic factors, such as hunting and habitat alteration. Unfortunately, accurately estimating abundance and demographic parameters is challenging, because direct monitoring of animals is often costly and impractical, particularly in densely forested regions or for animals that occur at low densities. Given these difficulties, most abundance estimates have relied on methods that are limited to small geographical areas or sample sizes, including track surveys [1], analysis of camera traps [2], and telemetry data [3]. Each of these methods by themselves do not provide a cost-effective means of estimating abundance and other demographic parameters across larger spatial scales at which most management occurs.

Statistical population reconstruction using integrated population models (IPMs) has emerged as a flexible framework for combining information from multiple studies using various, disparate datasets (e.g., aerial surveys, radio-collared individuals, age-at-harvest), and even from different parts of a state or region, to provide a more robust and cost-effective means of estimating species abundance and demographics across large spatial scales [4, 5]. This method simultaneously estimates multiple demographic parameters (e.g., annual abundance, recruitment, and survival) and their uncertainties throughout time, and can be used to provide separate estimates for different sexes and age classes. Such models have previously been used to estimate abundance and trends of wildlife species, such as American marten (Martes americana), black bears (Ursus americanus), and mountain lions (Puma concolor) [6–8].

Accurately estimating the abundance and trajectory of the moose (Alces alces) population in northeastern Minnesota (MN) is of current interest due to a recent and abrupt decline that was detected via aerial surveys between 2010 and 2013 [9]. At its nadir in 2013, this population estimate was 69% lower than when at its peak in 2006 (2,760 versus 8,840), but it appeared to have stabilized during 2012–2020 as estimated by aerial surveys [9, 10]. A study of demographics of the northeastern population in 2002–2008 predicted a slow reduction in numbers (long-term stochastic annual growth rate [λ] of 0.85,) with modeled adult and calf survival rates of 0.74–0.85 and 0.24–0.56, respectively [11]. However, the abrupt decline in northeastern MN was not detected by the annual aerial surveys until 2010 [11–13], which illustrated that demographic modeling may reveal population trajectories before they are reflected in total population estimates by aerial survey.

In response to the rapid decline of moose in the northeastern population, the MN Department of Natural Resources (MNDNR), Grand Portage Band of Lake Superior Chippewa, and Voyageurs National Park all independently initiated studies of adult and calf survival and cause-specific mortality (Fig 1). These studies built upon previous research [11, 14], but aimed to better understand causes of mortality [15, 16]. The more recent research employed global positioning system (GPS) collars and other remote monitoring techniques (e.g., internal temperature monitors, movement analyses) to track survival, habitat use, causes of mortality, physiological condition, and disease transmission dynamics [17–27].

[Figure omitted. See PDF.]

Primary moose range in Minnesota (red outline) that is surveyed annually by cooperators Minnesota Department of Natural Resources (MNDNR), Fond du Lac Band of Lake Superior Chippewa (FDL), and 1854 Treaty Authority; and 4 study areas that contained collared moose: Voyageurs National Park (A), Grand Portage Indian Reservation (B), MNDNR study (2012–2016; C), and MNDNR-FDL-1854 Treaty Authority study (2005–2008; D).

Our goal was to integrate these multiple data streams into a unified model that would accurately describe past population dynamics and future projections of the northeastern MN moose population. Specifically, we used statistical population reconstruction to estimate population abundance, recruitment, and survival rates using all available data. We also examined the sensitivity of model estimates to fluctuations in adult survival and per-capita recruitment (number of calves that survived to 1 year per adult female alive during the previous year) to determine which may be more important in predicting population growth and used time series analysis to project population estimates 10 years into the future to inform management and conservation concerns. Given recent declines in moose abundance occurring broadly across North America [28, 29], our study demonstrates the utility of statistical population reconstruction for understanding moose population dynamics.

Materials and methods

Study area

Our study occurred in northeastern MN, near the southern limit of the distributional range of moose (Fig 1) [11, 28]. Our study area was a mosaic of the Superior National Forest and various Tribal, state, county, and private lands (Fig 1), as well as the federal lands of Voyageurs National Park (VNP). Moose are a subsistence food used by the Anishinaabeg (people) of the Grand Portage Band of Lake Superior Chippewa historically and presently. The Grand Portage Band is a federally recognized Indian tribe in extreme northeastern MN and proudly exercises its rights to food sovereignty through subsistence hunting and fishing. Voyageurs National Park is just west of primary moose range, which is delineated by MNDNR Section of Wildlife field and research staff (Fig 1). Moose occur outside of primary range, but at low densities. Statewide moose harvest was closed during 1922–1971, because of low moose numbers, and then reopened in the northwestern and northeastern portions of the state with limited permits issued [30]. Harvest was stopped in the northwest in 1997, but continued in the northeast. In 2007, hunters were restricted to harvesting antlered adult males only [30]. Moose harvests were then suspended in MN from 2013 until 2016, when a tribal subsistence harvest was resumed [31–33]. Moose harvests do not occur in VNP.

Our study area is part of the Northern Superior Upland within the Laurentian mixed forest province [34]. The vegetative cover is a mosaic of wetlands, stands of northern white cedar (Thuja occidentalis), black spruce (Picea mariana), tamarack (Larix laricina), and upland stands of balsam fir (Abies balsamea), jack pine (Pinus banksiana), eastern white pine (P. strobus), and red pine (P. resinosa), intermixed with quaking aspen (Populus tremuloides) and paper birch (Betula papyrifera).

Moose range in this region overlapped with gray wolves (Canis lupus) and American black bears, both of which prey upon adult and calf moose [14, 20, 25, 35, 36]. Adult and calf moose hair was present in relatively few wolf scats from VNP (0–4% occurrence) compared to scats from other areas of moose range in MN (7–22% occurrence) [37, 38]. The moose population in northeastern MN were afflicted by various parasites and disease, including infestation by winter ticks (Dermacentor albipictus) and infection by meningeal worm (Parelaphostrongylus tenuis) and giant liver fluke (Fascioloides magna) [24, 26, 39].

Aerial surveys

As part of the ongoing monitoring and management of moose in northern MN that have taken place since the 1960s, the MNDNR, in cooperation with Fond du Lac Band of Lake Superior Chippewa (FDL) and 1854 Treaty Authority, conducted an aerial survey of the northeastern moose population each winter using an updated and standardized approach since 2005 [9]; however, a survey was not conducted in 2021 due to the COVID-19 pandemic. Timing of surveys was typically during the first two weeks of January; however, insufficient snow depth postponed the 2012 survey until 26 January to 9 February [40]. The surveys were conducted using helicopters over a total area of approximately 15,500 km². This area was divided into 436 rectangular survey plots of approximately 36 km² each, 36 to 52 of which were selected each year using a stratified random sampling protocol based on moose density (low, medium, high). Moose density strata were classified collaboratively by MNDNR, FDL, and 1854 Treaty Authority staff and are reevaluated every 5 years based on expert knowledge and previous survey results. Each sighted moose was classified as either a calf, adult female, or adult male based on body size and presence of vulva patch and/or antlers; uncorrected estimates (without a sightability correction) adjusted for sampling were then used to calculate adult male:female and calf:adult female ratios at the population level [9, 41]. A sightability model was then used to estimate overall abundance. Visual obstruction was calculated as the proportion of area within a 10-m radius surrounding the first moose observed in a group that was not visible and used to adjust each estimate and corresponding 90% confidence intervals (CI; Table 1) [9, 41]. We used the estimated annual abundance of calves, adult females, and adult males derived from the aerial surveys in the IPM below. We scaled the variance of the overall point count on the proportion of calves to obtain variance estimates for calf abundance.

[Figure omitted. See PDF.]

Age-class-specific aerial survey data with corresponding annual totals and 90% confidence intervals for moose in northeastern Minnesota, USA, 2005–2020 [9]. Total abundance is corrected for sightability, abundance of calves, adult females, and adult males is derived from reported calf:adult female and adult male:adult female ratios.

Adult survival rates

In addition to aerial survey data, we used adult moose survival data collected via telemetry from 2005 to 2019 by four different studies throughout northeastern MN (Fig 1). We excluded animals with collar failures from the data in the year of collar failure (i.e., right-censoring), animals that died as a result of capture, and young-of-the-year from any further analysis. Collar failure was assumed to be independent of moose fate. The remaining animals in each study were pooled together to determine annual mortality and associated at-risk counts as a measure of adult survival rates (Table 2).

[Figure omitted. See PDF.]

Telemetry data from four different studies of annual mortality (v) and associated at-risk counts (n) for yearling and adult moose in northeastern Minnesota, USA, 2005–2019.

We used annual adult survival rates from two previous studies by MNDNR, FDL, and the 1854 Treaty Authority [11, 14, 24, 42]. The earlier MNDNR-FDL-1854 Treaty Authority study used 150 adult moose (95 F/55 M) collared during 2002–2008 [11]; however, we only used survival rates that coincided with the aerial survey (2005–2007). We used pooled adult survival estimates, because there was no difference in survival between males and females [11, 14]. The more recent MNDNR study was conducted from 2013 to 2016 and used 173 adult moose (123 F/50 M) [24]. Differences in survival between males and females were not reported, so we used the pooled adult survival estimates from this study. Details of animal capture, handling, collaring, and monitoring can be found in the source publications [11, 14, 24, 42].

We used 2 additional sources of adult moose survival data from study sites that are adjacent to the aerial survey area (Fig 1). Voyageurs National Park collared 21 moose (14 F/7 M) to study moose survival from 2010 to 2017. Grand Portage Indian Reservation collared 99 adult moose (76 F/23M) between 2010 and 2019. All capture and handling protocols were conducted in accordance with requirements of the University of MN Institutional Animal Care and Use Committee (protocols 1803-35736A and 0192A75532) and the guidelines of the American Society of Mammalogists [25, 43, 44]. We calculated Kaplan-Meier survival estimates using the “survival” package in Program R [45, 46]. Because adult moose captures typically occurred in mid-winter (Jan–Mar), we modeled annual survival using the calendar year (i.e., t₀ = 1 Jan) [47]. Collared moose that survived multiple years contributed an observation for each year they were alive, yielding 98 moose-years for Voyageurs National Park and 302 moose-years for Grand Portage. We used the “survdiff” function in the “survival” R package, which uses a log-rank test, to examine differences in overall survival between sexes [45, 46].

Population reconstruction of moose in MN

Population reconstruction typically begins by specifying a projection matrix to describe the change in the number of animals in each cohort over time. Consider a hypothetical population of moose divided into four classes (male and female, calves and adults) monitored over Y consecutive years, where N_ij is the abundance in winter of animals of class j in year i. Under this framework, all individuals born during the same year constitute a single cohort that is subsequently subjected to annual mortality from various causes. Previous reconstructions have then used an age-at-harvest matrix to represent each cohort [48–50]; however, with the exception of tribal subsistence harvest averaging about 40 moose per year [32, 33], moose are not regularly harvested in MN. As such, we did not explicitly model the impacts of harvest mortality. In lieu of these data, we used aerial survey data to represent each cohort as a separate diagonal, where the observed counts, a_ij, are a function of the initial abundance of the corresponding cohort and the annual survival rate (to be estimated as parameters). Simulation studies have demonstrated that statistical reconstruction provides an unbiased estimate of population abundance [49]. Due to the difficulty associated with identifying sex of moose calves during aerial surveys, and the assumptions of a 50:50 sex-ratio of calves at birth with no sex differences in first year survival, we pooled male and female calves into a single cohort, for a total of A = 3 classes (calves, adult females, and adult males; Table 1). We defined adults as moose >1.5 years old, as they are classified in the aerial survey.

An objective function or estimator was then used to determine which set of model parameters best describes the observed data. We used a chi-square objective function to model the difference between the observed and predicted number of animals in each cohort and the joint difference for the entire matrix aswhere is the cell-specific chi-square calculation [7, 51]. The difference for the cell represented by the total number of adult females in year 2 (i.e., N₂₂), for example, can be written as follows:where a₂₂ is the number of adult females in year 1 observed via aerial survey, N₁₁ and N₁₂ are the initial calf and adult female cohort abundance in year 1, S₁ is the annual survival rate in year 1 (which we assumed to be constant for males and females but different between years), and 0.5 represents the sex-at-birth ratio to separate calves into adult females and males after the first year of life [52].

In addition to aerial survey data, we used information from collared individuals with known fates to help estimate annual survival by comparing the observed number of mortalities each year to that expected under the model parameterization as follows:where S_i is again the annual survival rate in year i, n_i is the number of collared animals alive at the beginning of year i, and v_i is the number of collared adult moose that died in year i.

We then used a spectral projected gradient method using the “spg” function in the BB package in Program R [53] to numerically solve for the minimum chi-square estimate. This allowed us to directly estimate annual survival (i.e., S_i), initial cohort abundances in year 1 (i.e., N₁₁, N₁₂, N₁₃), and recruitment in subsequent years (i.e., N₂₁, N₃₁,…,N_Y1). All other female and male adult abundances were estimated based on the invariance property:

We calculated standard errors (SEs) for the minimum chi-square estimates using a numerical estimate of the inverse Hessian [48, 49, 54] using the “numDeriv” package in Program R [55]. Because reconstruction models consistently underestimate uncertainty [50], we inflated all standard errors by the goodness-of-fit scale parameter suggested by previous research [56]:where the statistic is based on the observed aerial survey data (a_ij) and their expected values under the reconstruction (N_ij). The degrees of freedom (df) are equal to A×Y−K, where K is the number of parameters estimated by the reconstruction. We then used these inflated standard errors to construct 90% confidence intervals for the model-derived estimates of annual population abundance and recruitment for moose in MN.

Sensitivity analysis of reconstruction estimates

Given the rapid decline in animals seen during aerial surveys between 2009 and 2013 (64.8% in five years), we investigated the sensitivity of reconstructed population estimates during these years by incrementally increasing either adult survival or recruitment by 0.1%, while holding the other constant, until the population decline was reversed (i.e., population abundance in 2013 was within 10% of that in 2009).

Population projection using reconstruction estimates

We projected our estimates of per-capita recruitment (number of calves that survived to one year per adult female alive during the previous year) and adult survival for an additional 10 years using the “forecast” package in Program R [57]. We then used the reconstructed estimates of calf, adult female, and adult male cohort abundance in 2020 as a starting point from which to predict cohort abundance from 2021 to 2030 using a stochastic version of the projection matrix approach described above [58].

Results

Survival estimates from collared moose

We did not detect a difference in overall survival in VNP by sex (χ²₁ = 0.20, P = 0.70). Adult annual survival estimates in years 2011 to 2017 for Voyageurs National Park ranged from 0.741 (95% CI 0.484–1.00) in 2015 to 1.00 in 2010 and 2013, with a mean annual survival of 0.893 (95% CI 0.833–0.958; Table 3). In 2010 and 2013, no collared moose mortalities occurred in VNP, precluding an estimate of variation in survival in those years. We did not detect a difference in overall survival in Grand Portage by sex (χ²₁ = 0.60, P = 0.40). Grand Portage adult annual survival in years 2010 to 2019 ranged from 0.591 (95% CI 0.417–0.837) in 2013 to 1.00 in 2012, with a mean annual survival of 0.833 (95% CI 0.794–0.874; Table 3). Because no collared moose mortalities occurred in Grand Portage Indian Reservation in 2012, we were precluded from estimating variation in survival.

[Figure omitted. See PDF.]

Estimates of annual survival and sex-ratios of collared adult moose in Voyageurs National Park and Grand Portage Indian Reservation, MN, USA, 2010–2021.

Population reconstruction of moose in MN

Using statistical population reconstruction with available aerial survey and telemetry data, we estimated fluctuations in adult survival, ranging from a maximum of 0.902 (SE = 0.043) in 2005 to a minimum of 0.690 (SE = 0.061) in 2019 (Fig 2). Per-capita recruitment (number of calves that survived to 1 year per adult female alive during the previous year) followed a similar cyclical pattern as adult survival, decreasing slightly from 0.386 (SE = 0.030) in 2005 to 0.303 (SE = 0.051) in 2019 (Fig 2). Winter moose abundance estimates showed a slow decline from an estimated 8,304 (90% CI = 7,797–8,788) animals in 2005 to 7,841 (90% CI = 6,702–8,933) in 2009 (Fig 3). This was followed by a sharp decline to 3,386 (90% CI = 2,681–4,243) animals in 2013, but remained steady afterwards to an estimated 3,163 (90% CI = 2,403–3,718) in 2020 (Fig 3). Annual recruitment followed a similar pattern and varied from a high of 1,683 (90% CI = 1,380–1,943) animals in 2005 to 502 (90% CI = 343–647) in 2020 (Fig 3).

[Figure omitted. See PDF.]

Estimated trends in annual survival (top) and per-capita recruitment (number of calves that survived to 1 year per adult female alive during the previous year) for moose in Minnesota (thick solid lines) between 2005 and 2019 based on statistical population reconstruction using integrated population models (IPMs), along with associated standard errors (error bars).

[Figure omitted. See PDF.]

Estimated trends in abundance (top) and calf recruitment (bottom) into the winter population of moose in Minnesota (thick solid lines) between 2005 and 2020 based on statistical population reconstruction using integrated population models (IPMs), along with associated 90% confidence intervals (shaded regions).

Sensitivity analysis of reconstruction estimates

Abundance estimates during the rapid decline from 2009 to 2013 were more sensitive to changes in adult survival than in recruitment. A 27.0% change in survival during the four years, while holding recruitment constant, resulted in a 2013 population abundance that was just 10% lower than that in 2009. To achieve a similar result while holding survival constant required an increase of 248.6% in recruitment during the four years.

Population projection using reconstruction estimates

Stochastic projections using forecasted fecundity and survival estimates resulted in a slowly decreasing population from a high of 3,244 (90% CI = 2,936–3,461) in 2021 to a low of 2,680 (90% CI = 1,298–4,550) in 2030, and a corresponding annual growth rate of 0.984 (90% CI = 0.940–1.020; Fig 4).

[Figure omitted. See PDF.]

Stochastic population projection of moose in Minnesota from 2020 to 2030 using forecasted estimates of annual survival and per-capita recruitment (number of calves that survived to 1 year per adult female alive during the previous year). Shaded regions represent 90% confidence intervals from 1,000 individual simulations.

Discussion

Statistical population reconstruction was consistent with a substantial decline in the northeastern MN moose population between 2009 and 2013, as was indicated by the original aerial surveys conducted throughout the region. Reconstruction estimates indicated that during this time, the number of moose in primary moose range in MN decreased substantially from about 7,800 animals in 2009 to about 3,400 in 2013, corresponding to a >50% decline over just four years. Since 2013, however, the population largely stabilized and displayed an oscillatory pattern with a slight overall decrease of approximately 6.6% over the next seven years to an estimated 3,163 (90% CI = 2,403–3,718) animals in 2020. Stochastic projections using forecasted demographic rates indicated that this trend is likely to continue for the next 10 years to an estimated 2,680 (90% CI = 1,298–4,550) animals in 2030, yet the 90% CI of λ included 1. This estimate closely matches simulated populations under a constant harvest of 150 adult males each year, but is less than populations under low harvest (40–80 adult males/yr; ~4,000 moose) [59].

Our results demonstrate the utility of using statistical population reconstruction to monitor moose population trends throughout northeastern MN and other parts of their North American range. When compared to estimates derived from aerial surveys, reconstruction estimates produced substantially narrower confidence intervals around similarly sized abundance estimates. For example, both aerial surveys and population reconstruction estimated similar abundances of 8,161 and 8,304 animals in 2005, respectively. However, the confidence interval around this reconstructed point estimate was approximately 20% of the confidence intervals around the aerial survey point estimate, a five-fold increase in precision. Although the increase in precision gained from reconstruction was substantially lower during many of the other years, reconstruction nonetheless provided a consistent improvement in precision when compared to estimates derived from aerial surveys (Fig 2). Moose sightability at the time of aerial surveys, owing to individual moose behavior, habitat use, and weather, can contribute added variability affecting point estimates. Using statistical population reconstruction also eliminated the biologically unrealistic fluctuations in population abundance observed in the original aerial survey estimates. For example, aerial survey estimates indicated that moose abundance rebounded from 2,760 (2,160–3,650) animals in 2013 to 4,350 (3,220–6,210) in 2014, representing an increase of 57.1% in just one year. Given moose reproductive patterns, such a steep increase over such a short period of time is biologically impossible [52, 60]. Reconstruction estimates during the same time period, on the other hand, indicate an increase of only 13.4%, from 3,386 (90% CI = 2,681–4,243) to 3,840 (90% CI = 3,146–4,650), which is a reasonable increase given reported moose reproduction estimates [52, 60]. Conversely, the minimum modeled adult survival rate of 0.69 in 2019 may warrant caution, because such a low estimate is not biologically consistent with observed population trends and calf:adult female ratios [59].

An additional benefit of using statistical population reconstruction to monitor moose throughout the northeastern region is that it can retroactively provide abundance estimates during years when aerial surveys are not conducted. The COVID-19 pandemic prevented the statewide aerial survey for moose in 2021. After aerial survey and telemetry data are collected in subsequent years, statistical population reconstruction can be used to impute the missing number of calves, adult females, and adult males in 2021. Similarly, estimates of annual survival derived from telemetry studies often include years where no animals were monitored. In the present study, to our knowledge there were no published telemetry data collected in 2008 and 2009, precluding a direct estimate of survival during those years. However, with the use of statistical population reconstruction, we were able to estimate survival during those years.

Annual survival of moose in MN appears to follow a pattern of years of high survival followed by years of low survival (Fig 3). However, there was a consistent period of low survival between 2009 and 2013, corresponding to the observed and subsequently confirmed population decline of moose during this time [24]. Combined with the results of the sensitivity analysis, which indicated that population growth is more sensitive to fluctuations in adult survival than in per capita recruitment, these results suggest that the observed decline in population abundance was most likely caused by lower adult survival from 2009 to 2013. The two lowest collared moose survival rates measured on the Grand Portage Indian Reservation also occurred during this time period. Additionally, opportunistically collected free-ranging moose that were necropsied showed health and disease issues were common during this same period [39]. Subsequent research on cause-specific mortality of adult moose in northeastern MN further highlighted the significant effect of disease and parasites, such as winter tick and meningeal worm, on adult moose survival [24, 26, 39]. Additionally, wolf populations may have been subsidized by white-tailed deer in areas of moose range in MN, leading to declines in moose numbers via apparent competition and inverse-density-dependent predation [61]. Moose population dynamics, like those of many other large herbivores, are more impacted by variation in adult survival compared to juvenile survival [62, 63]. Adult survival typically varies little [63], but in populations exhibiting low and variable adult survival, populations decline [64]. Conversely, increases in adult survival can improve population performance [65, 66].

We believe our general approach was useful for a more comprehensive assessment of moose population dynamics of northeastern MN based on the integration of several different sources of information (i.e., aerial surveys and four separate telemetry studies). Future research should build upon this foundation to explore how the incorporation of other supporting data can improve reconstruction estimates and help to estimate additional model parameters not considered here. Data on annual pregnancy rates, calf survival, and twinning rates, for example, could be used to separate the effects of reproductive success from calf mortality, thereby allowing us to better identify the driving forces behind observed trends in annual recruitment. Additional finer-scale studies, such as those ongoing at Grand Portage (S. A. Moore, unpublished data), that incorporate predator density, experimental manipulations of predator density, and effects of alternate (non-moose) prey of predators will be useful in teasing apart factors driving recruitment and mortality.

Conclusion

Statistical population reconstructions confirmed that moose abundance in northeastern MN declined rapidly from 2009 to 2013 but has remained relatively stable during 2013–2020. Our results suggest that this decline was due primarily to low adult survival during those years. Our approach increased precision of population estimates gained from the state’s annual aerial survey and can further be used to impute missing values when surveys cannot be conducted, such as occurred in 2021 due to the COVID-19 pandemic. Continued monitoring of vital rates of collared moose through the use of telemetry, such as that continuing to be undertaken by Grand Portage Band of Lake Superior Chippewa on Grand Portage Indian Reservation and in ceded territory in Superior National Forest, will aid in refining future estimates of population trends and projections and contribute to more precise knowledge of the population across time. As of publication, a moratorium on state-permitted collaring of moose is still in effect; this order does not restrict tribal activities (Executive Order 15–10, 28 Apr 2015). Without additional data streams to inform the aerial survey estimates, projections are less useful to managers of moose populations, especially when explicit mechanisms driving the trends are unknown.

Acknowledgments

This material is based upon work supported in part by the Center for Applied Mathematics (CAM) Summer Research Program at the University of St. Thomas. We thank T. Arnold for inspiring this work, J. Fieberg for providing parameters from Lenarz et al. (2010), the MNDNR Wildlife Health Program, and dozens of volunteers and technicians for all moose capture and monitoring work. Thank you to T. Jung and an anonymous reviewer for comments that improved the clarity of the manuscript.

Federally recognized Indian tribes in northeastern Minnesota, USA, including the Grand Portage Band of Lake Superior Chippewa, Bois Forte Band of Chippewa, and Fond du Lac Band of Lake Superior Chippewa, proudly exercise their rights to food sovereignty through subsistence hunting and fishing. Moose are a primary subsistence food used by the Anishinaabeg (people) historically and presently. Management for and research on maintaining this moose population as a vital subsistence species thus sets the context for this paper examining the population trends of this culturally important resource.

Citation: Severud WJ, Berg SS, Ernst CA, DelGiudice GD, Moore SA, Windels SK, et al. (2022) Statistical population reconstruction of moose (Alces alces) in northeastern Minnesota using integrated population models. PLoS ONE 17(9): e0270615. https://doi.org/10.1371/journal.pone.0270615

About the Authors:

William J. Severud

Contributed equally to this work with: William J. Severud, Sergey S. Berg

Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Software, Supervision, Visualization, Writing – original draft, Writing – review & editing

E-mail: [email protected]

Affiliation: Department of Veterinary Population Medicine, University of Minnesota, Saint Paul, Minnesota, United States of America

https://orcid.org/0000-0003-0150-5986

Sergey S. Berg

Contributed equally to this work with: William J. Severud, Sergey S. Berg

Roles: Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Department of Computer and Information Sciences, University of St. Thomas, Saint Paul, Minnesota, United States of America

https://orcid.org/0000-0002-2707-294X

Connor A. Ernst

Roles: Conceptualization, Formal analysis, Methodology, Writing – original draft

Affiliation: Department of Mathematics, University of St. Thomas, Saint Paul, Minnesota, United States of America

Glenn D. DelGiudice

Roles: Data curation, Funding acquisition, Investigation, Writing – review & editing

Affiliation: Forest Wildlife Populations and Research Group, Minnesota Department of Natural Resources, Forest Lake, Minnesota, United States of America

Seth A. Moore

Roles: Conceptualization, Data curation, Funding acquisition, Writing – review & editing

Affiliation: Department of Biology and Environment, Grand Portage Band of Lake Superior Chippewa, Grand Portage, Minnesota, United States of America

Steve K. Windels

Roles: Data curation, Funding acquisition, Writing – review & editing

Affiliation: Voyageurs National Park, International Falls, Minnesota, United States of America

Ron A. Moen

Roles: Data curation, Funding acquisition, Methodology, Writing – review & editing

Affiliation: Center for Water and the Environment, University of Minnesota, Duluth, Minnesota, United States of America

Edmund J. Isaac

Roles: Data curation, Investigation, Writing – review & editing

Affiliation: Department of Biology and Environment, Grand Portage Band of Lake Superior Chippewa, Grand Portage, Minnesota, United States of America

Tiffany M. Wolf

Roles: Data curation, Funding acquisition, Writing – review & editing

Affiliation: Department of Veterinary Population Medicine, University of Minnesota, Saint Paul, Minnesota, United States of America

https://orcid.org/0000-0003-1740-4442

References

1. Zielinski WJ, Slauson KM, Carroll CR, Kent CJ, Kudrna DG. Status of American martens in coastal forests of the Pacific states. Journal of Mammalogy. 2001;82: 478–490.

2. Heilbrun RD, Silvy NJ, Peterson MJ, Tewes ME. Estimating bobcat abundance using automatically triggered cameras. Wildlife Society Bulletin. 2006;34: 69–73.

3. Rivest L-P, Couturier S, Crépeau H. Statistical methods for estimating caribou abundance using postcalving aggregations detected by radio telemetry. Biometrics. 1998;54: 865–876.

4. Arnold TW, Clark RG, Koons DN, Schaub M. Integrated population models facilitate ecological understanding and improved management decisions: Integrated Population Models. Jour Wild Mgmt. 2018;82: 266–274.

5. Zipkin EF, Saunders SP. Synthesizing multiple data types for biological conservation using integrated population models. Biological Conservation. 2018;217: 240–250.

6. Berg SS, Erb JD, Fieberg JR, Forester JD. Utility of radio-telemetry data for improving statistical population reconstruction. Journal of Wildlife Management. 2017;81: 535–544.

7. Fieberg JR, Shertzer KW, Conn PB, Noyce KV, Garshelis DL. Integrated population modeling of Black Bears in Minnesota: Implications for Monitoring and Management. PLOS ONE. 2010;5: e12114. pmid:20711344

8. Howard AL, Clement MJ, Peck FR, Rubin ES. Estimating mountain lion abundance in Arizona using statistical population reconstruction. Journal of Wildlife Management. 2020;84: 85–95.

9. DelGiudice GD. 2020 aerial moose survey. Saint Paul, USA: Minnesota Department of Natural Resources; 2020 pp. 1–8.

10. Severud WJ, DelGiudice GD, Bump JK. Comparing survey and multiple recruitment–mortality models to assess growth rates and population projections. Ecol Evol. 2019;9: 12613–12622. pmid:31788201

11. Lenarz MS, Fieberg J, Schrage MW, Edwards AJ. Living on the edge: Viability of moose in northeastern Minnesota. Journal of Wildlife Management. 2010;74: 1013–1023.

12. ArchMiller AA, Dorazio RM, St. Clair K, Fieberg JR. Time series sightability modeling of animal populations. PLoS ONE. 2018;13: e0190706. pmid:29329309

13. DelGiudice GD. 2013 aerial moose survey final results. Saint Paul, USA: Minnesota Department of Natural Resources; 2013 pp. 1–6.

14. Lenarz MS, Nelson ME, Schrage MW, Edwards AJ. Temperature mediated moose survival in northeastern Minnesota. Journal of Wildlife Management. 2009;73: 503–510.

15. Butler E, Carstensen M, DelGiudice GD. Proposal title: Determining causes of death in declining moose population. Environmental and Natural Resources Trust Fund-LCCMR. Saint Paul: Minnesota Department of Natural Resources; 2011.

16. DelGiudice GD, Severud WJ, Butler E, Carstensen M, Moen R. Proposal title: Determining temporal and spatial survival and cause-specific mortality of moose (Alces alces) calves in northeastern Minnesota. Saint Paul: Minnesota Department of Natural Resources; 2012.

17. DelGiudice GD, Severud WJ, Obermoller TR, Wright RG, Enright TA, St-Louis V. Monitoring movement behavior enhances recognition and understanding of capture-induced abandonment of moose neonates. Journal of Mammalogy. 2015;96: 1005–1016.

18. DelGiudice GD, Severud WJ, Obermoller TR, St-Louis V. Gaining a deeper understanding of capture-induced abandonment of moose neonates. Journal of Wildlife Management. 2018;82: 287–298.

19. DelGiudice GD, Ahmadkhani M, St-Louis V, Severud WJ, Obermoller TR. Exploring the role of parental proximity in the maternal–neonate bond and parental investment in moose (Alces alces) through postcapture movement dynamics. Ecology and Evolution. 2021;11: 8699–8712.

20. Severud WJ, DelGiudice GD, Obermoller TR, Enright TA, Wright RG, Forester JD. Using GPS collars to determine parturition and cause-specific mortality of moose calves. Wildlife Society Bulletin. 2015;39: 616–625.

21. Severud WJ, DelGiudice GD, Obermoller TR. Survival, cause-specific mortality, and space use of moose calves in a declining population. In: Cornicelli L, Carstensen M, D’Angelo G, Larson MA, Lawrence JS, editors. Summaries of Wildlife Research Findings 2015. Saint Paul, USA: Minnesota Department of Natural Resources; 2017. pp. 110–123.

22. Herberg AM, St-Louis V, Carstensen M, Fieberg J, Thompson DP, Crouse JA, et al. Calibration of a rumen bolus to measure continuous internal body temperature in moose. Wildlife Society Bulletin. 2018;42: 328–337.

23. Obermoller TR. Using movement behavior of adult female moose to estimate survival and cause-specific mortality of calves in a declining population. Thesis, University of Minnesota. 2017.

24. Carstensen M, Hildebrand EC, Plattner D, Dexter M, St-Louis V, Jennelle C, et al. Determining cause-specific mortality of adult moose in northeast Minnesota, February 2013–July 2017. In: Cornicelli L, Carstensen M, Larson MA, Davros N, Davis B, editors. Summaries of Wildlife Research Findings 2016. Saint Paul, USA: Minnesota Department of Natural Resources; 2018. pp. 236–247.

25. Wolf TM, Chenaux-Ibrahim YM, Isaac EJ, Moore SA. Neonate health and calf mortality in a declining population of North American moose (Alces alces americanus). Journal of Wildlife Diseases. 2021;57: 40–50. pmid:33635990

26. Oliveira-Santos LGR, Moore SA, Severud WJ, Forester JD, Isaac EJ, Chenaux-Ibrahim Y, et al. Spatial compartmentalization: A nonlethal predator mechanism to reduce parasite transmission between prey species. Science Advances. 2021;7: eabj5944. pmid:34936450

27. Van de Vuurst P, Moore SA, Isaac EJ, Chenaux-Ibrahim Y, Wolf TM, Escobar LE. Reconstructing landscapes of ungulate parturition and predation using vegetation phenology. Wang Y-P, editor. Current Zoology. 2021; zoab058. pmid:35592351

28. Timmermann HR, Rodgers AR. The status and management of moose in North America-circa 2015. Alces. 2017;53: 1–22.

29. Jensen WF, Smith JR, Carstensen M, Penner CE, Hosek BM, Maskey JJ. Expanding GIS analyses to monitor and assess North American moose distribution and density. Alces. 2018;54: 45–54.

30. Minnesota Department of Natural Resources [MNDNR]. Minnesota moose research and management plan. Saint Paul, USA: Division of Fish and Wildlife; 2012 p. 51.

31. DelGiudice GD. 2012 Minnesota moose harvest. Saint Paul, USA: Minnesota Department of Natural Resources; 2012.

32. Edwards AJ. 1994–2017 big game and furbearer harvest report. Duluth, USA: 1854 Treaty Authority, Resource Management Division; 2018 p. 12.

33. Schrage M. 2016 wildlife harvest record for the Fond du Lac Band of Lake Superior Chippewa. Cloquet, USA: Fond du Lac Resource Management Division; 2018 p. 12.

34. MNDNR. Ecological classification system. Minnesota Department of Natural Resources, Saint Paul, Minnesota, USA. 2015 [cited 30 Nov 2016]. Available: http://www.dnr.state.mn.us/ecs/index.html

35. Patterson BR, Benson JF, Middel KR, Mills KJ, Silver A, Obbard ME. Moose calf mortality in central Ontario, Canada. Journal of Wildlife Management. 2013;77: 832–841.

36. Severud WJ, DelGiudice GD, Obermoller TR, Ryan RJ, Smith BD. An alternate method to determine moose calving and cause-specific mortality of calves in northeastern Minnesota. In: Cornicelli L, Carstensen M, Grund MD, Larson MA, Lawrence JS, editors. Summaries of Wildlife Research Findings 2014. Saint Paul, USA: Minnesota Department of Natural Resources; 2015. pp. 93–108.

37. Chenaux-Ibrahim Y. Seasonal diet composition of gray wolves (Canis lupus) in northeastern Minnesota determined by scat analysis. Thesis, University of Minnesota, Duluth. 2015.

38. Gable TD, Windels SK, Bruggink JG, Barber-Meyer SM. Weekly summer diet of gray wolves (Canis lupus) in northeastern Minnesota. American Midland Naturalist. 2018;179: 15–27.

39. Wünschmann A, Armien AG, Butler E, Schrage M, Stromberg B, Bender JB, et al. Necropsy findings in 62 opportunistically collected free-ranging moose (Alces alces) from Minnesota, USA (2003–2013). Journal of Wildlife Diseases. 2015;51: 157–165. pmid:25390764

40. Lenarz MS. 2012 aerial moose survey. Saint Paul, USA: Minnesota Department of Natural Resources; 2013.

41. Giudice J. Analysis report: 2020 MNDNR Aerial Moose Survey. Saint Paul, USA: Biometrics Unit, Section of Wildlife, Minnesota Department of Natural Resources; 2020.

42. Lenarz MS, Nelson ME, Schrage MW, Edwards AJ. Moose population dynamics in northeastern Minnesota. In: Wingate PJ, Kimmel RO, Lawrence JS, Lenarz MS, editors. Summaries of Wildlife Research Findings 2005. Saint Paul, USA: Minnesota Department of Natural Resources; 2006. pp. 34–38.

43. Sikes RS, the Animal Care and Use Committee of the American Society of Mammalogists. 2016 Guidelines of the American Society of Mammalogists for the use of wild mammals in research and education. Journal of Mammalogy. 2016;97: 663–688. pmid:29692469

44. Windels SK. 2014 Voyageurs National Park moose population survey report. Fort Collins, CO, USA: National Park Service; 2014.

45. R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2021. Available: https://www.R-project.org/

46. Therneau T. A package for survival analysis in R. 2021. Available: https://CRAN.R-project.org/package=survival

47. Fieberg J, DelGiudice GD. What time is it? Choice of time origin and scale in extended proportional hazards models. Ecology. 2009;90: 1687–1697. pmid:19569383

48. Gast CM, Skalski JR, Isabelle JL, Clawson MV. Random effects models and multistage estimation procedures for statistical population reconstruction of small game populations. PLOS ONE. 2013;8: e65244. pmid:23755199

49. Gast C, Skalski JR, Beyer DE. Evaluation of fixed- and random-effects models and multistage estimation procedures in statistical population reconstruction. Journal of Wildlife Management. 2013;77: 1258–1270.

50. Gast CM. Fixed and random effects models and multistage estimation procedures for statistical population reconstructions. Dissertation, University of Washington. 2012. Available: https://www.proquest.com/docview/1011325940/abstract/3616EC7A498B42A7PQ/1

51. Berg SS, Palmer LL. A comparison of multinomial likelihood and chi-square approaches to statistical population reconstruction. J Biol Syst. 2021;29: 543–559.

52. Boer AH. Fecundity of North American moose (Alces alces): A review. Alces Supplement. 1992;1: 1–10.

53. Varadhan R, Gilbert P. BB: An R package for solving a large system of nonlinear equations and for optimizing a high-dimensional nonlinear objective function. Journal of Statistical Software. 2010;32: 1–26.

54. Skalski JR, Millspaugh JJ, Clawson MV. Comparison of statistical population reconstruction using full and pooled adult age-class data. PLOS ONE. 2012;7: e33910. pmid:22470491

55. Gilbert P, Varadhan R. numDeriv: Accurate numerical derivatives. R package version 2016.8–1.1. 2019. Available: https://CRAN.R-project.org/package=numDeriv

56. Skalski JR, Townsend RL, Gilbert BA. Calibrating statistical population reconstruction models using catch-effort and index data. Journal of Wildlife Management. 2007;71: 1309–1316.

57. Hyndman RJ, Khandakar Y. Automatic time series forecasting: The forecast package for R. Journal of Statistical Software. 2008;27: 1–22.

58. Berg S. Using demographically stochastic modeling to study the effects of cub survival on Amur leopard population trends. 2nd Stochastic Modeling Techniques and Data Analysis International Conference. 2012. Available: https://ir.stthomas.edu/cas_stats_pub/11

59. Moen RA. Modeling harvest scenarios for the northeastern Minnesota moose population. University of Minnesota Duluth: Natural Resources Research Institute; 2021 p. 20. Report No.: NRRI/TRC-2021/04, Release 1.0.

60. Ruprecht JS, Hersey KR, Hafen K, Monteith KL, DeCesare NJ, Kauffman MJ, et al. Reproduction in moose at their southern range limit. Journal of Mammalogy. 2016;97: 1355–1365.

61. Barber-Meyer SM, Mech LD. White-tailed deer (Odocoileus virginianus) subsidize gray wolves (Canis lupus) during a moose (Alces americanus) decline: A case of apparent competition? Canadian Field-Naturalist. 2016;130: 308–314.

62. Gaillard JM, Festa-Bianchet M, Yoccoz NG, Loison A, Toigo C. Temporal variation in fitness components and population dynamics of large herbivores. Annual Review of Ecology and Systematics. 2000;31: 367–393.

63. Gaillard J-M, Festa-Bianchet M, Yoccoz NG. Population dynamics of large herbivores: variable recruitment with constant adult survival. Trends in Ecology and Evolution. 1998;13: 58–63. pmid:21238201

64. Murray DL, Cox EW, Ballard WB, Whitlaw HA, Lenarz MS, Custer TW, et al. Pathogens, nutritional deficiency, and climate influences on a declining moose population. Wildlife Monographs. 2006;166: 1–30.

65. McNay RS, Lamb CT, Giguere L, Williams SH, Martin H, Sutherland GD, et al. Demographic responses of nearly extirpated endangered mountain caribou to recovery actions in Central British Columbia. Ecological Applications. 2022;32. pmid:35319129

66. Lamb CT, Willson R, Richter C, Owens-Beek N, Napoleon J, Muir B, et al. Indigenous-led conservation: Pathways to recovery for the nearly extirpated Klinse-Za mountain caribou. Ecological Applications. 2022;32: e2581. pmid:35319140