Geographical gradients of marine herbivorous fishes: patterns and processes

We present new data and the first rigorous analysis of latitudinal and thermal gradients of diversity, density and biomass of marine herbivorous fishes and review proposed explanatory mechanisms. Consistently, negative relationships between latitude, and positive relationships between sea surface temperature (SST), and relative richness and relative abundance of herbivorous fishes were found worldwide. Significant differences in the strength of gradients of richness and abundance with latitude and SST between tropical and extratropical zones were found consistently across ocean basins. Standardized sampling along the western Atlantic also showed negative relationships between latitude and total density and biomass. The trends, however, are driven by different components of the fish assemblages (i.e., scarids in the Caribbean and acanthurids in Brazil). Patterns of abundance along thermal gradients, generally associated with extensive latitudinal gradients, also were found at the local scale. Feeding rate of the ocean surgeonfish Acanthurus bahianus decreases with temperature more rapidly than the mean metabolic rate of teleost fishes. This relationship suggests a temperature-related physiological constraint. From the new standardized and comparative data presented and the review of the explanatory hypotheses, we conclude that temperature-related feeding and digestive processes are most likely involved in the distribution patterns of herbivorous fishes. Electronic Supplementary Material Supplementary material is available for this article at     

Secondary succession dynamics in estuarine marshes across landscape-scale salinity gradients

Secondary succession plays a critical role in driving community structure in natural communities, yet how succession dynamics vary with environmental context is generally unknown. We examined the importance of seedling and vegetative recruitment in the secondary succession of coastal marsh vegetation across a landscape-scale environmental stress gradient. Replicate bare patches were initiated in salt, brackish, and oligohaline marshes in Narragansett Bay, Rhode Island, USA, and allowed to recover unmanipulated or with colonizing seedlings or vegetative runners removed for three years. Seed dispersal and seed bank studies were conducted at the same sites. We found that rates of recovery were 3-10 times faster in brackish and oligohaline marshes than in salt marshes. The fast pace of recovery in oligohaline marshes was driven by seedling colonization, while recovery was dominated by vegetative runners in brackish marshes and by both seedlings and runners in salt marshes. Seed and seedling availability was much greater in oligohaline marshes with up to 24 times the seed bank density compared with salt marshes. In contrast to the facilitated succession generally found in salt marshes, oligohaline marshes follow the tolerance model of succession where numerous species colonize from seed and are slowly displaced by clonal grasses whose recovery is slowed by preemptive competition from seedlings, contributing to the higher species diversity of oligohaline marshes. These findings reveal fundamental differences in the dynamics and assembly of marsh plant communities along estuarine salinity gradients that are important for conceptually understanding wetlands and for guiding the management and restoration of various types of coastal marshes.     

Identification of phosphorus dynamics in Kootenay Lake,Canada

This paper regards the use of system theory techniques in the development of models describing the phosphorus kinetics in fresh waters with the objective to understand the observed kinetics better in order to obtain indications on the design of future field experiments. Identification techniques have been applied to the data regarding 49 experiments on phosphorus cycling in Kootenay Lake, British Columbia, Canada, in order to determine the minimal number of compartments required to model the phosphorus dynamical behaviour. Canonical state-space models have then been identified both from raw and smoothed data. The parameters in these models do not have a direct physical meaning; they can, however, be identified from the measured data and the obtained models describe with good accuracy the considered kinetics. The identification of physical models has also been performed; in this case, however, only two compartments have been considered according to the results of the identifiability analysis of the complete model parameters. Results give further insight on the phosphorus processes and show that the complete phosphorus cycle through the four compartments described by Lean can be identified only when phosphorus cycles rapidly through the system.     

Optimal control of Atlantic population Canada geese

