Seasonal variations in the diversity and benthic community structure of subtidal artificial oyster reefs adjacent to the Luanhe River Estuary, Bohai Sea

Artificial oyster reefs provide important spawning and nursery grounds for a variety of fishes and large mobile crustaceans. Between July 2016 and May 2017, seasonal surveys of species composition and community structure were performed in the artificial oyster reef area and control area adjacent to the Luanhe River Estuary in China. During the survey year, 56 species belonging to 50 genera, 45 families, and 19 orders were recorded. The dominant economically important fish and mobile crustaceans were Hexagrammos otakii, Pholis fangi, Sebastes schlegelii, Charybdis japonica, and Oratosquilla oratoria. Resident fishes belonged to the Cynoglossidae, Paralichthyidae, Pleuronectidae, and Gobiidae families. Seasonally important fish species included Lateolabrax japonicus, Konosirus punctatus, Thryssa kammalensis, Hexagrammos agrammus, and Acanthopagrus schlegelii. The ranges of H' values among stations were 1.18–2.16, 0.65–1.75, 1.18–2.06, and 0.62–1.92 in spring, summer, autumn, and winter, respectively. The benthic organisms present in the community of artificial oyster reef areas can be classified into groups according to month and season. The abundance biomass curves showed that the oyster reef area in spring, autumn, and winter experienced low disturbance, whereas the community structure in summer was subject to large variations from external disturbance. We also found that as the age of the oyster reefs increased, the percentage of oysters in the low shell height group (< 40 mm) decreased. The oyster density was 324 ind/m2 for the reef created in 2016, 724 ind/m2 for the reef created in 2015, and 364 ind/m2 for the reef created in 2013. These findings can be used to develop suitable management strategies for the sustainable maintenance of artificial oyster reef ecosystems.


Study stations and sampling
The Luanhe River, which is 1200-km long, arises at the foot of the Yanshan Mountains and flows into Bohai Bay.The study area was adjacent to Xiangyun Cove, located near the Luanhe River Estuary, at the northeast part of Bohai Bay (Fig. 1).The study area (39° 10″14.78″-39°10″53.86″N, 118° 59″ 30.57″-119° 1″ 48.72″ E) is surrounded by a 4-km long, 8-m high breakwater.Since 2013, a series of stone and concrete artificial reef  July 27-30 2017.Some of the crab pots at the stations were lost during the study, so no data were available for these stations.The length of a single crab pot was 8 m, and five were connected together at each station.
Additionally, oysters were collected by SCUBA divers, who scraped them off 0.5 × 0.5 m 2 quadrats in July 2016.Both valves of C. gigas are large and concave with rippled shell layers 38 .Torigoe (1981) described them as 'oval to spatulate' 38 , but Wakiya (1929) described one adult shell as 'extremely elongated' 39 .We also used a ship's crane to raise two artificial oyster reefs in May 2017 and used shovels to scrape off all sessile attached organisms.All individuals in each oyster reef were enumerated.After the surveys, fishes and crustaceans and other organisms including oysters were identified to the lowest possible taxonomic level, counted, and weighed to the nearest 0.1 g of wet weight in the laboratory.
We used a YSI multi-parameter water quality analysis measurer (EXO-2, YSI, Yellow Springs, OH, USA) to measure the water temperature (°C), salinity (ppt), pH, dissolved oxygen content (mg L −1 ), total dissolved solids content (g L −1 ), chlorophyll content (μg L −1 ), and turbidity (NTU) among the survey stations (St.1-St.12,St.A1-St.A2) in 9:00-12:00, May 30 and August 11, 2017.In addition, we used an acoustic doppler velocimeter (Linquest, USA) to measure the velocity in the study area before and after the reef construction.Before the reef construction (July 30-31 and August 6-7 2007), the current in the study area was dominated by tidal current and was the characteristics of reciprocating flow with a weak tidal residual current.The measured maximum velocity of spring tide was 0.86 m s −1 with the flow direction of 252°, and that of neap tide was 0.66 m s −1 with the flow direction of 64°.After the reef construction (September 2017), the average velocity of the spring tide and neap tide among survey stations in the study area was 0.28-0.37m s −1 .

