# Response of blacktip reef sharks Carcharhinus melanopterus to shark bite mitigation products

Madeline Thiele, Johann Mourier, Yannis Papastamatiou, Laurent Ballesta, Eric Chateauminois, Charlie Huveneers

► To cite this version:

Madeline Thiele, Johann Mourier, Yannis Papastamatiou, Laurent Ballesta, Eric Chateauminois, et al.. Response of blacktip reef sharks Carcharhinus melanopterus to shark bite mitigation products. Scientific Reports, 2020, 10 (1), pp.3563. 10.1038/s41598-020-60062-x. hal-03411075

HAL Id: hal-03411075

https://hal.umontpellier.fr/hal-03411075

Submitted on 29 Nov 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
OPEN

# Response of blacktip reef sharks *Carcharhinus melanopterus* to shark bite mitigation products

Madeline Thiele<fnref n="1" />, Johann Mourier²,³,⁴, Yannis Papastamatiou⁵, Laurent Ballesta⁶, Eric Chateauminois⁷ & Charlie Huveneers<fnref n="1" />*

Globally, the frequency of shark bites is rising, resulting in an increasing demand for shark deterrents and measures to lessen the impact of shark bites on humans. Most existing shark protection measures are designed to reduce the probability of a bite, but fabrics that minimise injuries when a shark bite occurs can also be used as mitigation devices. Here, we assessed the ability of the Ocean Guardian Scuba7 and Kevlar material to reduce the likelihood of blacktip reef sharks, *Carcharhinus melanopterus*, from feeding, and to minimise injuries from shark bites. Sharks were enticed to consume a small piece of local reef fish (bait) placed between the two Scuba7 electrodes with the deterrents randomly being turned on or kept off. In the second experiment, the bait was attached to a small pouch made of either standard neoprene or neoprene with a protective layer of Kevlar around it. The Scuba7 reduced the proportion of baits being taken by 67%, (from 100% during control trials to 33%). Sharks also took more time to take the bait when the device was active (165 ± 20.40 s vs. 38.9 ± 3.35 s), approached at a greater distance (80.98 ± 1.72 cm vs. 38.88 ± 3.20 cm) and made a greater number of approaches per trial (19.38 ± 2.29 vs. 3.62 ± 0.53) than when the Scuba7 was inactive. The sizes of punctures from shark bites were significantly smaller on neoprene with Kevlar compared to standard neoprene (3.64 ± 0.26 mm vs. 5.88 ± 0.29 mm). The number of punctures was also fewer when Kevlar was used (14.92 ± 3.16 vs. 74.1 ± 12.44). Overall, the Ocean Guardian Scuba7 and Kevlar reduced the impact of blacktip reef shark bites. These findings may help consumers make informed decisions when purchasing shark deterring and protective products.

Although shark-human interactions remain rare and unlikely events, their frequency has been increasing globally since the 1980s<fnref n="1" />,². The primary cause of this increase is unknown, but factors typically listed include growth in human population and increased ocean use³, habitat destruction and modification, degrading water quality, changing weather patterns and climate change, a change in the abundance and distribution of prey (potentially shifting shark distributions into more human populated areas), and increasing shark abundance⁴,⁵. Whilst increasing, the irregular occurrence of shark bites hinders our ability to evaluate the relative importance of any of these potential causes (but see Afonso, *et al.*⁴ and Meyer, *et al.*⁶). However, as the number of people undertaking aquatic activities continues to rise and shark conservation actions improve, it is likely that shark bite incidents will increase⁷. Whilst irregular, the gruesome nature and potential lethality of shark bites is a focal point for traditional and social media, and entertainment (e.g., films and cartoons)⁸,⁹. This exaggerates the public's perceived risk of shark bites¹⁰, leading to public pressure for the development and implementation of shark bite mitigation measures.

Shark bite mitigation measures are generally categorised in two broad approaches: area protection and personal protection¹¹,¹². The former aims to decrease the spatio-temporal overlap between humans and sharks, i.e. how frequently humans and sharks are in close proximity (e.g., the shark spotters program¹³). However, limitations of area protection strategies, such as weather conditions impacting effectiveness¹⁴ or inapplicability in some

<footnote n="1">¹Southern Shark Ecology Group, College of Science and Engineering, Flinders University, Bedford Park, SA, 5042, Australia. ²MARBEC, Univ Montpellier, CNRS, IFREMER, IRD, Sète, France. ³EPHE, PSL Research University, CRIOBE USR3278 EPHECNRS-UPVD, 66860, Perpignan, France. ⁴Labex Corail, CRIOBE, 98729, Moorea, French Polynesia. ⁵Department of Biological Sciences, Florida International University, North Miami, Florida, USA. ⁶Andromede Oceanologie, Place Cassan, 34280, Carnon, France. ⁷Centre Securite Requin, 25 F avenue des artisans, zone Artisanale de la Pointe des Chateaux, 97436, Saint-Leu, Reunion Island, France. *email: Charlie.huveneers@flinders.edu.au</footnote>
regions, led to a range of personal shark deterrent devices being developed to deter sharks from approaching humans¹². The increasing demand for personal shark deterrents has led to the development of several new products, but few have been scientifically tested. Out of those tested, fewer have shown to significantly affect shark behaviour and reduce probability of shark bites (but see¹²,¹⁵,¹⁶). Most studies have also focused on white sharks, *Carcharodon carcharias*, because it is the species responsible for the most fatal shark bites⁵. However, at least 12 other species are responsible for shark bites⁵, and while bites from most of these other species are not fatal, they can still cause severe injuries and have negative implications to public perception of sharks. Mitigation of shark-human conflict and testing of shark deterrents should therefore not only focus on those species responsible for fatal bites, but should expand to other species.

When shark bites cannot be prevented by area protection or personal deterrents, the ability to reduce the severity of injuries can help in alleviating negative shark-human interactions. One of the most notable pieces of equipment to reduce wounds from shark bites is the chainmail suit developed in the 1980s. The suit seemingly reduced wounds from shark bites and is worn by divers undertaking shark feeding during wildlife tourism¹⁷. However, its heavy weight and lack of flexibility make it impractical in many situations (e.g., surfing, snorkeling) and only usable in limited diving situations. Recent material development has enabled the incorporation of puncture-resistant fabrics with neoprene, e.g., Kevlar. While these fabrics have been shown to reduce cuts¹⁸,¹⁹, their ability to withstand shark bites and reduce injuries is unknown.

