NUTRIENT ENRICHMENT AND ITS EFFECT ON THE PHYTOPLANKTON COMMUNITY OF HRAZDAN RIVER IN THE YEREVAN DISTRICT, REPUBLIC OF ARMENIA Текст научной статьи по специальности «Биологические науки»

^äV.-: vi Трансформация экосистем issn 2619 оэ4х print

. ISSN 2619-0931 Online

/L •;.• Ecosystem Transformation www.ecosysttrans.com

Article

Nutrient enrichment and its effect on the phytoplankton community of Hrazdan River in the Yerevan District, Republic of Armenia

Lilit G. Stepanyan* , Evelina Kh. Ghukasyan

Institute of Hydroecology and Ichthyology of the Scientific Center of Zoology and Hydroecology, National Academy of Sciences of the Republic of Armenia, P. Sevak str. 7, 0014, Yerevan, Armenia

*[email protected]

Received: 29.03.2021 Revised: 17.05.2021 Accepted: 27.05.2021 Published online: 17.08.2021

DOI: 10.23859/estr-210329 UDC 574.52

Abstract. The phytoplankton community and nutrient enrichment of Hrazdan River have been studied in the Yerevan District, Republic of Armenia. Water was sampled in spring and summer of 2019. The content of phosphate, ammonium, and nitrite ions were significantly higher downstream of the Yerevan District compared to the stations located upstream. The latter were characterized by higher abundance of diatoms (Melosira vanans and Rhoico-sphenia curvata), while Oscillatoria limnetica, O. chlorina, and Anabaena sphaerica dominated downstream, reflecting different eutrophic state of the river parts.

Key words: Hrazdan River, phytoplankton community, nutrient enrichment.

To cite this article. Stepanyan, L.G., Ghukasyan, E.Kh., 2021. Nutrient enrichment and its effect on the phytoplankton community of Hrazdan River in the Yerevan District, Republic of Armenia. Ecosystem Transformation 4 (3), 3-12. https://doi.org/10.23859/estr-210329

Introduction

Rivers passing through megacities are exposed to high threat by various anthropogenic activities. Particularly, industrial and household wastewater discharge is the common source for organic load into the aquatic ecosystems (Adeyemo et al., 2008). Such a load poses huge risks for all aquatic organisms and, in particular, for phytoplankton communities. The distribution of different algae species depends on interactions between environmental abiotic and biotic factors. Deterioration of environmental conditions due to increased organic matter income leads to obvious changes in phytoplankton community, promoting succession from diatoms to Cyanophyta. Therefore, the use of phytoplankton as an indicator for ecological status of surface water bodies is widely

accepted (Reynolds et al., 2002; Swaminathan, 2003).

Hrazdan River is one of the major tributaries of Araks River on the territory of the Republic of Armenia. It flows from Lake Sevan through the Sevan, Hrazdan, Charentsavan, Abovyan, and Yerevan cities; its length is 141 km. However, negative effect of Yerevan city on ecological status of the river exceeds the cumulative negative effect of all the settlements upstream. The water of Hrazdan River is used for irrigation, energetic, industrial, and other purposes (Chilingaryan et al., 2002, Stepanyan, 2009).

The study aims to assess the nutrient enrichment and corresponding changes in phytoplankton community of Hrazdan River on the territory of Yerevan District.

Table 1. Brief characteristics of sampling sites.

Station no. Latitude Longitude

1 N 40.280873° E 44.589203°

2 N 40.171° E 44.499669°

3 N 40.159089° E 44.1389°

4 N 40.148791° E 44.589203°

Description of sampling site

Getamej village, upstream Yerevan City

City center of Yerevan

Upstream Yerevanyan Lich Reservoir, downstream the confluence site of Djhrvejh and Getar tributaries

Verin Charbakh district, downstream of Yerevanyan Lich Reservoir

Table 2. The hydrophysical parameters and water pH in Hrazdan River at the territory of Yerevan District. T - water temperature, V - flow velocity.

Station

Parameters 12 3 4

Spring Summer Spring Summer Spring Summer Spring Summer

T, °C 13 16 14 21 15 20 15 25

V, m/c 0.3 0.6 0.7 0.5 0.2 0.4 1 0.6

pH 7.5 7.2 8.2 8.5 8.1 7.9 8.0 8.2

Material and methods