Statistical analysis
The Shannon-Weaver diversity index H' 40 , Pielou's evenness index J 41 , Margalef richness index D 42 , and Pinkas relative importance index IRI 43 were used to analyze the community diversity, and a cluster analysis was performed to analyze the community structure in different seasons.The calculation formulae were as follows: where p i is the proportion of total samples number (N) belonging to ith species; S is the total number of species; N i is the ratio of the species s number to the total sample number; W i is the ratio of the species i's biomass to total samples biomass; F i is the frequency of occurrence of species i at each station.In the cases of H' > 3, 1 ≤ H' ≤ 3, and 0 < H' < 1, the biological community was in a situation of undisturbed, moderately disturbed, and highly disturbed, respectively 44 .The cases of IRI > 1000, 100 ≤ IRI < 1000, 10 ≤ IRI < 100, and IRI < 10 indicated dominant, important, common, and rare species, respectively 45 .
Multivariate analysis was performed through hierarchical cluster analysis, which is used to delineate groups of related objects within nonmetric multidimensional scaling (nMDS) ordination space, and nMDS, which is an ordination procedure based on ranked similarity data.Regarding the cluster analysis, the square root of biomass at each station was calculated to stabilize the impact of dominant species in the system.In calculating the Bray-Curtis similarity coefficient matrix, two-dimensional scaling ranking and hierarchical clustering was obtained through the application of nMDS.The confidence of nMDS calculations was evaluated according to the stress coefficient.Stress coefficients show how well the multivariate pattern is represented within lower dimensional space 46,47 .
Finally, abundance biomass curves (ABCs) were used to understand the biological community disturbance degree according to the relative position of the abundance curve and biomass dominance curve in a same coordinate system.If the biomass curve was above the abundance curve, it indicated an "undisturbed" community; if the position of curves were in contrast to each other, it indicated a "serious interference"; if two curves were generally intersecting each other, it meant a "moderate interference".The ABCs were drawn based on the proportion of the biomass and abundance of each species in the total biomass and abundance of marine organisms in different seasons.The formula was expressed as follows 48 : where S is the species number, B i is the biomass cumulative percentage of the species I, A i is the abundance cumulative percentage of the species i.The cases of W > 0, W < 0, and W ≈ 0, indicate an undisturbed, disturbed, and moderately disturbed community, respectively.The closer the W value is to 1, the richness of each species is closer to the similarity, and vice versa for when the W value is closer to − 1.
The formulae of the niche width index 40 and niche overlap index 49 were as follows: www.nature.com/scientificreports/where P ij is the ratio of the ith species abundance at station j to the total abundance at station j; P kj is the ratio of the kth species abundance at station j to the total abundance at station j; B i is the niche width of species i; and R is the total station number.The larger the value of B i , the larger the species niche width.Q ik is the ecological overlap index of the species i and k with the value range of 0 and 1.The larger the Q ik , the stronger the similarity in resource utilization in the species i and k.The cases of B i ≥ 2.0, 1.0 ≤ B i < 2.0, and B i < 1.0 indicate a wide, medium, and narrow niche species, respectively.The cases of Q ik > 0.7, 0.4 < Q ik < 0.7, and Q ik < 0.4, within the value range of 0.0 to 1.0, indicate a high, median, and low niche overlap degree, respectively 50 .Ecological niche width refers to the sum of the varieties of resources that can be exploited by the organisms.It is an index of resource diversity exploited by the organism.When a species' ecological niche width is wider, the specialization of this species is smaller (i.e., the species has a wide ecological niche).A narrow niche width indicates that a species is specialized.Niche overlap means that two or more species with a similar niche inhabit the same space and share or compete for the same resource.These statistical analyses were conducted in the software PRIMER 5.0 (Plymouth Institute of Oceanography, Plymouth, UK) and the DPS platform 51 .