This study aimed to test a personal shark deterrent (Ocean Guardian Scuba7; hereafter referred to as Scuba7) and to test a new material proposed to reduced injuries from shark bites (Kevlar). The Scuba7 is one of several products from Ocean Guardian aimed at deterring sharks by overwhelming their ampullae of Lorenzini¹². The Scuba7 differs from other products in the size and configuration of electrodes which are placed on the ankle and tank of the diver (vs the Freedom7 where one of the electrode streams ~1.5m behind the diver or snorkeler). Previous studies assessing the ability of electric deterrents to reduce the probability of shark bites have focused on white sharks, with testing predominantly occurring when the number of sharks present was relatively small, e.g., 1–4¹²,¹⁵,¹⁶. The response to an electric deterrent can, however, vary between species. In addition, competition can change behaviour and hierarchy in sharks²⁰ and decrease the efficacy of magnets as shark deterrents²¹. Testing the effectiveness of electric deterrents on a different species and with a large number of sharks present, potentially affecting shark motivation through competition, will provide a better understanding of electric deterrent efficacy. Kevlar is a material derived from the aromatic polyamide family, which makes it strong, lightweight, flexible, and highly resistant to punctures and lacerations. Previous studies have shown the effectiveness and strength of Kevlar in hypervelocity and high impact situations, resulting in this material being used in bumper shields and bulletproof clothing²²,²³. These properties have the potential to provide better protection for ocean-users when a layer of Kevlar is added onto a standard neoprene wetsuit. Here, we tested the ability of the Scuba7 to prevent a medium-bodied shark species from feeding on a piece of fish and whether the addition of Kevlar on neoprene can reduce injuries such as those occurring when fishers or scientists handle sharks, or during scuba-diving and wildlife tourism activities.

## Methods

**Study location and species.** The blacktip reef shark is a medium-sized shark, up to 1.5 m total length, and is widely distributed on Indo-West Pacific coral reefs²⁴. The species mostly inhabits shallow reefs and sand-flats of atolls and high islands²⁵⁻²⁸ and occasionally non-reef environments²⁹. Blacktip reef sharks are the subject of wildlife tourism, for example in Mo'orea, French Polynesia³⁰, and they are often encountered by scuba-divers and snorkelers and during wildlife tourism operations targeting other species (e.g., bull sharks in Fiji³¹; lemon sharks in French Polynesia³²; fish feeding in Palau; M. Thiele personal observation). These operations often use provisioning to attract sharks in close vicinity of tourists, which has resulted in accidental bites³³,³⁴.

Testing was undertaken on the top of the reef crest at the South pass of Fakarava Atoll, in the Tuamotu archipelago of French Polynesia (−16°52'S, −145°46'W), holding a population of ~100 blacktip reef sharks³⁵. At this site, blacktip reef sharks are fed pieces of meat and fish daily when food is prepared for guests of the Tetamanu village, resulting in up to 30 blacktip reef sharks simultaneously frequenting and being fed next to the food preparation area.

**Field testing.** We conducted tests on two products to assess the ability of mitigating bites from blacktip reef sharks. First, we tested the ability of an electric field-based deterrent developed for divers to repel sharks (i.e., the Scuba7 manufactured by Ocean Guardian). Second, we tested the ability of a layer of Kevlar glued to neoprene wetsuit material (hereafter referred to as Kevlar-neoprene) to reduce punctures or tears.

*Scuba7 testing.* The Scuba7 was deployed on the seabed (80% sand and 20% rocks and broken corals) in ~75 cm of water between 0900 and 1200 hours during seven consecutive days in June 2017. A coloured rope was placed 50 cm around the Scuba7 to facilitate estimates of the distance between sharks and the electrodes (Fig. 1).

Each day, we conducted 12 trials; six with the Scuba7 active and six with the Scuba7 inactive (control). Trials were block randomised so that the status of the Scuba7 (on or off) was randomly selected before each set of two trials. Status of the Scuba7 was confirmed by a green LED light present when turned on. This indicator was not covered during the trials so that the status of the device could be confirmed during the trials. The LED light was not expected to affect sharks as trials occurred during daylight reducing the relative intensity of the LED light, i.e. the intensity of the green light was not strong compared to normal sunlight. During each trial a small piece of local reef fish (bait) was placed between the two electrodes (85 cm from each electrode) using a long pole. Experimenters were not in the water during the trials and stood on a platform ~3 m from the experimental equipment and therefore did not influence shark behaviour. Each trial lasted three minutes or until a shark took the bait. Trial duration was short because of the large number of blacktip reef sharks present and their conditioning
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Diagram showing the experimental set-up and the two dominant current directions during the trials. Image drawn by the authors using SketchUp (www.sketchup.com).</figcaption></figure>

to being fed at that location resulting in sharks quickly taking bait. Preliminary tests showed that a bait lasts less than one minute before being consumed by a shark when the Scuba7 is absent. Trials were filmed by a GoPro Hero3+ video camera and footage was viewed using a VLC player. Coding of the video footage was analysed in a single-blind fashion because the coder did not participate in the trials and had no prior knowledge of the status of the Scuba7 when coding videos. Individual sharks could not be identified due to the lack of distinguishable markings or scars and because most sharks were of similar length. Instead, we recorded the maximum number of sharks during each trial within 30 m of the Scuba7. We used the following terminology to describe shark behaviour and code the videos, modified from previous studies¹²,¹⁵:

**Approach:** a directed swim towards the experimental set-up (each time a shark veered away from the bait and swam back, we classified it as a new approach);

**Approach direction:** whether sharks approach the Scuba7 from parallel to the dipole axis or perpendicular to the dipole axis;

**Distance to electrode:** distance from the sharks' snout, where the shark's sensory organs targeted by the deterrents (ampullae of Lorenzini) are located, to the closest electrode. Distance was visually estimated to the nearest 10 cm using a rope placed at a known length (50 cm) around the Scuba7 (Fig. 1); and

**Reaction:** a behavioural reaction from an individual shark towards the experimental set-up (e.g., tail flick, muscle spasm, head shake, fast direction change).

