Predation, habitat complexity, and variation in density-dependent mortality of temperate reef fishes

Density dependence in demographic rates can strongly affect the dynamics of populations. However, the mechanisms generating density dependence (e.g., predation) are also dynamic processes and may be influenced by local conditions. Understanding the manner in which local habitat features affect the occurrence and/or strength of density dependence will increase our understanding of population dynamics in heterogeneous environments. In this study I conducted two separate field experiments to investigate how local predator density and habitat complexity affect the occurrence and form of density-dependent mortality of juvenile rockfishes (Sebastes spp.). I also used yearly censuses of rockfish populations on nearshore reefs throughout central California to evaluate mortality of juvenile rockfish at large spatial scales. Manipulations of predators (juvenile bocaccio, S. paucispinus) and prey (kelp, gopher, and black-and-yellow [KGB] rockfish, Sebastes spp.) demonstrated that increasing the density of predators altered their functional response and thus altered patterns of density dependence in mortality of their prey. At low densities of predators, the number of prey consumed per predator was a decelerating function, and mortality of prey was inversely density dependent. However, at high densities of predators, the number of prey killed per predator became an accelerating response, and prey mortality was directly density dependent. Results of field experiments and large-scale surveys both indicated that the strength of density-dependent mortality may also be affected by the structural complexity of the habitat. In small-scale field experiments, increased habitat complexity increased the strength of density-dependent mortality. However, at large scales, increasing complexity resulted in a decrease in the strength of density dependence. I suggest that these differences resulted from scale-dependent changes in the predatory response that generated mortality. Whether increased habitat complexity leads to an increase or a decrease in the strength of density-dependent mortality may depend on how specific predatory responses (e.g., functional or aggregative) are altered by habitat complexity. Overall, the findings of this study suggest that rates of demographic density dependence and the resulting dynamics of local populations may largely depend upon attributes of the local habitat.     

Synergistic influences of phase, density, and climatic variation on the dynamics of fluctuating populations

Although ecologists have long recognized that certain mammalian species exhibit high-amplitude, often multiannual, fluctuations in abundance, their causes have remained poorly understood and the subject of intense debate. A key contention has been the relative role of density-dependent and density-independent processes in governing population dynamics. We applied capture-mark-recapture analysis to 25 years of monthly trapping data from a fluctuating prairie vole Microtus ochrogaster population in Illinois, USA, to estimate realized population growth rates and associated vital rates (survival and recruitment) and modeled them as a function of vole density and density-independent climatic variation. We also tested for phase dependence and seasonality in the effects of the above processes. Variation in the realized population growth rate was best explained by phase-specific changes in vole density lagged by one month and mean monthly temperatures with no time lags. The underlying vital rates, survival and recruitment, were influenced by the additive and interactive effects of phase, vole density, and mean monthly temperatures. Our results are consistent with the observation that large-scale population fluctuations are characterized by phase-specific changes in demographic and physiological characteristics. Our findings also support the growing realization that the interaction between climatic variables and density-dependent factors may be a widespread phenomenon, and they suggest that the direction and magnitude of such interactive effects may be phase specific. We conclude that density-dependent and density-independent climatic variables work in tandem during each phase of density fluctuations to drive the dynamics of fluctuating populations.     

The interplay of sex ratio, male success and density-independent mortality affects population dynamics

Environmental constraints can limit a population to a certain size, which is usually called the carrying capacity of a habitat. Besides to this ‘external’ factor, which is mainly determined by the limitation of resources, we investigate here another set of population-intrinsic factors that can limit a population size significantly below the maximum sustainable size. Firstly, density-independent mortality is a prominent factor in all organisms that show age-related and/or accidental death. Secondly, in sexually reproducing organisms the sex ratio and the success of pairing is important for finding reproductive partners. Using a simple model, we demonstrate how sex ratio, mating success and gender-specific mortality can strongly affect the speed of population growth and the maximum population size. In addition, we demonstrate that density-independent mortality, which is often neglected in population models, adds a very important feature to the system: it strongly enhances the negative influence of unbiased sex ratios and inefficient pairing to the maximum sustainable population size. A decrease of the maximum population size significantly affects a population's survival chance in inter-specific competition. Thus, we conclude that the inclusion of density-independent mortality is crucial, especially for models of species that reproduce sexually. We show that density-independent mortality, together with biased sex ratios, can significantly lower the abilities of a population to survive in conditions of strong inter-specific competition and due to the Allee effect. We emphasize that population models should incorporate the sex ratio, male success and density-independent mortality to make plausible predictions of the population dynamics in a gender-structured population. We show that the population size is limited by these intrinsic factors. This is of high ecological significance, because it means that there will always be resources available in any habitat that allows other species (e.g., invaders) to use these resources and settle successfully, if they are sufficiently adapted.     

