Vol. 13: 301–310, 2021 ENVIRONMENT INTERACTIONS Published August 12 https://doi.org/10.3354/aei00409 Aquacult Environ Interact

OPENPEN ACCESSCCESS Ecological co-benefits from sea cucumber farming: Holothuria scabra increases growth rate of

Jessica Arnull1,*, A. Meriwether W. Wilson1, Kitty Brayne3, Kyle Dexter1,2, Angelo G. Donah4, Charlotte L. A. Gough3, Timothy Klückow4, Bryne Ngwenya1, Alexander Tudhope1

1School of GeoSciences, University of Edinburgh, James Hutton Road, Edinburgh EH9 3FE, UK 2Royal Botanic Garden Edinburgh, 20a Inverleith Row, Edinburgh EH3 5LR, UK 3Blue Ventures, The Old Library, Trinity Road, Bristol BS2 0NW, UK 4Blue Ventures, Villa Huguette, Lot II U 86, Cité Planton, Ampahibe, 101 Antananarivo,

ABSTRACT: Sea cucumber aquaculture is increasing in extent and importance throughout the Indo-Pacific region, supplying a luxury seafood market in Asia. In this context, the grow-out of hatchery-bred juveniles in community-farmed pens is proving to be a viable model, providing increased income security and alternative livelihood options to resource-limited communities. Here, we report a study of the impacts of such sea cucumber farming on the growth of seagrass (a favourable habitat for the animals) at a village-scale aquaculture site in southwest Madagascar. Using experiments, we found that the presence of the hatchery-bred sea cucumber Holothuria scabra (sandfish), at stocking densities of 300 g m−2 (similar to the density used in the farmed pens, but relatively high for natural populations), resulted in a large (~30%), statistically significant increase in the leaf extension rate of the locally dominant seagrass species Thalassia hemprichii. However, the other dominant seagrass species, Cymodocea serrulata, did not significantly change its leaf extension rate in the presence of H. scabra. Since seagrass is a globally important coastal habitat, supporting high biodiversity, carbon sequestration, shoreline stability and nursery grounds for commercial and small-scale fisheries, the positive effect of H. scabra farming on the growth rate of at least one dominant seagrass species implies potential important ecological co- benefits. These co-benefits of H. scabra farming are likely to be relevant across the tropical Indo- Pacific coastlines, where this species is cultured.

KEY WORDS: Holothuria scabra · Sea cucumber aquaculture · Seagrass · Thalassia hemprichii · Cymodocea serrulata · Madagascar · LMMA · Ecosystem services

1. INTRODUCTION 1833 (‘sandfish’ or known in local dialect in SW Madagascar as ‘zanga fotsy’), which inhabit shallow Holothuroids (sea cucumbers) are benthic organ- coastal waters such as seagrass beds across the Indo- isms belonging to Phylum Echinodermata. They Pacific (Mercier et al. 2000, Hamel et al. 2001, Purcell have been harvested for centuries to supply Asian et al. 2012a,b). Rising demand in the 1980s led to the seafood markets, primarily as a luxury dried food depletion and over-extraction of sea cucumber fish- item known as bêche-de-mer or trepang (Purcell eries and natural stocks (Mercier et al. 1999). In some 2014, Juinio-Meñez et al. 2017). One of the most cases, exploitation became so severe that fisheries valuable of these species is Holothuria scabra Jaeger, governance and regulatory measures alone have

© The authors 2021. Open Access under Creative Commons by *Corresponding author: [email protected] Attribution Licence. Use, distribution and reproduction are un- restricted. Authors and original publication must be credited. Publisher: Inter-Research · www.int-res.com 302 Aquacult Environ Interact 13: 301–310, 2021