**Kevlar-neoprene testing.** We tested the Kevlar-neoprene between 0900 and 1200 hours on four days in June 2018. The Kevlar-neoprene consisted of ~1 mm Kevlar knitted sheets (19.5% Kevlar, 70% nylon, 10.5% spandex) glued on both sides of standard neoprene. Sections of Kevlar-neoprene and standard neoprene (without Kevlar) were sawed into 15 × 30 cm pouches on which pieces of meat or local reef fish were attached to entice blacktip reef sharks to bite the neoprene. We also had two neoprene thicknesses for each neoprene type: 5 mm and 3 mm. Neoprene pouches were left in the water until a shark bit the neoprene. Although bite intensity was not recorded, we ensured homogeneity of bites across the material tested by leaving the neoprene until a bite with similar intensity occurred, i.e. head shakes with sharks dragging the pouch in the water, as determined by the senior author (CH). Damages on the neoprene was quantified by a coder who was not present during the trials.

**Statistical analysis.** *Scuba7 testing.* We investigated the effect of the Scuba7 on the following response variables: whether the bait was taken, time taken to take the bait, distance to electrode, number of approaches, and whether a reaction was observed. Following the method in Huveneers *et al.*¹⁵, we minimised potential temporal correlation by testing the effects of the Scuba7 on all response variables using a generalised linear mixed-effects model (GLMM) with *Day* coded as a random effect, the natural logarithm of the maximum number of sharks as offset, and the status of the Scuba7 (on or off) and the approach direction as fixed effects. Current direction was initially included as a fixed factor but was later removed from the model due to strong correlation with approach direction. Dominant current direction was perpendicular to the dipole axis on 4 of the 7 trial days and parallel to the dipole axis on the remaining days, ensuring a relatively balanced current direction across the study period. The error structure of GLMM corrects for non-independence of statistical units due to shared temporal structure and permits the random effects variance explained at different levels of clustering to be decomposed. The inclusion of *Day* as a random effect enabled the analysis to account for the lack of independence in behaviour within each Day. We determined the most appropriate statistical family and error distribution for each analysis by examining the distribution of the response variable and visually inspecting the residuals for the saturated models. Models were fitted with a gaussian distribution for all parameters tested, except to test the effect of the Scuba7 on the likelihood of the bait being taken for which a binomial distribution with logit link function was used. The number of approaches was scaled and centred prior to fitting the GLMM. The full model was dredged, creating a set of new models containing all possible combinations of the factors. These were ranked on decreasing model fit, using Akaike's information criterion corrected for small sample size (AIC<sub>c</sub>)³⁶. The bias-corrected relative weight
of evidence for each model, given the data and the suite of candidate models considered, was the AICc weight; the smaller the weight, the lower its contribution to parameter estimates³⁶.

**Kevlar-neoprene testing.** We also used a GLMM to test whether the following response variables were affected by the Kevlar-neoprene: number of punctures (regardless of whether they went through the neoprene or not); whether a puncture went through the neoprene (i.e., hole); and the length of the hole. Day was again coded as a random effect, with the neoprene type and thickness included as fixed effects.

All models were undertaken in R (v.3.5.0) using the lme and glmer function in the nlme and lme4 packages³⁷,³⁸.

**Electric field modelling.** The effects of undertaking the trials in shallow water was investigated by modelling the electric field propagation across different situations, e.g. in shallow or open water. The electric field was modelled as a dielectric dipole, where each conductive electrode was considered a point load with a charge of q and −q at 1.7 m from each other. Each charge delivers a potential scalar field $V_q$ (Eq. 1):

$$V_q(d) = \frac{q}{4\pi\epsilon_0\epsilon_r} \frac{1}{d} \qquad (1)$$

where $\epsilon_0$ is the permittivity of the empty space; $\epsilon_r$ the relative permittivity of the considered material; and d the distance from the charge q.

The total electric field delivered by the dipole [q, −q] is then calculated by adding the potential scalar field of each electrode, and applying the classic electromagnetism results derived from the Poisson equation (Eqs. 2 and 3):

$$\vec{E}(r) = -\text{grad}(V_{tot}) \qquad (2)$$

$$\vec{E} = E_r \vec{U}_r + E_\theta \vec{U}_\theta \qquad (3)$$

$$E_r = K \left[ \frac{r - a \cos\theta}{(a^2 + r^2 - 2ar \cos\theta)^{3/2}} - \frac{r + a \cos\theta}{(a^2 + r^2 + 2ar \cos\theta)^{3/2}} \right]$$

$$E_{\theta} = K \left[ \frac{a \sin\theta}{(a^2 + r^2 - 2ar \cos\theta)^{3/2}} + \frac{a \sin\theta}{(a^2 + r^2 + 2ar \cos\theta)^{3/2}} \right]$$

with $K = \frac{q}{4\pi\epsilon_0\epsilon_r}$

where $V_{tot}$ is the total scalar field; r the distance from mid-point between the electrode; $\vec{U}_r$ and $\vec{U}_\theta$ the polar projection axes; and a the distance between the electrode and mid-point between the electrodes. This provides the electric field delivered by the Scuba7 in open water.

This electric field was then modified to account for the possible effects of the surface and substrate using the method of image charges. This method states that, for a given distribution of charges, the solution of the Poisson equation is unique. Using this principle, we can add some fictive charges to simulate the presence of a non-conductive field (air or seabed) and compare their effect on the electric field. However, this method does not allow us to model the presence of the surface and seabed concurrently. The electric field propagation was therefore modelled for the following situations:

(1) Mid-water: Scuba7 away from the surface or the seabed;

(2) Surface: Scuba7 was 70 cm from the surface which is considered to be non-conductive. The effect of the seabed was ignored;

(3) Bottom non-conductive: Scuba7 was 5 cm from the seabed which is considered to be non-conductive. The effect of the surface was ignored; and

(4) Bottom conductive: Scuba7 was 5 cm from the seabed which is considered to be conductive. The effect of the surface was ignored.

Cases (3) and (4) give us information about the effect of the conductive parameter of the seabed.

The selected distances to the surface and seabed replicated the position of the Scuba7 in the trials undertaken whereby water depth was ~75 cm with the electrodes being ~5 cm above the seabed. For each situation, the electric field propagation is represented along the horizontal plane positioned 20 cm above the seabed (Z = −55 cm, which is the average swimming depth of sharks during the trials) and in the vertical plane going through the electrodes.

