Density dependence in marine fish populations revealed at small and large spatial scales

Experimental manipulation of population density has frequently been used to demonstrate demographic density dependence. However, such studies are usually small scale and typically provide evidence of spatial (within-generation) density dependence. It is often unclear whether small-scale, experimental tests of spatial density dependence will accurately describe temporal (between-generation) density dependence required for population regulation. Understanding the mechanisms generating density dependence may provide a link between spatial experiments and temporal regulation of populations. In this study, I manipulated the density of recently settled kelp rockfish (Sebastes atrovirens) in both the presence and absence of predators to test for density-dependent mortality and whether predation was the mechanism responsible. I also examined mortality of rockfish cohorts within kelp beds throughout central California to evaluate temporal (between-generation) density dependence in mortality. Experiments suggested that short-term behavioral responses of predators and/or a shortage of prey refuges caused spatial density dependence. Temporal density dependence in mortality was also detected at larger spatial scales for several species of rockfish. It is likely that short-term responses of predators generated both spatial and temporal density dependence in mortality. Spatial experiments that describe the causal mechanisms generating density dependence may therefore be valuable in describing temporal density dependence and population regulation.     

Cleaner-client interactions on a Caribbean reef: influence of correlates of parasitism

We investigated the influence of known correlates of parasitism, namely fish density, body size and social behaviour, on three highly variable aspects of the interactions between cleaning gobies (Elacatinus spp.) and their clients, on a Barbadian coral reef. We specifically considered (1) variability in client visit rate to cleaning stations, (2) cleaning goby preference for specific clients and (3) variation in the time spent cleaning by cleaners. Using phylogenetically independent contrasts, we found that client species that were abundant on the reef visited cleaning stations more often than rarer client species. This could be due to the positive relationship between ectoparasite load and client density, or alternatively may simply reflect the frequency of contact between cleaners and clients. Cleaning gobies spent more time cleaning large-bodied clients, which usually have higher ectoparasite loads, although cleaning goby preference for clients was influenced by none of the correlates of client parasitism. Overall, factors assumed to correlate with ectoparasite load had a limited influence on the variability observed in the interactions between cleaning gobies and their clients. Received: 27 October 1999 / Received in revised form: 11 January 2000 / Accepted: 24 January 2000     

Widespread density-dependent seedling mortality promotes species coexistence in a highly diverse Amazonian rain forest

Negative density-dependent mortality can promote species coexistence through a spacing mechanism that prevents species from becoming too locally abundant. Negative density-dependent seedling mortality can be caused by interactions among seedlings or between seedlings and neighboring adults if the density of neighbors affects the strength of competition or facilitates the attack of natural enemies. We investigated the effects of seedling and adult neighborhoods on the survival of newly recruited seedlings for multiple cohorts of known age from 163 species in Yasuni National Park, Ecuador, an ever-wet, hyper-diverse lowland Amazonian rain forest. At local scales, we found a strong negative impact on first-year survival of conspecific seedling densities and adult abundance in multiple neighborhood sizes and a beneficial effect of a local tree neighborhood that is distantly related to the focal seedling. Once seedlings have survived their first year, they also benefit from a more phylogenetically dispersed seedling neighborhood. Across species, we did not find evidence that rare species have an advantage relative to more common species, or a community compensatory trend. These results suggest that the local biotic neighborhood is a strong influence on early seedling survival for species that range widely in their abundance and life history. These patterns in seedling survival demonstrate the role of density-dependent seedling dynamics in promoting and maintaining diversity in understory seedling assemblages. The assemblage-wide impacts of species abundance distributions may multiply with repeated cycles of recruitment and density-dependent seedling mortality and impact forest diversity or the abundance of individual species over longer time scales.     

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.     

Neighborhood and community interactions determine the spatial pattern of tropical tree seedling survival

