Spawning aggregations of checkered snapper (Lutjanus decussatus) and blackspot snapper (L. fulviflamma): seasonality, lunar-phase periodicity and spatial distribution within spawning ground

Snappers (family Lutjanidae) are important fisheries target species and some species are known to form spawning aggregations at particular spawning grounds. The present study investigated the ecological characteristics of fish aggregations of two snapper species (checkered snapper Lutjanus decussatus and blackspot snapper L. fulviflamma) that form at a particular site. Specifically, the aims were to clarify (1) seasonality and lunar-phase periodicity of fish aggregation formation, (2) fine-scale spatial distribution of fish density (spatial variations of fish density at intervals of several-tens meters) within the aggregation site, (3) size and age frequency distributions of fishes in the aggregation site, (4) gonad development, (5) to compare fish abundance between inside and outside the aggregation site, and (6) to verify that fish aggregations of the two snapper species were spawning aggregation. Underwater observations using a 600 m × 5 m transect revealed that greater fish abundance of Lutjanus decussatus was found monthly between May and October, and clear positive peaks in the fish abundance were found only around the last-quarter moon. This lunar-related periodicity in the increase of fish abundance was confirmed by a time-series analysis (correlogram). Within the aggregation site, L. decussatus showed a relatively uniform distribution. In contrast, greater fish abundance of L. fulviflamma was found monthly between April and October, and clear positive peaks in the fish abundance were found around the last-quarter moon (April, May, June and October) or new moon (July, August and September). This lunar-related periodicity was also confirmed by correlogram. Lutjanus fulviflamma showed a relatively clumped distribution within the aggregation site. Most females of the two species in the aggregation site had hydrated eggs, indicating that the two species form aggregations for reproduction. The two species, although occurring simultaneously, are considered to form aggregations of conspecifics only. For L. decussatus, average fork length and age of males and females were 229.2 mm and 243.9 mm and 9.4 years and 8.1 years, respectively. For L. fulviflamma, average fork length and age of males and females were 233.9 mm and 246.9 mm and 6.8 years and 8.1 years, respectively. Fish abundance inside the aggregation site was 266.8-fold and 141557.1-fold greater than those outside the aggregation site for L. decussatus and L. fulviflamma, respectively. These results showed that (1) fish aggregation formation of the two snapper species was predictably repeated in particular months and lunar-phase, (2) it was predictably found at the particular site, (3) the fish abundance in the aggregation site markedly exceeded the fish abundance outside the aggregation site, and (4) the two species form aggregations for reproduction. Therefore, it is suggested that the fish aggregations for the two species can be regarded as spawning aggregations.

of the two snapper species was predictably repeated in particular months and lunarphase, (2) it was predictably found at the particular site, (3) the fish abundance in the aggregation site markedly exceeded the fish abundance outside the aggregation site, and (4) the two species form aggregations for reproduction.Therefore, it is suggested that the fish aggregations for the two species can be regarded as spawning aggregations.

