Contrasting Structural Changes in Two Benthic Invertebrate Communities After an Extraordinary Rainfall Event in Liuqiu Island, Taiwan

© 2025 Academia Sinica, Taiwan Open Access Contrasting Structural Changes in Two Benthic Invertebrate Communities After an Extraordinary Rainfall Event in Liuqiu Island, Taiwan Dun-Ru Kang1, Meng-Ying Kuo2, Shui-Kai Chang3,4 , Jing-Ying Wu1, Te-Yu Liao1, Mei-Fang Lin5,6 , Shyh-Min Chao7, and Li-Lian Liu1,4,* 1Department of Oceanography, National Sun Yat-Sen University, Kaohsiung 804, Taiwan. *Correspondence: E-mail: [email protected] (Liu) E-mail: [email protected] (Kang); [email protected] (Wu); [email protected] (Liao) 2Taiwan Ocean Research Institute, National Applied Research Laboratories, Kaohsiung 852, Taiwan. E-mail: [email protected] (Kuo) 3Graduate Institute of Marine Affairs, National Sun Yat-sen University, Kaohsiung 804, Taiwan. E-mail: [email protected] (Chang) 4Sustainable Ocean Governance Center, National Sun Yat-sen University, Kaohsiung 804, Taiwan 5Department of Marine Biotechnology and Resources, National Sun Yat-Sen University, Kaohsiung 804, Taiwan. E-mail: [email protected] (Lin) 6Doctoral Degree Program in Marine Biotechnology, National Sun Yat-sen University, Kaohsiung 804, Taiwan 7National Museum of Natural Science, Taichung 404, Taiwan. E-mail: [email protected] (Chao) Received 1 October 2024 / Accepted 16 April 2025 / Published 23 July 2025 Communicated by Benny K.K. Chan This study compared the changes in benthic community structures of the Shanfu (SF) and Duozaiping (DU) coral reef intertidal zones of Liuqiu Island between 2020 and 2022 after a prolonged precipitation event between late July and early August 2021. The trend of extreme rainfall on the island was also explored using the Standardized Precipitation Index (SPI) for coastal management use. The adjacent study sites of SF and DU have areas of 16,203 and 31,930 m2, respectively. Variations in the diversity indexes between sites (SF and DU) and years (2020, 2021, and 2022) were determined by two-way ANOVA tests. SF and DU differed significantly in species density, richness, dominance, diversity, and evenness indexes. Site factor had a more significant effect than year on these indices. The benthic community at DU was distinctly clustered into before, after, and 1-year after rainfall groups, but this was not the case for SF. Echinodermata was the taxon primarily responsible for the differences. Based on the meteorological data from 1997 to 2022, the SPIs varied from -1.83 (severe drought) to +1.90 (severe precipitation). Temporally, precipitation intensity and frequency increased during 2016–2021 compared with 2000–2006, with a peak precipitation event of 1,168 mm rainfall over 14 days in 2021. The broader reef flat at DU may have resulted in thermal and salinity stresses lasting longer for stenohaline echinoderms than at SF. Results from this island highlighted the differential vulnerability of the benthic communities from different sites to natural disturbances for the first time, allowing policymakers and stakeholders to formulate effective regional management strategies to minimize the impact of extreme climate events. Key words: Diversity index, Standardized Precipitation Index, SPI, Climate change, Low salinity, Echinoderms Citation: Kang DR, Kuo MY, Chang SK, Wu JY, Liao TY, Lin MF, Chao SM, Liu LL. 2025. Contrasting structural changes in two benthic invertebrate communities after an extraordinary rainfall event in Liuqiu Island, Taiwan. Zool Stud 64:23. doi:10.6620/ZS.2025.64-23. BACKGROUND Climate-induced disturbances to marine ecosystems are expected to increase over the next few decades, posing threats such as elevated temperature, ocean acidification, sea level rise, and increasing storm frequency and severity (IPCC 2022). These effects can have shortor long-term consequences and vary depending on the ecosystem. In 2021, the Emergency Event Database recorded 432 natural disasters Zoological Studies 64:23 (2025) doi:10.6620/ZS.2025.64-23 1 © 2025 Academia Sinica, Taiwan worldwide (CRED 2022), with 223 of these events (51.6%) being classified as flooding. Taiwan’s climate tends to be warm and humid, with a southwestern wind in summer and cool and dry with a northeastern prevailing wind in winter (Wu and Kuo 1999; Henny et al. 2021). Rainfall occurs mainly between May and September, depending on the frontal systems of the Mei-Yu season (Chen and Chen 2003) in May and June, with typhoons in July and August. In the Mei-Yu and typhoon seasons, rainfall is weighty over the southwest mountain slopes (Henny et al. 2021) especially between May and September. In this study, Taiwan’s extreme rainfall (ER. Studies analyzing extreme precipitation trends have shown that extraordinary rainfall events in Taiwan have increased markedly between 1960 and 2015, particularly since 2000, in all seasons except fall. Moreover, studies have shown that rainfall has become more intense and unpredictable from 2003 onwards (Tung et al. 2022). In southwestern Taiwan, there has been an overall increasing trend in both the intensity and frequency of extreme rainfall alongside total annual precipitation. Short bursts of excessively heavy rainfall generally result in flooding and landslides which are considered Taiwan’s most significant natural hazard. Extreme climate events, such as heatwaves, cold snaps, droughts, and intense rainfall, impact various marine ecosystems and adversely affect Taiwan’s fishing industry (Hsieh et al. 2008; MOA 2020). In Penghu, an extreme cold in the coral reef ecosystem was related to a two-week low water temperature (minimum 11.7℃) in 2008 (Hsieh et al. 2008). Records show that 183 fish species in 58 families suffered significant mortality. Deaths of macroinvertebrates, including echinoderms, crustaceans, mollusks, and corals, were mentioned, but no further qualitative or quantitative descriptions were available of these dead animals. Another study reported losses of more than 73 tons of wild fish and an 80% mortality rate among caged aquaculture fish (Chang et al. 2013). It took 53 months for the wild fish stocks to recover (Chen et al. 2020). A 96-hour rainfall of 2,058–2,888 mm was recorded in southern Taiwan when Typhoon Morakot hit Taiwan in early August 2009 (Central Weather Administration (CWA), Taiwan, ROC). This record reached the average annual precipitation level of 2,500 mm in the plains of Taiwan. Pond culture and coastal fishery losses were estimated at NT$4.7 billion, as reported by the Ministry of Agriculture (MOA 2020). However, no assessment was made of the impact on wild marine fauna. Although heavy rainfall or flooding is rare, their impacts on coral reefs can be severe and last more than a decade (Jokiel et al. 1993; van Woesik et al. 1995; Jones and Berkelmans 2014). The responses of benthic fauna to hyposaline conditions vary by location and taxa. For example, exceptional rainfall with heavy flooding occurred in the inner Gulf of Thailand between July and December 2011 (Sangmanee et al. 2012). At Khang Khao Island, the density of the sea urchin Diadema setosum decreased from 9.5 in July 2010 to 5.2 individuals/m2 in August 2011. Low salinity also killed more than 60% of corals on the island. Conversely, the mortality rate of the thorny oyster Spondylus cf. versicolor and the sea cucumber Holothuria leucospilota was less than 5%. In the Urbinu Lagoon, France, exceptional rainfall occurred between October 31 and November 1, 1993 (Fernandez et al. 2006). The salinity dropped to 23 psu over three days, and the sea urchin Paracentrotus lividus population was reduced by 50% between 1990 and 1994. Low salinity, turbidity, and siltation might have caused the population crash, and it would take six years to recover completely. Floodwaters associated with tropical cyclones affecting coral communities were reported at Keppel Island, Australia, in 1990 and 2011 (van Woesik et al. 1995; Jones and Berkelmans 2014) and in Hervey Bay in 2010 and 2013 (Bulter et al. 2015) and Hawaii (Banner 1968; Jokiel et al. 1993). Regarding Hawaii, the re-occurrence of freshwater “killing” events were documented in Kaneohe Bay in May 1965, January 1988, July 2014, and April 2018 (Banner 1968; Jokiel et al. 1993; Bahr et al. 2015a b; Rodgers et al. 2021). Mass mortality events from the intertidal zone to 2 m depth were reported, including corals, sea cucumbers, crabs, and fishes (Jokiel et al. 1993). As smaller organisms rotted away quickly, evidence of the “killing event” was most apparent in corals and larger invertebrates. The highest mortality rates were near the estuary mouth, then declined to zero within 1 km. In contrast to the subtidal zones, however, there were no mortality events or adverse effects on fish and other organisms. For example, littorinid snails, periwinkle snails, hermit crabs, limpets, and barnacles remained in good condition. In 2018, at Pila‘a and Ha‘ena on the island of Kaua‘i, a flooding event caused an increase in coral bleaching (55%). It decreased in sea urchin and fish populations (34–48% and 58%, respectively, Rodgers et al. 2021) and changed the community structure of corals. In 2021, many extreme climate events occurred worldwide, including extraordinarily high temperatures caused by heatwaves across the Americas and the Mediterranean, which were accompanied by devastating wildfires. Extraordinary levels of rainfall were experienced in China and parts of Europe, causing severe flooding, numerous casualties, and economic losses (WMO 2021). In early August 2021, within 123 page 2 of 13Zoological Studies 64:23 (2025) © 2025 Academia Sinica, Taiwan hours, Tropical Storm Lupit brought in 1,218–1,638 mm of rainfall in southern Taiwan (CWA, Taiwan, ROC). Currently, our knowledge of the influence of extreme climate events on coral reef ecosystems is mainly restricted to corals and fish (Jones and Berkelmans 2014; Bahr et al. 2015a b; Butler et al. 2015; Rodgers et al. 2021; Xie et al. 2020). Excessive precipitation is also expected to impact intertidal benthic invertebrate communities on coral reefs. However, studies of these impacts remain limited. This study used data from a regular monitoring program on benthic invertebrates at two coral reef sites on Liuqiu Island in southern Taiwan. The responses of the benthic invertebrate communities at Shanfu (SF) and Duozaiping (DU) were compared before and after the August 2021 rainfall event over a three-year investigation (2020–2022). MATERIALS AND METHODS Study area and sampling methodology Liuqiu Island, also known as Xiao Liuqiu, Siaoliouciou, and Sio Liu-khiu, is located 14 km off southwestern Taiwan and has approximately 13,000 residents. It is a tropical island fringed with coral reefs, with a land area of 6.8 km2 (Tourism Administration, Ministry of Transportation and Communications, Taiwan, ROC) and 12 km of coastline. The coast is a popular tourist attraction because of its scenic views and marine biodiversity. The monthly mean sea surface temperatures of the island range from 25.0 to 30.0℃ (2004–2023) (CWA, Taiwan, ROC). We investigated the benthic invertebrate communities between 2020 and 2022 during low tide at the adjacent sites of SF and DU (Fig. 1). The reef flats of the surveyed sites were 16,203 and 31,930 m2, respectively. We set up two transect lines at each site with a 1 × 1 m2 quadrat at 5-m intervals between the high and low tide levels. We counted the total number of invertebrates (> 5 mm) in the quadrats, which included crustaceans, mollusks, annelids, and echinoderms. Specimens were identified in the field or the lab using a magnifier or a microscope. Species identification was undertaken using reference studies for invertebrates in general (Lo 2013), crustaceans (Miyake 1982; Chen and Lo 2014), annelids (Fauchald 1977; Ge et al. 2018), sipunculas (Cutler 1994), mollusks (Okutani 2000; Chen et al. 2012; Chiu and Su 2019 2020; Jie Fig. 1. The study area and sampling sites in Liuqiu Island, Taiwan. N page 3 of 13Zoological Studies 64:23 (2025) © 2025 Academia Sinica, Taiwan 2019; Chiu et al. 2021), echinoderms (Chao 1998 2018; Chao and Su 2009), and sponges (Hooper et al. 2002). The designation of species name follows the WoRMS taxonomic system (http://www.marinespecies.org/). Data analysis We analyzed the data regarding K-dominance curves, density, and diversity indexes to obtain quantitative estimates of the spatial and temporal variability of the benthic communities at SF and DU. A K-dominance plot was applied to visualize the cumulative species abundance distribution pattern. Diversity indexes measure species richness (total number of species), evenness (relative abundance of each species), and composition (the identity of species) within a specific area. We used Margalef’s Richness index (R), Shannon’s Diversity index (Ho), Pielou’s Evenness index (J0), and Simpson’s Dominance index (D) (Clarke and Warwick 2001). We used two-way analyses of variance (ANOVAs) to separately evaluate the effects of site (SF and DU) and year (2020, 2021, and 2022) on the density and diversity indexes of the benthic communities. The Kolmogorov-Smirnov test was used to check data normality before the analyses. Whenever the ANOVA showed significant differences, the Tukey post hoc test was conducted to explore the differences in the interested factors. Temporal variations of the benthic communities at SF and DU were analyzed separately by cluster analyses of the Bray-Curtis similarity (BCS) indices (Primer 6.0) to investigate the influence of precipitation (Clarke and Warwick 2001). The contribution of each species to the cluster was further determined by principal component analysis (PCA). All statistical analyses were run on the Primer 6.0 software package (Primer-E Ltd., Plymouth, UK). Rainfall records of Liuqiu Island The heaviest rainfall on Liuqiu Island occurred between July 28 and August 10, 2021. The rainfall trend was investigated using the Standardized Precipitation Index (SPI) based on precipitation records from the CWA, ROC. The data included monthly total precipitation amount (mm), annual maximum consecutive precipitation days, and annual precipitation amount (mm). The SPI was calculated at a time scale of 12 months using precipitation data from 1997 to 2022 to explore the climate change trend on Liuqiu Island. The SPI was developed by McKee et al. (1993) to define and monitor drought events at multiple time scales (1-, 3-, 6-, 12-, and 24-months). The World Meteorological Organization (WMO) recommended the SPI as an index for describing meteorological droughts (Svoboda et al. 2012). It categorizes climate conditions ranging from extreme wetness ≥ +2 SPI to extreme drought ≤ -2 SPI; -1 to +1 as a normal regime; +1 to +1.5 and +1.5 to +2.0 as moderate and severe precipitation, respectively (McKee et al. 1993; Seiler et al. 2002). Conversely, moderate and severe drought are characterized by the same SPI ranges but with negative values. RESULTS Species composition and diversity indexes of the benthic communities at SF and DU The total number of benthic invertebrate species at SF and DU varied among sampling months in 2020–2022, i.e., 39–69 vs. 27–71, respectively (Table 1). The shape of the cumulative species dominance plots (K-dominance curves) had gentle slopes at SF contrasting at DU (Fig. 2). Echinoderms comprised 25.8–49.3 individuals/m2 at DU, and 3.8 to 19.7 individuals/m2 at SF (Fig. 3). Mean densities, Margalef’s Richness index (R), Shannon’s Diversity index (Ho), Pielou’s Evenness index (J0), and Simpson’s Dominance index (D) at SF and DU were 17–41 vs. 31–59 (individuals/m2), 7.8–16.9 vs. 17.6–27.8, 2.2–3.1 vs. 2.0–2.6, 0.66–0.83 vs. 0.54–0.70, and 0.08–0.25 vs. 0.14–0.23, respectively (Table 1). The two-way ANOVA tests showed that the site had a much more significant effect than year on density and diversity indexes, except for dominance (Table 2). Moreover, density and all the diversity indexes differed significantly between sites (2-way ANOVAs; d.f. = 1; F = 43.8, P < 0.001 for density; F = 73.1; P < 0.001 for richness index). DU had greater values than SF for density, richness, and dominance. By contrast, DU’s diversity and evenness values were lower than SF’s. There were no significant differences in density among 2020, 2021, and 2022. Diversity and Evenness indexes differed yearly in descending order of 2022 > 2021 > 2020 (2-way ANOVAs; d.f. = 2; F = 4.4, P < 0.001 for diversity index), while the richness and dominance indexes increased from 2020 to 2022 (2-way ANOVAs; d.f. = 2; F = 4.7, P < 0.05 for richness index). Spatial-temporal patterns of community structure at SF and DU The Bray-Curtis similarity analyses showed that the benthic community at DU was separated into three groups: before, after, and 1-year after rainfall (Fig. 4). At DU, the five principal components (PC1–PC5) page 4 of 13Zoological Studies 64:23 (2025) © 2025 Academia Sinica, Taiwan accounted for 53.3, 21.4, 11.1, 6.6, and 4.1% of the total variance, respectively. The separation was mainly contributed by PC1, including species of Breviturma brevipes (-0.594), Echinometra mathaei (-0.568), Ophiocoma scolopendrina (0.488), Ophiocoma sp. (-0.242), and Holothuria atra (-0.109), respectively (Fig. 5). In contrast, the clustered groups at SF did not correlate with the level of rainfall (Fig. 4). PC1−PC5 for SF comprised 58.8, 25.5, 5.7, 3.3, and 2.0% of the total variance, respectively. The major contributors of PC1 were E. mathaei (-0.931), Calcinus latens (0.171), Nerita albicilla (-0.142), Stomopneustes variolaris (-0.138), and Porites lobate (-0.103) (Fig. 5). The DU cluster demonstrated a year-long influence of rainfall on the benthic community. Fig. 2. K-dominance curves of the benthic communities at (a) Shanfu (SF) and (b) Duozaiping (DU) in 2020–2022. (n): the number of species. page 5 of 13Zoological Studies 64:23 (2025) © 2025 Academia Sinica, Taiwan Rainfall records of Liuqiu Island The precipitation records for Liuqiu Island between 1997 and 2022 showed that the lowest annual precipitation occurred in 2002 (833.5 mm) (Fig. 6), and the highest occurred in 2016 (2,682 mm). The yearly maximum consecutive precipitation days were 14, with four records that fell in August−September 2001 (225 mm), June 2016 (404.5 mm), June 2021 (644 mm), and July−August 2021 (1168.5 mm). In 2021, the extreme rainfall lasted from July 28 to August 10, with the most rain falling on August 1 (272.5 mm). The 12-month SPI values are shown in figure 6, along with the monthly total precipitation records for Liuqiu Island. Extreme precipitation and drought as ≥ +2 and ≤ -2 SPI were not observed. Severe rainfall occurred in 2016 and 2021 with SPI values of +1.81 and +1.90; conversely, severe drought occurred in 2002 (-1.83). The 2021 extreme rainfall matched the highest value of +1.90. Most of the consecutive positive SPI values of > +1 occurred after 2016. There seems to be a tendency for wetter than average periods since then. Fig. 3. The composition and density of invertebrate taxa at (a) Shanfu (SF) and (b) Duozaiping (DU) in 2020–2022. page 6 of 13Zoological Studies 64:23 (2025) © 2025 Academia Sinica, Taiwan DISCUSSION A 14-day extreme precipitation event (1,168 mm) was recorded in July/August 2021 (CWB, Taiwan, ROC), the highest level since 1997. Additionally, SPI values have shown that Liuqiu Island has experienced increasingly severe precipitation since 2016. The total number of species in the two benthic communities varied from 39 to 69 at SF and 27 to 71 at DU between 2020 and 2022. Site factor (SF and DU) had a more significant effect than year (2020, 2021, and 2022) on indexes of density, richness, dominance, diversity, Table 1. Diversity indices listed for the benthic communities of Shanfu (SF) and Duozaiping (DU). Density (no./m2); R: Margalef’s Richness index; Ho: Shannon’s Diversity index; J0: Pielou’s Evenness index; D: Simpson’s Dominance index Shanfu 2020 2021 2022 Month369389359 Species 43 61 60 62 46 39 56 59 69 Density 22.7 ± 3.3 40.7 ± 28.5 26.1 ± 6.4 22.4 ± 6.5 17.4 ± 5.2 20.8 ± 13.9 17.8 ± 6.7 18.3 ± 0.9 33.1 ± 7.8 R 11.7 ± 1.6 16.9 ± 10.0 11.4 ± 2.5 10.0 ± 2.0 9.0 ± 1.5 11.1 ± 4.6 7.8 ± 2.5 8.7 ± 0.1 13.1 ± 2.8 Ho 2.24 ± 0.44 2.61 ± 0.18 2.63 ± 0.29 2.76 ± 0.18 2.62 ± 0.41 2.19 ± 0.10 2.92 ± 0.24 2.93 ± 0.02 3.10 ± 0.01 J00.66 ± 0.12 0.71 ± 0.02 0.71 ± 0.02 0.75 ± 0.00 0.77 ± 0.03 0.68 ± 0.12 0.83 ± 0.03 0.81 ± 0.00 0.79 ± 0.01 D 0.25 ± 0.11 0.14 ± 0.03 0.15 ± 0.05 0.12 ± 0.02 0.14 ± 0.06 0.22 ± 0.09 0.09 ±0.03 0.09 ± 0.01 0.08 ± 0.01 Duozaiping 2020 2021 2022 Month3679389358 Species 42 52 52 49 52 64 71 43 51 27 Density 55.0 ± 8.5 48.9 ± 13.9 59.3 ± 14.0 57.8 ± 15.3 57.7 ± 9.9 57.8 ± 28.3 49.1 ± 10.4 49.2 ± 11.5 39.8 ± 1.1 30.7 ± 13.2 R 27.8 ± 4.2 24.2 ± 5.7 27.6 ± 5.3 25.9 ± 4.9 27.8 ± 5.7 25.0 ± 10.2 21.5 ± 7.1 22.0 ± 4.6 17.7 ± 0.5 17.6 ± 3.5 Ho 1.96 ± 0.08 2.03 ± 0.18 2.16 ± 0.01 2.27 ± 0.24 1.96 ± 0.21 2.48 ± 0.05 2.57 ± 0.50 2.34 ± 0.04 2.41 ± 0.08 1.99 ± 0.06 J00.57 ± 0.01 0.56 ± 0.07 0.59 ± 0.03 0.62 ± 0.03 0.54 ± 0.05 0.65 ± 0.02 0.65 ± 0.11 0.69 ± 0.03 0.69 ± 0.01 0.70 ± 0.09 D 0.22 ± 0.03 0.23 ± 0.04 0.18 ± 0.01 0.18 ± 0.03 0.23 ± 0.04 0.14 ± 0.00 0.15 ± 0.07 0.16 ± 0.01 0.15 ± 0.01 0.20 ± 0.01 Table 2. Summary of two-way ANOVA results for the variables of density and diversity indices Variable Factor DF MS F Sig. Density Site 1 6177.41 43.81 < 0.001 DUa, SFb Year 2 395.89 2.81 > 0.05 2020a, 2021a, 2022a Site * Year 2 238.63 1.69 > 0.05 R Site 1 1433.96 73.05 < 0.001 DUa, SFb Year 2 92.28 4.70 < 0.02 2020a, 2021ab, 2022b Site * Year 2 24.26 1.24 > 0.05 Ho Site 1 1.81 24.86 < 0.001 SFa, DUb Year 2 0.32 4.35 < 0.001 2022a, 2021ab, 2020b Site * Year 2 0.23 3.13 > 0.05 J0 Site 1 0.12 43.06 < 0.001 SFa, DUb Year 2 0.04 15.34 < 0.001 2022a, 2021b, 2020b Site * Year 2 0 0.04 > 0.05 D Site 1 0.01 5.54 < 0.03 DUa, SFb Year 2 0.01 5.18 < 0.01 2020a, 2021ab, 2022b Site * Year 2 0 1.71 > 0.05 Density (no./m2); R: Margalef’s Richness index; Ho: Shannon’s Diversity index; J0: Pielou’s Evenness index; D: Simpson’s Dominance index. Different letters indicate that the values significantly differ (Tukey’s test; p < 0.05). page 7 of 13Zoological Studies 64:23 (2025) © 2025 Academia Sinica, Taiwan and evenness. The benthic communities at DU were characterized by before, after, and 1-year after rainfall groups, chiefly comprising brittle stars (B. brevipes, O. scolopendrina, Ophiocoma sp.), sea urchins (E. mathaei), and sea cucumbers (Holothuria atra). Extreme precipitation in coastal areas had a year-long effect on benthic communities at DU but not at SF. The broader reef flat at DU, which contained more stenohaline echinoderms, suffered greater thermal and salinity stresses during low tide than SF, which may explain the difference. Abiotic factors, such as flow direction and Fig. 4. Dendrograms of the hierarchical clustering of the benthic communities based on Bray-Curtis similarities at (a) Shanfu (SF) and (b) Duozaiping (DU) in 2020–2022. page 8 of 13Zoological Studies 64:23 (2025) © 2025 Academia Sinica, Taiwan nutrient load, also govern the vulnerability of benthic communities to natural disturbances. Regarding global synthesis and meta-analysis of the contributions of coral reefs to coastal hazard risk reduction, Ferrario et al. (2014) reported that wave energy and wave height were attenuated significantly across more expansive reef flats. Conversely, as in the case of our study at DU and SF, large reef flats retain desiccation, thermal, and salinity stresses longer. Regarding salinity tolerance of benthic invertebrates, Echinodermata is a stenohaline group compared with other invertebrates. Based on the experimental and observational estimates in the North Sea (Geburzi et al. 2022), the salinity tolerances of benthic crustaceans (13 spp.), echinoderms (3 spp.), molluscs (3 spp.), and polychaetes (1 sp.) were reported as 0.1–55, 18–36, 4–60, and 0.5–45 psu, respectively. Another study investigated the effect of salinity fluctuations on echinoderm and mollusc communities at Caye Cochons (CC) and Cap du Corbeau (CdC) in the northern Gulf of St. Lawrence (Drouin et al. 1985). Sites of CC and CdC differed markedly in the amplitude of salinity fluctuations, i.e., 10–30 vs. 24–30 psu, respectively. In the intertidal area to 2 m deep, five echinoderms distributed at CdC, i.e., an urchin (Strongylocentrotus droebachiens), two ophiuroids (Ophiopholis aculeata and Ophiura robusta), two starfishes (Asterias vulgaris and Leptasterias polaris). However, at CC, O. aculeata and O. robusta were less abundant, and there was no A. vulgaris. The situation was different concerning mollusc species. The abundance of five molluscs was similar at both sites, including a chiton (Tonicella marmorea), a limpet (Acmaea testudinalis), a snail (Buccinum undatum), and two clams (Hiatella arctica and Mytilus edulis). Meanwhile, there was less chiton of Tonicella rubra at CC. There was a striking difference in the echinoderm diversity between CC and CdC but not molluscs. This discrepancy was likely due to the lower salinity tolerance of echinoderms than molluscs. Consistent with these studies, our results showed a higher proportion of stenohaline echinoderms in the benthic invertebrate communities at DU (72–88%) than at SF (22–68%). Heavy rainfall or flooding is rare, and their impacts on benthic coral reef communities primarily focus on corals and fish. For example, the flooding associated with Tropical Cyclone Joy affected nearshore coral communities in the Keppel Islands of Australia from December 23, 1990, to January 7, 1991 (van Woesik et al. 1995). Coral mortality reached 85% on leeward reefs, particularly affecting Acropora spp., Pocillopora damicornis, and Seriatopora hystrix. Subsequently, the community composition shifted to more flood-resistant species. Moreover, Tropical Cyclone Tasha caused flooding in Keppel Bay in 2011 (Jones and Berkelmans 2014). This time, coral mortality reached 100% to 2 m depth. Based on the previous flood in 1991, coral recovery is expected to take 10–15 years. In another study in Hervey Bay, Australia, flooding in 2010 and 2013 caused coral mortality of 28% and 40%, respectively (Bulter et al. 2015), while changes in water quality parameters and salinity in submarine groundwater lasted for six and ten months, respectively. A study in Hong Kong showed that a minimum fiveFig. 5. The contribution of each species in the benthic communities of (a) Shanfu (SF) and (b) Duozaiping (DU) to the first and second principal component scores in 2020–2022. page 9 of 13Zoological Studies 64:23 (2025)