Species composition, dominant species, and biological community diversity
The  3).In addition, the IRI of A. amurensis in May (spring) was twice that in January (winter) (2155.25 vs. 1161.64).The IRI of C. japonica was similar between summer (9377.69 and 7737.86 for July and August) and autumn (8148.62),and these values were double that in May (4631.09)(Table 3).
Regarding community diversity, the mean H' values were 1. 73

Charybdis japonica
Japanese swimming crab
The analysis of seasonal variations revealed that H' and J were more similar among the months of March, May, July, September, January, and July in the control area (from 1.41 ± 0.22 to 1.97 ± 0.12 and from 0.63 ± 0.06 to 0.77 ± 0.10, respectively) compared with August (0.99 ± 0.27 and 0.48 ± 0.06) and December (0.84 ± 0.20 and 0.44 ± 0.10).The values of D were similar in March and July in the control area (2.61 ± 0.64 vs. 2.86 ± 0.32) (Table 4).www.nature.com/scientificreports/

Analysis of community structure and ABCs
The biological community of the artificial oyster reef area was classified into six groups according to months and seasons: group I (control area), group II (July, August, September), group III (May), group IV (December), group V (March), and group VI (January).These groupings were consistent with the results of nMDS and the clustering analysis (R = 0.824, P < 0.01) (Figs. 2, 3).In addition, the stress coefficient of 0.15 shown by the twodimensional scale sorting of each group indicated a certain explanatory significance for the clustering results of all the stations in different seasons (Figs. 2, 3).The dominance curves of the biomass and abundance all intersected and overlapped in different seasons, indicating that the artificial oyster reef ecosystem experienced "moderate interference" during the survey year.The biomass curves in the reef areas were above the abundance curves in May (spring), September (autumn), and January (winter) as well as in July in the control area, but not in the oyster reef areas in summer, indicating that the oyster reef areas in spring, autumn, and winter experienced low disturbance.In summer, the W value was negative (-0.001 in July) and lower than that in the other three seasons, and the curves intersected in several situations, indicating large variations due to external disturbance (Fig. 4).

Niche width and niche overlap
The niche width ranges of the marine organisms were 0.06-1.58    www.nature.com/scientificreports/65 (38.01%), and 26 (15.20%) in spring, summer, autumn, winter, and July in the control area, respectively.The overlap index value of A. ommaturus to other species was highest in summer and winter compared to in spring and autumn.The species A. ommaturus, H. otakii, and P. gravieri occurred at a high frequency and had a higher niche overlap degree compared with other species in winter.The index value of C. japonica was the highest, except for in winter, and that of S. schlegelii was the highest in summer compared with the other seasons (Table 5).

Analysis of oyster reefs and environmental factors
In the oyster reefs, 14 species of sessile organisms belonging to 14 genera, 13 families, and 10 orders were found.Among them, 10 and 13 species were found at St.1 (10 genera, 9 families, and 9 orders) and St.2 (13 general, 12 families, and 9 orders), respectively.They had higher biomass and abundance at St.2 compared to St.1 (Table 6).The density of oysters was highest for the oyster reef created in 2015, and then it decreased with increasing reef age (Fig. 5).The percentage of oysters in shell height group < 40 mm decreased with increasing reef age, whereas the percentage of oysters in the > 100 mm group increased with increasing reef age, indicating oyster growth.The value range of water temperature in May and August 2017 were 16.34-18.67°C and 25.63-27.83°C respectively in the reef area, and the value difference of mean values among stations in the reef area in May and August was 2 °C.In August, the average value of water temperature in the reef area was 1.07 °C lower than that www.nature.com/scientificreports/outside the reef area (26.90 °C against 28.07 °C).In terms of salinity, the average values in May and August were closed (32.52 ppt against 32.25 ppt); the average value in the reef area was a little higher than the outside (32.25 against 31.88 in August).The average value of pH in May and August was similar (8.87 against 8.81), and there was similar in and outside the reef area in August (8.81against 8.84).The dissolved oxygen value in May was higher than that of August in the reef area (7.91 against 6.03 mg L −1 ); the value of dissolved oxygen outside the reef area was higher than that of reef area in August (6.41against 6.03 mg L −1 ).The average and range value of total dissolved solids in May was 31.73 and 31.43-31.90g L −1 .The average value of chlorophyll in May and August was 6.54 and 5.03 μg L −1 respectively with the range of 2.41-9.68 and 3.48-7.38μg L −1 ; the chlorophyll value outside the reef was higher than that of reef area (5.52 against 5.03) in August.In terms of turbidity, the value in May was far lower than that of August (27.49against 76.08 NTU), and the value in the reef area was a little lower than outside the reef area (76.08 against 81.26 NTU) (Table 7).

