Homes for the Holothurians: Optimizing sea cucumber management strategies through biophysical modeling

By. Mia N. Judicpa & Gabrielle A. Abrahan

As the gastronomic demand for our marine resources grows, our natural stock can barely keep up. Luckily, ecology experts work with modellers to explore different approaches in preserving the biodiversity of our oceans. In this feature, we take a look at how computational science can help us achieve sustainable trade of the famous bêche-de-mer, one of the most sought-after delicacies in the world.


The sea cucumber craze has been around since the 1300s.1 Only the affluent Chinese were able to sample its unique mouthfeel back in the day, but the appetite for it has since spread to other parts of the world. Nowadays, multiple kiosks and restaurants in East and Southeast Asia sell this curious delicacy.

Despite their chewy meat, sea cucumbers lack flavor. Chefs compensate for this by preparing them in sauces or broths that pack a punch.
Photo: Michelin Magazine

Sold fresh or dried, sea cucumbers are used to embellish a dish or added as the main protein source. Despite being a prized ingredient, they actually have little flavor and are usually slow-cooked or braised with other savory ingredients. With market prices reaching tens of thousands of pesos,2 sea cucumbers are some of the most expensive marine treats: the more eccentric the appearance, the higher the value.3

As the demand for sea cucumbers increases and their population in the intertidal zone is depleted, harvesters risk diving into deeper portions of the ocean to obtain them. Their numbers are dwindling, and sea cucumber experts around the world are scrambling to find ways to save them.

Earthworms of the sea

Sea cucumbers or holothurians (Class Holothuroidea) belong to Phylum Echinodermata, alongside commonly known but equally bizarre-looking marine invertebrates like the starfish, sea urchin, sand dollar, and sea daisy. The name comes from the Greek words echino (“spiky”) and derma (“skin”) — the spiky skin of sea cucumbers is part of the reason why they are such a prized delicacy.

Echinoderms are known for their radial symmetries.
Graphics: Roy C. Corvera

Around 1,250 species4 of sea cucumbers have been discovered worldwide. Most are found in the benthic zones of the Indian Ocean and the western Pacific region, sprawled anywhere from the seashore to seafloors as deep as six kilometers.

Due to their sedentary lives in the benthos, most sea cucumbers feed on detritus, bacteria, and biomass that settle from the waters above. They break these sediments down into even smaller particles that bacteria can consume, resulting in a process that is similar to how earthworms digest and recycle organic material.

Sea cucumbers creep around the ocean floor, feeding on organic matter and excreting inorganic nitrogen and phosphorus.
Illustration: Roy C. Corvera

Through the mechanisms of their digestion, sea cucumbers participate in the decomposition of organic matter and in nutrient cycling, both key elements to a healthy marine ecosystem. Research has shown how the absence of sandfish (Holothuria scabra) causes shallow coastal ecosystems to become stressed, as the sediments that sea cucumbers recycle typically act as buffers to sudden changes in organic matter levels.5 As these ecosystems become stressed, their function and productivity are lowered. Another study has suggested that the removal of sea cucumbers from coral reefs increases the susceptibility of some corals to sediment-borne pathogens.6

The elongated bodies of sea cucumbers hold tube feet that resemble tentacles. These are used for feeding and locomotion.

Pricy cucumbers

Sea cucumber trade is lucrative; most of the global stock is directed to China, Hong Kong, and Taiwan.1 Its dried form, the bêche-de-mer, is the most popular, with China dominating the top spots for both manufacturing and consumption of the delicacy. In the Philippines, the sea cucumber (“balat” or “balatan”) can be added as the base meat to two beloved local dishes, kinilaw and adobo.

The value of sea cucumbers is not limited to the exciting food experience; Chinese records dating back to the Ming dynasty talk about the medicinal properties of these so-called “ginseng of the sea.” With the Western market opening its doors to Eastern nutraceuticals, researchers have begun studying potential natural products and have identified several bioactive compounds with anticoagulant, antioxidant, anticancer, antimicrobial, or antiviral properties.

