Skip to main content

Vertical segregation of holoplanktonic molluscs in the epipelagic layer, southern Gulf of Mexico


The small-scale vertical distribution of five genera of holoplanktonic molluscs inhabiting the upper oceanic layer (0–105 m) was analysed to test whether evident distribution patterns could be attributed to habitat partitioning. Zooplankton samples were taken over neritic waters during the day and night using a 505-μm multiple closing net at five levels (0–6, 6–12, 12–18, 45–55 and 95–105 m) of the water column. Flowmeters were placed in each net to estimate the amount of filtered water. From the 55,654 identifiable specimens of holoplanktonic molluscs, only 6.2 % were in juvenile stage, and Cavolinia, Diacavolinia, Diacria, Oxygyrus and Clio were the most abundant genera. Adults of these genera comprised 3 % of the total abundance. The abundance data were examined using a spatial segregation index, and its significance was tested with null model methods based on Monte Carlo randomizations. Results indicated that adults showed a much more overlapped distribution in the water column than the in the earlier stages. Indeed, juveniles of Cavolinia, Diacavolinia and Oxygyrus were mainly found in the 0–18-m layer, whereas Diacria and Clio were recorded in the 45–105-m stratum. Null model arguments revealed a significant vertical segregation among them. Potential ecological factors involve preferential spawning areas of adults, avoidance of competition for feeding and spatial resources, use of visual capabilities in searching prey, brood protection and avoidance of strong turbulence conditions.


The oceanic plankton distribution is subject to many spatiotemporal scales of variability nested in a hierarchic cascade pattern (Haury et al. 1978; Allen and Starr 1982). From meso- to large scales, the distribution of planktonic organisms is influenced mostly by physical factors, whereas at fine scale, ecological and biological processes become more important (Wiens 1989; Sanvicente-Añorve et al. 2006). In the horizontal plane, planktonic organisms have a limited capability to define their position in space; however, in the vertical axis, many species are well capable of vertical migrations related to their response to light, the availability of food and/or avoidance of predation (Basedow et al. 2008; Cohen and Forward 2009).

The zooplankton distribution throughout the water column displays different patterns of vertical structure for a variety of groups (Criales-Hernández et al. 2008). Field and theoretical studies suggested that this heterogeneous use of space can be mainly attributed to physical processes, seasonal, ontogenetic or diel vertical migrations, body size, morphologies, feeding resources, reproductive success, predation risk and/or competitive abilities (Jeffries and Lawton 1984; Fossheim and Primicerio 2008; Marrari et al. 2011; Sato et al. 2011). The vertical movements of zooplankters exhibit high variability in amplitude, which increases with the developmental stage of the organisms (Dale and Kaartvedt 2000; Cohen and Forward 2009). For instance, in Norwegian waters, the youngest stages of a copepod species were generally confined to the upper 30-m layer, whereas the older specimens were distributed progressively deeper according to their increasing developmental stage (Dale and Kaartvedt 2000).

In the pelagic environment, molluscs represent only a small fraction of the zooplankton biomass. Nevertheless, these animals may have a significant impact on the ecology of epipelagic marine communities and on the global cycle of carbon dioxide. Although many of the planktonic molluscs are rare, they are nonetheless of interest in illustrating the evolution of morphological adaptations to the pelagic environment and unique lifestyles in the ocean (Lalli and Gilmer 1989). Compared to crustaceans or fishes, studies dealing with the vertical distribution or behavioural migration of molluscs are very scarce (Solis and von Westernhagen 1978; Seapy 1990, 2008) and have seldom dealt with their earlier developmental stages. In the Arctic Ocean, Kobayashi (1974) has evidenced several developmental stages of Limacina helicina at different depths of the water column, but the causes remain poorly understood. Based on the assumption that the youngest stages of pelagic molluscs have a limited vertical motility in the water column, we propose that their position in the upper oceanic layer is therefore related to their ecological requirements and/or life history characteristics.

Materials and methods

Field survey and laboratory analyses

