118 Bulgarian Journal of Agricultural Science, 19 (Supplement 1) 2013, 118 125 Agricultural Academy AGE COMPOSITION AND GROWTH RATE OF THE SPAWNING PART OF THE POPULATION OF PONTIC SHAD ALOSA IMMACULATA (BENNETT, 1835) IN THE BULGARIAN SECTOR OF DANUBE RIVER D. ROZDINA, G. RAIKOVA-PETROVA and P. MIRTCHEVA Sofi a University, Faculty of Biology, Department of Common and Applied Hydrobiology, BG 1164 Sofi a, Bulgaria Abstract ROZDINA, D., G. RAIKOVA-PETROVA and P. MIRTCHEVA, 2013. Age composition and growth rate of the spawning part of the population of pontic shad Alosa immaculata (Bennett, 1835) in the Bulgarian sector of Danube River. Bulg. J. Agric. Sci., Supplement 1: The age structure and the growth ratе of the population of Pontic shad Alosa immaculata (Bennet, 1835) in the Bulgarian sector of the Danube River are studied. The catchments have been done in April 2010 and April 2012 during the spawning migration of the species. Altogether 239 specimens have been examined. The age of the individuals varied between II and V. The most abundant are II and III age groups representing 53.2% and 37.9% of the spawning part of the population respectively. The standard length ranged from 15.4 cm to 37 cm and the weight ranged from 37 g to 351 g. The equation representing the relation between the standard length and the scale radius in the population is L = 0.3955S 0.996; r = 0.8. The length weight relationship is described by the equation W = 0.0234L 2.7315, r = 0.95. The von Bertalanffy`s length and weight growth equations are L t = 35.749 [1 e -0.493209 (t-0.34115) ] and W t = 436.228 [1 e -0.319987 (t 1.04376) ] 2.7315 respectively. Key words: Pontic shad, Alosa immaculata, Danube River, age, growth Abbreviations: NAFA National Agency of Fisheries and aquaculture Introduction The Pontic shad (Alosa immaculata) is anadromous fish species, belonging to the family Clupeidae. The species is native for Bulgaria, Georgia, Moldova, Romania, Russia, Serbia, Turkey and Ukraine. It occurs in Black sea and Sea of Azov and for spawning, migrates in Danube, Dnepr, Dniester, Don, Bug etc. Since 2001, the species is established in the Sea of Marmara (Eryilmaz, 2001). The species has been regularly studied on the territory of Bulgaria (Kolarov, 1958a; 1958b; 1960a, 1960b; 1961; 1963; 1964; 1965; 1978; 1980; 1982; 1983; 1985; 1989; Ivanov and Kolarov, 1979 and Prodanov and Kolarov, 1983) but in the recent years the studies are rare (Yankova et al., 2011). The structure of the Pontic shad in Danube River has been also studied by Ciolac and Patriche (2004). The species is vulnerable according to IUCN and Bulgarian Red Data Book (http://e-ecodb.bas.bg/rdb/bg/ vol2/alpontic.html; http://www.iucnredlist.org/apps/redlist/ details/907/0). It is also included in Annex 2 and 4 of the Bulgarian Biodiversity Act (http://www.biodiversity. bg/files/file/zak_bg_biodiv.pdf). A. immaculata is a commercial fish species. Its catchments in the last decade vary significantly with a maximum in 2004 (21.95 t) and minimum in 2007 (0.2 t) (according to data from NAFA). This catchment vari- E-mail: d_rozdina@abv.bg, galerida@abv.bg, petia.mircheva@abv.bg
Age Composition and Growth Rate of the Spawning Part of the Population of Pontic Shad... 119 ability focused our attention on studying some basic biological parameters of the population. The aim of the article is to assess the age and size structure, growth rate and condition factor of the reproductive part of the population of A. immaculata in Bulgarian sector of Danube River. Material and Methods The specimens of A. immaculata were collected in April 2010 and April 2011 during the spawning migration of the species. Spawning sites were located in the area between Lom and Vidin. The geographical coordinates of the sampling points are as follows: Sampling point 1: 43 54 17.96 N and 22 50 27.45 E; Sampling point 2: 43 47 40.46 N and 23 4 50.41 E; Sampling point 3: 43 50 45.61 N and 23 18 44.17 E. Altogether 239 specimens were collected (134 individuals in 2010 and 159 in 2011 respectively). Gill nets were used with the size of the eye 32 88 mm. Each specimen was measured the standard length (L) to the nearest 1 mm, the total weight (W) and the gutted weight (w) to the nearest 1g. The age was determined by the scales at magnification of 17.5x with Projector Dokumator, Lasergeret (Carl Zeiss, Jena). Age and size structure was studied according to (Chugunova, 1959). Length and weight growth rates were determined by back calculation using the following equations: L = a + bs (L fish length, in cm and S scale radius, in units) and W = a.l b (W fish weight, in g) (Le Cren, 1951). Back-calculated lengths and weights at age were used to calculate the von Bertalanffy s growth parameters, using the software Growth II. Linear and weight growths were described using the von Bertalanffy (1938) growth equations: L t = L (1-e -k(t-to) ), where L t (in cm) is the length of fish at age t, L (in cm) is the asymptotic length, k (in yr -1 ) is the growth coefficient and t o (in yr) represents the age a fish would have at zero length. W t =W [1-e -k(t-to) ] b, where W t (in g) is the weight of fish at age t, W (in g) is the asymptotic weight, k (in yr -1 ) is the growth coefficient and t o (in yr) represents the age a fish would have at zero length. To asses if the growth potential of the population was well used the coefficient of Hohendorf was applied (Hohendorf, 1966). Condition of the population has been studied in three ways: By the equation of Fulton: k f = (W/L 3 ).100, where W is the fish weight (in g) and L is the fish length (in cm) By the equation of Fulton but instead the exponent 3, the exponent b from the length-weight relationship of the population was used (k b ). Calculating fish weight at the length of 5 cm, 10 cm, 15 cm, 20 cm and 25 cm with the use of lengthweight relationship. The population with higher weight at the same length of fish has better condition (Goldspink, 1979; De Silva, 1985; Basami and Grove, 1985; Raikova-Petrova, 1992 and Zivkov, 1999). Results and Discussion Age structure of Pontic shad in Danube River was represented from 4 age groups (from II to V) (Table 1). The average age of the population was 2.6 years. The most abundant were two and three year old individuals. The oldest specimens were two 5 years old males, which represent the less abundant age group. The absence of individuals older than 5 years was due to the commercial fishing of the species in the Danube Delta where selectively the biggest and the oldest individuals were caught. The established age structure shows, that in the spawning part of the population of A. immaculata replenishment dominates over the residue. Kolarov (1965) reported similar results. Kolarov (1960a) reported similar age structure for the Black sea catchments along the Bulgarian shore. Ciolac and Patriche (2004) report six age groups (from 2 to 7 years old individuals) for Pontic shad in Danube River. In a study the most abundant were three (41.2%) and four (30.3%) years old individuals. Kolarov (1965; 1980) reported for the Bulgarian sector of Danube River age structure of 1 to 6 years old individuals of Pontic shad.
120 D. Rozdina, G. Raikova-Petrova and P. Mirtcheva Table 1 Age-size composition of Alosa immaculata catchments in Bulgarian sector of Danube River Size class Age 2 3 4 5 % n 15 15.9 3 1 1.4 4 16 16.9 12 1 4.4 13 17 17.9 34 12 15.7 46 18 18.9 20 20 1 14.0 41 19 19.9 15 8 7.8 23 20 20.9 19 4 7.8 23 21 21.9 7 6 1 4.8 14 22 22.9 4 5 1 3.4 10 23 23.9 4 1.4 4 24 24.9 2 2 2 2.0 6 25 25.9 2 5 3 1 3.8 11 26 26.9 5 18 8 10.6 31 27 27.9 11 9 5 1 8.9 26 28 28.9 10 15 2 9.2 27 29 29.9 6 2 1 3.1 9 30 30.9 2 1 1.0 3 31 31.9 1 0.3 1 37 37.9 1 0.3 1 Total number 156 111 24 2 100 293 % 53.2 37.9 8.2 0.7 Size structure of A. immaculata was represented with 18 size classes (Table. 1). The most abundant were the individuals with size range 17 17.9; 18 18.9 and 26 26.9. The biggest individual was male with length 37 cm, weight 225 g and was 3 years old. The specimen with the highest weight (W = 351 g) was two years old male with a length of 29 cm. The individual with the smallest size had a length of 15.4 cm, weight 37 g, and was 3 years old. The size classes with the highest age range were 25 25.9 and 27 27.9. The III age group had the highest size range and included 17 size classes. The V age group included only two size classes. Kolarov (1978) has established 19 size classes in the area of Silistra and Svishtov. The most abundant were 21 22 and 23 34 size classes for the both area respectively. Other authors describe smaller number of size classes (Pavlov, 1953; Kolarov, 1964; 1980). During our study, we established a size structure with the smallest average size of the individuals in comparison with previous studies. The main reason is the high commercial catchments rates in Danube Delta where selectively the biggest individuals are caught. Temp of growth The equation representing the relation between the fish length (L) and the scale radius (S) for the population of Pontic shad is L = 0.3955S 0.996; r = 0.8. The relation is represented on Figure 1. On Table 2 are represented the back calculated lengths at age and the annual increments. 40 35 30 25 20 15 10 5 0 0 10 20 30 40 50 60 70 80 90 S Fig. 1. Relation between the fish length (L, cm) and scale radius (S, units of ocular-micrometer) in Pontic shad (Alosa immaculata) in the Bulgarian sector of Danube River L, cm W, g 600 500 400 300 200 100 0 W = 0,0234L 2,7315 r = 0.95 0 5 10 15 20 25 30 35 40 L, cm Fig. 2. Relation between the length (L, cm) and weight (W, g) of Pontic shad (Alosa immaculata) in the Bulgarian sector of Danube River
Age Composition and Growth Rate of the Spawning Part of the Population of Pontic Shad... 121 Table 2 Back calculated length (L, cm) at age and annual length increments for Pontic shad (Alosa immaculata) in the Bulgarian sector of Danube River Generation (year) Age group Back calculated average lengths (L, cm) at the end of each vegetation period Number of fish L1' L2' L3' L4' L5' 2009 II 11.58 17.94 100 2008 II III 11.62 19.43 32.45 111 2007 III IV 10.67 19.45 26.21 43.10 60 2006 IV 10.51 19.41 25.23 30.98 20 2005 V 8.10 15.22 20.76 24.91 30.05 2 Average length 10 496 18.29 26.20 32 997 30.05 293 Increments 10 496 7.8 7.9 6.8 2.9 Table 3 Back calculated weight (L, cm) at age and annual weight increments for Pontic shad (Alosa immaculata) in the Bulgarian sector of Danube River Generation (year) Age group Back calculated average weights (W. g) at the end of each vegetation period Number of fish W1' W2' W3' W4' W5' 2009 II 18.83 62.24 100 2008 II III 19.01 77.41 314.06 111 2007 III IV 15.06 77.62 175.37 682.11 60 2006 IV 14.44 77.19 157.89 276.77 20 2005 V 7.09 39.71 92.69 152.54 254.67 2 Average weight 14.89 66.83 185.00 370.47 254.67 293 Increments 14.89 51.9 118.2 185.5 115.8 The length-weight relationship for the population of A. immaculata in Bulgarian sector of Danube River is W = 0.0234L 2.7315, r = 0.95 (Figure 2). The back calculated weight at age is represented on Table 3. The average length and weight increases with increasing the fish age. Similar results were obtained for other populations of A. immaculata (Kolarov, 1960а; 1964; 1965; 1978; 1980; 1983; 1985; Ciolac and Patriche, 2004 and Ergüden et al., 2007). Generation 2005 showed the slowest temp of growth while the fastest growth rate was found for generations 2007 and 2008. The nature of linear and weight growth more precisely is described by the change of the annual increments. The length increments decrease with increasing fish age, while the weight increments increase with increasing fish weight until the fifth year. The yearlings show highest values of length increments. On the fifth year negative length and weight increments were established (Tables 2 and 3). The negative increments are due to the small number of 5 year old fish represented in the excerption as well as due to the commercial fishing in Danube s Delta when the biggest individuals are caught. Kolarov (1978) has established decreasing the average weight of the Pontic shad from Danube s mouth to Silistra and Svishtov. In comparison to other populations (Pavlov, 1953; Kolarov, 1964; 1965; 1978; 1980; 1983; 1985; Ciolac and Patriche, 2004; Ergüden et al., 2007) of the species the studied by us showed the lower linear increments during the first and fifth year of life, while two, three and four year old fish showed higher increments. The higher increments for the two and three year old indi-
122 D. Rozdina, G. Raikova-Petrova and P. Mirtcheva Observed Predicted Fig. 4. von Bertalanfy s growth equation curve describing the linear growth of Pontic shad (Alosa immaculata) from Bulgarian sector of Danube River Observed Predicted Fig. 3. Ford-Walford curve for the linear growth of Pontic shad (Alosa immaculata) from Bulgarian sector of Danube River viduals are due to the faster growth rates of the younger age groups and the earlier maturation of the population. For some of the populations the changes in the annual increments are uneven (Pavlov, 1953; 1965; 1978; 1980; 1985; Ciolac and Patriche, 2004 and Ergüden et al., 2007). For other populations there is a clear tendency of decreasing the increment with the fish age (Kolarov, 1964; 1978; 1983, 1985 based on data by Miklashevskaya, 1953). Observed Predicted Fig. 5. Ford-Walford curve for the weight growth of Pontic shad (Alosa immaculata) from Bulgarian sector of Danube River Concerning weight increments in comparison to other Pontic shad populations the studied by us showed highest values during the third and fourth year. As for the linear increments changes in the weight increments are uneven for some of the populations (Kolarov, 1964;
Age Composition and Growth Rate of the Spawning Part of the Population of Pontic Shad... 123 Table 4 Von Bertalanffy s growth parameters for the linear and weight growth of A. immaculata Linear growth Water body and autor L, cm k t o Danube River (our data) 35 749 0.493209 0.34115 Danube River (Kolarov, 1980) 57.38 0.1067 1727 Danube River, Black sea (Kolarov, 1983) 40.43 0.2705-0.218 Black sea (Prodanov, Kolarov, 1983) 40.43 0.27-0.218 Weight growth W k t o Danube River (our data) 436 228 0.319987 104 376 Danube River, Black sea (Kolarov, 1983) 859.7 0.2416 0.405 Black Sea (Prodanov, Kolarov, 1983) 860.58 0.241 0.421 Table 5 Variations between the average weights of Pontic shad (Alosa immaculata) from different water bodies calculated at the same length by the use of length-weight relationship Water body and author Population lengthweight Average weights calculated at one and the same length equation W 5 W 10 W 15 W 20 W 25 Danube River our data W = 0.0234L2.7315 1.9 12.6 38.2 83.7 154.1 Kolarov (1983) W = 0.0629L2.574 4.0 23.6 67.0 140.4 249.4 Sea of Marmara, Ergüden et al. (2007) W = 0.0163L2.8511 1.6 11.6 36.8 83.5 157.7 Black Sea, Kalayci et al. (2007) W = 0.0046L3.1237 0.7 6.1 21.7 53.3 107.0 Black Sea, Yılmaz, Polat (2011) W = 0.0032L3.285 0.6 6.2 23.4 60.1 125.1 Black Sea, Yankova et al. (2011) W = 0.071L2.488 3.9 21.8 59.9 122.5 213.5 1978; 1980; 1983; 1985; Pavlov, 1953 and Ergüden et al., 2007). For other populations weight increments increase with increasing the fish age (Kolarov, 1978; Ciolac and Patriche, 2004). There is a linear regression between the average fish length at age t years (L t ) and the average length in one year (L t+1 ) (Figure 3). Based on this the linear growth of A. immaculata can be described with the von Bertalanffy s equation (1938) (Figure 4). reservoir (Raikova-Petrova and Zivkov, 1992). There is a linear regression between the average fish weight in the end of the t year (W t ) and the weight one year later (W t+1 ) (Figure 5). Based on this the weight growth of A. immaculata can be described with the von Bertalanffy s equation (1938) for the weight growth (Figure 6). L t = 35.749 [1 e -0.493209 (t-0.34115) ] The received value for the asymptotic growth was 35.749 cm. This value was lower than the maximal one obtained by us (37 cm) due to the fast temp of growth of the young fish and absence of adult individuals in the excerption. Similar result has been established for the zander (Stizostedion lucioperca) in Ovcharitsa cooling Observed Predicted Fig. 6. von Bertalanfy s growth equation curve describing the weight growth of Pontic shad (Alosa immaculata) from the Bulgarian sector of Danube River
124 D. Rozdina, G. Raikova-Petrova and P. Mirtcheva W t = 436.228 [1 e -0.319987 (t 1.04376) ] 2.7315 The obtained asymptotic weight for the conditions in Danube River was 436.228 g. The value of the coefficient of Hohendorf (1966) for the linear growth was 1.03 and for the weight growth was 0.8. The values close to one mean that the population of Pontic shad in Danube River best uses its growth potential. Comparing the temp of growth with the coefficient k from the von Bertalanffy s equation the population studied by us shows the fastest growth rates (Table. 4). Lowest temp of linear growth had the population studied by Kolarov (1980) in Danube River and lowest weight growth rates had the population studied by Prodanov and Kolarov (1983). Condition factor The received value for the coefficient of Fulton was 1.22 and for k b it was 2.8. Comparing to other populations the value we received for the coefficient of Fulton was one of the highest (Kolarov, 1965; 1978; 1980). K b for the population studied by us had higher value (2.8) than the one for the population in Marmara Sea (1.88) Ergüden et al. (2007b) and lower values than the population in Black Sea and Danube River (6.8) Kolarov (1983a) and the one in Black Sea (7.26) (Yankova et al., 2011) Calculated fish weight at the length of 5 g, 10 g, 15 g, 20 g and 25 g with the use of length-weight relationship is represented on Table. 5. This index had average values for the population studied by us. The populations studied by Kalayci et al. (2007) and Yilmaz and Polat (2011) had lower condition values while the highest condition was established from Kolarov (1983). Conclusions The studied population of A. immaculata is young with average age of 2.6 years. The most abundant are the fishes with size range 17 19 cm. The age structure is represented from 4 age groups (2 to 5 years old fish). According to von Bertalanffy s growth model in the Danube River conditions L = 35.749 cm and W = 436.228 g. Highest length and weight increments are established for the fish from the middle age groups. The values of the coefficient of Hohendorf close to one show that the population of A. immaculata in Danube River uses maximally well its growth potential. References Basami, R. and D. Grove, 1985. Studies on feeding, growth and production of a recruited inshore populations of Pleuronectes platessa (L.) at East Anglesey, North Wales. J. Fish Biol., 27: 765 783. Bertalanffy, L., 1938. A quantitative theory of organic growth. (Inquiries on growth laws. II). Human Biol., 10 (2): 182 213. Chugunova, N., 1959: Handbook on studying fish age and growth. Moscow, ASUSSR 156 p (Ru). Ciolac, A. and N. Patriche, 2004. Structure of danube shad (Alosa pontica Eichwald, 1838) spawner flocks migrating for reproduction in Danube River. Applied Ecology and Environmental Research, 2 (2): 53 58. De Silva, S., 1985: Body condition and nutritional ecology of Orefchromis mossambicus (Pisces, Cichlidae) populations of man-made lakes in Sri Lanka. J. Fish Biol., 27: 621 633. Ergüden, D., C. Turan and C. Çevik, 2007. The growth features of Pontic Shad Alosa pontica (Eichwald, 1838) in the Sea of Marmara, Turkey. Journal of Biol. Sc., 7 (4): 685 688. Eryilmaz, L., 2001. A study on the bony fishes caught in the south of the Sea of Marmara by bottom trawling and their morphologies. Turk J Zool., 25: 323 342. Goldspind, C., 1979. The population density, growth rate and production of roach, Rutilus rutilus, in Tjeukemeer, the Netherlands. J. Fish Biol., 15: 473 498. Hohendorf, K., 1966. Eine Diskussion der Bertalanffy Funktionen und ihre Anwendung zur Charakterisierung des Wachstums von Fishen. Kieler Meeresforsschungen, 22 (1): 70 90. Ivanov, L. and P. Kolarov, 1979. Relation between the catchments of Pontic shad (Alosa kessleri pontica, Eichw) and solar activity. Societas internationalis limnologiae SIL, XIX, Jubiläumstagung donauforschung, pp. 389 396 (Ru). Klayci, F., N. Samsun, S. Bilgin and O. Samsun, 2007. Length-weight relationship of 10 fish species caught by bottom trawl and midwater trawl from the middle Black Sea, Turkey. TrJFAS, 7: 33 36. Kolarov, P., 1958a. On the reasons for depletion of Pontic shad catchments on our territory. Ribno Stopanstvo, 1: 6 8 (Bg).
Age Composition and Growth Rate of the Spawning Part of the Population of Pontic Shad... 125 Kolarov, P., 1958b. On Pontic shad (Alosa kessleri pontica Eichw., Pisces) movements along Bulgarian Black Sea shore. Priroda, 1: 70 72 (Bg). Kolarov P., 1960a. Status and prospects of Pontic shad catchments along Bulgarian shore. Priroda, 4: 54 58 (Bg). Kolarov, P., 1960b. Some biological observations on Pontic shad (Alosa kessleri pontica Eichw.). SRIFFI, 2, pp. 13 34 (Bg). Kolarov, P., 1961. Regarding to the oncoming of the Pontic shad to the Bulgarian Black Sea shore. Ribno Stopanstvo, 6: 6 7 (Bg). Kolarov, P., 1963. Biological notes on Pontic shad caught along Bulgarian Danube River shore in 1962. Ribno Stopanstvo, 8 (Bg). Kolarov, P., 1964. Some features on age structure changes of Pontic shad (Alosa kessleri pontica Eichw.). Reports of the Institute of Aquaculture and Fisheries 5: 93 116 (Bg). Kolarov, P., 1965. On the biological characteristics of Alosa kessleri pontica Eichw. from the Bulgarian sector of Danube River. Reports of SIFO, 6: 87 87 (Bg). Kolarov, P., 1978. Some changes in Alosa kessleri pontica Eichw population during the spawning migration. Digest Limnologia na Balgarskiya Sector na Reka Dunav, pp. 238 340 (Bg). Kolarov, P., 1980. Particularities of Pontic shad (Alosa kessleri pontica Eichw.) in 1979 in Bulgarian aquatory. Fisheries, 27 (4): 17 19 (Bg). Kolarov, P., 1982. Dynamics of Pontic Shad (Alosa kessleri pontica Eichw.) Stocks. Reports of IFR, vol. XIX, Zemizdat, pp. 35 55 (Bg). Kolarov, P., 1983. Some basic parameters of the Pontic shad (Alosa kesleri pontica Eichw.) population. Hydrobiologiya, 19: 60 69 (Bg). Kolarov, P., 1985. Biological particularities and population dynamics of the anadromous fish species. Dissertation to obtain the Doctor of Science degree, 300 pp. (Bg). Kolarov, P., 1989. State and stocks of some Bulgarian Black sea anadromous fish species. Hydrobiologiya, 34: 50 58 (Bg). Le Cren, E. D., 1951. The length-weight relationship and seasonal cycle in gonad weight and condition in the Perch (Perca fluviatilis). J. Anim. Ecol., 20: 210 218. Pavlov, P. I., 1953. Biological and commercial characteristics of the spawning schools of the Pontic shad. Tr. Inst Hydrobiology AN USSR, vol. 28 (Ru). Prodanov, K. and P. Kolarov, 1983. On the question of the rational exploitation of the fish populations. Reports of IFR Varna, 20: 47 70 (Bg). Raikova-Petrova, G., 1992. Comparative population biology of the pike-perch (Stizostedion lucioperca L.) in Batak and Ovcharitsa dams. Dissertation to obtain the degree Doctor, 181 pp. (Bg). Raikova-Petrova, G. and M. Zivkov, 1993. Age and growth rate of pikeperch (Stizostedion lucioperca L.) in Ovcharitsa cooling reservoir. Hydrobiologiya, 38: 67 81 (Bg). Yankova, M., D. Pavlov, V. Raykov, V. Mihneva and G. Radu, 2011. Lenght-weight relationships of ten fish species from the Bulgarian Black Sea waters. Turk J. Zoo.l, 35 (2): 265 270. Yilmaz, S. and N. Polat, 2011. Length-weight relationship and condition factor of pontic shad, Alosa immaculata (Pisces: Clupeidae) from the southern Black Sea. Research Journal of Fisheries and Hydrobiology, 6 (2): 49 53. Zivkov, М., 1999. Factors, patterns and methodological importance of the populational-biology variability of the freshwater fishes. Dissertation to obtain the degree Doctor of Science, 406 pp. (Bg). http://e-ecodb.bas.bg/rdb/bg/vol2/alpontic.html; http://www.iucnredlist.org/apps/redlist/details/907/0 http://www.biodiversity.bg/files/file/zak_bg_biodiv.pdf