Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Grouper (Epinephelidae) spawning aggregations affect activity space of grey reef sharks, Carcharhinus amblyrhynchos, in Pohnpei, Micronesia

  • Kevin L. Rhodes ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Writing – original draft, Writing – review & editing

    klrhodes_grouper@yahoo.com

    Affiliations MarAlliance, San Francisco, CA, United States of America, Pacific Marine Science and Conservation, Grass Valley, CA, United States of America

  • Ivy Baremore,

    Roles Formal analysis, Methodology, Writing – original draft, Writing – review & editing

    Affiliation MarAlliance, San Francisco, CA, United States of America

  • Rachel T. Graham

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Writing – original draft, Writing – review & editing

    Affiliation MarAlliance, San Francisco, CA, United States of America

Abstract

Fish spawning aggregations (FSA) act as biological hotspots that concentrate food and nutrients across a broad trophic spectrum. In Pohnpei (Federated States of Micronesia), 20 female grey reef sharks (Carcharhinus amblyrhynchos) were acoustically tagged at two multi-species grouper (Epinephelidae) FSA to examine the likelihood that these mesopredators utilize FSA as a seasonal food source. Both FSA sites are within small-scale MPAs, thus providing a secondary opportunity to examine their conservation potential during these ephemeral events. Shark movement and residency was gauged against known spatial and temporal grouper reproductive patterns using an array of 15 and 50 acoustic receivers at Ant Atoll and Pohnpei (Island), respectively. Activity space was investigated using Kernel Density estimates of individual sharks, and residency indices (RI) were analyzed based on daily and monthly occurrence at the array. Three distinct residency patterns were identified: transient, semi-transient, or resident (Daily RI <0.40, >0.40 to 0.80, or >0.80, respectively). Generalized linear mixed models (GLMMs) were used to identify biological and environmental factors influencing shark activity space, including month, temperature, shark size, spawning month, and residency pattern. Findings revealed significant changes in average monthly residency indices and kernel densities during spawning months in support of an opportunistic foraging strategy around FSA. Monthly residency was higher during spawning months among semi-resident and transient sharks, while average monthly activity space was concentrated around FSA. Best-fit models for the GLMM indicated that activity spaces were most influenced by month and grouper spawning month. Seven of 20 sharks demonstrated inter-island movement and wide variations in individual movement and spatial requirements were shown. The concentration of sharks and groupers at unprotected FSA sites increases their vulnerability to fishing and supports the need for combined area and non-area management measures to effectively protect these species.

Introduction

Fish spawning aggregations (FSA) represent critical life history phases for a number of coral reef fishes [1] and serve as biological hotspots that provide food and nutrients to marine organisms across a wide trophic spectrum [27]. In many locales, FSA are multi-species and may be comprised of 100s or 1000s of individuals [810], thus representing a rapid and substantial increase in biomass and nutrient flow within areas used as spawning sites. Since FSA are spatially and temporally predictable, they are highly attractive to both fishers who can extract high catch volumes over brief time periods [1113] and to marine organisms that gain from the elevated food and nutrients available during these events [14]. For fisheries, the benefits of FSA fishing are well known [13, 15] as are the impacts that targeted FSA fishing can have on stocks [16, 17], particularly when FSA remain unprotected. In extreme instances, aggregation fishing can fully extirpate FSA [18], with concomitant losses in food and nutrients to the organisms reliant on them.

Among the many species benefitting from FSA through foraging are sharks and rays, with perhaps the most spectacular foraging events associated with whale sharks (Rhincodon typus) [4, 1922] and reef-associated mantas (e.g. Alfred manta, Mobula alfredi) [23]. Less observed, but equally spectacular are the foraging events on FSA by subtropical and tropical reef-associated requiem sharks (Carcharhinidae), including lemon sharks (Negaprion brevirostris; Poey 1868) [24], grey reef sharks (Carcharhinus amblyrhynchos; Bleeker 1856) [14, 25], and blacktip reef sharks (Carcharhinus melanopterus; Quoy and Gaimard 1824) [26]. During these periods, sharks may become vulnerable to the same types of gear used for target aggregating fish (Rhodes pers. observ.), thus protective measures that target spawning fishes may also provide a temporary respite from fishing for sharks.

In the tropics, there is a drive to establish marine protected areas (MPAs) as a means to protect reef-associated marine resources. Similarly, the International Marine Protected Area Congress seeks to designate 10% of the world’s oceans as MPAs by 2020 [27]. Thus, the need to examine the utility and optimal design of these set-asides for highly mobile and migratory organisms, such as transient spawners and reef-associated sharks is imperative. For transient FSA-forming species, a number of conventional and acoustic telemetry studies has been conducted to examine seasonal and monthly movement patterns, residency, visitation frequency and fishing vulnerability [2833]. Those studies largely confirm the utility of MPAs in protecting aggregating fish, but point to the need for additional management measures to protect migrating fish and fish using aggregation sites. Similarly, numerous tagging studies have focused on reef-associated and pelagic shark movement, both to assess MPA effectiveness and identify other appropriate conservation measures [24, 3438].

Among coastal sharks, grey reef sharks are one of the most prominent tropical reef-associated Indo-Pacific species. Several studies have investigated their movement [3438] and the potential benefits of small-scale MPAs (herein ≤1 to 5 km2). In general, findings have shown sex-specific differences in habitat preferences and movement patterns and distances, with females generally less likely to undergo long-distance movement than males. Grey reef sharks also tend to be more resident when reefs are continuous or closely spaced [37]. Where reefs are semi-isolated, grey reef sharks appear to be largely resident with some long-distance movement possible [3840]. In contrast, in isolated reef settings, grey reef shark residency can be high [36], but movements away from isolated habitats may also be substantial, covering hundreds of kilometers [41].

