Fisheries and Aquatic Sciences
The Korean Society of Fisheries and Aquatic Science
RESEARCH ARTICLE

Diversity and abundance of the fish community in Elala and Gereb Tsedo streams, Tekeze Sub-Basin: Northern Ethiopia

Solomon Tesfay1,*https://orcid.org/0000-0001-7760-0045, Mekonen Teferi1,*https://orcid.org/0000-0002-1868-3319, Tsegazeabe Hadush Haileselasie1https://orcid.org/0000-0003-4567-620X
1Department of Biology, CNCS, Mekelle University, Mekelle 231, Tigray, Ethiopia
*Corresponding author: Mekonen Teferi, Department of Biology, CNCS, Mekelle University, Mekelle 231, Tigray, Ethiopia, Tel: +251-914708783, E-mail:mekonen.teferi@mu.edu.et
*Corresponding author: Solomon Tesfay, Department of Biology, CNCS, Mekelle University, Mekelle 231, Tigray, Ethiopia, Tel: +251-914110672, E-mail:Solomon232000@yahoo.com

Copyright © 2024 The Korean Society of Fisheries and Aquatic Science. This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Jan 15, 2024; Revised: Mar 11, 2024; Accepted: May 31, 2024

Published Online: Aug 31, 2024

Abstract

Fishes are found in different lentic and lotic aquatic ecosystems. Determining the composition of fish communities is a crucial component of water quality evaluations and conservation strategies in flowing rivers. The present study is carried out to investigate fish species diversity, richness, and abundance in two second-order streams: Gereb Tsedo (GTS) and Elala in the Tekeze sub-basin. We collected fish with backpack electrofishing during August and December 2013 and March 2014. Fish species richness and the total number of species per stream and per microhabitat were calculated for each stream. Species diversity for each site was calculated using the Shannon index of diversity using the standard formula H’ = –Σpi ln (pi). Diversity data among habitats and streams were compared using one-way analysis of variance (ANOVA) and Tukey post hoc comparisons, using the statistical software package STATISTICA 11. A total of 6,554 specimens representing four species: Garra blanfordii, Grra ignestii, Garra dembecha and Garra aethiopica, were caught. They are all indigenous riverine fish species of the family Cyprinidae. Fish abundance was significantly higher in GTS than Elala stream (t-test, t = 1.444, df = 3, p < 0.05). Across the two streams, the overall relative fish abundance was higher in pools, with 53%, followed by runs and riffles with 35% and 12%, respectively. Garra blanfordii was dominant in pools while G. aethiopica was the least dominant. Stream and habitat types were likely the more important driving factors behind the observed patterns of diversity, although the impact of other, not evaluated variables cannot be excluded. Most likely, the observed patterns do not represent a signature of fish dispersal limitation given that almost all species were widely distributed across the study area.

Keywords: Abundance; Cyprinidae; Diversity; Garra; Stream

Introduction

Fishes are the most diverse of any vertebrate group and teleosts diversity in particular constitutes an assemblage that is larger than all other vertebrates put together. Fish, in a sense, do practically everything. They have experimented in an evolutionary sense with nearly every kind of adaptation that can be seen modified in more derived vertebrates (Kumar, 2013). These adaptations may be structural, physiological, as well as behavioral (Kumar, 2013; Morrow & Fischenich, 2000).

Maintaining the species diversity of fish that live in streams is one of the streams’ ecological goals. Identifying patterns of variation in stream-dwelling fish assemblages, as well as potential contributory reasons, is a major issue in stream ecology (Godinho et al., 2000; Koel & Peterka, 2003). The ability to properly estimate ensemble attributes (e.g., species richness, composition, relative abundance) has crucial consequences for the management, and conservation of ecological assemblages (Godinho et al., 2000).

Changes in species richness, relative abundance, and species composition in aquatic environments have been thought to imply a loss of biodiversity or biotic integrity (Cheng et al., 2012), and fish assemblages have been recognized as reliable indicators in reflecting aquatic ecosystem health (Ibarra et al., 2003; Rashleigh, 2004) or habitat quality such as flow regulation, habitat degradation (like physical habitat alteration and fragmentation), environmental contamination (e.g., eutrophication, acidification, and chemical pollution) (Vila-Gispert et al., 2002). In addition, in fisheries research, there is often a need to sample fish from natural environments. This may be for monitoring (e.g., traditional fisheries assessments including determining species abundance, diversity, and population structure) or to collect fish for research (e.g., for implantation of telemetry tags, collection of fish for lab studies). As a result, fishes have received considerable attention in biomonitoring, characteristics of their environment, and other forms of aquatic life (Vila-Gispert et al., 2002).

In addition, the characterization of fish community composition is an important aspect of conservation plans and water quality assessments in flowing waters. Fish community composition and abundance are sensitive to a variety of environmental perturbations. Efforts to quantify these parameters involve the creation of appropriate sampling procedures or standards which themselves are often heavily influenced by considerations of manpower and budget (Kimmel & Argent, 2006).

It is a common view that freshwater fish live not in random groupings but in structured communities held together by favorable abiotic and biotic mechanisms. In some cases, a small number of local environmental variables such as temperature, dissolved oxygen, pH, electrical conductivity, etc., seem to exercise a strong influence on community structure, while in others it is related to a wider range of factors such as geographical factors like altitude, surface area, etc. (Edds, 1993). The structure of fish assemblages may also be a result of stochastic factors such as disturbance or deterministic interactions like competition and predation (Grossman & Freeman, 1987).

