A correction has been applied to this article in:
Corrigendum: Modelling the Distribution of a Commercial NE-Atlantic Sea Cucumber, Holothuria mammata: Demographic and Abundance Spatio-Temporal Patterns
Pedro M. Félix1*
Ana Pombo2
Francisco Azevedo e Silva1
Tomás Simões1
Tiago A. Marques3,4,5
Ricardo Melo1,6Cristina Rocha3,7João Sousa2
Eliana Venâncio2
José L. Costa1,5
Ana C. Brito1,6- 1MARE – Marine and Environmental Sciences Centre, Faculdade de Ciências, Universidade de Lisboa, Lisbon, Portugal
- 2MARE – Marine and Environmental Sciences Centre, ESTM, Polytechnic Institute of Leiria, Leiria, Portugal
- 3CEAUL – Centro de Estatística e Aplicações, Faculdade de Ciências, Universidade de Lisboa, Lisbon, Portugal
- 4Centre for Research into Ecological and Environmental Modelling, The Observatory, University of St Andrews, St Andrews, United Kingdom
- 5Departamento de Biologia Animal, Faculdade de Ciências, Universidade de Lisboa, Lisbon, Portugal
- 6Departamento de Biologia Vegetal, Faculdade de Ciências, Universidade de Lisboa, Lisbon, Portugal
- 7Departamento de Estatística e Investigação Operacional, Faculdade de Ciências, Universidade de Lisboa, Lisbon, Portugal
There is an increasing demand for sea cucumbers, for human consumption, mainly from Asian markets and, as a consequence, NE-Atlantic species are now new targets for exploitation and exportation. Holothuria mammata is one of the most valuable species in Europe. However, the lack of historical economic interest in this species in most European countries has also led to a lack of studies concerning biological and ecological aspects on wild populations and this is a major issue for stock management. This study aims to determine the temporal and spatial patterns of distribution of H. mammata, considering its abundance and demographic structure in a NE-Atlantic area, SW Portugal, as a function of environmental conditions. For that, a population from a marine protected area was followed for 1 year at 1.5-month intervals. Throughout the coastal area, six sites were selected and at each sampling campaign three random transects per site and substrate (rock and sand) in which all H. mammata individuals were counted and measured. For each site and survey several environmental parameters of interest, from the water column, the sediment and substrate cover, were also measured. Generalized Linear Models were used to model the spatial and temporal distribution of the species according to environmental conditions, to determine the species’ habitat preferences. The distribution models indicate that abiotic and biotic parameters of the water column are not the main drivers shaping the distribution of H. mammata. The species has a patchy distribution, and its habitat preferences depend on environmental stability, the presence of shelter and habitat complexity, which is more important for smaller, more vulnerable, individuals, while bigger size classes tend to venture more into less stable environments in an opportunistic fashion. The knowledge of these population traits is determinant to develop stock management measures, which are now urgent to prevent the depletion of commercial sea cucumber populations in the NE-Atlantic. Sustainable fisheries policies should be developed and start by considering to delimit fishing areas and periods, considering the species spatial and temporal distribution patterns.
Introduction
Sea cucumbers are a group of echinoderms with important ecological functions. Deposit-feeders compose a large group of these holothurians, which play a critical role in habitat structuring through bioturbation, recycling and redistributing nutrients (e.g., MacTavish et al., 2012; Purcell et al., 2016). Their feeding behaviour, ingesting sediments and digesting its organic content, help to maintain or improve the physico-chemical processes of benthic habitats, by decomposing sediment organic matter and controlling populations of bacteria, fungi, and phytobenthos (e.g., Uthicke, 1999; Michio et al., 2003; Namukose et al., 2016). Sea cucumbers are a highly priced commodity in China and have been harvested for more than three centuries (Schwerdtner Máñez and Ferse, 2010). As a consequence of high fishing pressure and intrinsic biological and ecological characteristics, such as low mobility, density dependent breeding, slow growth of many species, or shallow water occurrence for others (Hamel and Mercier, 1996a; Kinch et al., 2008; Mercier and Hamel, 2009; Purcell et al., 2013), this group is vulnerable to overfishing. The majority of exploited sea cucumber species are, in fact, deposit-feeders (most of the formerly known aspidochirotids) (Purcell et al., 2016), and several of these are currently overexploited or depleted (Purcell et al., 2013). This can create a set of inevitable ecological consequences that can result in a curtailment in primary productivity, decrease in benthic diversity and biomass, decrease of buffering capacity of surrounding seawater acidity, preclude symbiotic associations, impair trophic relationships, etc. (Moriarty et al., 1985; Uthicke, 1999; Michio et al., 2003; Solan, 2004; So et al., 2010; Wolkenhauer et al., 2010; Schneider et al., 2011, 2013; Purcell et al., 2016). Considering the ecosystem services provided by these species, the overexploitation of multiple populations has the potential to create severe cascading effects in the ecosystems.
Setting limits to fisheries is now urgent but, for that to be done efficiently, improved knowledge on the ecological traits of exploited species is of the utmost importance. Distribution patterns and their drivers are a fundamental set of information for stock management, e.g., establishing geographical limits to harvest and avoid excessive capture of breeders or specific size-classes, particularly due to a common patchy distribution of sea cucumbers, where high density areas can be followed by adjacent zero or near-zero abundances (e.g., Džeroski and Drumm, 2003; Mendes et al., 2006; Eckert, 2007; Shears and Babcock, 2007; Dissanayake and Stefansson, 2010; Domínguez-Godino and González-Wangüemert, 2020). The study of these traits takes on a new importance considering that this is not only species-specific but varies according to habitat, i.e., the distribution patterns of a species may depend on the habitat it occurs in Sloan and von Bodungen (1980), Navarro et al. (2013), Domínguez-Godino and González-Wangüemert (2020). Depending on habitat type, densities of deposit-feeders may be higher in shallow-waters, but may vary in depth, for reasons like reduced exposure from UV radiation (Domínguez-Godino and González-Wangüemert, 2020). A vegetated cover can be the main driver for aggregation, as it provides shelter (Dissanayake and Stefansson, 2012; Navarro et al., 2014; Siegenthaler et al., 2015, 2017) or the opposite (Tuya et al., 2006), but in some cases the preference can change according to life-cycle stage (Mercier et al., 2000; Hamel et al., 2001; Eriksson et al., 2012; Aydin, 2019a). Rocky reefs can also be preferred when compared to adjacent areas without rocky substrate, offering more complex habitats, regarding food sources or shelter (Zhou and Shirley, 1996; Džeroski and Drumm, 2003; Mendes et al., 2006). Generally, density is not related to organic content of the sediment, although it can be related to its granulometry and, consequently, to diet quality through benthic biodiversity (Sloan and von Bodungen, 1980; Mercier et al., 2000; Dissanayake and Stefansson, 2012; Navarro et al., 2014; Domínguez-Godino and González-Wangüemert, 2020; Viyakarn et al., 2020). Muddy sediments, with more organic content, have, typically, less microbenthic diversity than sandy sea beds (Underwood and Barnett, 2006). Often, the main habitat preferences stem on lower hydrodynamic characteristics of the area and the presence of shelter (Sloan and von Bodungen, 1980; Mendes et al., 2006; Al-Rashdi et al., 2007; Morgan, 2011).
However, and in spite of the ecological relevance of these species and risk of overexploitation, these traits are fairly described for the most relevant species in the trade market (e.g., Mercier et al., 2000; Shiell, 2007; Morgan, 2011; Dissanayake and Stefansson, 2012; Purcell et al., 2012, 2016; Kashio et al., 2016) but for others are sparse and scarce. Emergent target species (from the Mediterranean and NE-Atlantic) are less studied, or not at all, or have a strong regional or habitat focus (e.g., Simunovic et al., 2000; Navarro et al., 2014; Siegenthaler et al., 2015, 2017; Marquet et al., 2017; Aydin, 2019b; Boncagni et al., 2019; Domínguez-Godino and González-Wangüemert, 2020). The introduction of these species in the international trade markets is the result of a stock depletion of traditional commercial species, mainly from the Indo-Pacific, increasing demand and consequent fisheries expansion to new areas (Purcell et al., 2013; Conand, 2017; González-Wangüemert et al., 2018; Dereli and Aydın, 2021). Holothuria mammata Grube, 1840, is currently one of the most relevant commercial species in the NE-Atlantic and Mediterranean (González-Wangüemert et al., 2014, 2018; Dereli and Aydın, 2021). Its geographical distribution covers the NE-Atlantic where it has its northern limits at the Bay of Biscay and southern limits at the meridional Macaronesian islands including the insular regions of Madeira, Azores and Canary islands, the north-African coast, occurring also in the Mediterranean sea (Tortonese, 1965; Borrero-Pérez et al., 2009). Still, the knowledge on ecological traits for this species is scarce, limited and regional, and insufficient to establish stock management policies, although relevant efforts have been put into it (Navarro et al., 2013; González-Wangüemert et al., 2014, 2016, 2018; Aydin and Erkan, 2015; Mustapha and Hattour, 2017; Siegenthaler et al., 2017; Azevedo E Silva et al., 2018; Félix et al., 2018; Simões et al., 2018; Aydin, 2019a,b, 2020; Sousa et al., 2020; Venâncio et al., 2021; Azevedo e Silva et al., 2021). There are several ecological details missing, either to allow management and control of harvest activities or for future mitigation actions involving population recovery or replenishment. Hence, there is a need for further studies on reproductive parameters (across the species distribution range), growth, recruitment, population structure and dynamics, meta-population variability, and other life-history traits.
This study aims to determine the temporal and spatial patterns of distribution of H. mammata, considering its abundance and demographic structure in a NE-Atlantic area, SW Portugal. Individual size and density will be modelled as a function of the site’s environmental characteristics to describe habitat preferences.
Materials and Methodology
Study Area
The study was carried out at a NE-Atlantic coastal area in the southwest of Portugal, within a marine protected area (MPA), the Arrábida Marine Park (38°26′50.4″N; 9°01′58.7″W). The study area is dominated by a rocky reef that ends in a sandy sea floor in a gradual transition around 3–10 m in depth, depending on site, as depth varies greatly throughout this coastline. This subtidal area displays a variety of habitats in complex diversity of macro- and microhabitats, which supports a high diversity of algae, invertebrates and fish that benefit from the hydro- and geomorphological conditions of the area. The productivity of these coastal waters increases in the summer, as a consequence of upwelling events (Wooster et al., 1976; Costa et al., 2013). Sampling was extended offshore up to the 20-meter bathymetric mark and the Sado estuary to delimit the actual range of occurrence of H. mammata (Figure 1). This is a currently unexploited population. Fisheries are mostly banned in this MPA, where local marine authorities display a regular surveillance for illegal activities, and sea cucumbers, in particular, are not yet a fishing target.
FIGURE 1
Figure 1. Map of the study area, with sampling sites, grouped according to habitat type, classified as; estuary (yellow); coastal (red); and offshore (blue).
Sampling
Considering the lack of knowledge on the distribution of H. mammata extensive preliminary surveys were conducted considering habitat type, depth and area to identify the distribution of the species in the region (Figure 1 – all sites). These surveys, that allowed to determine which areas should be considered for the bulk study, continued throughout a 1-year period, at each of its four seasons (in January, May, August, and November), to account for temporal variability and guarantee the sites with no occurrences were true zeroes. The majority of estuarine sites go from muddy to sandy sediment, without cover, except for sites SE7 and SE8 that were set within subtidal patches of Zostera marina, Linnaeus. Sites SE9 and SE10 were also exceptions regarding the spatial distribution of sampling sites. They are at a close distance but have markedly different properties: SE10 is the only site within the estuarine area with a subtidal isolated rocky outcrop that resembles the coastal rocky reefs; and SE9 is an open area without shelter or biotic cover of the substrate. The marine area (offshore and southern coast) is dominated by a sandy substrate regardless of depth. Sites SC2 to SC6 differ from the remaining sites by the presence of a rocky reef (Supplementary Table 1).
Based on these surveys, each site with a recorded presence of H. mammata was included in the distribution range of the species. This distribution area included the sites SC2 to SC6 and SE10, which were then specifically set based on spatial heterogeneity considering hydromorphological characteristics (exposure to tidal currents), available substrates (sand, rock, density of algae cover in each of the previous), shelter frequency (available crevices in rock and in sand/rock transition areas) species densities and distance to the estuary. As this area of occurrence was the one used to assess the distribution patterns of the species, sampling effort was increased to 1.5-month intervals, between the months of January 2018 and March 2019 (1.5 sampling intervals were set to match the seasonal samplings that occur every 3 months). Thus, for sites SC2 to SC6 and SE10, 9 sampling campaigns were performed, while for all other sites, only four sampling campaigns were done, in total.
To ensure comparability across the study area sampling campaigns were always based on a catch per unit of effort (cpue) of number of individuals per unit of area, always carried out between 9 to 12 h (a.m.) and with three replicates per site. Sampling was primarily done by scuba-diving (visual census on random transects of 30 m × 3 m, parallell to the coastline), except for sites at depths higher than 12 m, with high currents or in upper estuarine areas, where visibility did not allow a visual census. In these sites, sampling was carried by beam-trawl with an opening of 2 m and a 2 mm mesh size at the cod end. Depending on location, tows of 5–10 min, also with three replicates per site, were performed. Typically, longer tows were possible offshore (higher depths and sandy bottoms) and shorter tows inside the estuary, with muddy bottoms and debris that promote a quick trawl clogging (Supplementary Table 1). Sites SC2 to SC6 and SE10 (species distribution area), unlike every other site, had two bottom types (rock and sand). For these, six random transects were done (three per substrate type) to represent the full heterogeneity of the site. The first transects were based on the deployment location and the following to the left and right of that location, parallell to the coastline. Transect areas were not repeated along the study period.
All H. mammata individuals detected in each transect were counted, and measured (nearest mm ±1), in situ, without manipulation to prevent muscle contraction and minimise size variability (using a flexible measuring tape). Since rocky substrates have sandy patches and vice-versa, the surface on which each individual occurred was also recorded (e.g., presence of an individual in sand substrate/patch of a rock transect), as was the cryptic condition of each individual (within macroalgal cover or sheltered in crevices and seams in rocks).
For each site (coastal, estuarine, and offshore) and each survey several environmental parameters of interest were measured, using a calibrated multiparametric sonde YSI-EXO2: water temperature, pH, salinity, dissolved oxygen (ODO), dissolved solids, turbidity, chlorophyll a, and depth. Site current was measured with a Doppler Current Sensor 4100. The type of rock surface for each site (plain or presenting crevices and seams) and the type of biological cover of the substrate (seagrass, macroalgae, or none) were also recorded. Grab samples were collected at each site for determination of sediment grain size and total organic matter (TOM) of sediments. The later was obtained by loss on ignition (480°C) and the former determined using 63 μm to 2 mm sieves to separate, respectively, the silt, sand and gravel fractions. Each fraction was then dried and weighed, and sediments assorted according to their percentages (Blott and Pye, 2001) and then classified with Shepard diagrams (Shepard, 1954).
Data Analysis
The data analysis was divided in two stages: (1) the exploratory assessment of occurrence areas for H. mammata and the identification of space and time differences in the distribution of the population; and (2) the use of empirical models to evaluate which environmental factors influence population density and size distribution of H. mammata. All statistical analyses were carried out in R software (R Core Team, 2020) and results considered significant at a p-value < 0.05.
To assess the environmental characteristics of the study area and to help explain the distribution range of H. mammata, a Principal Components Analysis (PCA) was used to group sites by their environmental parameters, considering three main habitats, classified as: estuary, coast, and offshore. The existence of spatial and temporal differences in the distribution of H. mammata, based on density and mean size, was assessed with a Kruskal–Wallis rank test. The selection of non-parametric techniques was justified by a non-normal distribution of density and a positively skewed distribution of mean size. Multiple comparisons were then conducted with the Hochberg method, for an adjusted p-value (Hochberg, 1988). An ANOVA was used to determine differences between sheltered and unsheltered behaviour (e.g., individuals in crevices or algae cover), based on individual sizes (normally distributed).
The second stage of the data analysis relied on empirical modelling to explain distribution patterns (density and mean size) according to environmental factors. To avoid misleading inferences due to similar effects of correlated variables (Morrissey and Ruxton, 2018), the collinearity between predictors to be included in the empirical models was assessed by applying a stepwise trait selection based on variance inflation factors (VIF), using the function vifstep() in the R package usdm (Naimi et al., 2014), with a VIF threshold set at th = 8, and checked with correlation plots. The selected predictors are described in Table 1. Apart from the variables excluded due to collinearity issues (dissolved solids, fine and medium sand fractions, and gravel), the two following variables were also left out. Salinity in the coastal area is more conservative than in estuarine environments. Since there is only one site at the estuary mouth (SE10) this could result in local salinity variations based on tidal regimes. Sea cucumbers are low mobility animals and salinity variations that are tidal dependent make a poor predictor for distribution patterns. Thus, other variables were chosen to represent a proximity to the estuary (e.g., longitude, turbidity, and ODO). Besides, longitude is not influenced by tidal regimes and salinity is, which, in turn, varies according to sampling hour and day (sampling logistics did not allow a standardised sampling regarding the tide). Shelter (presence or absence of crevices or seams where sea cucumbers can find shelter in), a binary variable representing a response to physical environmental stressors, like current, or predation, was also removed, because it was highly correlated to site current and caused convergence issues in the models. Sites without shelter were, simultaneously, the sites with the highest exposure to hydrodynamic influence (high current).
TABLE 1
Table 1. Selected variables based on variance inflation factors to be included in empirical models assessing habitat preferences for Holothuria mammata.
The Generalized Linear Models (GLM) were implemented for two response variables: (1) density and (2) mean size, with Tweedie and Gamma distributions, respectively, and log-link functions (Tweedie, 1984). The ΔAIC was the measure derived from the AIC (Akaike information criterion) that was used as the model selection criterion and significant interactions were assessed considering the signs of the coefficients for the interaction and individual variables and with surface plots (Feld et al., 2016). A thorough description on the used models and the analytical procedures that supported the interpretation of the results are detailed in Supplementary Material 2. For relevant predictors in the models that showed a variations in time, a causality test, based on the ccf() cross-correlation function, was used to determine relationships between two time series: the relevant predictor variable and the response variable, across the sampling period. The autocorrelation function (ACF) plotted the results with all tested time-lags.
Results
Space-Time Distribution Patterns
The distribution range of H. mammata in the sampled area, like the other two commercial holothurians (Holothuria forskali and Holothuria arguinensis), was restricted to rocky reefs. Marine open areas without hard substrate had no occurrences and at rocky reefs no individuals were found beyond a maximum distance of 7 m from the transition border between rock and sand substrates. The estuarine area showed a wider heterogeneity concerning its environmental features and the dissimilarities from the adjacent marine environment were, mainly, related to the granulometric characteristics of the sediment, its organic content, turbidity, dissolved oxygen and temperature (Figure 2). The only estuarine site with H. mammata (or any of the aforementioned) was SE10, a site near the estuary mouth distinguished by the presence of a rocky outcrop and, in more than a year of sampling, no occurrence was observed at the two low depth nearby sites, SE9 and SE8, the latter in a seagrass meadow patch.
FIGURE 2
Figure 2. Principal Component Analysis (PCA) diagram of environmental data by sampled site, grouped according to habitat, classified as coast, estuary, and offshore.
Within the area of distribution of H. mammata 2114 individuals, ranging from 119 to 510 mm, were observed and counted in a total sampled area of 29,160 m2. Density per replicate ranged from 0 to 120 ind/100 m2 (or 12,000 ind/ha – a maximum value recorded at SC2) (Figure 3). The majority of individual counts were made on rock transects (94%), when compared to sand transects. Of those, 68% were settled on a sandy patch. In sand transects, the large majority was, in fact, settled on sand (82%) and remaining on small rocky patches. However, it should be considered that sand transects have very few rock outcrops and in rock transects sandy patches are common. Only a small number of individuals preferred a macroalgal cover, as only 5% occurred in the algal cover of rock transects and none in the algal cover of sand transects. Of all specimens, 34% were found to be sheltered in crevices at the time of sampling (daytime) and these individuals were overall smaller in size (ANOVA: p = 0.0086). When comparing sizes between substrates of settlement, individuals occurring on rock were also smaller than those on sand (ANOVA: p < 0.0001) and so were those occurring within an algal cover (ANOVA: p = 0.039).
FIGURE 3
Figure 3. Temporal density data for Holothuria mammata at Arrábida, Setúbal, during the nine month sampling period. Red dots at each month connected by a line represent the evolution of the mean.
Both density and mean size showed a significant spatial variation (KW: p = 2.474e-14 and p = 1.408e-13, respectively), but not temporal, although overall average density increased from spring up until early autumn (Figure 3). However, spatial differences in density were not related to site proximity (Figure 4). The most similar pair in density was SC5 and SE10 (KW, Hochberg: adjusted p = 1) and the most dissimilar pair was SC2 and SC3 (KW, Hochberg, adjusted p = 1.003e-10). Mean size, however, tended to increase towards the estuary, considering the four most distant sites (Figure 5). Size-classes presented a unimodal distribution at all sites. Mean size and density have a negative a linear relationship, with higher densities corresponding to lower size classes (R2 = 0.12, F1,155 = 20.4, p = 1.24e-05).
FIGURE 4
Figure 4. Density of Holothuria mammata by sampling site, arranged according to distance to estuary (left to right). Red dotted line represents the overall mean value. * represents the mean of each sampling station.
FIGURE 5
Figure 5. Size-class distribution of Holothuria mammata by sampling site, arranged according to distance to estuary (bottom to top). Vertical lines represent the site mean.
Density Distribution Models
The GLM used to explain the density distribution of H. mammata, rendered only one final model with a ΔAIC > 2 and adjusted R2 of 62% and an explained deviance of 68.6% (Table 2). Density, consistently higher at rock transects, decreased with increasing current, with pH (between the values of 7.71 and 8.54), with proximity to the estuary and increased with depth (Figure 6). The pH did not show a spatial variation, but a temporal one, increasing in the rainy months. The cross-correlation function showed a significant negative correlation at lag = 0 (months) (Figure 7A), with the peak of pH corresponding to the lowest average density value (Figure 7B). The significant interaction between depth and current, with a positive coefficient contrasting with the negative coefficients of its individual variables, showed an antagonistic effect, i.e., depth cancelled the negative effect of current on density. The interaction plot (Figure 8), detailing the effect, shows that between the depth values of 6 and 8 m, there was little variation in density, as opposed to lower depth values where the current had an effect on density, which decreased as current increased. In other words, the intensity of current only had a negative influence on density at low depths. At higher depths the hydrological characteristics of the site had less impact on density, as it even increased with increasing current at the two deepest sites.
TABLE 2
Table 2. Fixed-effect GLM model results with a Tweedie distribution (R2 = 62%), explaining density distribution patterns of Holothuria mammata at Arrábida, Setúbal.
FIGURE 6
Figure 6. Partial plots of the density model for all significant predictors – Fixed-effect GLM with a Tweedie distribution – explaining distribution patterns of Holothuria mammata at Arrábida, Setúbal.
FIGURE 7
Figure 7. Cross-correlation analysis for the time series of pH of the water column and mean density of Holothuria mammata, at Arrábida, Setúbal, during the nine month sampling period: (A) ACF plot depicting correlations (vertical lines) and significance threshold (horizontal dotted lines); (B) average monthly values for pH and density.
FIGURE 8
Figure 8. Surface plot depicting the interaction behaviour between current and depth in the density model, explaining distribution patterns of Holothuria mammata at Arrábida, Setúbal.
Size Distribution Models
The GLM used to explain the size distribution rendered four top models with a ΔAIC < 2. Model averaging for these models (full and conditional) revealed different levels of importance for several predictors, with substrate and distance to estuary (longitude) as the most relevant variables (Table 3 – full average), and coarse sand, ODO, TOM, Chla, Turbidity and the interaction depth:current as less important in explaining the variance of mean size for this species (Table 3 – conditional average). Considering the sum of weights of each variable in all four models, ODO and depth still presented a relatively relevant contribution in explaining the variance of the response variable. Coarse sand, TOM, the interaction depth:current and turbidity were relatively similar in contribution to the final model, but Chla was only represented in one model (Table 4). These results were corroborated by the assessment of the confidence intervals of the predictor’s coefficients (Figure 9). Nonetheless, of all the predictors regarded as less important, coarse sand and the interaction depth:current provided a more precise estimate of effect (narrower confidence intervals). Depth, although with a high importance (sum of weights), crosses zero and shows a wider confidence interval and, hence, not significant in either averaged model (full or conditional). Considering the interpretation of all model validation methods, the key predictors explaining size were substrate (larger individuals on sand) and longitude (larger individuals closer to the estuary). To a lesser extent, coarse sand (larger individuals associated to the higher granulometric fraction of sand), TOM (larger individuals at sites with more organic content in sediment: values between 0.3 and 8.5%), and turbidity (larger individuals at lower turbidity sites) also seem relevant. The relationship of larger individuals with more productive sites (Chla) was the least important association. The interaction depth:current was not assessed as to its type due to the unsignificant relevance of its individual predictors in the model. Of the predictors, ODO, TOM, and Chla, only the latter showed an obvious pattern, although like pH, it was not spatial, but temporal, with higher average values in the summer months. The cross-correlation function showed a significant positive correlation at lag = −1 (months) between Chla and mean size (Figure 10A), depicting an increase in local productivity preceded by a peak in mean length (Figure 10B).
TABLE 3
Table 3. Fixed-effect GLM averaged model results with a Gamma distribution (averaged R2 = 49%, as the average of all models weighted by their model weights), explaining size distribution patterns of Holothuria mammata at Arrábida, Setúbal.
TABLE 4
Table 4. Model importance values extracted from the GLM averaged model, considering the four top models with a ΔAIC < 2.
FIGURE 9
Figure 9. Confidence intervals for the predictor coefficients of the model – Fixed-effect GLM averaged model with a Gamma distribution – explaining the size distribution patterns of Holothuria mammata at Arrábida, Setúbal.
FIGURE 10
Figure 10. Cross-correlation analysis for the time series of Chla and mean total length of Holothuria mammata, at Arrábida, Setúbal, during the nine month sampling period: (A) ACF plot depicting correlations (vertical lines) and significance threshold (horizontal dotted lines); (B) average monthly values for Chla and mean total length.
Discussion
The population of H. mammata at this NE-Atlantic rocky reef was the healthiest studied up to this moment. It showed the largest individuals reported so far, with a maximum total length of 51 cm, densities up to 120 ind/100 m2 and unimodal size-class distributions. This reflects the good environmental conditions, successful recruitment and an unexploited population. The Mediterranean and other NE-Atlantic regions show populations of this species with lower maximum sizes, the majority between 25 and 38 cm (Navarro et al., 2013; González-Wangüemert et al., 2014, 2018; Marquet et al., 2017; Mustapha and Hattour, 2017; Siegenthaler et al., 2017). Previous records show that H. mammata in the NE-Atlantic attain bigger sizes than in the Mediterranean and the largest in the Portuguese coast, with Peniche, a highly productive marine area, presenting the previous largest documented specimen, with 43 cm (González-Wangüemert et al., 2016). Densities found in the present study were also the highest, compared to a maximum density of 55 ind/100 m2 found in the Mediterranean, Aegean Sea (Aydin, 2019a), 16 ind/100 m2 reported in the Canary Islands (Navarro et al., 2013), and particular low densities observed at Ria Formosa (south Portugal) of 1.2 ind/100 m2 (Siegenthaler et al., 2017). However, densities tend to be much lower in habitats without rocky bottoms (Navarro et al., 2013) and Ria Formosa is a shallow mesotidal lagoon lacking the typical H. mammata habitat (Siegenthaler et al., 2017). Like at Peniche, the Arrábida coast is also a highly productive area, benefiting from upwelling events in the summer, nutrient inputs from the estuary and from a complex subtidal habitat (Wooster et al., 1976; Cabeçadas et al., 1999; Costa et al., 2013) that can have positive effects on the condition of the population. The high productivity together with the subtidal distribution of the species, unlike intertidal distributions (Siegenthaler et al., 2017), provide more stable environmental conditions, which may be the reason for higher abundances and lower variation in time.
This population was not only restricted to the marine environment but it preferred rocky reefs as opposed to other complex environments such as stands of macroalgae and seagrass meadows. Siegenthaler et al. (2017) had already reported that the preference of H. mammata for seagrass meadows at Ria Formosa did not reflect the species’ habitat preferences, but a distribution driven by the lack of rocky substrata and the consequent search for shelter and a more complex habitat (González-Wangüemert et al., 2014, 2016). Despite the higher occurrence found on rocky substrate, smaller animals were the ones sheltered in crevices or macroalgae, contrary to those on sand. There are a number of sea cucumber predators, with sea stars as the most common, particularly at the juvenile stage (Francour, 1997). Hence, those larger that venture more into open spaces (for feeding activities) are less prone to be preyed upon, as predation risk declines with growth (Shiell and Knott, 2008; Purcell, 2010; Ceccarelli et al., 2018). It is important to stress that while sandy patches in rocky areas have shelter close by, sand areas do not, which may be the cause for the proportion of individuals on sand to be higher when considering the individual presence, rather than the transects density (sandy areas vs. rock). So, this may be either a defence strategy for smaller individuals, that use sandy patches to feed, and an opportunistic behaviour exclusive of larger animals that move more into sandy areas, or instead, a diel feeding behaviour of lower size-classes that may be more active at night, as this is a species with reported nocturnal activity (Borrero-Pérez et al., 2010; Navarro et al., 2013; Siegenthaler et al., 2017). Nighttime was not covered in the present study, but there was a large proportion of active individuals at this NE-Atlantic rocky reef during the day. This was not only assessed by the number of individuals on sand, because H. mammata feeds exclusively on sand in this area (Azevedo e Silva et al., 2021) and needs to move from the rock substrate to a sandy area or patch to feed, but also demonstrated by the observation of feeding behaviour at sampling time for several individuals (footage obtained at the study area between 10 and 12 h a.m. at https://seacucumber.eu/en/o-pepino-do-mar).
Like other sea cucumber species (Džeroski and Drumm, 2003; Mendes et al., 2006; Eckert, 2007; Shears and Babcock, 2007; Dissanayake and Stefansson, 2010; Domínguez-Godino and González-Wangüemert, 2020), the population studied in the present work showed a patchy distribution, with areas of high abundance followed by zero or near zero sites. However, this patchiness is size dependent. Lower size-classes aggregated in higher density groups in more stable environments (away from the estuary).
The models results suggest that abiotic and biotic parameters of the water column are not the main drivers shaping the distribution of H. mammata in the area. Similarly to other sea cucumber species (e.g., Sloan and von Bodungen, 1980; Sonnenholzner, 2003; Mendes et al., 2006; Morgan, 2011; Domínguez-Godino and González-Wangüemert, 2020), its habitat preferences depend on environmental stability, as opposed to the estuarine variability, low hydrodynamic conditions and the presence of rocky substrate that may offer shelter and habitat complexity, more likely to increase food quality. Depth plays an important role by providing a quieter refuge in areas that are more susceptible to tidal influence. Individuals tend to use bathymetry to their advantage in areas where current is stronger, occupying deeper areas. The significant physical-chemical predictor of the water column explaining density was pH. However, pH shows little variation in space, but is highly variable in time. Its increase in the autumn and winter months is most likely related to continental run-off from the limestone ridge of Arrábida that characterises the entire coastline (Kullberg et al., 2012). This temporal variation showed a negative correlation with the mean density that, within this time-frame, represents an aggregation phenomenon that increases capturability during the months that coincide with the pH decrease in this coastal area. In echinoderms, pH may affect abundance, distribution, and also reproductive success, as a consequence of induced physiological stress (Lawrence and Herrera, 2000; Walag and Canencia, 2016). Coincidently, the highest densities found occurred during spring/summer months, which correspond to the species reproductive season (unpublished data for the area; Marquet et al., 2017), and may reflect a relationship between reproduction and suitable environmental conditions, where pH can play a role, promoting an aggregation behaviour during this period (Hamel and Mercier, 1996b; Mercier and Hamel, 2009; Leite-Castro et al., 2016). The model results corroborate the relationship between density and size. Despite the lower densities on sand, these were represented by the largest individuals. Smaller individuals were distributed in rocky substrate and at an increasing distance from the estuary, where, by comparison, there is less environmental variability. The relationships with H. mammata sizes and dissolved oxygen, organic content in the sediment and Chla were, in general, difficult to assess due to the low importance of these predictors in the model as well as their inconspicuous patterns. Only Chla showed an evident temporal variation, increasing in the summer, which is a common phenomenon in the area, due to seasonal upwelling cycles (Wooster et al., 1976). Similarly to pH and density, this rising local productivity is associated an aggregation behaviour of larger size-classes in the spring/summer months. Hence, there is a strong suggestion of an aggregation pattern of breeders in a specific period, which relates to environmental conditions that potentially favour the physiological process of gonad development. The species’ preference for coarser sand fractions may be related to a more diverse diet and less related to organic sediment content (Underwood and Barnett, 2006). So, these associations may be explained by the displacement ability of larger individuals towards more favourable conditions, here expressed by coarser sediment fractions, higher productivity, less turbidity and more organic content of the sediment (never higher than 8.5%), at a lower risk of predation or due to a higher resilience, both privileged by size.
Conclusion
In general, H. mammata prefers a more stable environment, which is more important for smaller individuals, while bigger size classes tend to move to less stable environments, probably in an opportunistic fashion due to their size related resilience. This behaviour could be related to feeding conditions associated to lower predation susceptibility of bigger animals, or with a reproductive behaviour leading the broodstock, as they reach maturity size, towards more favourable areas for reproduction, since sea cucumber larvae are pelagic (Venâncio et al., 2021) and require currents for dispersion (Pedrotti and Fenaux, 1992).
Conducting research in an unexploited area provided important advantages in the study of biological and ecological traits, closer to a pristine condition, in comparison to populations that are subject to harvest, several of which already showing signs of exploitation by a reduction of larger size classes, thus creating multimodal distributions (e.g., González-Wangüemert et al., 2014, 2018; Aydin and Erkan, 2015). Although, this work was conducted in a limited area of an Atlantic rocky-reef, there is a consistency in the drivers that shape the distribution of different detritivorous sea cucumbers, when compared to other studies (e.g., Sloan and von Bodungen, 1980; Sonnenholzner, 2003; Mendes et al., 2006; Morgan, 2011; Domínguez-Godino and González-Wangüemert, 2020). These species tend to prefer the available shelter, a more complex habitat for feeding and areas with less environmental variability. This provides guidance for future stock management measures in a time when only a small proportion of the species’ distribution range is studied. At the studied rocky coastal area, the habitat preferences of H. mammata showed that, at this point, even lacking important information like growth parameters or size at first maturity to allow the establishment of recovery periods for the population and determine, e.g., fishing quotas for each size-class, some red flags can be raised as to the potential susceptibility of H. mammata to fisheries. The patchy distribution and the high density areas should be regarded as sensitive areas, as these are more attractive for fishermen, allowing a reduced fishing effort and increased yields, which can lead to a significant reduction in spawning biomass. Access to high density of breeders should also be restricted to allow a spillover effect and help sustain the adjacent populations (Purcell and Kirby, 2006). On an undesirable note, and to the best of our knowledge, no study, so far, covers the distribution of wild juveniles of H. mammata under 100 mm. The understanding of juvenile habitat preferences, distribution and the integration of these size-classes in growth models is key knowledge and, currently, a large gap. During the sampling campaigns no small juveniles were observed, which suggests a different habitat preference or, most likely, a cryptic or diel behaviour, as observed for other sea cucumber species (e.g., Sloan, 1979; Shiell, 2004; Purcell, 2010; Soliman et al., 2019).
Although, currently, the exploitation of NE-Atlantic sea cucumbers is restricted in space and to a relatively low number of illegal fishermen, the high-income potential that these animals yield will soon widen the interest that easily leads to a high fishing pressure. This study represents a fundamental step towards a better understanding of the populations’ dynamics and can contribute with vital knowledge support decision-making processes on stock management.
Data Availability Statement
The datasets presented in this article are not readily available because these are under an embargo period. Requests to access the datasets should be directed to PF, pmfelix@fc.ul.pt.
Author Contributions
PF, AB, and JC designed the study. FA, TS, and PF carried out the field work. PF, TM, and CR developed the data analysis. AP, RM, JS, and EV provided the contributions to several aspects of the study and the manuscript. All co-authors revised the manuscript.
Funding
This study was financed by the Operational Programme Mar2020, MAR-02.01.01-FEAMP-0052, “Newcumber – Avanços para o cultivo sustentável de pepinos do mar”. It received further financial support from Fundação para a Ciência e Tecnologia (projects UIDB/04292/2020, UIDB/00006/2020, CoastNet – PINFRA/22128/2016, AB with the Scientific Stimulus Programme – CEECIND/00095/2017 and FA with the individual research grant 2020.09563.BD). This publication was financed by Portuguese national funds through FCT – Fundação IP under project reference UIDB/04292/2020, and by the European Union’s Horizon 2020 Research and Innovation Program under grant agreement N810139: Project Portugal Twinning for Innovation and Excellence in Marine Science and Earth Observation – PORTWIMS.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors would like to thank Ricardo Lima for developing the map for the study area and Pedro Segurado for his crucial help with the interaction plots.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2021.675330/full#supplementary-material
References
Al-Rashdi, K. M., Claereboudt, M. R., and Al-Busaidi, S. S. (2007). Density and size distribution of the sea cucumber, Holothuria scabra (Jaeger, 1935), at Six Exploited Sites in Mahout Bay, Sultanate of Oman. J. Agric. Mar. Sci. 12, 43–51. doi: 10.24200/jams.vol12iss0pp43-51
Aydin, M. (2019a). Biometry, density and the biomass of the commercial sea cucumber population of the Aegean Sea. Turkish J. Fish. Aquat. Sci. 19, 463–474. doi: 10.4194/1303-2712-v19_6_02
Aydin, M. (2019b). Density and biomass of commercial sea Cucumber species relative to Depth in the Northern Aegean Sea. Thalassas 35, 541–550. doi: 10.1007/s41208-019-00144-4
Aydin, M. (2020). Length-weight relationships and condition factor of four different sea cucumber species in the Aegean Sea. J. Anatol. Environ. Anim. Sci. 5, 80–85. doi: 10.35229/jaes.677940
Aydin, M., and Erkan, S. (2015). Identification and some biological characteristics of commercial sea cucumber in the Turkey coast waters. Int. J. Fish. Aquat. Stud. 3, 260–265.
Azevedo e Silva, F., Brito, A. C., Simões, T., Pombo, A., Marques, T. A., Rocha, C., et al. (2021). Allometric relationships to assess ontogenetic adaptative changes in three NE Atlantic commercial sea cucumbers (Echinodermata, Holothuroidea). Aquat. Ecol. 55, 711–720. doi: 10.1007/s10452-021-09856-3
Azevedo E Silva, F., Simões, T., Raposo, A., Sousa, J., Venâncio, E., Pombo, A., et al. (2018). Reproductive biology of north-eastern Atlantic sea cucumber species: a preliminary assessment. Front. Mar. Sci. 5:66. doi: 10.3389/conf.FMARS.2018.06.00066
Blott, S. J., and Pye, K. (2001). GRADISTAT: a grain size distribution and statistics package for the analysis of unconsolidated sediments. Earth Surf. Process. Landforms 26, 1237–1248. doi: 10.1002/esp.261
Boncagni, P., Rakaj, A., Fianchini, A., and Vizzini, S. (2019). Preferential assimilation of seagrass detritus by two coexisting Mediterranean sea cucumbers: Holothuria polii and Holothuria tubulosa. Estuar. Coast. Shelf Sci. 231:106464. doi: 10.1016/j.ecss.2019.106464
Borrero-Pérez, G. H., Gómez-Zurita, J., González-Wangüemert, M., Marcos, C., and Pérez-Ruzafa, A. (2010). Molecular systematics of the genus Holothuria in the Mediterranean and Northeastern Atlantic and a molecular clock for the diversification of the Holothuriidae (Echinodermata: Holothuroidea). Mol. Phylogenet. Evol. 57, 899–906. doi: 10.1016/j.ympev.2010.08.019
Borrero-Pérez, G. H., Pérez-Ruzafa, A., Marcos, C., and González-Wangüemert, M. (2009). The taxonomic status of some Atlanto-Mediterranean species in the subgenus Holothuria (Echinodermata: Holothuroidea: Holothuriidae) based on molecular evidence. Zool. J. Linn. Soc. 157, 51–69. doi: 10.1111/j.1096-3642.2009.00529.x
Cabeçadas, G., Nogueira, M., and Brogueira, M. (1999). Nutrient dynamics and productivity in three european estuaries. Mar. Pollut. Bull. 38, 1092–1096. doi: 10.1016/s0025-326x(99)00111-3
Ceccarelli, D., Logan, M., and Purcell, S. (2018). Analysis of optimal habitat for captive release of the sea cucumber Holothuria scabra. Mar. Ecol. Prog. Ser. 588, 85–100. doi: 10.3354/meps12444
Conand, C. (2017). Expansion of global sea cucumber fisheries buoys exports. Rev. Biol. Trop. 65, S1–S10. doi: 10.15517/rbt.v65i1-1.31661
Costa, B. H. E., Erzini, K., Caselle, J. E., Folhas, H., and Gonçalves, E. J. (2013). Reserve effect within a temperate marine protected area in the north-eastern Atlantic (Arrábida Marine Park, Portugal). Mar. Ecol. Prog. Ser. 481, 11–24. doi: 10.3354/meps10204
Dereli, H., and Aydın, M. (2021). Sea cucumber fishery in Turkey: management regulations and their efficiency. Reg. Stud. Mar. Sci. 41:101551. doi: 10.1016/j.rsma.2020.101551
Dissanayake, D. C. T., and Stefansson, G. (2010). Abundance and distribution of commercial sea cucumber species in the coastal waters of Sri Lanka. Aquat. Living Resour. 23, 303–313. doi: 10.1051/alr/2010031
Dissanayake, D. C. T., and Stefansson, G. (2012). Habitat preference of sea cucumbers: Holothuria atra and Holothuria edulis in the coastal waters of Sri Lanka. J. Mar. Biol. Assoc. U.K. 92, 581–590. doi: 10.1017/S0025315411000051
Domínguez-Godino, J. A., and González-Wangüemert, M. (2020). Habitat associations and seasonal abundance patterns of the sea cucumber Holothuria arguinensis at Ria Formosa coastal lagoon (South Portugal). Aquat. Ecol. 54, 337–354. doi: 10.1007/s10452-020-09746-0
Džeroski, S., and Drumm, D. (2003). Using regression trees to identify the habitat preference of the sea cucumber (Holothuria leucospilota) on Rarotonga, Cook Islands. Ecol. Modell. 170, 219–226. doi: 10.1016/s0304-3800(03)00229-1
Eckert, G. L. (2007). Spatial patchiness in the sea cucumber Pachythyone rubra in the California Channel Islands. J. Exp. Mar. Biol. Ecol. 348, 121–132. doi: 10.1016/j.jembe.2007.04.004
Eriksson, H., Jamon, A., and Wickel, J. (2012). Observations on habitat utilization by the sea cucumber Stichopus chloronotus. SPC Beche-de-mer Inf. Bull. 32, 39–42.
Feld, C. K., Segurado, P., and Gutiérrez-Cánovas, C. (2016). Analysing the impact of multiple stressors in aquatic biomonitoring data: a ‘cookbook’ with applications in R. Sci. Total Environ. 573, 1320–1339. doi: 10.1016/j.scitotenv.2016.06.243
Félix, P., Pombo, A., Azevedo, E., Silva, F., Simões, T., Sousa, J., et al. (2018). Newcumber: advances towards the sustainable rearing of sea cucumbers. Front. Mar. Sci. 5:60. doi: 10.3389/conf.FMARS.2018.06.00060
Francour, P. (1997). Predation on holothurians: a literature review. Invertebr. Biol. 116, 52–60. doi: 10.2307/3226924
González-Wangüemert, M., Aydin, M., and Conand, C. (2014). Assessment of sea cucumber populations from the Aegean Sea (Turkey): first insights to sustainable management of new fisheries. Ocean Coast. Manag. 92, 87–94. doi: 10.1016/j.ocecoaman.2014.02.014
González-Wangüemert, M., Domínguez-Godino, J. A., and Cánovas, F. (2018). The fast development of sea cucumber fisheries in the Mediterranean and NE Atlantic waters: from a new marine resource to its over-exploitation. Ocean Coast. Manag. 151, 165–177. doi: 10.1016/j.ocecoaman.2017.10.002
González-Wangüemert, M., Valente, S., Henriques, F., Domínguez-Godino, J. A., and Serrão, E. A. (2016). Setting preliminary biometric baselines for new target sea cucumbers species of the NE Atlantic and Mediterranean fisheries. Fish. Res. 179, 57–66. doi: 10.1016/j.fishres.2016.02.008
Hamel, J.-F., Conand, C., Pawson, D. L., and Mercier, A. (2001). The sea cucumber Holothuria scabra (Holothuroidea: Echinodermata): its biology and exploitation as Beche-de-mer. Adv. Mar. Biol. 41, 129–223. doi: 10.1016/s0065-2881(01)41003-0
Hamel, J.-F., and Mercier, A. (1996a). Early development, settlement, growth, and spatial distribution of the sea cucumber Cucumaria frondosa (Echinodermata: Holothuroidea). Can. J. Fish. Aquat. Sci. 53, 253–271. doi: 10.1139/f95-186
Hamel, J.-F., and Mercier, A. (1996b). Evidence of Chemical Communication During the Gametogenesis of Holothuroids. Ecology 77, 1600–1616. doi: 10.2307/2265555
Hochberg, Y. (1988). A sharper bonferroni procedure for multiple tests of significance. Biometrika 75, 800–803. doi: 10.2307/2336325
Kashio, S., Yamana, Y., Furukawa, N., Uekusa, R., and Goshima, S. (2016). Seasonal microhabitat use patterns in Japanese sea cucumber Apostichopus japonicus in Funka Bay, Hokkaido, northern Japan. Aquac. Sci. 64, 371–378. doi: 10.11233/aquaculturesci.64.371
Kinch, J. P., Purcell, S., Uthicke, S., and Friedman, K. (2008). “Papua New Guinea: a hotspot of sea cucumber fisheries in the Western Central Pacific,” in Sea Cucumbers. A Global Review of Fisheries and Trade, eds V. Toral-Granda, A. Lovatelli, and M. Vasconcellos (Rome: FAO), 57–77.
Kullberg, J. C., Pais, J., Almeida, J. A., and Mateus, O. (2012). Contributo do Património Geológico e Geomorfológico na Candidatura da Arrábida (Portugal) a Património Mundial Misto. Santos.
Lawrence, J. M., and Herrera, J. (2000). Stress and deviant reproduction in echinoderms. Zool. Stud. 39, 151–171.
Leite-Castro, L. V., de Souza Junior, J., Salmito-Vanderley, C. S. B., Nunes, J. F., Hamel, J.-F., and Mercier, A. (2016). Reproductive biology of the sea cucumber Holothuria grisea in Brazil: importance of social and environmental factors in breeding coordination. Mar. Biol. 163:67. doi: 10.1007/s00227-016-2842-x
MacTavish, T., Stenton-Dozey, J., Vopel, K., and Savage, C. (2012). Deposit-Feeding Sea cucumbers enhance mineralization and nutrient cycling in organically-enriched coastal sediments. PLoS One 7:e50031. doi: 10.1371/journal.pone.0050031
Marquet, N., Conand, C., Power, D. M., Canário, A. V. M., and González-Wangüemert, M. (2017). Sea cucumbers, Holothuria arguinensis and H. mammata, from the southern Iberian Peninsula: variation in reproductive activity between populations from different habitats. Fish. Res. 191, 120–130. doi: 10.1016/j.fishres.2017.03.007
Mendes, F., Marenzi, A., and Domenico, M. (2006). Population patterns and seasonal observations on density and distribution of Holothuria grisea (Holothuroidea: Aspidochirotida) on the Santa Catarina Coast, Brazil. SPC Beche-de-mer Inf. Bull. 23, 5–10.
Mercier, A., Battaglene, S. C., and Hamel, J.-F. (2000). Periodic movement, recruitment and size-related distribution of the sea cucumber Holothuria scabra in Solomon Islands. Hydrobiologia 440, 81–100. doi: 10.1007/978-94-017-1982-7_8
Mercier, A., and Hamel, J. (2009). Endogenous and exogenous control of gametogenesis and spawning in echinoderms. Adv. Mar. Biol. 55, 1–302.
Michio, K., Kengo, K., Yasunori, K., Hitoshi, M., Takayuki, Y., Hideaki, Y., et al. (2003). Effects of deposit feeder Stichopus japonicus on algal bloom and organic matter contents of bottom sediments of the enclosed sea. Mar. Pollut. Bull. 47, 118–125. doi: 10.1016/s0025-326x(02)00411-3
Morgan, A. D. (2011). Patterns of distribution and abundance of the temperate sea cucumber Australostichopus mollis on a rocky subtidal reef. New Zeal. J. Zool. 38, 195–206. doi: 10.1080/03014223.2011.563426
Moriarty, D. J. W., Pollard, P. C., Hunt, W. G., Moriarty, C. M., and Wassenberg, T. J. (1985). Productivity of bacteria and microalgae and the effect of grazing by holothurians in sediments on a coral reef flat. Mar. Biol. 85, 293–300. doi: 10.1007/BF00393250
Morrissey, M. B., and Ruxton, G. D. (2018). Multiple regression is not multiple regressions: the meaning of multiple regression and the non-problem of Collinearity. Philos. Theory Pract. Biol. 10:3. doi: 10.3998/ptpbio.16039257.0010.003
Mustapha, K., and Hattour, A. (2017). Estimation de la population des holothuries des régions Jarzouna Metline, El Hawar (Ras Jbel) et la lagune de Bizerte. Bull. l’Institut Natl. Sci. Technol. Salammbô 44, 147–163.
Naimi, B., Hamm, N. A. S., Groen, T. A., Skidmore, A. K., and Toxopeus, A. G. (2014). Where is positional uncertainty a problem for species distribution modelling? Ecography 37, 191–203. doi: 10.1111/j.1600-0587.2013.00205.x
Namukose, M., Msuya, F. E., Ferse, S. C. A., Slater, M. J., and Kunzmann, A. (2016). Growth performance of the sea cucumber Holothuria scabra and the seaweed Eucheuma denticulatum: integrated mariculture and effects on sediment organic characteristics. Aquac. Environ. Interact. 8, 179–189. doi: 10.3354/aei00172
Navarro, P. G., García-Sanz, S., and Tuya, F. (2013). Patrones de abundancia y talla de Holothuria sanctori, Holothuria mammata y Holothuria arguinensis (Echinodermata: Holoturoidea) en la isla de Gran Canaria, Atlántico oriental. Rev. Biol. Mar. Oceanogr. 48, 273–284. doi: 10.4067/S0718-19572013000200007
Navarro, P. G., García-Sanz, S., and Tuya, F. (2014). Contrasting displacement of the sea cucumber Holothuria arguinensis between adjacent nearshore habitats. J. Exp. Mar. Biol. Ecol. 453, 123–130. doi: 10.1016/j.jembe.2014.01.008
Pedrotti, M., and Fenaux, L. (1992). Dispersal of echinoderm larvae in a geographical area marked by upwelling (Ligurian Sea, NW Mediterranean). Mar. Ecol. Prog. Ser. 87, 217–227. doi: 10.3354/meps087217
Purcell, S., Conand, C., Uthicke, S., and Byrne, M. (2016). Ecological roles of exploited sea cucumbers. Oceanogr. Mar. Biol. Annu. Rev. 54, 367–386. doi: 10.1201/9781315368597-8
Purcell, S., Samyn, Y., and Conand, C. (2012). Commercially Important Sea Cucumbers of the World. Rome: FAO Food and Agriculture Organization.
Purcell, S. W. (2010). Diel burying by the tropical sea cucumber Holothuria scabra: effects of environmental stimuli, handling and ontogeny. Mar. Biol. 157, 663–671. doi: 10.1007/s00227-009-1351-6
Purcell, S. W., and Kirby, D. S. (2006). Restocking the sea cucumber Holothuria scabra: sizing no-take zones through individual-based movement modelling. Fish. Res. 80, 53–61. doi: 10.1016/j.fishres.2006.03.020
Purcell, S. W., Mercier, A., Conand, C., Hamel, J.-F., Toral-Granda, M. V., Lovatelli, A., et al. (2013). Sea cucumber fisheries: global analysis of stocks, management measures and drivers of overfishing. Fish Fish. 14, 34–59. doi: 10.1111/j.1467-2979.2011.00443.x
R Core Team (2020). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Schneider, K., Silverman, J., Kravitz, B., Rivlin, T., Schneider-Mor, A., Barbosa, S., et al. (2013). Inorganic carbon turnover caused by digestion of carbonate sands and metabolic activity of holothurians. Estuar. Coast. Shelf Sci. 133, 217–223. doi: 10.1016/j.ecss.2013.08.029
Schneider, K., Silverman, J., Woolsey, E., Eriksson, H., Byrne, M., and Caldeira, K. (2011). Potential influence of sea cucumbers on coral reef CaCO 3 budget: a case study at One Tree Reef. J. Geophys. Res. 116:G04032. doi: 10.1029/2011JG001755
Schwerdtner Máñez, K., and Ferse, S. C. A. (2010). The history of Makassan Trepang fishing and trade. PLoS One 5:e11346. doi: 10.1371/journal.pone.0011346
Shears, N. T., and Babcock, R. C. (2007). Quantitative Description Mainland New Zealand’s Shallow Subtidal Reef Communities. New Zealand. Available online at: https://www.doc.govt.nz/globalassets/documents/science-and-technical/sfc280entire.pdf
Shepard, F. P. (1954). Nomenclature based on sand-silt-clay ratios. SEPM J. Sediment. Res. 24, 151–158. doi: 10.1306/d4269774-2b26-11d7-8648000102c1865d
Shiell, G. (2004). Field observations of juvenile sea cucumbers. SPC Beche-de-mer Inf. Bull. 20, 6–11.
Shiell, G. R. (2007). Spatial distribution and temporal shifts in the biology of the commercial sea cucumber Holothuria whitmaei [Echinodermata: Holothuroidea], Ningaloo Reef, Western Australia. SPC Beche-de-mer Inf. Bull. 26, 35–36.
Shiell, G. R., and Knott, B. (2008). Diurnal observations of sheltering behaviour in the coral reef sea cucumber Holothuria whitmaei. Fish. Res. 91, 112–117. doi: 10.1016/j.fishres.2007.12.010
Siegenthaler, A., Cánovas, F., and González-Wangüemert, M. (2015). Spatial distribution patterns and movements of Holothuria arguinensis in the Ria Formosa (Portugal). J. Sea Res. 102, 33–40. doi: 10.1016/j.seares.2015.04.003
Siegenthaler, A., Canovas, F., and Wanguemert, M. G. (2017). Outlanders in an Unusual Habitat: Holothuria mammata (Grube, 1840) Behaviour On Seagrass Meadows from Ria Formosa (S Portugal). Turkish J. Fish. Aquat. Sci. 17, 1031–1038. doi: 10.4194/1303-2712-v17_5_19
Simões, T., Azevedo, E., Silva, F., Sousa, J., Brito, A., Pombo, A., et al. (2018). Ecological traits of sea cucumbers with commercial relevance from the north-eastern Atlantic coast. Front. Mar. Sci. 5:147. doi: 10.3389/conf.fmars.2018.06.00147
Simunovic, A., Piccinetti, C., Bartulovic, M., and Grubelic, I. (2000). Distribution and abundance of the species Holothuria tubulosa GMELIN, 1788 and Holothurla forskall DELLE CHIAJE, 1823 (Holothurla Echinodermata) in the Adriatic Sea. Acta Adriat. 41, 3–16.
Sloan, N., and von Bodungen, B. (1980). Distribution and Feeding of the Sea Cucumber Isostichopus badionotus in Relation to Shelter and Sediment Criteria of the Bermuda Platform. Mar. Ecol. Prog. Ser. 2, 257–264. doi: 10.3354/meps002257
Sloan, N. A. (1979). Microhabitat and resource utilization in cryptic rocky intertidal echinoderms at Aldabra Atoll, Seychelles. Mar. Biol. 54, 269–279. doi: 10.1007/BF00395789
So, J. J., Hamel, J.-F., and Mercier, A. (2010). Habitat utilisation, growth and predation of Cucumaria frondosa: implications for an emerging sea cucumber fishery. Fish. Manag. Ecol. 17, 473–484. doi: 10.1111/j.1365-2400.2010.00747.x
Solan, M. (2004). Extinction and ecosystem function in the marine benthos. Science 306, 1177–1180. doi: 10.1126/science.1103960
Soliman, T., Reimer, J. D., Kawamura, I., van der Meij, S. E. T., Reijnen, B. T., and Paulay, G. (2019). Description of the juvenile form of the sea cucumber Thelenota anax H. L. Clark, 1921. Mar. Biodivers. 49, 547–554. doi: 10.1007/s12526-017-0820-2
Sonnenholzner, J. (2003). Seasonal variation in the food composition of Holothuria theeli (Holothuroidea: Aspidochirotida) with observations on density and distribution patterns at the central coast of Ecuador. Bull. Mar. Sci. 73, 527–543.
Sousa, J., Félix, P. M., Brito, A. C., Venâncio, E., Azevedo e Silva, F., Simões, T., et al. (2020). The effects of stocking density on physiological traits in Holothuria forskali broodstock. Aquac. Res. 52, 1–13. doi: 10.1111/are.15007
Tortonese, E. (1965). Fauna d’Italia: Vol. VI, Echinodermata. Calderini. doi: 10.1002/iroh.19680530111
Tuya, F., Hernández, J. C., and Clemente, S. (2006). Is there a link between the type of habitat and the patterns of abundance of holothurians in shallow rocky reefs? Hydrobiologia 571, 191–199. doi: 10.1007/s10750-006-0240-y
Tweedie, M. C. K. (1984). “An index which distinguishes between some important exponential families,” in Statistics: Applications and New Directions, eds J. K. Ghosh and J. Roy (Calcutta: Indian Statistical Institute), 579–604.
Underwood, G. J. C., and Barnett, M. (2006). “What determines species composition in microphytobenthic biofilms?,” in Functioning of Microphytobenthos in Estuaries. Microphytobenthos Symposium, ed. J. Kromkamp (Amsterdam: Royal Netherlands Academy of Arts and Sciences), 121–138.
Uthicke, S. (1999). Sediment bioturbation and impact of feeding activity of Holothuria (Halodeima) atra and Stichopus chloronotus, two sediment feeding holothurians, at Lizard Island, great barrier reef. Bull. Mar. Sci. 64, 129–141.
Venâncio, E., Félix, P. M., Brito, A. C., Sousa, J., Azevedo e Silva, F., Simões, T., et al. (2021). Do broodstock diets influence viability and larval development of Holothuria mammata? Aquaculture 536:736431. doi: 10.1016/j.aquaculture.2021.736431
Viyakarn, V., Chavanich, S., Heery, E., and Raksasab, C. (2020). Distribution of sea cucumbers, Holothuria atra, on reefs in the upper Gulf of Thailand and the effect of their population densities on sediment microalgal productivity. Estuar. Coast. Shelf Sci. 235:106514. doi: 10.1016/j.ecss.2019.106514
Walag, A. M. P., and Canencia, M. O. P. (2016). Physico-chemical parameters and macrobenthic invertebrates of the intertidal zone of Gusa, Cagayan de Oro City, Philippines. AES Bioflux 8, 71–82. doi: 10.6084/M9.FIGSHARE.2057994.V1
Wolkenhauer, S.-M., Uthicke, S., Burridge, C., Skewes, T., and Pitcher, R. (2010). The ecological role of Holothuria scabra (Echinodermata: Holothuroidea) within subtropical seagrass beds. J. Mar. Biol. Assoc. U.K. 90, 215–223. doi: 10.1017/S0025315409990518
Wooster, W. S., Bakun, A., and McLain, D. R. (1976). Seasonal upwelling cycle along the eastern boundary of the North Atlantic. J. Mar. Res. 34, 131–141.
Keywords: Holothuroidea, density, size-class distribution, habitat preferences, environmental parameters, rocky reef, Portugal
Citation: Félix PM, Pombo A, Azevedo e Silva F, Simões T, Marques TA, Melo R, Rocha C, Sousa J, Venâncio E, Costa JL and Brito AC (2021) Modelling the Distribution of a Commercial NE-Atlantic Sea Cucumber, Holothuria mammata: Demographic and Abundance Spatio-Temporal Patterns. Front. Mar. Sci. 8:675330. doi: 10.3389/fmars.2021.675330
Received: 02 March 2021; Accepted: 02 July 2021;
Published: 04 August 2021.
Edited by:
Donata Melaku Canu, National Institute of Oceanography and Experimental Geophysics (OGS), ItalyReviewed by:
Vinko Bandelj, National Institute of Oceanography and Applied Geophysics, ItalyNina Bednarsek, Southern California Coastal Water Research Project, United States
Copyright © 2021 Félix, Pombo, Azevedo e Silva, Simões, Marques, Melo, Rocha, Sousa, Venâncio, Costa and Brito. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Pedro M. Félix, pmfelix@fc.ul.pt