In the western and central Pacific, the Micronesia Challenge was initiated in 2005 with the primary goal of incorporating 30% of available nearshore marine habitat into marine reserves (http://micronesiachallenge.org). While there has been success in achieving this goal, most reserves are small (<10 km2), nearly all are unenforced and only a few have been placed in areas of biological significance, i.e. critical habitat or areas of high biodiversity. Indeed, a recent assessment showed a limited efficacy of MPAs in the region [42]. Currently, it is unclear whether the intended benefits of these reserves will be realized, particularly for mobile animals, such as reef-associated sharks and aggregating groupers [43]. Questions about the efficacy of MPAs for sharks have arisen based on empirical evidence. For example, in the Great Barrier Reef grey reef sharks have continued to experience large-scale declines (17% yr-1) in spite of relatively strict enforcement in four types of management zones that include no-entry reserves [44]. In contrast, subsequent fishery-dependent findings [45] reported that current management for sharks was effective, although the majority of shark by-catch was grey reef sharks.

Similar to the use of MPAs for sharks, the results of small-scale MPAs for FSA have been mixed. In Pohnpei, Micronesia, year-over-year declines in fish density have continued inside a small-scale MPA meant to protect spawning grouper [10]. In contrast, improvements in fish density have been shown for MPA-protected FSA in Papua New Guinea [46]. In the Caribbean, improvements in FSA abundance and size were demonstrated following combined area and harvest measures for red hind (Epinephelus guttatus) [47]. Clearly, there is a need for greater understanding of habitat use, movement and vulnerability of these mobile organisms.

Although grey reef sharks are among the most common reef-associated mesopredators, they are also among the most vulnerable to small-scale fisheries. Grey reef sharks were listed by the International Union for the Conservation of Nature (IUCN) as Near Threatened (NT) during the last global evaluation in 2005. Population trends at that time were listed as unknown [48]. In recent years, declines in grey reef sharks have been associated with targeted and non-targeted fisheries that are drawn to the commercial value of dried fins, meat and teeth. In many small-scale commercial fisheries in the tropics, grey reef sharks are commonly taken as by-catch [49] and often killed and discarded because of the perception by some fishers that they impede fishing success (Rhodes pers. observ.). Among larger scale commercial fisheries, passive sexual selection of sharks has been reported from longline fishing [5051], due in part to segregation by depth in some shark species [52]. Moreover, high site fidelity and philopatry have also been identified among coastal shark species [5355], suggesting localized depletions are possible where sharks are concentrated and remain unprotected. As a result of these various shark fisheries and observed sex-specific behavioral differences [5657], populations in some locales have been severely depleted [58]. Even with fisheries management in place, grey reef sharks may be vulnerable to fishing mortality as they have a moderate rebound potential, with gestation over 9-months, litter size ranging from 3–6 pups, and a mean of 4.1 pups (e.g. Hawaii) [59].

In Pohnpei (Federated States of Micronesia, FSM hereafter), populations of all nearshore shark species are anecdotally reported as declining (Rhodes, pers. observ.). Local fisheries for sharks include the capture of juveniles for teeth and jaws for the curio trade, however on the main island and atolls immediately surrounding Pohnpei (Island), there is no targeted food fishery for sharks. Small-scale commercial and subsistence fishers view sharks of all kinds as a nuisance and are known to remove the tail and discard the carcass when landed (Rhodes, pers observ.). Some small-scale shark finning is known to occur within the artisanal fishery. Further investigation is needed to gauge the scale and impact of the existing fishery, notably in light of the ban on shark fishing in Pohnpei’s territorial waters enacted by the FSM Government in 2013. The national ban was expanded as part of the declaration of the Regional Micronesia Shark Sanctuary in 2015, which encompasses 6.5 million square kilometers, and where the region’s governments prohibit the commercial fishing of sharks, retention of sharks caught as by-catch, and the trade, possession, and sale of shark products. The impact of the sanctuary on a reversal in declines of shark populations has not been evaluated.

Within the state, there are currently 16 small-scale no-take marine reserves ranging in coverage from 0.34 km2 to 41.4 km2 (mean = 7.4 km2; median = 2.1 km2); however, most are poorly enforced [10, 60], including the Kehpara Marine Sanctuary (KMS) (Fig 1). The KMS was established in 1999 to protect seasonal spawning aggregations of camouflage grouper, Epinephelus polyphekadion (Bleeker 1849; VU A2bd), brown-marbled grouper, Epinephelus fuscoguttatus (Forsskål 1775; VU A2bd+4bd), and squaretail coralgrouper, Plectropomus areolatus (Rüppell 1830; VU A2bd) (http://www.iucnredlist.org).

thumbnail
Fig 1. Map of the study site.

Pohnpei and Ant Atoll, Federated States of Micronesia showing receiver locations. Boxes with yellow shading indicate the boundaries of no-take marine reserves and land masses are shaded grey. Orange boxes denote fish spawning aggregation (FSA) sites.

https://doi.org/10.1371/journal.pone.0221589.g001

Four sections of nearby Ant Atoll were also designated for marine area protection by Pohnpei State in 2007 (Fig 1) and in 2008 the entirety of Ant (34 km2) became a UNESCO Man and Biosphere Reserve (http://www.unesco.org). The lagoon and channel are now regularly patrolled by private enforcement. Within the KMS, decadal declines in the FSAs of all three aggregating groupers have been associated with fishing, following a relaxation in enforcement [10]. There has been no baseline assessment of aggregating groupers at Ant or of sharks at any location in Pohnpei.

Within both the KMS and the main channel of Ant Atoll, seasonal (January-May) FSA of these aforementioned three groupers form c. two weeks prior to each full moon [35]. Aggregation abundance and density increases daily for all species prior to spawning and the arrival time and duration of individuals varies by sex [6163]. Spawning occurs on the last day or days of the FSA for each respective species, after which time individuals disperse to home reef areas. Each aggregation typically numbers in the 100s to 1,000s of individuals. At Ant, a number of other teleosts also aggregate during these months, creating an additional attraction to fishers and possibly sharks. Sharks have always been present at the Kehpara and Ant FSA during past research studies [63] and were often captured and released, however it is unclear whether this presence simply represented resident grey reef sharks or whether there was an attraction to the site for trophic opportunism.

The objectives of the current study were to:

  1. examine monthly activity spaces, as kernel densities, residency and timing of grey reef shark movement in order to discern whether there are significant differences between grouper non-spawning and spawning periods to imply attraction to the site for trophic interactions, and;
  2. identify what level of protection the small-scale MPAs that characterize much of the developing Pacific, including the existing MPAs at Kehpara and Ant Atoll, might provide for grey reef sharks times when they are concentrated, e.g. during grouper reproductive periods.

Materials and methods

Site description

Between 22 and 26 January 2010 and 31 January to 16 March 2012, grey reef sharks were captured at two separate seasonal (transient) grouper spawning aggregation sites each located within small-scale MPAs (Fig 1). In Pohnpei, the KMS was established by a private landowner as a permanent no-take zone in 1995 and formalized by Pohnpei State in 1999. The KMS encompasses a well-studied multi-species grouper FSA site that has been described in detail in previous reports [35, 3941]. The KMS is ca. 1.5 km2 and includes both inner and outer reef habitat along the southwestern portion of Pohnpei (Kitti Municipality). Within the KMS, groupers aggregate monthly on the outer reef wall over an approximate 2-week period prior to full moon between January and May: P. areolatus—January-May; E. fuscoguttatus—February-May; E. polyphekadion—February-March or March-April. The initiation of each species’ spawning season and the month of peak density are all dependent on the timing of full moon relative to winter solstice [10]. In this way, the reproductive season for E. polyphekadion, for example, can be either February-March or March-April. For P. areolatus, there are indications the FSA in January represents an all-male aggregation that starts around new moon and remains in place over a 6-week period [62]. In subsequent months, male P. areolatus depart along with females following spawning around the full moon. All three species have been shown or observed utilizing common reproductive migratory corridors to reach and depart the FSA site and at last some individuals travel up to 25 km from the FSA to reach home reefs immediately after spawning [6264]. As previously indicated, the species-specific duration of FSA varies, with P. areolatus and E. fuscoguttatus maintaining aggregations over c. 2 weeks during each reproductive month, while the E. polyphekadion FSA currently extends over 5 days or less [10]. This truncated aggregation period is considered an effect of overfishing and represents a decline from two decades earlier when these FSA formed and persisted for 10–12 days [61]. Depending on the timing of seasonal aggregation formation, February, March and April are peak aggregation months in terms of abundance, and months where either two or all three species can be present.

Prior acoustic telemetry studies and fishery-independent surveys for grouper at the KMS showed that grey reef sharks inhabit the area and readily prey on line-caught groupers during retrieval. Enforcement at KMS occurred from the time of its designation as an MPA until 2005 when the private landowner died. Since that time, poaching has been ongoing, with year-over-year declines in grouper abundance and density [10]. Another FSA site, located in the channel of neighboring Ant Atoll was chosen as the second location for grey reef shark tagging. Unlike Kehpara, the Ant FSA was largely unprotected and actively fished until 2015, particularly during grouper spawning times. Ant Sanctuary was established as a no-take zone in 2007 and is comprised of four sections totally 34 km2, however only about one-quarter of the total area is dedicated to the channel where the FSA forms (Fig 1). The site is currently being enforced by private individuals residing on the atoll.

Tagging

For tagging, grey reef sharks were captured by traditional fishers fishing on snorkel. Fishing was conducted in open water seaward of the reef. Fishers used standard 80-lb clear monofilament tuna line fitted with braided 1/16” stainless steel wire leader and a swivel attached to a 16/0 gauge circle hook. Fishing was conducted during daylight hours using various baits, including grouper, skipjack (Katsuwonus pelamis) or kawakawa (Euthynnus affinis). Once captured, grey reef sharks were brought to the surface, held boatside on the fishing line to allow water movement across the gills, immobilized using a tail rope and subsequently inverted to induce tonic immobility and facilitate tag implantation. Each captured shark was then measured across the body length in a curved line to the nearest cm pre-caudal (PCL) and total length (TL) [65], and sexed visually prior to surgery. For implantation of acoustic tags (Vemco V16-6L coded acoustic transmitter, 69 kHz, 90-sec blanking interval, estimated tag life = 1877 d) (Amirix Systems, Nova Scotia, Canada), a c. 4–5 cm incision into the peritoneal cavity was made using a sterile #10 surgical scalpel blade. Prior to insertion, tags were coated with over-the-counter triple antibiotic ointment. Following insertion, incisions were closed using Ethicon 5–0 cutting edge needle and braided silk sutures. Tagging procedures typically lasted c. 8 to 10 min. All fish were checked to ensure they were healthy and active prior to release at or near the point of capture. Of the nine tags deployed in 2010, one was recovered from an animal previously tagged in a separate study and had a maximum of 1411 days remaining battery life. Twelve new tags were implanted in grey reef sharks in 2012, bringing the study total to 21 animals. Pohnpei State Office of Fisheries and Aquaculture reviewed and approved the protocols and provided oversight of the research.

Acoustic monitoring

Between January 2010 and October 2013, monitoring was conducted on tagged grey reef sharks using 65 Vemco VR2W receivers that were moored along the outer barrier reefs and channels leading into Ant and Pohnpei lagoons (Fig 1). A concentration of receivers was placed within and adjacent to the KMS as part of a separate, ongoing grouper tagging study. Otherwise, receivers were placed north and south of the FSA at c. 4-km intervals. The receivers were placed to simultaneously research grey reef sharks and aggregating groupers [63] and covered c. 80 km of combined reef area. All receivers were suspended in the water column at c. 15–20 m using a hard plastic sub-surface float. All receivers were downloaded and maintained at 1-year intervals or less, depending on opportunity and need. At Ant Atoll, 15 receivers remained deployed along 45 km of outer barrier reef, inclusive of two receivers along the main channel leading into the lagoon (Fig 1). Prior range testing showed that detections of 250 m were common in unobstructed areas [64], similar to that which characterizes the outer reef mooring sites.

Data analysis

Total detections at each receiver within the array were examined for the study period. Detections at Pohnpei and Ant Atoll were examined by day and month to examine changes in shark presence over time and determine activity spaces. Daily detections (shark detected at the array at least twice during one 24-hr period) of individual sharks were calculated to monitor presence of the sharks within the receiver array and at the FSA, and to assess changes in detection frequency during grouper spawning months. Residency of sharks (residency index = RI) was calculated as the total number of days detected divided by the total number of days the shark was expected at the array if it were detected every day (date array removed-date tagged) [66]. RI was also calculated for days spent within the FSA and the KMS, and during new and full moon cycles. Sharks were considered “resident” when overall RI > 0.80, “semi-transient” when RI > 0.40 and < 0.80, and “transient” when RI < 0.40, similar to the ‘cut-off’ values used by other recent researchers [14]. A monthly RI was also calculated for each shark to determine how individuals’ residency patterns changed over the course of the year. Monthly RI was plotted separately for the resident groups to assess whether these patterns differed among transient, semi-transient, and resident sharks at Ant Atoll and Pohnpei. Differences in RI by month and transient group were examined using Analysis of Variance (ANOVA) if data were determined to be normally distributed and homogenous, otherwise a Kruskal-Wallis test was applied [67]. Monthly RI differences between grouper spawning months and non-spawning months were tested using a Student’s t-test. Differences in mean daily RI between new and full moon periods at the FSA were also tested using a Student’s t-test. Lunar phases were determined using the R package Lunar [68].

The utilization distribution (UD) of individual sharks was examined using fixed Kernel Density (KD) analyses of detections for each individual at 30-minute time steps using a bivariate normal kernel function [69]. The KD calculation uses a kernel method to estimate the KD using relocation (detection) data [70], and a smoothing parameter was chosen based on visual choice after several trials. Core activity space (50% KD) and the activity space extent (95% KD) were examined for changes in habitat use over the course of 46 months. Sequential monthly KD values were calculated for each individual; in instances where the number of detections was insufficient to calculate KD, the individual was omitted from analysis for that month. All calculations and KD were made using the adehabitatHR and sp packages in R [69, 7172], and visualizations for KD were made in QGIS.

Generalized linear mixed models (GLMMs) were used to examine the monthly changes in activity spaces of sharks in relation to seasonal (January—May) grouper spawning aggregations and other factors, such as water temperature, shark size (TL), and residency status of individuals (transient, semi-transient, resident). Individual (transmitter code) was included in the model as a random factor to account for repeated-measures [73]. Estimates of KD were tested for normality and square root transformed, if necessary [74], and models were checked for multicollinearity by calculating variance inflation factors (VIF). Final model selection was based on Akaike Information Criteria (AIC) [73, 75], and the dredge function [76] (package MuMIn) was used to generate a suite of models for comparison [77]. Models were built with all combinations of factors for 50% and 95% KD (global model: KD ~ month + spawning month + temperature + TL + residency status) and were tested against the null model using maximum likelihood. Model averaging on the top candidate models was conducted: parameter estimates were averaged for a subset of models where deltaAIC < 2, and cumulative AIC weights (wi) were used to determine the influence of each factor on activity space, with wi > 0.5 considered significant drivers of activity space [73, 77]. Month was a sequential factor and spawning months were binary (spawning months = January through the end of May of each year), based on previously published information [10]. Temperature was a calculated monthly average based on data loggers that were concurrently placed at the study area.

Differences in monthly estimated 50% and 95% KD were analyzed using a non-parametric Kruskal-Wallis test because the data were non-normal and non-homogenous [67]. In order to assess individual variation in monthly activity space use, Kruskal-Wallis tests were also run separately for monthly KD estimates by residency status (transient, semi-transient, and resident) and plotted. A post-hoc pair-wise Wilcox test was run if differences were significant to determine where differences occurred. A student’s t-test was run to assess differences in activity spaces between grouper spawning months and non-spawning months.

Results

Spatial habitat use and residency

Between 22 and 26 January 2010 and 31 January to 16 March 2012, a total of 21 grey reef sharks were captured and fitted internally with an acoustic transmitter, which included four sharks inside Ant Channel and the remainder of sharks within the KMS (Table 1). All individuals were adult females as determined by a combination of size and the absence of claspers. Adult females ranged in size from 101 to 160 cm TL (mean = 147.0 cm TL). The smallest grey reef shark (101 cm TL) tagged at Ant (Transmitter 51334) was detected at the array for only a single day and was removed from the analysis. Among the 20 remaining individuals, total detection days ranged from 250–1347 d (Table 1).

thumbnail
Table 1. Tagging information, residency indices (RI), and Kernel Density (KD) estimates1.

https://doi.org/10.1371/journal.pone.0221589.t001

Residency indices at the array by daily occurrence were high for all sharks combined, with an average (±SD, hereafter) RI of 0.64±0.29 (Table 1). Of the 20 sharks used in the analysis, six individuals were considered transient, eight were semi-transient and six were highly resident (Table 1). All sharks tagged at Ant were considered semi-transients (Table 1).

The number of days sharks were detected within the KMS was similar to those spent outside the KMS (8,103 days inside vs. 7,439 outside; 64% vs. 71% of total days each shark was detected at the array). Two of the 20 sharks were never detected at the KMS, both of which were initially tagged at Ant. Residency within the KMS was highly variable, with an average overall RI of 0.48±0.36. Sharks that were considered transient at the array spent the least amount of time within the KMS (RI = 0.22±0.11), followed by semi-transient sharks (RI = 0.34±0.29) (Table 1). Of six sharks considered to be resident to the array, five had a daily RI ≥ 0.90±0.17 within the KMS (Table 1). The 18 sharks detected at the FSA throughout the study had an average daily RI of 0.47±0.35. The daily RI was not significantly different at the FSA between full moon and new moon periods (Student’s t-test, p = 0.37). Average monthly RI at the array for all sharks was not significantly different across all months (ANOVA, p>0.05), with an average of 0.80±0.30 (Fig 2). There were insignificant trends (ANOVA, p>0.05) in decreasing average monthly RI among semi-transient and transient sharks, while resident sharks maintained similarly high RI across all months (Fig 3A). Monthly average RI was significantly higher at the FSA during spawning months versus non-spawning months for all residency types (0.59±0.38 vs. 0.52±0.37, respectively; Student’s t-test, p<0.05).

thumbnail
Fig 2. Kernel density and residency index estimates.

Monthly 50% and 95% kernel density utilization (KD) and residency index (RI) estimates for all female grey reef sharks Carcharhinus amblyrhynchos monitored by acoustic array at Pohnpei and Ant Atoll from January 2010-October 2013. Error bars represent standard error. The outlined box represents the aggregation months.

https://doi.org/10.1371/journal.pone.0221589.g002

thumbnail
Fig 3. Monthly residency index and kernel density.

(A) Monthly residency index (number of days detected divided by the number of days in each month) and (B) Monthly KD for transient, semi-transient, and resident female grey reef sharks tagged at Pohnpei and Ant Atoll, Micronesia during two tagging events in 2010 and 2012. The black box outlines aggregation months (January–May), while the shaded area represents months when either two or all three groupers are present (February–April).

https://doi.org/10.1371/journal.pone.0221589.g003

Activity space

Large variations in activity space (KD) were observed among individuals. By residency type, semi-transient and transient sharks showed a much wider range of movement and activity space outside of the aggregation periods (Fig 3B). Activity spaces were small overall, with only one shark utilizing both Pohnpei and Ant Atoll regularly (Transmitter 48010, 95% KD = 109.0 km2). The extent of activity spaces of all other sharks was < 45 km2 (95% KD; Fig 4; Table 1) and most were under 20 km2. The semi-transient sharks had the largest average activity spaces (mean 50% KD = 2.1 km2, 95% KD = 20.6 km2), while the resident sharks had the smallest (Table 1). The six resident sharks at the array had an average activity space core of just 0.32 km2 (50% KD) and an extent of 2.1 km2 (95% KD). All resident sharks spent a high number of days within the KMS (Table 1).

thumbnail
Fig 4. Shark activity space.

Activity spaces for 20 female grey reef sharks calculated as 95% Kernel density (KD) utilization for A) resident sharks overall; B) resident sharks during grouper spawning months; C) semi-transient sharks overall; D) semi-transient sharks during grouper spawning months; E) transient sharks overall; and F) transient sharks during grouper spawning months.