Knowledge of the diversity, distribution, and relative abundance of fish in the streams is of vital importance to ensure better water quality management and conservation activities. The Ethiopian highland streams and rivers seem to be a center for fish diversity on the African continent (Stiassny & Getahun, 2007). Few taxonomic studies in these water systems showed that Garra was the dominant genera (Stiassny & Getahun, 2007). Moreover, artificial micro dams constructed in Tigray were also inhabited mainly by Garra species (Asmelash et al., 2007; Dejenie et al., 2008; Teferi et al., 2013, 2014). As to our current knowledge, there is no study on the diversity, distribution, and abundance of fish communities on the highland streams or rivers of Tigray. Therefore, the present study is an attempt: (1) To document the freshwater fish species abundance and relative abundance in the study streams, (2) To compare the freshwater fish species richness and diversity in the study streams, (3) To compare the fish abundance, richness, and diversity of Elala and Gereb Tsedo (GTS) streams with their confluence sites.

Materials and Methods

Description of study area

The current study was conducted in two streams Elala and GTS that pass through Mekelle city (Fig. 1). The Tigray National Regional State, located in northern Ethiopia, has Mekelle as its capital. The city’s North Western and South Western sections are traversed by the sporadic streams GTS and Elala, respectively. People use these stream waters for washing clothes and for livestock drinking. Moreover, these streams provide several services to the community such as household use, irrigation, production of building materials (bricks), and carwash activities along their gradients from the headwater as they pass down through the city. These two streams join together and make Mariam Dahan stream. This common stream further continues down to make Romanat stream (Fig. 1).

fas-27-8-501-g1
Fig. 1. Map of study area showing Mekelle City and sampling sites indicated with triangular shapes along the streams. A, Aragure; B, May weyni primary school; C, Mariam Dahan(Confluence); D, Mariam Dahan stream; RF, Romanat fall (Stream).
Download Original Figure

Elala’s headwaters are located in Aragure, about 16 kilometers from Quiha (Fig. 1). The second stream (GTS) is known by several names as it flows through the city. The initial source and headwater of this stream rises from the compounds of May Weini complete primary school (Fig. 1) and progresses to Catholic and Adventist schools, where it is known as May Tsaeda Egam. After that it continues to run downhill on the west side of Abraha Castle (CAS site) and named May Gifaf before it reaches Kebele 11. Then further continues downward and is named GTS. It continues its way with some small tributaries joining it with wastewater discharged from the town and finally joins Elala stream at Mariam Dahan (Fig. 1).

Data collection

Fish specimens were collected using standardized backpack electrofishing: depletion sampling techniques as described in Zippin (1958) and (“Zippin” method) (SFCC, 2007; Zippin, 1958). Electrofishing from the two study streams and Mariam Dahan, where the two streams meet, was carried out during the third week of August, December 2013, and March 2014. Before their confluence, each stream had six sampling locations, and once they united, they each had two common sites. Using global positioning system (GPS) data, sites were identified and situated at an average distance of 0.65 km between two successive sampling locations, taking into account the streams’ wadeability and accessibility. The stretch of each sampling site (having 3 habitats) was approximately 50 to 100 meters. The time for each run was 5 minutes, and the time gap between the runs was 20 minutes. A consistent sampling design for each site in each sampling period was applied to avoid biased results (Jha et al., 2005). Samples were collected in a particular stream at the scheduled time for all study periods along the same reach, which included all major microhabitats of the stream: pools, runs, and riffles (Jayaratne & Surasinghe, 2010; Jones et al., 2003). Pools are deep with slow water. Runs are deep with fast water and little or no turbulence. Riffles are shallow with fast, turbulent water flowing over rocks (Cushing & David Allan, 2001). Each of the two study streams had six sites with three habitats (pools, runs, and riffles) and two sample runs (Jayaratne & Surasinghe, 2010; Jones et al., 2003). In addition, two sites were visited downstream of the confluence, each containing three habitats and two streams. The total number of sampling sites upstream and downstream of the confluence of both streams was 14 (6 upstream of the confluence and two downstream).

Fish specimens were collected and identified using keys and descriptions of Stiassny & Getahun (2007). Specimens of every species were tallied. Five to ten specimens of unknown species were sent to the lab for species validation and confirmation after being stored in 10% buffered formalin.

Fish species richness and the total number of species per stream and per microhabitat were calculated for each stream (Magurran, 1988). Species diversity for each site was calculated using the Shannon index of diversity using the standard formulae: (Kwak & Peterson, 2007; Shannon & Weaver, 1949). H’ = –Σpi ln(pi), where pi = ni/N; ni is the number of individuals of ‘i’th species and N = Σni. Then these biodiversity indices were used to compare the pattern of fish composition among the study sites and between microhabitats within the stream.

At each sampling site, the following physicochemical parameters such parameters measured: water flow (velocity, m/s), water depth (cm), water temperature (°C), pH, conductivity (µS/cm), and dissolved oxygen (mg/L). Stream velocity was measured using the floatation method of Gordon et al. (1992). Conductivity, pH, oxygen concentration, and water temperatures were measured in-situ using a conductivity meter (SX713, SANXINS, Shenyang, China), pH meter (pH-013, HINOTEK, Ningbo, China), and oxygen meter (HQ40d multimeter, Hach, Loveland, CO, USA), respectively. Besides, other parameters like water transparency are measured using Snell’s tube (diameter 6 cm) (Sovell et al., 2000; Tesfay et al., 2019; Van de Meutter, 2005; Van de Meutter et al. ,2007). Turbidity and chlorophyll-a concentrations (as a proxy of phytoplankton biomass) were measured using a fluorometer (AquaFluor 8000-001, Turner Designs, San Jose, CA, USA) in the field. In addition, macrophyte cover (%) and shading by riparian (%) were estimated following the protocol of van de Meutter (2005).

Data analysis