The study area comprised neritic waters of the southern Gulf of Mexico, between 18°–22°N and 91°–95°W. Zooplankton sampling was carried out during two periods, May and November 1995, over a sampling grid including 28 oceanographic stations (Fig. 1). Five vertical strata were sampled (0–6, 6–12, 12–18, 45–55 and 95–105 m), and a total of 187 samples were collected during the day and night using a multiple opening–closing net equipped with 75-cm-diameter and 505-μm-mesh-size nets. Day (night) period ranged from 7 am to 6 pm (7 pm to 6 am). A flowmeter was placed on the mouth of each net to estimate the amount of water filtered by nets. Samples were preserved in a 4 % formaldehyde seawater solution neutralized with sodium borate. Salinity and temperature measurements were also taken with a CTD probe at each sampling station. In the laboratory, zooplankton biomass, taken as a measure of food availability, was estimated as displaced volume and standardized to 100 m3 of water (ml 100 m−3). All pelagic molluscs were separated from samples and identified to the genus level. Abundance values of the five genera considered in this study were expressed as ind 100 m−3.

Fig. 1
figure 1

Location of sampling stations in neritic waters of the southern Gulf of Mexico

Data analysis

Mean monthly wind speeds (3.8 m s−1 for May and 6.1 m s−1 for November; Anonymous 1999) were used to estimate the turbulent kinetic energy (k), the dissipation rate of kinetic energy \( \left( \varepsilon \right) \) and the length (η) and velocity (v t ) Kolmogorov scales. Wind speed values were used for the calculations of the significant wave height (H 1/3) and the period (T) of waves, as the JONSWAP Spectrum stated (Hasselmann et al. 1973):

$$ H_{1/3} = 4\left[ {1.67 \times 10^{ - 7} \frac{F}{g}} \right]^{1/2} $$
$$ \omega = \frac{2\pi }{T} = 22\left( {\frac{{g^{2} }}{{U_{10} F}}} \right)^{1/3} $$

where F = fetch length, here taken as 300 km; g = acceleration of gravity; U 10 = wind speed at a height of 10 m; ω = wave frequency.

The values of H 1/3 and T were used to estimate the wave length (L) and the horizontal and vertical velocities (u and w) of the water particles (SPM 1977):

$$ L = \frac{{gT^{2} }}{2\pi }\tanh \left( {\frac{2\pi d}{L}} \right) $$
$$ u = \frac{{H_{{{1 \mathord{\left/ {\vphantom {1 3}} \right. \kern-0pt} 3}}} gT}}{2L}\frac{{\cosh \left( {{\raise0.7ex\hbox{${2\pi \left( {z + d} \right)}$} \!\mathord{\left/ {\vphantom {{2\pi \left( {z + d} \right)} L}}\right.\kern-0pt} \!\lower0.7ex\hbox{$L$}}} \right)}}{{\cosh \left( {{\raise0.7ex\hbox{${2\pi }$} \!\mathord{\left/ {\vphantom {{2\pi } L}}\right.\kern-0pt} \!\lower0.7ex\hbox{$L$}}} \right)}}\cos \left( {\omega t} \right) $$
$$ w = \frac{{H_{{{1 \mathord{\left/ {\vphantom {1 3}} \right. \kern-0pt} 3}}} gT}}{2L}\frac{{\sinh \left( {{\raise0.7ex\hbox{${2\pi \left( {z + d} \right)}$} \!\mathord{\left/ {\vphantom {{2\pi \left( {z + d} \right)} L}}\right.\kern-0pt} \!\lower0.7ex\hbox{$L$}}} \right)}}{{\cosh \left( {{\raise0.7ex\hbox{${2\pi }$} \!\mathord{\left/ {\vphantom {{2\pi } L}}\right.\kern-0pt} \!\lower0.7ex\hbox{$L$}}} \right)}}\sin \left( {\omega t} \right) $$

where H 1/3 = significant wave height (the average height of the largest one-third of the waves); t = time (varying from 0 to T); z = water column depth (negative); and d = bottom depth (positive, here equal to 100 m).

Based on the previous calculations, the fluctuations from the horizontal and vertical mean velocities were computed:

$$ \tilde{u}^{2} = \frac{4}{T}\int\limits_{0}^{T/4} {\left( {\bar{u} - u} \right)^{2} } {\text{d}}t\quad \tilde{w}^{2} = \frac{4}{T}\int\limits_{0}^{T/4}{\left( {\bar{w} - w} \right)^{2} } {\text{d}}t $$