been incapable of restoring populations (Mercier et 2000, Hamel et al. 2001, Wolkenhauer et al. 2010, al. 1999, Anderson et al. 2011, Friedman et al. 2011, Purcell 2012, Hair et al. 2016, Ceccarelli et al. 2018). Purcell et al. 2012a, Conand 2018), leading to the Whilst seagrass habitats provide benefits for H. development of the sea cucumber aquaculture sector scabra, seagrass may also benefit from the presence (Eriksson et al. 2012, Purcell et al. 2012b). Typically of holothurians (Wolkenhauer et al. 2010). Through selling for up to US$165 kg−1 (for specimens <8 cm) their diel burying cycle, species like H. scabra bury and more than US$840 kg−1 (for specimens >12 cm) into the sediment, potentially having a greater in Hong Kong or Guangzhou, H. scabra is now the impact on surface and sub-surface sediment dis- most extensively cultured tropical species of sea placement and mixing than other holothurians cucumber (Purcell 2014, Juinio-Meñez et al. 2017). (Yamanouchi 1939, 1956, Mercier et al. 1999, Purcell In the southwest region of Madagascar, sea cucum- et al. 2016). This burying has been hypothesised to ber fishing in shallow coastal lagoons is an integral increase benthic primary production by irrigating aspect of the livelihoods of the Vezo people who and oxygenating deeper sediment layers (Mercier et inhabit this region (Rasolofonirina 2007, Robinson & al. 1999, Purcell et al. 2016), allowing for greater Pascal 2009). Since 2009, the British-based conserva- belowground seagrass growth (Wolkenhauer et al. tion group (BV) has been making con- 2010, Rougier et al. 2013). This behaviour may (1) siderable pro gress working with communities to transfer or displace buried sedimentary organic mat- develop an innovative sea cucumber aquaculture ter to the sediment surface, encouraging aerobic livelihood model in the village of Tampolove, a small remineralisation and (2) release nutrients trapped in Vezo fishing village situated in the protected waters interstitial waters, thereby increasing nutrient efflux of the Velondriake locally marine managed area to both the water and sediment columns for uptake (LMMA) (Cripps & Harris 2009). Here, BV and the by seagrass (Massin 1982, Uthicke 1999, Grall & Velondriake Association have been working with the Chauvaud 2002, de Witt 2009). Similarly, the inges- village to establish a viable community-led sea tion and excretion of large volumes of sediment by cucumber aquaculture business that allows a steady holothurians has been shown to enrich dissolved monthly income through the sale of mature adults nutrients in the surrounding sediment (Webb et al. whilst reducing wild-catch fishing pressure and 1977, Conde et al. 1991, Uthicke 2001, Costa et al. potentially supporting the recovery of wild sea 2014), which may also stimulate seagrass growth cucumber populations (Gardner et al. 2020). (Hughes et al. 2004). In support of this hypothesis, a Although the last 2 decades have seen increased ef- 2 yr experiment carried out by Wolkenhauer et al. forts in sea cucumber aquaculture (e.g. Mercier et al. (2010) found that when wild H. scabra were 1999, Lavitra et al. 2010, Purcell 2010, Hair et al. 2016, excluded from experimental pens, seagrass produc- 2020, Altamirano et al. 2017, Sinsona & Juinio-Meñez tivity and biomass were diminished in some (but not 2018), research addressing the basic biology and all) experimental periods. ecology of hatchery-bred H. scabra remains limited. If this ecological mutualism does exist between H. Studies have identified substantial ecological impacts scabra and seagrass, it is likely that under certain cir- when benthic marine organisms are re moved from a cumstances a greater density of H. scabra could be particular coastal environment (Solan et al. 2004), and positively correlated with greater seagrass growth, given the importance of sea cucumbers in recycling as suggested by Rougier et al. (2013), noting that and remineralising nutrients and organic matter there is a recognised maturity-dependency on sand- through feeding, excretion and bioturbation activities fish preference for seagrass (e.g. Ceccarelli et al. (Purcell et al. 2016), it is reasonable to assume that 2018). As seagrass meadows are considered highly both the removal and the introduction of these ani- productive ecosystems (Duarte & Chiscano 1999, mals — especially at high densities — could generate Hemminga & Duarte 2000), a positive mutualism considerable ecological impacts on the surrounding could enhance wider ecosystem benefits. Seagrass coastal ecosystem (Thomson 2017, Ceccarelli et al. meadows provide a number of key ecological func- 2018, Lee et al. 2018). tions and services, including shoreline stabilisation, Of particular importance is the relationship that H. coastline protection, nutrient cycling, commercial scabra have with seagrass — a favourable habitat for and small-scale fisheries support, biodiversity the animals. Sea cucumber larvae rely on seagrass enhancement, nursery grounds and habitats for for their settling cues, and the seagrass then provides invertebrates, fish and large predators (Hemminga & the juveniles with protection from predation and a Duarte 2000, Orth et al. 2006, Heck et al. 2008, suitable substrate on which to grow (Mercier et al. Belshe et al. 2017) and livelihood opportunities for Arnull et al.: H. scabra increases growth rate of seagrass 303

many of the world’s poorest communities (Unsworth Assassins; a broad and shallow bay located ~150 km et al. 2010). As of yet, no studies have addressed the north of Toliara in southwest Madagascar (22° 13.3’ S, potentially important role of burying holothurians on 43° 16.0’ E; Fig. 1). The bay is well-flushed, with a seagrass ecosystems in the context of an aquaculture water depth ranging to ~10 m in the central area, a farm, where population densities tend to be higher tidal range of ~2.9 m during spring tides and ~1.5 m than in natural local settings. during neap tides and no major freshwater inputs. Here, we report a study of a village-scale, com - On average, this region has a hot, semi-arid climate, munity-led aquaculture site in southwest Madagas- with open water sea surface temperature (SST) typi- car conducted in natural seagrass habitat. The goals cally ranging from about 24.1°C in August to 29.1°C were to ascertain whether the presence of hatchery- in February (mean monthly SST 1982−2019, IGOSS- bred sandfish affected growth of the dominant sea- nmc; Reynolds et al. 2002). The 4 ha farm is located grass species in the study area, adjacent to a commu- ~750 m offshore from the village (Fig. 1) and operates nity H. scabra farming operation. We hypothesised as a no-take zone (with the exception of the principal that the addition of H. scabra would increase sea- sandfish predator, the blue swimming crab Thala- grass growth through increased bioturbating and mita crenata, which are removed when encoun- burying activities, improving sediment oxygenation tered). The boundaries are delineated by a grid of 41 and enhancing organic matter remineralisation, farmed pens (30 × 30 m; Fig. 2), and are used to grow thereby increasing primary productivity. out hatchery-supplied juvenile Holothuria scabra until they reach marketable size. The farm has an average water depth of 0.5 m relative to low water 2. MATERIALS AND METHODS spring tide level and sits atop a thin veneer of sedi- ment, ranging in thickness from 0.01−1.12 m. The 2.1. Study site seabed at the site is mostly covered with seagrass with some bare sandy patches and is dominated by Tampolove is a small village (population ~500 peo- poorly−very poorly sorted fine sand with an average ple) situated on the southern shore of the Bay of carbonate content of around 45−50%. In total, 7 spe-

Fig. 1. Location of Tampolove within the Bay of Assassins in southwest Madagascar 304 Aquacult Environ Interact 13: 301–310, 2021

current stock intake, the experiment did not attempt to simulate the natural variability in the size of the animals usually found in the commercially farmed pens, which range from <15 g to >500 g. Furthermore, as the experi- ment was based on biomass (g m−2) rather than ind. m−2, the number of an- imals per pen varied between and Fig. 2. (a) Tampolove site map (not to scale) with location of experimental sites within the treatments, although care indicated. M: mirador watch tower. (b) View of the aquaculture farm and was taken to ensure the size and num- mirador ber of H. scabra was roughly consistent between pens (Table S1). cies of seagrass have been identified and are listed here in order of decreasing prevalence across the site: Thalassia hemprichii, Cymodocea serrulata, Ha - 2.3. Pen characterisation lodule uninervis, Thalassodemdron ciliatum, Halo - phila ovalis, H. decipiens and C. rotundata. Prior to the experiment, data were collected for several descriptor variables for each experimental pen through visual surveys. Seagrass cover (%) and 2.2. Experimental design species composition (%) were calculated using stan- dardized estimates from Seagrass Watch (McKenzie The experiment was conducted over 2 tidal cycles et al. 2001) to visually estimate seagrass cover on a during July and August 2019 and was designed to scale of 1−100%. Sediment thickness (m) and sea- identify any differences in seagrass growth rates asso- and rockbed elevation were calculated by measur- ciated with the presence of H. scabra. Due to the na- ing sediment and water depth and then calibrating ture of the site and experiment, all sampling occurred against a local tidal datum (spring low tide level on foot during spring low tides. The experiment was observed during the experiment). The number and undertaken in 4 m2 pens that were constructed (see type of holothurian and thalassinidean shrimp sedi- Text S1 in the Supplement at www. int-res. com/ ment mounds (species unknown and both of which articles/ suppl/ q013 p301 _supp. pdf for detail) in sea- may have an impact on sediment oxygenation) grass beds at 4 (unfarmed) sites (3 pens site−1) on the present within each pen was recorded. Sediment periphery of the farm (n = 12) in areas void of wild sea samples were collected from the top 5 cm of sedi- cucumber populations. Site location (and pen posi- ment prior to the experiment for subsequent analy- tions therein) were chosen to ensure comparable sea- sis of sediment grain size and sediment organic car- grass coverage, and the position of each treatment bon content (Corg %). pen was haphazardly chosen within each site. All 4 sites (identified A−D; Fig. 2) were densely covered with mixed beds of either T. hemprichii and H. unin- 2.4. Seagrass growth ervis (dominating Sites A and B), or C. serrulata and H. uninervis (dominating Sites C and D). At each site, Seagrass growth was measured by determining 3 different treatments were applied: control pens, the vertical growth in seagrass blades over the void of H. scabra; medium-density pens, where pens course of 3 spring tides. Growth was measured in T. were stocked with H. scabra at a density of 150 g m−2; hemprichii blades at Sites A and B, and in C. serru- and high-density pens, where pens were stocked with lata blades at Sites C and D (even where T. H. scabra at a density of 300 g m−2. The high-density hemprichii dominated cover in pens C1, C2 and D1; treatment was based on the farm’s optimal stocking see Table 1). Following the ‘hole-punch’ method density, which is determined by the critical biomass detailed by Zieman (1974), seagrass blades were value of the sediment (Lavitra et al. 2010, Plotieau et stapled at the base of the leaf, left to grow and har- al. 2013). Animals used to stock the experimental vested after 14 d. The distance from the base of the pens weighed between 30 and 150 g and were sup- leaf to the staple was measured wet (cm), allowing plied by the hatchery used to stock the commercial for leaf extension rates to be expressed as cm pens. Whilst these sandfish were representative of the shoot−1 d−1. Arnull et al.: H. scabra increases growth rate of seagrass 305

To account for any between- or within-site varia- checked for normality and homoscedasticity with tion in environmental conditions, seagrass extension respect to fitted values and explanatory variables. To rates were first calculated without the presence of H. evaluate the importance of holothurian density, we scabra over the first tidal cycle (pre-stocking). In compared the values of Akaike’s information crite- 2 each experimental pen, a small area of 20 cm was rion corrected for small sample size (AICC) of a model haphazardly marked (ensuring sufficient coverage of with holothurian density as a fixed effect versus a the corresponding seagrass species at the respective null model without. A difference in AICC values sites) using a quadrat and small metal rods. Within greater than 2 can be taken as evidence that there is these 20 cm2 areas, a minimum of 25 seagrass blades a statistically detectable effect (Burnham & Anderson were marked individually with a staple and their 2004). We conducted additional tests on the pairs of number recorded, although due to the heterogeneity sites with similar dominating seagrass species (i.e. in seagrass cover, this number varied between pens; just Sites A and B, and just Sites C and D). Finally, to 14 d later, all stapled blades were harvested and assess if any significant effects of holothurian density stored in water for subsequent measurement. Pens were present at medium densities or just at high den- were then stocked with the sea cucumbers (at the sities, we constructed similar models but grouped the respective densities) and the process was repeated control with medium-density plots and alternatively for another 20 cm2 area within the pens over the next grouped the medium-density and high-density plots. tidal cycle (post-stocking). The dynamic coastal set- The AICC values of these models were compared to ting meant that the blades naturally suffered varying the model with 3 separate density categories. Signif- degrees of loss (e.g. through grazing or disturbance), icance was accepted at p ≤ 0.05, and analyses were meaning not all blades that were stapled survived conducted in R v.4.0.2. the experiment. The average number of stapled blades retrieved at the end of the experiment was 25 but ranged from 16−32 among pens. 3. RESULTS

3.1. Pen characterisation 2.5. Elemental and sediment grain size analysis Environmental conditions were similar among

Sediment samples were analysed for total Corg (%) pens within the same site, though there were some content using an organic elemental analyser (Carlo noteworthy between-site differences. Site A was Erba NA 2500, Eager 300 software). Sediment grain most anomalous, with considerably lower seagrass size was measured using a Beckman-Coulter LS 230 cover and a thinner sediment layer than the other 3 laser particle size analyser equipped with a fluid sites. Sites A and B had greater average water depth module (see Text S2). (lower seabed elevations) than Sites C and D. Sea- grass species composition within Sites C and D were also considerably more varied than those from Sites 2.6. Statistical analyses A and B (Table 1), with pens C1, C2 and D1 domi- nated by Thalassia hemprichii. Despite these differ- All descriptive statistics are presented as means ± ences, linear regression analyses indicated no signif- SE unless stated otherwise. Given non-normality of icant relationships between pen characteristics and data and small sample sizes, non-parametric Mann- seagrass leaf extension rates. Whitney U-tests were used to identify any significant within-site and between-site differences for pre- stocking data as well as any significant differences 3.2. Seagrass growth between pre-stocking and post-stocking control data. Linear regression was used to determine rela- Due to the physiological differences in the seagrass tionships between pen characteristics and seagrass species, analyses were carried out separately for the leaf growth rates. To assess the impact of stocking 2 pairs of sites (A−B and C−D), allowing differences density, we constructed mixed effects models with in growth rates between species to be measured. the growth of individual grass blades (cm shoot −1 During pre-stocking, no statistically significant dif- d−1) as the response variable, holothurian density as a ferences were detected between pens within a single fixed effect (control, medium and high) and site as a site. When control-pen data were compared between random effect (A, B, C and D). Model residuals were the pre- and post-stocking growing phases (thus 306 Aquacult Environ Interact 13: 301–310, 2021

Table 1. Descriptor variables characterising the environmental conditions of each individual aquaculture pen. Corg: organic carbon. Th: Thalassia hemprichii; Cs: Cymodocea serrulata; Hu: Halodule uninervis; S: shrimp mound; H: holothurian mound. Seabed and rock elevation data are relative to low water spring tide of 0.35 m (17 May 2018 at 14:15 h)

Pen Seagrass Species Number Total Corg Sediment grain Sediment Seabed Rock cover (%) composition of (%) size (median depth (m; mean elevation elevation (%) mounds value; μm) ± SE, n = 3) (m) (m)

A1 78 94:6 (Th:Hu)1 (S) 0.5 131 0.17 ± 0.003 0.08 −0.05 A2 75 95:5 (Th:Hu)1 (H) 0.6 155 0.14 ± 0.003 A3 65 92:8 (Th:Hu) 0 0.7 144 0.11 ± 0.002 B1 95 95:5 (Th:Hu)2 (S) 0.7 132 0.35 ± 0.006 0.05 −0.14 B2 94 90:10 (Th:Hu) 0 0.8 125 0.28 ± 0.004 B3 94 92:8 (Th:Hu)1 (H) 0.7 124 0.31 ± 0.002 C1 93 85:10:5 (Th:Cs:Hu)2 (S) 0.9 131 0.30 ± 0.003 0.25 −0.10 C2 80 70:20:10 (Th:Cs:Hu)1 (S) 1.0 113 0.23 ± 0.003 C3 95 35:50:15 (Th:Cs:Hu)1 (H) 0.9 121 0.20 ± 0.003 D1 92 65:5:30 (Th:Cs:Hu) 0 0.8 131 0.27 ± 0.002 0.23 0.01 D2 75 10:65:25 (Th:Cs:Hu)2 (S) 0.7 115 0.29 ± 0.003 D3 95 20:75:5 (Th:Cs:Hu) 0 0.8 106 0.22 ± 0.003 accounting for external factors that may have influ- post-stocking phase to the addition of sea cucum- enced growth rates between the 2 sampling periods), bers. During post-stocking, in both the medium- and there were no statistically significant differences in high-density pens, leaf extension rates increased for seagrass growth rates between the 2 growing phases T. hemprichii but not C. serrulata (Fig. 3). The high- for either species (T. hemprichii: U = 1467, p = 0.670, density pens displayed the highest growth rates for n = 111; Cymodocea serrulata: U = 924.5, p = 0.295, T. hemprichii (0.307 ± 0.023 cm shoot−1 d−1), which n = 92). These observations improve confidence in was ~33% faster on average than that observed in attributing any observed between-treatment differ- the controls (0.231 ± 0.021 cm shoot−1 d−1). The addi- ences in seagrass growth rate within sites during the tion of Holothuria scabra, irrespective of density, pro-

Fig. 3. Mean (±SE) seagrass leaf extension rates for the 2 most dominant seagrass species. Stocking densities (g Holothuria scabra m−2) used in the 3 treatments are shown below each bar for the 2 growing phases: pre-stocking and post-stocking. Each bar represents the mean growth rate of all individual seagrass blades from Sites A and B for Thalassia hemprichii and from Sites C and D for Cymodocea serrulata for each treatment (number of observations are shown underneath error bars) Arnull et al.: H. scabra increases growth rate of seagrass 307

duced a negligible effect on growth rates for C. ser- which may have led to different growth rate re- rulata, with the high-density pens producing some- sponses that were not fully captured within this what (but not statistically significantly) lower growth study. As there are many different aspects to plant rates than both the control and medium pens. expansion and maturation, it is difficult to capture a Our mixed effects statistical models indicate that holistic picture of whole plant growth with just one holothurian stocking density did have a statistically method. Future experiments would benefit from detectable effect on the growth of T. hemprichii measuring belowground biomass growth through Δ (Table 2; AICC of 2.5 between null model and full rhizome marking and plastochrone interval. Ranges model), but not C. serrulata. This effect seems to be in seagrass growth observed through the experimen- driven by higher seagrass growth in the high-density tal period correspond reasonably well with maxi- plots since grouping the medium-density plots with mum growth rates found in other studies for T. the control plots into a single category improved the hemprichii (e.g. Ogden & Ogden 1982, Duarte 1991, model (i.e. lowered the AICC values; Table 2). Con- Erftemeijer & Herman 1994, Lin & Shao 1998), but for versely, when we conducted analyses across all 4 Cymodocea serrulata, our growth ranges were more groups of pens (i.e. pens dominated by T. hemprichii variable compared to other studies (e.g. Udy & Den- and those dominated other seagrass species), we did nison 1997, Kamal et al. 1999, de Boer 2000, Yamada not detect any statistical effect of holothurian density et al. 2018).

(null model had a lower AICC value by 0.7 units). To Bioturbators are a natural component of seagrass emphasise this result, when the analysis was re- ecosystems, providing essential services within sea- stricted to just C. serrulata growth, the null model grass meadows, without which they would struggle outperformed the full model by 4.1 AICC units. to exist (Kristensen et al. 2012, Thomson 2017). Whilst the role of sea cucumbers as bioturbators is less well understood, these results indicate that some 4. DISCUSSION species like H. scabra could have greater ecosystem value than previously thought through their alter- The experimental data indicate that Holothuria ation of sediment biology, geochemistry and physical scabra have a net overall positive effect on seagrass structure. For , sustained functioning of leaf growth but that this effect is conditional on 2 roots and rhizomes depends on an adequate supply considerations. Firstly, the data suggest that there is of oxygen (Duarte et al. 1998). Therefore, it is likely some threshold effect of H. scabra density on sea- that burying activities carried out by H. scabra desta- grass growth as only the high-density pens experi- bilise the stratified sediment causing increased aera- enced a large (~30%) and statistically significant tion, thus increasing sediment oxygen supply and increase in growth rates of Thalassia hemprichii. This contributing to plant growth. These actions may also may imply that only at certain density will H. scabra potentially release organic material and nutrients make a notable difference on the local environment. trapped in sediment porewaters back into the water Secondly, it suggests that this relationship may only column (Massin 1982, Uthicke 1999), thereby provid- hold true for certain species of seagrass. This may be ing additional nutrients for seagrass. Likewise, the due to differing functional traits of the 2 species, ingestion of large volumes of organic compounds contained within subsurface sediments and subse- quent excretion of this material on the seabed results Table 2. Model performance for mixed effect models to eval- in the recycling of nutrients (particularly ammonium; uate the effect of holothurian stocking density on seagrass growth for Sites A and B only (Thalassia hemprichii). Dif- cf. Uthicke 2001 for Holothuria atra) back into the ferences in Akaike’s information criterion corrected for water column and sediments (Webb et al. 1977, Δ small sample size ( AICC) values are given relative to the Massin 1982, Wainright 1990, Conde et al. 1991, Δ best-performing model, which itself has a AICC value of 0 Uthicke 2001), thereby stimulating seagrass growth (p < 0.05). C: control (no Holothuria scabra); M: medium stocking density; H: high stocking density (Hughes et al. 2004, de Witt 2009, Costa et al. 2014). To more fully constrain and understand the impact of farmed H. scabra on seagrass growth and produc- Model AIC ΔAIC C C tivity, future studies would benefit from considering Null model (only random effects) −1128.0 3.7 3 additional factors that were beyond the scope of the Model with C, M and H separate −1130.5 1.2 current study. Firstly, estimation of changes in root Model grouping M and H vs. C −1129.9 1.8 and rhizome biomass and production to give a more Model grouping C and M vs. H −1131.7 0 complete picture of the impacts of total seagrass pro- 308 Aquacult Environ Interact 13: 301–310, 2021

ductivity. Secondly, undertaking multi-year experi- Anderson SC, Flemming JM, Watson R, Lotze HK (2011) ments to assess long-term impacts. And, finally, Serial exploitation of global sea cucumber fisheries. Fish Fish 12: 317−339 greater replication (including more experimental Belshe EF, Mateo MA, Gillis L, Zimmer M, Teichberg M pens) to increase statistical confidence and help rule (2017) Muddy waters: Unintentional consequences of out the possibility of ‘non-demonic’ intrusion (sensu research obscure our understanding of Hurlbert 1984) influencing the interpretation of organic carbon dynamics in seagrass ecosystems. Front Mar Sci 4: 125 observed differences in growth rates. Burnham KP, Anderson DR (2004) Multimodel inference: understanding AIC and BIC in model selection. Sociol Methods Res 33: 261−304 5. CONCLUSIONS Ceccarelli DM, Logan M, Purcell SW (2018) Analysis of opti- mal habitat for captive release of the sea cucumber Holothuria scabra. Mar Ecol Prog Ser 588: 85−100 This study found that Holothuria scabra, stocked at Conand C (2018) Tropical sea cucumber fisheries: changes the relatively high population density used in aqua- during the last decade. Mar Pollut Bull 133:590−594 culture pens, significantly increased the growth rate Conde JE, Diaz H, Sambrano A (1991) Disintegration of of a common seagrass species. Since seagrass eco- holothurian fecal pellets in beds of the seagrass Thalas- sia testudinum. J Coast Res 7: 853−862 systems are considered to be amongst the most pro- Costa V, Mazzola A, Vizzini S (2014) Holothuria tubulosa ductive and diverse on earth (Belshe et al. 2017, Gmelin 1791 (Holothuroidea, Echinodermata) enhances Unsworth et al. 2010) but are also threatened glob- organic matter recycling in Posidonia oceanica mead- ally by human activities (Wylie et al. 2016), the po- ows. J Exp Mar Biol Ecol 461: 226−232 sitive impact of H. scabra aquaculture on growth Cripps G, Harris A (2009) Community creation and manage- ment of the Velondriake marine protected area. Blue rates raises the possibility of important ecological co- Ventures Conservation Report. Blue Ventures Conserva- benefits. The sea cucumber farm at Tampolove oper- tion, London ates under a strict ‘no-take-zone’ policy and can be de Boer WF (2000) Biomass dynamics of seagrasses and the viewed as ‘non-invasive’ (i.e. the farm requires no role of and seagrass vegetation as different nutrient sources for an intertidal ecosystem. Aquat Bot inputs such as feed and uses a native and endan- 66: 225−239 gered species that were previously highly abundant de Witt TH (2009) The effects of bioturbation and bioirriga- in this area). Therefore, not only is the farm offering tion on seagrasses. In: Nelson WG (ed) Seagrasses and an alternative livelihood, but the high density of protective criteria: a review and assessment of research status. US Environmental Protection Agency, Newport, farmed H. scabra may simultaneously protect, con- OR, p 176−205 serve and even ‘boost’ the seagrass habitat, with Duarte CM (1991) Allometric scaling of seagrass form and potential knock-on effects to recovery of local fish productivity. Mar Ecol Prog Ser 77: 289−300 stocks and wild sandfish populations. This finding Duarte CM, Chiscano CL (1999) Seagrass biomass and pro- has relevance to sandfish farming across many duction: a reassessment. Aquat Bot 65: 159−174 Duarte CM, Merino M, Agawin NSR, Uri J and others (1998) regions of the tropical Indo-Pacific. Root production and belowground seagrass biomass. Mar Ecol Prog Ser 171: 97−108 Erftemeijer PL, Herman PM (1994) Seasonal changes in Acknowledgements. This research was conducted as part of environmental variables, biomass, production and nutri- a joint research project between the University of Edinburgh ent contents in two contrasting tropical intertidal sea- and Blue Ventures and was co-funded by the Prince Albert II grass beds in South Sulawesi, Indonesia. Oecologia 99: of Monaco Foundation, Project No. 2505: Community-Based 45−59 Aquaculture as a Catalyst for Locally Managed Marine Eriksson H, Robinson G, Slater M, Troell M (2012) Sea Areas: Developing a Scalable Framework for Economic and Environmental Sustainability. We are very grateful to the cucumber aquaculture in the western Indian Ocean: Blue Ventures staff, whose guidance in the field was invalu- challenges for sustainable livelihood and stock improve- able (in particular to Godard Diome, Feliny Bino, Alexandre ment. Ambio 41: 109−121 Razfimandimby and Marceline ‘Kokoly’ Narindra) and also Friedman K, Eriksson H, Tardy E, Pakoa K (2011) Manage- to Amber Carter, Kathryn Taylor, Mario Ray and Seb ment of sea cucumber stocks: patterns of vulnerability Hennige for field assistance and Gavin Sim for assistance and recovery of sea cucumber stocks impacted by fish- with sediment analyses. We are also greatly appreciative to ing. Fish Fish 12: 75−93 the reviewers, whose constructive comments and suggestions Gardner CJ, Cripps G, Day LP, Dewar K and others (2020) A helped us significantly improve the manuscript. decade and a half of learning from Madagascar’s first locally managed marine area. Conserv Sci Prac 2: e298 Grall J, Chauvaud L (2002) Marine eutrophication and ben- LITERATURE CITED thos: the need for new approaches and concepts. Glob Change Biol 8: 813−830 Altamirano JP, Recente CP, Rodriguez JCJ (2017) Substrate Hair C, Mills DJ, McIntyre R, Southgate PC (2016) Optimis- preference for burying and feeding of sandfish ing methods for community-based sea cucumber ranch- Holothuria scabra juveniles. Fish Res 186: 514−523 ing: experimental releases of cultured juvenile Holo - Arnull et al.: H. scabra increases growth rate of seagrass 309

thuria scabra into seagrass meadows in Papua New Ogden JC, Ogden NB (1982) A preliminary study of two Guinea. Aquacult Rep 3:198−208 representative seagrass communities in Palau, Western Hair C, Foale S, Daniels N, Minimulu P, Aini J, Southgate Caroline Islands (Micronesia). Aquat Bot 12: 229−244 PC (2020) Social and economic challenges to commu- Orth RJ, Carruthers TJ, Dennison WC, Duarte CM and oth- nity-based sea cucumber mariculture development in ers (2006) A global crisis for seagrass ecosystems. Bio- New Ireland Province, Papua New Guinea. Mar Policy science 56: 987−996 117: 103940 Plotieau T, Baele JM, Vaucher R, Hasler A, Kounad D, Eeck- Hamel JF, Conand C, Pawson DL, Mercier A (2001) The sea haut I (2013) Analysis of the impact of Holothuria scabra cucumber Holothuria scabra (Holothuroidea: Echinoder- intensive farming on sediment. Cah Biol Mar 54: 703−711 mata): its biology and exploitation as beche-de-mer. Adv Purcell SW (2010) Diel burying by the tropical sea cucumber Mar Biol 41: 129−223 Holothuria scabra: effects of environmental stimuli, han- Heck KL Jr, Carruthers TJB, Duarte CM, Hughes AR, dling and ontogeny. Mar Biol 157:663−671 Kendrick G, Orth RJ, Williams SW (2008) Trophic trans- Purcell SW (2012) Principles and science of stocking marine fers from seagrass meadows subsidize diverse marine areas with sea cucumbers. In: Hair CA, Pickering TD, and terrestrial consumers. Ecosystems 11:1198−1210 Mills DJ (eds) Asia−Pacific tropical sea cucumber aqua- Hemminga MA, Duarte CM (2000) Light, carbon and nutri- culture. ACIAR Proceedings No. 136. Australian Centre ents. In: Hemminga MA, Duarte CM (eds) Seagrass ecol- for International Agricultural Research, Canberra, ogy. Cambridge University Press, Cambridge, p 99−145 p 92−103 Hughes AR, Bando KJ, Rodriguez LF, Williams SL (2004) Purcell SW (2014) Value, market preferences and trade of Relative effects of grazers and nutrients on seagrasses: a bêche-de-mer from Pacific Island sea cucumbers. PLOS meta-analysis approach. Mar Ecol Prog Ser 282: 87−99 ONE 9: e95075 Hurlbert SH (1984) Pseudoreplication and the design of eco- Purcell SW, Hair CA, Mills DJ (2012a) Sea cucumber cul- logical field experiments. Ecol Monogr 54: 187−211 ture, farming and sea ranching in the tropics: pro - Juinio-Meñez MA, Tech ED, Ticao IP, Gorospe JR, Edul- gress, problems and opportunities. Aquaculture 368- lantes CMA, Rioja RAV (2017) Adaptive and integrated 369: 68−81 culture production systems for the tropical sea cucumber Purcell SW, Samyn Y, Conand C (2012b) Commercially Holothuria scabra. Fish Res 186: 502−513 important sea cucumbers of the world. FAO Species Cat- Kamal A, Misri K, Japar S, Hishamuddin O, Hidir H (1999) alogue for Fishery Purposes No. 6. FAO, Rome Leaf growth, production and ecological aspects of toothed Purcell SW, Conand C, Uthicke S, Byrne M (2016) Ecologi- seagrass Cymodocea serrulata (R. Br.) Aschers. et Magnus cal roles of exploited sea cucumbers. Oceanogr Mar Biol at Port Dickson, Negeri Sembilan, Malaysia. In: Interna- Annu Rev 54: 367−386 tional Conference on the International Oceanographic Rasolofonirina R (2007) Sea cucumbers in Madagascar. In: Data and Information Exchange in the Western Pacific Conand C, Muthiga N (eds) Commercial sea cucumbers: (IODE−WESTPAC—ICIWP ’99), Lang kawi (Malaysia), a review for the Western Indian Ocean. Western Indian 1−4 November 1999. Intergovern mental Oceanographic Ocean Marine Science Association (WIOMSA), Nairobi, Commission of UNESCO, Paris, p 173−184 p 30−38 Kristensen E, Penha-Lopes G, Delefosse M, Valdemarsen T, Reynolds RW, Rayner NA, Smith TM, Stokes DC, Wang W Quintana CO, Banta GT (2012) What is bioturbation? (2002) An improved in situ and satellite SST analysis for The need for a precise definition for fauna in aquatic sci- climate. J Clim 15: 1609−1625 ences. Mar Ecol Prog Ser 446:285−302 Robinson G, Pascal B (2009) From hatchery to community— Lavitra T, Rasolofonirina R, Eeckhaut I (2010) The effect of Madagascar’s first village-based holothurian mariculture sediment quality and stocking density on survival and programme. SPC Beche-de-mer Inf Bull 29:38−43 growth of the sea cucumber Holothuria scabra reared in Rougier A, Ateweberhan M, Harris A (2013) Strategies for nursery ponds and sea pens. West Indian Ocean J Mar improving survivorship of hatchery-reared juvenile Sci 9: 153−164 Holothuria scabra in community-managed sea cucumber Lee S, Ford AK, Mangubhai S, Wild C, Ferse SC (2018) farms. SPC Beche-de-mer Inf Bull 33: 14−22 Effects of zanga (Holothuria scabra) removal on shallow- Sinsona MJ, Juinio-Meñez MA (2018) Effects of sediment water sediments in Fiji. PeerJ 6: 47−73 enrichment with macroalgae, Sargassum spp., on the Lin HJ, Shao KT (1998) Temporal changes in the abundance behavior, growth, and survival of juvenile sandfish, and growth of intertidal Thalassia hemprichii seagrass Holothuria scabra. Aquacult Rep 12: 56−63 beds in southern Taiwan. Bot Bull Acad Sin 39:191−198 Solan M, Cardinale BJ, Downing AL, Engelhardt KAM, Massin C (1982) Effects of feeding on the environment: Ruesink JL (2004) Extinction and ecosystem function in Holothuroidea. In: Jangoux M, Lawrence JM (eds) Echi- the marine benthos. Science 306: 1177−1180 noderm nutrition. AA Balkema, Rotterdam, p 493−497 Thomson ACG (2017) The role of bioturbators in seagrass McKenzie LJ, Campbell SJ, Roder CA (2001) Seagrass- blue carbon dynamics. PhD dissertation, University of Watch: manual for mapping & monitoring seagrass Technology, Sydney resources by community (citizen) volunteers. QFS, NFC, Udy JW, Dennison WC (1997) Growth and physiological Cairns responses of three seagrass species to elevated sediment Mercier A, Battaglene SC, Hamel JF (1999) Daily burrowing nutrients in Moreton Bay, Australia. J Exp Mar Biol Ecol cycle and feeding activity of juvenile sea cucumbers 217: 253−277 Holothuria scabra in response to environmental factors. J Unsworth RK, Cullen LC, Pretty JN, Smith DJ, Bell JJ (2010) Exp Mar Biol Ecol 239: 125−156 Economic and subsistence values of the standing stocks Mercier A, Battaglene SC, Hamel JF (2000) Settlement pref- of seagrass fisheries: potential benefits of no-fishing mar- erences and early migration of the tropical sea cucumber ine protected area management. Ocean Coast Manage Holothuria scabra. J Exp Mar Biol Ecol 249: 89−110 53: 218−224 310 Aquacult Environ Interact 13: 301–310, 2021

Uthicke S (1999) Sediment bioturbation and impact of feed- odermata: Holothuroidea) within subtropical seagrass ing activity of Holothuria (Halodeima) atra and Stichopus beds. J Mar Biol Assoc UK 90: 215−223 chloronotus, two sediment feeding holothurians, at Wylie L, Sutton-Grier AE, Moore A (2016) Keys to successful Lizard Island, Great Barrier Reef. Bull Mar Sci 64: blue carbon projects: lessons learned from global case 129−141 studies. Mar Policy 65: 76−84 Uthicke S (2001) Nutrient regeneration by abundant coral Yamada H, Nakamoto K, Hayakawa J, Kawamura T, Kon K, reef holothurians. J Exp Mar Biol Ecol 265:153−170 Shimabukuro H, Fukuoka K (2018) Seasonal variations Wainright SC (1990) Sediment-to-water fluxes of particulate in leaf growth of Cymodocea serrulata in subtropical material and microbes by resuspension and their contri- seagrass meadows. Fish Sci 84: 461−468 bution to the planktonic food web. Mar Ecol Prog Ser 62: Yamanouchi T (1939) Ecological and physiological studies 271−281 on the holothurians in the of Palao Islands. Webb KL, D’Elia CF, Dupaul WD (1977) Biomass and Palao Trop Biol Stn Stud 25: 603−634 nutrition flux measurements on Holothuria atra pop - Yamanouchi T (1956) The daily activity rhythms of the ulations on windward reef flats at Enewetak, Mar- holothurians in the coral reef of Palao Islands. Publ Seto shall Islands. Proc 3rd Int Coral Reef Symp, Miami 1: Mar Biol Lab 5: 347−362 410−415 Zieman JC (1974) Methods for the study of the growth and Wolkenhauer SM, Uthicke S, Burridge C, Skewes T, Pitcher production of turtle grass, Thalassia testudinum König. R (2010) The ecological role of Holothuria scabra (Echin- Aquaculture 4: 139−143

Editorial responsibility: Catriona MacLeod, Submitted: December 3, 2020 Hobart, Tasmania, Australia Accepted: June 1, 2021 Reviewed by: 2 anonymous referees Proofs received from author(s): August 6, 2021