Ontogenetic and sex-specific differences in density-dependent habitat selection of a marine fish population

The spatial dynamics of species are the result of complex interactions between density-independent and density-dependent sources of variability. Disentangling these two sources of variability has challenged ecologists working in both terrestrial and aquatic ecosystems. Using a novel spatially explicit statistical model, we tested for the presence of density-independent and density-dependent habitat selection in yellowfin sole (Limanda aspera) in the eastern Bering Sea. We found specificities in the density-dependent processes operating across ontogeny and particularly with gender. Density-dependent habitat expansion occurred primarily in females, and to a lesser degree in males. These patterns were especially evident in adult stages, while juvenile stages of both sexes exhibited a mix of different dynamics. Association of yellowfin sole with substrate type also varied by sex and to a lesser degree with size, with large females distributed over a wider range of substrates than males. Moreover, yellowfin sole expanded northward as cold subsurface waters retracted in summer, suggesting high sensitivity to arctic warming. Our findings illustrate how marginal habitats can play an important role in buffering density-dependent habitat expansion, with direct implications for resource management. Our spatially explicit modeling approach is effective in evaluating density-dependent spatial dynamics, and can easily be used to test similar hypotheses from a variety of aquatic and terrestrial ecosystems.     

Vertebrate herbivory impacts seedling recruitment more than niche partitioning or density-dependent mortality

In tropical forests, resource-based niches and density-dependent mortality are mutually compatible mechanisms that can act simultaneously to limit seedling populations. Differences in the strengths of these mechanisms will determine their roles in maintaining species coexistence. In the first assessment of these mechanisms in a Congo Basin forest, we quantified their relative strengths and tested the extent to which density-dependent mortality is driven by the distance-dependent behavior of seed and seedling predators predicted by the Janzen-Connell hypothesis. We conducted a large-scale seed addition experiment for five randomly selected tropical tree species, caging a subset of seed addition quadrats against vertebrate predators. We then developed models to assess the mechanisms that determine seedling emergence (three months after seed addition) and survival (two years after seed addition). As predicted, both niche differentiation and density-dependent mortality limited seedling recruitment, but predation had the strongest effects on seedling emergence and survival. Seedling species responded differently to naturally occurring environmental variation among sites, including variation in light levels and soil characteristics, supporting predictions of niche-based theories of tropical tree species coexistence. The addition of higher densities of seeds into quadrats initially led to greater seedling emergence, but survival to two years decreased with seed density. Seed and seedling predation reduced recruitment below levels maintained by density-dependent mortality, an indication that predators largely determine the population size of tree seedlings. Seedling recruitment was unrelated to the distance to or density of conspecific adult trees, suggesting that recruitment patterns are generated by generalist vertebrate herbivores rather than the specialized predators predicted by the Janzen-Connell hypothesis. If the role of seed and seedling predation in limiting seedling recruitment is a general phenomenon, then the relative abundances of tree species might largely depend on species-specific adaptations to avoid, survive, and recover from damage induced by vertebrate herbivores. Likewise, population declines of herbivorous vertebrate species (many of which are large and hunted) may trigger shifts in species composition of tropical forests.     

Consequences of patch reef spacing for density-dependent mortality of coral-reef fishes