https://doi.org/10.1371/journal.pone.0221589.g004

In general, activity spaces increased over the course of the year, although there were not significant differences detected among months (averaged by calendar month) (Kruskal-Wallis p > 0.05, Fig 2). Overall, activity spaces for resident sharks remained constant, while non-significant increases were observed among semi-transient and transient sharks for the core (50% KD) (Fig 3B) and extent (95% KD) of activity spaces. Average monthly 50% KD and 95% KD estimates were significantly lower during grouper spawning months than non-spawning months (2.5 vs. 7.3 km2 and 7.8 vs. 38.8 km2, respectively; Student’s t-test, p < 0.05) (Figs 2 and 4).

Results for mixed models indicated that month and grouper spawning period influenced the activity space of sharks at the study site (Tables 2 and 3). The best-fit model for core activity space (50% KD) (Table 2) included month and grouper spawning period, and the best-fit model for the extent of activity space (95% KD) (Table 3) included only month. All candidate models were significantly different than the null model (p < 0.001). Model results were similar between the 50% and 95% KD estimates, with the top three models for both scenarios including month, grouper spawning period, and residency status (Tables 2 and 3). AIC values were similar for the top candidate models, indicating weak support.

thumbnail
Table 2. The 15 top-ranked GLMMs for the 50% Kernel Density (KD) activity spaces.