Discussion
Oysters successfully create vertical complex relief structures that provide microhabitats 52 for a variety of marine organisms.Fifty-six species were found in the artificial oyster reef areas of this study.Oyster reefs are critical stable habitats for a variety of resident and seasonal resident fishes such as S. ommaturus, C. stigmatias, P. fangi,   56 suggested that the structured habitat constructed by oyster reefs can support a diversity of taxa, including macroalgae, sessile and mobile invertebrate epifauna, infauna, fish, and birds that may be present at decreased numbers or absent in adjacent unvegetated soft-sediment habitats .Dumbauld et al. 57 found the juvenile Dungeness crab (Cancer magister) numbers were enhanced in the oyster shell habitat and thereby compensated for habitat loss caused by dredging.Additionally, Tolley and Volety (2005) found that crab and fish density, biomass, and diversity were all greater on reefs compared with an unstructured sand bottom at a Florida location 58 .Humphries et al. 59 reported that the nekton assemblage at oyster reef sites had greater diversity, biomass, and abundance compared with mud bottom sites in Louisiana .It is well-recognized that artificial structures such as artificial oyster reefs can be quite different from those in adjacent rocky areas 60 and can comprise a diverse assemblage of macroalgae and filter-feeding invertebrates 61 .Powers et al. 62 verified that the emergent habitat provided by mesh bags led to increasing densities of mobile invertebrates and juvenile fish in the case of on-ground clam culture in the United States compared with the adjacent sand flats and natural seagrass areas .However, in the current study, the total species number and diversity H' of the control area were greater than those of the oyster reef area throughout the year-long study from spring to winter (27 vs. 12-26 and 0.84 ± 0.20 to 1.73 ± 0.30 vs. 1.97 ± 0.12).We speculate that because the control area is very close to the oyster reef area, the boundary between reefs and barren areas can attract more fishes and crustaceans than inside the area.Fish species such as gray snapper and crested goby juveniles feed outside the oyster reef area, exhibiting little diet overlap 54 .The restored oyster reef area might provide a corridor between sheltering and foraging grounds 63 ; therefore, the reefs in barren areas have a high heterogeneity, and there are complex interspecies relationships 64 .
The community structure of marine organisms in the artificial oyster reefs in the estuary was more easily influenced by disturbances throughout the year.All of the ABCs intersected and overlapped during the seasons, showing various degrees of disturbances (the shift between wet and dry seasons).In our study area, the natural annual runoff of Luanhe was reduced to zero during the summer (June to August), resulting in salinity changes 65 .In the Loxahatchee Estuary in the USA, the optimal salinity zone for oysters has shifted upstream, resulting in the death and subsequent burial of old oyster reefs 66 .Freshwater releases can flush oyster larvae downstream to locations that have unsuitable substrate and create unfavorable salinity conditions for larval survival in estuarine areas 67 .Oysters in southwest Florida spawn continuously, with peak recruitment occurring from May to November.At that time, large freshwater flows from Lake Okeechobee and conditions during the summer expose oyster larvae and associated organisms to low salinities and flush the larvae downstream to locations with substrates that are not suitable for settlement 67 .Variations in physico-chemical characteristics beneath the oyster habitat could lead to a displacement of large-bodied macrofauna (e.g., large bivalves, heart urchins, brittle stars) and the enhancement of small-bodied disturbance-tolerant opportunistic species (e.g., marine worms and capitellid polychaetes) 61 .
We also detected large seasonal and spatial variations in species composition and community structure.The biological community of the artificial oyster reef areas can be divided according to months.For example, the abundances of transient finfish averaged across all of the reefs over time revealed a strong seasonal pattern, but no obvious interannual pattern related to oyster abundance was reported 68 .Our study area is an important spawning and nursery ground for S. schlegelii from spring to summer.They migrate out of the reef areas in autumn when water temperature decreases.This fish's niche width in summer was lower than that in winter (1.44 vs. 1.82) because it mainly preys on the crustacean species A. chinensis in summer 69 .Grabowski reported that for reefs constructed in summer, the development of fish and mobile crustaceans was completed during the next spring-summer season 70 .We also found that the values of H', J, and D in summer were lower than those in the other three seasons.This result shows that the large biomass and abundance of several dominant species such as S. ommaturus, S. schlegelii, and C. japonica had larger impacts on the community structure in summer.
Artificial oyster reefs provide spatial refuges from predation and can alleviate food limitation 71,72 .For example, the range of niche overlap index values among seasons was 0.0-1.0,showing high heterogeneity among seasons and stations.Crustaceans such as L. vittata (117.99 g y −1 and 190 ind.y −1 ) in spring and P. gravieri (225.57g y −1 and 315 ind.y −1 ) and C. affinis (96.77 g y −1 and 101 ind.y −1 ) in winter may have a large impact on the community structure due to their high abundances.Oysters can promote pelagic fauna by preventing primary production from entering microbial loops, thus allowing it to pass up the food chain to bottom-feeding crustaceans 11 .The www.nature.com/scientificreports/resident crustacean species C. japonica, with large biomass and abundance (24,341.89g y −1 and 700 ind.y −1 ), had high niche overlap values with other organisms during the year from spring to autumn.Artificial oyster reefs on mudflats can increase the amount of habitat for these crustaceans, thus further increasing the secondary productivity of the estuary and prey organism foraging 73,74 .There are positive correlations between the diversity and abundance of reef-associated species and oyster shell height and biomass 75 .Reef macroinvertebrate communities respond positively to habitat restoration 76 .The dominance of the lower trophic level consumers might be attributed to the high productivity of the environment 77 .Furthermore, the economically important species A. fangsiao (217.42 g y −1 and 12 ind.y −1 ), O. variabilis (914.00 g y −1 and 6 ind.y −1 ), and R. venosa (9447.73g y −1 and 85 ind.y −1 ) had large biomass and abundance in this artificial oyster reef, highlighting the value of oyster reef habitats for maintaining high densities of resident species such as mollusks.
Hard substrates on the seabed such as live and dead oysters may provide novel habitats for fouling organisms and associated mobile biota 78 .The aggregation of various fish and crustacean species around the artificial structures, including artificial oyster reefs, is also well-recognized 79,80 .Xu et al. 81 investigated the community structure of polychaetes in reef and non-reef areas and found that the abundance of carnivore species in the reef area was higher than that in the non-reef area, which might be due to the halo effects of the oyster shell reef.Leguerrier et al. (2004) suggested that oyster aquaculture could increase the food supply to various fishes, which was predicted to be true as the result of increasing meiofaunal production 82 .However, there has been some discussion about the potential for direct negative effects of cultured oysters and mussels on fish population, primarily due to predation on fish eggs and larvae 83 .Dumbauld et al. 25 did not detect an overall increase in fish richness or abundance adjacent to oysters.
Regarding the development of policies for oyster reef protection, 60-80 fishing boats operate in the open sea of Tianjin Hangu Dashentang each year.Chinese fishermen locally harvested Volachlamys hirasei in the 1970s, Ranapa venosa in the 1980s, and oyster species such as Ostrea talienwhanensis, Ostrea denselamellosa, and C. gigas in the 1990s-2000s 16 .According to the fishing statistics of Dashentang village, about 100,000 t of living oysters were caught from 1999 to 2006, and 90% of them were C. gigas and O. talienwhanensis.This fishing behavior led to the rapid decline of the height of living oyster reefs and to the annual increase in the rate of empty shells.This was exacerbated by the sale of attachment substrates required for the survival of oyster larvae at a very low price 16 .Thus, we suggest that it is necessary to set up protection areas for artificial and wild oyster reefs and to implement actions such as quota catch and limits on catch size and period.