The spatial configuration of habitat patches can profoundly affect a number of ecological interactions, including those between predators and prey. I examined the effects of reef spacing on predator-prey interactions within coral-reef fish assemblages in the Bahamas. Using manipulative field experiments, I determined that reef spacing influences whether and how density-dependent predation occurs. Mortality rates of juveniles of two ecologically dissimilar species (beaugregory damselfish and yellowhead wrasse) were similarly affected by reef spacing; for both species, mortality was density dependent on reef patches that were spatially isolated (separated by 50 m), and density independent on reef patches that were aggregated (separated by 5 m). A subsequent experiment with the damselfish demonstrated that a common resident predator (coney) caused a substantial proportion of the observed mortality, independent of reef spacing. Compared to isolated reefs, aggregated reefs were much more likely to be visited by transient predators (mostly yellowtail snappers), regardless of prey density, and on these reefs, mortality rates approached 100% for both prey species. Transient predators exhibited neither an aggregative response nor a type 3 functional response, and consequently were not the source of density dependence observed on the isolated reefs. These patterns suggest that resident predators caused density-dependent mortality in their prey through type 3 functional responses on all reefs, but on aggregated reefs, this density dependence was overwhelmed by high, density-independent mortality caused by transient predators. Thus, the spatial configuration of reef habitat affected both the magnitude of total predation and the existence of density-dependent mortality. The combined effects of the increasing fragmentation of coral reef habitats at numerous scales and global declines in predatory fish may have important consequences for the regulation of resident fish populations.     

Abundances of surgeonfishes on patch-reefs in Caribbean Panamá: due to settlement,or post-settlement events?

I examined whether patterns of settlement of planktonic juveniles determined the size and composition of adult populations of three Caribbean surgeonfishes (A canthurus bahianus, A. chirurgus and A. coeruleus) on six isolated patch-reefs in Panamá. I monitored settlement monthly on six reefs (0.14 to 0.67 ha) for 8 yr (1978–1986), and censused adults on those reefs annually. Year-to-year fluctuations in the absolute and relative abundances of settlers arriving on each reef were not extreme. Densities of adults in 1986 were not correlated with the densities of settlers that arrived between 1979 and 1984. Survival of settlers, although apparently unrelated to the density of adults, was negatively related to the density of settlers per unit of adult habitat (but not per unit of settlement habitat) for each species. Changes in adult populations of each species between 1978 and 1986 were not correlated with the abundance of settlers relative to the initial number of adults. The relative abundances of adults of each species in 1986 did not reflect the relative abundances of settlers (1979–1984); these abundances were due to both spatial variation and interspecific variation in survival. In some cases, relocation of adults between isolated patch-reefs (as much as 100 m from other reefs) substantially altered the size and composition of surgeonfish populations. Thus, either post-settlement mortality or relocation overrode effects of settlement patterns in determining adult abundances on the six reefs. Relocation is important in organizing reef-fish communities and must be further examined, particularly in large habitat mosaics.     

Modeling the Effects of Fishing and Implications for the Design of Marine Protected Areas: Juvenile Fish Responses to Variations in Seafloor Habitat

Abstract: A number of recent studies have linked post-settlement survivorship of Atlantic cod (  Gadus morhua ) with the complexity of the seafloor to which fish settle. Survivorship is greater in habitats of higher complexity (e.g., pebble-cobble substratum with emergent epifauna> pebble-cobble> sand), where cover provides shelter from predators. Fishing with mobile gear such as bottom trawls and dredges reduces the complexity of seafloor habitats. We used a dynamic model to (1) link patterns in habitat-mediated survivorship of post-settlement juvenile cod with spatial variations in habitat complexity, (2) simulate habitat change based on fishing activities, and (3) determine the role of marine protected areas in enhancing recruitment success. Density-dependent natural mortality was specified as three alternative functional response curves to assess the influence of different predator foraging strategies on juvenile survivorship during the first 12 months of demersal existence. We applied the model to a theoretical patch of hard-bottom substrata and to a case study based on seafloor habitat distributions at Stellwagen Bank National Marine Sanctuary (Gulf of Maine, Northwest Atlantic). Our results demonstrate that patterns in the shape of response surfaces that show the relationship between juvenile cod survivorship and density as well as movement rate were similar regardless of functional response type, that juvenile cod movement rates and post-settlement density were critical for predicting the effects of marine protected-area size on survivorship, and that habitat change caused by fishing has significant negative effects on juvenile cod survivorship and use of marine protected areas can ameliorate such effects.     