Management of Canada geese (Branta canadensis) can be a balance between providing sustained harvest opportunity while not allowing populations to become overabundant and cause damage. In this paper, we focus on the Atlantic population of Canada geese and use stochastic dynamic programming to determine the optimal harvest strategy over a range of plausible models for population dynamics. There is evidence to suggest that the population exhibits significant age structure, and it is possible to reconstruct age structure from surveys. Consequently the harvest strategy is a function of the age composition, as well as the abundance, of the population. The objective is to maximize harvest while maintaining the number of breeding adults in the population between specified upper and lower limits. In addition, the total harvest capacity is limited and there is uncertainty about the strength of density-dependence. We find that under a density-independent model, harvest is maximized by maintaining the breeding population at the highest acceptable abundance. However if harvest capacity is limited, then the optimal long-term breeding population size is lower than the highest acceptable level, to reduce the risk of the population growing to an unacceptably large size. Under the proposed density-dependent model, harvest is maximized by maintaining the breeding population at an intermediate level between the bounds on acceptable population size; limits to harvest capacity have little effect on the optimal long-term population size. It is clear that the strength of density-dependence and constraints on harvest significantly affect the optimal harvest strategy for this population. Model discrimination might be achieved in the long term, while continuing to meet management goals, by adopting an adaptive management strategy.     

Irruptive population dynamics in Yellowstone pronghorn.

Irruptive population dynamics appear to be widespread in large herbivore populations, but there are few empirical examples from long time series with small measurement error and minimal harvests. We analyzed an 89-year time series of counts and known removals for pronghorn (Antilocapra americana) in Yellowstone National Park of the western United States during 1918-2006 using a suite of density-dependent, density-independent, and irruptive models to determine if the population exhibited irruptive dynamics. Information-theoretic model comparison techniques strongly supported irruptive population dynamics (Leopold model) and density dependence during 1918-1946, with the growth rate slowing after counts exceeded 600 animals. Concerns about sagebrush (Artemisia spp.) degradation led to removals of >1100 pronghorn during 1947-1966, and counts decreased from approximately 700 to 150. The best models for this period (Gompertz, Ricker) suggested that culls replaced intrinsic density-dependent mechanisms. Contrary to expectations, the population did not exhibit enhanced demographic vigor soon after the termination of the harvest program, with counts remaining between 100 and 190 animals during 1967 1981. However, the population irrupted (Caughley model with a one-year lag) to a peak abundance of approximately 600 pronghorn during 1982-1991, with a slowing in growth rate as counts exceeded 500. Numbers crashed to 235 pronghorn during 1992-1995, perhaps because important food resources (e.g., sagebrush) on the winter range were severely diminished by high densities of browsing elk, mule deer, and pronghorn. Pronghorn numbers remained relatively constant during 1996-2006, at a level (196-235) lower than peak abundance, but higher than numbers following the release from culling. The dynamics of this population supported the paradigm that irruption is a fundamental pattern of growth in many populations of large herbivores with high fecundity and delayed density-dependent effects on recruitment when forage and weather conditions become favorable after range expansion or release from harvesting. Incorporating known removals into population models that can describe a wide range of dynamics can greatly improve our interpretation of observed dynamics in intensively managed populations.     

Higher trophic level prey does not represent a higher quality diet in a threatened seabird: implications for relating population dynamics to diet shifts inferred from stable isotopes

Diet quality is a key determinant of population dynamics. If a higher trophic level, more fish-based diet is of higher quality for marine predators, then individuals with a higher trophic level diet should have a greater body mass than those feeding at a lower trophic level. We examined this hypothesis using stable isotope analysis to infer dietary trophic level and foraging habitat over three years in eastern rockhopper penguins Eudyptes chrysocome filholi on sub-Antarctic Campbell Island, New Zealand. Rockhopper penguins are ‘Vulnerable’ to extinction because of widespread and dramatic population declines, perhaps related to nutritional stress caused by a climate-induced shift to a lower trophic level, lower quality diet. We related the stable nitrogen (δ¹⁵N) and carbon (δ¹³C) isotope values of blood from 70 chicks, 55 adult females, and 55 adult males to their body masses in the 2010, 2011, and 2012 breeding seasons and examined year, stage, age, and sex differences. Opposite to predictions, heavier males consumed a lower trophic level diet during incubation in 2011, and average chick mass was heavier in 2011 when chicks were fed a more zooplankton-based, pelagic/offshore diet than in 2012. Contrary to the suggested importance of a fish-based diet, our results support the alternative hypothesis that rockhopper penguin populations are likely to be most successful when abundant zooplankton prey are available. We caution that historic shifts to lower trophic level prey should not be assumed to reflect nutritional stress and a cause of population declines.     

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.     