Study was carried out during spring (May) and summer (July) seasons of 2019. Water was sampled at four stations for hydrochemical and hydrobiologic analyses (Table 1).

The phytoplankton community was analyzed by the standard methods (Abakumov, 1983). A 1-L water sample was taken from each site and immediately fixed with 40% formaldehyde solution (0.4% final concentration) and stored in a dark place until the total sedimentation of seston. The laboratory analyses were performed according to V.A. Abakumov (1983).

Quantitative and qualitative analyses of phyto-plankton were performed under a microscope in the 0.01-mL Nageotte chamber. Phytoplankton biomass was calculated by the stereometric (cell volume) method.

The species of planktonic algae were identified by the taxonomic keys and the manuals for freshwater ecosystems (Hambaryan and Shahazizyan, 2014; Proshkina-Lavrenko and Makarova, 1968; Streble and Krauter, 2001; Tsarenko, 1990).

Water temperature and velocity were measured at the sampling site. Hydrochemical parameters such as DO (dissolved oxygen) and pH were measured in parallel with the hydrophysical parameters by Milwaukee (PH 56 PRO) Waterproof pH meter and Milwaukee (MW 600) Dissolved Oxygen meter. The nutrients' concentrations were measured in the laboratory

according to the ISO methods by spectrophotometer. Ammonium nitrogen concentration was determined according to the ISO 5664:2006, nitrite nitrogen, ISO 6777:1984, nitrate nitrogen, ISO 6777:1984, and phosphate phosphorus, ISO 6878.

Results and discussion

Abiotic parameters

The hydrophysical parameters are presented in Table 2. Water temperature ranged from 13 to 25 °C, flow velocity, from 0.2 to 1.0 m/s.

According to A.M. Nikanorov (2001), the optimal pH range for sustainable aquatic life varies from 6.5 to 8.5. Thus, optimal pH conditions have been registered at all stations of the studied part of Hrazdan River (Table 2).

Generally, the dissolved oxygen concentration (DO) was lower in summer than in spring. Moreover, DO was significantly higher upstream Yerevan city than at the other stations. The lowest DO concentration (2.9 mgO2/L) was registered at the station no. 3 in summer, partly due to the urban sewage discharge brought to this site by the Getar River, a tributary of the Hrazdan River (Fig. 1).

Similarly to 2003-2006 (Stepanyan, 2009), significant nutrient enrichment has recorded in the Hrazdan River in 2019. Nitrate content was high at the stations nos. 1 and 3, which was still within the

Fig. 1. Dissolved oxygen concentration at the sampling sites. Hereinafter, the background environmental values are indicated according to State Standards of Republic of Armenia.

mgN/l

3.5

2.5

1.5

0.5

lll.il-Spring Summer Spring Summer Spring Summer Spring Summer 12 3 4

Ecological norm

Fig. 2. Nitrate concentrations at the sampling sites.

ecological norm standards for the river ecosystem health and biodiversity conservation1 (Fig. 2).

The concentrations of nitrite and ammonium ions in the riverine waters usually varied from hundredths to tenths of a milligram per liter (Figs. 3, 4). Within our

1 Government Decision No. 75-N "On defining the standards for water quality of each water basin management area depending on local characteristics". Web page. URL: https://www.e-gov.am/u_ files/file/decrees/kar/2011/02/11_0075.pdf (accessed: 27.01.2021)

study, nitrite content varied from 0.01 to 0.32 mgN/L. In summer, nitrite content significantly increased at the stations located within the Yerevan city (stations nos. 2, 3, and 4).

Ammonium concentration also exceeded ecological norm at the stations nos. 3 and 4 at both seasons (Fig. 4). High nitrite and ammonium content corresponds to the eutrophic status of the waters (Nikanorov, 2001).

Phosphate concentrations were exceptionally high in summer almost at all sampling sites except

mgN/l

0.35 0.3 0.25 0.2 0.15 0.1 0.05 0

_ _ __ H __

Spring Summer 1 Spring Summer 2 Spring Summer 3 Spring Summer 4

Fig. 3. Nitrite concentrations at the sampling sites.

■Ecological norm

Spring Summer 1

Spring Summer 2

Spring Summer 3

Spring Summer 4

■Ecological norm

Fig. 4. Ammonium concentrations at the sampling sites.

station no. 4 (Fig. 5). Such pattern observed downstream Yerevanyan Lich Reservoir could be explained by the algal bloom in the reservoir itself. As the algae consume phosphates, a bloom may minimize the phosphate content in the outflowing water.