Diversity data among habitats and streams were compared using one-way analysis of variance (ANOVA) and Tukey post hoc comparisons, using the statistical software package STATISTICA 11 (Statsoft, Tulsa, OK, USA). The normality and homogeneity were tested using the Sharpiro-Wilk and Levene’s tests, respectively. In addition, parametric t-tests were used to compare variables like species diversity, richness, evenness, total fish abundance, and fish species abundance within the stream and between the streams.

Composition of the fish community based on catch per unit of effort, that is to say, individuals/run (King, 1995) of each species of overall streams as well as from each separate stream was analyzed using species abundance and total fish abundance data. Two-way ANOVA was used to test the effect of season, habitat, and their interactions on fish abundance in the streams for temporal (season) and habitat differences.

Results

Physico-chemical and environmental characteristics

The physicochemical factors of the three streams tested (Table 1) show that dissolved oxygen concentration, pH, Chlorophyll a, and macrophyte coverage (%) were relatively high in Elala compared to the other two streams. Furthermore, in December, Elala stream site E2 had a maximum dissolved oxygen concentration of 14.6 mg/L and GTS stream site had a minimum dissolved oxygen concentration of 1.7 mg/L. During the study periods, the pH concentration fluctuated from 6.3 to 8.13. A minimum of 6.3 was reported at GTS stream site GTS-Rn (GTS run), and a maximum of 8.13 was observed at Elala (E2). In addition, the surface water temperature varied from 12.4°C to 26.8°C in Elala, 17°C to 25.1°C in GTS, and 19.3°C to 26.7°C during the study period.

Table 1. Physicochemical and environmental (habitat) characteristics of the streams Gereb Tsedo, Elala and the common downstream of the two streams, Mariam Dahan and Romant fall
Variable GerebTsedo (Ave. ± SD) Elala (Ave. ± SD) Mariam Dahan (Ave. ± SD)
Altitude (m) 2,071.6 ± 31.1 2,011.5 ± 8.01 1,899.5 ± 79.9
Water depth (cm) 39.0 ± 30.4 44.3 ± 34.47 40.1 ± 29.7
River bed width (m) 2.7 ± 1.6 3.5 ± 2.0 4.1 ± 3.0
Water transparency (Snell, cm) 20.8 ± 10.2 19.73 ± 9.73 19.5 ± 10.3
pH 7.4 ± 0.4 7.72 ± 0.28 7.5 ± 0.4
Conductivity (µS/cm) 142.6 ± 17.4 115.1 ± 22.8 106.3 ± 4.1
Chlorophyll a (µg/L) 57.3 ± 20.8 73.3 ± 26.2 62.2 ± 21.04
Dissolved oxygen (mg/L) 5.4 ± 1.9 7.1 ± 2.2 5.3 ± 1.74
Temperature (°C) 20.8 ± 1.7 21.5 ± 3.5 21.3 ± 2.6
Water flow velocity (m/s) 0.2 ± 0.1 0.127 ± 0.11 48.4 ± 34.0
% Macrophyte coverage 14.76 ± 16.17 18.36 ± 23.43 13.3 ± 8.1
% Riparian coverage 58.7 ± 29.0 33.4 ± 31.5 48.4 ± 34.0

Physicochemical and environmental characteristics were taken from 14 study sites; each site with 3 habitats (42 microhabitats × 3 seasons) of the streams and the average value and SD is presented.

Download Excel Table
Fish abundance between and within streams and among habitats

We collected a total of 8,410 individuals that belong to eight different fish species of the family Cyprinidae: Garra aethiopica, Garra blanfordii, Garra dembecha, Garra ignestii, Varicorhinus beso, Labeobarbus intermedius, Labeo niloticus, and Labeo forskalii (Table 3; Fig. 2). Among these, the first four Garra species, except G. aethiopica were recorded from the two study streams. However, G. aethiopica was detected only in GTS. Both streams account for about 6,554 (77.9%) of the total fish population or individuals collected during the study periods.

fas-27-8-501-g2
Fig. 2. Fish species recorded from the study streams. GB, Garra blanfordii; GI, Garra ignestii; GD, Garra dembecha; GA, Garra aethopica; LBI, Labeobarbus intermedius; LN, Labeo niloticus; LF, Labeo forskalii; VB, Varicorhinus beso.
Download Original Figure

In GTS the total fish abundance (4,870; 74.3%) was significantly higher than Elala stream (1,684; 25.7%) (t-test, t = 1.444, df = 3, p < 0.05; Table 2). The total catch of G. blanfordii and G. ignestii in GTS were also significantly higher than that of Elala (Fig. 3; Table 3). There were no significant differences in the abundances of each of the other two species (G. aethiopica and G. dembecha) between the two streams (Fig. 3; Table 3). Comparison of fish species abundance within stream showed that G. blanfordii are dominant in GTS stream (Fig. 4A), while there were no significant differences among the species of Elala stream (Fig. 4B). There were no catch of G. aethiopica in Elala and lower catch in only two sites of GTS stream (Figs. 3 and 4; Table 2). Moreover, this species was not recorded from their confluence point (first site of Mariam Dahan stream) at Mariam Dahan church (Table 4). However, G. aethiopica was also recorded from the second site of Mariam Dahan which is from Romant stream (particularly below Romant Fall).

Table 2. t-test results comparing each fish species abundance between Gereb Tsedo (GTS) and Elala streams
Response variable GTS (Mean ± SE) Elala (Mean ± SE) t-value df p-value
Garra aethiopica 7.17 ± 4.93 0.00 ± 0 1.45 10 0.177
Garra blanfordii 601.00 ± 260.27 200.50 ± 42.09 1.52 10 0.039*
Garra dembecha 75.17 ± 30.73 60.33 ± 23.80 0.38 10 0.710
Garra ignestii 128.33 ± 50.07 19.83 ± 15.96 2.06 10 0.043*
Total fish abundance 811.67 ± 309.55 280.67 ± 62.77 3.20 10 0.0084**

