Characterizing cumulative impacts using a brook trout population dynamics model

An individuals-based modelling framework is used to characterize the nature of exploitation and toxaphene stressors acting simultaneously on a population of brook trout (Salvelinus fontinalis) in terms of age 0 + and adult abundance, survivorship and population size-structure. A no-stressor control case was estimated against which exploitation-only, toxaphene-only and cumulative exploitation and toxaphene stressor cases were compared to determine the extent and significance of impacts. Single stressor case results were used to predict cumulative impacts by assuming additivity and predictions compared to modelled cumulative stressor results. Comparisons indicated the inadequacy of assuming additivity in predicting cumulative impacts. A factorial experimental design was used to estimate the size and significance of interactive effects. Effects are substantial and underscore the necessity of interpreting probable impacts of increases in a single stressor in conjunction with knowledge of other stressors acting on a population. A positive functional relationship between variability in population abundance and stress was estimated and is suggested as a potentially useful means of characterizing risks posed to populations by increases in, or additions to, population stressors. Multiple stressors were also demonstrated to effectively eliminate the significance of density-dependent mortality processes in determining age 0 + and adult abundance. Taken together, results indicate the inappropriateness of attempting to predict additional perturbation impacts without considering the sum of population stressors and their associated interactions.     

SPAnDX: a process-based population dynamics model to explore management and climate change impacts on an invasive alien plant,Acacia nilotica

This paper describes a process-based metapopulation dynamics and phenology model of prickly acacia, Acacia nilotica, an invasive alien species in Australia. The model, SPAnDX, describes the interactions between riparian and upland sub-populations of A. nilotica within livestock paddocks, including the effects of extrinsic factors such as temperature, soil moisture availability and atmospheric concentrations of carbon dioxide. The model includes the effects of management events such as changing the livestock species or stocking rate, applying fire, and herbicide application. The predicted population behaviour of A. nilotica was sensitive to climate. Using 35 years daily weather datasets for five representative sites spanning the range of conditions that A. nilotica is found in Australia, the model predicted biomass levels that closely accord with expected values at each site. SPAnDX can be used as a decision-support tool in integrated weed management, and to explore the sensitivity of cultural management practices to climate change throughout the range of A. nilotica. The cohort-based DYMEX modelling package used to build and run SPAnDX provided several advantages over more traditional population modelling approaches (e.g. an appropriate specific formalism (discrete time, cohort-based, process-oriented), user-friendly graphical environment, extensible library of reusable components, and useful and flexible input/output support framework).     

Recolonizing wolves and mesopredator suppression of coyotes: impacts on pronghorn population dynamics.

Food web theory predicts that the loss of large carnivores may contribute to elevated predation rates and, hence, declining prey populations, through the process of mesopredator release. However, opportunities to test predictions of the mesopredator release hypothesis are rare, and the extent to which changes in predation rates influence prey population dynamics may not be clear due to a lack of demographic information on the prey population of interest. We utilized spatial and seasonal heterogeneity in wolf distribution and abundance to evaluate whether mesopredator release of coyotes (Canis latrans), resulting from the extirpation of wolves (Canis lupus) throughout much of the United States, contributes to high rates of neonatal mortality in ungulates. To test this hypothesis, we contrasted causes of mortality and survival rates of pronghorn (Antilocapra americana) neonates captured at wolf-free and wolf-abundant sites in western Wyoming, USA, between 2002 and 2004. We then used these data to parameterize stochastic population models to heuristically assess the impact of wolves on pronghorn population dynamics due to changes in neonatal survival. Coyote predation was the primary cause of mortality at all sites, but mortality due to coyotes was 34% lower in areas utilized by wolves (P < 0.001). Based on simulation modeling, the realized population growth rate was 0.92 based on fawn survival in the absence of wolves, and 1.06 at sites utilized by wolves. Thus, wolf restoration is predicted to shift the trajectory of the pronghorn population from a declining to an increasing trend. Our results suggest that reintroductions of large carnivores may influence biodiversity through effects on prey populations mediated by mesopredator suppression. In addition, our approach, which combines empirical data on the population of interest with information from other data sources, demonstrates the utility of using simulation modeling to more fully evaluate ecological theories by moving beyond estimating changes in vital rates to analyses of population-level impacts.     