Using these models, we calculated the field strength that sharks were exposed to when crossing through the dipole axis along three trajectories going through the electrode mid-point: parallel, perpendicular, and diagonally. This allowed us to compare the field strengths that sharks were exposed to accounting for the effects of the surface and substrate as well as the effect of approach direction.
<table><thead><tr><th>Level</th><th>β</th><th>SE</th><th>z</th><th>P</th></tr></thead><tbody><tr><td><em>Distance</em></td><td></td><td></td><td></td><td></td></tr><tr><td>intercept</td><td>3.46</td><td>0.51</td><td>6.76</td><td></td></tr><tr><td>Deterrent</td><td>3.55</td><td>0.40</td><td>8.88</td><td></td></tr><tr><td>Approach direction</td><td>3.14</td><td>0.46</td><td>6.79</td><td></td></tr><tr><td>Deterrent* Approach direction</td><td>-1.09</td><td>0.53</td><td>-2.05</td><td></td></tr><tr><td><em>Time</em></td><td></td><td></td><td></td><td></td></tr><tr><td>Intercept</td><td>-0.82</td><td>0.09</td><td>-9.28</td><td></td></tr><tr><td>Deterrent</td><td>1.64</td><td>0.12</td><td>13.12</td><td></td></tr><tr><td><em>Number of approaches</em></td><td></td><td></td><td></td><td></td></tr><tr><td>Intercept</td><td>-1.18</td><td>0.19</td><td>-6.25</td><td></td></tr><tr><td>Deterrent</td><td>1.77</td><td>0.17</td><td>10.23</td><td></td></tr><tr><td><em>Bait consumption</em></td><td></td><td></td><td></td><td></td></tr><tr><td>Intercept</td><td>20.40</td><td>496.71</td><td>0.04</td><td></td></tr><tr><td>Deterrent</td><td>-23.91</td><td>496.71</td><td>-0.05</td><td></td></tr><tr><td><em>Reaction</em></td><td></td><td></td><td></td><td></td></tr><tr><td>Intercept</td><td>-1.53</td><td>0.46</td><td>-3.35</td><td>0.001</td></tr><tr><td>Deterrent</td><td>6.03</td><td>0.67</td><td>9.05</td><td>&lt;0.001</td></tr><tr><td>Approach distance</td><td>-0.07</td><td>0.01</td><td>-5.94</td><td>&lt;0.001</td></tr><tr><td>Approach direction</td><td>0.71</td><td>0.69</td><td>1.03</td><td>0.304</td></tr><tr><td>Deterrent* Approach direction</td><td>-1.99</td><td>0.78</td><td>-2.55</td><td>0.011</td></tr></tbody></table>

**Table 1.** Estimated deterrent level coefficients (β) and their standard errors (SE), z-values of factors included in the top-ranked model (indicated for each variable), and the individual coefficient Type I error estimate (P) when relevant. All models include *Day* as a random effect. Baselines are Scuba7 off and approaches in-between electrodes.

## Results

**Scuba7 testing.** We conducted a total of 84 trials (42 with the Scuba7 active and 42 with the Scuba7 inactive) during which up to 24 blacktip reef sharks were sighted within 30 m of the bait. We recorded 966 approaches, with the mean (±standard error) distance between sharks and an electrode being 74.36 ± 1.61 cm (range: 0–220 cm).

The status of the Scuba7 significantly affected all response variables investigated: whether the bait was taken, time taken to take the bait, distance to electrode, number of approaches, and whether a reaction was observed (Table 1; Fig. 2a–f; Tables S1–5). When the Scuba7 was inactive, sharks took the bait during all trials in an average of 38.9 ± 3.35 s. In contrast, sharks took the bait on 33% of the trials and took about four times longer (165 ± 20.40 s) when the Scuba7 was active (Fig. 2a,b). This resulted in a greater number of approaches per trial (3.62 ± 0.53 vs. 19.38 ± 2.29 when the Scuba7 was turned off or on respectively), which also increased with the number of sharks present during the trials (Fig. 2c; Table S3). The Scuba7 also affected the minimum distance between sharks and electrode, which doubled from 38.88 ± 3.20 cm to 80.98 ± 1.72 cm. The minimum distance was also affected by approach direction, with sharks getting closer when approaching perpendicular to the dipole axis in-between the electrodes (40.28 ± 2.62 cm) compared to in line with the dipole axis (84.95 ± 1.78 cm) (Fig. 2d). Each time a shark consumed the bait with the Scuba7 active, approach direction was perpendicular to the dipole axis and between the electrodes.

We observed reactions 322 times at an average distance of 40.68 ± 1.84 cm (range: 0–160 cm) from an electrode, with most reactions (95%) occurring when the Scuba7 was active. The likelihood of reaction was affected by the combination of the deterrent status and approach direction (Table S5). When the Scuba7 was active, a reaction was more likely to occur when approaching in line with the dipole axis compared to perpendicular (66% vs. 34% respectively), but this effect was no longer apparent when the Scuba7 was inactive (47% vs 53%) (Fig. 2e). The likelihood of a reaction increased as sharks got closer to the electrode, but this was mostly driven by trials when the Scuba7 was active (Table S5; Fig. 2f).

**Kevlar-neoprene testing.** We obtained 22 bites, 12 on Kevlar-neoprene (six on 5 mm, six on 3 mm) and 10 bites on standard neoprene (six on 5 mm, four on 3 mm). The size of holes was significantly smaller when Kevlar-neoprene was used (3.64 ± 0.26 mm) compared to standard neoprene (5.88 ± 0.29 mm; t-value: 3.68; P < 0.0001), but was not affected by neoprene thickness (t-value: -0.29; P = 0.77) or the interaction between neoprene type and thickness (t-value: -1.12; P = 0.26) (Table 2; Fig. 3a), even though both these were included in the top-ranked model (Table S6). In the case of the number of holes, neoprene type was the only factor included in the top-ranked model (wAIC<sub>c</sub> = 0.76; Table S7). The number of holes was significantly less when Kevlar-neoprene was used (14.92 ± 3.16) compared to standard neoprene (74.1 ± 12.44) (Table 2; Fig. 3b). The proportion of holes that went through the fabric was the only variable affected by neoprene thickness, with the interaction between neoprene type and thickness being included in the top-ranked model and having the largest coefficient (Table 2; wAIC<sub>c</sub> = 0.99; Table S8). The proportion of holes through the thick fabric was similar between the standard and
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