p < 0.05;

p <0.01.

Download Excel Table
Table 3. Total fish abundance recorded from Gereb Tsedo, Elala and Mariam Dahan streams (Site wise), Tekeze sub-basin, northern Ethiopia
Fish species GTS ELA MD
MW MTE CAS MG K11 GTS E1 E2 E3 E4 E5 E6 MD RF Total
GA 0 14 0 0 29 0 0 0 0 0 0 0 0 33 76
GB 103 456 1,804 462 713 68 212 329 50 137 178 297 566 634 6,009
GD 34 40 47 120 207 3 167 90 19 31 26 29 145 132 1,090
GI 52 77 259 56 308 18 99 0 9 0 11 0 122 140 1,151
VB 0 0 0 0 0 0 0 0 0 0 0 0 0 28 28
LBI 0 0 0 0 0 0 0 0 0 0 0 0 0 25 25
LN 0 0 0 0 0 0 0 0 0 0 0 0 0 12 12
LF 0 0 0 0 0 0 0 0 0 0 0 0 0 19 19
Total 189 587 2,110 638 1,257 89 478 419 78 168 215 326 833 1,023 8,410

The sites are arranged according to their increasing distance from the source.

GA, Garra aethiopica; GB, Garra blanfordii; GD, Garra dembecha; GI, Garra ignestii; VB, Varicorhinus beso; LBI, Labeobarbus intermedius; LN, Labeo niloticus; LF, Labeo forskalii; GTS, Gereb Tsedo; ELA, Elala; MD, Mariam Dahan; MW, May Weini; MTE, May Tsaeda Egam; CAS, Castle; MG, May Gifaf; K11, Kebele 11; GTS, Gereb Tsedo; E1, Feleg Daero; E2, Elala2; E3, Elala3; E4, Elala4; E5, Elala5; E6, Elala6; RF, Romanat Fall.

Download Excel Table
fas-27-8-501-g3
Fig. 3. Total number of fish individuals per species caught in Gereb Tsedo and Elala streams. Error bars indicate standard errors. Bars labeled by a different letter are significantly different, p < 0.05; these statistical comparisons were only made between the abundances of the same species from two streams.
Download Original Figure
fas-27-8-501-g4
Fig. 4. Fish species abundance. Total number of fish species caught in Gereb Tsedo (a), Elala stream (b). Error bars indicate standard errors. Bars labeled by a different letter are significantly different (p < 0.05). These statistical comparisons were made among species abundances within the same stream.
Download Original Figure

In the two streams, the highest total relative fish abundance was recorded in pools followed by runs and riffles (Table 4). Species-wise, G. blanfordii was dominant which account for about 71% of the pool fish community recorded in the two streams followed by G. ignestii and G. dembecha with 15% and 13% of the pool fish community, respectively. Where the least abundant in pool was G. aethiopica only accounted for 1% of the pool fish community (Table 5). Overall, G. blanfordii showed significantly higher numerical abundance and dominated the fish community in pools (Table 5).

Table 4. Fish species abundance and assemblage attributes (total abundance, richness, diversity and evenness) for different habitats sampled in streams Elala and GerebTsedo
Fish species Pool Run Riffle Total RSA (%)
Garra aethiopica 41 2 0 43 0.7
Garra blanfordii 2,464 1,772 573 4,809 73.4
Garra dembecha 444 197 172 813 12.4
Garra ignestii 493 341 55 889 13.5
Total 3,442 2,312 800 6,554 100
RA (%) 53 35 12 100 NA
Richness (S) 4 4 3 4 NA
Shannon (H’) 0.835 0.702 0.754 0.79 NA
Evenness (J) 0.576 0.505 0.708 0.551 NA

RSA, relative species abundance; RA, relative abundance; NA, not applicable.

Download Excel Table
Table 5. Fish species diversity indices (Shannon and Weaver species Diversity) (H’), relative abundance (RA), species richness (S) and species evenness (J) of each sampling site in Gereb Tsedo and Elala streams, Mekelle city, northern Ethiopia
Stream Code Pool Run Riffle RA
H' S J H' S J H' S J
GTS MW 0.895 3 0.816 1.004 3 0.915 0.990 3 0.903 2.3
MTE 0.890 4 0.609 0.567 3 0.588 0.728 3 0.691 7.0
CAS 0.437 3 0.516 0.484 3 0.541 0.435 3 0.515 25.1
MG 0.762 3 0.714 0.709 3 0.677 0.839 3 0.771 7.6
K11 1.08 4 0.736 0.766 4 0.538 0.942 3 0.855 15.0
GTS 0.480 2 0.808 0.652 2 0.960 0.451 3 0.785 1.1
Mean (SD) 0.757 (0.25) 3.2 (0.75) 0.641 (0.12) 0.697 (0.18) 3.0 (0.63) 0.703 (0.19) 0.731 (0.24) 3.0 (0.00) 0.753 (0.14) 9.7 (9.01)
Elala E1 0.720 3 0.685 0.777 3 0.725 0.673 2 0.980 5.7
E2 0.356 2 0.714 0.651 2 0.959 0.679 2 0.986 5.0
E3 0.815 3 0.753 0.731 3 0.692 0.668 2 0.975 0.93
E4 0.418 2 0.759 0.500 2 0.825 0.691 2 0.998 2.0
E5 0.456 3 0.526 0.586 2 0.898 0.621 2 0.931 2.6
E6 0.171 2 0.593 0.447 2 0.782 0.530 2 0.850 3.9
Mean (SD) 0.489 (0.24) 2.5 (0.55) 0.672 (0.09) 0.615 (0.13) 2.33 (0.52) 0.814 (0.10) 0.644 (0.06) 2.0 (0.00) 0.953 (0.06) 3.4 (1.83)
Below junction MD 0.695 3 0.668 0.543 3 0.574 0.988 3 0.895 9.9
RF 1.401 8 0.507 1.046 8 0.356 0.930 5 0.507 12.2
Mean (SD) 1.048 (0.50) 5.500 (0.52) 0.588 (0.11) 0.795 (0.36) 5.500 (3.56) 0.465 (0.15) 0.959 (0.04) 4.0 (1.41) 0.701 (0.27) 11.05 (1.63)