Seasonal timing of first rain storms affects rare plant population dynamics

A major challenge in forecasting the ecological consequences of climate change is understanding the relative importance of changes to mean conditions vs. changes to discrete climatic events, such as storms, frosts, or droughts. Here we show that the first major storm of the growing season strongly influences the population dynamics of three rare and endangered annual plant species in a coastal California (USA) ecosystem. In a field experiment we used moisture barriers and water addition to manipulate the timing and temperature associated with first major rains of the season. The three focal species showed two- to fivefold variation in per capita population growth rates between the different storm treatments, comparable to variation found in a prior experiment imposing eightfold differences in season-long precipitation. Variation in germination was a major demographic driver of how two of three species responded to the first rains. For one of these species, the timing of the storm was the most critical determinant of its germination, while the other showed enhanced germination with colder storm temperatures. The role of temperature was further supported by laboratory trials showing enhanced germination in cooler treatments. Our work suggests that, because of species-specific cues for demographic transitions such as germination, changes to discrete climate events may be as, if not more, important than changes to season-long variables.     

Native and exotic plant species exhibit similar population dynamics during succession

A growing body of literature has led to the debate in invasion biology whether exotic species perform within communities differently than native taxa due to inherent advantages. To address this issue, the population dynamics of native and exotic plant species were assessed from a 48-year record of permanent plot data from the Hutcheson Memorial Forest Center (New Jersey, USA) to determine rate of increase, lag time, maximum frequency, and the year of peak frequency. Overall, native and exotic species exhibited very similar population dynamics. Rates of increase and length of lag times were similar between native and exotic taxa but were strongly influenced by plant life form. Short-lived species were characterized by rapid population growth rates and short lag times. Growth rates decreased and lag times increased with species longevity. Overall, correlations between population metrics were the same in native and exotic taxa, suggesting similar trade-offs in life history patterns. The one difference observed was that, in native species, peak frequency was negatively associated with the year of peak frequency (i.e., early-successional species tended to become more abundant), while there was no relationship in exotic species. These analyses show that exotic species behave in essentially the same way as native taxa within dynamic communities. This suggests that abundant native and exotic plant species are exploiting the same range of ecological strategies resulting in similar roles within communities.     

Modeling the impacts of hunting on the population dynamics of red howler monkeys (Alouatta seniculus)

Overexploitation of wildlife populations occurs across the humid tropics and is a significant threat to the long-term survival of large-bodied primates. To investigate the impacts of hunting on primates and ways to mitigate them, we developed a spatially explicit, individual-based model for a landscape that included hunted and un-hunted areas. We used the large-bodied neotropical red howler monkey (Alouatta seniculus) as our case study species because its life history characteristics make it vulnerable to hunting. We modeled the influence of different rates of harvest and proportions of landscape dedicated to un-hunted reserves on population persistence, population size, social dynamics, and hunting yields of red howler monkeys. In most scenarios, the un-hunted populations maintained a constant density regardless of hunting pressure elsewhere, and allowed the overall population to persist. Therefore, the overall population was quite resilient to extinction; only in scenarios without any un-hunted areas did the population go extinct. However, the total and hunted populations did experience large declines over 100 years under moderate and high hunting pressure. In addition, when reserve area decreased, population losses and losses per unit area increased disproportionately. Furthermore, hunting disrupted the social structure of troops. The number of male turnovers and infanticides increased in hunted populations, while birth rates decreased and exacerbated population losses due to hunting. Finally, our results indicated that when more than 55% of the landscape was harvested at high (30%) rates, hunting yields, as measured by kilograms of biomass, were less than those obtained from moderate harvest rates. Additionally, hunting yields, expressed as the number of individuals hunted/year/km², increased in proximity to un-hunted areas, and suggested that dispersal from un-hunted areas may have contributed to hunting sustainability. These results indicate that un-hunted areas serve to enhance hunting yields, population size, and population persistence in hunted landscapes. Therefore, spatial regulation of hunting via a reserve system may be an effective management strategy for sustainable hunting, and we recommend it because it may also be more feasible to implement than harvest quotas or restrictions on season length.     