**Figure 2.** Effects of the Scuba7 on blacktip reef sharks (*Carcharhinus melanopterus*) on (a) time taken to take the bait; (b) percentage of baits taken; (c) number of approaches per trial; (d) distance to the electrode; (e) body reactions (white = approaches with no reaction, light grey = reaction during outside approaches, darker grey = reaction during in-between approaches); and (f) logistic regression of the likelihood of a reaction occurring as distance from electrode increases for trials when the Scuba7 was off (top) and on (bottom). Mean values are indicated by the crosses; median by the horizontal bar; the length of the box is the inter-quartile range; whiskers represent quartiles; and circles are extreme values.

Kevlar-neoprene (2.8 ± 1.2% and 0.7 ± 0.7% respectively), but differed in the thin fabric, where the standard neoprene (20.0 ± 4.9%) had more holes through the fabric than the Kevlar-neoprene (1.7 ± 1.7%) (Fig. 3c).

**Electric field modelling.** The models show that the electric field propagation was affected by the surface and seabed (Fig. 4). In both planes, the electric fields propagated slightly further when the Scuba7 was close to the surface and the seabed compared to a mid-water situation. For example, 1 Vm⁻¹ was modelled ~1.1 m horizontally from the electrodes in mid-water, but ~1.5 m when close to the non-conductive seabed. However, this difference was not as pronounced close to the electrodes, with the field strength of the Scuba7 in mid-water and close to the surface being similar at ~40 cm and that of the Scuba7 close to the non-conductive seabed propagating only ~15 cm further. The propagation of the electric field was most affected by the conductivity of the seabed, with field strength at 50 cm being ~1 Vm⁻¹ vs. ~10 Vm⁻¹ for the model with non-conductive seabed.
<table><caption>**Table 2.** Estimated deterrent level coefficients (β) and their standard errors (SE), z-values of factors included in the top-ranked model (indicated for each variable), and the individual coefficient Type I error estimate (P). Baselines are Kevlar-neoprene and thick neoprene.</caption><thead><tr><th>Level</th><th>β</th><th>SE</th><th>z</th><th>P</th></tr></thead><tbody><tr><td colspan="5"><em>Size of holes</em></td></tr><tr><td>intercept</td><td>3.77</td><td>0.70</td><td>5.39</td><td>&lt;0.001</td></tr><tr><td>Fabric type</td><td>2.90</td><td>0.79</td><td>3.68</td><td>&lt;0.001</td></tr><tr><td>Fabric thickness</td><td>-0.33</td><td>1.12</td><td>-0.29</td><td>0.770</td></tr><tr><td>Fabric*Thickness</td><td>-1.40</td><td>1.24</td><td>-1.12</td><td>0.261</td></tr><tr><td colspan="5"><em>Number of holes</em></td></tr><tr><td>intercept</td><td>14.92</td><td>7.97</td><td>1.87</td><td>0.076</td></tr><tr><td>Fabric type</td><td>59.18</td><td>11.83</td><td>5.01</td><td>&lt;0.001</td></tr><tr><td colspan="5"><em>Proportion of holes punctured</em></td></tr><tr><td>intercept</td><td>0.01</td><td>0.02</td><td>0.33</td><td>0.744</td></tr><tr><td>Fabric type</td><td>0.02</td><td>0.03</td><td>0.76</td><td>0.455</td></tr><tr><td>Fabric thickness</td><td>0.01</td><td>0.03</td><td>0.35</td><td>0.729</td></tr><tr><td>Fabric*Thickness</td><td>0.16</td><td>0.04</td><td>3.80</td><td>0.001</td></tr></tbody></table>

The field strength that sharks were exposed to varied substantially by the approach direction (Fig. 5). Sharks were exposed to a field <8 Vm⁻¹ when approaching perpendicularly or diagonally to the dipole axis and passing in-between the electrodes. However, electric field exposure increased to 80 Vm⁻¹ as sharks swam on top of the electrodes. The increased electric field exposure was more pronounced when the Scuba7 was close to the seabed and non-conductive, nearly reaching an exposure of 115 Vm⁻¹. The field strength only reached 60 Vm⁻¹ when the seabed was conductive.

Discussion

The global rise in shark attacks has prompted the development of new methods to prevent or reduce the risk and injuries from shark bites. While several studies have previously tested the effects of electric shark deterrents¹²,¹⁵,¹⁶,³⁹, this is the first focused on a small species, the blacktip reef sharks, with many sharks present simultaneously, and undertaken in shallow water (<1 m). Blacktip reef sharks exhibited an obvious behavioural change when exposed to an active Scuba7, with sharks remaining further away from the bait, which was not consumed as often. This resulted in an increased number of attempts to consume the bait. However, blacktip reef sharks were still able to come within close proximity of the active deterrent (80.98 ± 1.72 cm) and consumed the bait in 33% of the trials. This only occurred when sharks approached the bait from in-between the electrodes. Theoretical modelling of the electric field supported these findings with sharks being exposed to weak electrical fields when approaching perpendicularly and in-between electrodes. Testing of the Kevlar-neoprene showed that in situations when personal deterrents do not prevent bites, injuries can be reduced through puncture-resistant fabric. The number and size of holes were all reduced by the Kevlar layer. Combined, these results show that the risk and injuries from blacktip reef sharks can be reduced by the use of personal deterrents and puncture-resistant fabric.

Whilst most studies have tested deterrents on individual sharks¹²,¹⁵,¹⁶, this study tested an electric deterrent around up to 30 sharks at any given time. Previous studies have shown that competition can change behaviour and hierarchy in sharks²⁰ and that high shark densities can result in non-selective mob rushes towards the bait and diminish the effectiveness of shark deterrents²¹. However, the number of sharks present in our study only affected the number of approaches, i.e. more sharks led to more approaches, but did not otherwise hamper the efficacy of the Scuba7. This either suggests that competition does not change the effects of electric deterrents on blacktip reef sharks or that competition did not occur in the situation tested. As testing was conducted in an area where sharks are frequently fed by locals, the regular food supply may have also reduced the need for competition.

