State-space models for the dynamics of wild animal populations

Identifying historical and future global change drivers that place species recovery at risk

Ecosystem management in the face of global change requires understanding how co-occurring threats affect species and communities. Such an understanding allows for effective management strategies to be identified and implemented. An important component of this is differentiating between factors that are within (e.g. invasive predators) or outside (e.g. drought, large wildfires) of a local manager's control. In the global biodiversity hotspot of south-western Australia, small- and medium-sized mammal species are severely affected by anthropogenic threats and environmental disturbances, including invasive predators, fire, and declining rainfall. However, the relative importance of different drivers has not been quantified. We used data from a long-term monitoring program to fit Bayesian state-space models that estimated spatial and temporal changes in the relative abundance of four threatened mammal species: the woylie (Bettongia penicillata), chuditch (Dasyurus geoffroii), koomal (Trichosurus vulpecula) and quenda (Isoodon fusciventor). We then use Bayesian structural equation modelling to identify the direct and indirect drivers of population changes, and scenario analysis to forecast population responses to future environmental change. We found that habitat loss or conversion and reduced primary productivity (caused by rainfall declines) had greater effects on species' spatial and temporal population change than the range of fire and invasive predator (the red fox Vulpes vulpes) management actions observed in the study area. Scenario analysis revealed that a greater extent of severe fire and further rainfall declines predicted under climate change, operating in concert are likely to further reduce the abundance of these species, but may be mitigated partially by invasive predator control. Considering both historical and future drivers of population change is necessary to identify the factors that risk species recovery. Given that both anthropogenic pressures and environmental disturbances can undermine conservation efforts, managers must consider how the relative benefit of conservation actions will be shaped by ongoing global change.

BIRDIE: A data pipeline to inform wetland and waterbird conservation at multiple scales

Introduction

Efforts to collect ecological data have intensified over the last decade. This is especially true for freshwater habitats, which are among the most impacted by human activity and yet lagging behind in terms of data availability. Now, to support conservation programmes and management decisions, these data need to be analyzed and interpreted; a process that can be complex and time consuming. The South African Biodiversity Data Pipeline for Wetlands and Waterbirds (BIRDIE) aims to help fast and efficient information uptake, bridging the gap between raw ecological datasets and the information final users need.

Methods

BIRDIE is a full data pipeline that takes up raw data, and estimates indicators related to waterbird populations, while keeping track of their associated uncertainty. At present, we focus on the assessment of species abundance and distribution in South Africa using two citizen-science bird monitoring datasets, namely: the African Bird Atlas Project and the Coordinated Waterbird Counts. These data are analyzed with occupancy and state-space models, respectively. In addition, a suite of environmental layers help contextualize waterbird population indicators, and link these to the ecological condition of the supporting wetlands. Both data and estimated indicators are accessible to end users through an online portal and web services.

Results and discussion

We have designed a modular system that includes tasks, such as: data cleaning, statistical analysis, diagnostics, and computation of indicators. Envisioned users of BIRDIE include government officials, conservation managers, researchers and the general public, all of whom have been engaged throughout the project. Acknowledging that conservation programmes run at multiple spatial and temporal scales, we have developed a granular framework in which indicators are estimated at small scales, and then these are aggregated to compute similar indicators at broader scales. Thus, the online portal is designed to provide spatial and temporal visualization of the indicators using maps, time series and pre-compiled reports for species, sites and conservation programmes. In the future, we aim to expand the geographical coverage of the pipeline to other African countries, and develop more indicators specific to the ecological structure and function of wetlands.

Integrated Population Models: Achieving Their Potential

Precise and accurate estimates of abundance and demographic rates are primary quantities of interest within wildlife conservation and management. Such quantities provide insight into population trends over time and the associated underlying ecological drivers of the systems. This information is fundamental in managing ecosystems, assessing species conservation status and developing and implementing effective conservation policy. Observational monitoring data are typically collected on wildlife populations using an array of different survey protocols, dependent on the primary questions of interest. For each of these survey designs, a range of advanced statistical techniques have been developed which are typically well understood. However, often multiple types of data may exist for the same population under study. Analyzing each data set separately implicitly discards the common information contained in the other data sets. An alternative approach that aims to optimize the shared information contained within multiple data sets is to use a “model-based data integration” approach, or more commonly referred to as an “integrated model.” This integrated modeling approach simultaneously analyzes all the available data within a single, and robust, statistical framework. This paper provides a statistical overview of ecological integrated models, with a focus on integrated population models (IPMs) which include abundance and demographic rates as quantities of interest. Four main challenges within this area are discussed, namely model specification, computational aspects, model assessment and forecasting. This should encourage researchers to explore further and develop new practical tools to ensure that full utility can be made of IPMs for future studies.

