Study Area and Sampling
Sıddıklı Küçükboğaz Dam Lake (formerly called Karababa Dam) is a zonal dam on Körpeli Boğaz Creek at the border of the Province of Kırşehir located from Central Anatolian region of Türkiye. This dam lake is composed of clay and rock, with a central core. Dam construction began in 1991 and was finished in 2002. The lake dam is made of clay and filled with rock. It has a 1.62 km2 surface area with an active water level of 25.3 hm3 (Akkan et al. 2018). Pike gill nets were simultaneously used to collect fish samples monthly from September 2015 to August 2016 at Sıddıklı Küçükboğaz Dam Lake. The nets were set at dusk, left in the water for 12-hour min, and hauled at 08:00–09:00 h. The nets were composed of pelagic gill nets with bar mesh sizes (knot to knot) of 20, 25, 30, 35, 40, 45, 55, 65, and 80 mm. All operations on fish capture and dead fish studied in the laboratory were carried out in accordance with animal health and welfare ethical rules. This study was approved by the animal experiments local ethics committee (document no: 68429034/05). Additionally, our study complies with ARRIVE 2.0 guidelines. Some physical and chemical parameters of the lake water, such as dissolved oxygen (DO), temperature, salinity, pH, conductivity, total dissolved solid (TDS) were measured monthly from September 2015 to August 2016 (Yazıcı 2018) and environmental variables were given Table 1.
Table 1
Monthly changes in some physicochemical parameters of surface water in Sıddıklı Dam Lake from September 2015 to August 2016
Months | Water Parameters |
Temperature | pH | DO | Salinity | TDS | Conductivity |
September 2015 | 22.91 | 8.22 | 9.38 | 0.48 | 6.63 | 940.00 |
October | 15.77 | 7.51 | 14.90 | 0.52 | 6.73 | 850.00 |
November | 10.40 | 8.41 | 7.16 | 0.78 | 9.95 | 1100.00 |
December | 3.68 | 8.41 | 11.03 | 0.63 | 8.29 | 750.00 |
January 2016 | 3.10 | 8.21 | 11.79 | 0.40 | 5.35 | 820.00 |
February | 7.40 | 8.15 | 11.76 | 0.44 | 5.77 | 890.00 |
March | 8.78 | 8.25 | 9.24 | 0.35 | 4.55 | 700.00 |
April | 14.93 | 8.16 | 8.79 | 0.34 | 4.48 | 690.00 |
May | 18.10 | 8.31 | 7.18 | 0.34 | 4.52 | 700.00 |
June | 20.85 | 8.40 | 6.86 | 0.32 | 4.28 | 650.00 |
July | 24.48 | 8.31 | 6.02 | 0.33 | 4.42 | 680.00 |
August | 23.05 | 8.37 | 5.97 | 0.33 | 4.72 | 700.00 |
Laboratory Methods and Stomach Content Analysis
A total of one hundred thirty-three (133) samples were examined for stomach analysis. In the laboratory, all fish samples were measured in cm (total length) and weighed in grams. The stomachs were removed by dissection from each specimen and preserved in a 4% formaldehyde solution for afterward analysis. The stomachs were opened during the examination, and the prey was identified, weighed in grams, sorted, and classified to the lowest taxonomic level before being preserved in 70% ethanol. Also, full and empty stomach weights were measured with a precision of 0.01 g. When a prey item was mostly digested, identification of prey fishes was based on scales, pharyngeal bones (cyprinids), opercular bones, vertebrae, and the location of the mouth and eyes (Pavlović et al. 2015).
To compare the change of feeding intensity between seasons the fullness index (FI = weight of stomach content/ weight of fish *100) and the vacuity index (VI%= the number of empty stomachs/total number of the examined stomachs* 100) were calculated (De Santis and Volta 2021). Low feeding activity is considered when high vacuity index is observed (Martinho et al., 2012). Kruskal-Wallis test (K-W test) was used to analyze whether seasons affected the fullness index (FI) in northern pike. Spearman’s rank correlation was used to determine the relationship between fullness index (FI) and physicochemical parameters of surface water. Also, A chi-square test (χ2) was applied to determine the vacuity index (%VI) changes between the seasons.
Traditional methods such as percentage frequency of occurrence (FO%= number of stomachs containing prey i/number of stomachs with any food item* 100), numerical percentage (N%= number of prey i/total number of all prey items* 100), and percentage by weight (W%= weight of prey i/total weight of all prey items* 100) of dietary analysis were used to determine feeding features (Hyslop 1980). The main food items were identified using index of relative importance (IRI) of Pinkas et al. (1971), as modified by Hacunda (1981).
$$\:IRI=\left(\%N+\%W\right)\times\:\%FO$$
This index has been expressed as the percentage of each prey item;
$$\:\%IRI=\left(IRI/\sum\:IRI\right)x\:100$$
For computation of the relative amounts of intraspecific competition between seasons, simplified Morisita-Horn index (Ch) based on %N data was used (Horn 1966):
$$\:{C}_{h}=\frac{2\:\left(\sum\:{p}_{ij}{p}_{ik}\right)}{\sum\:{p}_{ij}^{2}+\sum\:{p}_{ik}^{2}}$$
where Ch is the Morisita-Horn index of diet overlap between different seasons pij is the proportion of food type “i” f the total food quantity by seasons “j,” pik is the proportion of food type “i” of the total food used by seasons “k” and, n is the total number of food types. The degree of overlap was classified as low (0.0–0.29), moderate (0.30–0.59), and high (0.60–1.00) (Langton 1982).
The selectivity of prey categories in the diet was statistically tested with x2-test, utilizing Pearre's C index of prey selection. The index value (Va) varies from − 1 (prey avoidance) to + 1 (prey selection), with 0 indicating random prey selection (Pearre 1982).
$$\:Va=\frac{({a}_{d}\times\:{b}_{e})-({a}_{e}\times\:{b}_{d})}{\sqrt{a\times\:b\times\:d\times\:e}}$$
where Va is Pearre’s index for pike selection of prey a, ad is the abundance of prey a in the diet, be is the abundance of all other prey in the environment, bd is the abundance of all other prey in the diet, and ae is the abundance of prey a in the environment.
Values without subscripts are expressed as follows:
The statistical significance of the selection index value (Va) was tested using the chi-squared test.
$$\:{x}^{2}=n\times\:{C}^{2}$$
The value of relative abundance (ae) used in the prey selection index for each fish species inhabiting Sıddıklı Dam Lake was obtained from Yazıcı (2018).
Feeding strategy was determined from the plot of percentages of prey-specific abundance (Pi%) against frequency of occurrence (FO). Prey specific abundance, the percent numerical abundance of a prey item averaged over the stomach samples in which it occurs, was calculated using the methodology in Amundsen et al. (1996).
$$\:Pi=\left(\frac{\sum\:Si}{\sum\:Sti}\right)*100$$
where Pi = Prey specific abundance of prey i; Si = Abundance of prey in stomachs and Sti = Total abundance of prey in predators that contain prey i. For specialist feeding, prey items appear in the upper part of the plot, while generalists have all prey points in the lower part.