Our results confirm a previous study showing that although an electric pulse can deter sharks from consuming a bait or breaching on an intended prey, it cannot prevent all bites¹²,¹⁵,¹⁶. In the present study, blacktip reef sharks only consumed the bait after approaching perpendicular and in-between to the dipole axis. This only occurred when the water current was running perpendicular to the axis of the electrodes and sharks followed the scent in that direction. This direction of approach allowed sharks to come closer to the bait without swimming directly over the electrodes. When sharks approached from a different direction, they were more likely to swim over the electrodes to reach the bait and be deterred when in close proximity to the electrodes. The theoretical models support this, with sharks only being exposed to a strong electrical field (>80 Vm⁻¹) when swimming over the electrodes, while not being exposed to fields >8 Vm⁻¹ when approaching perpendicularly or diagonally and in-between the electrodes. This is likely due to electric field strength rapidly decreasing with distance⁴⁰, with the strength of electric field diminishing from ~115 Vm⁻¹ at the electrodes to 7–16 Vm⁻¹ when 80 cm away from the electrodes¹⁶,⁴¹. These values were obtained by measuring electric fields produced by various Ocean Guardian products (i.e. Freedom7, Surf+ Freedom7) placed at the surface in a 4 m deep sheltered bay¹⁶ and in a 1.2 m deep tank⁴¹, and are similar to values obtained from our theoretical models. While our models show that proximity to the surface and seabed can affect electric field propagation, these effects are relatively small close to the electrodes where sharks are most deterred. This highlights the importance of adequate electrode placement as the ability to deter sharks will depend on the strength of the electric field which is most strongly affected by the distance to the electrodes¹². Overall, this study shows that electric fields have a similar ability to reduce bait consumption
<figure><img src="image_9.png" /><figcaption>**Figure 3.** Damage of blacktip reef shark (*Carcharhinus melanopterus*) bites on standard and Kevlar-neoprene, based on (a) size of the holes; (b) number of holes per trial; and (c) percentage of holes that went through the fabric.</figcaption></figure>

in blacktip reef sharks and white sharks. The similar behavioural response across the two species was observed regardless of the water depth difference between the two studies and proximity of the seabed in the present study. This might be expected as many elasmobranch species have a similar electrosensory threshold sensitivity⁴². However, detection threshold and behavioural response to electric fields can vary extensively between species⁴³⁻⁴⁵. For example, maximum voltage gradient tolerated by *Sphyrna lewini* is 18.5 Vm⁻¹ and is higher than the deterrent threshold suggested for *Triakis semifasciata* (9.6 Vm⁻¹), *C. obscurus* (7–10 Vm⁻¹), and *C. leucas* (3 Vm⁻¹)⁴⁴,⁴⁶. Results from this study should, therefore, not be extrapolated to other species or situations and further testing of electric field-based shark deterrents on other potentially dangerous species (e.g., tiger shark, *Galeocerdo cuvier*, bull shark *C. leucas*) is recommended.

The Kevlar layer glued to the neoprene reduced damages to the neoprene and therefore, Kevlar-neoprene has the potential to minimise injuries and blood loss from shark bites. While the reduced damaged is likely due to the Kevlar layer, it is possible that the method of attachment, i.e. glued on the neoprene, contributed to the reduction of punctures by decreasing neoprene flexibility. Our results also found no differences in the number of punctures and puncture size between thick (5 mm) and thin (3 mm) neoprene. This suggests that a thin layer of Kevlar could be used for protection from shark bites without unnecessarily using a thick layer of neoprene, thus promoting useability and wearability of the Kevlar-neoprene material in a range of applications. Most aquatic activities using wetsuits (e.g., surfing, snorkelling, diving) requires the material to be light and flexible. Although the Kevlar-neoprene is lighter and more practical than the chainmail suit, its suitability across various activities
<figure><img src="image_10.png" /><figcaption>**Figure 4.** Theoretical electric field strength overlayed for electrodes in Mid-water (full lines), close to the surface (dashed lines), close to a non-conductive seabed (dashed-dot lines), and close to a conductive seabed (dotted lines) in (a) horizontal plane (Z = -55 cm) and (b) vertical plane Y = 0 cm). Coloured lines represent the magnitude of the electric field (red: 1 Vm⁻¹; cyan: 10 Vm⁻¹); white lines represent the vector field lines for (b) and a projected vector field line on the horizontal plane for (a).</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 5.** Field strength that sharks were exposed to when crossing through the dipole axis along (a) three different approach directions (parallel = green; perpendicular = red; diagonal = blue) and for (b) four model situations (Mid-water = full line; Surface = dashed line; non-conductive seabed = dashed-dot line; conductive seabed = dotted line). Note that dashed lines are difficult to see because they overlay on top of the full lines.</figcaption></figure>
remains to be tested, but the use of a puncture-resistant layer glued or incorporated to standard neoprene material might bridge the gap between currently used wetsuits and the heavy and relatively impractical chainmail suit.

We are cognisant that the Kevlar-neoprene will not prevent or lessen bone breakages and crushing injuries. However, the primary cause of death in fatal shark bites is blood loss caused from vascular lacerations⁴⁷⁻⁵⁰. The Kevlar-neoprene could therefore still prove effective in reducing haemorrhaging by providing additional time for a shark bite victim to be attended to by emergency services, thus reducing the likelihood of death⁴⁸,⁴⁹. Furthermore, many shark bites around the world are inflicted by small shark species (e.g., blacktip sharks *C. limbatus*; spinner sharks *C. brevipinna*; Caribbean reef sharks *C. perezii*; silky sharks *C. falciformis*)⁵¹ that are less likely to cause crushing injuries but are capable of inflicting damages to tissues, ligaments, and nerves⁵². The ability of Kevlar-neoprene to reduce injuries will differ between species and will be affected by species size, jaw strength, teeth morphology, and bite kinematics. Similarly to the need to test electric field-based deterrents on a range of species, testing Kevlar-neoprene on other potentially dangerous species (e.g., white shark, *Carcharodon carcharias*; tiger shark; bull shark) is warranted.

Whilst many studies focus on bites from large predatory sharks, this study is relevant to smaller sharks, which are also responsible for many shark bites worldwide¹. For example, Florida has the highest number of incidents, with most bites likely to be from blacktip or spinner sharks³. Our findings are also relevant to wildlife tourism and shark feeders who often use gloves or to spearfishers who are also victims of bites from small sharks. We also acknowledge that we did not test damages to skin, muscle, or other material mimicking those, which would normally be under the neoprene. We instead infer that reduced punctures or tear in the Kevlar-neoprene compared to standard neoprene is likely to result in reduced injuries to a person wearing this material. Future studies should validate this assumption.