Population dynamics of the barnacle <Emphasis Type="Italic">Chthamalus montagui</Emphasis> at two spatial and temporal scales in northern Spain

This study describes spatial and temporal patterns of variability in population parameters in the barnacle Chthamalus montagui Southward in three localities of northern Spain and evaluates whether density-dependent settlement may regulate population dynamics. The sampling design considered two spatial scales, localities and sites within localities, and two temporal scales, years and six month intervals. Density, amount of free space, mortality, growth rate and magnitude of settlement (both absolute and per unit of free space) were obtained from photographs of permanent quadrats and from direct counts in the field. The number of settlers in scraped and untouched quadrats was used to estimate the importance of the presence of conspecifics in settlement. Significant variation at the two spatial and temporal scales was found for most parameters. Large spatial and temporal variations in adult mortality rate, density, and settlement were observed. Patterns of mortality were not consistent with differences in density among localities. Differences in settlement among localities were maintained through time. We suggest that magnitude of settlement is regulated by persistent features such as topography or local water circulation. We assume that early post-settlement mortality does not differ among localities. In the absence of differential mortality, settlement determines average population density within localities. Within localities, settlement was independent of density and free space. No consistent evidence was found on preferential settlement at the vicinity of conspecifics. The main conclusion is that density-dependent settlement is not relevant for the regulation of the populations of C. montagui in the northern Spain. Regulation might occur by density-dependent processes within the adult fraction of the population and/or the larval phase before settlement.     

Density-dependent settlement and mortality structure the earliest life phases of a coral population

The local densities of heterospecifics and conspecifics are known to have profound effects on the dynamics of many benthic species, including rates of settlement and early post-settlement survivorship. We described the early life history of the Caribbean coral, Siderastrea radians by tracking the population dynamics from recently settled planulae to juveniles. Through three years of observation, settlement correlated with the abundance of other benthic organisms, principally turf algae (negatively) and crustose coralline algae (positively). In addition, adult density showed independent effects on coral settlement and early post-settlement survivorship. Settlement rates increased across low levels of adult cover and saturated at a maximum around 10% cover. Early post-settlement survivorship decreased with adult cover, revealing structuring density dependence in coral settlers. The earliest life stages of corals are defined by low survivorship, with survivorship increasing appreciably with colony size. However, recent settlers (one-polyp individuals, < 1-year-old) are more likely to grow into two-polyp juveniles than older single polyps (> 1-year-old) that were delayed in their development. The early benthic phase of corals is defined by a severe demographic bottleneck for S. radians, with appreciable density-dependent and density-independent effects on survivorship. For effective management and restoration of globally imperiled coral reefs, we must focus more attention on this little studied, but dynamic, early life history period of corals.     

Local population size in a flightless insect: importance of patch structure-dependent mortality

In spatially heterogeneous systems, utilizing population models to integrate the effects of multiple population rates can yield powerful insights into the relative importance of the component rates. The relative importance of demographic rates and dispersal in shaping the distribution of the western tussock moth (Orgyia vetusta) among patches of its host plant was explored using stage-structured population models. Tussock moth dispersal occurs passively in first-instar larvae and is poor or absent in all other life stages. Spatial surveys suggested, however, that moth distribution is not well explained by passive dispersal; moth populations were greater on small patches and on isolated ones. Further analysis showed that several local demographic rates varied significantly with patch characteristics. Two mortality factors in particular may explain the observed patterns. First, crawler mortality both increased with patch size and was density-dependent. A single-patch difference equation model showed mortality related to patch size is strong enough to overcome the homogenizing effect of density dependence; greater equilibrium densities were predicted for smaller patches. Second, although three rates were found to vary with local patch density, only pupal parasitism by a chalcid wasp could potentially account for higher moth abundances on isolated patches. A spatially explicit simulation model of the multiple-patch system showed that spatial variation in pupal parasitism is indeed strong enough to generate such a pattern. These results demonstrate that habitat spatial structure can affect multiple population processes simultaneously, and even relatively low attack rates imposed on a reproductively valuable life stage of the host can have a dominant effect on population distribution among habitat patches.     