GTS, Gereb Tsedo; MW, May Weini; MTE, May Tsaeda Egam; CAS, Castle; MG, May Gifaf; K11, Kebele 11; E1, Feleg Daero; E2, Elala 1; E3, Elala 2; E4, Elala3; E5, Elala4; E6, Elala 5; MD, Mariam Dahan; RF, Romanat Fall.

Download Excel Table

There was no significant seasonal difference in fish abundance but a pattern of increase in total relative fish abundance was observed between sampling periods (Two-way ANOVA).

Fish species diversity and richness between streams and among habitats

Fish species richness varied from 3–4 in the streams with the highest being in GTS stream (Tables 2 and 3; Fig. 5A and 5B) where the number of sample sites was the same in both streams. Species richness streams’ habitat ranges from 2–4. The highest fish species richness was recorded in pools. We have observed a significant difference in fish species richness between riffles of the two streams (t-test, t = 7, df = 5, p < 0.001, Fig. 5B).

fas-27-8-501-g5
Fig. 5. Fish species richness. (A) Fish species richness (S) of total number of fish species between the two streams; comparison was done. (B) Fish species richness (S) of total number of fish species among habitat types between the two streams; comparison was done. Error bars indicate standard errors. Bars labeled by a different letter are significantly different (p < 0.05). These statistical comparisons were made between similar habitats of the streams species.
Download Original Figure

The relative occurrence of the fish varies between sampling sites and among streams. For example, G. aethiopica, was found in two sites, while G. dembecha, G. blanfordii and G. ignestii are very widespread in GTS. Although widespread, G. dembecha and G. ignestii were more frequently found at relatively low abundances (Table 3; Fig. 2). Meanwhile G. blanfordii and G. dembecha were observed in all the study sites of Elala but G. ignestii occurred in 50% of the sampling sites of Elala stream (Table 3; Fig. 2).

There were no statistically significant differences in species diversity (Shannon-Weaver Diversity, H’) between streams as well as between habitat types (t-test, t = 1.095, df = 5, p < 0.323; Fig. 6A, Table 5 and between habitats types; pools t = 1.475, df = 5, p < 0.200, runs t = 0.403, df = 5, p < 0.704, riffles t = 0.964, df = 5, p < 0.379, Fig. 6B, Table 5).

fas-27-8-501-g6
Fig. 6. Fish species diversity. (A) Shannon-Weaver species Diversity (H') of total number of fish species between the two streams; comparison was done. (B) Shannon-Weaver species Diversity (H') of total number of fish among habitat types between the two streams; comparison was done. Bars labeled by the same letter are not significantly different, p < 0.05; these statistical comparisons were made between similar habitats of the streams.
Download Original Figure
Fish abundance, species diversity, and richness of the confluence point in comparison to the two streams

Mariam Dahan stream is the confluence of Elala and GTS streams. This site has two sampling sites Mariam Dahan and Romant. The fish species in the first site (Mariam Dahan) were similar to the two streams. The cyprinids of the genus Garra namely; G. blanfordii, G. dembecha and G. ignestii (Fig. 2) were recorded from this site. The total fish abundance of the three fish species was 566,145, and 122, respectively. Similar to most sites of the two streams the pools of this site retain the highest total fish abundances and riffles were the lowest. The fish species richness of the three microhabitats was equal; three species each (Table 5). Whereas the second site (Romant stream) retains a maximum of eight species where four of the fish species were not recorded from the 13 upper sites (headwater streams). These fish species were all cyprinids namely, V. beso, L. intermedius, L. niloticus, and L. forskalii (Table 3; Fig. 2).

Discussion

From the above observations, it is obvious that all the physicochemical parameters and habitat characteristics like water depth, velocity, and river bed width are some of the major factors for the distribution, abundance, and richness of species in different habitats. Habitat wise and/or between streams parametric t-test and other statistical comparisons showed the variability of the above metrics. Gorman & Karr (1978), Arunachalam (2000), Johal et al. (2002) and Negi et al. (2007) made similar observations.

To the best of our knowledge, this study represents the first study of riverine fish diversity, richness, and abundance in the streams: GTS and Elala. We investigate the patterns of fish population composition along longitudinal (habitat) gradients. We showed that the two streams and their confluence point (Mariam Dahan) are inhabited generally by small-sized noncommercial fish Garra species and there was no longitudinal variation in species composition. However, at the second site of Mariam Dahan in Romanat particularly below the long waterfall (Romanat Fall) the species richness doubled. Among these species, three of them were commercially important.

Understanding the diversity, richness, and abundance of fish in stream and river ecosystems is among the central goals of tropical ecological research (Herder & Freyhof, 2006). Dominance by the ecologically diverse Cyprinid family is common in East African and Southwest Asian freshwater systems (Beamish et al., 2006; Stiassny & Getahun, 2007) which are also evidenced in the results of the current study. This is because cyprinids have evolved partially through highly adapted body forms and mouth structures to occupy all habitats (particularly lakes) throughout their distributions (Ward-Campbell et al., 2005).