https://doi.org/10.1371/journal.pone.0221589.t002

thumbnail
Table 3. The 15 top-ranked GLMMs for the 95% Kernel Density (KD) activity spaces.

https://doi.org/10.1371/journal.pone.0221589.t003

Results from the 15 top-ranked GLMMs for the 95% Kernel Density (KD) activity spaces of female grey reef sharks Carcharhinus amblyrhynchos at Pohpei and Ant Atoll from January 2010-October 2013. The best-fit model with the lowest AIC is shown in bold. Month: sequential month from the start of the study; Spawning: grouper spawning period (binary); Resident: residency status (transient, semi-transient, resident); Temp: temperature; TL: total length.

Model averaging of cumulative AIC weights supported the model selection, and parameters in the top candidate models for both 50% and 95% KD estimates were shown to influence (wi > 0.5) activity spaces of grey reef sharks tagged at Pohnpei (Table 4). Month was the most influential parameter for the core and extent of activity space, followed by grouper spawning period, with activity spaces contracting during these periods (Fig 4 and Table 4).

thumbnail
Table 4. Cumulative AIC weights (w), parameter estimates, standard error (SE), and p-values where deltaAIC < 2.

https://doi.org/10.1371/journal.pone.0221589.t004

During aggregation months, RI was lowest and most variable during January when only P. areolatus is present in the FSA site. Previous work at the FSA site suggests that it is during this period when P. areolatus aggregations are highly or exclusively composed of males.