Density regulation amplifies environmentally induced population fluctuations

Background

Density-dependent regulation is ubiquitous in population dynamics, and its potential interaction with environmental stochasticity complicates the characterization of the random component of population dynamics. Yet, this issue has not received attention commensurate with its relevance for descriptive and predictive modeling of population dynamics. Here we use a Bayesian modeling approach to investigate the contribution of density regulation to population variability in stochastic environments.

Methods

We analytically derive a formula linking the stationary variance of population abundance/density under Gompertz regulation in a stochastic environment with constant variance to the environmental variance and the strength of density feedback, to investigate whether and how density regulation affects the stationary variance. We examine through simulations whether the relationship between stationary variance and density regulation inferred analytically under the Gompertz model carries over to the Ricker model, widely used in population dynamics modeling.

Results

The analytical decomposition of the stationary variance under stochastic Gompertz dynamics implies higher variability for strongly regulated populations. Simulation results demonstrate that the pattern of increasing population variability with increasing density feedback found under the Gompertz model holds for the Ricker model as well, and is expected to be a general phenomenon with stochastic population models. We also analytically established and empirically validated that the square of the autoregressive parameter of the Gompertz model in AR(1) form represents the proportion of stationary variance due to density dependence.

Discussion

Our results suggest that neither environmental stochasticity nor density regulation can alone explain the patterns of population variability in stochastic environments, as these two components of temporal variation interact, with a tendency for density regulation to amplify the magnitude of environmentally induced population fluctuations. This finding has far-reaching implications for population viability. It implies that intense intra-specific resource competition increases the risk of environment-driven population collapse at high density, making opportune harvesting a sensible practice for improving the resistance of managed populations such as fish stocks to environmental perturbations. The separation of density-dependent and density-independent processes will help improve population dynamics modeling, while providing a basis for evaluating the relative importance of these two categories of processes that remains a topic of long-standing controversy among ecologists.

Quantifying the checks and balances of collaborative governance systems for adaptive carnivore management

Recovering or threatened carnivore populations are often harvested to minimise their impact on human activities, such as livestock farming or game hunting. Increasingly, harvest quota decisions involve a set of scientific, administrative and political institutions operating at national and sub‐national levels whose interactions and collective decision‐making aim to increase the legitimacy of management and ensure population targets are met. In practice, however, assessments of how quota decisions change between these different actors and what consequences these changes have on population trends are rare.

We combine a state‐space population modelling approach with an analysis of quota decisions taken at both regional and national levels between 2007 and 2018 to build a set of decision‐making models that together predict annual harvest quota values for Eurasian lynx (Lynx lynx) in Norway.

We reveal a tendency for administrative decision‐makers to compensate for consistent quota increases by political actors, particularly when the lynx population size estimate is above the regional target. Using population forecasts based on the ensemble of decision‐making models, we show that such buffering of political biases ensures lynx population size remains close to regional and national targets in the long term.

Our results go beyond the usual qualitative assessment of collaborative governance systems for carnivore management, revealing a system of checks and balances that, in the case of lynx in Norway, ensures both multi‐stakeholder participation and sustainable harvest quotas. Nevertheless, we highlight important inter‐regional differences in decision‐making and population forecasts, the socio‐ecological drivers of which need to be better understood to prevent future population declines.

Synthesis and applications. Our work analyses the sequence of decisions leading to yearly quotas for lynx harvest in Norway, highlighting the collaborative and structural processes that together shape harvest sustainability. In doing so, we provide a predictive framework to evaluate participatory decision‐making processes in wildlife management, paving the way for scientists and decision‐makers to collaborate more widely in identifying where decision biases might lie and how institutional arrangements can be optimised to minimise them. We emphasise, however, that this is only possible if wildlife management decisions are documented and transparent.