Differences in stream types and habitat types each explained a significant proportion of the variation in the numerical abundance of fishes. Fish abundance strongly differed between streams and among habitat types. GTS harbors the highest fish abundance which might be due to higher food availability from organic matter and waste material inputs from the town. Yet, large debris of riparian vegetation and waste disposal from the town was observed at various points via the river course of GTS. Similar studies elsewhere showed a positive effect of the input of plant debris into streams on fish abundance (Wright & Flecker, 2004).

Our findings also showed that, among the habitats, the pool has the largest fish abundance, which is consistent with previous research (Freeman & Marcinek, 2006). The large and significant amount of variation explained by habitat type may indicate that habitat use is indeed potentially an important factor in determining total fish abundance (Tesfay et al., 2019).

The current study’s environmental variables (such as temperature, pH, and dissolved oxygen) had rather modest or weak associations with fish abundance. For the prediction of fish diversity and abundance, it was found that stream type and habitat type were the best predictors for fish abundance and species richness, compared to physicochemical variables (Tesfay et al., 2019).

Moreover, the two studied streams GTS and Elala including their junction area Mariam Dahan are generally intermittent and very shallow streams. But at some intervals they consist of some deep pools, important habitats for the Garra fish serving as a refuge and feeding habitat during the severe dry season due to the formation of patches. This mainly happened as result of pumping of the stream water for irrigation and construction issues such as bricks production.

Physical habitat structure is widely recognized as a major determinant of the distribution, abundance, and diversity of stream fishes (Pusey et al., 1993; Rice, 2005). The study streams are quite small and sporadic, with relatively shallow runs and riffles that are vulnerable to assaults from predators. Because of this, the research streams and pools had greater species diversity, abundance, and richness. Additionally, among the habitat types it was demonstrated that the pools are the lotic habitats which also support the highest mean species richness followed by runs and riffles.

Pool and run habitats were architecturally more complex and always deeper than riffles, as shown from the research streams. As a result, they offer fish a better habitat than riffles. Riffles often had a simple structural design, were shallow, and had little fish diversity or richness. Thompson & Larsen (2004) made a comparable observation and evaluation. Additionally, pools sustain the greatest diversity of fish because they have higher depths and slower currents than riffles and runs, as noted by Negi & Negi (2010) in their explanation of the Shannon-Weaver diversity index.

In this study, among the Garra fishes G. blanfordii dominates the abundance data, as they occupy all possible habitats. This might be due to their high adaptive variability, similar results were revealed by Johnson & Arunachalam (2009) in the study of cyprinids. Jayaram (1999) and others have also reported similar results on the cyprinids like Garra mullya. The wide distribution and large number of these species suggest that they can tolerate a wide range of environmental conditions (Pusey et al., 1993).

After the confluence point particularly in the second site of Mariam Dahan below the Romanat Fall, the fish species richness increased. The Garra species recorded in GTS and Elala streams were also found in Romant stream including G. aethiopica which was limited to GTS stream. Similar to GTS stream it was found in the microhabitats such as pools and runs. Although fish assemblage structure is usually higher in lower regions (lower altitudes) of streams and rivers in general as reported by different authors (e.g., Habteselassie, 2012; Paller, 1994; Rahel & Hubert, 1991), the barriers created by the natural waterfall might also affect the variability in fish species assemblage, distribution, and richness between the up and down steam of the waterfall (Liu et al., 2021). Literatures indicate that waterfalls are an important feature of freshwater ecosystems, affecting not only the physical environment but also the biological communities inhabiting them (Evelyn Hutchinson, 1957; Torrente-Vilara et al., 2011). Studies have shown that waterfalls can significantly impact the fish assemblage structure of rivers and streams, with different species adapting to the varying conditions created by waterfalls (Barbosa et al., 2015; Kano et al., 2012; Liu et al., 2021; Torrente‐Vilara et al., 2011). The effects of waterway barriers such as natural waterfalls on fish movements are expected to produce differing assemblage structures in riverine ecosystems (Rahel, 2007). In addition, such natural barriers may cause faunal discontinuities and increase dissimilarities in the ichthyofauna, as changes in landscape characteristics cause habitat alterations and increase species turnover along the longitudinal gradient (Rahel & Hubert, 1991).

Conclusion

In this study, eight fish species of family Cyprinidae representing four genera were identified. We observed that the fish abundance differed between streams and among habitat types. GTS harbors the highest fish abundance and among the habitats, the pool has the largest fish abundance. In general, the fish abundance in the study streams of GTS and Elala exhibited clear patterns among the gradients in stream type and habitat type. Most likely, habitat type and food resource differences between the streams, possibly combined with other unmeasured environmental gradients, were the most important driving factors behind variation among fish abundances. We do not interpret the observed patterns as signatures of dispersal limitation given that almost all species were widely distributed across the study area, and because habitat and stream type were found to be strongly collinear with each other. The wide spread of fish species in most of the study sites indicates a relatively high degree of mobility of species across the stream courses via the network of stream connections. Hence the effect of dispersal limitation may not be the case except for the large and commercially important fish’s limited upward dispersal by the waterfall.

Competing interests

No potential conflict of interest relevant to this article was reported.

Funding sources

This work was funded by the College of Natural and Computational Sciences (CNCS), Mekelle University, from its recurrent budgets.

Acknowledgements

The authors would like to acknowledge the Department of Biology, College of Natural and Computational Sciences, Mekelle University, for laboratory facilities and financial support. We are grateful to Belay Gebreyohannes and to all Aquatic research team members for their help during field sample collection and technical support.