Sea cucumber harvest is a crucial alternative source of income for many Philippine island and coastal communities, so experts place high importance on sustainable production and consumption. The ever-increasing demand for sea cucumbers gives rise to serial exploitation, whereby previously untouched resources can be exploited if the resources in one location are depleted.7 This has caused their populations to collapse in many parts of the world, including the Philippines.

Sandfish in one of the Bolinao Marine Laboratory (BML) hatchery tubs: Some sandfish have burrowed under the sand, hiding from either light or cold temperatures, depending on their age.
Photo: Roy C. Corvera

Ranches provide a way to grow and farm sea cucumbers without threatening their populations in the wild. But how do we choose where to set up these ranches?

At the sea ranch in Victory, Bolinao, cultured sandfish juveniles are released in ocean nursery pens and then to the unenclosed area in the 5-hectare sea ranch where they are free-ranging. The sea ranch is guarded and managed by local community partners of the UP MSI. The sea cucumbers are harvested when they attain a marketable size of > 320 g.

Photo: Olivier Josh Caasi

Juvenile sandfish are laid out on a white slate, ready to be weighed and measured.

The sandfish of the BML are routinely weighed by researchers from the Inverts lab.

Photos: Roy C. Corvera

Dr. Marie Antonette Juinio-Meñez and other marine scientists worked with Josephine Dianne Deauna, a graduate of the MS Marine Science program of the UP Marine Science Institute (currently a PhD student in biological oceanography at the University of Hawai’i), to simulate the movement of sea cucumber larvae and determine which locations in western central Philippines are best for housing broodstock. Their paper, “Identification of Priority Sites to Support Management of Commercially Important Sea Cucumber Species by Applying Infomap and Habitat Filters to Larval Dispersal Data,” demonstrates a mathematical approach to conservation research, and how this method can optimize the conduct of field work.

A researcher carefully scrapes off parts of the skin of a Stichopus horrens specimen. This skin sample will be submitted for DNA analysis.

Tomi Garpa and Josh Caasi, two researchers who specialize in marine invertebrates, hold out some of the juvenile sandfish under their care.

Photos: Roy C. Corvera

Simulating larval travel

Most sea cucumbers propagate through sexual reproduction. Both males and females release gametes into the surrounding water, and the sperm cells fertilize the egg cells upon encounter. The larvae formed by this spawning process are swept around by ocean currents until they acquire enough nutrients to grow and metamorphose into a mature state that can settle on the ocean floor (settlement). The sea cucumbers that manage to survive untimely death post-settlement are called “recruits.”

Due to the sheer number of larvae produced during spawning season and the many factors affecting their travel in the water column, it is difficult to predict where matured sea cucumbers are located. Computational marine scientists can bypass this drawback by using models that mimic their ecological behavior. In this study, Deauna used a technique called biophysical modeling, through which she ran a simulation of larval travel with known ocean current and biological patterns.

A biophysical model takes the physical properties of an ecosystem or defined location in the ocean and uses those properties to form a virtual diorama with the help of math translated into lines of code. Scientists can then observe how the model responds to changes in physical and/or biological factors.

Pursuing particles

The team modeled the larval dispersal of two economically viable sea cucumber species in the Philippines, Holothuria scabra and Stichopus horrens. To effectively visualize the traversal of larvae, Deauna treated them as particles.

Stichopus horrens and Holothuria scabra are two species of sea cucumber that are being studied and farmed in the BML in Bolinao, Pangasinan.

The biophysical model was designed using the Connectivity Modeling System (CMS), an open-source software that tracks particles as they move randomly in the ocean. She considered these key factors: the recruitment behavior of sea cucumbers, ocean circulation, and geographical location.

The study’s domain, consisting of the coasts of Palawan, Panay, and Mindoro, was divided into 480 7 km2 squares called grid cells. Twenty random locations were identified within each square as a release point for 1,000 particles, for a total of 20,000 particles released. This number was doubled during the months of April, May, and June to account for the spawning season. These grid cells were then grouped spatially to fit a 10 km radius (the standard size of a marine protected area or MPA), resulting in 127 locations or nodes that were considered for further analysis.

Daily surface current data was obtained from a regional Hybrid Coordinate Ocean Model (HYCOM) model run of the Philippine archipelago.10 Since ocean circulation is highly affected by the seasonal climate, Deauna ran the model under four time periods (northeast monsoon, dry-transition, southwest monsoon, and early northeast monsoon) and observed which oceanographic features would dominate the domain during each period.