Phytoplankton community

Earlier studies of phytoplankton carried out in 2004-2006 in Hrazdan River at the territory of Yerevan District revealed 108 species belonging to six phytoplankton groups (Stepanyan, 2009).

During present study of phytoplankton community, 55 species belonging to six phytoplankton groups have been registered (Table 3). As found generally, diatoms dominated, but Cyanophyta were characterized by higher abundance than usually observed. Although some species of Dinophyta and Euglenophyta, new groups of algae for Hrazdan River, were recorded in 2019, total species diversity of phytoplankton decreased comparing to previous studies (Badalyan et al., 2005; Stepanyan et al., 2005; Stepanyan, 2009).

Table 3. List of phytoplankton species at the sampling sites. "+" - species is present; "-" - species is absent. Chorological types (Geo): c - cosmopolite; b - boreal, a-a - Arctic-Alpine. Halobity (Hal): mh - mesohalobe; i - oligohalobe-indifferent; hl - oligohalobe-halophilous; hb - oligohalobe-halophobous. Habitat types (Hab): B - benthic; P - planktonic; P-B - planktonic-benthic. Saprobity (S): o - oligosaprobic; o-p - oligo-beta-mesosaprobic; p - beta-mesosaprobic; p-o - beta-oligomesosaprobic; p-a - beta-alfa-mesosaprobic; a-p - alfa-beta-mesosaprobic; x - xenosaprobic; x-o - xeno-oligosaprobic; x-p - xeno-beta-mesosaprobic; o-x - oligo-xenosaprobic, o-a - oligo-alfa-mesosaprobic; p - polysaprobic; "-" - no data (Barinova et al, 2006).

12 3 4

(spring/ (spring/ (spring/ (spring/ Geo Hal summer) summer) summer) summer)

Hab S

Cyanophyta

Anabaena sphaerica Bornet & Flahault

Aphanothece clathrata West & G.S. West

Microcystis aeruginosa Kutzing

M. wessenbergii Komarek in Joosen

Oscillatoria chlorina Kützing ex Gomont

Oscillatoria limnetica Lemmermann

Bacillariophyta

Achnanthes taeniata Grunow

Amphora ovalis (Kützing) Kützing

Ceratoneis arcus (Ehrenberg) Kützing

Cocconeis placentula Ehrenberg

C. pediculus Ehrenberg Cyclotella comta Kützing

C. stelligera (Cleve & Grunow) Van Heurck

Cymatopleura solea (Brebisson) W. Smith

Cymbella prostrata (Berkeley) Cleve

C. ventricosa C. Agardh

Diatoma hiemale var. hiemale (Roth) Heib

D. vulgaris Bory

Fragilaria capucina Desmazieres

F. construens (Ehrenberg) Grunow

F. crotonensis Kitton

-/-

+/+ +/+ +/+

-/--/-/+ -/+ -/-

+/+ -/+ +/-

+/-/--/-+/+ +/+ +/-+/+ +/+ +/-

-/+/--/+ +/+ -/--/-

-/--/--/-/+ -/--/-

-/--/--/--/+/-/-/+ -/+/-

-/-+/+ +/+ +/+ -/+ -/+

+/+ +/-+/-

-/+ -/--/-

-/-/+ -/--/-+/+ -/-+/+ -/-+/+

-/+ +/+ +/+ +/+ +/+ -/+

+/--/+ -/-

+/+ -/+ -/-

-/--/-/+ -/+ +/+ -/-+/+ -/+ +/+

ci ci c hl c-

- hl ci a-a i

ci ci ci

P

hb

o-ß ß

o-a

o-a p

o-ß

P P P P

P-B P-B

B -

B a-ß

B o-x

P-B o-ß

B o-a

ß-o

hl

P-B x

P-B o

B o-ß

B o-a

P-B ß-o

P-B ß Bo

P-B o

P a-ß

c

c

c

c

c

c

c

c

c

1 (spring/ summer) 2 (spring/ summer) 3 (spring/ summer) 4 (spring/ Geo Hal summer) Hab S

Gomphonema constrictum Ehrenberg in Kützing -/+ -/- -/+ -/- c i B o

Gomphonema olivaceum (Hornemann) Brebisson +/- -/- +/- -/- c i B ß-a

Melosira varians C. Agardh +/+ +/+ +/- +/+ c hl P-B a-ß