A process-based population dynamics model to explore target and non-target impacts of a biological control agent

The risks and benefits associated with efforts to control invasive alien species using classical biological control are being subjected to increasing scrutiny. A process-based population dynamics model was developed to explore the interactions between a folivorous biological control agent, Cleopus japonicus, and its plant host Buddleja davidii. The model revealed that climate could have a significant impact upon the interactions between B. davidii and C. japonicus. At the coolest sites, the impact of C. japonicus on B. davidii was slowed, but it was still eventually capable of controlling populations of B. davidii. At the warmer sites where both B. davidii and C. japonicus grew faster, B. davidii succumbed rapidly to weevil damage. We hypothesise that barring an encounter with a natural enemy, C. japonicus will eventually be able to provide sustained control B. davidii throughout the North Island of New Zealand. The model scenarios illustrate the potential for the C. japonicus population to attain high densities rapidly, and to defoliate patches of B. davidii, creating the potential for spill-over feeding on non-target plants. The potential magnitude of this threat will depend partly on the climate suitability for C. japonicus, the pattern by which it migrates in response to a reduction in the available leaf resource, and the suitability of non-target plants as hosts. In all migration scenarios considered, the pattern of population growth and resource consumption by C. japonicus was exponential, with a strong tendency toward complete utilisation of resource patches more quickly at the warmer compared to colder sites. In addition to providing some useful hypotheses about the effects of climate on the biological control system, and the non-target risks, it also provides some insight into the mechanisms by which climate affects the system.     

Estimating population dynamics without population data

We develop a biologically correct cost system for production systems facing invasive pests that allows the estimation of population dynamics without a priori knowledge of their true values. We apply that model to a data set for olive producers in Crete and derive from it predictions about the underlying population dynamics. Those dynamics are compared to information on population dynamics obtained from pest sampling with extremely favorable results.     

Estimating population dynamics without population data

We develop a biologically correct cost system for production systems facing invasive pests that allows the estimation of population dynamics without a priori knowledge of their true values. We apply that model to a data set for olive producers in Crete and derive from it predictions about the underlying population dynamics. Those dynamics are compared to information on population dynamics obtained from pest sampling with extremely favorable results.     

Temperature-associated population diversity in salmon confers benefits to mobile consumers

Habitat heterogeneity can generate intraspecific diversity through local adaptation of populations. While it is becoming increasingly clear that population diversity can increase stability in species abundance, less is known about how population diversity can benefit consumers that can integrate across population diversity in their prey. Here we demonstrate cascading effects of thermal heterogeneity on trout-salmon interactions in streams where rainbow trout rely heavily on the seasonal availability of anadromous salmon eggs. Water temperature in an Alaskan stream varied spatially from 5 degrees C to 17.5 degrees C, and spawning sockeye salmon showed population differentiation associated with this thermal heterogeneity. Individuals that spawned early in cool regions of the 5 km long stream were genetically differentiated from those spawning in warmer regions later in the season. Sockeye salmon spawning generates a pulsed resource subsidy that supports the majority of seasonal growth in stream-dwelling rainbow trout. The spatial and temporal structuring of sockeye salmon spawn timing in our focal stream extended the duration of the pulsed subsidy compared to a thermally homogeneous stream with a single population of salmon. Further, rainbow trout adopted movement strategies that exploited the multiple pulses of egg subsidies in the thermally heterogeneous stream. Fish that moved to track the resource pulse grew at rates about 2.5 times higher than those that remained stationary or trout in the reference stream with a single seasonal pulse of eggs. Our results demonstrate that habitat heterogeneity can have important effects on the population diversity of dominant species, and in turn, influence their value to species that prey upon them. Therefore, habitat homogenization may have farther-reaching ecological effects than previously considered.     

Effects of nanoparticles on Daphnia magna population dynamics

