Fish Sci (2014) 80:21–41 DOI 10.1007/s12562-013-0682-x
ORIGINAL ARTICLE
Biology
Differences in fish assemblage structure between vegetated and unvegetated microhabitats in relation to food abundance patterns in a mangrove creek Kusuto Nanjo • Hiroyoshi Kohno • Yohei Nakamura Masahiro Horinouchi • Mitsuhiko Sano
•
Received: 31 August 2013 / Accepted: 20 October 2013 / Published online: 11 December 2013 Ó The Japanese Society of Fisheries Science 2013
Abstract In order to clarify the mechanisms determining fish distribution patterns in a mangrove system on Iriomote Island, in southern Japan, fish assemblage structures were determined by visual observation, along with food abundance and environmental factors, in an area of mangrove roots on the banks, and a bare sand area at the center, within downstream, midstream and upstream portions of a branch creek from the Urauchi River. The fish assemblage structures differed significantly between the area types, with the mangrove-root area supporting a more diverse and abundant fish fauna. A canonical correspondence analysis revealed that the relationships between fish distribution and their food abundance differed among trophic groups. Benthic crustacean or plant feeders were positively associated with their prey i.e. crabs and macroalgae—in other words,
these trophic groups were abundant in downstream and/ or midstream mangrove-root areas in which their prey were also particularly abundant. However, zooplankton feeders did not show such relationships, their abundance being positively associated with fine sediment particles (characteristic of areas with weak water movement). These results suggested that food availability is a major factor determining the distribution patterns of benthic crustacean feeders and plant feeders, whereas for zooplankton feeders other factors, such as sheltering effects against water current and/or predators, may be more significant.
K. Nanjo (&) Atmosphere and Ocean Research Institute, The University of Tokyo, Kashiwa, Chiba 277-8564, Japan e-mail:
[email protected]
Introduction
H. Kohno Okinawa Regional Research Center, Tokai University, Uehara, Taketomi, Okinawa 907-1541, Japan Y. Nakamura Graduate School of Kuroshio Science, Kochi University, 200 Monobe, Nankoku, Kochi 783-8502, Japan M. Horinouchi Research Center for Coastal Lagoon Environments, Shimane University, 1060, Nishikawatsu-cho, Matsue, Shimane 690-8504, Japan M. Sano Graduate School of Agricultural and Life Sciences, The University of Tokyo, Bunkyo, Tokyo 113-8657, Japan
Keywords Distribution Fish assemblage Food abundance Mangrove estuary Microhabitat Trophic group
Tropical and subtropical mangrove estuaries often support large numbers of fish species and individuals, including some that are commercially important, compared to nearby open bare sand/mud flats, which are usually characterized by fewer species and fewer individuals [1–5]. Habitat complexity structured by mangrove vegetation has been considered a major factor responsible for such a difference due to greater habitat complexity providing shelter from predation and/or strong water movement, increased microhabitat availability, and abundant food [2, 6]. For example, the structural complexity of mangrove prop roots and pneumatophores could reduce predation risks [7, 8], such being especially important for small-sized fishes, such as juveniles, because of their vulnerability to predation [9]; thus, the latter sometimes
123
22
select mangrove areas in order to reduce predation risks [10]. In addition, the complex structure of mangrove roots enhances sedimentation rates, trapping suspended particles and leaf litter [11], and also harbor abundant epiphytic algae on their large surface area [12, 13], thereby enabling large numbers of small invertebrates (food resources for mangrove fishes), which feed upon such organic materials [14, 15], to reside within the mangrove system [7]. A greater abundance of food provides potential advantages, such as faster growth [16] to a size that is relatively free from predation [9, 17, 18]. Accordingly, small-sized fishes sometimes prefer mangrove areas to unvegetated, bare areas [7]. Several recent studies of fish assemblages in mangrove habitats, however, have suggested that fish species richness and/or abundance are not always greater in mangrove areas [19–21]. For example, in temperate Australia, Smith and Hindell [19] examined fish assemblage structures in the Barwon River Estuary, using fyke, gill and pop nets, and found that fishes were less abundant in the mangroveforested areas than in the unvegetated channels. In addition, Tse et al. [22] reported that fish species richness and gut fullness did not differ between a mangrove estuary and unvegetated mudflat within Tolo Harbour in southern China, suggesting that mangrove habitats may not necessarily be a better feeding area for fishes compared to mudflats. These findings are contrary to the general notion that mangrove-vegetated areas support a large number of fishes by harboring rich prey resources. Accordingly, the present study re-evaluated the relationships between fish distribution and prey abundance in a mangrove system. We conducted the study on a small spatial scale (i.e., vegetated/unvegetated microhabitats within a single mangrove creek), since the abiotic environmental factors (e.g., temperature) associated with such may not differ strongly between microhabitats, whereas they may affect fish distribution on a broader spatial scale due to significant differences among microhabitats, which may lead to the misinterpretation of results [23]. In order to re-evaluate the relationships between fish distribution and prey abundance in a mangrove system, we examined and compared fish assemblage structures and their prey abundance in mangrove-vegetated areas and adjacent open areas within a single mangrove creek.
Materials and methods Study site The study was conducted in a creek of the Urauchi River (24°240 N, 123°460 E), situated on the northern side of Iriomote Island, Ryukyu Islands, Japan (Fig. 1), in August
123
Fish Sci (2014) 80:21–41
(summer) and November (autumn) of 2008, February (winter) and May (spring) of 2009, August (summer) and November (autumn) of 2009, and February (winter) and May (spring) of 2010. The former 4 months were regarded arbitrarily as the first year, and the latter 4 months as the second year. The study creek (1.1-km-long, 10–80 m in width until 900 m from the creek mouth), branching from the main river at a point about 2 km upstream from the river mouth, was fringed by dense, undisturbed mature mangrove forest, comprising mainly the red mangrove Rhizophora stylosa. Sampling stations were established at the lower- (0–250 m from the creek mouth), middle- (350–550 m from the creek mouth) and upper-portions (650–850 m from the creek mouth) of the study creek (hereafter referred to as ‘‘downstream’’, ‘‘midstream’’ and ‘‘upstream’’, respectively. See Fig. 1). Two microhabitat types within each station were targeted: (1) an offshore edge portion of the sand/mud area with mangrove roots on the creek bank (hereafter referred to as ‘‘mangrove-root area’’), and (2) an unvegetated sand area in the center of the creek (hereafter referred to as ‘‘bare sand area’’). The mean creek widths were 45.0 ± 12.3 m (standard deviation) downstream, 29.0 ± 6.1 m midstream and 23.4 ± 2.9 m upstream (n = 5). The mean root densities of R. stylosa in the mangrove-root area were 66.6 ± 14.8/ m2 downstream, 63.0 ± 12.6/m2 midstream and 59.0 ± 11.6/m2 upstream (n = 5). The tidal range within the creek was approximately 1.5 m, prop roots being inundated at high tide and partially exposed at low tide. Water temperature and salinity were measured at 10-min intervals by a small data logger (Compact-CT, JFE Advantech Co, Ltd.) anchored to a prop root, whereas water turbidity was measured using a portable water quality analyzer (Model U-20XD, Horiba Ltd.) at high tide, with ten replicates at each station in each season. The preliminary study conducted from 2007 through 2008 confirmed no difference in the three environmental factors between the microhabitats, although they varied among stations and/or seasons (three-way analysis of variance, water temperature: f = 0.014, p [ 0.05 for microhabitats; f = 1.563, p [ 0.05 for stations; f = 266.996, p \ 0.05 for seasons; salinity: f = 0.215, p [ 0.05 for microhabitats; f = 4.006, p \ 0.05 for stations, f = 2.814, p \ 0.05 for seasons; water turbidity: f = 0.038, p [ 0.05 for microhabitats; f = 27.015, p \ 0.05 for stations, f = 2.792, p \ 0.05 for seasons) (Nanjo K, unpubl. data). Sediment samples were taken from five points located randomly in each microhabitat at each station in each season. Dry sediment samples were sieved through mesh trays of decreasing size (2,000, 1,000, 500, 250, 180, 125, 63 lm) stacked vertically on a vibratory sieve shaker (Analysette 3 Pro, Fritsch Co.). The sediment retained in each sieve was weighed to the nearest 0.001 g and mean median grain sizes were calculated.
Fish Sci (2014) 80:21–41
23
each transect was positioned within 5–10 m from the outer border of the mangrove-root area. All fishes within the transect areas were counted by an observer (Nanjo K) using a snorkel, and their total lengths were estimated in size classes of 1 cm. Individual fishes were identified to species levels following Nakabo [24].
East China Sea Japan 123˚ 46’ E
Food abundance
China Taiwan
Iriomote Island Ryukyu Islands 24˚ 24’ N
0
1 km
ORRC
Urauchi Bay
Urauchi River
N E
W Downstream
Midstream
S
Upstream
Fig. 1 Map of the Urauchi River mangrove estuary, Iriomote Island, Ryukyu Islands, Japan. Shaded areas indicate mangrove forests. Circles show each station (downstream, midstream and upstream). filled square, Okinawa Regional Research Center, Tokai University (ORRC)
Sampling design Fishes In order to clarify fish assemblage structures, visual censusing was conducted in both area types at high tide (ca. water depth 1.2 m) between 09:00 and 17:00 h at each station in each season in both years. For each census, sixbelt transects 1 m wide and 20 m long, parallel to the creek flow direction and separated from each other by at least 10 m, were established randomly in each microhabitat at depths of ca. 1.2 m (downstream), 1.4 m (midstream) and 0.8 m (upstream), using a scaled measure. For the mangrove-root area census, each transect was placed within the outer border of that area. For the bare sand area census,
Zooplankton, crabs, gammaridean amphipods, polychaetes, macroalgae and detritus were chosen as major food resources, since these food items were the most important (especially for the dominant species) for the mangrove fish assemblage at the present study site [25]. The six food types were collected at each station in each microhabitat in each season in both years. Zooplankton was collected at high tide, while the other food items were sampled at low tide between 09:00 and 17:00 h. Zooplankton, such as calanoid and cyclopoid copepods and their nauplii, was collected using a plankton net (20cm mouth diameter, 75-cm-long and 100-lm mesh size) equipped with a flow meter. The net was held in front of the collector (Nanjo K), who walked against the water flow for 20 m. Subsequently, the net contents were retrieved and the volume of water filtered was recorded. In the mangrove-root area, the collector maintained a position as close as possible to the prop roots. Five replicates were obtained for each microhabitat. Zooplankton densities were expressed as individual numbers per m3. In order to collect macroalgae, such as the red algae Bostrychia spp., attached to prop roots, mangrove root-bark samples were taken from roots some 30 cm above the substrate (ca. 20 cm2/root), with ten replicates collected from each mangrove-root area. The lack of above-ground structures in the bare sand areas precluded the collection of macroalgal samples. Epiphytic macroalgae on root-bark were removed (scraped off) using a small knife and weighed to the nearest 0.001 g (wet weight), with the algal biomass expressed as wet weight per surface area of each root-bark sample (g/cm2). Crabs on the sediment surface were collected from five quadrats (1 9 1 m) established randomly within each microhabitat. In the mangrove-root areas, crabs on the roots were also captured. Fiddler crabs (Uca spp.) were excluded from the crab samples because they are not taken by fishes in this mangrove estuary [see 25]. Gammaridean amphipods and polychaetes were sampled using a cylindrical core sampler (18 cm in diameter), collecting 500-cm3 of the sediment from five points located randomly in each microhabitat. Immediately after collection, the samples were preserved in 10 % buffered formalin, and the zooplankton, gammaridean amphipods, and polychaetes within each sample were later sorted and counted under a binocular microscope.
123
24
The amount of detritus in the sediment was expressed by ignition loss. Five (ca. 10 g) sediment samples were collected to a depth of 1 cm from the surface using a small spoon, dried to a constant weight and then combusted at 560 °C for 8 h, and reweighed. The amount of detritus was expressed as the percentage of the dry weight reduced by combustion to the total dry sediment weight. Data analysis Species richness and abundance of fishes in each microhabitat were expressed as mean species and individual numbers per transect (20 m2), respectively. In addition, mangrove fishes were divided into seven trophic categories based on published dietary data [e.g., 25], including three well-represented groups—benthic crustacean feeders (preying mostly on crabs, shrimps and gammaridean amphipods), zooplankton feeders, and plant feeders (taking mainly macroalgae)—in addition to polychaete feeders, detritus feeders, fish feeders and insect feeders. The latter trophic categories were represented mostly by a few lowdensity species (see Appendix) and were not treated specifically in the analyses described below. Since the objective of the present study was to determine differences among microhabitats, stations and seasons, rather than differences between years, statistical analyses were made in each year. Three-way analyses of variance (ANOVAs) were employed to compare fish assemblage properties (i.e., mean species and individual numbers per transect of the overall fish assemblage and of each dominant trophic group) and abundances of food resources, i.e., zooplankton (ind/m3), crabs (ind/m2), gammaridean amphipods (ind/500 cm3), polychaetes (ind/ 500 cm3) and detritus (%), among microhabitats, stations and seasons. Since neither plant nor zooplankton feeders occurred in the bare sand area, mean species and individual numbers of such groups in the mangrove-root areas were compared among stations and seasons, using two-way ANOVA. For algal biomass (g/cm2), two-way ANOVA was employed for the same reason. When significant differences were recognized (p \ 0.05), post hoc Games– Howell tests were applied. Because the first-order interactions (microhabitat 9 station, microhabitat 9 season or station 9 season) were often significant (see Table 2), oneway ANOVAs and Games–Howell tests were conducted to compare the values between microhabitats at each station in each season, among stations in each microhabitat in each season, and among seasons in each microhabitat at each station. Some ANOVA results are not presented here due to space limitations. Prior to ANOVA, the detritus data (i.e., ignition loss) were transformed to arcsine, whereas other data were square-root (x ? 1) transformed, to improve homogeneity of variances. Such transformation, however,
123
Fish Sci (2014) 80:21–41
did not produce homogeneous variances in some cases. In such cases, to compensate for the increased likelihood of Type 1 errors (the rejection of a true null hypothesis), the significance level was set at 0.01 [26]. In order to classify the assemblages into groups with similar species occurrence patterns, a cluster analysis was conducted. Initially, the degree of similarity of the assemblages between microhabitats at each station in each season in each year was calculated, using the Bray–Curtis similarity index, based on the individual numbers of each species. Subsequently, the resultant similarity matrix was subjected to an average linkage clustering method. Prior to the cluster analysis, data were square-root (x ? 1) transformed. Species contributing significant variation between the groups were identified using the SIMPER (similarity percentage) subroutine. In order to assess the relationships among fish distributions and biotic (food abundance) and abiotic (environmental) factors, a canonical correspondence analysis (CCA) [27] was conducted. For the CCA, data for ten dominant fish species from the study site were used to express species variables (see Table 1), in order to reduce the effect of rare species distorting the ordination of CCA [28]. The abundances of six food resources (zooplankton, crabs, gammaridean amphipods, polychaetes, detritus, and macroalgae) and environmental factors (water temperature, salinity, water turbidity, and median grain size of sediment) were used as explanatory variables.
Results Fish assemblage structure A total of 31,055 individuals, representing 29 families and 85 species, were observed during the study period (Table 1 and Appendix). In terms of individual numbers, the brackish cardinal fish Apogon amboinensis (8,539 individuals), hover goby Parioglossus raoi (8,499), P. palustris (3,375), freshwater demoiselle Neopomacentrus taeniurus (1,194), brackish damselfish Pomacentrus taeniometopon (1,137), mangrove snapper Lutjanus argentimaculatus (971), blacktail snapper L. fulvus (923), Ambasis miops (880), Parioglossus rainfordi (848) and Pomadasys argenteus (746) were dominant, accounting for 87.3 % of the total. In the present study, most fishes, including dominant species, were observed in the mangrove-root area at all stations in all seasons, although only a few low-density species, such as mojarra Gerres erythrourus, thornfish Terapon jarbua and sand goby Favonigobius reichei, occurred in the bare sand area (see Fig. 2; Table 1 and Appendix). One-way ANOVAs conducted following threeway ANOVAs (see explanation in Data analysis; Table 2)
B B
Apogon amboinensis (A. amb)
Lutjanus argentimaculatus (L. arg)
Lutjanus fulvus (L. ful)
Lutjanidae
Z
Parioglossus raoi (P. rao)
B B
Apogon amboinensis (A. amb)
Lutjanus argentimaculatus (L. arg)
Lutjanus fulvus (L. ful)
Lutjanidae
2–4
Z Z
Parioglossus rainfordi (P. rai)
Parioglossus raoi (P. rao)
6.3
19.2
21.2
4.5
20.5
4.3
23.8
5.3
4.8
25.2
16.0
6.3
14.3
4.2
33.3
23.0
4.3
1.8
0.7
13.2
13.7
3.0
9.8
9.0
186.0
2.8
19.2
7.3
5.5
9.5
10.3
4.3
6.0
7.8
135.3
S
30.7
1.0
4.7
1.2
1.7
2.5
24.5
44.8
6.5
1.0
0.5
1.0
14.2
8.3
M
S
2.7
1.3
13.5
24.8
6.0
21.3
4.8
51.3
13.2
2.7
17.5
6.5
15.0
13.8
3.0
45.7
0.2
M
S
M
M
S
Downstream
Midstream
Downstream
Upstream
Autumn
Summer
16.2
0.5
17.8
9.7
12.3
2.7
5.5
27.3
136.2
0.2
14.3
4.3
4.0
7.3
12.0
11.5
2.8
28.5
125.2
15.5
M
S
Midstream
14.8
5.5
46.7
0.3
2.0
13.5
16.3
6.3
44.2
0.2
0.8
5.2
5.0
M
S
Upstream
29.3
4.2
10.0
10.7
0.2
16.5
8.5
17.5
68.7
13.8
9.2
11.7
9.0
2.2
8.0
6.0
24.5
M
S
Downstream
Winter
6.7
81.8
10.0
44.8
5.8
10.3
8.7
3.5
14.8
107.2
6.7
81.8
3.5
17.7
6.7
3.0
9.8
3.5
11.7
120.5
M
S
Midstream
TC trophic category, TL total length (cm), M mangrove-root area, S bare sand area. Trophic category: B benthic crustacean feeders, Z zooplankton feeders, Pl plant feeders
1–5
1–4
Z
1–14
1–8
10–31
3–20
3–50
3–7
1–6
1–5
2–4
1–4
1–14
1–8
10–31
2–24
3–55
3–7
1–5
TL
Parioglossus palustris (P. pal)
Pl
Pomacentrus taeniometopon (P. tae)
Ptereleotridae
Pl
Neopomacentrus taeniurus (N. tae)
Pomacentridae
B
Pomadasys argenteus (P. arg)
Haemulidae
B
Ambassis miops (A. mio)
Apogonidae
Z
Z
Parioglossus rainfordi (P. rai)
Ambassidae
Second year
Z
Parioglossus palustris (P. pal)
Pl
Pomacentrus taeniometopon (P. tae)
Ptereleotridae
Pl
Neopomacentrus taeniurus (N. tae)
Pomacentridae
B
Pomadasys argenteus (P. arg)
Haemulidae
B
Ambassis miops (A. mio)
Z
TC
Apogonidae
Species (Abbreviation)
Ambassidae
First year
Family
Table 1 Mean individual number of dominant fish species per transect (20 m2, n = 6) at the study site in the first and second years
60.7
14.2
53.2
0.3
0.2
1.3
13.0
78.5
14.3
64.0
0.2
1.5
0.8
M
S
Upstream
51.0
7.0
5.5
13.8
33.3
15.0
12.5
9.3
21.7
45.8
12.0
10.8
10.0
10.3
12.2
11.3
4.8
40.0
16.7
M
S
Downstream
Spring
82.0
3.3
29.8
9.0
6.3
4.0
4.2
3.0
115.3
31.7
99.2
19.2
41.3
7.8
1.0
19.0
2.5
4.2
82.5
10.0
M
S
Midstream
0.3
0.8
101.2
1.7
62.7
0.5
0.2
0.3
1.7
6.2
121.8
1.7
64.3
M
S
Upstream
Fish Sci (2014) 80:21–41 25
123
26
Fish Sci (2014) 80:21–41
revealed that the mean total species and individual numbers were significantly higher in the mangrove-root area than in the bare sand area at all stations in all seasons in both years (p \ 0.001). In the mangrove-root area, mean total midstream species and individual numbers were higher compared to downstream and/or upstream in each season in both years, although in some cases the differences were marginal (for species number in winter in the second year, p = 0.034; for individual number in winter in the first year, p = 0.021; for other cases, all p \ 0.01). In contrast, the mean total species and individual numbers in the bare sand areas seldom differed among stations and seasons. The trophic guild structure of the fish assemblage is summarized in Fig. 3. Because benthic crustacean feeders, plant feeders, and zooplankton feeders were dominant in the present study (see Fig. 3, Appendix), these three categories were applied to the analyses. One-way ANOVAs conducted First year
25
Mean species number/20 m 2
following three-way ANOVAs (see explanation in Data analysis; Table 2) indicated that the mean species and individual numbers of benthic crustacean feeders were higher in the mangrove-root area than in the bare sand area at each station in each season in both years, although in some cases the differences were marginal (for mean upstream individual number in winter and spring in the second year, p = 0.011; for other cases, all p \ 0.01). In the mangrove-root area, the mean individual number was greater downstream and/or midstream compared to upstream in each season in both years (for winter and spring in the second year, p \ 0.05; for other cases, all p \ 0.01). Since the plant feeders and zooplankton feeders never occurred in the bare sand area (see Fig. 3; Table 1 and Appendix), the mean species and individual numbers of these groups occurring in the mangrove-root area were compared among stations and seasons by two-way Second year
20
15
10
5
0
Mean individual number/20 m2
400
300
200
100 5
0 DS MS US
DS MS US
DS MS US
DS MS US
DS MS US
DS MS US
DS MS US
DS MSUS
Summer
Autumn
Winter
Spring
Summer
Autumn
Winter
Spring
Fig. 2 Mean fish species and individual numbers per transect (20 m2, n = 6) in the mangrove-root (filled square) and bare sand (open square) areas at each station in each season in the first and second years. Bars indicate standard deviation. DS downstream, MS midstream, US upstream
123
3.000
120
Residuals
3.000
6
Microhabitat 9 Station 9 Season
Individuals
Species
3 6 6 120
Station 9 Season
Microhabitat 9 Station 9 Season
Residuals
2
Microhabitat 9 Station
Microhabitat 9 Season
3
6 120
Microhabitat 9 Station 9 Season
Residuals 2
6
Station 9 Season
Season
3
Microhabitat 9 Season
Station
2
Microhabitat 9 Station
1
3
Season
Microhabitat
1 2
Microhabitat Station
18.00
3.000
1.600
1.000
7.200
203.5
6.100
212.6
1,774
0.110
0.030
0.120
0.480
0.830
0.100
105.8 1.710
2.000
3 6
Microhabitat 9 Season
Station 9 Season
88.00
18.00
3 2
Season
Microhabitat 9 Station
87.00
2
5,615
0.120
Station
120
Residuals
0.080
0.090
0.820
0.380 2.370
1
6
Microhabitat 9 Station 9 Season
240.5 4.120
Microhabitat
3 6
Microhabitat 9 Season
3 2
Season Microhabitat 9 Station
Station 9 Season
1 2
Microhabitat
Station
Benthic crustacean feeders
Individuals
Species
Overall fish assemblage 2,004
0.539
0.318
2.389
67.94
2.047
70.98
592.1
0.303
1.085
4.497
7.752
0.895
983.4 15.92
0.979
0.781
6.367
31.59
6.489
31.50
2,024
0.692
0.709
6.873
3.148 19.74
34.33
3
3
0.778
0.926
120
6
6
3
2
\0.001** 0.072
3
2
\0.001** 0.111
1
120
6
6
\0.001**
0.934
0.375
3
2
\0.001** 0.005*
3
1 2
120
6
0.446
\0.001** \0.001**
0.443
6
\0.001** 0.587
3 2
\0.001**
2
\0.001** \0.001**
1
120
6
\0.001**
0.656
6
\0.001** 0.643
3 2
0.028 \0.001**
1 2
\0.001**
Games–Howell test
\0.001**
p
3.900
1.900
1.700
12.70
171.9
13.90
177.2
2,119
0.090
0.010
0.030
0.380
0.490
0.590
120.0 0.730
3.000
3.000
3.000
17.00
121.0
16.00
124.0
5,516
0.090
0.050
0.080
0.370
0.630 1.590
1.980
246.1
MS
df
f
df
MS
Second year
First year
0.480
0.446
3.264
44.17
3.578
45.53
544.6
0.114
0.387
4.343
5.572
6.701
1,364 8.341
1.005
1.122
6.073
43.47
5.596
44.45
1,982
0.536
0.927
4.251
7.219 18.20
22.64
2,816
f
0.822
0.846
0.024
\0.001**
0.016
\0.001**
\0.001**
0.995
0.886
0.006*
0.005*
\0.001**
\0.001** \0.001**
0.426
0.356
\0.001**
\0.001**
0.001*
\0.001**
\0.001**
0.780
0.478
\0.001**
\0.001** \0.001**
\0.001**
\0.001**
p
Games–Howell test
Table 2 Results of three-way ANOVA and two-way ANOVA testing differences in the mean numbers of fish species and individuals per transect (20 m2, n = 6) among microhabitats, stations and seasons in the first year and second year
Fish Sci (2014) 80:21–41 27
123
123 p
Games–Howell test
3 6 60
Season
Station 9 Season
Residuals
60
Residuals 2
6
Station 9 Season
Station
2 3
60
Residuals
Station
6
Station 9 Season
Season
3
60
Residuals
Season
6
Station 9 Season 2
3
Season
Station
2
Station
2.760
0.710
16.98
59.63
0.077
0.039
1.175
3.724
7.960
5.380
207.1
88.10
0.064
0.314
0.120
0.777
0.258
6.157
21.62
0.514
15.33
48.58
0.677
26.03
11.07
4.902
1.874
12.15
Sum = Aut \ Win = Spr
\0.001**
0.954
0.001*
\0.001**
US \ DS = MS
Win = Spr \ Sum
\0.001** 0.795
US \ DS = MS
\0.001**
0.669
DS \ US
* p \ 0.01, ** p \ 0.001
60
6
3
2
60
6
3
2
60
6
3
2
60
6
\0.001**
3
0.144
2
\0.001**
\0.001**
2.360
1.260
6.350
113.5
0.070
0.076
0.449
3.628
7.890
14.20
253.5
173.6
0.067
0.110
0.532
0.882
MS
df
f
df
MS
Second year
First year
DS downstream, MS midstream, US upstream, Sum summer, Aut Autumn, Win winter, Spr spring
Individuals
Species
Plant feeders
Individuals
Species
Zooplankton feeders
Table 2 continued
0.534
2.689
48.11
1.085
6.385
51.54
1.800
32.14
22.01
1.638
7.907
13.11
f
\0.001**
0.780
0.054
\0.001**
0.382
\0.001**
\0.001**
0.114
\0.001**
US \ DS = MS
US \ DS = MS
Sum \ Win \ Spr, Aut \ Spr
DS \ MS = US
Sum \ Win \ Spr, Aut \ Spr
\0.001** 0.152
DS \ MS = US
Games–Howell test
\0.001**
p
28 Fish Sci (2014) 80:21–41
Fish Sci (2014) 80:21–41 Benthic crustacean feeders
Zooplankton feeders
First year
25
Plant feeders
Others
Second year
20
Mean species number/20 m2
Fig. 3 Mean fish species and individual numbers per transect (20 m2, n = 6) for each trophic group in the mangrove-root and bare sand areas at each station in each season in the first and second years. DS downstream, MS midstream, US upstream
29
15
10
5
0 350
Mean individual number/20 m2
300 250
200
150 100 50
0 DS MS US DS MS US Summer
Autumn
ANOVA. The mean species and individual numbers of plant feeders were greater downstream and midstream than upstream in both years (Table 2), and the mean species number also varied significantly among seasons in both years, being higher in summer than in winter and spring of both, although a post hoc test did not detect a significant difference (Table 2). The mean species number of zooplankton feeders tended to be higher midstream and/or upstream than downstream, and highest in spring and lowest in summer in both years, although the differences between pairs were not always statistically significant in the first year (Table 2). The mean individual number showed similar tendencies, being higher in midstream and/or upstream than downstream, and highest in winter and spring (lowest in summer and autumn) (Table 2). The cluster analyses revealed that fish assemblages in each microhabitat fell into six spatial groups (Groups 1-1–1-6) in the first year and five (Groups 2-1–2-5) in the
DS MS US DS MS US Winter
Spring
DS MS US DS MS US Summer
Autumn
DS MS US DS MS US Winter
Spring
second, at a similarity level of 45 (Fig. 4). In both years, fish assemblages in one microhabitat type sometimes grouped together, but never with assemblages in the alternative type. The mangrove-root area assemblages could be subdivided into downstream and midstream (Groups 1-1 and 2-1), and upstream (Groups 1-2 and 2-2) assemblages, whereas such among-station differences were not observed for the bare sand area assemblages. In Groups 1-1 and 2-1, benthic crustacean feeders, such as Apogon amboinensis, Lutjanus argentimaculatus and L. fulvus, plant feeders (Neopomacentrus taeniurus and Pomacentrus taeniometopon) and zooplankton feeders (Parioglossus raoi and P. rainfordi) were dominant species. SIMPER analysis showed that those species contributed greatly to similarities within the groups (Table 3). In Group 1-2 and 2-2, zooplanktivous fish P. raoi and P. palustris were abundant and primarily responsible for the similarity within the groups (Table 3). Within the bare sand area
123
30
Fish Sci (2014) 80:21–41
Spring Spring Winter May Autumn Winter Autumn Summer Winter Summer Autumn Spring Winter Autumn Summer Spring Winter Autumn Summer Summer Winter Autumn Spring Summer
MS DS MS US MS DS DS DS US US US US US US US MS MS MS DS MS DS DS DS MS
Microhabitat S S 1–6 S S S 1–5 S S 1–4 S S S 1–3 S M M 1–2 M M M M M M 1–1 M M M M S
Winter Autumn Winter Autumn Spring Spring Winter Summer Spring Summer Summer Autumn Autumn Summer Spring Winter Spring Winter Autumn Summer Spring Winter Summer Autumn
MS DS DS MS US MS US US DS MS DS US US US US US DS DS DS DS MS MS MS MS
S S S S S S S S S S S S M M M M M M M M M M M M
Season
(a) First year
Station
*
(b) Second year
0
20
40
60
80
2–5
2–4
2–3
2–2
Food abundance and environmental factors The spatial and seasonal abundance patterns of the studied food resources are summarized in Fig. 5. Food item abundance sometimes differed among microhabitats, stations and/or seasons. Statistical values for some comparisons are shown in Table 4 (some values omitted due to space limitations). Environmental factors, shown in Table 5, are not dealt with specifically in the text for the same reason, although some brief descriptions are provided below. Of the main food items, crabs, polychaetes and detritus were more abundant in the areas with mangrove roots than in the bare sand areas, whereas zooplankton tended to be more abundant in the bare sand areas, although the differences were only marginally significant (Fig. 5. See also Table 4). Gammaridean amphipod density did not differ between the microhabitats. Macroalgae were attached to the prop roots in the mangrove-root area at each station, but were lacking in the adjacent open microhabitats, which consisted of soft substrata lacking above-ground, hard structures. The abundance of macroalgae, crabs, gammaridean amphipods and detritus differed among stations, although not so for the remaining items (Table 4). Macroalgae were significantly more abundant downstream and midstream than upstream (Table 4), whereas crabs and detritus were mostly abundant midstream and upstream, respectively (Fig. 5), although post hoc tests did not always detect significant differences between station-pairs (Table 4). Zooplankton density showed seasonal differences in both years (Table 4), being highest in summer and lowest in winter (Fig. 5. See also Table 4). Polychaete density showed a similar seasonal fluctuation pattern in the first year, although not in the second (Table 4). No seasonal differences in abundance were apparent for the remaining items (Table 4).
2–1
100
Bray-Curtis similarity index Fig. 4 Dendrograms of the cluster analyses showing similarities of fish assemblages based on the density of each fish species in the mangrove-root and bare sand areas at each station in each season in the first and second years. DS downstream, MS midstream, US upstream, M mangrove-root area, S bare sand area. At a Bray-Curtis similarity index level of 45, the assemblages were divided into six groups in the first year (1-1–1-6) and five groups in the second year (2-1–2-5). Asterisk indicates an assemblage not belonging to any group
123
groups, Gerres erythrourus, Terapon jarbua, Bathygobius fuscus and Favonigobius reichei were the dominant species.
The relationships between fishes, food resources and abiotic factors The first two axes of the CCA ordination explained 87.1 % of the variances of species-explanatory variable bi-plots in the first year (axis 1, 67.0 %; axis 2, 20.1 %) and 83.9 % of the variance in the second year (axis 1, 62.0 %; axis 2, 21.9 %) (Fig. 6). In the first year, the vectors of detritus and water turbidity, which tended to be more abundant or higher upstream, were on the right (positive) side of axis 1, whereas the vectors of other factors, such as macroalgae, crabs and salinity, which tended to be more abundant or higher downstream and/or midstream, were on the left
Fish Sci (2014) 80:21–41
31
Table 3 Results of a SIMPER analysis showing contributing rates of fish species to similarity of each group in the first year and second year First year Group/species
Second year TC
Mean ind.
Contribution (%)
Apogon amboinensis
B
75.9
17.2
Apogon amboinensis
B
82.4
17.7
Parioglossus raoi
Z
43.4
12.0
Neopomacentrus taeniurus
Pl
16.6
9.7
Group 1-1
Group/species
TC
Mean ind.
Contribution (%)
Group 2-1
Pomacentrus taeniometopon
Pl
12.0
7.5
Parioglossus raoi
Z
34.2
9.4
Pomadasys argenteus
B
10.0
6.2
Pomacentrus taeniometopon
Pl
11.8
8.9
Parioglossus rainfordi
Z
8.5
5.9
Lutjanus fulvus
B
11.8
7.9
Neopomacentrus taeniurus Lutjanus fulvus
Pl B
8.5 7.8
5.7 5.4
Lutjanus argentimaculatus Pomadasys argenteus
B B
10.1 5.5
6.9 4.4
Parioglossus palustris
Z
12.9
3.3
Lutjanus argentimaculatus
B
8.8
5.3
Parioglossus palustris
Z
11.1
3.5
Parioglossus raoi
Z
51.9
20.4
Parioglossus raoi
Z
65.4
29.1
Parioglossus palustris
Z
41.8
18.2
Parioglossus palustris
Z
44.8
22.9
Apogon amboinensis
B
14.3
12.4
Group 1-2
Group 2-2
Acentrogobius moloanus
B
2.5
6.6
Acentrogobius moloanus
B
3.8
5.6
Lutjanus argentimaculatus
B
1.0
4.5
Parioglossus rainfordi
Z
5.6
5.3
Terapon jarbua
B
1.5
4.4
Lutjanus argentimaculatus
B
1.9
5.0
Parioglossus rainfordi
Z
5.6
4.0
Terapon jarbua
B
1.9
4.1
B
0.2
100
B
0.5
100
Favonigobius reichei
B
0.7
43.1
B
0.2
100
Yongeichthys criniger Gerres erythrourus
B B
0.3 0.3
39.8 17.1
B
0.2
100
B
0.2
100
Favonigobius reichei
B
0.5
57.6
Bathygobius fuscus
B
0.2
42.4
Group 1-3 Terapon jarbua
Group 2-3
Group 1-4
Terapon jarbua Group 2-4 Favonigobius reichei Group 2-5 Bathygobius fuscus
Group 1-5 Favonigobius reichei Group 1-6
TC trophic category, Mean ind. mean individual number per transect (20 m2, n = 6) in each group. Trophic category: B benthic crustacean feeders, Z zooplankton feeders, Pl plant feeders
(negative) side (Fig. 6a). Macroalgae (correlation coefficient, r = -0.81), detritus (r = 0.77), salinity (r = -0.64) and crabs (r = -0.60) were highly correlated with axis 1, and water turbidity (r = 0.50) and crabs (r = 0.48) were highly correlated with axis 2. Regarding the relationships of fish trophic groups and environmental and food resource vectors in the ordination, the dominant species of benthic crustacean feeders, such as A. amboinensis, P. argenteus, L. argentimaculatus and L. fulvus, and plant feeders (N. taeniurus and P. taeniometopon), located on the left side in the bi-plots, were positively associated with their main food items (i.e., crabs and macroalgae) and salinity, which tended to be higher midstream and downstream (Fig. 6a). In contrast, most zooplankton feeders (P. raoi, P. palustris and P. rainfordi) fell on the right side, being positively associated with detritus and water turbidity, but negatively with prey abundance.
Ambassis miops, occupying the center of axis 1 and upper part of the ordination, was positively associated with water turbidity (Fig. 6a). Similar results were obtained for the second year (Fig. 6b).
Discussion The present study revealed that fish assemblage structures differed very considerably between microhabitats, such differences being consistent across stations and seasons. The mangrove-root area consistently supported greater numbers of fish species and individuals compared to the unvegetated substrata, contrary to the findings of several previous studies [19–21], even though supported the general notion that mangrove-vegetated areas harbored large numbers of fish species and individuals [e.g., 1].
123
32
Fish Sci (2014) 80:21–41 First year Mean biomass (g/cm2 )
0.15
Second year
Macroalgae
0.10
0.05
Mean individual number/m 2
Mean individual number/m 3
0 30000
Zooplankton
20000
10000
0
Crabs 20
10
Mean individual number/500 cm 3
Mean individual number/500 cm3
0
Gammaridean amphipods 400
200
0
400
200
0 20
Ignition loss (%)
Polychaetes
Detritus
10
0
DS MS US DS MS US Summer
Autumn
DS MS US
DS MS US
DS MS US
DS MS US
DS MS US
DS MS US
Winter
Spring
Summer
Autumn
Winter
Spring
Fig. 5 Mean abundance of each food resource in the mangrove-root (filled square) and bare sand (open square) areas at each station in each season in the first and second years. Abundances of each food resource expressed as: mean biomass (g/cm2, n = 5) for macroalgae, mean individual numbers for zooplankton (/m3, n = 5), crabs (/m2,
123
n = 5), gammaridean amphipods (/500 cm3, n = 5), polychaetes (/500 cm3, n = 5) and ignition loss (%) of sediment (n = 5) for detritus. Bars indicate standard deviation. DS downstream, MS midstream, US upstream
3 6 108
Season
Station 9 Season
Residuals
3 2 3 6 6 96
Station
Season
Microhabitat 9 Station
Microhabitat 9 Season
Station 9 Season
Microhabitat 9 Station 9 Season
Residuals
3 2 3 6 6 96
Season
Microhabitat 9 Station
Microhabitat 9 Season
Station 9 Season
Microhabitat 9 Station 9 Season
Residuals
3 2 3 6 6 96
Station
Season
Microhabitat 9 Station
Microhabitat 9 Season
Station 9 Season
Microhabitat 9 Station 9 Season
Residuals 1 2 3
Microhabitat
Station
Season
Polychaetes
1 2
Microhabitat
Gammaridean amphipods
1 2
Microhabitat
Station
Crabs
1 2
Microhabitat
Zooplankton
2
Station
Macroalgae
38.90
6.800
537.1
13.73
9.060
20.63
12.45
73.44
18.06
102.3
4.920
0.340
0.190
0.080
0.150
6.500
0.350
5.410
138.2
461.0
98.00
167.0
169.0
262.0
31,186
1,561
2,188
0.0002
0.0001
0.0007
0.0067
9.561
1.682
132.1
0.660
1.503
0.906
5.349
1.315
7.448
0.358
0.551
0.242
0.448
19.25
1.024
16.02
409.6
0.213
0.363
0.366
0.568
67.62
3.386
4.745
0.596
3.710
37.42
2 3
\0.001**
1
96
6
6
3
2
3
0.191
\0.001**
0.682
0.186
0.441
0.006*
0.274
1 2
0.551
96
6
6
3
\0.001**
0.768
0.961
0.720
3 2
0.385 \0.001**
1 2
\0.001**
96
6
6
3
\0.001**
0.972
0.901
0.778
2
3
0.569
\0.001**
1
108
6
3
2
2 Win \ Aut \ Spr \ Sum
US \ DS = MS
0.038
0.032
0.733
0.014
\0.001**
Games–Howell test
12.80
22.10
653.9
20.43
5.700
6.100
44.45
85.62
38.45
188.4
18.09
0.022
0.004
0.005
0.028
0.408
0.011
0.347
9.590
885.0
88.00
156.0
481.0
75.00
47,203
2,138
4,117
0.0004
0.0002
0.0001
0.0100
MS
df
p
Second year f
df
MS
First year
1.974
3.397
100.6
0.279
0.299
2.176
4.192
1.883
9.225
0.885
0.168
0.211
1.266
18.27
0.488
15.56
429.7
0.099
0.176
0.543
0.085
53.31
2.414
4.649
0.670
0.292
28.19
f
0.123
0.038
\0.001**
0.945
0.936
0.096
0.018
0.138
\0.001**
0.349
0.985
0.972
0.291
\0.001**
0.691
\0.001**
\0.001**
0.996
0.983
0.654
0.918
\0.001**
0.095
0.034
0.674
0.831
\0.001**
p
MS \ DS = US
Win \ Aut \ Spr = Sum
US \ DS = MS
Games–Howell test
Table 4 Results of three-way ANOVA and two-way ANOVA testing differences in the mean abundances of each food resource among microhabitats, stations and seasons in the first year and second year
Fish Sci (2014) 80:21–41 33
123
Fish Sci (2014) 80:21–41
Table 4 continued
Greater food availability may be a factor responsible for greater fish species richness and abundance in the mangrove-root area. Mangrove root structures trap suspended organic particles [11] and harbor abundant epiphytic algae [12, 13], such organic materials attracting many small invertebrates, which in turn serve as important food resources for fish [e.g., 7]. In fact, the present study demonstrated that crabs and polychaetes, along with macroalgae and detritus, were most abundant in the mangrove-root areas, fishes utilizing such food items also frequenting the mangrove-root areas but seldom or never the sand areas. This strongly supported the importance of food availability in their habitat choice. The abundance patterns of two dominant trophic groups (benthic crustacean feeders and plant feeders) from downstream to upstream matched those of their food items. Benthic crustacean feeders were most abundant in the midstream mangrove-root area where their prey crabs (mainly grapsids, sesarmids and ocypodids) were also most abundant. Plant feeders also frequented the mangrove-root areas, grazing actively upon macroalgae on the roots (Nanjo K, pers. obs.), but never occurred in the sand areas lacking such macroalgae. They were especially abundant downstream and midstream where macroalgae flourished densely. Additionally, in the CCA ordination, most component species of the dominant trophic groups (e.g., L. argentimaculatus and P. taeniometopon) showed a positive correlation with food abundance, further evidence for the importance of food availability in determining fish distributions in a mangrove system. Because a greater food supply provides a potential advantage for fishes, such as faster growth [16], many mangrove fishes prefer mangrove-root areas with abundant food [see also 13]. The latter conducted field experiments at the Caribbean Island Curac¸ao using artificial mangrove-root units and found that herbivorous fish species occurred in the experimental units with both mangrove roots and fouling algae but not in units without algae, thereby concluding that food was important for fish distribution. The importance of food availability for fish distribution has been pointed out elsewhere [29–32]. Contrary to the findings for the above two trophic groups, the distribution patterns of zooplankton feeders could not be explained by food availability. Although they were restricted entirely to the mangrove-root area, zooplankton tended to be more abundant in the bare sand area. Additionally, in the CCA ordination, most component species did not show a positive correlation with their prey. A similar phenomenon, the abundance patterns of planktivorous fishes not coinciding with those of pelagic copepods in Alligator Creek, northern Australia, previously reported by Robertson et al. [33], adds weight to the possibility that factors other than food availability may be
0.988 0.152 0.080
0.010 6
96
0.276
0.833
0.261
\0.001** 11.21 2 \0.001**
0.920
1.309 0.110 3 0.356
0.465
1.357 0.110
0.040 6
3
\0.001**
\0.001** 14.02
1,261 103.6 1
2
\0.001**
\0.001**
1.150
6.500 96
0.967
0.738
0.688
0.895 0.372 2.400 6
0.227 1.500 6
0.422
0.376 2.400
2.700 3
2
df df
123
* p \ 0.01, ** p \ 0.001
DS downstream, MS midstream, US upstream, Sum summer, Aut Autumn, Win winter, Spr spring
0.450 0.970 Residuals
0.070 6
96
Microhabitat 9 Station 9 Season
0.080
25.33 2 Station 9 Season
1.910
1.092 3 Microhabitat 9 Season
0.080
0.358
0.159 1.588
1.089
6 Microhabitat 9 Station
0.080 3 Season
0.120
15.33
1,179 88.72 1
2
Microhabitat
Station
Detritus
4.100 96 Residuals
1.150
0.803 0.505 6 Microhabitat 9 Station 9 Season
2.100
0.302 1.221 6 Station 9 Season
5.000
0.006*
0.841 0.174
4.389 17.80 3 Microhabitat 9 Season
0.700 2 Microhabitat 9 Station
MS
f
p
Games–Howell test
Second year First year
MS
f
p
Games–Howell test
34
4.9 ± 4.6
2.0 ± 1.0
1.7 ± 1.1
Water turbidity (NTU)
3.8 ± 1.5
M mangrove-root area, S bare sand area
157.6 ± 10
Water turbidity (NTU)
Grain size (lm)
22.3 ± 2.2
24.8 ± 9.4
Water temperature (°C)
Salinity
Second year
150.6 ± 14
22.9 ± 1.9
25.3 ± 8.2
200.4 ± 23
218.8 ± 24
22.1 ± 2.2
162.6 ± 5
8.1 ± 2.8
21.9 ± 8.5
21.1 ± 2.5
158.2 ± 4
3.3 ± 1.2
26.2 ± 2.3
227.2 ± 21
236.8 ± 18
134.8 ± 20
12.6 ± 7.9
15.6 ± 10.2
20.5 ± 2.5
130.4 ± 20
5.7 ± 2.5
25.9 ± 2.7
21.4 ± 2.8
M
Upstream
183.6 ± 18
181.0 ± 10
S
178.8 ± 8
25.1 ± 1.5
150.2 ± 7
3.1 ± 1.7
22.1 ± 4.8
157.0 ± 8
6.7 ± 4.4
27.2 ± 4.6
26.3 ± 2.2
151.8 ± 13
5.1 ± 1.5
27.7 ± 7.3
25.6 ± 2.0
M
M
S
7.7 ± 3.5 122.6 ± 15
3.7 ± 1.8 142.8 ± 24
M
229.0 ± 30
30.3 ± 2.4 24.8 ± 2.8
184.0 ± 15
Downstream
Water temperature (°C)
Grain size (lm)
5.3 ± 1.0 154.4 ± 17
9.2 ± 7.1 109.6 ± 13
Midstream
Salinity
First year
30.0 ± 2.1 27.4 ± 3.2
215.2 ± 23
23.8 ± 9.4
23.6 ± 2.3
M
Downstream S
5.3 ± 2.7 146.0 ± 7.2
29.9 ± 1.6 29.4 ± 3.1
S
Spring
196.8 ± 8
189.8 ± 15
30.1 ± 1.7 30.9 ± 2.2
M
Winter
142.6 ± 18
Water turbidity (NTU)
Grain size (lm)
29.8 ± 1.7
27.5 ± 4.5
Water temperature (°C)
Salinity
Second year
134.4 ± 23
Water turbidity (NTU)
Grain size (lm)
29.6 ± 1.3
31.4 ± 2.6
Water temperature (°C)
Salinity
First year
S
M
M
S
Downstream
Midstream
Downstream
Upstream
Autumn
Summer
192.8 ± 9
193.0 ± 15
S
200.4 ± 23
182.0 ± 8
S
173.6 ± 8
8.6 ± 3.9
27.0 ± 4.5
26.7 ± 2.4
155.0 ± 17
5.8 ± 1.7
25.8 ± 7.1
26.2 ± 2.3
M
Midstream
162.6 ± 5
6.1 ± 3.6
20.7 ± 5.3
24.5 ± 1.7
150.2 ± 8
5.0 ± 1.2
21.9 ± 8.4
24.3 ± 3.1
M
Midstream
237.2 ± 16
225.2 ± 18
S
227.2 ± 21
225.2 ± 9
S
134.0 ± 5
12.3 ± 3.7
21.8 ± 7.5
26.3 ± 2.7
134.4 ± 12
6.0 ± 1.7
19.9 ± 7.5
25.9 ± 2.7
M
Upstream
134.8 ± 20
7.4 ± 3.3
18.4 ± 7.7
24.2 ± 1.7
121.4 ± 15.6
8.7 ± 1.5
16.0 ± 8.1
23.7 ± 3.3
M
Upstream
181.2 ± 12
188.8 ± 7
S
183.6 ± 18
185.6 ± 8
S
Table 5 Summary of the environmental factors (mean ± standard deviation) in each microhabitat in each station in each season in the first year and second year in a study creek in the Urauchi River mangrove estuary
Fish Sci (2014) 80:21–41 35
123
36
Fish Sci (2014) 80:21–41
-1.0
(a) First year
A. mio
0.5
Turb Crab
P. arg P. tae Gs N. tae L. ful
L. arg P. pal P. rao
Amphi
P. rai
-1.0
-0.5
0
De
Po
Sal A. amb
Algae
CCA 2
Temp Zp
-1.0
-0.5
0
1.0
0.5
CCA 1
(b) Second year 1.0
A. mio
0.5
Gs
Crab Sal
0
CCA 2
Algae
Zp
P. arg L. arg N. tae P. tae Temp L. ful A. amb
Turb
Amphi
P. rai
De P. pal
P. rao
1.0
-0.5
Po
-1.0
-0.5
0
0.5
1.0
CCA 1
Fig. 6 Canonical correspondence analysis (CCA) ordination diagrams based on the densities of ten dominant fish species in the mangrove-root and bare sand areas at each station in each season in the first and second years; first axis is horizontal, second axis is vertical. Environmental factors and food resources represented by vectors. Temp water temperature, Sal salinity, Turb water turbidity, Gs median grain size, Algae macroalgae, Zp zooplankton, Crab crabs, Amphi gammaridean amphipods, Po polychaetes, De detritus. Black symbols, benthic crustacean feeders; gray symbols, zooplankton feeders; white symbols, plant feeders. Species abbreviations given in Table 1
more significant in the habitat choice of small planktivorous fishes within a mangrove system. The sheltering effect of habitat structural complexity against strong water movement and/or predators is sometimes considered as an important factor influencing fish distribution [7, 8, 34, 35]. In the present study, plankton
123
feeders comprised mainly hover gobies, such as Parioglossus raoi and P. palustris, hovering in and around the root structure, into which they usually retreated (indicating a sheltering effect) when the water current became strong or predatory fishes approached (Nanjo K, pers. obs.). Because hover gobies are small (up to ca. 40 mm TL) and are likely to have relatively weak swimming ability, therefore being highly vulnerable to predation [9], such a sheltering effect of root structure is likely to be essential. The CCA results supported this conclusion. In the ordination, hover gobies showed a positive correlation with detritus and water turbidity, and negative with sediment particle size. A greater amount of detritus, higher water turbidity and smaller sediment particles are characteristic of areas with weak water movement [36]. In fact, hover gobies were most abundant in the upstream area of the present study site, where water movement was attenuated. The sheltering effect may therefore be the most important determinant of their distribution pattern, compensating for the disadvantage of a lower food supply. In summary, the present fish assemblages mainly comprised mangrove-root associated fishes, resulting in greater fish species richness and abundance in mangrove-root areas. Food availability may be responsible for the distribution patterns of some mangrove-associated fishes, such as benthic crustacean feeders and plant feeders, which were restricted mostly to the mangrove-root areas, but not so for other types of mangrove fishes, such as zooplankton feeders. For the latter, other possible factors, such as the sheltering effects of mangrove root structure against water current and/or predators, may be more significant. Accordingly, the fish response to factors provided by mangrove vegetation, such as food availability and sheltering effects, may be species-specific, not to mention region-specific, since fish assemblage structures (e.g., species compositions and trophic structures) may differ greatly among regions [23]. Therefore, a small-scale approach, as employed here, should contribute to a better understanding of the relationships between mangrove habitat and resident fishes, in addition to being useful for designing fishery management strategies or conservation policies suitable to a targeted region. Acknowledgments We are grateful to Ken Sakihara, Akira Mizutani and the Okinawa Regional Research Center, Tokai University, for assistance with the fieldwork. Constructive comments on the manuscript from Graham Hardy and two anonymous reviewers were much appreciated. This study was funded by a Grant-in-Aid for Scientific Research (B) from the Japan Society for the Promotion of Science (No. 21380121).
Appendix See Table 6.
F B B
Ambassis sp.
Epinephelus ongus
Caranx ignobilis
Caranx melampygus
Caranx papuensis
Caranx sexfasciatus
Ambassidae
Serranidae
Carangidae
Po Z
Chaetodon auriga
Heniochus acuminatus
Chaetodontidae
Gobiidae B B
Acentrogobius janthinopterus
Acentrogobius moloanus
B
B
Ophiocara porocephala
B
Eleotris melanosoma
Pl
Omox sp.
Butis amboinensis
Pl
Omox biporos
Eleotridae
Pl
Omobranchus ferox
Blenniidae
B
Terapon jarbua
Teraponidae
I
Toxotes jaculatrix
Toxotidae
B
B
Monodactylus argenteus
B
Parupeneus spilurus
B
Sparidae sp.
Parupeneus indicus
B
Acanthopagrus sivicolus
Po
Gerres oyena B
B
Gerres filamentosus
Acanthopagrus pacificus
B
B
Lutjanus russellii
Gerres erythrourus
B
Lutjanus fulviflamma
B
Z
F
Monodactylidae
Mullidae
Sparidae
Gerreidae
Lutjanidae
B
Pterois volitans
Scorpaenidae
F
Strongylura incisa
Belonidae
I
D
Liza melinopterus
Zenarchopterus dunckeri
D
F
Liza macrolepis
F
Evenchelys macrurus
TC
Echidna rhodochilus
Species
Hemiramphidae
Mugilidae
Muraenidae
First year
Family
3–9
6–11
20–22
5
10
3
3–5
2–5
9–23
10
4–5
12–20
1–10
10
8–14
20
11–33
15–30
5–18
8
4–20
6–18
8–20
13–15
8–18
9–15
8–16
12–16
5
15
30
4–10
8–21
8–20
150–200
30
TL
0.8
1.5
0.2
0.8
0.8
1.0
3.3
0.5
0.7
0.2
2.3
2.0
0.2
0.5
1.0
0.2
0.2
0.2
0.2
0.2
1.8
0.2
0.3
0.3
1.3
0.2
0.7
1.0
0.2
0.2
0.2
0.2
S
1.8
0.3
0.3
0.2
0.2
0.2
0.8
0.8
M
0.2
S
0.7
0.2
0.2
0.3
0.8
7.3
0.7
7.7
1.3
1.2
2.8
0.7
0.7
0.2
0.8
1.8
0.2
M
0.3
S
M
M
S
Downstream
Midstream
Downstream
Upstream
Autumn
Summer
0.3
0.2
0.7
0.7
9.0
0.3
1.2
0.3
1.3
2.2
0.2
0.2
0.3
0.8
6.8
0.3
M
S
Midstream
3.8
2.2
0.2
0.2
0.2
2.0
0.5
0.5
0.2
M
0.2
S
Upstream
1.0
0.7
0.2
0.5
1.5
0.2
2.0
0.2
2.2
5.3
3.3
1.3
0.3
0.2
M
0.3
S
Downstream
Winter
0.3
0.3
3.0
0.2
1.5
5.5
0.5
0.7
0.2
0.3
M
0.2
S
Midstream
1.7
0.7
0.5
2.0
0.2
0.2
0.3
0.2
M
0.2
S
Upstream
0.2
0.2
0.7
0.2
0.2
0.7
2.8
0.7
1.3
0.7
0.2
M
S
Downstream
Spring
0.2
0.5
0.3
5.3
0.2
0.2
0.3
0.2
0.5
0.3
0.2
M
S
Midstream
2.7
1.0
1.5
0.3
2.8
2.0
M
S
Upstream
Table 6 Mean individual number of fish species (except ten dominant species) per transect (20 m2, n = 6) in the mangrove-root area and bare sand area at each station in each season in the first and second years
Fish Sci (2014) 80:21–41 37
123
123 Pl B B B B B B D B B B
Oligolepis acutipennis
Oligolepis stomias
Oxyurichthys cornutus
Oxyurichthys ophthalmonema
Oxyurichthys visayanus
Pandaka lidwilli
Pandaka trimaculata
Pseudogobius javanicus
Redigobius balteatus
Vanderhorstia lanceolata
Yongeichthys criniger
B B
Arothron reticularis
Chelonodon patoca
Epinephelus ongus
Caranx melampygus
Serranidae
Carangidae
D
F
B
I
Liza melinopterus
Zenarchopterus dunckeri
Hemiramphidae
D
Liza macrolepis
Mugilidae
F
Evenchelys macrurus
Muraenidae
Second year
B
Arothron manilensis
Tetraodontidae
F
Sphyraena barracuda
Sphyraenidae
Pl
Pl
Siganus guttatus
Acanthurus xanthopterus
Pl
Siganus fuscescens
Z
B
Glossogobius circumspectus
Parioglossus interruptus
F
Glossogobius biocellatus
Z
B
Favonigobius reichei
Parioglossus dotui
B
Exyrias puntang
4–7
B B
B B
Bathygobius fuscus
Cryptocentrus singapurensis
Drombus sp.
3–6 5–8
Pl
Amblygobius phalaena
Eviota prasina
4–10
Pl
Amblygobius linki
9–15
12–16
4–10
8–21
8–20
150–230
12
10
6
18–30
4–10
5–20
5
2–3
1–5
2–11
8
5
2–5
1
1–2
8
5–9
4–10
10
2–10
4–8
1–10
2–7
4–13
6
1–6
4–6
B
Acentrogobius multifasciatus
TL
TC
Species
Acanthuridae
Siganidae
Ptereleotridae
Family
Table 6 continued
0.8
0.7
0.8
0.2
0.2
0.2
2.7
3.5
0.2
5.5
0.2
0.2
0.2
3.2
3.3
0.3
0.7
0.5
0.7
0.5
1.0
0.3
0.3
1.7
0.8
0.3
1.5
3.7
1.5
15.2
1.3
0.3
0.5
2.7
0.2
0.8
0.2
0.3
1.0
0.3
0.2
0.2
0.8
4.8
0.2
0.2
0.3
1.3
7.8
0.8
M
S
0.2
0.2
1.7
0.2
2.0
2.0
1.7
1.7
0.2
0.3
0.2
2.2
0.5
M
0.5
0.2
0.2
S
Downstream
Upstream
M
M
S
Midstream
Downstream S
Autumn
Summer
0.7
0.8
8.5
0.2
0.8
0.2
0.2
1.7
1.2
0.5
2.7
0.2
M
0.2
S
Midstream
0.2
0.2
2.2
0.5
0.3
0.2
4.3
0.5
0.2
0.5
S
Upstream M
0.2
0.2
0.3
0.5
0.2
0.5
0.2
1.2
0.2
M
0.2
0.5
0.3
S
Downstream
Winter
0.2
0.3
0.3
0.7
0.5
0.2
0.3
2.8
0.3
0.2
0.2
S
Midstream M
0.3
0.3
1.0
0.7
1.0
0.5
0.2
0.2
S
Upstream M
0.5
0.2
0.2
3.0
1.7
19.7
0.2
1.0
0.2
0.2
0.2
M
0.7
0.2
S
Downstream
Spring
0.3
0.2
0.3
0.2
0.3
0.7
0.3
23.3
0.7
0.2
0.3
2.2
2.0
0.3
3.8 0.7
0.2
0.2
S
Midstream M
0.2
0.2
0.7
0.5
1.7
1.7
0.7
0.5
0.2
0.2
S
Upstream M
38 Fish Sci (2014) 80:21–41
B
Acanthopagrus sivicolus
Po Z
Chaetodon auriga
Heniochus acuminatus
Chaetodontidae
B B
Eleotris melanosoma
Ophiocara porocephala
Gobiidae
B
Butis amboinensis
Eleotridae
B Pl B B B B F B B Pl B
Amblygobius phalaena
Bathygobius fuscus
Cryptocentrus singapurensis
Exyrias puntang
Favonigobius reichei
Glossogobius biocellatus
Glossogobius circumspectus
Myersina macrostoma
Oligolepis acutipennis
Oxyurichthys cornutus
B
Acentrogobius multifasciatus Pl
B
Acentrogobius moloanus
Amblygobius hectori
B
Acentrogobius janthinopterus
Amblygobius linki
B
Acentrogobius audax
Pl
Omox biporos
Blenniidae
B
Terapon jarbua
Teraponidae
I
Toxotes jaculatrix
Toxotidae
B B
Parupeneus spilurus
Monodactylus argenteus
B
Parupeneus indicus
Monodactylidae
Mullidae
Acanthopagrus pacificus
Sparidae
B B
Lethrinidae sp.
Lethrinidae
B
Po
Plectorhinchus albovittatus
B
B
Lutjanus sp.
Gerres oyena
B
Gerres erythrourus
B
B
Caranx sexfasciatus
Lutjanus russellii
B
Caranx papuensis
Lutjanus fulviflamma
TC
Species
Haemulidae
Gerreidae
Lutjanidae
Family
Table 6 continued
4–10
2–10
4
4–7
1–10
2–7
4–13
7
4
4–10
1–6
6
5
3–9
6–11
4
20–25
5
10
3–5
9–23
10
4–5
12–20
1–10
10
8–14
11–33
15–30
6
12
5–22
4–20
8
6–18
8–20
13–15
8–18
TL
2.2
0.3
0.8
0.7
0.5
0.5
0.7
0.5
1.2
0.2
0.2
0.8
0.8
0.5
3.3
0.5
0.5
2.2
3.7
0.7
0.2
3.2
0.7
0.2
0.2
0.8
0.3
0.2
1.8
0.2
0.7
0.7
0.2
1.8
4.7
0.7
1.8
1.5
0.8
4.3
0.3
0.2
1.7
0.2
1.2
0.2
0.2
0.7
M
0.2
0.2
S
1.2
0.2
1.3
0.8
0.2
1.8
7.3
0.3
3.5
0.2
1.8
0.2
0.2
M
0.2
S
Downstream
Upstream
M
M
S
Midstream
Downstream S
Autumn
Summer
0.2
0.2
0.5
0.3
0.2
0.3
1.8
9.0
1.7
0.3
2.0
2.0
3.5
M
0.2
S
Midstream
0.5
1.0
3.5
0.5
7.2
1.5
0.2
0.2
0.2
1.8
0.5
0.3
0.7
0.2
S
Upstream M
0.2
0.2
0.2
0.2
0.8
0.2
2.7
1.5
1.2
5.2
10.5
3.5
1.2
M
0.2
S
Downstream
Winter
0.3
0.2
1.8
0.2
0.3
0.7
0.3
3.0
0.2
1.5
2.3
1.3
0.3
0.2
0.2
0.2
S
Midstream M
0.2
0.7
0.7
0.2
1.0
0.3
0.2
2.7
1.7
0.3
0.5
0.2
S
Upstream M
0.2
0.5
0.2
0.2
0.2
0.5
0.2
2.8
0.2
0.2
0.2
1.0
0.3
1.2
1.8
M
0.3
0.2
S
Downstream
Spring
1.5
0.5
0.5
2.5
1.3
0.2
5.3
0.2
0.2
0.3
0.2
1.7
0.2
0.2
S
Midstream M
0.7
1.3
0.7
0.2
2.7
0.8
1.3
0.7
1.0
0.2
S
Upstream M
Fish Sci (2014) 80:21–41 39
123
123 D B B B
Pseudogobius javanicus
Redigobius balteatus
Yongeichthys criniger
Gobiidae sp.
B B
Arothron manilensis
Chelonodon patoca
Tetraodontidae 14
6
18–30
4–10
5–20
5
8
2–3
1–5
4
2–11
3
2–5
1–2
8
5–10
TL
0.2
1.8
0.8
0.2
0.2
5.3
0.2
0.2
2.0
0.3
1.0
2.2
3.3
0.7
3.2
7.8
2.2
M
S
0.2
0.2
1.2
3.5
0.2
1.5
M
S
Downstream
Upstream
M
M
S
Midstream
Downstream S
Autumn
Summer
0.2
0.7
0.2
0.2
0.2
M
S
Midstream
0.2
5.3
0.5
0.5
S
Upstream M
2.5
1.7
1.2
0.2
M
S
Downstream
Winter
0.5
6.7
0.3
0.2
S
Midstream M
2.8
0.2
0.3
0.5
1.0
S
Upstream M
0.3
1.7
3.8
0.5
0.2
M
S
Downstream
Spring
0.7
0.2
0.3
46.7
0.8
0.2
1.3
S
Midstream M
0.7
0.7
4.2
2.0
0.7
S
Upstream M
TC trophic category, TL total length (cm), M mangrove-root area, S bare sand area. Trophic category; B benthic crustacean feeders, Z zooplankton feeders, Pl plant feeders, Po polychaete feeders, F fish feeders, D detritus feeders, I insect feeders
F
Sphyraena barracuda
Sphyraenidae
Pl
Acanthurus xanthopterus
Pl
Siganus guttatus
Acanthuridae
Pl
Siganus fuscescens
D
Scatophagus argus
Siganidae
Z
B
Pandaka lidwilli
Z
B
Oxyurichthys visayanus
Parioglossus interruptus
B
Oxyurichthys ophthalmonema
Parioglossus dotui
TC
Species
Scatophagidae
Ptereleotridae
Family
Table 6 continued
40 Fish Sci (2014) 80:21–41
Fish Sci (2014) 80:21–41
References 1. Robertson AI, Blaber SJM (1992) Plankton, epibenthos and fish communities. In: Robertson AI, Alongi DM (eds) Tropical mangrove ecosystems. American Geophysical Union, Washington DC, pp 173–224 2. Blaber SJM (2000) Tropical estuarine fishes: ecology, exploitation, and conservation. Blackwell, Oxford 3. Kathiresan K, Bingham BL (2001) Biology of mangroves and mangrove ecosystems. Adv Mar Biol 40:81–251 4. Nagelkerken I, Blaber SJM, Bouillon S, Green P, Haywood M, Kirton LG, Meynecke JO, Pawlik J, Penrose HM, Sasekumar A, Somerfield PJ (2008) The habitat function of mangroves for terrestrial and marine fauna: a review. Aquat Bot 89:155–185 5. Nagelkerken I (2009) Evaluation of nursery function of mangroves and seagrass beds for tropical decapods and reef fishes: patterns and underlying mechanisms. In: Nagelkerken I (ed) Ecological connectivity among tropical coastal ecosystems. Springer, New York, pp 357–399 6. Manson FJ, Loneragan NR, Skilleter GA, Phinn SR (2005) An evaluation of the evidence for linkages between mangroves and fisheries: a synthesis of the literature and identification of research directions. Oceanogr Mar Biol Annu Rev 43:483–513 7. Laegdsgaard P, Johnson C (2001) Why do juvenile fish utilise mangrove habitats? J Exp Mar Biol Ecol 257:229–253 8. Nanjo K, Nakamura Y, Horinouchi M, Kohno H, Sano M (2011) Predation risks for juvenile fishes in a mangrove estuary: a comparison of vegetated and unvegetated microhabitats by tethering experiments. J Exp Mar Biol Ecol 405:53–58 9. Sogard SM (1992) Variability in growth rates of juvenile fishes in different estuarine habitats. Mar Ecol Prog Ser 85:35–53 10. Grol MGG, Nagelkerken I, Rypel AL, Layman CA (2011) Simple ecological trade-offs give rise to emergent cross-ecosystem distributions of a coral reef fish. Oecologia 165:79–88 11. Woodroffe C (1992) Mangrove sediments and geomorphology. In: Robertson AI, Alongi DM (eds) Tropical mangrove ecosystems. American Geophysical Union, Washington DC, pp 7–41 12. Zuccarello GC, Yeates PH, Wright JT, Bartlett J (2001) Population structure and physiological differentiation of haplotypes of Caloglossa leprieurii (Rhodophyta) in a mangrove intertidal zone. J Phycol 37:235–244 13. Verweij MC, Nagelkerken I, de Graaff D, Peeters M, Bakker EJ, van der Velde G (2006) Structure, food and shade attract juvenile coral reef fish to mangrove and seagrass habitats: a field experiment. Mar Ecol Prog Ser 306:257–268 14. Camilleri JC (1992) Leaf-litter processing by invertebrates in a mangrove forest in Queensland. Mar Biol 114:139–145 15. Kon K, Kurokura H, Hayashizaki K (2007) Role of microhabitats in food webs of benthic communities in a mangrove forest. Mar Ecol Prog Ser 340:55–62 16. Jones GP (1986) Food availability affects growth in a coral reef fish. Oecologia 70:136–139 17. Wilson DT, Meekan MG (2002) Growth-related advantages for survival to the point of replenishment in the coral reef fish Stegastes partitus (Pomacentridae). Mar Ecol Prog Ser 231:247–260 18. Meekan MG, von Kuerthy C, McCormick MI, Radford B (2010) Behavioural mediation of the costs and benefits of fast growth in a marine fish. Anim Behav 79:803–809
41 19. Smith TM, Hindell JS (2005) Assessing effects of diel period, gear selectivity and predation on patterns of microhabitat use by fish in a mangrove dominated system in SE Australia. Mar Ecol Prog Ser 294:257–270 20. Payne NL, Gillanders BM (2009) Assemblages of fish along a mangrove-mudflat gradient in temperate Australia. Mar Freshw Res 60:1–13 21. Wang M, Huang Z, Shi F, Wang W (2009) Are vegetated areas of mangroves attractive to juvenile and small fish? The case of Dongzhaigang Bay, Hainan Island, China. Estuar Coast Shelf Sci 85:208–216 22. Tse P, Nip THM, Wong CK (2008) Nursery function of mangrove: a comparison with mudflat in terms of fish species composition and fish diet. Estuar Coast Shelf Sci 80:235–242 23. Faunce CH, Serafy JE (2006) Mangroves as fish habitat: 50 years of field studies. Mar Ecol Prog Ser 318:1–18 24. Nakabo T (2002) Fishes of Japan with Pictorial Keys to the Species, English edn. Tokai University Press, Tokyo 25. Nanjo K, Kohno H, Sano M (2008) Food habits of fishes in the mangrove estuary of Urauchi River, Iriomote Island, southern Japan. Fish Sci 74:1024–1033 26. Underwood AJ (1981) Techniques of analysis of variance in experimental marine biology and ecology. Oceanogr Mar Biol Annu Rev 19:513–605 27. terBraak CFJ (1986) Canonical correspondence analysis: a new eigenvector technique for multivariate direct gradient analysis. Ecology 67:1167–1179 28. terBraak CFJ, Verdonschot PFM (1995) Cannonical correspondence analysis and related multivariate methods in aquatic ecology. Aquat Sci 57:255–289 29. Connolly RM (1994) Removal of seagrass canopy: effects on small fish and their prey. J Exp Mar Biol Ecol 184:99–110 30. Horinouchi M (2007) Distribution patterns of benthic juvenile gobies in and around seagrass habitats: effectiveness of seagrass shelter against predators. Estuar Coast Shelf Sci 72:657–664 31. Wen CKC, Pratchett MS, Almany GR, Jones GP (2013) Role of prey availability in microhabitat preferences of juvenile coral trout (Plectropomus: Serranidae). J Exp Mar Biol Ecol 443:39–45 32. Horinouchi M, Mizuno N, Jo Y, Fujita M, Suzuki Y, Aranishi F, Sano M (2013) Habitat preference rather than predation risk determines the distribution patterns of filefish Rudarius ercodes in and around seagrass habitats. Mar Ecol Prog Ser 488:255–266 33. Robertson AI, Dixon P, Daniel PA (1988) Zooplankton dynamics in mangrove and other nearshore habitats in tropical Australia. Mar Ecol Prog Ser 43:139–150 34. Almany GR (2004) Does increased habitat complexity reduce predation and competition in coral reef fish assemblages? Oikos 106:275–284 35. Horinouchi M (2007) Review of the effects of within-patch scale structural complexity on seagrass fishes. J Exp Mar Biol Ecol 350:111–129 36. Chambers PA, Prepas EE, Hamilton HR, Bothwell ML (1991) Current velocity and its effect on aquatic macrophytes in flowing waters. Ecol Appl 1:249–257
123