Navicula cryptocephala Kutz. +/- -/- +/+ +/+ ci B x-o

N. menisculus Schumann -/- +/- -/- -/- ci B x-ß

N. pupula Kützing +/- -/- -/- -/- c hl B x-o

N. pygmaea Kützing -/- -/- -/+ -/- c mh - ß-o

N. radiosa Kutzing -/- -/- +/- -/- ci B o

N. rhyncocephala Kutz -/- -/- +/- -/- c hl B ß

Nitszchia amphibia Grunow -/- -/+ -/- -/- ci P-B o

N. dissipata (Kützing) Rabenhorst -/- -/- +/- -/- ci B x

N. linearis W. Smith -/- -/- -/- +/- ci B x

Pinnularia virdis Ehrb. +/+ +/+ +/+ -/+ ci P-B o-x

P. leptosoma (Grunow) Cleve +/+ -/- +/- -/- b i B o

Rhoicophenia curvata (Kutzing) Grunow +/- +/+ -/+ -/+ ci P-B x-o

Stephanodiscus astraea (Kützing) Grunow +/- +/+ -/+ +/+ ci P ß

S. hantzschii Grunow in Cleve & Grunow -/- -/- +/- -/- ci P a-ß

Surirella biseriataBrebisson in Brebisson & Godey +/- -/- -/- -/- ci - o-ß

S. ovata Kützing +/- -/- -/+ -/- ci B o-a

Surirella sp. -/- +/- -/- -/- -- - -

S. tenera var. tenera Greg. -/- -/- -/+ -/- ci P-B o

Chlorophyta

Ankistrodesmus falcutus (Corda) Ralfs -/- -/- -/- +/- c hb P-B ß

Binuclearia lauterbornii (Schmidle) Proschkina-Lavrenko -/- -/+ -/- -/- -- - -

Characium nasutum Rabenhorst +/- -/- -/- -/-

Chroococcus sp. -/- +/- -/- -/- -- - -

12 3 4

(spring/ (spring/ (spring/ (spring/ Geo Hal Hab S summer) summer) summer) summer)

Dictyosphaerium pulchellum H.C. Wood

Pandorina morum (Mull) Bory

Scenedesmus armatus var. armatus (Chod) G.M. Smith

Scenedesmus opolensis var. opolensis P. Richt.

Selenastrum gracile Reinsch

Trebouxia humicola G.S. West & F.E. Fritsch

Euglenophyta

Trachelamonas volvocina (Ehrenberg) Ehrenberg

Dinophyta

Peridinium willei Huitfeldt-Kaas Xanthophyta

Tribonema monochloron Pascher & Geitler in Pascher

-/--/--/--/-

-/--/+/-/--/-

-/-/+ -/--/-/+ -/-

-/--/--/-

-/--/--/--/--/--/-

-/--/--/-

-/+ -/-/+ -/+ -/+/-

-/+

-/+

i P

i P

ß

P-B o-a P-B ß P-B o-a

+/+ c i B ß

hl P o-ß

c

ß

c

c

Microalgal communities of the Hrazdan River are represented by benthic, planktonic, and planktonic-benthic algae, mostly by benthic algae (Table 3). Forty-five phytoplankton species are indicators of organic pollution, presented mostly by p-mesosaprobic species. Four ecological groups of salinity indicators have been identified, the representatives of oligohalobe-indifferent group predominate. Most species are eurybionts.

In spring, diatoms strongly prevailed at the station no. 1, comprising 96% of phytoplankton community. Melosira varans dominated both by abundance (45%) and biomass (53%).

Downstream the station no. 1, there was a decrease in the quantitative parameters of phytoplankton, the dominant groups of algae changed as well. Cyanophyta prevailed at the stations nos. 3 and 4. Correspondingly, Microcystis aeruginosa and Aphanothece clathrata were the most abundant species at these stations. Representatives of Chlorophyta and Euglenophyta were also revealed in spring at the study sites.

Compared to spring, both the quantitative and qualitative parameters of microalgae increased in summer (Figs. 6, 7). However, diatoms still dominated at the stations nos. 1 and 2 due to Melosira varians and Rhoicosphenia curvata, respectively.

Significant increase of Cyanophyta abundance was registered at the stations nos. 3 and 4 following the nutrients' dynamics. Bloom of Oscillatoria species (O. limnetica and O. chlorina) was observed; they comprised 77% of the total abundance and 59% of the total biomass at the station no. 3; 75% and 59% at the station no. 4, respectively.