Fish population dynamics in a seasonally varying wetland

Small fishes in seasonally flooded environments such as the Everglades are capable of spreading into newly flooded areas and building up substantial biomass. Passive drift cannot account for the rapidity of observed population expansions. To test the ‘reaction-diffusion’ mechanism for spread of the fish, we estimated their diffusion coefficient and applied a reaction-diffusion model. This mechanism was also too weak to account for the spatial dynamics. Two other hypotheses were tested through modeling. The first—the ‘refuge mechanism’—hypothesizes that small remnant populations of small fishes survive the dry season in small permanent bodies of water (refugia), sites where the water level is otherwise below the surface. The second mechanism, which we call the ‘dynamic ideal free distribution mechanism’ is that consumption by the fish creates a prey density gradient and that fish taxis along this gradient can lead to rapid population expansion in space. We examined the two alternatives and concluded that although refugia may play an important role in recolonization by the fish population during reflooding, only the second, taxis in the direction of the flooding front, seems capable of matching empirical observations. This study has important implications for management of wetlands, as fish biomass is an essential support of higher trophic levels.     

Long-term population dynamics of seeded plants in invaded grasslands

In recent decades, dozens of studies have involved attempts to introduce native and desirable nonnative plant species into grasslands dominated by invasive weeds. The newly introduced plants have proved capable of establishing, but because they are rarely monitored for more than four years, it is unknown if they have a high likelihood of persisting and suppressing invaders for the long-term. Beyond invaded grasslands, this lack of long-term monitoring is a general problem plaguing efforts to reintroduce a range of taxa into a range of ecosystems. We introduced species from seed and then periodically measured plant abundances for nine years at one site and 15 years at a second site. To our knowledge, our 15-year data are the longest to date from a seeding experiment in invaded, never-cultivated grassland. At one site, three seeded grasses maintained high densities for three or more years, but then all or nearly all individuals died. At the second site, one grass performed similarly, but two other grasses proliferated and at least one greatly suppressed the dominant invader (Centaurea maculosa). In one study, our point estimate suggests that the seeded grass Thinopyrum intermedium reduced C. maculosa biomass by 93% 15 years after seeding. In some cases, data from three and fewer years after seeding falsely suggested that seeded species were capable of persisting within the invaded grassland. In other cases, data from as late as nine years after seeding falsely suggested seeded populations would not become large enough to suppress the invader. These results show that seeded species sometimes persist and suppress invaders for long periods, but short-term data cannot predict if, when, or where this will occur. Because short-term data are not predictive of long-term seeded species performances, additional long-term data are needed to identify effective practices, traits, and species for revegetating invaded grasslands.     

Composite forces shape population dynamics of copepod crustaceans

Understanding the processes that control species abundance and distribution is a major challenge in ecology, yet for a large number of potentially important organisms, we know little about the biotic and abiotic factors that influence population size. One group of aquatic organisms that defies traditional demographic analyses is the Crustacea, particularly those with complex life cycles. We used likelihood techniques and information theoretics to evaluate a suite of models representing alternative hypotheses on factors controlling the abundance of two copepod crustaceans in a small, tropical floodplain lake. Quantitative zooplankton samples were collected at three stations in a Venezuelan floodplain lake from June through December 1984; the average sampling interval was two days. We constructed a series of models with stage structure that incorporated six biotic and abiotic covariates in various combinations to account for temporal changes in abundance of these target species and in their population growth rates. Our analysis produced several novel insights into copepod population dynamics. We found that multiple forces affected the abundance of particular stages, that these factors differed between species as well as among stages within each species, and that biotic processes had the largest effects on copepod population dynamics. Density dependence had a large effect on the survival of Oithona amazonica copepodites and on population growth rate of Diaptomus negrensis.     