Availability of data and materials

Upon reasonable request, the datasets of this study can be available from the corresponding author.

Ethics approval and consent to participate

It is to declare that we have all the ethical approval and consent to take participate in research. The study compiles all regulations for ethical clearance. This study does not involve the use of human subjects (human-related experimentation) but involves field sampling of fish from non-protected streams. Permission to carry out field observation was obtained from Mekelle University, College of Natural and Computational Sciences, along with project registration CNCS/MSc/RB/12/2013.

References

1.

Arunachalam M. Assemblage structure of stream fishes in the Western Ghats (India). Hydrobiologia. 2000; 430:1-31

2.

Asmelash T, Dejenie T, Decleerck S, Nyssen J, van der Guch K, Risch S, et al. Ecological atlas of reservoirs in Tigray, Northern Ethiopia. Tigray livelihood papers No. 4 VLIR – Mekelle University IUC Programme. 2007.

3.

Barbosa TAP, Benone NL, Begot TOR, Gonçalves A, Sousa L, Giarrizzo T, et al. Effect of waterfalls and the flood pulse on the structure of fish assemblages of the middle Xingu river in the eastern Amazon basin. Braz J Biol. 2015; 75:78-94

4.

Beamish FWH, Sa-ardrit P, Tongnunui S. Habitat characteristics of the cyprinidae in small rivers in Central Thailand. Environ Biol Fish. 2006; 76:237-53

5.

Cheng L, Lek S, Lek-Ang S, Li Z. Predicting fish assemblages and diversity in shallow lakes in the Yangtze River basin. Limnologica. 2012; 42:127-36

6.

Cushing CE, David Allan J. Streams: their ecology and life. Cambridge, MA: Academic Press. 2001

7.

Dejenie T, Asmelash T, De Meester L, Mulugeta A, Gebrekidan A, Risch S, et al. Limnological and ecological characteristics of tropical highland reservoirs in Tigray, Northern Ethiopia. Hydrobiologia. 2008; 610:193-209

8.

Edds DR. Fish assemblage structure and environmental correlates in Nepal’s Gandaki river. Copeia. 1993; 1993:48-60

9.

Evelyn Hutchinson G. Concluding remarks. Cold Spring Harbor symposia on quantitative biology. Cold Spring Harb Symp Quant. 1957; 22:415-27

10.

Freeman MC, Marcinek PA. Fish assemblage responses to water withdrawals and water supply reservoirs in piedmont streams. Environ Manag. 2006; 38:435-50

11.

Godinho FN, Ferreira MT, Santos JM. Variation in fish community composition along an Iberian river basin from low to high discharge: relative contributions of environmental and temporal variables. Ecol Freshw Fish. 2000; 9:22-9

12.

Gordon ND, McMahon TA, Finlayson BL. Stream hydrology: an introduction for ecologists. Chichester: John Wiley & Sons. 1992.

13.

Gorman OT, Karr JR. Habitat structure and stream fish community. Ecology. 1978; 59:507-15

14.

Grossman GD, Freeman MC. Microhabitat use in a stream fish assemblage. J Zool. 1987; 212:151-76

15.

Habteselassie R. Fishes of Ethiopia: annotated checklist with pictorial identification guide. Addis Ababa: EFASA. 2012; p p. 250.

16.

Herder F, Freyhof J. Resource partitioning in a tropical stream fish assemblage. J Fish Biol. 2006; 69:571-89

17.

Ibarra AA, Gevrey M, Park YS, Lim P, Lek S. Modeling the factors that influence fish guilds composition using a back-propagation network: assessment of metrics for indices of biotic integrity. Ecol Model. 2003; 160:281-90

18.

Jayaram KC. The freshwater fishes of the Indian region. New Delhi: Narendra Publishing House. 1999.

19.

Jayaratne R, Surasinghe T. General ecology and habitat selectivity of fresh water fishes of the Rawan Oya, Kandy, Sri Lanka. Sabaramuwa Univ J. 2010; 9:11-43

20.

Jha BR, Waidbacher H, Sharma S, Straif M. Length-weight relationship of sucker head, Garra Gotyla Gotyla (Gray, 1830) in different rivers of Nepal and the influence of monsoon. Int J Environ Sci Technol. 2005; 2:147-53

21.

Johal M, Krishan-Tandon K, Kumar-Tyor A, Kumar-Rawal Y. Fish diversity in different habitats in the streams of lower Middle Western Himalayas. Pol J Ecol. 2002; 50:45-56.

22.

Johnson JA, Arunachalam M. Diversity, distribution and assemblage structure of fishes in streams of southern Western Ghats, India. J Threat Taxa. 2009; 1:507-13

23.

Jones NE, Tonn WM, Scrimgeour GJ, Katopodis C. Productive capacity of an artificial stream in the Canadian Arctic: assessing the effectiveness of fish habitat compensation. Can J Fish Aquat Sci. 2003; 60:849-63

24.

Kano Y, Nishida S, Nakajima J. Waterfalls drive parallel evolution in a freshwater goby. Ecol Evol. 2012; 2:1805-17

25.

Kimmel WG, Argent DG. Efficacy of two-pass electrofishing employing multiple units to assess stream fish species richness. Fish Res. 2006; 82:14-8

26.

King M. Fisheries biology, assessment and management. Oxford: Fishing News Books. 1995.

27.

Koel TM, Peterka JJ. Stream fish communities and environmental correlates in the Red River of North, Minnesota and North Dakota. Environ Biol Fish. 2003; 67:137-55

28.

Kumar R. Biology of fishes. Form and function, modes of reproduction. FISH/BIOL 311 [Internet] 2013.[cited 2013 Oct 22]https://files.polsri.ac.id/Magazines/All%20Animals/Fish//Quentin_Bone,_Richard_Moore-Biology_of_Fishes_(English_and_French_Edition)-Taylor_&_Francis(2008).pdf.

