1. Introduction
Most lotic ecosystems in the world are subjected to various anthropogenic impacts, which together with the global warming processes disrupt the integrity of lotic ecosystems and adversely affect macrozoobenthic communities and their habitats [
8,
25]. As key communities in rivers, through which the main flow of matter and energy passes, they have proven their resilience and ability to quickly recover in harsh environments [
8,
11,
31]. The presence and abundance of macroinvertebrates are determined by their location and their ability to tolerate the specific environmental conditions in which they live [
32]. The species structure of the benthic communities is used as an indicator for their status and for the aquatic ecosystem health [
7,
11,
31].
Several studies on benthic macrofauna in Maritsa River and some of its tributaries have been published to date [
5,
24,
33,
34,
35,
37]. The taxonomic compositions and dominant structures of the macrozoobenthos in Maritsa River and estuarine zones of its main tributaries (Chepelarska, Sazliyka, Stryama, and Topolnitsa) recorded in 2020 were published by Park et al. [
20]. Varadinova et al. [
36] analysed the macroinvertebrate communities at the study sites in relation to various environmental factors and determined the ecological status of the studied aquatic ecosystems through the macrozoobenthos in 2021 and 2022.
Research on benthic macroinvertebrates in five major River Watershed Regions (RWRs: Han, Geum, Nakdong, Yeongsan, and Seomgin) has been published [
2,
6,
8,
13,
17]. The taxonomic composition and species diversity of benthic macroinvertebrate communities within the RWRs from 2010 to 2020 were identified, and their similarities as well as correlations with physicochemical variables were analysed [
14,
15]. Kwak et al. [
8] identified distributional characteristics based on the frequency of benthic macroinvertebrate taxa in the Han RWR (the South, North, and mainstream of Han River, and Anseong stream) during the period from 2009 to 2016, and analysed their relationship with altitude, current velocity, and substrate composition.
The aim of this study is to compare the composition and diversity of the macrozoobenthos in two selected rivers with similar latitude and hydromorphological features but located on different continents with distinct geographical characteristics: Maritsa River in Bulgaria and Han River in South Korea, considering different ecological situations and various anthropogenic impacts.
2. Material and Methods
Maritsa River, located in the East Aegean Sea River Basin District in Bulgaria and Han River, situated in the Han RWR in South Korea, were selected for this study. A total of 15 study sites from different sections of the main rivers and their tributary systems were sampled to ensure the representativeness [
20]. These were both unaffected reference sites and sites affected by diverse human activities. Some characteristics of the sampling sites were published by Park et al. [
19,
20]. Two sampling sites were located downstream in 2021 compared to their locations in 2020: site 13 in Maritsa River (4.9 km) and site 3 in Han River (17.4 km).
Benthic invertebrates from the different sites in Maritsa River were collected according to EU standards and national water legislation (standards BDS EN ISO 10870:2012 and EN 16150:2012) from 11 August to 17 September 2020 and on 23 and 24 July 2021. The macrozoobenthos surveys of Han River were conducted twice a year (spring and autumn) as part of the Water Environment Monitoring Network system, following the Guideline for Aquatic Ecosystem Status survey and health assessment method (2019-52) [
9]. The data is publicly available on the Water Environment Information System. Data from September and October (autumn data) in both 2020 and 2021 were used in this study.
The similarity in the taxonomic composition of the macrozoobenthos in the two rivers, was presented by the Sørensen similarity index [
29]. Cluster analysis using the Complete linkage method and Bray-Curtis similarity in the Statistical Software Package Primer 6 was also conducted. Margalef's index of species richness (d) [
12], Simpson’s index of dominance (c) [
28], Pielou’s index of evenness (е) [
21] and Shannon-Weaver diversity index (H) [
27] were used for the species diversity analysis.
4. Discussion
The numbers of identified taxa have decreased over the last several decades in both rivers. Even in the relatively short period of this study, a considerable reduction in species richness was observed in Maritsa River and a very slight decrease in Han River. A tendency of reducing the number of species recorded in Maritsa River have been observed for the last 65 years: from 305 species in 1955 [
21,
22] to 284 in 1966 [
23] and 229 in 1976-1977 [
35]. In 2021, the species richness (110 taxa) reduced by almost a third when compared to 1955. Even if we compare the total number of taxa recorded in Maritsa River in both 2020 and 2021 (192) the reduction is more than 1,5 times. This trend has been related to the deteriorating conditions of the aquatic environment, particularly to the increased pollution of the river over the years [
22,
23,
33,
34,
35], but some hydromorphological changes as well as the global warming can also contribute to that. In a study of Han River, conducted at similar locations, 120 species were recorded in 1994 [
1], but the number has been reduced twice for the last 30 years. Between 1986 and 2002, the number of taxa in the lower stretches of the Han RWR increased from 68 to 114 species, attributed to the Ecological Stream Restoration Project initiated in 1987 [
2,
19]. However, despite such efforts, the overall reduction in species richness couldn't be avoided, due to the rapidly expanding urbanization and increasing population.
In this study, the most diverse and dominant groups are Ephemeroptera and Trichoptera, in the composition of which there are both indicative of unaffected conditions species and species, which are tolerant to adverse influence (
Figure 1). In Bulgarian rivers, Diptera, Ephemeroptera, Oligochaeta, Plecoptera and Trichoptera have been commonly found [
3,
11,
18,
26]. However, in Maritsa River, the richness of the order Ephemeroptera, historically the most dominant group, was reduced considerably, from 51 species in 1955-1966 [
22,
23] to 16 species in 2021. In the main RWRs of South Korea, the phylum Arthropoda accounted for over 80% of benthic macroinvertebrates between 2009 and 2016: family Chironomidae had the largest number of taxa (20.8%), followed by the families Hydropsychidae (17.1%), Baetidae (12.6%) and Tubificidae (10.3%) [
8,
17]. The Han RWR exhibited the highest number of species and abundance among the main RWRs, with Chironomidae accounting for 18.6% (26 taxa), Hydropsychidae for 18.2% (25 taxa), and Baetidae for 14.8% (20 taxa). From 1986 to 2002, the most dominant groups in the Han RWR were Ephemeroptera and Trichoptera upstream, and Mollusca and Annelida downstream [
2]. In the recent five years, Chironomidae (31.95%), Ephemeroptera (30.12%), and Trichoptera (11.66%) were the most abundant and dominant groups in the Han RWR [
17]. In line with previous studies, the absence of Plecoptera species, known for their high sensitivity to pollution, is likely attributed to the anthropogenic impacts on the lotic ecosystems of Han River [
1,
2].
The low levels of similarity in the taxonomic composition of macrozoobenthos between the years and rivers (
Figure 3) might be due to variations in the hydrological characteristics recorded in the two years, leading to different water environmental conditions. It is notable that the reference sites 1 and 2 in both Maritsa and Han rivers formed independent but distinct clusters for the two lotic ecosystems. The reference sites, characterized by conditions closely resembling natural environments, exhibited among the highest species diversity in the rivers is the cause of the distinct separation.
Different benthic groups had the greatest species richness and abundance in Maritsa River in both years (
Figure 1 & 4). That fact together with the substantial decrease in the total number of taxa and the total abundance of macrozoobenthos in 2021 testified to a considerable change in the water environmental conditions in Maritsa River compared to 2020. The species richness and abundance of macrozoobenthic communities decreased at all sampling sites except for sites 8 and 13 where the reduced number of taxa, especially at site 13 (only 6 tolerant species remained) increased considerably their abundance. That resulted in strong dominance of one freshwater snail (class Gastropoda.) at site 8 and of two taxa of aquatic worms (class Oligochaeta) at site 13, and very low species diversity compared to the other sites. Sites 8 and 13 were located at the mouths of Vacha and Sazliyka rivers tributaries of Maritsa River [
19,
20,
36]. Their upstream stretches were characterized with good ecological status, however, downstream, they accumulated various anthropogenic pollutants, including toxins and domestic [
38,
39]. Despite self-purification processes, the lower stretches of the rivers remained polluted, especially when the river flows were reduced. Very strong dominance (
Figure 7) and very low species diversity were well exhibited also at site 6 in 2021 (
Figure 9).
The causes of the changes in the macrozoobenthic communities are difficult to precisely define due to the multiple factors and their correlations [
7]. Two main reasons could explain the differences between the two years. Firstly, the values of the most physicochemical parameters of the water environment were changed in 2021 [
36] and the macrozoobenthic communities responded to them with changes in the taxonomic composition and abundance. The indicative response of pollution-sensitive macrozoobenthic taxa has also been proven in other studies of Bulgarian rivers [
4,
30,
31]. Secondly, there was a difference in the sampling periods with a shorter sampling duration of two days in 2021 compared to a month-prolonged sampling in 2020. Macrozoobenthic communities change in dominance through the year by their specific life cycles and environmental conditions [
18]. The sampling period in 2021 might be before the blooming of certain macrozoobenthic species.
Although the number of taxa in Han River in 2020 and 2021 was almost the same, and the order Ephemeroptera dominated in terms of species richness and abundance in both years, the conditions in the water environment were not stable. Тhis was confirmed by the greatly changed taxonomic composition of the macrozoobenthos in 2021. The total abundance increased in that year, but mainly upstream (2-3 times) and at sites 13, 11 and 14 (22.5, 12.3 and 11.5 times, respectively). The physicochemical and biological parameters remained unchanged upstream (at sites 1 to 3). Downstream (at sites 12 to 15), the status based on physicochemical parameters got worse by one class, but the one assessed on macrozoobenthos got better at sites 11 to 13 and 15, and only had worsened at site 14 [
17]. The number of taxa as well as the abundance of the macrozoobenthic communities were very low at all sites from 6 to 15 in both years (except for the abundance at site 13 in 2021). Firstly, these stretches are characterized by relatively uniform hydro-morphological features, artificial structural elements such as concrete blocks and impermeable pavement along the riverbanks and riverbed. These elements isolate the river ecosystem from its surroundings, thereby creating an environment unsuitable to many benthic organisms [
10,
17]. Secondly, the downstream region is stagnant due to a slow and deep flow rate, which consequently limits the available habitat for macrozoobenthos [
16]. Thirdly, the inflow of pollutants from the central metropolitan area, particularly affecting sites 11 to 15, contributes to the dominance of contamination indicator taxa such as Chironomidae sp. and
Limnodrilus gotoi [
10,
16]. Historically, the ecological status downstream has consistently been 'poor' or 'very poor' [
17]. The small number of taxa, represented in limited quantities, had led to low species diversity and relatively close values of dominance and evenness as observed at sites 11 and 3 in 2021. In fact, some of the communities were almost destroyed completely.