Increasing drought favors nonnative fishes in a dryland river: evidence from a multispecies demographic model

Matrix community models for ecology and evolution

Ecological communities are shaped by biotic interactions as well as environmental forces, and both must be incorporated to obtain models capable of forecasting realistic community dynamics. Many community models first specify pairwise biotic interactions and then secondarily examine how extrinsic factors such as abiotic conditions affect species abundances. A disadvantage of this approach is that the species interactions themselves are often environment and context specific, making parameterization difficult. We propose an alternative approach, matrix community models (MCMs), which are sets of matrix population models linked by an assumption of aggregate density dependence. MCMs incorporate detailed species autecology but are neutral with respect to pairwise species interactions, instead allowing interactions to be revealed within the model structure. These model-revealed species interactions, including competitive exclusion, facilitation, and interference competition, shape the distribution and abundance of species within communities and generate empirically testable predictions about species interactions. We develop a framework for building MCMs using vital rates in a stochastic, multispecies framework. Single-species matrix population models are connected via an assumption of aggregate density dependence, pairwise species interactions are estimated with sensitivity analysis, and community trajectories are analyzed under different environmental regimes using standard statistical tools and network analysis. MCMs have the advantage that pairwise species interactions need not be specified a priori, and that mechanistic demographic-environment linkages permit forecasting of community dynamics under novel, non-stationary environmental regimes. A challenge is that species' autecological vital rates, such as fecundity, growth and survivorship, must be measured under a diverse range of environmental conditions to parameterize the models. We illustrate the approach with examples and discuss prospects for future theoretical and empirical developments.

Simple statistical models can be sufficient for testing hypotheses with population time-series data

Time-series data offer wide-ranging opportunities to test hypotheses about the physical and biological factors that influence species abundances. Although sophisticated models have been developed and applied to analyze abundance time series, they require information about species detectability that is often unavailable. We propose that in many cases, simpler models are adequate for testing hypotheses. We consider three relatively simple regression models for time series, using simulated and empirical (fish and mammal) datasets. Model A is a conventional generalized linear model of abundance, model B adds a temporal autoregressive term, and model C uses an estimate of population growth rate as a response variable, with the option of including a term for density dependence. All models can be fit using Bayesian and non-Bayesian methods. Simulation results demonstrated that model C tended to have greater support for long-lived, lower-fecundity organisms (K life-history strategists), while model A, the simplest, tended to be supported for shorter-lived, high-fecundity organisms (r life-history strategists). Analysis of real-world fish and mammal datasets found that models A, B, and C each enjoyed support for at least some species, but sometimes yielded different insights. In particular, model C indicated effects of predictor variables that were not evident in analyses with models A and B. Bayesian and frequentist models yielded similar parameter estimates and performance. We conclude that relatively simple models are useful for testing hypotheses about the factors that influence abundance in time-series data, and can be appropriate choices for datasets that lack the information needed to fit more complicated models. When feasible, we advise fitting datasets with multiple models because they can provide complementary information.

Life and death in a dynamic environment: Invasive trout, floods, and intraspecific drivers of translocated populations

Understanding the relative strengths of intrinsic and extrinsic factors regulating populations is a long-standing focus of ecology and critical to advancing conservation programs for imperiled species. Conservation could benefit from an increased understanding of factors influencing vital rates (somatic growth, recruitment, survival) in small, translocated populations, which is lacking owing to difficulties in long-term monitoring of rare species. Translocations, here defined as the transfer of wild-captured individuals from source populations to new habitats, are widely used for species conservation, but outcomes are often minimally monitored, and translocations that are monitored often fail. To improve our understanding of how translocated populations respond to environmental variation, we developed and tested hypotheses related to intrinsic (density dependent) and extrinsic (introduced rainbow trout Oncorhynchus mykiss, stream flow and temperature regime) causes of vital rate variation in endangered humpback chub (Gila cypha) populations translocated to Colorado River tributaries in the Grand Canyon (GC), USA. Using biannual recapture data from translocated populations over 10 years, we tested hypotheses related to seasonal somatic growth, and recruitment and population growth rates with linear mixed-effects models and temporal symmetry mark-recapture models. We combined data from recaptures and resights of dispersed fish (both physical captures and continuously recorded antenna detections) from throughout GC to test survival hypotheses, while accounting for site fidelity, using joint live-recapture/live-resight models. While recruitment only occurred in one site, which also drove population growth (relative to survival), evidence supported hypotheses related to density dependence in growth, survival, and recruitment, and somatic growth and recruitment were further limited by introduced trout. Mixed-effects models explained between 67% and 86% of the variation in somatic growth, which showed increased growth rates with greater flood-pulse frequency during monsoon season. Monthly survival was 0.56-0.99 and 0.80-0.99 in the two populations, with lower survival during periods of higher intraspecific abundance and low flood frequency. Our results suggest translocations can contribute toward the recovery of large-river fishes, but continued suppression of invasive fishes to enhance recruitment may be required to ensure population resilience. Furthermore, we demonstrate the importance of flooding to population demographics in food-depauperate, dynamic, invaded systems.