Factors affecting survival and recruitment of 3531 individually mapped seedlings of Myristicaceae were examined over three years in a highly diverse neotropical rain forest, at spatial scales of 1-9 m and 25 ha. We found convincing evidence of a community compensatory trend (CCT) in seedling survival (i.e., more abundant species had higher seedling mortality at the 25-ha scale), which suggests that density-dependent mortality may contribute to the spatial dynamics of seedling recruitment. Unlike previous studies, we demonstrate that the CCT was not caused by differences in microhabitat preferences or life history strategy among the study species. In local neighborhood analyses, the spatial autocorrelation of seedling survival was important at small spatial scales (1-5 m) but decayed rapidly with increasing distance. Relative seedling height had the greatest effect on seedling survival. Conspecific seedling density had a more negative effect on survival than heterospecific seedling density and was stronger and extended farther in rare species than in common species. Taken together, the CCT and neighborhood analyses suggest that seedling mortality is coupled more strongly to the landscape-scale abundance of conspecific large trees in common species and the local density of conspecific seedlings in rare species. We conclude that negative density dependence could promote species coexistence in this rain forest community but that the scale dependence of interactions differs between rare and common species.     

The scale and cause of spatial heterogeneity in strength of temporal density dependence

The importance of density dependence in natural communities continues to spark much debate because it is fundamental to population regulation. We used temporal manipulations of density to explore potentially stabilizing density dependence in early survivorship among six local populations of a tropical damselfish (Dascyllus flavicaudus). Specifically, we tested the premise that spatial heterogeneity in the strength of temporal density dependence would reflect variation in density of predators, the agent of mortality. Our field manipulations revealed that mortality among successive cohorts of young fishes was density dependent at each reef, but that its strength varied by approximately 1.5 orders of magnitude. This spatial heterogeneity was well predicted by variation among the six reefs in the density of predatory fishes that consume juvenile damselfishes. Because density dependence arose from competition for enemy-free space within a shelter coral, the mortality consequence of the competition depended on the neighborhood density of predators. Thus, the scale of heterogeneity in the density dependence largely reflected attributes of the environment that shaped the local abundance of predators. These results have important implications for how ecologists explore regulatory processes in nature. Failure to account for spatial variation could frequently yield misleading conclusions regarding density dependence as a stabilizing process, obscure underlying mechanisms influencing its strength, and provide no insight into the spatial scale of the heterogeneity. Further, models of population dynamics will be improved when experimental approaches better estimate the magnitude and causes of variation in strength of stabilizing density dependence.     

Strength of evidence for density dependence in abundance time series of 1198 species