Inter-reef and long-distance movement

In general, grey reef sharks tended to be highly resident to the island or atoll where they were captured, however six females tagged at Pohnpei (35%) were detected at Ant Atoll. Of the four sharks tagged at Ant Atoll, Transmitter 51334 was detected only for one day and subsequently removed from the analysis. Two females were never detected at any station at Pohnpei and the fourth (Transmitter 48010) was detected at Pohnpei for 39 days out of the 318 total detection days, with the remainder of detections at Ant.

Discussion

In Pohnpei, Micronesia, grouper spawning events appear to influence the movement patterns and activity spaces of semi-transient and transient female grey reef sharks. During the grouper spawning season, grey reef shark activity spaces contracted around FSA sites, daily detections increased and the number of tagged sharks present at the FSA rose, in support of prior research that demonstrated the use of FSA for foraging by elasmobranchs. While the scale of movement and activity space for all sharks extended well beyond MPA boundaries (95% KD range = 2.1–109.0 km2; average = 12.0±25.3 km2) during non-spawning months, most sharks tagged were detected within the KMS on more than 50% of monitored days. For resident sharks, the small activity space and extensive time within the MPA provides proof that small-scale MPAs have benefits to at least some individuals and offers partial benefits to others. These combined findings highlight that, although limited spatially and temporally, small-scale MPAs have utility in reducing the vulnerability of sharks and groupers to fishing during periods when they are concentrated. Conversely, the large-scale movements and associated activity spaces of sharks outside of grouper spawning times, and the use of common reproductive migratory corridors by groupers, supports a combined approach to management that limits fishing impacts outside of reproductive periods and away from MPAs.

Previous published reports have shown the direct trophic benefits to elasmobranchs from FSA foraging [24, 14, 2021, 2325, 35]. For example, spawning aggregations of E. polyphekadion in Fakarava, French Polynesia, and presumably elsewhere, create the conditions for an inverted biomass pyramid that allows grey reef (and possibly other) sharks to maximize foraging success and minimize energetic, long-distance movements during aggregation periods [14]. Although typically less abundant than Fakarava, E. polyphekadion FSA of similar densities [33, 47] are known throughout the Indo-Pacific [32, 78], many in association with the same (or congeneric) grouper species found in Pohnpei. While predation success rates were recorded as low in Fakarava, the concentration of prey nonetheless increased opportunities for foraging relative to non-reproductive periods [14]. The FSA also appeared to attract sharks from a substantially larger catchment area. Similar foraging events have been observed by whale sharks on cubera snapper (Lutjanus cyanopterus) (Cuvier 1828) FSA [4] and whale sharks and bull sharks (Carcharhinus leucas) (Muller & Henle 1839) on mutton snapper (Lutjanus analis) (Cuvier 1828) spawn and spawners, respectively [79]. Alfred mantas have been observed schooling and feeding on eggs produced within surgeonfish FSA [23]. For other reef and nearshore species, oophagy has been widely reported, while defecation of digested materials provide nutrients to a wide range of pelagic and benthic organisms [27]. Clearly, FSA serve as productivity hotspots and have important benefits to ecosystem health and trophodynamics at multiple scales [57]. These benefits are expected to increase as FSA biomass and reproductive output increase. In recognition of these food web linkages, FSA protection provides multiple benefits to ecosystems, including sharks and rays.

In Pohnpei, acoustic telemetry surveys clearly suggest that small-scale MPAs can provide ephemeral protection to female C. amblyrhynchos, with the level of protection varying by shark residency type. Specifically, resident sharks would be expected to receive greater benefits both spatially and temporally from these small set-asides, while for transient and semi-transient sharks, the protections more closely resemble that of aggregating grouper, given the extent of movement and use of common reproductive migratory corridors by these transient spawners. For aggregating groupers, at least some individuals have been detected or captured up to 27 kms from the FSA site [6264], with pathways typically along outer reef areas to reach home reefs. The scale of movement suggests that catchment areas for groupers are on the order of 100s of km2 [63], thus similar to non-resident sharks. For coastal sharks as well as aggregating grouper, large-scale no-entry or no-take management options appear to be more optimal solutions than small MPAs to long-term population persistence. However, the degree of area protection needed is both politically and socially untenable in PICTs where coastal fisheries represent important sources of protein and income, and fishing areas are inherently limited.

Although there is no formal enforcement, the Federated States of Micronesia provided blanket protection for sharks from fishing within the 3-million km2 economic exclusion zone (EEZ) in 2015 through the passage of Public Law No. 18–108. Owing to the limited reliance on sharks for either food or income, the law received little opposition. In addition to the FSM shark law, in 2017 Pohnpei State extended its own exclusion zone from 12 to 24 miles (38.6 km) (Public Law 19–167) preventing foreign commercial fishing from nearing its coastlines. Such large-scale MPAs could provide the protection required by these animals, however it is currently unclear how effective enforcement of such large areas can be accomplished.

For grouper FSA, Pohnpei modified its two-month (1 March– 30 April) grouper sales ban to extend protection throughout the spawning season (1 January to 31 May). The law was further amended to extend the sales ban to include catch and possession of groupers during this period. Unfortunately, the law provided a loophole allowing individual fishers to take 10 groupers per person per day, which greatly offsets any gains that would have otherwise been realized. Moreover, extended seasonal bans tend to exert greater fishing pressure on other vulnerable species when fishers’ primary goals are catch volume and not species composition [80]. In contrast, an inclusion of a catch ban to the grouper sales ban may lessen the attractiveness of FSA fishing and thus reduce vulnerability on both groupers and sharks at the aggregation site. Ultimately, an ecosystem-based management plan would benefit these various species and, in particular, be more suitable to both resident and transient sharks.