INTRODUCTION
Coral reef fishes show diverse reproductive behaviors (Thresher, 1984).Among these diverse fish species, some fish species form aggregations with greater densities during restricted seasons and lunar phases at particular spawning grounds (Nemeth, 2009).These fish aggregations are known as spawning aggregations (Sadovy de Mitcheson & Colin, 2012).Domeier (2012) defined spawning aggregations as conspecific individuals gather at specific sites in a specific period.Over 80 fish species are regarded to form spawning aggregations (Sadovy de Mitcheson & Colin, 2012), although it is suggested that the actual number of species would be greater than 160 (Claydon, 2004).Two types of spawning aggregations have been reported: resident and transient spawning aggregations (Domeier & Colin, 1997).Resident spawning aggregations are characterized by smaller-sized species (e.g., parrotfishes, surgeonfishes and wrasses), shorter migration distance (within several kilometer scale), shorter duration of the spawning event (several hours) and almost daily aggregation.On the other hand, transient spawning aggregations are characterized by larger-sized species (e.g., emperorfishes, groupers and snappers), longer migration distances (scales of several to several-hundreds of kilometers), longer duration of the spawning event (several days to several weeks) and almost monthly or annual aggregations (Nemeth, 2009).
Two snapper species, checkered snapper Lutjanus decussatus (Cuvier, 1828) and blackspot snapper L. fulviflamma (Forsskål, 1775), are important fisheries target species in Okinawan coral reefs (Akita et al., 2016).Nanami et al. (2010) indicated the possibility that L. decussatus forms spawning aggregations in this region by clarifying its reproductive biology.The main spawning season was between June and October and clear lunarsynchronized fluctuations in the gonadsomatic index (GSI) of females were found.The highest GSI values were found around the last-quarter moon phase during the main spawning season.However, direct evidence (fish aggregation at a particular site) has not yet been reported.In contrast, Shimose & Nanami (2015) showed that the main spawning season of L. fulviflamma occurred between April and August.However, lunarsynchronized reproductive activity as well as the possibility of spawning aggregation formation of L. fulviflamma have not been examined yet.Although there has been no ecological information from the local communities about spawning aggregations of the two species, fish aggregations of the two snapper species were found in an Okinawan coral reef (Fig. 1, Videos S1 and S2).This finding suggests that the fish aggregation might be a spawning aggregation and should be appropriately protected since these two snapper species are target species of commercial catch and recreational fishing in the region.Thus, understanding precise ecological characteristics of fish aggregations of the two snapper species would be useful for effective management including the necessary spatial scale and duration of spawning ground protection.
The aims of the present study were to clarify ecological aspects of the fish aggregations of Lutjanus decussatus and L. fulviflamma.Specifically, the aims were to clarify: (1) seasonality and lunar-phase periodicity of fish aggregation formation, (2) fine-scale spatial distribution of fish density (spatial variations of fish density at intervals of several-tens of meters) at the aggregation site, (3) size and age frequency distributions of fishes in the aggregation site, (4) gonad development of fish individuals that were captured in the aggregation site, (5) to compare fish abundance between inside and outside the aggregation site, and (6) to verify that fish aggregations of the two snapper species were spawning aggregation.Since the present study is the first examination about fish aggregation of the two snapper species, the results would contribute to a more comprehensive understanding of ecological aspects about spawning aggregations of snappers in coral reefs.

MATERIALS AND METHODS
Main method of this study was field observation.Some fish individuals were caught as samples to examine the gonad development.By placing on ice, these fish individuals were immediately killed to minimize pain.Okinawa prefectural government fisheries coordination regulation No. 37 approved the sampling procedure.This regulation permits capture of fishes for scientific purposes on Okinawan coral reefs.

Study site
This study was conducted at Sekisei lagoon in the Yaeyama Islands, Okinawa, located in the southernmost part of Japan (Fig. 1).One fish aggregation site was located in the Sekisei lagoon, although the precise location is not shown.This is because the study site is not protected during the spawning periods of the two snapper species yet and there is concern that showing the precise location might lead to overfishing the fish aggregations of the two snapper species.

Temporal variations in the fish abundance within the aggregation site
To clarify the monthly and weekly variation in the fish abundance at the aggregation site, daytime underwater observations were conducted.After preliminary surveys between September and October 2020, one line transect (600 m ×5 m) was set to extend over the whole aggregation site between November 2020 and December 2021.Underwater observations on the transect were conducted using SCUBA.During the main spawning season of the two snapper species (between May and October), weekly observations were conducted (Table 1).The observation days were adjusted to be carried out during each of the four lunar phases whenever possible (new moon, first-quarter moon, full moon and last-quarter moon).When rough sea conditions prevented such adjustments, the observation days were set as several days earlier or later of the four lunar phases (Table 1).During non-spawning seasons (between November 2020 and April 2021, and between November and December 2021), underwater observations were conducted around two lunar phases (last-quarter moon and new moon).This was because: (1) preliminary surveys between September and October 2020 revealed the fish aggregations of the two snapper species were found around the last-quarter moon and new moon and (2) Nanami et al. (2010) showed that Lutjanus decussatus showed a greater peak of ovary development during the last-quarter moon.
The number of individuals of the two snapper species on the above-mentioned 600 m ×5 m line transect was counted per every 1 min.During the observations, a portable GPS receiver, sealed in a waterproof case and attached to a buoy, was towed and the course and distance of the tracks were obtained in accordance with the protocol detailed in Nanami et al. (2017).The water depth range at which the underwater observations were conducted was c.a. 10 m to 15 m.