## Conclusion

The global increase in number of shark bites has led to the development of a range of shark-mitigation measures and personal deterrents rapidly becoming commercially available. Our study unambiguously shows that electric field-based deterrents can change the behaviour of sharks and reduce shark bite risk. In cases when a shark bite still occur, a Kevlar layer glued to neoprene reduces punctures and might therefore reduce injuries and blood loss. Our results will enable private and government agencies to make informed decisions about these devices and to implement shark bite mitigation strategies. Our findings will also inform ocean users and allow them to decide which product might be most suitable.

## Data availability

The datasets generated and analysed during the current study are available from the corresponding author on reasonable request.

Received: 30 July 2019; Accepted: 31 January 2020;

Published online: 27 February 2020

## References

1. Chapman, B. K. & McPhee, D. Global shark attack hotspots: Identifying underlying factors behind increased unprovoked shark bite incidence. *Ocean Coast. Manage.* **133**, 72–84 (2016).

2. McPhee, D. Unprovoked shark bites: Are they becoming more prevalent? *Coast. Manage.* **42**, 478–492 (2014).

3. Burgess, G., Buch, R., Carvalho, F., Garner, B. & Walker, C. *Factors contributing to shark attacks on humans: A Volusia County, Florida, case study.* **1** 541–565 (CRC Press, 2010).

4. Afonso, A. S., Niella, Y. V. & Hazin, F. H. Inferring trends and linkages between shark abundance and shark bites on humans for shark-hazard mitigation. *Marine and Freshwater Research* **68**, 1354–1365 (2017).

5. West, J. G. Changing patterns of shark attacks in Australian waters. *Marine and Freshwater Research* **62**, 744–754, https://doi.org/10.1071/MF10181 (2011).

6. Meyer, C. G. et al. Habitat geography around Hawaii's oceanic islands influences tiger shark (Galeocerdo cuvier) spatial behaviour and shark bite risk at ocean recreation sites. *Scientific reports* **8**, 4945 (2018).

7. Merchant, J. In *AGU Fall Meeting Abstracts*. 3–13 (Advancing Earth and Space Science) (2017).

8. Muter, B. A., Gore, M. L., Gledhill, K. S., Lamont, C. & Huveneers, C. Australian and US news media portrayal of sharks and their conservation. *Conservation Biology* **27**, 187–196, https://doi.org/10.1111/j.1523-1739.2012.01952.x (2013).

9. Sabatier, E. & Huveneers, C. Changes in Media Portrayal of Human-wildlife Conflict During Successive Fatal Shark Bites. *Conservation and Society* **16**, 338–350 (2018).

10. Crossley, R., Collins, C. M., Sutton, S. G. & Huveneers, C. Public perception and understanding of shark attack mitigation measures in Australia. *Human dimensions of wildlife* **19**, 154–165 (2014).

11. Gray, G. M. & Gray, C. A. Beach-User Attitudes to Shark Bite Mitigation Strategies on Coastal Beaches; Sydney, Australia. *Human Dimensions of Wildlife* **22**, 282–290 (2017).

12. Huveneers, C. et al. Effects of an electric field on white sharks: *in situ* testing of an electric deterrent. *PloS one* **8**, e62730, https://doi.org/10.1371/journal.pone.0062730 (2013).

13. Engelbrecht, T., Kock, A., Waries, S. & O'Riain, M. J. Shark Spotters: Successfully reducing spatial overlap between white sharks (*Carcharodon carcharias*) and recreational water users in False Bay, South Africa. *PloS one* **12**, e0185335, https://doi.org/10.1371/journal.pone.0185335 (2017).

14. Curtis, T. H. et al. Responding to the risk of White Shark attack. *Global Perspectives on the Biology and Life History of the White Shark*. CRC Press, 477–510 (2012).

15. Huveneers, C. et al. Effectiveness of five personal shark-bite deterrents for surfers. *PeerJ* **6**, e5554 (2018).

16. Kempster, R. M. et al. How Close is too Close? The Effect of a Non-Lethal Electric Shark Deterrent on White Shark Behaviour. *PloS One* **11**, e0157717, https://doi.org/10.1371/journal.pone.0157717 (2016).

17. Warner, A. J. *Relative tensile strengths of chainmail weaves*, Massachusetts Institute of Technology, (2015).

18. Bao, L., Sato, S., Morioka, H. & Soma, S. Improving stab-resistant textile materials with a non-woven fabric structure. *Journal of Textile Engineering* **62**, 37–42 (2016).

19. Zhu, D., Mobasher, B. & Rajan, S. In *Dynamic Behavior of Materials, Volume 1* 147–152 (Springer, 2011).

20. Brena, P. F., Mourier, J., Planes, S. & Clua, E. E. Concede or clash? Solitary sharks competing for food assess rivals to decide. *Proceedings of the Royal Society B: Biological Sciences* **285**, 20180006 (2018).
21. Robbins, W., Peddemors, V. & Kennelly, S. Assessment of permanent magnets and electropositive metals to reduce the line-based capture of Galapagos sharks, Carcharhinus galapagensis. *Fisheries Research* **109**, 100–106, https://doi.org/10.1016/j.fishres.2011.01.023 (2011).

22. Hayhurst, C. J., Hiermaier, S. J., Clegg, R. A., Riedel, W. & Lambert, M. Development of material models for Nextel and Kevlar-epoxy for high pressures and strain rates. *International Journal of Impact Engineering* **23**, 365–376 (1999).

23. Tham, C., Tan, V. & Lee, H.-P. Ballistic impact of a KEVLAR® helmet: Experiment and simulations. *International Journal of Impact Engineering* **35**, 304–318 (2008).

24. Compagno, L., Dando, M. & Fowler, S. A field guide to the sharks of the world (2005).

25. Hobson, E. S. Feeding behavior in three species of sharks (1963).

26. Mourier, J., Vercelloni, J. & Planes, S. Evidence of social communities in a spatially structured network of a free-ranging shark species. *Anim. Behav.* **83**, 389–401 (2012).

27. Papastamatiou, Y., Caselle, J., Friedlander, A. & Lowe, C. Distribution, size frequency, and sex ratios of blacktip reef sharks Carcharhinus melanopterus at Palmyra Atoll: a predator-dominated ecosystem. *Journal of Fish Biology* **75**, 647–654 (2009).