To understand the connections among geography, current patterns, and sea cucumber reproduction, the nodes were grouped into coherent regions using the Infomap algorithm. This technique identifies networks of nodes that have the highest probabilities of larval exchange (i.e., larvae released within a coherent region are not likely to settle outside of it.)

Ultimately, the domain was divided into three regions: southeastern Palawan (Group I); eastern Mindoro, western Panay, Romblon, Marinduque, and Magsaysay and Cuyo Islands (Group II); and western and northern Palawan, western Mindoro, and Agutayan (Group III). These groupings are a result of the seasonal fluctuations of the Mindoro Jet, the dominant hydrodynamic feature in the domain. It is strongest during the early northeast monsoon and weakens in this order: southwest monsoon (SWM), northeast monsoon (NEM), and dry transition season. Interestingly, the Mindoro Jet is weakest during the peak spawning season, allowing larvae to drift farther distances.

The 127 nodes in this study’s domain (comprising Palawan, Panay, and Mindoro) were grouped into “coherent regions” — each region having the highest probability of larval exchange within (i.e., larvae released within a coherent region are not likely to settle outside of it.) Three groups were identified as a result. Habitat types present in the nodes were detected using satellite data. Nodes marked with an upright triangle had coral, mangrove, and seagrass; upside-down triangles represent corals and seagrass; diamonds represent mangroves and seagrass; the rectangle represents only seagrass; while circles represent nodes with none of the 3 habitat types.

Generally, there is a high degree of larval exchange among nodes and even across the coherent regions. Sites were classified as either strong sources (where larvae originate) or sinks (where larvae settle). Sites in Group II acted as strong larval sources during the NEM, while those in Group III were strong sources all year round. Most particles from Group II settled in sites within Group I and III. While all locations exhibited source and/or sink characteristics, Group II locations were the main sources and Group I locations, the main sinks. These patterns make sense if you look at the flow of the Mindoro Jet.

The dispersal behavior of the larval particles was mainly influenced by the seasonal fluctuations of the Mindoro Jet, the dominant hydrodynamic feature in the domain.

Not all nodes are equal

The spawning and settlement processes for sandfish typically require seagrass and mangroves, while those for S. horrens need seagrass and corals. Within the domain, 41 nodes featured all three habitat types (seagrass, mangrove, and coral). Twenty-three had seagrass and coral, two had seagrass and mangrove, one had only seagrass, and the remaining 60 had none of the three habitat types. This means that, at a total node count of 64, more locations were suitable for S. horrens. After the simulation, it was also found that S. horrens was distributed throughout a larger area than the sandfish, which had only 43 nodes.

The nodes with suitable habitats were then ranked according to four indices that quantify larval settlement and recruitment. Highly-ranked sites were those that supplied the most larvae (to other sites) while receiving larvae from many other sources. Those sites could be prioritized for protection and broodstock collection, since the spawning adults in those locations were projected to contribute the most larvae to the domain and be genetically different from each other.

“Ecologically meaningful” management areas

Simulating how larval particles travel helped Deauna and her co-authors pinpoint sites in the western central Philippine region that fisherfolk and concerned institutions can prioritize for conservation and farming purposes. They were able to classify sixteen sites around northern Palawan and western Panay as “high-priority.” Locating these sites is an initial step in maximizing strategies for the sustainable management and aquaculture of sea cucumbers. Commercial breeders of sea cucumbers can better manage their broodstock by taking into consideration areas where many sea cucumbers spawn, as well as where most larvae end up.

“Ideally, we would like to protect all 127 locations. Unfortunately, that’s not feasible, so we need to choose which ones to protect first,” Deauna said, emphasizing the importance of ranking the sites. “The locations highlighted in this study can be used by existing hatcheries to further optimize their farming, whether it be for releasing larvae or selecting broodstock.”

She states that the information gathered from the study are just the results of a simulation, and that actual groundwork in confirming its findings and identifying compatible sites is paramount. “It’s baseline information, and only serves as a guide for further studies.”