Aerial survey estimates of polar bears and their tracks in the Chukchi Sea

Polar bears are of international conservation concern due to climate change but are difficult to study because of low densities and an expansive, circumpolar distribution. In a collaborative U.S.-Russian effort in spring of 2016, we used aerial surveys to detect and estimate the abundance of polar bears on sea ice in the Chukchi Sea. Our surveys used a combination of thermal imagery, digital photography, and human observations. Using spatio-temporal statistical models that related bear and track densities to physiographic and biological covariates (e.g., sea ice extent, resource selection functions derived from satellite tags), we predicted abundance and spatial distribution throughout our study area. Estimates of 2016 abundance ([Formula: see text]) ranged from 3,435 (95% CI: 2,300-5,131) to 5,444 (95% CI: 3,636-8,152) depending on the proportion of bears assumed to be missed on the transect line during Russian surveys (g(0)). Our point estimates are larger than, but of similar magnitude to, a recent estimate for the period 2008-2016 ([Formula: see text]; 95% CI 1,522-5,944) derived from an integrated population model applied to a slightly smaller area. Although a number of factors (e.g., equipment issues, differing platforms, low sample sizes, size of the study area relative to sampling effort) required us to make a number of assumptions to generate estimates, it establishes a useful lower bound for abundance, and suggests high spring polar bear densities on sea ice in Russian waters south of Wrangell Island. With future improvements, we suggest that springtime aerial surveys may represent a plausible avenue for studying abundance and distribution of polar bears and their prey over large, remote areas.

Separating the effects of climate, bycatch, predation and harvesting on tītī (Ardenna grisea) population dynamics in New Zealand: A model-based assessment

A suite of factors may have contributed to declines in the tītī (sooty shearwater; Ardenna grisea) population in the New Zealand region since at least the 1960s. Recent estimation of the magnitude of most sources of non-natural mortality has presented the opportunity to quantitatively assess the relative importance of these factors. We fit a range of population dynamics models to a time-series of relative abundance data from 1976 until 2005, with the various sources of mortality being modelled at the appropriate part of the life-cycle. We present estimates of effects obtained from the best-fitting model and using model averaging. The best-fitting models explained much of the variation in the abundance index when survival and fecundity were linked to the Southern Oscillation Index, with strong decreases in adult survival, juvenile survival and fecundity being related to El Niño-Southern Oscillation (ENSO) events. Predation by introduced animals, harvesting by humans, and bycatch in fisheries also appear to have contributed to the population decline. It is envisioned that the best-fitting models will form the basis for quantitative assessments of competing management strategies. Our analysis suggests that sustainability of the New Zealand tītī population will be most influenced by climate, in particular by how climate change will affect the frequency and intensity of ENSO events in the future. Removal of the effects of both depredation by introduced predators and harvesting by humans is likely to have fewer benefits for the population than alleviating climate effects.

Uncovering ecological state dynamics with hidden Markov models

Ecological systems can often be characterised by changes among a finite set of underlying states pertaining to individuals, populations, communities or entire ecosystems through time. Owing to the inherent difficulty of empirical field studies, ecological state dynamics operating at any level of this hierarchy can often be unobservable or 'hidden'. Ecologists must therefore often contend with incomplete or indirect observations that are somehow related to these underlying processes. By formally disentangling state and observation processes based on simple yet powerful mathematical properties that can be used to describe many ecological phenomena, hidden Markov models (HMMs) can facilitate inferences about complex system state dynamics that might otherwise be intractable. However, HMMs have only recently begun to gain traction within the broader ecological community. We provide a gentle introduction to HMMs, establish some common terminology, review the immense scope of HMMs for applied ecological research and provide a tutorial on implementation and interpretation. By illustrating how practitioners can use a simple conceptual template to customise HMMs for their specific systems of interest, revealing methodological links between existing applications, and highlighting some practical considerations and limitations of these approaches, our goal is to help establish HMMs as a fundamental inferential tool for ecologists.

State-space models reveal a continuing elephant poaching problem in most of Africa