Oscillatoria species reach considerable abundance in naturally eutrophic and polluted streams and lakes (Babanazarova et al., 2013; Prowse, 1969). They are generally benthic by preferable habitat, but may develop actively in plankton communities. Both registered species of this genus refer to the indicators of organic pollution (Table 2) (Barinova, 2006).

Ability of Cyanophyta to produce toxins attracts much attention, especially when these species dominate in eutrophic waters, making concern to water quality (Codd, 2000). Oscillatoria species are the most common agents of freshwater toxic blooms. They produce hepatotoxin microcystin (Chia et al, 2018). Microcystin is reported to cause death of wild animals and agricultural livestock (Carmichael,

1988), it has been as well recognized as a potential threat to the human health in the countries, where water supplies are contaminated with cyanobacteria (Banerjee, 2020; Gkelis and Zaoutsos, 2014; Yu,

1989). Therefore, the bloom of Oscillatoria species

mgP/l

0.3

0.25 0.2 0.15 0.1 0.05 0

Spring Summer 1

Spring Summer 2

Spring Summer 3

Spring Summer 4

Ecological norm

Fig. 5. Phosphate concentrations at the sampling sites.

E

D>

V.

mß

Spring | Summer Spring | Summer 1 2

Spring | Summer 3

Spring | Summer 4

■ Cyanophyta DBacillariophyta DChlorophyta BOthers Fig. 6. Seasonal dynamics of phytoplankton abundance at the sampling sites.

2500000 2000000 S 1500000 1000000 H 500000 0

a> o

fe? S> ö OT Ö

Spring | Summer 1

Spring | Summer 2

Spring | Summer 3

Spring | Summer 4

■ Cyanophyta nil Bacillariophyta DChlorophyta B Others Fig. 7. Seasonal dynamics of phytoplankton biomass at the sampling sites.

could cause a serious harm to the economy of the agriculture and recreational capacities of the area, since these sectors of economics are the main water consumers.

Abundance of nitrogen-fixing species Anabaena sphaerica (148,000 cell/L) were observed at the station no. 4. Our observations performed in Yerevanyan Lich reservoir during summer of 2019 evidenced on the dominance of Oscillatoria chlorina and Anabaena sphaerica, proving the impact of the reservoir on the Hrazdan River (Annual Report..., 2019).

Conclusion

As a corollary, we conclude that high content of phosphates and nitrogen lead to increase of algae growth rates and ultimately reduce dissolved oxygen concentration in water. Nutrient enrichment at the stations nos. 3 and 4 is significantly higher than at the sites located upstream. The Oscillatoria bloom has been registered in the stations nos. 3 and 4 in summer. The influence of Yerevanyan Lich reservoir on phytoplankton community of these stations is obvious. Therefore, better management and strict control of the nutrient load to Hrazdan River from Yerevan city is urgently needed.

ORCID

L.G. Stepanyan 0000-0002-0960-7234

E.Kh. Gukasyan 0000-0001-8866-9785

References

Abakumov, V.A., 1983. Rukovodstvo po metodam gidrobiologicheskogo analiza poverkhnostnykh vod i donnykh otlojenij [Manual for hydrobiological analysis of surface water and sediments]. Gidrometeoizdat, Leningrad, USSR, 240 p. (In Russian).

Adeyemo, O.K., Adedokun, O.A., Yusuf, R.K., Adeleye, E.A., 2008. Seasonal changes in physico-chemical parameters and nutrient load of river sediments in Ibadan city, Nigeria. Global NEST Journal 10 (3), 326-336.

Annual Reports of the Institute of Hydroecology and Ichthyology of Scientific Center of Zoology and Hydroecology NAS RA 2019. Report on major research outcomes for 2019, issued by decision of the Presidium of NAS RA, N1 (1768) on 18 March of 2020. 66-69. (In Armenian).

Babanazarova, O., Sidelev, S., Schischeleva, S., 2013, The structure of winter phytoplankton in Lake Nero, Russia, a hypertrophic lake dominated by Planktothrix-like Cyanophyta. Aquatic Biosystems 9 (18). https://doi.org/10.1186/2046-9063-9-18