Summary
As an important foraging area and component of the estuarine landscape, artificial oyster reefs created in the barren areas near the Luanhe River Estuary provide important spawning and nursery grounds for a variety of fishes and large mobile crustaceans, and thus have a positive impact on biodiversity in the estuary area.In this study, the dominant species such as resident rocky fishes and large crustaceans appeared to have important impacts on the community structure and diversity of the artificial oyster reefs ecosystem.We found that as oyster reef age increased, the percentage of oysters in the low shell height group (< 40 mm) decreased and that of the high shell height group (> 100 mm) increased.The density of oysters was 324 ind/m 2 in the oyster reef created in 2016, 724 ind/m 2 in the oyster reef created in 2015, and 364 ind/m 2 in the oyster reef created in 2013.Future research should focus on the trade-off between the economic and ecological benefits of artificial oyster reef habitat restoration.A system of ecological assessment indicators for the sustainable management of artificial oyster reefs should be developed.The impact of global climate change with increasing water temperature brings uncertainties regarding the management of oyster reefs, especially in the summer.The species composition and community structure of the artificial oyster reefs within the temperate Luanhe River Estuary should also be compared with those in other temperate and subtropical zones globally.

Figure 2 .
Figure 2. Cluster analysis of the biological communities among sampling stations according to the Bray-Curtis similarity.