The most comprehensive data on poaching of African elephants comes from the Monitoring the Illegal Killing of Elephants (MIKE) program, which reports numbers of illegally killed carcasses encountered by rangers. Recent studies utilizing MIKE data have reported that poaching of African elephants peaked in 2011 and has been decreasing through 2018. Closer examination of these studies, however, raises questions about the conclusion that poaching is decreasing throughout the continent. To provide more accurate information on trends in elephant poaching, we analyzed MIKE data using state-space models. State-space models account for missing data and the error inherent when sampling carcasses. Using the state-space model, for 2011–2018, we found no significant temporal trends in rates of illegal killing for Southern, Central and Western Africa. Only in Eastern Africa have poaching rates decreased substantially since 2011. For Africa as a whole, poaching did decline for 2011–2018, but the decline was entirely due to Eastern African sites. Our results suggest that poaching for ivory has not diminished across most of Africa since 2011. Continued vigilance and anti-poaching efforts will be necessary to combat poaching and to conserve African elephants.

An improved understanding of ungulate population dynamics using count data: Insights from western Montana

Understanding the dynamics of ungulate populations is critical given their ecological and economic importance. In particular, the ability to evaluate the evidence for potential drivers of variation in population trajectories is important for informed management. However, the use of age ratio data (e.g., juveniles:adult females) as an index of variation in population dynamics is hindered by a lack of statistical power and difficult interpretation. Here, we show that the use of a population model based on count, classification and harvest data can dramatically improve the understanding of ungulate population dynamics by: 1) providing estimates of vital rates (e.g., per capita recruitment and population growth) that are easier to interpret and more useful to managers than age ratios and 2) increasing the power to assess potential sources of variation in key vital rates. We used a time series of elk (Cervus canadensis) spring count and classification data (2004 to 2016) and fall harvest data from hunting districts in western Montana to construct a population model to estimate vital rates and assess evidence for an association between a series of environmental covariates and indices of predator abundance on per capita recruitment rates of elk calves. Our results suggest that per capita recruitment rates were negatively associated with cold and wet springs, and severe winters, and positively associated with summer precipitation. In contrast, an analysis of the raw age ratio data failed to detect these relationships. Our approach based on a population model provided estimates of the region-wide mean per capita recruitment rate (mean = 0.25, 90% CI = 0.21, 0.29), temporal variation in hunting-district-specific recruitment rates (minimum = 0.09; 90% CI = [0.07, 0.11], maximum = 0.43; 90% CI = [0.38, 0.48]), and annual population growth rates (minimum = 0.83; 90% CI = [0.78, 0.87], maximum = 1.20; 90% CI = [1.11, 1.29]). We recommend using routinely collected population count and classification data and a population modeling approach rather than interpreting estimated age ratios as a substantial improvement in understanding population dynamics.

Data Integration for Large-Scale Models of Species Distributions

With the expansion in the quantity and types of biodiversity data being collected, there is a need to find ways to combine these different sources to provide cohesive summaries of species' potential and realized distributions in space and time. Recently, model-based data integration has emerged as a means to achieve this by combining datasets in ways that retain the strengths of each. We describe a flexible approach to data integration using point process models, which provide a convenient way to translate across ecological currencies. We highlight recent examples of large-scale ecological models based on data integration and outline the conceptual and technical challenges and opportunities that arise.

Combining capture-recapture data and known ages allows estimation of age-dependent survival rates

In many animal populations, demographic parameters such as survival and recruitment vary markedly with age, as do parameters related to sampling, such as capture probability. Failing to account for such variation can result in biased estimates of population-level rates. However, estimating age-dependent survival rates can be challenging because ages of individuals are rarely known unless tagging is done at birth. For many species, it is possible to infer age based on size. In capture-recapture studies of such species, it is possible to use a growth model to infer the age at first capture of individuals. We show how to build estimates of age-dependent survival into a capture-mark-recapture model based on data obtained in a capture-recapture study. We first show how estimates of age based on length increments closely match those based on definitive aging methods. In simulated analyses, we show that both individual ages and age-dependent survival rates estimated from simulated data closely match true values. With our approach, we are able to estimate the age-specific apparent survival rates of Murray and trout cod in the Murray River, Australia. Our model structure provides a flexible framework within which to investigate various aspects of how survival varies with age and will have extensions within a wide range of ecological studies of animals where age can be estimated based on size.