Analysis of periodicity in the increase of fish abundance in relation to lunar phase
To determine the statistical significance in periodicity of fish increasing at the aggregation site, a time-series analysis (correlogram) was applied.In this procedure, data obtained between May 3 and October 20 was used for the analysis (Table S1), since greater fish abundance at the aggregation site was found in the study period (see RESULTS).In the analysis, 24 time periods were established, in which each time period includes one lunar phase (last-quarter moon, new moon, first-quarter moon and full moon) (Table S1).
Nanami (2023), PeerJ, DOI 10.7717/peerj.159916/26 was determined by using 95% confidence intervals (CI) as follow: where T is the number of observations.It can be regarded that the auto-correlation coefficient is significant when the following equation was found: where p k is the value of auto-correlation at k th time lag, and | p k | is the absolute value of p k .If |p k | is greater than the 95% CI, auto-correlation at k th time lag was significant.
For example, auto-correlation coefficient at 4th time lag (p 4 ) was over 0.40, the p 4 is a significantly positive value.Namely, the periodic increase in fish abundance was found at the same particular lunar phase.Since data for July 22 (full moon phase) could not be collected due to a typhoon, the missing value was imputed.In this procedure, other fish abundance data obtained around full moon phase (May 26, June 25, August 24, September 21 and October 20) was averaged and the average value was applied to impute the missing value.

Fine-scale spatial variations in fish density at the aggregation site
The above-mentioned underwater observations revealed that fish abundance showed positive peaks around last-quarter moon or new moon (see RESULTS).At the two lunar phases, fine-scale spatial variations in fish abundance (variations in fish abundance at intervals of several-tens meters) within the aggregation site were examined.The abovementioned 600 m × 5 m line transect was divided into 1-minute sub-transects (average distance ± standard deviation = 20.5 m ± 4.1).Then, the number of individuals of the two snapper species and the distance for the 1-minute sub-transect were obtained.From the data, the number of individuals was converted to density (per 20 m × 5 m) for each 1-minute transect.Fish density on the 1-minute sub-transect was individually plotted by bubble plot along the 600 m × 5 m line transect.
The fine-scale spatial variation in fish density per 20 m × 5 m was shown as frequency data (histogram).Then, Kolmogorov-Smirnov test was applied to test the significant difference in fine-scale spatial variation in fish density between the two fish species.Domeier (2012) and de Mitcheson & Erisman (2012) have indicated that presence of females with matured eggs in the hydrated stage (hydrated eggs) in aggregation sites is one of the direct evidence that the fish aggregation can be regarded as a spawning aggregation.To clarify whether females have hydrated eggs in the aggregation site, individuals of the two species were caught by spear-gun just after the above-mentioned underwater observations around the last-quarter and new moon (Table 1).In total, 73 and 92 individuals were caught for L. decussatus and L. fulviflamma, respectively (Table 1).In the laboratory, fork length (FL), whole body weight and gonad weight were measured.The gonadsomatic index (GSI) was calculated by using the formula: GSI = Gonad weight (g)/[whole body weight (g) − gonad weight (g)] × 100.

Verification of spawning by ovarian development
To obtain histological observations, the gonads were preserved in 20% buffered formalin over 48 h and then kept in 70% ethanol baths.Embedded pieces of gonads were sectioned and stained with Mayer's hematoxylin-eosin.Under microscopic observations, developmental stages of ovaries were examined whether the gonads were sufficiently developed for spawning (hydrated eggs).The categorization of ovarian developmental stages followed Ohta & Ebisawa (2015) and Ohta et al. (2017).According to the categorization of Ohta & Ebisawa (2015), oocytes with migration nuclear stage, pre-maturation stage and maturation stage were defined as hydrated stage.