Spherical TiO₂ nanoparticles (npTiO₂) were prepared by controlled hydrolysis of tetraethoxy orthotitanate under a nitrogen atmosphere. ZnO nanoparticles (npZnO) were prepared using hydrothermal methods. The crystal structure, chemical, thermal and morphological properties of npZnO and npTiO₂ were characterised using Fourier Transform Infrared Spectrometer, enery-dispersive X-ray spectroscopy, X-ray diffraction, and scanning electron microscope techniques. The short- and long-term experiments were started with neonates taken from the same culture and laboratory condition. In the acute experiments, npTiO₂, npZnO, and cocktail concentrations were applied. 96h-LC₅₀ values were 1.8, 0.7, and 0.1?mg?L^?1, respectively (p?<?.05). For the chronic experiments, different npTiO₂ concentrations were performed. 21d-LC₅₀ value was 1.0?mgL^?1 (p?<?.05). Morphometry became progressively worse in concentrations of more than 1?mgL^?1 npTiO₂. Neonate and young individuals were more sensitive to death because of their low tolerance. This result was affected by population progeny and growth rates (p?<?.05). While control and 0.5?mgL^?1 npTiO₂ groups were determined as growing population, 1.5 and 2?mgL^?1 npTiO₂ groups had decreased population size as R₀ values. Consequently, the relationships between nanoparticle accumulation within Daphnia magna and its population structure and body morphometry for each concentration were important indicators. Its tolerance level to nanoparticles under laboratory conditions reflected its replacement and behaviour in the ecosystem.     

Processing of allochthonous macrophyte subsidies by sandy beach consumers: estimates of feeding rates and impacts on food resources

Allochthonous subsidies of organic material can profoundly influence population and community structure; however, the role of consumers in the processing of these inputs is less understood but may be closely linked to community and ecosystem function. Inputs of drift macrophytes subsidize sandy beach communities and food webs in many regions. We estimated feeding rates of dominant sandy beach consumers, the talitrid amphipods (Megalorchestia corniculata, in southern California, USA, and Talitrus saltator, in southern Galicia, Spain), and their impacts on drift macrophyte subsidies in field and laboratory experiments. Feeding rate varied with macrophyte type and, for T. saltator, air temperature. Size-specific feeding rates of talitrid amphipods were greatest on brown macroalgae (Macrocystis, Egregia, Saccorhiza and Fucus). Rates for large individuals of both species ranged from ∼40 mg wet wt individual⁻¹12 h⁻¹ on brown macroalgae to negligible feeding by M. corniculata on a vascular plant (surfgrass). Amphipod growth rates were also greatest on Macrocystis and lowest on surfgrass, Phyllospadix. For a Californian beach with substantial inputs of macrophyte wrack (>70 kg wet wt m⁻¹ month⁻¹ in summer), we estimated that the population of talitrid amphipods could process an average of 55% of the palatable Macrocystis input. Our results indicate that talitrid amphipod populations can have a significant impact on drift macrophyte processing and fate and that the quantity and composition of drift macrophytes could, in turn, limit populations of beach consumers.     

Uncertainty analysis of transient population dynamics

Two types of demographic analyses, perturbation analysis and uncertainty analysis, can be conducted to gain insights about matrix population models and guide population management. Perturbation analysis studies how the perturbation of demographic parameters (survival, growth, and reproduction parameters) may affect the population projection, while uncertainty analysis evaluates how much uncertainty there is in population dynamic predictions and where the uncertainty comes from. Previously, both perturbation analysis and uncertainty analysis were conducted on the long-term population growth rate. However, the population may not reach its equilibrium state, especially when there is management by harvesting or hunting. Recently, there has been an increased interest in short-term transient dynamics, which can differ from asymptotic long-term dynamics. There are currently techniques to conduct perturbation analyses of short-term transient dynamics, but no techniques have been proposed for uncertainty analysis of such dynamics. In this study, we introduced an uncertainty analysis technique, the general Fourier Amplitude Sensitivity Test (FAST), to study uncertainties in transient population dynamics. The general FAST is able to identify the amount of uncertainty in transient dynamics and contributions by different demographic parameters. We applied the general FAST to a mountain goat (Oreamnos americanus) matrix population model to give a clear illustration of how uncertainty analysis can be conducted for transient dynamics arising from matrix population models.     