The mean values (\( \bar{u} \) and \( \bar{w} \)) of these velocities were estimated, taken sin(ωt) and cos(ωt) equal to \( \frac{2}{\pi } \), since \( \frac{4}{T}\int_{0}^{T/4} {\sin \omega t\begin{array}{*{20}c} {{\text{d}}t = \frac{4}{T}\int_{0}^{T/4} {\cos \omega t{\text{d}}t = } \frac{2}{\pi }} \\ \end{array} } \) in Eqs. 4 and 5.

The turbulent kinetic energy (k) is a measure of the turbulence intensity. It is defined as the kinetic energy per unit mass (J kg−1) of the horizontal \( \left( {\tilde{u}^{2} } \right) \) and vertical \( \left( {\tilde{w}^{2} } \right) \) velocity variances (Burchard 2002):

$$ k = \frac{1}{2}\left( {\tilde{u}^{2} + \tilde{w}^{2} } \right) $$

Kolmogorov theory states that a turbulent flow is formed by eddies that decrease in diameter with depth. The dissipation rate of the turbulent kinetic energy through the spectrum of eddies in the water column is related to the turbulent velocities \( \left( {\sqrt k } \right) \) and the diameter of eddies \( \left( \ell \right) \), with units of velocity squared per second (m2 s−3) (Tennekes and Lumley 1972):

$$ \varepsilon = \frac{AU}{\ell } $$

where A = constant near to one; U = root-mean-square velocity fluctuations, equal to \( \sqrt k ;\ell \) = diameter of the gyre.

At the superficial waters, ℓ is equal to H 1/3; at deeper waters, ℓ is twice the amplitude \( \left( \zeta \right) \) of the waves, defined as:

$$ \zeta = \frac{{H_{{{1 \mathord{\left/ {\vphantom {1 3}} \right. \kern-0pt} 3}}} }}{2}\frac{{\sinh \left[ {{\raise0.7ex\hbox{${2\pi \left( {z + d} \right)}$} \!\mathord{\left/ {\vphantom {{2\pi \left( {z + d} \right)} L}}\right.\kern-0pt} \!\lower0.7ex\hbox{$L$}}} \right]}}{{\sinh \left( {{\raise0.7ex\hbox{${2\pi d}$} \!\mathord{\left/ {\vphantom {{2\pi d} L}}\right.\kern-0pt} \!\lower0.7ex\hbox{$L$}}} \right)}} $$

The dissipation of energy mainly occurs at the smallest turbulence scales, that is, at the so-called Kolmogorov microscales. At this point, the driving energy is finally overcome by the viscosity \( \left( {\upsilon \approx 1 0^{ - 6} \begin{array}{*{20}c} {{\text{m}}^{ 2} \, {\text{s}}^{ - 1} } \\ \end{array} } \right) \), and the size of the smallest eddy (η) and its characteristic velocity (v t ) are defined as (Tennekes and Lumley 1972) follows:

$$ \eta \equiv \left( {{\raise0.7ex\hbox{${\upsilon^{3} }$} \!\mathord{\left/ {\vphantom {{\upsilon^{3} } \varepsilon }}\right.\kern-0pt} \!\lower0.7ex\hbox{$\varepsilon $}}} \right)^{1/4}$$
$$ v_{t} \equiv \left( {\upsilon \varepsilon } \right)^{1/4}$$

These Kolmogorov microscales are crucial to understanding the interactions of planktonic organisms with their microdistribution patterns in the water column (Rothschild and Osborn 1988).

A spatial segregation index (D) was also calculated between all possible pairs formed by the most abundant genera. This index is (White 1983):

$$ D = \frac{1}{2}\sum\limits_{i = 1}^{n} {\left| {\frac{{N_{1i} }}{{N_{1} }} - \frac{{N_{2i} }}{{N_{2} }}} \right|} $$

where N 1i  = number of individuals of genus 1 at station i; N 2i  = number of individuals of genus 2 at station i; N 1 = total number of individuals of genus 1; N 2 = total number of individuals of genus 2.