Population limitation is a fundamental tenet of ecology, but the relative roles of exogenous and endogenous mechanisms remain unquantified for most species. Here we used multi-model inference (MMI), a form of model averaging, based on information theory (Akaike's Information Criterion) to evaluate the relative strength of evidence for density-dependent and density-independent population dynamical models in long-term abundance time series of 1198 species. We also compared the MMI results to more classic methods for detecting density dependence: Neyman-Pearson hypothesis-testing and best-model selection using the Bayesian Information Criterion or cross-validation. Using MMI on our large database, we show that density dependence is a pervasive feature of population dynamics (median MMI support for density dependence = 74.7-92.2%), and that this holds across widely different taxa. The weight of evidence for density dependence varied among species but increased consistently, with the number of generations monitored. Best-model selection methods yielded similar results to MMI (a density-dependent model was favored in 66.2-93.9% of species time series), while the hypothesis-testing methods detected density dependence less frequently (32.6-49.8%). There were no obvious differences in the prevalence of density dependence across major taxonomic groups under any of the statistical methods used. These results underscore the value of using multiple modes of analysis to quantify the relative empirical support for a set of working hypotheses that encompass a range of realistic population dynamical behaviors.     

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.     

Safety in numbers and the spatial scaling of density-dependent mortality in a coral reef fish

In coral reef fishes, density-dependent population regulation is commonly mediated via predation on juveniles that have recently settled from the plankton. All else being equal, strong density-dependent mortality should select against the formation of high-density aggregations, yet the juveniles of many reef fishes aggregate. In light of this apparent contradiction, we hypothesized that the form and intensity of density dependence vary with the spatial scale of measurement. Individual groups might enjoy safety in numbers, but predators could still produce density-dependent mortality at larger spatial scales. We investigated this possibility using recently settled juvenile bluehead wrasse, Thalassoma bifasciatum, a small, aggregating reef fish. An initial caging experiment demonstrated that juvenile bluehead wrasse settlers suffer high predation, and spatial settlement patterns indicated that bluehead wrasse juveniles preferentially settle in groups, although they are also found singly. We then monitored the mortality of recently settled juveniles at two spatial scales: microsites, occupied by individual fish or groups of fish and separated by centimeters, and sites, consisting of approximately 2400-m2 areas of reef and separated by kilometers. At the microsite scale, we measured group size and effective population density independently and found that per capita mortality decreased with group size but was not related to density. At the larger spatial scale, however, per capita mortality increased with settler density. This shift in the form of density dependence with spatial scale could reconcile the existence of small-scale aggregative behavior typical of many reef fishes with the population-scale density dependence that is essential to population stability and persistence.     

Tracking bioaccumulation in aquatic organisms: A dynamic model integrating life history characteristics and environmental change

We have developed a dynamic model to track the evolution of contaminant concentration in an aquatic organism as a function of season and ontogeny throughout its life cycle. We have focused our analysis on the round goby (Apollonia melanostoma), a globally distributed invasive forage fish. By integrating bioenergetics with a bioaccumulation model, we illustrate how life history characteristics interact to influence contaminant accumulation. We use uncertainty and sensitivity analyses to assess how the model output is affected by uncertainty and variability in model parameters. We then demonstrate the influence of important physiological characteristics on contaminant accumulation with two scenarios. First, we probe the influence of sexual dimorphism by comparing gender-specific accumulation of a standard polychlorinated biphenyl congener, PCB153, in male and female round gobies. We hypothesize that lipid loss in female gobies during spawning season leads to a decrease in the PCB body burden compared to male gobies. Second, we compare PCB accumulation in the round goby and in the mottled sculpin (Cottus bairdi), the native forage fish that the round goby displaced in southern Lake Michigan, to determine whether the invasive species has an intrinsically different bioaccumulation potential than the native one. Our non-intuitive findings from these simulations illustrate how the interaction of growth rate with other life history characteristics lead to unexpected bioaccumulation patterns. The model we present is a flexible tool that integrates complex and dynamic interactions among environmental parameters, thus providing a means to better assess the potential for chemical accumulation in human and wildlife populations, and aiding the development of ecological forecasts.     

Using an individual-based model to quantify scale transition in demographic rate functions: Deaths in a coral reef fish

Scientifically informed population management requires quantitatively accurate demographic rate functions that apply at the spatial scale at which populations are actually managed, but practical constraints confine most field measurements of such functions to small study plots. This paper employs an individual-based population growth model to extrapolate the death rate function in a well-studied coral reef fish, the bridled goby Coryphopterus glaucofraenum, from the scale of coral reef “cells” at which it was measured to the larger scale of an entire coral reef. Density dependence in the whole-reef function actually proves stronger than in the local function because high goby density occasionally arises in local patches with few refuges from predators, producing very high mortality there. This IBM-based approach extends the reach of scale transition theory by examining considerably more realistic models than standard analytical methods can presently handle.     

Seasonal density dependence, timing of mortality, and sustainable harvesting

Birth-pulse populations are often characterized with discrete-time models, that use a single function to relate the post-breeding population size to the post-breeding size of the previous year. Recently, models of seasonal density dependence have been constructed that emphasize interactions during shorter time periods also. Here, we study two very simple forms of density-dependent mortality, that lead to Ricker and Beverton-Holt type population dynamics when viewed over the whole year. We explore the consequences of harvest timing to equilibrium population sizes under such density dependence. Whether or not individual mortality compensates for the harvested quota, the timing of harvesting has a strong impact on the sustainability of a harvesting quota. Further, we show that careless discretization of a continuous mortality scheme may seriously underestimate the reduction in population size caused by hunting and overestimate the sustainable yield. We also introduce the concept of the demographic value of an individual, which reflects the expected contribution to population size over time in the presence of density dependence. Finally, we discuss the possibility of calculating demographic values as means of optimizing harvest strategies. Here, a Pareto optimal harvest strategy will minimize the loss of demographic value from the population for a given yield.     

Behavioral responses of territorial red squirrels to natural and experimental variation in population density

The relative scarcity of studies at the intersection of behavioral and population ecology is surprising given the presumed importance of behavior in density-dependent population regulation. Here we tested whether North American red squirrels (Tamiasciurus hudsonicus) adjust their behavior in response to local population density and whether they use rates of territorial vocalizations in their local neighborhood to assess density. We examined these relationships using 18?years of live trapping and 20?years of behavioral data that were collected across natural variation in local population density. To disentangle the effects of population density on behavior from those due to changes in per capita food abundance or changes in the frequency of antagonistic interactions with neighbors, we also experimentally manipulated population density with long-term food supplementation as well as perceived population density with long-term playbacks of territorial vocalizations. The frequency with which squirrels emitted territorial vocalizations was positively associated with local population density. In contrast, antagonistic physical interactions observed between squirrels and territorial intrusions were rare and the frequency of intrusions was weakly and negatively, not positively, associated with population density. Squirrels experiencing naturally and experimentally high density conditions spent less time in the nest and feeding but more time being vigilant. Similar density-dependent changes in behavior were observed in response to our manipulations of perceived population density, indicating that vocalization rates and not physical interactions or food abundance were the mechanism by which squirrels assessed and responded behaviorally to changes in local density.     

Density dependence and density independence during the early life stages of four marine fish stocks

Recruitment variability caused by density-dependent and density-independent processes is an important area within the study of fish dynamics. These processes can exhibit nonlinearities and nonadditive properties that may have profound dynamic effects. We investigate the importance of population density (i.e., density dependence) and environmental forcing (i.e., density independence) on the age-0 and age-1 abundance of capelin (Mallotus villosus), northeast Arctic cod (Gadus morhua), northeast Arctic haddock (Melanogrammus aeglefinus), and Norwegian spring spawning herring (Clupea harengus) in the Barents Sea. We use statistical methods that explicitly account for nonlinearities and nonadditive interactions between internal and external variables in the abundance of these two pre-recruitment stages. Our results indicate that, during their first five months of life, cod, haddock, and herring experience higher density-dependent survival than capelin. The abundance of age-0 cod depends on the mean age and biomass of the spawning stock, a result which has implications for the management of the entire cod stock. Temperature is another important factor influencing the abundance at age-0 and age-1 of all four species, except herring at age-1. Between age-0 and age-1, there is an attenuation of density-dependent survival for cod and herring, while haddock and capelin experience density dependence at high and low temperatures, respectively. Predation by subadult cod is important for both capelin and cod at age-1. We found strong indications for interactions among the studied species, pointing to the importance of viewing the problem of species recruitment variability as a community, rather than as a population phenomenon.     

Refuge sharing network predicts ectoparasite load in a lizard

Living in social groups facilitates cross-infection by parasites. However, empirical studies on indirect transmission within wildlife populations are scarce. We investigated whether asynchronous overnight refuge sharing among neighboring sleepy lizards, Tiliqua rugosa, facilitates indirect transmission of its ectoparasitic tick, Amblyomma limbatum. We fitted 18 neighboring lizards with GPS recorders, observed their overnight refuge use each night over 3 months, and counted their ticks every fortnight. We constructed a transmission network to estimate the cross-infection risk based on asynchronous refuge sharing frequencies among all lizards and the life history traits of the tick. Although self-infection was possible, the network provided a powerful predictor of measured tick loads. Highly connected lizards that frequently used their neighbors’ refuges were characterized by higher tick loads. Thus, indirect contact had a major influence on transmission pathways and parasite loads. Furthermore, lizards that used many different refuges had lower cross- and self-infection risks and lower tick loads than individuals that used relatively fewer refuges. Increasing the number of refuges used by a lizard may be an important defense mechanism against ectoparasite transmission in this species. Our study provides important empirical data to further understand how indirectly transmitted parasites move through host populations and influence individual parasite loads.     

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.     

Experimental test of parasitism hypothesis for population cycles of a forest lepidopteran

Population cycles of herbivores are thought to be driven by trophic interaction mechanisms, either between food plant and herbivore or between the herbivorous prey and its natural enemies. Observational data have indicated that hymenopteran parasitoids cause delayed density-dependent mortality in cyclic autumnal moth (Epirrita autumnata) populations. We experimentally tested the parasitism hypothesis of moth population cycles by establishing a four-year parasitoid-exclusion experiment, with parasitoid-proof exclosures, parasitoid-permeable exclosures, and control plots. The exclusion of parasitoids led to high autumnal moth abundances, while the declining abundance in both the parasitoid-permeable exclosures and the control plots paralleled the naturally declining density in the study area and could be explained by high rates of parasitism. Our results provide firm experimental support for the hypothesis that hymenopteran parasitoids have a causal relationship with the delayed density-dependent component required in the generation of autumnal moth population cycles.     

Spatial and temporal variability modify density dependence in populations of large herbivores

A central challenge in ecology is to understand the interplay of internal and external controls on the growth of populations. We examined the effects of temporal variation in weather and spatial variation in vegetation on the strength of density dependence in populations of large herbivores. We fit three subsets of the model ln(Nt) = a + (1 + b) x ln(N(t-1)) + c x ln(N(t-2)) to five time series of estimates (Nt) of abundance of ungulates in the Rocky Mountains, USA. The strength of density dependence was estimated by the magnitude of the coefficient b. We regressed the estimates of b on indices of temporal heterogeneity in weather and spatial heterogeneity in resources. The 95% posterior intervals of the slopes of these regressions showed that temporal heterogeneity strengthened density-dependent feedbacks to population growth, whereas spatial heterogeneity weakened them. This finding offers the first empirical evidence that density dependence responds in different ways to spatial heterogeneity and temporal heterogeneity.     

Nonlinear effects of climate and density in the dynamics of a fluctuating population of reindeer

Nonlinear and irregular population dynamics may arise as a result of phase dependence and coexistence of multiple attractors. Here we explore effects of climate and density in the dynamics of a highly fluctuating population of wild reindeer (Rangifer tarandus platyrhynchus) on Svalbard observed over a period of 29 years. Time series analyses revealed that density dependence and the effects of local climate (measured as the degree of ablation [melting] of snow during winter) on numbers were both highly nonlinear: direct negative density dependence was found when the population was growing (Rt > 0) and during phases of the North Atlantic Oscillation (NAO) characterized by winters with generally high (1979-1995) and low (1996-2007) indices, respectively. A growth-phase-dependent model explained the dynamics of the population best and revealed the influence of density-independent processes on numbers that a linear autoregressive model missed altogether. In particular, the abundance of reindeer was enhanced by ablation during phases of growth (Rt > 0), an observation that contrasts with the view that periods of mild weather in winter are normally deleterious for reindeer owing to icing of the snowpack. Analyses of vital rates corroborated the nonlinearity described in the population time series and showed that both starvation mortality in winter and fecundity were nonlinearly related to fluctuations in density and the level of ablation. The erratic pattern of growth of the population of reindeer in Adventdalen seems, therefore, to result from a combination of the effects of nonlinear density dependence, strong density-dependent mortality, and variable density independence related to ablation in winter.     

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.