Spatiotemporal dynamics in variable population interactions with density-dependent interaction coefficients

In this article, I present a two-patch metapopulation model with locally explicit dynamics to study the effect of spatial heterogeneity and dispersal upon population interactions with variable or conditional outcomes. These are interactions that may be either detrimental or beneficial for each species depending on the balance of the density-dependent costs and benefits involved. The local dynamics respond to density-dependent α-interaction functions that may change sign, thus yielding a diversity of possible local outcomes for the association in terms of type of interaction and in the number of stable solutions. The spatiotemporal model predicts that the fragmentation of space and dispersal between patches may cause further variation in these outcomes. First, the demographic performance of a species in the association is enhanced if migrations cause a proportional increase of individuals of its own species; being so, a victim may become a mutualist or an exploiter, an excluded species may invade, and a good competitor may overcome its own carrying capacity: the ‘enhancement effect of dispersal’; a sort of rescue effect in source-sink dynamics. The underlying mechanisms involve an interplay between density-dependent effects of dispersal per se and the relative local and global average α-interaction functions, which involve costs and benefits at both the local and regional level that may either counteract or reinforce each other; thus, localities and/or populations may change dynamically their sink or source role in the spatial dynamics. A significant insight arises herewith: in the context of variable or conditional interactions the concept of the role of a species does not make strict sense; it becomes a spatiotemporal dynamic quality. Second, regardless of which species disperses, bifurcation of equilibria may occur in those patches that receive the migrating individuals, and annihilation of equilibria in those from where migration leaves; thus, the number of equilibria increases or decreases accordingly.     

Habitat-specific impacts of multiple consumers on plant population dynamics

Multiple consumers often attack seeds, seedlings, and adult plants, but their population-level consequences remain uncertain. We examined how insect and small mammal consumers influenced the demography and abundance of the perennial shrub, bush lupine (Lupinus arboreus). In grassland and dune habitats we established replicate experimental lupine populations in 81-m2 plots that were either protected from, or exposed to, herbivorous voles and granivorous mice (via fencing) and/or root feeding insects (via insecticide treatment). Populations were initiated with transplanted seedlings in 1999 and 2000. We followed the demography of these cohorts, subsequent generations, and the seed bank for 5.5 years. Voles and insects killed many seedlings in dune (1999 only) and grassland (1999 and 2000) habitats. After 2000, insects and voles had minimal effects on seedling or adult survival. Seed predation by granivorous mice, however, greatly depressed seedling recruitment, resulting in lower adult lupine abundance in control plots vs. those protected from rodents. In grasslands, initial effects of voles and insects on seedling survival produced large differences among treatments in adult plant density and the cumulative number of seeds produced throughout the experiment. Differences among grassland populations in seed rain, however, had little influence on the magnitude of seedling recruitment into this habitat. Instead, recruitment out of a preexisting seed bank compensated for the lack of seed production in populations exposed to consumers. Shading by dense adults in plots protected from consumers limited seedling establishment within these populations. Although differences among populations in cumulative seed rain did not influence adult establishment, populations protected from consumers accumulated substantially larger seed banks than controls. These results illustrate how density dependence, habitat-specific seed dynamics, and particular demographic impacts of consumers interact to shape plant population responses to consumers.     

Interacting trophic forcing and the population dynamics of herring