The values of D range from 0, in the case of the maximum spatial co-occurrence between the two genera, to 1, which indicates a perfect spatial segregation. The level of significance of an observed index value was estimated through the use of null models (Manly 1991). This procedure consists of the following: (1) estimation of the observed index Do from two vectors, (2) random reallocation of one vector (genera 1 or 2), (3) the repetition of the second step a large number of times (1,000 in this case) to estimate new random values (Dr) in order to find the null distribution of D and (4) the comparison of the Do value within the null distribution: If Do represents a typical value within the D distribution, then the spatial pattern of the two genera seems to be random. In contrast, if Do is on one of the extremes of the distribution, the conclusion is that the two genera are spatially associated (lower tail) or segregated (upper tail). In this study we tested the level of spatial segregation among the genera, and the P value was considered to be the proportion of the Dr values higher than or equal to Do. A MATLAB procedure was developed for these calculations.


Environmental conditions

The hydrological data indicated that surface temperature conditions were warmer in May (~28 °C) than in November (24–26 °C); in deeper strata (45–105 m), the temperature varied from 18.5 to 24 °C in May, and from 21 to 26 °C in November due to a deeper thermocline. Surface salinity values were more variable in November (30–33.5) than in May (35.9–36.8), due to higher continental water discharges in November; in deeper waters (45–105 m), the salinity values were more stable than in surface waters, with values around 36.1 in both seasons.

The turbulence parameter profiles in the southern Gulf of Mexico were different in the two studied periods, especially in the 0–40-m stratum. During May, when calm meteorological conditions prevailed, winds of 3.1 m s−1 induced turbulent kinetic energy (k) values ranging from 1.51 × 10−2 to 4.42 × 10−10 J kg−1, and the associated dissipation rate of energy (ɛ) ranged from 2.26 × 10−3 to 1.14 × 10−10 m2 s−3 (Table 1). Moreover, the size of the smallest eddy was less than 1 mm at surface waters. During November, winds of 6.1 m s−1 induced k values that varied from 3.27 × 10−2 to 1.18 × 10−7 J kg−1 and ɛ values from 4.19 × 10−3 to 3.32 × 10−8 m2 s−3. The values of the Kolmogorov microscales were on the same order of magnitude as those observed in May (Table 1).

Table 1 Values of turbulent kinetic energy (k), dissipation rate of energy (ɛ) and length (η) and velocity (v t ) Kolmogorov scales at different discrete depths under different wind conditions (U 10) in the southern Gulf of Mexico

Vertical distribution of zooplankton biomass and holoplanktonic molluscs

The values of zooplankton biomass (ml 100 m−3) were highest in the upper three levels (0–18 m) during both seasons with no significant differences (t test, P > 0.05) between the comparable strata. In turn, pelagic molluscs’ densities were highest in May, especially in the upper three strata, where significant differences (t test, P < 0.05) were observed between months.

The pelagic mollusc community was mainly composed (92.5 %) of organisms belonging to the Thecosomata order. From the 55,654 identifiable specimens of holoplanktonic molluscs collected in the two seasons, only 6.2 % were in juvenile stage, with Cavolinia, Diacavolinia, Diacria, Oxygyrus and Clio being the most abundant genera; juveniles of Oxygyrus were only collected in May. Adults of these genera represented 3 % of the total abundance.

The patterns of vertical distribution revealed that juveniles of Cavolinia, Diacavolinia and Oxygyrus mainly inhabit the upper layer (0–18 m), while those of Diacria and Clio are more abundant in the deeper stratum (45–105 m) (Table 2), at both noon and midnight (Table 3). The highest spatial segregation values were found between the pairs formed by combinations of genera from the upper (Cavolinia, Diacavolinia and Oxygyrus) and lower strata (Diacria and Clio). Moreover, null models revealed significant P values in the pairs formed by Cavolinia or Diacavolinia with Diacria or Clio (Table 4). The distribution of Diacavolinia and Diacria in the five horizontal planes (Figs. 2, 3) illustrates well the genera inhabiting the upper and lower layers. Adults of the same genera exhibited less marked differences between the upper and deeper strata of the water column (Table 5).

Table 2 Mean abundance (ind 100 m−3 ± SD) of five genera of juvenile holoplanktonic molluscs in five strata (m) of the water column during two sampling periods
Table 3 Day (D) and night (N) mean abundance (ind 100 m−3) of five genera of juvenile holoplanktonic molluscs in five strata (m) of the water column during two sampling periods
Table 4 Values of spatial segregation index (Do) estimated from combinations of five genera of juvenile holoplanktonic molluscs collected in five strata of the water column during May (upper matrix) and November (lower matrix)
Fig. 2
figure 2