Badalyan, K., Stepanyan, L., Hambaryan, L., 2005. Study of some spaces of Cyanophyta in the phytoplankton of the River Hrazdan. Abstract's book of 3rd Young Medics' International Conference. Armenian medical association, Yerevan, 37.

Barinova, S.S., Medvedeva, L.A., Anisimova, O.V., 2006. Bioraznoobrazie vodoroslej - indikatorov okrughajushej sredy [Biodiversity of algal indicators in environmental assessment]. Pilies Studio, Tel-Aviv, Israel, 497 p. (In Russian).

Banerjee, A., 2020. Liver toxicity of chemical warfare agents. In: Gupta, R. (ed.), Handbook of Toxicology of Chemical Warfare Agents. Academic Press, Cambridge, USA, 659-671. https://doi. org/10.1016/B978-0-12-819090-6.00040-4

Carmichael, W.W., 1988. Toxins of freshwater algae. In: Tu, A.T. (ed.), Handbook of Natural Toxins, Marine Toxins and Venoms. Marcel Dekker, New York, USA, 121-157.

Chia, M.A., Jankowiak, J.G., Kramerb, B.J., Goleskib, J.A., Huang, I-Sh., Zimbac, P.V., Bittencourt-Oliveiraa, M. do Ca., Gobler, C.J., 2018. Succession and toxicity of Microcystis and Anabaena (Dolichospermum) blooms are controlled by nutrient-dependent allelopathic interactions. Harmful Algae 74, 67-77.

Chilingaryan, L.A., Mnacakanyan, B.P., Axababyan, K.A., Tokmajyan, H. V., 2002. Hydrography of rivers and lakes of Armenia. Agropress, Yerevan, Armenia, 49 p. (In Armenian).

Codd, G.A., 2000. Cyanophytal toxins, the perception ofwater quality, and the prioritisation of eutrophication control. Ecological Engineering 16, 51-60.

Gkelis, S., Zaoutsos, N., 2014. Cyanotoxin occurrence and potentially toxin producing Cyanophyta in freshwaters of Greece: a multidisciplinary approach. Toxicon 78, 1-9.

Hambaryan, L., Shahazizyan, I., 2014. Brief decisive, educational manual for genera of freshwater algae. YSU Press, Yerevan, Armenia, 61 p. (In Armenian).

Nikanorov,A.M., 2001. Gidrokhimia [Hydrochemistry]. Gidrometeoizdat, Saint-Petersburg, Russia, 447 p. (In Russian).

Proshkina-Lavrenko, A.I., Makarova, I.V., 1968. Vodorosli planktona Kaspijskogo morja [Planktonic algae of the Caspian Sea]. Nauka, Leningrad, USSR, 295 p. (In Russian).

Prowse, G.A., 1969. The role of cultured pond fish in the control of eutrophication in lakes and dams. Internationale Vereinigung für Theoretische und Angewandte Limnologie: Verhandlungen 17 (2), 714-718.

Reynolds, C.S., Huszar, V., Kruk, C., Naselli-Flores, L., Melo, S., 2002. Towards a functional classification of the freshwater phytoplankton. Journal of Plankton Research 24, 417-428.

Stepanyan, L.G., 2009. Hydrobiological and hydrochemical study of Hrazdan hydro ecosystem in the territory of Yerevan district. Biological sciences PhD thesis. Yerevan, Armenia, 142 p.(In Armenian).

Stepanyan, L.G, Hambaryan, L.R., Hovhannisyan, R.H., 2005. Study of phytoplankton community dynamics of the Hrazdan River in Yerevan area. Biological Journal of Armenia 57 (3-4), 257-264.

Streble, H., Krauter, D., 2001. Das Leben im Wassertropfen. Kosmos, Stuttgard, Germany, 415 p. (In German).

Swaminathan, M.S., 2003. Biodiversity: an effective safety net against environmental pollution. Environmental Pollution 126, 287-291.

Tsarenko, P.M., 1990. Kratkij opredelitel chlorokokovykh vodoroslej Ukrainskoj SSR [Short key of Chlorococcal Algae of Ukraine]. Naukova Dumka, Kiev, USSR, 106 p. (In Russian).

Yu, S.Z., Tang, Z.Y., Wu, M.C., Xia, S.S., 1989. Drinking water and primary liver cancer. In: Tang, Z.Y., Wu, M.C., Xia, S.S. (eds.), Primary Liver Cancer. China Academic Publisher, Beijing, China, 30-37.