Small pelagic fish occupy a central position in marine ecosystems worldwide, largely by determining the energy transfer from lower trophic levels to predators at the top of the food web, including humans. Population dynamics of small pelagic fish may therefore be regulated neither strictly bottom-up nor top-down, but rather through multiple external and internal drivers. While in many studies single drivers have been identified, potential synergies of multiple factors, as well as their relative importance in regulating population dynamics of small pelagic fish, is a largely unresolved issue. Using a statistical, age-structured modeling approach, we demonstrate the relative importance and influence of bottom-up (e.g., climate, zooplankton availability) and top-down (i.e., fishing and predation) factors on the population dynamics of Bothnian Sea herring (Clupea harengus) throughout its life cycle. Our results indicate significant bottom-up effects of zooplankton and interspecific competition from sprat (Sprattus sprattus), particularly on younger age classes of herring. Although top-down forcing through fishing and predation by grey seals (Halichoerus grypus) and Atlantic cod (Gadus morhua) also was evident, these factors were less important than resource availability and interspecific competition. Understanding key ecological processes and interactions is fundamental to ecosystem-based management practices necessary to promote sustainable exploitation of small pelagic fish.     

Importance of life cycle events in the population dynamics of Gonyaulax tamarensis

Life cycle changes that allow populations of the toxic dinoflagellate Gonyaulax tamarensis Lebour to inhabit the benthos and the plankton alternately are important factors regulationg the initiation and decline of blooms in restricted embavments. When the dynamics of these estuarine populations were monitored during “bloom” and “non bloom” years, it was shown that: (1) each year, germination of benthie cysts inoculated the overlying waters during the vernal warming period, but a large residual population remained in the sediments throughout the blooms; (2) the resulting planktonic population began growth under suboptimal temperature conditions; (3) the populations developed from this inoculum through asexual reproduction until sexuality (and cyst formation) were induced; (4) encystment was not linked to any obvious environmental cue and occurred under apparently optimal conditions; and (5) an increase in the number of non-mitotic swimming cells (planozygotes, the precursors to dormant cysts) accompanied the rapid decline of the planktonic population. Thus encystment, in combination with hypothesized losses due to advection and grazing, contributed substantiatly to the decline of the vegetative cell population. We conclude that the encystment/excystment cycle temporally restricts the occurrence of the vegetative population and may not be optimized for rapid or sustained vegetative growth and bloom formation in shallow embayments. The factors that distinguish “bloom” from “non-bloom” years thus appear to be operating on the growth of the planktonic population.     

Spatial scaling of avian population dynamics: population abundance, growth rate, and variability

Synchrony in population fluctuations has been identified as an important component of population dynamics. In a previous study, we determined that local-scale (<15-km) spatial synchrony of bird populations in New England was correlated with synchronous fluctuations in lepidopteran larvae abundance and with the North Atlantic Oscillation. Here we address five questions that extend the scope of our earlier study using North American Breeding Bird Survey data. First, do bird populations in eastern North America exhibit spatial synchrony in abundances at scales beyond those we have documented previously? Second, does spatial synchrony depend on what population metric is analyzed (e.g., abundance, growth rate, or variability)? Third, is there geographic concordance in where species exhibit synchrony? Fourth, for those species that exhibit significant geographic concordance, are there landscape and habitat variables that contribute to the observed patterns? Fifth, is spatial synchrony affected by a species' life history traits? Significant spatial synchrony was common and its magnitude was dependent on the population metric analyzed. Twenty-four of 29 species examined exhibited significant synchrony in population abundance: mean local autocorrelation (rho)= 0.15; mean spatial extent (mean distance where rho=0) = 420.7 km. Five of the 29 species exhibited significant synchrony in annual population growth rate (mean local autocorrelation = 0.06, mean distance = 457.8 km). Ten of the 29 species exhibited significant synchrony in population abundance variability (mean local autocorrelation = 0.49, mean distance = 413.8 km). Analyses of landscape structure indicated that habitat variables were infrequent contributors to spatial synchrony. Likewise, we detected no effects of life history traits on synchrony in population abundance or growth rate. However, short-distance migrants exhibited more spatially extensive synchrony in population variability than either year-round residents or long-distance migrants. The dissimilarity of the spatial extent of synchrony across species suggests that most populations are not regulated at similar spatial scales. The spatial scale of the population synchrony patterns we describe is likely larger than the actual scale of population regulation, and in turn, the scale of population regulation is undoubtedly larger than the scale of individual ecological requirements.     