Sharks are often sex segregated during non-reproductive periods [51, 53, 59] and are reportedly impacted differentially by various types of fishing gear. In Pohnpei, fishing in relatively shallow water (c. 30 m) along outer barrier reef walls and slopes predominantly yielded only large females. The paucity of males and juveniles in this (and other) study somewhat limits the conclusions that can be made regarding species-level conservation for grey reef sharks using small-scale MPAs, however it also highlights the need to conduct more robust investigations integrating all life history stages. The studies conducted to date also reveal the complexity associated with shark conservation using small-scale area protection. Reef connectivity, shark development stage, reproduction, parturition and sex clearly have a bearing on reserve design, including placement, area and depth [81]. Future efforts to describe the conservation benefits of small- and large-scale MPAs, especially those focused on FSA areas should, if possible, incorporate all sizes and sexes of target animals utilizing a wider capture area, and use a range of habitats to enable more generalized statements on the design and spatial habitat required for effective nearshore shark conservation. Finally, information on grey shark reproductive periods and locations would benefit our understanding of spatial habitat associations and conservation needs.

Acknowledgments

The Pohnpei State Department of Lands and Natural Resources and Office of Fisheries and Aquaculture provided support throughout the study. Taylor Paul, Trace Loren, Scotty Malakai, and Quincy Malakai captured the animals. Rich Adams provided housing during tagging operations and Anthony Nolan provided goods and services throughout the project. This paper is dedicated to Zeddy Seymour.