Human impacts, plant invasion, and imperiled plant species in California.

Invasive species are one of the fastest growing conservation problems. These species homogenize the world's flora and fauna, threaten rare and endemic species, and impose large economic costs. Here, we examine the distribution of 834 of the more than 1000 exotic plant taxa that have become established in California, USA. Total species richness increases with net primary productivity; however, the exotic flora is richest in low-lying coastal sites that harbor large numbers of imperiled species, while native diversity is highest in areas with high mean elevation. Weedy and invasive exotics are more tightly linked to the distribution of imperiled species than the overall pool of exotic species. Structural equation modeling suggests that while human activities, such as urbanization and agriculture, facilitate the initial invasion by exotic plants, exotics spread ahead of the front of human development into areas with high numbers of threatened native plants. The range sizes of exotic taxa are an order of magnitude smaller than for comparable native taxa. The current small range size of exotic species implies that California has a significant "invasion debt" that will be paid as exotic plants expand their range and spread throughout the state.     

Evaluating the impacts of multiple generalist fungal pathogens on temperate tree seedling survival

Host-specific mortality driven by natural enemies is a widely discussed mechanism for explaining plant diversity. In principle, populations of plant species can be regulated by distinct host-specific natural enemies that have weak or nonexistent effects on heterospecific competitors, preventing any single species from becoming dominant and thus promoting diversity. Two of the first steps in exploring the role of natural enemies in diversity regulation are to (1) identify potential enemies and (2) evaluate their levels of host specificity by determining if interactions between any one host and its enemy have equivalent survival impacts on co-occurring host species. We developed a bioinformatics framework to evaluate impacts of potential pathogens on seedling survival, for both single and multiple infections. Importantly, we consider scenarios not only if there are specialist pathogens for each plant, but also when generalist pathogens have differential effects on multiple host species, and when co-infection has species-specific effects. We then applied this analytical framework to a field experiment using molecular techniques to detect potential fungal pathogens on co-occurring tree seedling hosts. Combinatorial complexity created by 160 plant-fungus interactions was reduced to eight combinations that affect seedling survival. Potential fungal pathogens had broad host ranges, but seedling species were each regulated by different combinations of fungi or by generalist fungi that had differential effects on multiple plant species. Soil moisture can have the potential to shift the nature of the interactions in some plant-fungal combinations from neutral to detrimental. Reassessing the assumption of single-enemy-single-host interactions broadens the mechanisms through which natural enemies can influence plant diversity.     

Ungulate population dynamics and optimization models

Optimal harvesting strategies for an ungulate population are estimated using stochastic dynamic programming. Data on the Llano Basin white-tailed deer (Odocoileus virginianus) population were used to construct a 2-variable population dynamics model. The model provided the basis for estimating optimal harvesting strategies as a feedback function of the current values of the state variables (prefawning older deer and juveniles). Optimal harvest strategies were insensitive to assumptions about the probability distributions of the stochastic variable (rainfall). The response of the population components to harvesting and the returns obtained from applying optimal strategies were explored through simulation. Mean annual harvest is about 15% of the population. Simplified harvesting strategies based on age-ratios as well as a simplified version based on optimal strategies—but assuming persisting equilibrium juvenile deer density—were compared to optimal strategies through examining values of information. Simplified harvesting strategies lead to a lower harvest over a 50-year simulation period.     

A simulation model to evaluate the impacts of invasive earthworms on soil carbon dynamics