How population dynamics shape the functional response in a one-predator-two-prey system

The type III functional response has historically been associated with switching predators; when there is a choice of prey the predator favors the more abundant prey type. Although this functional response has been found in experiments where both prey densities are manipulated, in real world studies the type II functional response is more commonly found. In modeling, the type III functional response is often used in systems where the second prey type is, implicitly, assumed to be constant. Here we define a functional response that takes into account both prey densities. This causes the functional response to show both type II and type III behavior, dependent on the interaction between the two prey densities. If we take into account population dynamics, we find a type II functional response in most cases, because predation regulates the relative prey densities. This explains why type III functional responses are found in experiments where both prey densities are manipulated, but type II functional responses occur when the feedback of population dynamics on the functional response is important. Furthermore, the results show that switching can have a stabilizing or destabilizing effect and can even lead to predator extinction.     

Anticipation and tracking of pulsed resources drive population dynamics in eastern chipmunks

Pulsed systems are characterized by boom and bust cycles of resource production that are expected to cascade through multiple trophic levels. Many of the consumers within pulsed resource systems have specific adaptations to cope with these cycles that may serve to either amplify or dampen their community-wide consequences. We monitored a seed predator, the eastern chipmunk (Tamias striatus), in an American beech (Fagus grandifolia) dominated forest, and used capture-mark-recapture analyses to estimate chipmunk vital rates and relate them to interannual variation in beech seed production. The summer activity and reproduction of adults anticipated autumn beech production, with high activity and intense reproduction occurring in summers prior to beech masts. Chipmunks also reproduced every spring following a beech mast. However, adult survival was independent of beech production. In contrast, juvenile survival was lower in years of mast failure than in years of mast production, but their activity was consistently high and independent of beech production. Population growth was strongly affected by the number of juveniles and therefore by beech seed production, which explains nearly 70% of variation in population growth. Our results suggest that a combination of resource-dependent reproduction and variable activity levels associated with anticipation and response to resource pulses allows consumers to buffer potential deleterious effects of low food abundance on their survival.     

Energy flow and population dynamics of the barnacle Balanus glandula

The energetics and population dynamics of a barnacle (Balanus glandula Darwin) population in British Columbia, Canada, were studied. Consumption, energy flow, production and mortality were 6844.6, 6667.0, 2896.5 and 2522.8 Kcal m^-2 year^-1, respectively. These energy flow and production values are among the highest for animal populations reported, and therefore strongly suggest the functional importance of E. glandula in littoral systems. The young age groups of the first-year settlements were most important in contributing to the energy flow, production and reproduction of the entire population. Most of the assimilated energy in the older age groups was used in respiration.     

Complex dynamics of the population with a simple age structure

The effects of the following modes of density-dependent control of population growth: density-dependent birth rate, adult survival rate, juvenile survival rate are compared based on the mathematical model of population dynamics. It is shown that the most efficient mechanisms limiting population size are decreasing with the growth of the adult population birth rate and/or the decreasing survival rate of the offspring with the increase in their number. However, these same mechanisms are responsible for oscillations of the population size and its chaotic change. The density-dependence of the adult survival rate is not efficient in constraining the population growth, but it can substantially limit the magnitude of oscillations of the population size.     

Biology and population dynamics of the intertidal isopod Cirolana harfordi

The distribution of Cirolana harfordi (Lockington) populations is determined largely by the availability of loose boulders on sandy beaches. The isopods swim out from under the rocks at high tide to feed. Their diet consists primarily of minute polychaetes and crustaceans. In addition, the isopods locate and utilize any available dead animal matter in the surf zone. Breeding occurs throughout most of the year, except for a brief lull in the fall. Females produce 1 or 2 broods of 18 to 68 young during their 2-year life-span. Marsupial incubation lasts 3 to 4 months. The population size and structure remained relatively constant from one year to the next. Mortality rates estimated for juveniles, males and females showed that newly emergent young and post-reproductive females suffer the greatest losses (up to 75% mortality per month).