Vertical distribution of juveniles of Diacavolinia in the southern Gulf of Mexico during two seasons

Fig. 3
figure 3

Vertical distribution of juveniles of Diacria in the southern Gulf of Mexico during two seasons

Table 5 Values of spatial segregation index (Do) estimated from combinations of five genera of adult holoplanktonic molluscs collected in five strata of the water column during May (upper matrix) and November (lower matrix)


The mean density values of pelagic molluscs in the warmer period, May, almost doubled those recorded in November in the 0–18-m layer. Accordingly, Xu and Li (2005) stated that temperature is a major factor determining the seasonal changes in the molluscs’ total abundance. In the southern Gulf of Mexico, higher temperatures seem to induce reproductive processes of pelagic molluscs, as the total density values showed (Table 1).

The vertical distribution of the holoplanktonic molluscs was closely related to the distribution of the zooplankton biomass (Table 1), suggesting that food availability plays an important role in the position of these organisms in the water column. Pelagic molluscs have a wide variety of feeding mechanisms; they feed upon an ample spectrum of microplankton particles and may also exhibit a certain degree of dietary specialization (Lalli and Gilmer 1989). The thecosomatous pteropods, like most of the specimens encountered here, catch their prey using a mucous web which can be swallowed by such predators in less than 3 min (Gilmer and Harbison 1986). In this study, the lowest observed density of molluscs in the 0–18-m layer occurred at the first level (0–6 m) in November, when the highest turbulence values (k and ɛ) were registered (Table 1). In accordance with this finding, the in situ Arctic explorations performed by Gilmer and Harbison (1991) revealed that the pteropods avoid the ~4-m upper layer because the strong turbulence could disrupt their fragile mucous web. Trophic and other ecological interactions take place in the so-called Kolmogorov microscales. According to Rothschild and Osborn (1988), the ecological processes of plankton are coupled with physical phenomena of the same scales of space and time. For instance, the size of the microparticles (15–800 μm) ingested by the molluscs (Gilmer and Harbison 1991) has the same magnitude as the length Kolmogorov scales here calculated (Table 1). Furthermore, the swimming (~4–14 cm s−1), escape (~10–30 cm s−1) and sinking (~0.5 cm s−1) speeds of the pelagic molluscs (Gilmer and Harbison 1986; Harbison and Gilmer 1992) are within the range of the turbulence velocities (\( \sqrt k \) and v t ) estimated in this study (Table 1).

The vertical distribution of both adults and juveniles of targeted genera was quite different. While juveniles showed a significant vertical segregation, adults exhibited a much more vertical overlapped distribution (Tables 4, 5). In the upper 100-m layer of Philippine waters, Solis and von Westernhagen (1978) reported different patterns in the vertical distribution of adults and juveniles of Cavolinia, Diacria and Clio. Indeed, no adults of Clio were caught in their samples, because juveniles seemed to be distributed shallower in the water column (Furnestin 1961).

The vertical partitioning of Cavolinia/Diacavolinia/Oxygyrus and Diacria/Clio does not match the vertical conditions of temperature and salinity. Also, the differences in the day/night abundance data of Cavolinia and Diacavolinia were not consistent in the two sampling periods (Table 3). Hence, we believe that the effects of hydrological conditions and diel movements of organisms have a negligible role in explaining the observed vertical segregation. Therefore, we propose that ecological interactions and/or life history characteristics of species can explain this result:

Suitable areas for eggs deposition

Taking into account that the earlier stages of planktonic organisms have a limited motility in the water column, we think that the selection of suitable areas for adults’ oviposition is the main cause and the first step for juvenile habitat partitioning. In accordance with Resetarits (1996), the choice of proper sites for eggs deposition affects hatching success, larval performance and parental fitness. This author stressed that the choice of an area for oviposition must optimize locally survivorship of phenotypes. Thus, this life history strategy might avoid competition among organisms having similar ecological requirements, as we show below.

Competition avoidance