Figure 3 .
Figure 3. Non-metric multidimensional scaling ordination plot of the biological communities among sampling stations according to the Bray-Curtis similarity.The area is separated into six groups according to months.The stress coefficient value was 0.15 in this study.

Figure 4 .
Figure 4.The ABCs of cumulative dominance (unit: %) against species rank varying with months of the seasons.W is the index value of the ABC.The blue inverted triangles denote the biomass curve, and the green triangles denote the abundance curve.
www.nature.com/scientificreports/H. otakii, and S. schlegelii, economically important crustaceans such as C. japonica and O. oratoria in different seasons, and molluscs such as R. venosa and O. variabilis.Similarly, 15 fish, 10 crustacean, and 2 cephalopods species and important resident species including O. oratoria, C. joyneri, C. japonica, A. hexanema, K. punctatus, and L. beka were found in the subtidal wetland with sandy substrate in the Tianjin Dashentang oyster reefs 17 .Quan and Wang (2013) identified significantly greater densities and biomass of living natural intertidal oysters Crassostrea sikamea in oyster aquaculture gear than in intertidal oyster reefs in the Xiangshan Bay of Zhejiang Province, China.Greater Pielou evenness and Shannon index values were observed in the intertidal oyster reefs 53 .Yeager et al. (2011) reported that artificial oyster reefs provide a critical refuge for diverse communities of fishes and invertebrates such as blue crabs and striped bass 54 .Gregalis et al. (2009) found higher abundances of small demersal fish and sessile invertebrate species on restored oyster reefs compared with unstructured bottom sediment in coastal Alabama 55 .Ruesink et al.

Figure 5 .
Figure 5. Percentage of oysters collected in July 2016 in shell height groups (%) and oyster density (ind/m 2 ) in the reefs created in different years.The left y-axis indicates the height groups < than 40 mm (denoted by blue), 40-70 mm (denoted by reddish brown), 70-100 mm (denoted by gray), and > 100 mm (denoted by yellow).The right y-axis indicates the density variations denoted by hollow circle in the reefs created in 2013, 2014, 2015, and 2016.