28. Speed, C. W. et al. Spatial and temporal movement patterns of a multi-species coastal reef shark aggregation. *Marine Ecology Progress Series* **429**, 261–275 (2011).

29. Chin, A., Heupel, M., Simpfendorfer, C. & Tobin, A. Ontogenetic movements of juvenile blacktip reef sharks: evidence of dispersal and connectivity between coastal habitats and coral reefs. *Aquat. Conserv: Mar. Freshwat. Ecosyst.* **23**, 468–474 (2013).

30. Gallagher, A. J. & Hammerschlag, N. Global shark currency: the distribution, frequency, and economic value of shark ecotourism. *Current Issues in Tourism* **14**, 797–812 (2011).

31. Brunnschweiler, J. M. & Baensch, H. Seasonal and long-term changes in relative abundance of bull sharks from a tourist shark feeding site in Fiji. *PLoS One* **6**, e16597 (2011).

32. Clua, E., Buray, N., Legendre, P., Mourier, J. & Planes, S. Business partner or simple catch? The economic value of the sicklefin lemon shark in French Polynesia. *Marine and Freshwater Research* **62**, 764–770 (2011).

33. Clua, E., Buray, N., Legendre, P., Mourier, J. & Planes, S. Behavioural response of sicklefin lemon sharks *Negaprion acutidens* to underwater feeding for ecotourism purposes. *Marine Ecology Progress Series* **414**, 257–266 (2010).

34. Randall, J. E. & Helfman, G. S. Attacks on humans by the blacktip reef shark (*Carcharhinus melanopterus*) (1973).

35. Mourier, J. et al. Extreme inverted trophic pyramid of reef sharks supported by spawning groupers. *Curr. Biol.* **26**, 2011–2016 (2016).

36. Burnham, K. P. & Anderson, D. R. *Model selection and multimodel inference: a practical information-theoretic approach*. (Springer-Verlag, 2002).

37. Bates, D., Kliegl, R., Vasishth, S. & Baayen, H. Parsimonious mixed models. *arXiv preprint arXiv:1506.04967* (2015).

38. Pinheiro, J., Bates, D., DebRoy, S. & Sarkar, D. R. Development Core Team. nlme: linear and nonlinear mixed effects models. *R package version* **3**, 1–103 (2010).

39. Broad, A., Knott, N., Turon, X. & Davis, A. R. Effects of a shark repulsion device on rocky reef fishes: no shocking outcomes. *Marine Ecology Progress Series* **408**, 295–298, https://doi.org/10.3354/meps08591 (2010).

40. Kalmijn, A. J. The electric sense of sharks and rays. *Journal of Experimental Biology* **55**, 371–383 (1971).

41. Chateauminois, E., Hoarau, M. & Maillard, F. Innovative projects of the shark risk reduction resource and support center (CRA)-results of experimental tests on individual electrical impulse repulsion equipment-final report. *CRA-V* **2**, 1 (2019).

42. Jordan, L. K., Mandelman, J. W. & Kajiura, S. M. Behavioral responses to weak electric fields and a lanthanide metal in two shark species. *Journal of Experimental Biology and Ecology* **409**, 345–350 (2011).

43. Kajiura, S. M. & Holland, K. N. Electroreception in juvenile scalloped hammerhead and sandbar sharks. *Journal of Experimental Biology* **205**, 3609–3621 (2002).

44. Marcotte, M. M. & Lowe, C. G. Behavioral responses of two species of sharks to pulsed, direct current electrical fields: testing a potential shark deterrent. *Marine Technology Society Journal* **42**, 53–61 (2008).

45. O'Connell, C. P., Abel, D. C., Gruber, S. H., Stroud, E. M. & Rice, P. H. Response of juvenile lemon sharks, Negaprion brevirostris, to a magnetic barrier simulating a beach net. *Ocean Coast. Manage.* **54**, 225–230 (2011).

46. Smith, E. Electro-physiology of the electrical shark-repellant. *The Transactions of the Institute of Electrical Engineers* **65**, 1–20 (1974).

47. Lentz, A. K. et al. Mortality and Management of 96 Shark Attacks and Development of a Shark Bite Severity Scoring System. *The American Surgeon* **76**, 101–106 (2010).

48. Ballas, R., Saetta, G., Peuchot, C., Elkienbaum, P. & Poinsot, E. Clinical features of 27 shark attack cases on La Réunion Island. *Journal of Trauma and Acute Care Surgery* **82**, 952–955, https://doi.org/10.1097/ta.0000000000001399 (2017).

49. Woolgar, J. D., Cliff, G., Nair, R., Hafez, H. & Robbs, J. V. Shark attack: review of 86 consecutive cases. *Journal of Trauma and Acute Care Surgery* **50**, 887–891 (2001).

50. Clua, E. & Reid, D. In *Post Mortem Examination and Autopsy-Current Issues From Death to Laboratory Analysis* (InTech, 2018).

51. French, L. (Florida Museum of Natural History University of Florida Gainesville, FL 32611 USA, 2018).

52. Hazin, F. & Afonso, A. Response: A conservation approach to prevention of shark attacks off Recife, Brazil. *Animal Conservation* **17**, 301–302 (2014).

## Acknowledgements

The study was approved by Flinders University Animal Ethics Committee approval E446. This study was conducted under the authority of the Direction de environment de Polynésie française (DIREN) under the convention 4469/MCE/ENV signed on 11 June 2017. Funding was provided by New South Wales Department of Primary Industries Shark Management Strategy to C. Huveneers and Blancpain Ocean Commitment to L. Ballesta. Sponsors supporting the expedition include: Tahiti Tourisme, Air Tahiti Nui, Aqualung, Nikon, Seacam, Aqualung and ApDiving. We thank S. and A. Richemond, and the staff from Tetamanu lodging for hosting the research team and providing technical assistance during the field work.

## Author contributions

C.H. conceived the study; C.H., J.M., Y.P., L.B. undertook the study; M.T., C.H., E.C. analysed the data; M.T., C.H., E.C., J.M., Y.P. wrote the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-60062-x.

Correspondence and requests for materials should be addressed to C.H.

Reprints and permissions information is available at www.nature.com/reprints.
**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_12.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020