Recent studies have reported that earthworm invasions alter native communities and impact nutrient cycling in terrestrial ecosystems. We developed a simulation model to evaluate the potential impacts of earthworm invasions on carbon dynamics, taking into consideration earthworm feeding strategies and priming effects on the microorganisms through their casting activities. Responses of carbon stocks (forest litter, soil organic matter, microbial biomass and earthworm populations) and carbon fluxes (litter decomposition, earthworm consumption, and microbial respiration) were used to evaluate an earthworm invasion of a forest ecosystem. Data from a northern temperate forest (Arnot Forest, New York) were adapted for model calibration and evaluation. Simulation results suggest that the impact and outcome of earthworm invasions are affected by pre-invasion resource availability (litter and soil organic matter), invasive earthworm assemblages (particularly feeding strategy), and invasion history (associated with earthworm population dynamics). The abovementioned factors may also determine invasion progress of earthworm species. The accuracy of the model could be improved by the addition of environmental modules (e.g., soil water regimes), precise parameters accounting for individual species attributes under different environmental conditions (e.g. utilization ability of different types of food resources), as well as earthworm population dynamics (size and structure) and interactions with predators and other invasive/indigenous species during the invasion progress. Such an earthworm invasion model could provide valuable evaluation of the complicated responses of carbon dynamics to earthworm invasions in a range of forest ecosystems, particularly under global change scenarios.     

Estimating population impacts via dynamic occupancy analysis of Before-After Control-Impact studies

Estimating environmental impacts on populations is one of the main goals of wildlife monitoring programs, which are often conducted in conjunction with management actions or following natural disturbances. In this study we investigate the statistical power of dynamic occupancy models to detect changes in local survival and colonization from detection-nondetection data, while accounting for imperfect detection probability, in a Before-After Control-Impact (BACI) framework. We simulated impacts on local survival and/or detection probabilities, and asked questions related to: (1) costs and benefits of different analysis models, (2) confounding changes in detection with changes in local survival, (3) sampling design trade-offs, and (4) species with low vs. high rates of turnover. Estimating seasonal effects on local survival and colonization, as opposed to estimating Before-After effects, had little effect on the power to detect changes in local survival. Estimating a parameter that accounted for pretreatment differences in local survival between Control and Impact sites decreased power by 50%, but it was critical to include when such differences existed. When the experimental treatment had a negative impact on species detectability but analysis assumed constant detection, the Type I error rates were dramatically inflated (0.20 0.33). In general, there was low power (< 0.5) to detect a 50% decrease in local survival for all combinations of sites (N = 50 vs. 100), seasons sampled (8 vs. 12), and visits per site per season (4 vs. 6). Unbalanced designs performed worse than balanced designs, with the exception of the case of treatments being implemented in different seasons at different sites. Adding more control sites improved the ability to detect changes in local survival. Surveying more seasons after impact resulted in modest power gains, but at least three seasons before impact were required to successfully implement BACI occupancy studies. Turnover rates had a low impact on power. Occupancy studies conducted in a BACI design offer the opportunity to detect environmental impacts on wildlife populations without the costs of intensive studies. However, given the low power to detect small changes (20%) in local survival, these studies should be used when researchers are confident that major treatment impacts will occur or very large sample sizes are obtainable.     

Simulation of herbicide impacts on a plant community: comparing model predictions of the plant community model IBC-grass to empirical data

Background

Semi-natural plant communities such as field boundaries play an important ecological role in agricultural landscapes, e.g., provision of refuge for plant and other species, food web support or habitat connectivity. To prevent undesired effects of herbicide applications on these communities and their structure, the registration and application are regulated by risk assessment schemes in many industrialized countries. Standardized individual-level greenhouse experiments are conducted on a selection of crop and wild plant species to characterize the effects of herbicide loads potentially reaching off-field areas on non-target plants. Uncertainties regarding the protectiveness of such approaches to risk assessment might be addressed by assessment factors that are often under discussion. As an alternative approach, plant community models can be used to predict potential effects on plant communities of interest based on extrapolation of the individual-level effects measured in the standardized greenhouse experiments. In this study, we analyzed the reliability and adequacy of the plant community model IBC-grass (individual-based plant community model for grasslands) by comparing model predictions with empirically measured effects at the plant community level.

Results

We showed that the effects predicted by the model IBC-grass were in accordance with the empirical data. Based on the species-specific dose responses (calculated from empirical effects in monocultures measured 4 weeks after application), the model was able to realistically predict short-term herbicide impacts on communities when compared to empirical data.

Conclusion

The results presented in this study demonstrate an approach how the current standard greenhouse experiments—measuring herbicide impacts on individual-level—can be coupled with the model IBC-grass to estimate effects on plant community level. In this way, it can be used as a tool in ecological risk assessment.