Precision gain versus effort with joint models using detection/non-detection and banding data

Capture-recapture techniques provide valuable information, but are often more cost-prohibitive at large spatial and temporal scales than less-intensive sampling techniques. Model development combining multiple data sources to leverage data source strengths and for improved parameter precision has increased, but with limited discussion on precision gain versus effort. We present a general framework for evaluating trade-offs between precision gained and costs associated with acquiring multiple data sources, useful for designing future or new phases of current studies.We illustrated how Bayesian hierarchical joint models using detection/non-detection and banding data can improve abundance, survival, and recruitment inference, and quantified data source costs in a northern Arizona, USA, western bluebird (Sialia mexicana) population. We used an 8-year detection/non-detection (distributed across the landscape) and banding (subset of locations within landscape) data set to estimate parameters. We constructed separate models using detection/non-detection and banding data, and a joint model using both data types to evaluate parameter precision gain relative to effort.Joint model parameter estimates were more precise than single data model estimates, but parameter precision varied (apparent survival > abundance > recruitment). Banding provided greater apparent survival precision than detection/non-detection data. Therefore, little precision was gained when detection/non-detection data were added to banding data. Additional costs were minimal; however, additional spatial coverage and ability to estimate abundance and recruitment improved inference. Conversely, more precision was gained when adding banding to detection/non-detection data at higher cost. Spatial coverage was identical, yet survival and abundance estimates were more precise. Justification of increased costs associated with additional data types depends on project objectives.We illustrate a general framework for evaluating precision gain relative to effort, applicable to joint data models with any data type combination. This framework evaluates costs and benefits from and effort levels between multiple data types, thus improving population monitoring designs.

A Bayesian state-space model using age-at-harvest data for estimating the population of black bears (Ursus americanus) in Wisconsin

Population estimation is essential for the conservation and management of fish and wildlife, but accurate estimates are often difficult or expensive to obtain for cryptic species across large geographical scales. Accurate statistical models with manageable financial costs and field efforts are needed for hunted populations and using age-at-harvest data may be the most practical foundation for these models. Several rigorous statistical approaches that use age-at-harvest and other data to accurately estimate populations have recently been developed, but these are often dependent on (a) accurate prior knowledge about demographic parameters of the population, (b) auxiliary data, and (c) initial population size. We developed a two-stage state-space Bayesian model for a black bear (Ursus americanus) population with age-at-harvest data, but little demographic data and no auxiliary data available, to create a statewide population estimate and test the sensitivity of the model to bias in the prior distributions of parameters and initial population size. The posterior abundance estimate from our model was similar to an independent capture-recapture estimate from tetracycline sampling and the population trend was similar to the catch-per-unit-effort for the state. Our model was also robust to bias in the prior distributions for all parameters, including initial population size, except for reporting rate. Our state-space model created a precise estimate of the black bear population in Wisconsin based on age-at-harvest data and potentially improves on previous models by using little demographic data, no auxiliary data, and not being sensitive to initial population size.

Applying bayesian population assessment models to artisanal, multispecies fisheries in the Northern Mokran Sea, Iran

Small-scale fisheries substantially contribute to the reduction of poverty, local economies and food safety in many countries. However, limited and low-quality catches and effort data for small-scale fisheries complicate the stock assessment and management. Bayesian modelling has been advocated when assessing fisheries with limited data. Specifically, Bayesian models can incorporate information of the multiple sources, improve precision in the stock assessments and provide specific levels of uncertainty for estimating the relevant parameters. In this study, therefore, the state-space Bayesian generalised surplus production models will be used in order to estimate the stock status of fourteen Demersal fish species targeted by small-scale fisheries in Sistan and Baluchestan, Iran. The model was estimated using Markov chain Monte Carlo (MCMC) and Gibbs Sampling. Model parameter estimates were evaluated by the formal convergence and stationarity diagnostic tests, indicating convergence and accuracy. They were also aligned with existing parameter estimates for fourteen species of the other locations. This suggests model reliability and demonstrates the utility of Bayesian models. According to estimated fisheries’ management reference points, all assessed fish stocks appear to be overfished. Overfishing considered, the current fisheries management strategies for the small-scale fisheries may need some adjustments to warrant the long-term viability of the fisheries.