Information processing in models for density-dependent emigration: A comparison

Density-dependent emigration has been recognized as a fitness enhancing strategy. Yet, especially in the modelling literature there is no consensus about how density-dependent emigration should quantitatively be incorporated into metapopulation models. In this paper we compare the performance of five different dispersal strategies (defined by the functional link between density and emigration probability). Four of these strategies are based on published functional relationships between local population density and emigration probability, one assumes density-independent dispersal. We use individual-based simulations of time-discrete metapopulation dynamics and conduct evolution experiments for a broad range of values for dispersal mortality and environmental stochasticity. For each set of these conditions we analyze the evolution of emigration rates in ‘monoculture experiments’ (with only one type of dispersal strategy used by all individuals in the metapopulation) as well as in selection experiments that allow a pair-wise comparison of the performance of each functional type. We find that a single-parameter ‘asymptotic threshold’ strategy - derived from the marginal value theorem - with a decelerating increase of emigration rate with increasing population density, out-competes any other strategy, i.e. density-independent emigration, a ‘linear threshold’ strategy and a flexible three-parameter strategy. Only when environmental conditions select for extremely high emigration probabilities (close to one), strategies may perform approximately equally. A simple threshold strategy derived for the case of continuous population growth performs even worse than the density-independent strategy. As the functional type of the dispersal function implemented in metapopulation models may severely affect predictions concerning the survival of populations, range expansion, or community changes we clearly recommend to carefully select adequate functions to model density-dependent dispersal.     

Do resources or natural enemies drive bee population dynamics in fragmented habitats?

The relative importance of bottom-up or top-down forces has been mainly studied for herbivores but rarely for pollinators. Habitat fragmentation might change driving forces of population dynamics by reducing the area of resource-providing habitats, disrupting habitat connectivity, and affecting natural enemies more than their host species. We studied spatial and temporal population dynamics of the solitary bee Osmia rufa (Hymenoptera: Megachilidae) in 30 fragmented orchard meadows ranging in size from 0.08 to 5.8 ha in an agricultural landscape in central Germany. From 1998 to 2003, we monitored local bee population size, rate of parasitism, and rate of larval and pupal mortality in reed trap nests as an accessible and standardized nesting resource. Experimentally enhanced nest site availability resulted in a steady increase of mean local population size from 80 to 2740 brood cells between 1998 and 2002. Population size and species richness of natural enemies increased with habitat area, whereas rate of parasitism and mortality only varied among years. Inverse density-dependent parasitism in three study years with highest population size suggests rather destabilizing instead of regulating effects of top-down forces. Accordingly, an analysis of independent time series showed on average a negative impact of population size on population growth rates but provides no support for top-down regulation by natural enemies. We conclude that population dynamics of O. rufa are mainly driven by bottom-up forces, primarily nest site availability.     

The control of emigration and its consequences for the survival of populations

Dispersal is the key process enhancing the long-term persistence of metapopulations in heterogeneous and dynamic landscapes. However, any individual emigrating from a occupied patch also increases the risk of local population extinction. The consequences of this increase for metapopulation persistence likely depend on the control of emigration. In this paper, we present results of individual-based simulations to compare the consequences of density-independent (DIE) and density-dependent (DDE) emigration on the extinction risk of local populations and a two-patch metapopulation. (1) For completely isolated patches extinction risk increases linearly with realised emigration rates in the DIE scenario. (2) For the DDE scenario extinction risk is nearly insensitive to emigration as longs as emigration probabilities remain below ≈0.2. Survival chances are up to half an order of magnitude larger than for populations with DIE. (3) For low dispersal mortality both modes of emigration increase survival of a metapopulation by ca. one order of magnitude. (4) For high dispersal mortality only DDE can improve the global survival chances of the metapopulation. (5) With DDE individuals are only removed from a population at high population density and the risk of extinction due to demographic stochasticity is thus much smaller compared to the DIE scenario.With density-dependent emigration prospects of metapopulations survival may thus be much higher compared to a system with density-independent emigration. Consequently, the knowledge about the factors driving emigration may significantly affect our conclusions concerning the conservation status of species.     