were dominant species in the control area (Table
56species recorded during the survey year belonged to 50 genera, 45 families, and 19 orders, and Table1lists the numbers in these categories in each survey month, including control sampling.R. venosa was present in spring to autumn; other species such as Asterias amurensis and Octopus variabilis were present in spring, summer, and winter; and Chaeturichthys stigmatias and Alpheus japonicus were present in spring, autumn, and winter (Table2).Compared with the oyster reef area, Philyra platycheir and Johnius belangerii were found only in the control area, and the abundances of D. edwardsii (153 vs. 17) and C. stigmatias (190 vs. 27) were greater in the control area than in the oyster reef area (Table2).

2 kj Table 1 .
Total number and number by month of order, family, genus, and species in the oyster reef areas and control area in different seasons.

Table 2 .
Species composition classified by order and family among seasons and the biomass (unit: g) and abundance (unit: ind.) of each species in the artificial oyster reef areas and control area.The symbol • means present in this month and season.Values in parentheses represent the niche width index.-meansno records.B and N represent the biomass and number.St.11 and the lowest (0.51) at St.7 in March.The mean values in summer were 0.63 ± 0.06 in July and 0.48 ± 0.06 in August, with the highest (0.70) at St.4 in July and the lowest (0.39) at St.12 in August.The mean value in autumn was 0.71 ± 0.11, with the highest (0.89) at St.3 and the lowest (0.57) at St.7 in September.The mean values in winter were 0.44 ± 0.10 in December and 0.77 ± 0.10 in January, with the highest (0.90) at St.3 in January and the lowest (0.30) at St.10 in December.The mean value in the control area was 0.70 ± 0.03, with the highest (0.76) at St.A3 and the lowest (0.66) at St.B3.Regarding D, the mean values in spring were 2.61 ± 0.64 in March and 1.93 ± 0.63 in May, with the highest (3.49) at St.9 in March and the lowest (1.38) at St.11 in May.The mean values in summer were 1.68 ± 0.26 in July and 1.43 ± 0.49 in August, with the highest (2.17) at St.4 in July and the lowest (0.65) at St.5 in August.The mean value in autumn was 1.74 ± 0.44 in September, with the highest (2.46) at St.2 and the lowest (1.04) at St.3.The mean values in winter were 1.67 ± 0.23 in December and 1.72 ± 0.53 in January, with the highest (2.39) at St.7 and the lowest (0.83) at St.9 in January.The mean value in the control area was 2.86 ± 0.32, with the highest (3.34) at St.A1 and the lowest (2.39) at St.B2 (Table Vol:.(1234567890) Scientific Reports | (2023) 13:17650 | https://doi.org/10.1038/s41598-023-44176-6www.nature.com/scientificreports/(0.65) at St.5 in August.The mean H' value in autumn was 1.47 ± 0.27, with the highest (2.06) at St.2 and the lowest (1.18) at St.7.The mean H' values in winter were 0.84 ± 0.20 in December and 1.41 ± 0.34 in January, with the highest (1.92) at St.7 in January and the lowest (0.62) at St.10 in December.The mean H' value in the control area was 1.97 ± 0.12 in July, with the highest (2.12) at St.A1 and the lowest (1.80) at St.B2 (Table 2).The ranges of H' values among stations were 1.18-2.16,0.65-1.75,1.18-2.06,and 0.62-1.92 in spring, summer, autumn, and winter, respectively.For J, the mean values in spring were 0.74 ± 0.12 in March and 0.70 ± 0.07 in May, with the highest (0.86) at

Table 3 .
Pinkas relative importance index IRI value of different species varying with months of the seasons in the artificial oyster reef areas and control area.
in July in the control area, (Table1).The dominant species in each season had a wide niche width, indicating a positive relationship between niche width and IRI.The wide niche species were H. otakii (1.58)in March and C. japonica (1.57) in May; C. japonica (1.70) in July and S. ommaturus (1.59) in August; C. japonica (2.79) and S. ommaturus (2.68) in September; S. ommaturus (3.48) and P. gravieri (2.49) in January; and P. gravieri (2.55), C. japonica (2.46), and C. stigmatias (2.35) in July in the control area (Table1).The number of species with a wide niche width was greater in autumn and winter than in spring and summer.Regarding the niche overlap index, 171 pairs were detected.The number of pairs with a high, median, and low niche overlap degree were 35 (20.47%), 42 (24.56%),and 94 (54.97%) in spring, respectively; 44 (25.73%),

Table 6 .
Species composition classified by order and family and the biomass B (unit: g) and abundance N (unit: ind.) of each species in the oyster reefs in May 2017.