A state-space modelling approach to wildlife monitoring with application to flying-fox abundance

Monitoring flying-foxes is challenging as their extreme mobility produces highly dynamic population processes, considerable logistic difficulty, and variability in estimated population size. We report on methods for inferring population trend for the population of the spectacled flying-fox (Pteropus conspicillatus) in Australia. Monthly monitoring is conducted at all known roost sites across the species’ range in the Wet Tropics Region. The proportion of animals in camps varies seasonally and stochastic environmental events appear to be influential. We develop a state-space model that incorporates these processes and enables inference on total population trends and uses early warning analysis to identify the causes of population dynamics. The model suggests that population growth rate is stable in the absence of cyclones, however, cyclones appear to impact on both survival and reproduction. The population recovered after two cyclones but declined after a third. The modelling estimates a population decline over 15 years of c. 75% (mean r = − 0.12yr⁻¹ and belief of negative trend is c. 83%) suggesting that conservation action is warranted. Our work shows that a state-space modelling approach is a significant improvement on inference from raw counts from surveys and demonstrates that this approach is a workable alternative to other methods.

Range position and climate sensitivity: The structure of among-population demographic responses to climatic variation

Species' distributions will respond to climate change based on the relationship between local demographic processes and climate and how this relationship varies based on range position. A rarely tested demographic prediction is that populations at the extremes of a species' climate envelope (e.g., populations in areas with the highest mean annual temperature) will be most sensitive to local shifts in climate (i.e., warming). We tested this prediction using a dynamic species distribution model linking demographic rates to variation in temperature and precipitation for wood frogs (Lithobates sylvaticus) in North America. Using long-term monitoring data from 746 populations in 27 study areas, we determined how climatic variation affected population growth rates and how these relationships varied with respect to long-term climate. Some models supported the predicted pattern, with negative effects of extreme summer temperatures in hotter areas and positive effects on recruitment for summer water availability in drier areas. We also found evidence of interacting temperature and precipitation influencing population size, such as extreme heat having less of a negative effect in wetter areas. Other results were contrary to predictions, such as positive effects of summer water availability in wetter parts of the range and positive responses to winter warming especially in milder areas. In general, we found wood frogs were more sensitive to changes in temperature or temperature interacting with precipitation than to changes in precipitation alone. Our results suggest that sensitivity to changes in climate cannot be predicted simply by knowing locations within the species' climate envelope. Many climate processes did not affect population growth rates in the predicted direction based on range position. Processes such as species-interactions, local adaptation, and interactions with the physical landscape likely affect the responses we observed. Our work highlights the need to measure demographic responses to changing climate.

Estimating Density and Temperature Dependence of Juvenile Vital Rates Using a Hidden Markov Model

Organisms in the wild have cryptic life stages that are sensitive to changing environmental conditions and can be difficult to survey. In this study, I used mark-recapture methods to repeatedly survey Anaea aidea (Nymphalidae) caterpillars in nature, then modeled caterpillar demography as a hidden Markov process to assess if temporal variability in temperature and density influence the survival and growth of A. aidea over time. Individual encounter histories result from the joint likelihood of being alive and observed in a particular stage, and I have included hidden states by separating demography and observations into parallel and independent processes. I constructed a demographic matrix containing the probabilities of all possible fates for each stage, including hidden states, e.g., eggs and pupae. I observed both dead and live caterpillars with high probability. Peak caterpillar abundance attracted multiple predators, and survival of fifth instars declined as per capita predation rate increased through spring. A time lag between predator and prey abundance was likely the cause of improved fifth instar survival estimated at high density. Growth rates showed an increase with temperature, but the preferred model did not include temperature. This work illustrates how state-space models can include unobservable stages and hidden state processes to evaluate how environmental factors influence vital rates of cryptic life stages in the wild.

Population assessment using multivariate time-series analysis: A case study of rockfishes in Puget Sound