In niche theory, the importance of competition for limited feeding or space resources has been widely recognized in the spatial structure of communities. In ecologically similar taxa, such as the cavoliniids here encountered, interspecific competition may account for mutual exclusion. When feeding, the cavoliniids produce a spherical mucous web which can attain 200 mm diameter in Cavolinia, 40 mm in Diacavolinia (as Cavolinia longirostris), 100 mm in Diacria and 40 mm in Clio (Gilmer and Harbison 1986). Then, considering a one-dimensional, 1-m-diameter plankton patch, this means that only five individuals deploying a 200-mm-diameter web can coexist. Unfortunately, the spatial relationships among neighbouring individuals in pelagic communities are practically unknown. However, interspecific competition for food and space seems to be a strong reason for habitat partitioning between Cavolinia/Diavolinia and Diacria/Clio groups.

Visual predation

Within the gastropods, the eyes of heteropod molluscs are unusually large and complex (Land 1982). In Oxygyrus, a small epipelagic carnivorous heteropod, the eyes make smooth scanning movements to detect their prey in the surrounding waters, allowing a wide field of view (Land 1982). Blumer (1999), who examined the photoreceptor organs of different growth stages of a heteropod species, indicated that while the larvae have only two sensory cells, the eyes of both juvenile and adults are well adapted to be visual predators. This fact explains the presence of Oxygyrus juveniles in the upper layer (0–18 m), where light intensity magnifies visual abilities for searching prey.

Brood protection

In the epipelagic layer, visual or tactile predators can be common, while physical refuges are rare (Hamner 1995). A way to protect progeny from predation is providing parental cares. Most pelagic opisthobranch molluscs are protandrous hermaphrodites, laying free-floating mucous egg masses and releasing free-swimming veliger larvae (Lalli and Wells 1973; Lalli and Gilmer 1989), vulnerable to predation. However, some mesopelagic molluscs have developed particular reproductive strategies to guarantee the survival of their offspring. For instance, Clio retains embryos and larvae within the mantle cavity, enhancing the offspring success (Tesch 1946). Brood protection is more effective in the darkest part of the water column studied, where Clio was mainly found (Table 2).

Turbulence conditions

Wind-induced turbulence can strongly impact the vertical distribution of planktonic organisms with either beneficial or detrimental consequences for their populations (Margalef 1997). While turbulence promotes encounters between predator and prey, a relatively high velocity would make catching the prey more difficult (Lewis and Pedley 2001). In the case of Oxygyrus, the only visual predator here studied, turbulence may enhance encounters with its prey. In the cavoliniids Cavolinia and Diacavolinia, which deploy and suspend beneath themselves a fragile mucous feeding web, turbulence may induce more suspended particles to be trapped in the web. There are very few in situ observations of pteropod feeding webs (Gilmer and Harbison 1986; Gilmer 1990; Harbison and Gilmer 1992), and no quantifiable relationships between turbulence intensity and mucous webs resistance have been observed. Empirical studies suggest that turbulence might be a major factor in determining the position of molluscs and other zooplankters in the water column, and sometimes the main responsible of vertical partitioning of co-generic species (Mackas et al. 1993; Tsurumi et al. 2005). It seems that suitable turbulence conditions for Cavolinia and Diacavolinia correspond to the third level (12–18 m), as revealed by their abundance data (Table 2). Even when Diacria and Clio could also obtain a beneficial effect from turbulent levels similar to those in the upper layer, their presence in the lower stratum suggests that other ecological or biological phenomena have more influence on their vertical distribution.