Invasion dynamics of the varnish clam (Nuttallia obscurata): a matrix demographic modeling approach

Marine invaders have become a significant threat to native biodiversity and ecosystem function. In this study, the invasion of the varnish clam (Nuttallia obscurata) in British Columbia, Canada, is investigated using a matrix modeling approach to identify the life history characteristics most crucial for population growth and to investigate population differences. Mark-recapture analyses and field collections from 2003 to 2004 were used to determine individual growth, survival rates, and fecundity for two sites. A multi-state matrix model was used to determine population growth rates and to conduct sensitivity and elasticity analyses. A life table response experiment was also used to determine what life history stage contributed most to observed differences in population growth rates. Population survey data were used in conjunction with the matrix model to determine plausible recruitment levels and to investigate recruitment scenarios. Both populations are currently declining but are likely sustainable because of the pulsed nature of large recruitment events. Survival of larger clams (>40 mm) is the most important for population growth based on elasticity and sensitivity analyses. Adult survival also had the largest influence on observed differences between site-specific population growth rates. The two populations studied differed in recruitment dynamics; one experiencing annual recruitment with higher post-settlement mortality and the other, episodic recruitment and lower post-settlement mortality. The most influential factor for the successful invasion of the varnish clam appears to be survival of the larger size classes. Therefore, any process that decreases adult survival (e.g., predation, commercial harvest) will have the greatest impact on population growth.     

Habitat modification affects recruitment of abalone in central New Zealand

The habitat experienced during early life-history stages can determine the number and quality of individuals that recruit to adult populations. In a field experiment, biogenic habitat complexity was manipulated (presence or absence of foliose macroalgae) at two depths (2–3 m and 5–6 m) and the habitat-dependent effects on recruitment of the black foot abalone (Haliotis iris) were examined at three field sites along the south coast of Wellington, New Zealand (41°20′S, 174°47′E), between July and November 2005. Recruit density (<5 weeks post-settlement) was measured on cobbles covered with crustose coralline algae. Habitats of low complexity (barrens treatments) had consistently greater densities of recruits than habitats of high complexity (algae treatments). However, recruits in algae habitats were larger, and for deep habitats, there was greater survival in algae habitats compared with barrens habitats. While depth had no significant effect on early recruit (<2 weeks post-settlement) density, late recruit (<5 weeks post-settlement) density was greater in shallow habitats, and so it seems recruit survival was greater in shallow habitats. In this experiment, algal habitat complexity had strong effects on early recruit abundance, but habitat-dependent variations in recruit growth and survival may modify initial patterns of abundance and determine recruitment to adult abalone populations.     

Model selection analysis of temporal variation in benefit for an ant-tended treehopper

Recent studies of mutualism have emphasized both that the net benefit to participants depends on the ecological context and that the density-dependent pattern of benefit is key to understanding the population dynamics of mutualism. Indeed, changes in the ecological context are likely to drive changes in both the magnitude of benefit and the density-dependent pattern of benefit. Despite the close linkage between these two areas of research, however, few studies have addressed the factors underlying variation in the density-dependent pattern of benefit. Here I use model selection to evaluate how variation in the benefits of a mutualism drives temporal variation in the density-dependent pattern of net benefit for the ant-tended treehopper Publilia concava. In the interaction between ants and treehoppers in the genus Publilia, ants collect the sugary excretions of treehoppers as a food resource, and treehoppers benefit both directly (e.g., by feeding facilitation) and indirectly (e.g., by predator protection). Results presented here show that temporal changes in the relative magnitude of direct and indirect benefit components of ant tending, especially the effectiveness of predator protection by ants, qualitatively change the overall pattern of density-dependent benefit between years with maximum benefit shifting from treehoppers in small to large aggregations. These results emphasize the need for empirical studies that evaluate the long-term dynamics of mutualism and theoretical studies that consider the population dynamics consequences of variation in the density-dependent pattern of benefit.     