Life and death in a dynamic environment: Invasive trout, floods, and intraspecific drivers of translocated populations

Understanding the relative strengths of intrinsic and extrinsic factors regulating populations is a long-standing focus of ecology and critical to advancing conservation programs for imperiled species. Conservation could benefit from an increased understanding of factors influencing vital rates (somatic growth, recruitment, survival) in small, translocated populations, which is lacking owing to difficulties in long-term monitoring of rare species. Translocations, here defined as the transfer of wild-captured individuals from source populations to new habitats, are widely used for species conservation, but outcomes are often minimally monitored, and translocations that are monitored often fail. To improve our understanding of how translocated populations respond to environmental variation, we developed and tested hypotheses related to intrinsic (density dependent) and extrinsic (introduced rainbow trout Oncorhynchus mykiss, stream flow and temperature regime) causes of vital rate variation in endangered humpback chub (Gila cypha) populations translocated to Colorado River tributaries in the Grand Canyon (GC), USA. Using biannual recapture data from translocated populations over 10 years, we tested hypotheses related to seasonal somatic growth, and recruitment and population growth rates with linear mixed-effects models and temporal symmetry mark-recapture models. We combined data from recaptures and resights of dispersed fish (both physical captures and continuously recorded antenna detections) from throughout GC to test survival hypotheses, while accounting for site fidelity, using joint live-recapture/live-resight models. While recruitment only occurred in one site, which also drove population growth (relative to survival), evidence supported hypotheses related to density dependence in growth, survival, and recruitment, and somatic growth and recruitment were further limited by introduced trout. Mixed-effects models explained between 67% and 86% of the variation in somatic growth, which showed increased growth rates with greater flood-pulse frequency during monsoon season. Monthly survival was 0.56-0.99 and 0.80-0.99 in the two populations, with lower survival during periods of higher intraspecific abundance and low flood frequency. Our results suggest translocations can contribute toward the recovery of large-river fishes, but continued suppression of invasive fishes to enhance recruitment may be required to ensure population resilience. Furthermore, we demonstrate the importance of flooding to population demographics in food-depauperate, dynamic, invaded systems.

Modeling effects of disturbance across life history strategies of stream fishes

A central goal of population ecology is to establish linkages between life history strategy, disturbance, and population dynamics. Globally, disturbance events such as drought and invasive species have dramatically impacted stream fish populations and contributed to sharp declines in freshwater biodiversity. Here, we used RAMAS Metapop to construct stage-based demographic metapopulation models for stream fishes with periodic, opportunistic, and equilibrium life history strategies and assessed their responses to the effects of invasion (reduced carrying capacity), extended drought (reduced survival and fecundity), and the combined effects of both disturbances. Our models indicated that populations respond differentially to disturbance based on life history strategy. Equilibrium strategists were best able to deal with simulated invasion. Periodic strategists did well under lower levels of drought, whereas opportunistic strategists outperformed other life histories under extreme seasonal drought. When we modeled additive effects scenarios, these disturbances interacted synergistically, dramatically increasing terminal extinction risk for all three life history strategies. Modeling exercises that examine broad life history categories can help to answer fundamental ecological questions about the relationship between disturbance resilience and life history, as well as help managers to develop generalized conservation strategies when species-specific data are lacking. Our results indicate that life history strategy is a fundamental determinant of population trajectories, and that disturbances can interact synergistically to dramatically impact extinction outcomes.

Resistance and resilience of pelagic and littoral fishes to drought in the San Francisco Estuary

Many estuarine ecosystems and the fish communities that inhabit them have undergone substantial changes in the past several decades, largely due to multiple interacting stressors that are often of anthropogenic origin. Few are more impactful than droughts, which are predicted to increase in both frequency and severity with climate change. In this study, we examined over five decades of fish monitoring data from the San Francisco Estuary, California, USA, to evaluate the resistance and resilience of fish communities to disturbance from prolonged drought events. High resistance was defined by the lack of decline in species occurrence from a wet to a subsequent drought period, while high resilience was defined by the increase in species occurrence from a drought to a subsequent wet period. We found some unifying themes connecting the multiple drought events over the 50-yr period. Pelagic fishes consistently declined during droughts (low resistance), but exhibit a considerable amount of resiliency and often rebound in the subsequent wet years. However, full recovery does not occur in all wet years following droughts, leading to permanently lower baseline numbers for some pelagic fishes over time. In contrast, littoral fishes seem to be more resistant to drought and may even increase in occurrence during dry years. Based on the consistent detrimental effects of drought on pelagic fishes within the San Francisco Estuary and the inability of these fish populations to recover in some years, we conclude that freshwater flow remains a crucial but not sufficient management tool for the conservation of estuarine biodiversity.