29.

Kwak TJ, Peterson JT. 15: Community indices, parameters, and comparisons.In: In: Guy CS, Brown ML, editors.editors Analysis and interpretation of freshwater fisheries data. Bethesda, MDThe American Fisheries Society 2007; p p. 677-763

30.

Liu F, Li M, Wang J, Gong Z, Liu M, Liu H, et al. Species composition and longitudinal patterns of fish assemblages in the middle and lower Yarlung Zangbo River, Tibetan Plateau, China. Ecol Indic. 2021; 125:107542

31.

Magurran AE. Ecological diversity and its measurement. Cambridge: Cambridge University Press. 1988

32.

Morrow JV, Fischenich C. Habitat requirements for freshwater fishes. Vicksburg, MS: EMRRP. 2000; p p. 14.

33.

Negi RK, Joshi BD, Negi T, Chand P. A study on stream morphology of some selected streams hill streams of district Nainital with special reference to its biotic communities. In: Proceedings of National Seminar on Limnology. 2007 Jaipur, India. .

34.

Negi RK, Negi T. Assemblage structure of stream fishes in the Kumaon Himalaya of Uttarakhand state, India. Life Sci J. 2010; 7:9-14.

35.

Paller MH. Relationships between fish assemblage structure and stream order in South Carolina coastal plain streams. Transact Am Fish Soc. 1994; 123:150-61

36.

Pusey BJ, Arthington AH, Read MG. Spatial and temporal variation in fish assemblage structure in the Mary River, south-eastern Queensland: the influence of habitat structure. Environ Biol Fish. 1993; 37:355-80

37.

Rahel FJ. Biogeographic barriers, connectivity and homogenization of freshwater faunas: it’s a small world after all. Freshw Biol. 2007; 52:696-710

38.

Rahel FJ, Hubert WA. Fish assemblages and habitat gradients in a rocky Mountain–Great plains stream: biotic zonation and additive patterns of community change. Trans Am Fish Soc. 1991; 120:319-32

39.

Rashleigh B. Relation of environmental characteristics to fish assemblages in the upper French broad river basin, North Carolina. Environ Monit Assess. 2004; 93:139-56

40.

Rice JC. Understanding fish habitat ecology to achieve conservation. J Fish Biol. 2005; 67:1-22

41.

Scottish Fisheries Coordination Centre [SFCC]. Eletrofishing team leader training manual. Fisheries management SVQ level 3: managing electrofishing operations. Inverness: SFCC. 2007.

42.

Shannon CW, Weaver W. The mathematical theory of communication. Champaign, IL: University of Illinois Press. 1949.

43.

Sovell L, Heiskary S, Anderson J. Using the transparency tube in Minnesota’s new Citizen Stream Monitoring Program. In: Proceedings of the 1999 National Water Quality Monitoring Conference. 2000 Austin, TX, USA. .

44.

Stiassny MLJ, Getahun A. An overview of labeonin relationships and the phylogenetic placement of the Afro-Asian genus Garra Hamilton, 1922 (Teleostei: Cyprinidae), with the description of five new species of Garra from Ethiopia, and a key to all African species. Zool J Linn Soc. 2007; 150:41-83

45.

Teferi M, Declerck SAJ, De Bie T, Gebrekidan A, Asmelash T, Dejenie T, et al. The ecology of the riverine Garra species (Teleostei, Cypriniformes) in reservoirs of the semi-arid highlands of northern Ethiopia: temporal dynamics of feeding activity. Inland Waters. 2013; 3:331-40

46.

Teferi M, Declerck SAJ, De Bie T, Lemmens P, Gebrekidan A, Asmelash T, et al. Strong effects of occasional drying on subsequent water clarity and cyanobacterial blooms in cool tropical reservoirs. Freshw biol. 2014; 59:870-84

47.

Tesfay S, Teferi M, Tsegazeabe HH. Habitat selectivity of fresh water fishes of two second-order tropical streams in Tigray, Northern Ethiopia. J Ecol Environ. 2019; 43:1-11

48.

Thompson L, Larsen R. Fish habitat in freshwater stream. Irvine, CA: UCANR Publications. 2004

49.

Torrente-Vilara G, Zuanon J, Leprieur F, Oberdorff T, Tedesco PA. Effects of natural rapids and waterfalls on fish assemblage structure in the Madeira river (Amazon Basin). Ecol Freshw Fish. 2011; 20:588-97

50.

Van de Meutter F. Local and regional processes in macroinvertebrate communities in shallow lakes. Belgium: Katholieke Universiteit Leuven, Faculteit Weteschappen Departement Biologie, Laboratorium voor Aquatische Ecologie. 2005.

51.

Van de Meutter F, de Meester L, Stoks R. Metacommunity structure of pond macroinvertebrates: effects of dispersal mode and generation time. Ecology. 2007; 88:1687-95

52.

Vila-Gispert A, Moreno-Amich R, Garcia-Berthou E. Gradients of life-history variation: an intercontinental comparison of fishes. Rev Fish Biol Fish. 2002; 12:417-27

53.

Ward-Campbell BMS, Beamish FWH, Kongchaiya C. Morphological characteristics in relation to diet in five coexisting Thai fish species. J Fish Biol. 2005; 67:1266-79

54.

Wright JP, Flecker AS. Deforesting the riverscape: the effects of wood on fish diversity in a Venezuelan piedmont stream. Biol Conserv. 2004; 120:439-47

55.

Zippin C. The removal method of population estimation. J Wildl Manag. 1958; 22:82-90