References

  1. 1. Sadovy de Mitcheson Y, Colin PL (editors). Reef Fish Spawning Aggregations: Biology, research and management. Springer, Dortrecht. 2012.
  2. 2. Robertson DR. On the spawning behavior and spawning cycles of eight surgeonfishes (Acanthuridae) from the Indo-Pacific. Environ Biol Fish 1983; 9: 193–223.
  3. 3. Sancho G, Petersen CW, Lobel PS. Predator-prey relations at a spawning aggregation site of coral reef fishes. Mar Ecol Prog Ser 2000; 203: 275–288.
  4. 4. Heyman WD, Graham RT, Kjerfve B, Johannes, RE. Whale sharks Rhincodon typus aggregate to feed on fish spawn in Belize. Mar Ecol Prog Ser 2001; 215: 275–282.
  5. 5. Nemeth RS. Ecosystem aspects of species that aggregate to spawn. In: Sadovy de Mitcheson Y, Colin PL (editors). Reef fish spawning aggregations: Biology, research and management, Springer, Dortrecht. 2012. pp. 21–56.
  6. 6. Archer SK, Allgeier JE, Semmens BX, Heppell SA, Pattengill CV, Rosemond AD, Bush PG, McCoy CM, Johnson BC, Layman CA. Hot moments in spawning aggregations: implications for ecosystem-scale nutrient cycling. Coral Reefs 2015; 34: 19–23.
  7. 7. Erisman B, Heyman W, Kobara S, Ezer T, Pittman S, Aburto-Oropeza O, Nemeth R. Fish spawning aggregations: where well-placed management actions can yield big benefits for fisheries and conservation. Fish Fish 2015; 18: 128–144.
  8. 8. Hamilton RJ, Giningele M., Aswani S, Ecochard JL. Fishing in the dark- local knowledge, night spearfishing and spawning aggregations in the Western Solomon Islands. Biol Conserv 2012; 145: 246–257.
  9. 9. Heyman WD, Kjerfve B. Characterization of transient multi-species reef fish spawning aggregations at Gladden Spit. Bull Mar Sci 2008; 83: 531–551.
  10. 10. Rhodes KL, Nemeth RS, Kadison E, Joseph J. Spatial, temporal and environmental dynamics of a multi-species epinephelid spawning aggregation in Pohnpei, Micronesia. Coral Reefs 2014; 33: 765–775.
  11. 11. Warren-Rhodes K, Sadovy Y, Cesar H. Marine ecosystem appropriation in the Indo-Pacific: A case study of the life reef fish food trade. Ambio 2002; 32: 481–488.
  12. 12. Hamilton RJ, Matawai M. Live reef fish trade causes rapid declines in abundance of squaretail coralgrouper (Plectropomus areolatus) at a spawning aggregation in Manus, Papua New Guinea. SPC Live Reef Fish Information Bulletin 2006; 16: 13–18. https://coastfish.spc.int/en/publications/bulletins/live-reef-fish/97-live-reef-fish-information-bulletin-16
  13. 13. Erisman B, Aburto-Oropeza O, Gonzalez-Abraham C, Mascarenas-Osorio, Moreno-Baez M, Hastings PA. Spatio-temporal dynamics of a fish spawning aggregation and its fishery in the Gulf of California. Scientific Reports 2012; 2: 1–11. pmid:22359736
  14. 14. Mourier J, Maynard J, Parravicini V, Ballesta L, Clua E, Domeier ML, et al. Extreme inverted trophic pyramid of reef sharks supported by spawning groupers. Curr Biol 2016; 26: 1–6.
  15. 15. Sadovy Y, Domeier M. Are aggregation-fisheries sustainable? Reef fish fisheries as a case study. Coral Reefs 2005; 24: 254–262.
  16. 16. Olsen DA, LaPlace JA. A study of a Virgin Islands grouper fishery based on a breeding aggregation. Proc Gulf Caribb Fish Inst 1978; 31: 130–144.
  17. 17. Sadovy de Mitcheson Y, Craig MT, Bertoncini AA, Carpenter KE, Cheung WWL, Choat JH, et al. Fishing groupers towards extinction: A global assessment of threats and extinction risks in a billion dollar fishery. Fish Fish 2013; 14: 119–136.
  18. 18. Aguilar-Perera A. Disappearance of a Nassau grouper spawning aggregation off the southern Mexican coast. Mar Ecol Prog Ser 2006; 327: 289–296.
  19. 19. Graham RT, Roberts CM, Smart, JCR. Diving behaviour of whale sharks in relation to a predictable food pulse. J Royal Soc Interface 2005; 3: 109–116.
  20. 20. Hoffmayer ER, Franks JS, Driggers WB III, Oswald KJ, Quattro JM. Observations of a feeding aggregation of whale sharks, Rhincodon typus, in the north central Gulf of Mexico. Gulf Caribb Res 2007; 19: 69–73.
  21. 21. de la Parra Venegas R, Hueter R, Gonzaléz Cano J, Tyminski J, Remolina JG, Maslanka M, et al. An unprecedented aggregation of whale sharks, Rhincodon typus, in Mexican coastal waters of the Caribbean Sea. PLoS ONE 2011; 6(4): e18994. pmid:21559508
  22. 22. Robinson DP, Jaidah MY, Jabado RW, Lee-Brooks K, Nour El-Din NM, et al. Whale sharks, Rhincodon typus, aggregate around offshore platforms in Qatari waters. PLoS ONE 2013; 8(3): e58255. pmid:23516456
  23. 23. Hartup JA, Marshell A, Kottemair M. Carlson P. Manta alfredi multispecies surgeonfish spawning aggregations. Coral Reefs 2013; 32: 367.
  24. 24. Pickard AE, Vaudo JJ, Wetherbee BM, Nemeth RS, Blondeau JB, Kadison EA, et al. Comparative use of a Caribbean mesophotic coral ecosystem and association with fish spawning aggregations by three specie of shark. PLoS ONE 2016; 11(5): e011221. pmid:27144275
  25. 25. Etpison MT, Colin PL. Blue water spawning by Moorish idols and orangespine unicornfish in Palau: Is it a “suicide mission”? Aqua 2018; 23: 121–136.
  26. 26. Weideli OC, Mourier J, Planes S. A massive surgeonfish aggregation creates a unique opportunity for reef sharks. Coral Reefs 2015; 34: 835.
  27. 27. Laffoley D, Baxter J, Lebebvre C, Sévin M-A, Simard F. Building MPA networks by 2020: IMPAC3 achievements, future challenges and next steps. Anim Conserv Mar Freshw Ecosyst 2014; 24 (Suppl. 2): 238–245.
  28. 28. Beets J, Friedlander A. Evaluation of a conservation strategy: A spawning aggregation closure for red hind, Epinephelus guttatus, in the U.S. Virgin Islands. Environ Biol Fish 1999; 55:91–98.
  29. 29. Starr RM, Sala E, Ballesteros E, Zabala M. Spatial dynamics of the Nassau grouper Epinephelus striatus in a Caribbean atoll. Mar Ecol Prog Ser 2007; 343: 239–249.
  30. 30. Whaylen L, Pattengill-Semmens CV, Semmens BX, Bush PG, Boardman MR. Observations of a Nassau grouper, Epinephelus striatus, spawning aggregation site in Little Cayman, Cayman Islands, including multi-species spawning information. Environ Biol Fish 2004; 70: 305–313.
  31. 31. Nemeth RS, Blondeau J, Herzlieb S, Kadison E. Spatial and temporal patterns of movement and migration at spawning aggregations of red hind, Epinephelus guttatus, in the U.S. Virgin Islands. Environ Biol Fish 2007; 78: 365–381.
  32. 32. Waldie PA, Almany GR, Sinclair-Taylor TH, Hamilton RJ, Potuku T, Priest MA, et al. Restricted grouper reproductive migrations support community-based management. Roy Soc Open Sci 2016; 3: 150694; Available from: http://dx.doi.org/10.1098/rsos.150694.52.
  33. 33. Bijoux J, Dagorn L, Cowley PD, Simier M, Adam P-A, Robinson J. Spawning aggregation dynamics of brown-marbled grouper and camouflage grouper at a remote Indian Ocean atoll. Endanger Species Res 2013; 22: 145–157.
  34. 34. Green A, Maypa A, Almany G, Rhodes K, Weeks R, Mumby PJ, et al. Larval dispersal and movement patterns of coral reef fishes, and implications for marine reserve network design. Biol Rev 2014; 90: 1215–1247. pmid:25423947
  35. 35. Espinoza M, Cappo M, Heupel MR, Tobin AJ, Simpfendorfer CA. Quantifying Shark Distribution Patterns and Species-Habitat Associations: Implications of Marine Park Zoning. PLoS ONE 2014; 9, Available from: https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0106885
  36. 36. Barnett A, Abrantes KG, Seymour J, Fitzpatrick R. Residency and spatial use by reef sharks of an isolated seamount and its implications for conservation. PLoS ONE 2012; 7: Available from: https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0036574
  37. 37. Vianna GMS, Meekan MG, Meeuwig JJ, Speed CW. Environmental influences on patterns of vertical movement and site fidelity of grey reef sharks (Carcharinus amblyrhynchos) at aggregation sites. PLoS ONE 2012; 8: Available from:
  38. 38. McKibben JN, Nelson DR. Patterns of movement and grouping of grey reef sharks, Carcharhinus amblyrhynchos, at Enewetak, Marshall Islands. Bull Mar Sci 1986; 38: 89–110.
  39. 39. Heupel MR, Simpendorfer CA, Fitzpatrick R. Large-scale movement and reef fidelity of grey reef sharks. PLoS ONE 2010; 5: 3. Available from: https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0009650
  40. 40. Espinoza M, Heupel MR, Tobin AJ, Simpendorfer CA. Residency patterns and movements of grey reef sharks (Carcharhinus amblyrhynchos) in semi-isolated coral reef habitats. Mar Biol 2015; 162: 343–358.
  41. 41. White TD, Carlisle AB, Kroodsma DA, Block BA, Casagrandi R, De Leo GA, et al. Assessing the effectiveness of a large marine protected area for reef shark conservation. Biol Conserv 2017; 207: 64–71.
  42. 42. Houk P, Camacho R, Johnson S, McLean M, Maxin S, Anson J, et al. The Micronesia Challenge: Assessing the relative contribution of stressors on coral reefs to facilitate science-to-management feedback. PLoS ONE 2015; 10(6): e0130823. pmid:26087252
  43. 43. Graham NAJ, Spalding MD, Sheppard CRC. Reef shark declines in remote atolls highlight the need for multi-faceted action. Aquat Conserv: Mar Freshw Ecosyst 2010; 20:543–548.
  44. 44. Robbins WD, Hisano M, Connolly SR, Choat HJ. Ongoing collapse of coral-reef shark populations. Curr Biol 2006; 16: 2314–2319. pmid:17141612
  45. 45. Heupel MR, Williams AJ, Welch DJ, Ballaugh A, Mapstone BD, Carlos G et al. Effects of fishing on tropical reef associated shark populations on the Great Barrier Reef. Fish Res 2009; 95: 350–361.
  46. 46. Hamilton RJ, Potuku T, Montambault LR. Community-based conservation results in the recovery of reef fish spawning aggregations in the Coral Triangle. Biol Conserv 2011; 144: 1850–1858.
  47. 47. Nemeth RS. Population characteristics of a recovering US Virgin Islands red hind spawning aggregation following protection. Mar Ecol Prog Ser 2005; 286: 81–97. pmid:16612415
  48. 48. Smale M.J. 2009. Carcharhinus amblyrhynchos. The IUCN Red List of Threatened Species 2009: e.T39365A10216946. Available from: http://dx.doi.org/10.2305/IUCN.UK.2009-2.RLTS.T39365A10216946.en. Downloaded on 09 December 2017
  49. 49. Anderson RC, Sheppard C, Spalding M, Crosby R. Shortage of sharks at Chagos. Newsletter of the IUCN/SSC Shark Specialist Group, UK; 1998. Available from: https://www.researchgate.net/profile/Mark_Spalding2/publication/263705420_Shortage_of_Sharks_at_Chagos/links/54c8dd9a0cf289f0ced10570/Shortage-of-Sharks-at-Chagos.pdf
  50. 50. Mucientes GR, Queiroz N, Sousa LL, Tarroso P, Sims DW. Sexual segregation of pelagic sharks and the potential threat from fisheries. Biol Lett 2009; 5: 156–159. pmid:19324655
  51. 51. Beerkircher LR, Cortes E, Shivji M. Characteristics of shark bycatch observed in pelagic longlines off the southeastern Unites States, 1992–2000. Mar Fish Rev 2002; 64: 40–49.
  52. 52. Sims DW. Differences in habitat selection and reproductive strategies of male and female sharks. In: Ruckstuhl KE, Neuhaus P, editors. Segregation in vertebrates: Ecology of two sexes, Cambridge University Press; 2005, pp. 127–147.
  53. 53. Bond ME, Babcock EA, Pikitch EK, Abercrombie DL, Lamb NF, Chapman DD. Reef sharks exhibit site-fidelity and higher relative abundance in marine reserves on the Mesoamerican Barrier Reef. PLoSONE 2012; 7: Available from: https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0032983
  54. 54. Heuter RE, Heupel MR. Evidence of philopatry in sharks and implications for the management of shark fisheries. J Northw Atl Fish Sci 2004; 35: 239–247.
  55. 55. Hueter RE, Huepel MR, Heist EJ, Keeney DB. Evidence of philopatry and implications for the management of shark fisheries. J. Northw. Alt. Fish. Sci. 2004; 35: 239–247.
  56. 56. Strasburg DW. 1958. Distribution, abundance, and habits of pelagic sharks in the central Pacific Ocean. Washington, US Govt. Print Office, Fishery Bulletin 138. Available from: https://spo.nmfs.noaa.gov/sites/default/files/pdf-content/fish-bull/strasburg_0.pdf
  57. 57. Jacoby DMP, Busawon DS, Sims DW. Sex and social networking: the influence of male presence on social structure of female shark groups. Behav Ecol 2010; 2: 808–818. Available from: https://doi.org/10.1093/beheco/arq061
  58. 58. Smith SE, Au DW, Show C. Intrinsic rebound potentials of 26 species of Pacific sharks. Mar Freshw Res 1998; 49: 663–678.
  59. 59. Wetherbee BM, Crow GL, Lowe CG. Distribution, reproduction and diet of the gray reef shark, Carcharhinus amblyrhynchos in Hawaii. Mar Ecol Prog Ser 1997; 151: 181–189.
  60. 60. Rhodes KL, Taylor BM, McIlwain JL. Detailed demographic analysis of an Epinephelus polyphekadion spawning aggregation and fishery. Mar Ecol Prog Ser 2011; 421: 183–198.
  61. 61. Rhodes KL, Sadovy Y. Temporal and spatial trends in spawning aggregations of camouflage grouper, Epinephelus polyphekadion in Pohnpei, Micronesia. Environ Biol Fish 2002: 63: 27–39.
  62. 62. Rhodes KL, Tupper MH. The vulnerability of reproductively active squaretail coralgrouper (Plectropomus areolatus) to fishing. Fish Bull 2008; 106: 194–203.
  63. 63. Rhodes KL, McIlwain J, Joseph E, Nemeth RS. Reproductive movement, residency and fisheries vulnerability of brown-marbled grouper, Epinephelus fuscoguttatus (Forsskål, 1775). Coral Reefs 2012; 31: 443–453.
  64. 64. Hutchinson N, Rhodes KL. Spatial management needs of squaretail coralgrouper, Plectropomus areolatus (Rüppell 1830). Coral Reefs 2010; 29, 567–576.
  65. 65. Francis MP. Morphometric minefields–towards a measurement standard for chondrichthyan fishes. Environ Biol Fish 2007; 77: 407–421.
  66. 66. Simpfendorfer CA, Yeiser BG, Poulakis GR, Stevens PW, Heupel MR. Environmental influences on the spatial ecology of juvenile smalltooth sawfish Pristis pectinata: Results from acoustic monitoring. PLoS ONE 2011; 6: 2. Available from: https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0016918
  67. 67. Zar JH. Biostatistical analysis. Pearson Education India; 1999.
  68. 68. Lararidis E. lunar: Lunar, Phase & Distance, Seasons and Other Environmental Factors. 2014; (Version 0.1–04): Available from: http://statistics.lazaridis.eu
  69. 69. Calenge C. The package adehabitat for the R software: a tool for the analysis of space and habitat use by animals. Ecol Modell 2006; 197: 516–519.
  70. 70. Worton B. Kernel methods for estimating the utilization distribution in home-range studies. Ecol 1989; 70: 164–168.
  71. 71. R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2008: ISBN 3-900051-07-0, Available from: http://www.R-project.org.
  72. 72. Pebesma EJ, Bivand RS. Classes and methods for spatial data in R. R News 2005; 5: 2. Available from: https://cran.r-project.org/doc/Rnews/.
  73. 73. Zuur AF, Ieno EN, Walker N, et al. Mixed effects models and extensions in ecology with R: New York; Springer; 2009.
  74. 74. Heupel MR, Simpfendorfer CA. Importance of environmental and biological drivers in the presence and space use of a reef-associated shark. Mar Ecol Prog Ser 2014; 496: 47–57.
  75. 75. Akaike H. Information theory and an extension of the maximum likelihood principle. In: Proceedings of the 2nd International Symposium on Information Theory, Tsahkadsor, Armenia SSR; 1973. pp. 267–281.
  76. 76. Barton K. MuMIn: multi-model inference. 2013. R package version 1.9. 5. Available from: http://CRAN.R-project.org/package=Mu Min
  77. 77. Schlaff AM, Meupel MR, Udyawer V, Simpendorfer CA. Biological and environmental effects on activity space of a common reef shark on an inshore reef. Mar Ecol Prog Ser 2017; 571: 169–181.
  78. 78. Johannes RE, Squire L, Graham T, Sadovy Y, Renguul H (1999) Spawning aggregations of groupers (Serranidae) in Palau. The Nature Conservancy Marine Research Series Publication No. 1, August 1999, 144 p. Available from: http://www.conservationgateway.org/Files/Pages/spawning-aggregations-gro.aspx
  79. 79. Graham RT, Castellanos DW. Apex predators target mutton snapper spawning aggregation. Coral Reefs 2012; 31: 1017.
  80. 80. Rhodes KL, Tupper MH, Wichilmel CB. Characterization and management of the commercial sector of the Pohnpei coral reef fishery, Micronesia. Coral Reefs 2008; 27:443–454.
  81. 81. Speed CW, Field IC, Meekan MG, Bradshaw CJA. Complexities of coastal shark movements and their implications for management. Mar Ecol Prog Ser 2010; 408: 275–293.