Wide-spread genetic variability and the paradox of effective population size in the gag,Mycteroperca microlepis,along the West Florida Shelf

Wide-ranging marine species are often described as having a low effective population size (N _e) to census size (N) ratio. This genetic phenomenon is typically attributed to large variation among individuals in reproductive success because of the high mortality rates and unpredictable environments associated with larval dispersal. In this study, we examined patterns of genetic variation in gag (Mycteroperca microlepis) on the West Florida Shelf across year classes of post-settlement juveniles and spawning adults. With no significant genetic differentiation among year classes despite varying recruitment dynamics, little evidence for chaotic genetic patchiness, and no truncation of adult genetic diversity in subsequent juvenile cohorts, there was little support for large variation among individual in reproductive success contributing to a low N _e/N ratio. In fact, the consistent lack of significant differences in annual recruitment classes indicated that reproductive success among individuals was resistant to skewing. Among the various evolutionary forces that may be affecting N _e, changes to demography due to fishing pressure are posited as a likely mechanism affecting current levels of genetic variation.     

Dispersal per recruit: an efficient method for assessing sustainability in marine reserve networks.

Marine reserves are an increasingly important tool for the management of marine ecosystems around the world. However, the effects of proposed marine reserve configurations on sustainability and yield of populations are typically not estimated because of the computational intensity of direct simulation and uncertainty in larval dispersal and density-dependent recruitment. Here we develop a method for efficiently assessing a marine reserve configuration for persistence and yield of a population with sedentary adults and dispersing larvae. The method extends the familiar sustainability criteria of individual replacement for single populations based on eggs-per-recruit (EPR) to spatially distributed populations with sedentary adults, a dispersing larval phase, and limited carrying capacity in the settlement-recruit relationship. We refer to this approach as dispersal-per-recruit (DPR). In some cases, a single DPR calculation, based on the assumption that post-settlement habitat is saturated (i.e., at maximum recruitment), is sufficient to determine population persistence, while in other cases further iterative calculations are required. These additional calculations reach an equilibrium more rapidly than a full simulation of age- or size-structured populations. From the DPR result, fishery yield can be computed from yield-per-recruit (YPR) at each point. We assess the utility of DPR calculations by applying them to single reserves, uniformly distributed systems of reserves, and randomly sized and spaced systems of reserves on a linear coastline. We find that for low levels of EPR in fished areas (e.g., 10% or less of the natural, unfished EPR when post-settlement habitats are saturated by 35% of natural settlement), a single DPR calculation is sufficient to determine persistence of the population. We also show that, in uniform systems of reserves with finite reserve size, maximal fisheries yield occurs when the density of reserves is such that all post-settlement habitat is nearly saturated with settlers. Finally, we demonstrate the application of this approach to a realistic proposed marine reserve configuration.     

Spatial variability in habitat associations of pre- and post-settlement stages of coral reef fishes at Ishigaki Island,Japan

Successful settlement of pelagic fish larvae into benthic juvenile habitats may be enhanced by a shortened settlement period, since it limits larval exposure to predation in the new habitat. Because the spatial distribution of marine fish larvae immediately prior to settlement versus during settlement was unknown, field experiments were conducted at Ishigaki Island (Japan) using light trap sampling and underwater visual belt transect surveys to investigate the spatial distribution patterns of selected pre- and post-settlement fishes (Acanthuridae, Pomacentridae, Chaetodonidae and Lethrinidae) among four habitats (seagrass bed, coral rubble, branching coral and tabular coral). The results highlighted two patterns: patterns 1, pre- and post-settlement individuals showing a ubiquitous distribution among the four habitats (Acanthuridae) and pattern 2, pre-settlement individuals distributed in all habitats, but post-settlement individuals restricted to coral (most species of Pomacentridae and Chaetodontidae) or seagrass habitats (Lethrinidae). The first pattern minimizes the transition time between the larval pelagic stage and acquisition of a benthic reef habitat, the latter leading immediately to a juvenile lifestyle. In contrast, the second pattern is characterized by high settlement habitat selectivity by larvae and/or differential mortality immediately after settlement.