Size and age frequency distributions
To determine the size frequency distribution of fishes at the aggregation site, histograms of fork length for the above-mentioned fish samples was plotted.Male and female were separately plotted, and probability density of size frequency was analyzed by R statistical computing language (function ''density'': R Core Team, 2022).
To determine the age frequency distribution of fishes, sagittal otoliths of the abovementioned fish samples were extracted from each fish, and cleaned in water and dried.Then, one otolith was embedded in epoxy resin and transversely sectioned into 0.5 mm-thick sections using ISOMET low speed saw (Buehler) and attached on a glass slide.The sectioned otoliths were observed under a microscope with reflected light at ×4 magnification, and the number of opaque rings was counted.The procedure of opaque rings count followed Ohta et al. (2017).The number of opaque rings on each otolith was counted twice.If the two counts coincided, the number of rings was used.However, if the two counts did not coincide, the number of opaque rings was counted once more and any two coinciding counts were used.Since Nanami (2021) and Shimose & Nanami (2015) revealed that each opaque ring was formed annually for the two species, number of opaque rings can be considered as age (year).Age frequency distribution for male and female was separately plotted, and probability density was analyzed by R statistical computing language (function ''density'': R Core Team, 2022).
Mann-Whitney U-test test was applied to determine the significant differences in average fork length and average age between males and females.Domeier (2012) has proposed the definition of spawning aggregations, indicating that fish abundance in aggregation sites is at least 4-fold greater than that outside the aggregation site.In order to verify the definitions, 65 study sites outside the aggregation site were established in Sekisei Lagoon (Fig. S2) and the number of the two snapper species was counted.Above-mentioned 20-minutes underwater survey with a portable GPS receiver was conducted at each site (details about method was shown in Nanami (2020).Fish abundance per 600 m ×5 m was estimated by using the fish count data and the measured distance.The estimated fish abundance per 600 m ×5 m among the 65 sites were averaged and regarded as average fish abundance outside the aggregation site.Then, the fish abundance per 600 m ×5 m inside the aggregation site was compared with that outside aggregation site.

RESULTS
Temporal variations in the fish abundance and reproductive activity in relation to lunar phase Greater fish abundance (number of individuals per 600 m ×5 m) of Lutjanus decussatus was found between May and October (Fig. 2A), when the water temperature exceeded 25 • C.During the six months, clear peaks in the fish abundance were only found around the last-quarter moon (Fig. 2B).The peak fish abundance ranged from 529 (October) to 1565 (June).Average GSI values ± SD (standard deviation) of females that were caught at the aggregation site was 10.34 ± 4.95 (ranged from 0.42 to 24.89) (Fig. 2C).About 69% of individuals (34 out of 49) had hydrated stage oocytes (Figs.2C-2E, Table S3).
Greater fish abundance of L. fulviflamma was clearly found between April and October (Fig. 3A), when water the temperature exceeded 25 • C.During the seven months, clear peaks in the fish abundance were found around the last-quarter moon in April, May, June and October (Fig. 3B).In contrast, clear peaks were found around the new moon in July, August and September (Fig. 3B).The peak fish abundance ranged from 660 (September) to 6337 (June).Average GSI values ± SD of females were 6.33 ± 2.21 (ranged from 3.44 to 14.51) (Fig. 3C).About 90% of individuals (37 out of 41) had hydrated stage oocytes (Figs.3C-3E, Table S4).

Periodicity of the increase of fish density
Correlogram of L. decussatus revealed that significant positive auto-correlation coefficients were found when time lags were 4 and 8 (Fig. 4A).Correlogram of L. fulviflamma revealed that significant positive auto-correlation coefficient was found when time lag was 4 (Fig. 4B).

Fine-scale spatial variations in fish density at the aggregation site
Fine-scale spatial distributions revealed that L. decussatus showed relatively uniform distribution at the aggregation site, especially around the last-quarter moon between May and October (Fig. 5).Fish densities per 20 m ×5 m were generally less than 140 individuals in most cases (Fig. 6).
In contrast, greater fish density of L. fulviflamma per 20 m ×5 m was found at particular areas (south-eastern side) in the aggregation site (Fig. 7).This tendency was clearly found around last quarter moon in April, May and June, and around new moon in July and August.In some cases, fish density within 20 m ×5 m was over 300 (Fig. 8).
Kolmogorov-Smirnov test revealed that a significant difference in fish density frequency (i.e., patterns in fine-scale spatial distribution) was found between the two fish species from May to October (p < 0.05, Table 2)., 9D).For L. fulviflamma, fork length of most individuals was less than 250 mm for males whereas 280 mm for females (Figs.10A, 10B).Average fork length of females (246.9 mm FL ± 23.0 SD) was significantly greater than males (233.9 mm FL ± 17.8 SD: Mann-Whitney

Comparison of fish density between inside and outside the aggregation site
Fish abundance per 600 m ×5 m of L. decussatus inside the aggregation site was 148.2 to 438.4-fold (average = 266.8-fold)greater than that outside the aggregation site (Table 3).Fish abundance of L. fulviflamma inside the aggregation site was 33000.0 to 316850.0-fold(average = 141557.1-fold)greater than that outside the aggregation site (Table 3).

Verification of spawning aggregation of two snapper species
For the two snapper species Lutjanus decussatus and L. fulviflamma, this study was the first attempt to examine whether fish aggregations at a particular site can be regarded as spawning aggregations.The present study showed: (1) repeated greater fish abundance of the two species and this is predictable in time (particular month and lunar-phase) and space, (2) the fish abundance inside the aggregation site is over 4-fold greater than that outside the aggregation site, and (3) most females inside the aggregation site had hydrated eggs.Thus, the fish aggregations of the two snapper species can be regarded as spawning aggregations.Graphs at upper right side show the temporal changes in number of fish individuals (see Fig. 3A) and red arrows represent the day when the data was collected.Lunar phases are abbreviated as LQM, last-quarter moon; NewM, New moon.''+'' and ''-'' mean after and before the lunar phase (e.g., ''NewM-1'' means 1 day before new moon).

Spawning season and spawning day
The present study revealed that fish aggregations of the species were found between May and October.Nanami et al. (2010) showed similar results from fish samples by commercial catch, i.e., estimated main spawning season of L. decussatus was between June and October.
Since developed-or matured-oocytes (tertiary yolk stage, migration nuclear stage, prematuration stage and maturation stage) were obtained for most female individuals between May and October, it is suggested that the main spawning season of L. decussatus can be regarded as between May and October in the study region.Since peaks of fish abundance were found only around the last-quarter moon and the periodicity in the increase of fish abundance during each last-quarter moon phase was significant, it is suggested that spawning occurs around the last-quarter moon.
The present study revealed that spawning aggregations of the species were found and the most females had hydrated eggs between April and October.Shimose & Nanami (2015) showed similar results from fish samples by commercial catch, i.e., estimated main spawning season of L. fulviflamma is between April and August.Therefore, it is suggested that the spawning season of L. fulviflamma can be regarded as between April and October in the study region.The peak fish abundance at the spawning ground was found during the last-quarter moon (in April, May, June and October) and new moon (in July, August and September), suggesting that spawning occurs around the last-quarter moon and new moon.This lunar-related periodicity in the increase of fish abundance at during a particular lunar phase (last-quarter moon or new moon) was also supported by the correlogram in the present study.This evidence for lunar-related spawning is the first finding for L. fulviflamma.
Since western Atlantic snapper species (Lutjanus cyanopterus (Cuvier, 1828) and L. jocu (Bloch & Schneider, 1801)) spawned for a period of four to seven consecutive days (Heyman & Kjerfve, 2008), the two snapper species in this study might spawn during several consecutive days around last-quarter moon and new moon.Intensive daily observations around the two lunar phases would clarify more precise ecological aspects about aggregation formation of the two snapper species.
Water temperature is one of the main factors for gonad development and reproduction of marine fishes (Wang et al., 2010).The present study revealed that reproductive activity of the two snapper species is likely to occur when water temperature is over 25 • C. In addition, decline of temperature (beginning of October) might be one of the factors leading to reduced ovarian development.This trend is consistent with cubera snapper (Lutjanus cyanopterus) in the Caribbean (Heyman et al., 2005;Motta et al., 2022), showing that spawning aggregation of L. cyanopterus related to increasing water temperature in the summer.Notes.

Fine-scale spatial variation in fish density within spawning ground
Some previous studies have shown fine-scale spatial variations of fish density at intervals of several-tens meters within spawning grounds for groupers (e.g., Colin, 2012;Nanami et al., 2017;Sadovy de Mitcheson et al., 2020).These studies have shown species-specific spatial variations among multiple grouper species and each species showed species-specific core sites, in which a very high density was found in a limited area within the spawning ground (Nanami et al., 2017;Sadovy de Mitcheson et al., 2020).Spatial variations of fish density at intervals of 250 m was also observed for two snapper species (Biggs & Nemeth, 2016), showing fine scale movement of two snapper species (L.cyanopterus and Lutjanus jocu) within the aggregation site in relation to spawning time.
The present study revealed a significant difference in the fine-scale spatial distributions in fish density between the two snapper species.Although the reasons why the two species showed species-specific spatial variations within the spawning ground remain unknown, this might be related with preference to a particular environmental condition or mating behavior.

Size and age frequency of fish aggregation
The present study showed that the maximum fork length of fish individuals at the spawning ground was 304 mm and 295.5 mm for L. decussatus and L. fulviflamma, respectively.In contrast, previous studies revealed that maximum fork length was 316.5 mm and 347.0 mm for Lutjanus decussatus and L. fulviflamma, respectively (Shimose & Nanami, 2015;Nanami, 2021).Since average fork length of fish individuals at the spawning ground was less than 250 mm for both species, most of the fish individuals forming the spawning aggregations were relatively small sized individuals.
A similar trend was also found for age frequency.Maximum age was respectively 24 and 23 for Lutjanus decussatus and L. fulviflamma (Shimose & Nanami, 2015;Nanami, 2021), whereas maximum age of fish individuals at the spawning ground was respectively 21 and 17 for L. decussatus and L. fulviflamma.Since average age of fish individuals at the spawning ground was less than 10 for both species, most of the fish individuals that form the spawning aggregations were relatively young individuals.

CONCLUSIONS
Results of the present study revealed that fish aggregations of the two snapper species (Lutjanus decussatus and L. fulviflamma) can be regarded as spawning aggregations due to their spatially and temporally predictable formation of fish aggregation as well as the presence of hydrated eggs in females within the aggregation site.In the Caribbean and Palau, spawning aggregations for some species of snappers have been previously found by local communities (Hamilton, Sadovy de Mitcheson & Aguilar-Perera, 2012).Since there is almost no local ecological information about the spawning aggregations of L. decussatus and L. fulviflamma, the present study is probably the first finding of a spawning aggregation of the two snapper species.As the two species are major fishery target species in the study region, the results of the present study should be applied to consider when and where a marine protected area should be established for effective protection of the spawning aggregation of the two species.Namely, the aggregation site of the two snapper species should be protected during their main spawning season (between April and October) and the protected area should completely cover the fish aggregations.
Table 3 Comparison about number of fish individuals per 600 m ×5 m of two snapper species between inside and outside aggregation site.Average number of fish individuals outside the aggregation site were obtained from underwater observations at 65 sites (Fig. S2).

Figure 1
Figure 1 Study site and fish aggregations of Lutjanus decussatus and L. fulviflamma.Location of Yaeyama Islands (A), Sekisei lagoon (B) (enclosed area by a yellow dotted line), fish aggregations of L. decussatus (C) and L. fulviflamma (D).Aerial photograph in (B) was provided by International Coral Reef Research and Monitoring Center.Full-size DOI: 10.7717/peerj.15991/fig-1

Figure 2 Figure 3
Figure 2 Temporal changes in number of fish individuals and gonad-somatic index for Lutjanus decussatus.Temporal changes in number of fish individuals with water temperature (A).Data in months enclosed by a dotted line in (A) are re-plotted in (B), showing lunar-related temporal changes in the number of fish individuals since a greater number of individuals were observed during the study period.Temporal changes in gonad-somatic index between April and October (C), together with oocyte developmental stage (pie charts).Numerals in pie charts indicate the number of females with eggs in each developmental stage.Lunar phases are abbreviated as LQM, last-quarter moon; NewM, New moon.''+'' and ''-'' mean after and before the lunar phase (e.g., ''NewM-1'' means 1 day before new moon).Oocyte developmental stage was abbreviated as PYS, primary yolk stage; TYS, tertiary yolk stage; MN, migration nuclear stage; PMA, pre-maturation stage; MA, maturation.Fish individual with hydrated egg (D) and crosssection of ovaries having oocytes with maturation stage (E).In (E), MA represents maturation stage of oocytes.Data of water temperature in (A) was provided by International Coral Reef Research and Monitoring Center.Full-size DOI: 10.7717/peerj.15991/fig-2

Figure 4
Figure 4 Correlogram showing auto-correlation coefficient with time lag, which examines the periodicity in fish increasing or decreasing at a particular time lag (for time lag, see Tables S1, S2 and Fig. S1).Horizontal dotted line represent 95% confidence interval of the auto-correlation coefficient.Black and white arrows represent significant positive and negative values of auto-correlation coefficient, showing significant periodicity at the time lag.Full-size DOI: 10.7717/peerj.15991/fig-4

Figure 5 NanamiFigure 6 NanamiFigure 7 Nanami
Figure 5 Spatial distribution of Lutjanus decussatus on the 600 m × 5 m line transect in the aggregation site.The transect was divided into 1-minute sub-transects.Fish data are shown as bubble plots and each bubble represents the fish density (number of individuals in 20 m × 5 m area) on each 1-minute subtransect.Cross marks represent no fishes in the sub-transect.Graphs at upper right side show the temporal changes in number of fish individuals (see Fig.2A) and red arrows represent the day when the data was collected.Lunar phases are abbreviated as LQM, last-quarter moon; NewM, New moon.''+'' and ''-'' mean after and before the lunar phase (e.g., ''NewM-1'' means 1 day before new moon).Full-size DOI: 10.7717/peerj.15991/fig-5

Figure 8 Nanami
Figure 8 Fish density frequency of Lutjanus decussatus for 1-minute sub-transect in the aggregation site.Fish density represents the number of fish individuals in the 20 m × 5 m area.Number of individuals at 1-minute sub-transect was converted to fish density (20 m × 5 m).Graphs at upper right side show the temporal changes in number of fish individuals (see Fig.3A) and red arrows represent the day when the data was collected.Lunar phases are abbreviated as LQM: last-quarter moon; NewM: New moon.''+'' and ''-'' mean after and before the lunar phase (e.g., ''NewM-1'' means 1 day before new moon).For actual spatial distributions, see Fig.7.Full-size DOI: 10.7717/peerj.15991/fig-8 Figure 9 (A-D) Size and age frequency of Lutjanus decussatus individuals that were captured at the aggregation site.Solid lines represent the probability density function.*: Among the 49 individuals in (B) age of one individual could not be identified due to difficulty in counting of number of opaque rings on otolith.Thus, sample size was 48 in (D).Full-size DOI: 10.7717/peerj.15991/fig-9 Figure 10 (A-D) Size and age frequency of Lutjanus fulviflamma individuals that were captured at the aggregation site.Solid lines represent the probability density function.Full-size DOI: 10.7717/peerj.15991/fig-10