Estimating a population's growth rate and year-to-year variance is a key component of population viability analysis (PVA). However, standard PVA methods require time series of counts obtained using consistent survey methods over many years. In addition, it can be difficult to separate observation and process variance, which is critical for PVA. Time-series analysis performed with multivariate autoregressive state-space (MARSS) models is a flexible statistical framework that allows one to address many of these limitations. MARSS models allow one to combine surveys with different gears and across different sites for estimation of PVA parameters, and to implement replication, which reduces the variance-separation problem and maximizes informational input for mean trend estimation. Even data that are fragmented with unknown error levels can be accommodated. We present a practical case study that illustrates MARSS analysis steps: data choice, model set-up, model selection, and parameter estimation. Our case study is an analysis of the long-term trends of rockfish in Puget Sound, Washington, based on citizen science scuba surveys, a fishery-independent trawl survey, and recreational fishery surveys affected by bag-limit reductions. The best-supported models indicated that the recreational and trawl surveys tracked different, temporally independent assemblages that declined at similar rates (an average of -3.8% to -3.9% per year). The scuba survey tracked a separate increasing and temporally independent assemblage (an average of 4.1% per year). Three rockfishes (bocaccio, canary, and yelloweye) are listed in Puget Sound under the US Endangered Species Act (ESA). These species are associated with deep water, which the recreational and trawl surveys sample better than the scuba survey. All three ESA-listed rockfishes declined as a proportion of recreational catch between the 1970s and 2010s, suggesting that they experienced similar or more severe reductions in abundance than the 3.8-3.9% per year declines that were estimated for rockfish populations sampled by the recreational and trawl surveys.

Immigration ensures population survival in the Siberian flying squirrel

Linking dispersal to population growth remains a challenging task and is a major knowledge gap, for example, for conservation management. We studied relative roles of different demographic rates behind population growth in Siberian flying squirrels in two nest-box breeding populations in western Finland. Adults and offspring were captured and individually identifiable. We constructed an integrated population model, which estimated all relevant annual demographic rates (birth, local [apparent] survival, and immigration) as well as population growth rates. One population (studied 2002-2014) fluctuated around a steady-state equilibrium, whereas the other (studied 1995-2014) showed a numerical decline. Immigration was the demographic rate which showed clear correlations to annual population growth rates in both populations. Population growth rate was density dependent in both populations. None of the demographic rates nor the population growth rate correlated across the two study populations, despite their proximity suggesting that factors regulating the dynamics are determined locally. We conclude that flying squirrels may persist in a network of uncoupled subpopulations, where movement between subpopulations is of critical importance. Our study supports the view that dispersal has the key role in population survival of a small forest rodent.

Strong "bottom-up" influences on small mammal populations: State-space model analyses from long-term studies

“Bottom‐up” influences, that is, masting, plus population density and climate, commonly influence woodland rodent demography. However, “top‐down” influences (predation) also intervene. Here, we assess the impacts of masting, climate, and density on rodent populations placed in the context of what is known about “top‐down” influences. To explain between‐year variations in bank vole Myodes glareolus and wood mouse Apodemus sylvaticus population demography, we applied a state‐space model to 33 years of catch‐mark‐release live‐trapping, winter temperature, and precise mast‐collection data. Experimental mast additions aided interpretation. Rodent numbers in European ash Fraxinus excelsior woodland were estimated (May/June, November/December). December–March mean minimum daily temperature represented winter severity. Total marked adult mice/voles (and juveniles in May/June) provided density indices validated against a model‐generated population estimate; this allowed estimation of the structure of a time‐series model and the demographic impacts of the climatic/biological variables. During two winters of insignificant fruit‐fall, 6.79 g/m² sterilized ash seed (as fruit) was distributed over an equivalent woodland similarly live‐trapped. September–March fruit‐fall strongly increased bank vole spring reproductive rate and winter and summer population growth rates; colder winters weakly reduced winter population growth. September–March fruit‐fall and warmer winters marginally increased wood mouse spring reproductive rate and September–December fruit‐fall weakly elevated summer population growth. Density dependence significantly reduced both species' population growth. Fruit‐fall impacts on demography still appeared after a year. Experimental ash fruit addition confirmed its positive influence on bank vole winter population growth with probable moderation by colder temperatures. The models show the strong impact of masting as a “bottom‐up” influence on rodent demography, emphasizing independent masting and weather influences; delayed effects of masting; and the importance of density dependence and its interaction with masting. We conclude that these rodents show strong “bottom‐up” and density‐dependent influences on demography moderated by winter temperature. “Top‐down” influences appear weak and need further investigation.