This study provides a major contribution to niche requirements of five pelagic mollusc genera (Cavolinia, Diacavolinia, Diacria, Oxygyrus and Clio) inhabiting the ocean epipelagic layer. An examination of their vertical position in the water column, based on null model arguments, evidenced a small-scale spatial segregation among them. Several ecological mechanisms are involved in the observed spatial segregation, each of them having differential importance upon the distribution of each genus. Habitat partitioning of ecologically related genera, such as cavoliniids, suggests an avoidance competition. Thus, position in the water column of juvenile pelagic molluscs represents an equilibrium between the need to minimize the detrimental consequences of competition, the risk of predation and/or the effects of strong turbulence against the necessity of inhabiting suitable areas for feeding and growth.


  • Allen THF, Starr B (1982) Hierarchy: perspectives for ecological complexity. University of Chicago Press, Chicago

    Google Scholar 

  • Anonymous (1999) Datos Climáticos en Veracruz, Ver., México 1917–1998. Comisión Nacional del Agua. Gerencia Estatal en Veracruz. Centro de Previsión del Golfo de México, México DF

  • Basedow SL, Edvardsen A, Tande KS (2008) Vertical segregation of Calanus finmarchicus copepodites during the spring bloom. J Mar Syst 70:21–32

    Article  Google Scholar 

  • Blumer MJF (1999) Development of a unique eye: photoreceptors of the pelagic predator Atlanta peroni (Gastropoda, Heteropoda). Zoomorphology 119:81–91

    Article  Google Scholar 

  • Burchard H (2002) Applied turbulence modelling in marine waters. Lecture notes in earth sciences 100. Springer, Berlin

  • Cohen JH, Forward RB Jr (2009) Zooplankton diel vertical migration—a review of proximate control. Oceanogr Mar Biol Annu Rev 47:77–110

    Article  Google Scholar 

  • Criales-Hernández MI, Schwamborn R, Graco M, Ayón P, Hirche HJ, Wolff M (2008) Zooplankton vertical distribution and migration off central Peru in relation to the minimum layer. Helgol Mar Res 62(Sup 1):S85–S100

    Google Scholar 

  • Dale T, Kaartvedt S (2000) Diel patterns in stage-specific vertical migration of Calanus finmarchicus in habitats with midnight sun. ICES J Mar Sci 57:1800–1818

    Article  Google Scholar 

  • Fossheim M, Primicerio R (2008) Habitat choice by marine zooplankton in a high-latitude ecosystem. Mar Ecol Prog Ser 364:47–56

    Article  Google Scholar 

  • Furnestin ML (1961) Pteropodes et heteropodes du plancton marocain. Rev Trav Inst Pêches marit 25:293–326

    Google Scholar 

  • Gilmer RW (1990) In situ observations of feeding behavior of thecosome pteropod molluscs. Amer Malac Bull 8:53–59

    Google Scholar 

  • Gilmer RW, Harbison GR (1986) Morphology and field behavior of pteropod molluscs: feeding methods in the families Cavoliniidae, Limacinidae and Peraclidae (Gastropoda: Thecosomata). Mar Biol 91:47–57

    Article  Google Scholar 

  • Gilmer RW, Harbison GR (1991) Diet of Limacina helicina (Gastropoda: Thecosomata) in Arctic waters in midsummer. Mar Ecol Prog Ser 77:125–134

    Article  Google Scholar 

  • Hamner WM (1995) Predation, cover, and convergent evolution in pelagic oceans. Mar Freshw Behav Physiol 26:71–89

    Article  Google Scholar 

  • Harbison GR, Gilmer RW (1992) Swimming, buoyancy and feeding in shelled pteropods: a comparison of filed and laboratory observations. J Moll Stud 58:337–339

    Article  Google Scholar 

  • Hasselmann K, Barnett TP, Bouws E, Carlson H, Cartwright DE, Enke K, Ewing JA et al (1973) Measurements of wind-wave growth and swell decay during the Joint North Sea Wave Project (JONSWAP). Dtsch Hydrogr Z Suppl A8:1–95

    Google Scholar 

  • Haury LR, McGowan JS, Wiebe PH (1978) Patterns and processes in the time and space scales for plankton distribution. In: Steel JH (ed) Spatial patterns in plankton communities. Plenum Press, London, pp 237–327

    Google Scholar 

  • Jeffries MJ, Lawton JH (1984) Enemy free space and the structure of ecological communities. Biol J Linn Soc 23:269–286

    Article  Google Scholar 

  • Kobayashi HA (1974) Growth cycle and related vertical distribution of the Thecosomatous pteropod Spiratella (“Limacina”) helicina in the central Arctic Ocean. Mar Biol 26:295–301

    Article  Google Scholar 

  • Lalli CM, Gilmer RW (1989) Pelagic snails. The biology of holoplanktonic gastropod mollusks. Stanford University Press, Stanford

    Google Scholar 

  • Lalli CM, Wells FE (1973) Brood protection in an epipelagic thecosomatous pteropod, Spiratella (“Limacina”) inflata (D’Orbigny). Bull Mar Sci 23:933–941

    Google Scholar 

  • Land MF (1982) Scanning eye movements in a heteropod mollusc. J Exp Biol 96:427–430

    Google Scholar 

  • Lewis DM, Pedley TJ (2001) The influence of turbulence on plankton predation strategies. J Theor Biol 210:347–356

    Article  PubMed  CAS  Google Scholar 

  • Mackas DL, Sefton H, Miller CB, Raich A (1993) Vertical habitat partitioning by large calanoid copepod in the oceanic subarctic Pacific during spring. Prog Oceanogr 32:259–294

    Article  Google Scholar 

  • Manly BFJ (1991) Randomization and Monte Carlo methods in biology. Chapman and Hall, London

    Google Scholar 

  • Margalef R (1997) Turbulence and marine life. Sci Mar 61:109–123

    Google Scholar 

  • Marrari M, Daly KL, Timonin A, Semenova T (2011) The zooplankton of Marguerite Bay, western Atlantic Peninsula—part II: vertical distributions and habitat partitioning. Deep-Sea Res 58:1599–1613

    Article  CAS  Google Scholar 

  • Resetarits WJ Jr (1996) Oviposition site choice and life history evolution. Amer Zool 36:205–215

    Google Scholar 

  • Rothschild BJ, Osborn TR (1988) Small-scale turbulence and plankton contact rates. J Plankton Res 10:465–474

    Article  Google Scholar 

  • Sanvicente-Añorve L, Soto LA, Espinosa-Fuentes ML, Flores-Coto C (2006) Relationship patterns between ichthyoplankton and zooplankton: a conceptual model. Hydrobiologia 559:11–22. doi:10.1007/s10750-005-1323-x

    Google Scholar 

  • Sato K, Yamaguchi A, Ueno H, Ikeda T (2011) Vertical segregation within four grazing copepods in the Oyashio region during early spring. J Plankton Res 33:1230–1238. doi:10.1093/plankt/fbr018

    Google Scholar 

  • Seapy RR (1990) Patterns of vertical distribution in epipelagic heteropod mollusc off Hawaii. Mar Ecol Prog Ser 60:235–246

    Article  Google Scholar 

  • Seapy RR (2008) Offshore-inshore and vertical distributional patterns of heteropod mollusks off leeward Oahu, Hawaii. Mar Biol 15:234–246

    Google Scholar 

  • Solis NB, von Westernhagen H (1978) Vertical distribution of euthecosomatous pteropods in the upper 100 m of the Hilutangan Channel, Cebu, The Philippines. Mar Biol 48:78–87

    Article  Google Scholar 

  • SPM (1977) Shore protection manual, vol I. US Army Coastal Engineering Research Center, Washington

    Google Scholar 

  • Tennekes H, Lumley JL (1972) A first course in turbulence. The MIT Press, Cambridge

    Google Scholar 

  • Tesch JJ (1946) The thecosomatous pteropods. Dana Report No. 28

  • Tsurumi M, Mackas DL, Whitney FA, DiBacco C, Galbraith MD, Wong CS (2005) Pteropods, eddies, carbon flux, and climate variability in the Alaska Gyre. Deep-Sea Res 52:1037–1053

    Article  CAS  Google Scholar 

  • White MJ (1983) The measurement of spatial segregation. Am J Sociol 88:1008–1018

    Article  Google Scholar 

  • Wiens JA (1989) Spatial scaling in ecology. Funct Ecol 3:385–397

    Article  Google Scholar 

  • Xu ZL, Li C (2005) Horizontal distribution and dominant species of heteropods in the east China Sea. J Plankton Res 27:373–382

    Article  Google Scholar 

Download references


We thank F. Zavala-García and M. Hermoso-Salazar for their invaluable assistance in this research and the anonymous reviewers for their comments and suggestions to improve our manuscript. We are also grateful to C. Ramírez and O. Barrientos for their help in searching literature. This study was partially supported by the Programa de Apoyo a Proyectos de Investigación e Innovación Tecnológica, Universidad Nacional Autónoma de México (IN229010 Project).

Author information

Authors and Affiliations


Corresponding author

Correspondence to L. Sanvicente-Añorve.

Additional information

Communicated by H.-D. Franke.

Rights and permissions

Reprints and permissions

About this article

Cite this article

Sanvicente-Añorve, L., Lemus-Santana, E., Flores-Coto, C. et al. Vertical segregation of holoplanktonic molluscs in the epipelagic layer, southern Gulf of Mexico. Helgol Mar Res 67, 397–405 (2013).

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: