Next Article in Journal
Assessment of Surface Water Quality in the Baia Mare Area, Romania
Next Article in Special Issue
A Digital Light Microscopic Method for Diatom Surveys Using Embedded Acid-Cleaned Samples
Previous Article in Journal
Effects of Artificial Water Withdrawal on the Terrestrial Water Cycle in the Yangtze River Basin
Previous Article in Special Issue
The Impacts of Different Anthropogenic Disturbances on Macroinvertebrate Community Structure and Functional Traits of Glacier-Fed Streams in the Tianshan Mountains
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Odonata Assemblages in Anthropogenically Impacted Habitats in the Drava River—A Long-Term Study

by
Marina Vilenica
1,* and
Zlatko Mihaljević
2
1
Faculty of Teacher Education, University of Zagreb, Trg Matice Hrvatske 12, 44250 Petrinja, Croatia
2
Department of Biology, Faculty of Science, University of Zagreb, Rooseveltov Trg 6, 10000 Zagreb, Croatia
*
Author to whom correspondence should be addressed.
Water 2022, 14(19), 3119; https://doi.org/10.3390/w14193119
Submission received: 10 September 2022 / Revised: 28 September 2022 / Accepted: 29 September 2022 / Published: 3 October 2022

Abstract

:
Lotic freshwater ecosystems are among the most threatened ecosystems worldwide due to the effects of multiple stressors, such as intensive land use in their catchments, morphological alterations, flow regulation, pollution, and climate change. Odonata are often used as valuable indicators of ecological integrity and anthropogenic disturbance of freshwater habitats. Here, we present the results of a study on Odonata assemblages in anthropogenically impacted habitats (hydropower plant reservoirs, tailrace canals, drainage ditches, and old river channels) conducted over a nine-year period. The negative impacts of anthropogenic activities on inhabiting biota were confirmed—with only 11 species recorded, the Odonata assemblages were species-poor and had low population densities. Although most species recorded were generalists, some species of national conservation concern were detected. Among the physico-chemical water parameters, the concentrations of ammonium, orthophosphates, nitrates, and mineral oils in the water were found to be the most important determinants of Odonata assemblages. The preservation of near-natural sites in the vicinity of anthropogenically impacted and man-made habitats is important for maintaining the local Odonata fauna and for the preservation of rare species. Our results highlight the importance of long-term data for determining the occurrence of Odonata species and monitoring their population dynamics.

1. Introduction

Freshwater ecosystems cover only about 0.8% of the land’s surface and account for no more than 0.01% of the world’s water volume [1]. However, they are of enormous importance for life on Earth and are characterized by impressive biodiversity [2,3], harboring more than 100,000 species [1]. They also provide various ecosystem services to humans, i.e., they are used as a source of drinking water, for irrigation in agriculture, for industrial purposes, for energy production, for recreation, etc. [4,5]. Nevertheless, freshwater ecosystems are among the most degraded natural systems in the world [6,7]. They are threatened by many anthropogenic stressors, such as pollution, hydromorphological alterations, habitat fragmentation, dam construction, flow regime alterations, introduction of invasive species, urbanization, and the overexploitation of freshwater resources due to population growth [1,8,9,10]. In addition, the impacts of climate change on freshwater ecosystems are increasing significantly, and due to rapid changes in air temperatures and precipitation patterns, they have already heavily modified the hydrological regimes of rivers worldwide [8,11]. Some regions are affected by severe droughts, while others are facing heavy floods [12]. Such stressors negatively impact aquatic communities, leading to their homogenization and shifts in taxonomic and functional species composition, resulting in a dramatic loss of freshwater biodiversity [13,14,15,16,17,18,19,20,21].
Many interacting factors influence subtle ecological processes that occur over a relatively long period of time. Therefore, these processes, such as population dynamics, life cycles, and reproductive patterns, as well as the influences of various disturbances on communities, require long-term studies to be understood [22]. Long-term data could also explain the limiting factors for different species, leading to a better interpretation of their past, present, and future distributions and understanding species’ resilience and responses to environmental stressors [23]. Therefore, long-term studies could serve as a useful tool to assess the impact of global change on biota [14]. Although the number of long-term studies on European benthic macroinvertebrates has increased in the 21st century, the availability of long-term data is still quite limited, including those on Odonata [14,24].
Odonata is an order of aquatic insects with an amphibious life cycle that includes an aquatic nymph phase and a terrestrial adult phase. Due to their predatory feeding habits, both nymphs and adults are at high levels in the trophic web [25]. In addition, the life history traits of different Odonata species vary widely, such as their dispersal ability, habitat and microhabitat preferences, and tolerance to water pollution and habitat degradation [26]. For these reasons, they are widely used as valuable indicators of ecological integrity and the anthropogenic disturbance of freshwater habitats [20,21,27,28]. The bottom substrate and submerged aquatic vegetation composition, as well as water clarity, are the most important determinants of Odonata nymph occurrence in a given habitat, while adults select habitats primarily based on structural heterogeneity, i.e., aquatic and riparian vegetation structure and shading [25,29,30]. The main objectives of this study were: (i) to determine Odonata species richness and their spatial and temporal distributions, and (ii) to determine the main environmental factors affecting Odonata occurrence in anthropogenically altered and man-made habitats in the Drava River system over a nine-year period. Species-rich Odonata assemblages were expected to be found in less anthropogenically impacted habitats.

2. Materials and Methods

2.1. Study Area

With a length of 719 km, the Drava River is the fourth-largest and longest tributary of the Danube. The river originates in Italy and flows through Austria, Slovenia, Hungary, and Croatia, encompassing a wide range of natural habitats along its middle and lower reaches that harbor unique animal and plant assemblages [31,32,33]. For these reasons, the Drava River has been a part of the Transboundary Biosphere Reserve Mura–Drava–Danube since 2012 [34]. However, the river has been heavily modified along its entire length by the construction of hydropower plants (HPPs), the last three of which, namely Varaždin, Čakovec, and Dubrava, are located in the north-western part of Croatia. After the construction of the HPPs, new man-made habitats were created (e.g., reservoirs, tailrace canals, and drainage ditches) that strongly influenced the water regime of the river [35].
This study was conducted at 17 sites along the entire hydropower system (Figure 1). Over a length of approximately 60 km, these sites belong to all of the major habitat types of the system (reservoir, tailrace canal, drainage ditch, and old river channel). Four sites belong to the Varaždin HPP, seven to the Čakovec HPP, and six to the Dubrava HPP (Figure 1).

2.2. Odonata Sampling

Odonata nymphs were sampled seasonally in April, July, October, and December over a nine-year period (2012–2020), together with other macroinvertebrates. Reservoirs (sites V1, V2, C1, C1A, C2, D1, and D2) were sampled using an Ekman grab with a sample area of 225 cm2, with three subsamples collected. The sampled macroinvertebrates were sorted out from the sediment by sieving with a benthos hand net (mesh size of 475 μm). At other sites, which included tailrace canals (sites V3, C5, and D5), drainage ditches (sites C3, C4, D3, and D4), and old river channels (sites V4, C6, and D6), macroinvertebrates were sampled using a benthos hand net (25 cm × 25 cm; mesh size = 500 µm). Ten subsamples were collected at each site and at each sampling event, and pooled into one composite sample. Odonata nymph abundance is presented as the number of individuals per square meter.
Nymphs were identified using the relevant identification keys [36,37,38]. Very young or damaged individuals were identified to the family level. All voucher specimens were deposited at the Department of Biology, Faculty of Science, Zagreb, Croatia.

2.3. Environmental Variables

At each study site, the following environmental parameters were measured at the time of macroinvertebrate sampling: water temperature, pH (using a WTW pH 330 pH meter), dissolved oxygen concentration and oxygen saturation (using a WTW Oxi 330/SET oximeter), and conductivity (with a WTW LF 330 conductivity meter). The rest of the environmental parameters (biological oxygen demand, chemical oxygen demand, ammonium, nitrite, nitrate, total nitrogen, total phosphorous, orthophosphate, pesticide, phenol, and mineral oil concentrations in water, as well as the number of fecal coliform bacteria in 1000 mL) were analyzed according to the standard methods [39] (Table 1).

2.4. Data Analyses

Cluster analysis was performed to reveal similarities in the Odonata assemblages among the study sites. The analysis was based on the Bray–Curtis similarity index. Study sites without Odonata records were excluded from the analysis.
To determine the relationship between Odonata assemblages and environmental variables, a redundancy analysis (RDA) with the Monte Carlo permutation test for significance (with 499 permutations) was performed. This analysis was conducted using data for 13 taxa, 153 sampling points, and 15 environmental variables (water temperature, pH, conductivity, dissolved oxygen concentration, biological oxygen demand, chemical oxygen demand, ammonium, nitrate, nitrite, total nitrogen, orthophosphate, total phosphorus concentrations, concentrations of pesticides, mineral oils, and counts of fecal coliform bacteria in water).
Bray–Curtis similarity index and cluster analyses were performed using the Primer 6 software package [40], and RDA analysis was performed using CANOCO for Windows (ver. 4.02) [41]. Species data were log(x + 1)-transformed prior to analyses.

3. Results

3.1. Environmental Variables

The mean values of environmental parameters slightly varied between habitat types (Table 1). The mean water temperature, ammonium, orthophosphate, and phenol concentrations were slightly higher in the old river channels than in the other habitat types. The mean dissolved oxygen concentration, oxygen saturation, biological and chemical oxygen demand, nitrate, total nitrogen, total phosphorous, and pesticide concentrations were slightly higher in the reservoirs than in the other habitat types. The mean conductivity, nitrite concentration, and fecal coliform bacteria counts in water were higher in the drainage ditches, while the mineral oil concentrations were higher in the tailrace canals than in the other habitat types. The mean pH was higher in the reservoirs and tailrace canals than in the other two habitat types.

3.2. Odonata Assemblages

During the study period, a total of 12 taxa of Odonata (11 species) was recorded (Table 2). Platycnemis pennipes and Onychogomphus forcipatus were the most widespread species, detected at eight and six sites, respectively (Table 2). Old river channels had the highest Odonata abundance (number of individuals per m2) and species richness (number of species) (sites C6, D6, and V4). No Odonata species were detected at several sites, including four reservoirs (C2, D2, V1, and V2) and two drainage ditches (C3 and C4) (Table 2).
The Odonata species richness and abundance varied between study years and habitat types (Figure 2 and Figure 3). Very low Odonata occurrence was observed at reservoirs, tailrace canals, and drainage ditches, where Odonata were detected in low abundances and only in some years, especially in the second half of the study period (after 2017) (Figure 2 and Figure 3). On the other hand, at the old river channels, Odonata were detected during almost the entire study period. The exception was the old river channel next to the Dubrava HPP (D6), where no Odonata were detected between 2013 and 2017 (Figure 2).
In the cluster analysis, Odonata assemblages were mainly not grouped by habitat type—only old river channels partially clustered together (Figure 4).

3.3. Odonata and Environmental Variables

The results of the ordination of species and environmental data from the RDA are shown in an F1 × F2 ordination plot (Figure 5). The distribution of Odonata assemblages was highly related to axes 1 and 2; the eigenvalues were 0.19 and 0.04, respectively, yielding a species–environment correlation of 22.5%. A Monte Carlo permutation test indicated that the species–environment ordination was significant (first axis: F-ratio = 31.71, p = 0.02; overall: trace = 0.24, p = 0.03), suggesting that taxa were significantly related to the tested set of environmental variables. Axis 1 was related to the orthophosphate (R = 0.368) and ammonium concentrations (R = 0.342) and axis 2 was associated with the nitrate (R = −0.100) and mineral oil concentrations in water (R = −0.097), suggesting that they were the most important parameters in explaining the patterns of Odonata assemblages in the habitats studied. Odonata taxa richness and/or abundance were higher at sites and in years with higher concentrations of ammonium and orthophosphates in water and lower concentrations of nitrates and mineral oils (Figure 6 and Figure 7).

4. Discussion

With a total of 11 Odonata species recorded, representing 16% of the Croatian and 8% of the European Odonata fauna [42,43], our results indicate low dragonfly species richness in anthropogenically impacted and man-made habitats in the Drava River system, confirming the results of previous studies [20,21,44]. Although a rather low species richness was expected, one should consider the limitations of the sampling method used, which was not designed as a model for the study of Odonata. Sampling Odonata along with other benthic macroinvertebrates could result in the under-collection of their nymphs because only a relatively small area of habitat is sampled with a benthos hand net [45,46]. Odonata nymphs are relatively large compared with most other benthic macroinvertebrates, and their abundance in benthos samples tends to be low [20,21,44]. In addition, many authors have already emphasized the need to sample all stages of the Odonata life cycle, i.e., including adults and exuviae, in order to obtain a complete list of species at a particular site [45,47,48,49]. The absence of Odonata at several sites throughout the study period could be explained by the combination of the sampling method used here and the low ecological quality of these habitats. Most likely, these sites were not suitable for Odonata nymph development due to the poor water quality, habitat homogenization, i.e., lack of suitable microhabitats (such as aquatic vegetation), and intense water level or discharge fluctuations [20,21,50,51]. Because Odonata nymphs were not collected throughout the study period, our results confirm the importance of long-term data collection to determine Odonata species richness and population dynamics. Previous long-term studies have already emphasized the disadvantages of short-term studies for capturing the true assemblage composition, pointing out that conclusions from short-term ecological studies can be misleading and reduce the relevance of management activities [52]. In the study by Donlý et al. [53], the authors observed different qualitative (species richness and specialization) and quantitative (abundance and dominance) results in long-term studies and long-term monitoring. To cover 95% of the species richness of Odonata, they needed 16 years. Based on their results, a sampling period of at least 10 years is suggested as the ideal period to identify more than 80% of the species.
The Odonata assemblages consisted mainly of generalists, such as Platycnemis pennipes, Erythromma lindenii, Ischnura elegans, and Calopteryx splendens, species that are relatively tolerant of the environmental conditions in their habitats [26,53,54], similar to previous studies on man-made and anthropogenically impacted habitats [20,21,44]. Nevertheless, there was an obvious difference in species richness and abundance between sites with significant anthropogenic impacts (reservoirs, tailrace canals, and drainage ditches) and sites with low human influences (i.e., near-natural old river channels). As expected, old river channels were more suitable for the development of a greater number of Odonata species, as such habitats were characterized by sections of slow-flowing water, lentic sections, and rich aquatic and riparian vegetation. Such habitat heterogeneity allows for the completion of the life cycle of both riverine species, such as Gomphus vulgatissimus and Onychogomphus forcipatus, and standing water Odonata species, such as Enallagma cyathigerum and Erythroma lindenii [54,55]. In addition, our results suggest that habitats with low anthropogenic influence, such as the old river channels associated with HPP reservoirs, may also provide adequate conditions for some species of conservation concern. There we detected Erythromma najas, a species characteristic of habitats with rich aquatic vegetation and standing or slow-flowing water, such as ponds, lakes, canals, and old oxbows [37,42,54]. The species is considered near-threatened in Croatia due to its low dispersal potential, introduction of invasive fish species, degradation of large lowland river floodplains, and construction of HPP reservoirs [42]. Several previous studies have shown that man-made and anthropogenically impacted habitats can harbor protected Odonata species if they are characterized by high habitat heterogeneity, i.e., have well-developed aquatic and riparian vegetation structures [20,21,44,56,57,58].
Many studies have already shown a close relationship between Odonata and environmental variables that characterize their habitats, such as water velocity, water temperature, pH, oxygen concentration, and substrate composition [20,21,28,50,55,56]. Anthropogenic interventions in freshwater ecosystems alter the morphology, hydrology, and physico-chemical water parameters, which also affect the inhabiting biota, including dragonflies [20,21,44,59,60,61]. Of the measured physico-chemical water parameters in our study area, Odonata assemblages were most significantly affected by the concentrations of ammonium, orthophosphates, nitrates, and mineral oils in water. Sites and years with higher concentrations of ammonium and orthophosphates in the water had higher numbers of Odonata species or their abundance. Some studies have shown that slightly elevated concentrations of ammonium from animal waste disposal and orthophosphates from pesticides in water can promote aquatic macrophyte growth [62,63]. Aquatic macrophytes are used by Odonata nymphs for hiding from predators and lurking for prey, and by adults for resting, perching, and oviposition [25,64,65,66,67]. Higher nutrient concentrations may also have resulted in higher population densities of Odonata prey, such as mayflies, which should be investigated in future analyses. Lower Odonata species richness, abundance, or even their absence was recorded at sites and in years with slightly elevated nitrate and mineral oil concentrations. Nitrates in the Drava River system originate mainly from agricultural activities and wastewater. Sites with slightly elevated nitrate levels, such as reservoirs, were most likely highly unsuitable habitats for Odonata, as the water levels fluctuated greatly, preventing the development of aquatic vegetation. An exception is the very shallow Varaždin HPP reservoir, with its dense Elodea canadensis growth, but also dense populations of predatory fish, which most likely also drastically reduce Odonata populations there. Mineral oils are commonly used in agriculture to prevent pest infestations and plant diseases, and negatively affect many terrestrial and some aquatic insects, including some Odonata [67,68].

5. Conclusions

Lotic freshwater ecosystems are among the most threatened ecosystems worldwide, affected by multiple stressors, such as intensive land use in their catchments, morphological alterations, flow regulation, pollution, and climate change. The current study confirmed the negative effects of such stressors on inhabiting biota; more precisely, species-poor Odonata assemblages were recorded with low population densities due to the effects of morphological modifications and water pollution. On the other hand, the preservation of near-natural sites in the vicinity of anthropogenically impacted and man-made habitats is important for maintaining the local Odonata fauna, as well as for the preservation of species of conservation concern. Our results highlight the importance of long-term data for identifying the occurrence of Odonata species and monitoring their population dynamics, and can be used for planning of management and protection activities of lowland rivers and their assemblages according to the requirements of the European Water Framework Directive.

Author Contributions

Conceptualization, M.V. and Z.M.; methodology, M.V. and Z.M.; validation, M.V. and Z.M.; formal analysis, M.V.; investigation, M.V. and Z.M.; resources, Z.M.; data curation, M.V. and Z.M.; writing—original draft preparation, M.V.; writing—review and editing, M.V. and Z.M.; visualization, M.V.; supervision, M.V. and Z.M.; project administration, Z.M.; funding acquisition, Z.M. Both authors have read and agreed to the published version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available from the corresponding author upon request.

Acknowledgments

This study was conducted as part of the long-term monitoring of the Drava River, “Biological assessments of surface waters at the Varaždin HPP, the Čakovec HPP and the Dubrava HPP”. We thank Perica Mustafić and colleagues from the University of Zagreb, Faculty of Science, for their help in sampling, sorting of collected samples, and measurement of environmental variables. We thank two anonymous reviewers and the academic editor for their useful comments and suggestions, which significantly improved this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cazzolla Gatti, R. Freshwater biodiversity: A review of local and global threats. Int. J. Environ. Stud. 2016, 73, 887–904. [Google Scholar] [CrossRef]
  2. Ivković, M.; Plant, A. Aquatic insects in the Dinarides: Identifying hotspots of endemism and species richness shaped by geological and hydrological history using Empididae (Diptera). Insect Conserv. Divers. 2015, 8, 302–312. [Google Scholar] [CrossRef]
  3. Previšić, A.; Graf, W.; Vitecek, S.; Kučinić, M.; Bálint, M.; Keresztes, L.; Pauls, S.U.; Waringer, J. Cryptic diversity of caddisflies in the Balkans: The curious case of Ecclisopteryx species (Trichoptera: Limnephilidae). Arthropod. Syst. Phylogeny 2014, 72, 309–329. [Google Scholar] [PubMed]
  4. Aylward, B.; Bandyopadhyay, J.; Belausteguigotia, J.-C. Freshwater Ecosystem Services. Ecosyst. Hum. Well-Being Policy Responses 2005, 3, 213–256. [Google Scholar]
  5. Naiman, R.J.; Magnuson, J.J.; McKnight, D.M.; Stanford, J.A.; Karr, J.R. Freshwater ecosystems and management: A national initiative. Science 1995, 270, 584–585. [Google Scholar] [CrossRef]
  6. Sala, O.E.; Chapin, F.S.; Armesto, J.J.; Berlow, E.; Bloomfield, J.; Dirzo, R.; Huber-Sanwald, E.; Huenneke, L.F.; Jackson, R.B.; Kinzig, A.; et al. Global Biodiversity Scenarios for the Year 2100. Science 2000, 287, 1770–1774. [Google Scholar] [CrossRef]
  7. Vörösmarty, C.J.; McIntyre, P.B.; Gessner, M.O.; Dudgeon, D.; Prusevich, A.; Green, P.; Glidden, S.; Bunn, S.E.; Sullivan, C.A.; Liermann, C.R.; et al. Global threats to human water security and river biodiversity. Nature 2010, 467, 555–561. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Däll, P.; Zhang, J. Impact of climate change on freshwater ecosystems: A global-scale analysis of ecologically relevant river flow alterations. Hydrol. Earth Syst. Sci. 2010, 14, 783–799. [Google Scholar] [CrossRef] [Green Version]
  9. Pham, H.V.; Torresan, S.; Critto, A.; Marcomini, A. Alteration of freshwater ecosystem services under global change—A review focusing on the Po River basin (Italy) and the Red River basin (Vietnam). Sci. Total Environ. 2019, 652, 1347–1365. [Google Scholar] [CrossRef] [PubMed]
  10. Žganec, K.; Lajtner, J.; Ćuk, R.; Crnčan, P.; Pušić, I.; Atanacković, A.; Kralj, T.; Valić, D.; Jelić, M.; Maguire, I. Alien macroinvertebrates in Croatian freshwaters. Aquat. Invasions 2020, 15, 593–615. [Google Scholar] [CrossRef]
  11. Carpenter, S.R.; Stanley, E.H.; Vander Zanden, M.J. State of the world’s freshwater ecosystems: Physical, chemical, and biological changes. Annu. Rev. Environ. Resour. 2011, 36, 75–99. [Google Scholar] [CrossRef]
  12. Giang, P.Q.; Toshiki, K.; Kunikane, S.; Sakata, M. Climate change challenges transboundary water resources management: Drawing from the case of Vietnam. In Proceedings of the 3rd International Conference on Chemical, Biological and Environmental Sciences, Kuala Lumpur, Malaysia, 8–9 January 2013. [Google Scholar]
  13. Arribas, P.; Abellán, P.; Velasco, J.; Millán, A.; Sánchez-Fernández, D. Conservation of insects in the face of global climate change. In Global Climate Change and Terrestrial Invertebrates; Johnson, S.N., Jones, T.H., Eds.; Wiley-Blackwell: Chichester, UK, 2017; p. 416. [Google Scholar]
  14. Assandri, G. Anthropogenic-driven transformations of dragonfly (Insecta: Odonata) communities of low elevation mountain wetlands during the last century. Insect Conserv. Divers. 2020, 14, 26–39. [Google Scholar] [CrossRef]
  15. Cerini, F.; Stellati, L.; Luiselli, L.; Vignoli, L. Long-term shifts in the communities of Odonata: Effect of chance or climate change? North-West J. Zool. 2020, 16, 1–6. [Google Scholar]
  16. Dudgeon, D.; Arthington, A.H.; Gessner, M.O.; Kawabata, Z.; Knowler, D.; Lévêque, C.; Naiman, R.J.; Prieur-Richard, A.-H.; Soto, D.; Stiassny, M.L.J.; et al. Freshwater biodiversity: Importance, status, and conservation challenges. Biol. Rev. 2006, 81, 163–182. [Google Scholar] [CrossRef]
  17. Marten, M. Environmental monitoring in Baden-Württemberg with special reference to biocoenotic trend-monitoring of macrozoobenthos in rivers and methodical requirements for evaluation of long-term biocoenotic changes. Aquat. Ecol. 2001, 35, 159–171. [Google Scholar] [CrossRef]
  18. Outhwaite, C.L.; Gregory, R.D.; Chandler, R.E.; Collen, B.; Isaac, N.J.B. Complex long-term biodiversity change among invertebrates, bryophytes and lichens. Nat. Ecol. Evol. 2020, 4, 384–392. [Google Scholar] [CrossRef]
  19. Vilenica, M.; Previšić, A.; Ivković, M.; Popijač, A.; Vučković, I.; Kučinić, M.; Kerovec, M.; Gattolliat, J.-L.; Sartori, M.; Mihaljević, Z. Mayfly (Insecta: Ephemeroptera) assemblages of a regulated perennial Mediterranean river system in the Western Balkans. Biologia 2016, 71, 1038–1048. [Google Scholar] [CrossRef]
  20. Vilenica, M.; Kerovec, M.; Pozojević, I.; Mihaljević, Z. Odonata assemblages in anthropogenically impacted lotic habitats. J. Limnol. 2020, 80, 1–9. [Google Scholar] [CrossRef]
  21. Vilenica, M.; Pozojević, I.; Vučković, N.; Mihaljević, Z. How suitable are man-made water bodies as habitats for Odonata? Knowl. Manag. Aquat. Ecosyst. 2020, 421, 1–10. [Google Scholar] [CrossRef] [Green Version]
  22. Franklin, J.F. Importance and Justification of Long-Term Studies in Ecology. In Long-Term Studies in Ecology; Likens, G.E., Ed.; Springer: New York, NY, USA, 1989; pp. 3–19. [Google Scholar] [CrossRef]
  23. Magurran, A.E.; Baillie, S.R.; Buckland, S.T.; Dick, J.M.; Elston, D.A.; Scott, E.M.; Smith, R.I.; Somerfield, P.J.; Watt, A.D. Long-term datasets in biodiversity research and monitoring: Assessing change in ecological communities through time. Trends Ecol. Evol. 2010, 25, 574–582. [Google Scholar] [CrossRef]
  24. Monzó, J.C.; Verdú, J.R. Effects of restoration and management of Mediterranean traditional water systems on Odonata alpha diversity: A long-term monitoring survey. Biodivers. Conserv. 2022, 31, 227–243. [Google Scholar] [CrossRef]
  25. Corbet, P.; Brooks, S. Dragonflies (Collins New Naturalist Library, Book 106); Collins: London, UK, 2008. [Google Scholar]
  26. Schmidt-Kloiber, A.; Hering, D. www.freshwaterecology.info—An online tool that unifies, standardises and codifies more than 20,000 European freshwater organisms and their ecological preferences. Ecol. Indic. 2015, 53, 271–282. [Google Scholar] [CrossRef]
  27. Kutcher, T.E.; Bried, J.T. Adult Odonata Conservatism as an Indicator of Freshwater Wetland Condition. Ecol. Indic. 2014, 38, 31–39. [Google Scholar] [CrossRef]
  28. Simaika, J.P.; Samways, M.J. An Easy-to-Use Index of Ecological Integrity for Prioritizing Freshwater Sites and for Assessing Habitat Quality. Biodivers. Conserv. 2009, 18, 1171–1185. [Google Scholar] [CrossRef] [Green Version]
  29. Vilenica, M. Ecological traits of dragonfly (Odonata) assemblages along an oligotrophic Dinaric karst hydrosystem. Ann. Limnol. Int. J. Limnol. 2017, 53, 377–389. [Google Scholar] [CrossRef] [Green Version]
  30. Vilenica, M.; Rebrina, F.; Matoničkin Kepčija, R.; Šegota, V.; Rumišek, M.; Ružanović, L.; Brigić, A. Aquatic macrophyte vegetation promotes taxonomic and functional diversity of Odonata assemblages in intermittent karst rivers in the Mediterranean. Diversity 2022, 14, 31. [Google Scholar] [CrossRef]
  31. Maguire, I.; Gottstein-Matočec, S. The distribution pattern of freshwater crayfish in Croatia. Crustaceana 2004, 77, 25–47. [Google Scholar] [CrossRef]
  32. Previšić, A.; Mihaljević, Z.; Kerovec, M. Caddisfly (Insecta: Trichoptera) fauna of altered and man-made habitats in the Drava River, NW Croatia. Nat. Croat. 2007, 16, 181–187. [Google Scholar]
  33. Csiky, J.; Purger, D. Monitoring of plant species along the Drava river and in Baranja (Croatia). In Biodiversity Studies along the Drava River; Purger, J.J., Ed.; University of Pécs: Pécs, Hungary, 2008; p. 328. [Google Scholar]
  34. UNESCO. 2012. Available online: https://en.unesco.org/biosphere/eu-na/mura-drava-danube (accessed on 3 September 2022).
  35. Bonacci, O.; Tadić, Z.; Trninić, D. Effects of dams and reservoirs on the hydrological characteristics of the lower Drava River. Regul. Rivers Res. Manag. 1992, 7, 349–357. [Google Scholar] [CrossRef]
  36. Gerken, B.; Sternberg, K. Die Exuvien Europaïscher Libellen—The Exuviae of European Dragonflies (Insecta, Odonata); Arnika & Eisvogel/Huxaria Druckerei GmbH: Höxter/Jena, Germany, 1999; p. 354. [Google Scholar]
  37. Askew, R.R. The Dragonflies of Europe, 2nd ed.; Harley Books: Essex, UK, 2004; p. 308. [Google Scholar]
  38. Brochard, C.; Groendijk, D.; van der Ploeg, E.; Termaat, T. Fotogids Larvenhuidjes van Libellen; KNNV Publishers: Zeist, The Netherlands, 2012; p. 320. [Google Scholar]
  39. APHA. Standard Methods for the Examination of Water and Wastewater, 18th ed.; American Public Health Association: Washington, DC, USA, 1992. [Google Scholar]
  40. Clarke, K.R.; Gorley, R.N. Primer V6: User Manual/Tutorial; Plymouth: Auckland, New Zealand, 2006. [Google Scholar]
  41. Ter Braak, C.J.F.; Šmilauer, P. Canoco Reference Manual and User’s Guide: Software for Ordination, Version 5.0; Microcomputer Power: Ithaca, NY, USA, 2012. [Google Scholar]
  42. Belančić, A.; Bogdanović, T.; Franković, M.; Ljuština, M.; Mihoković, N.; Vitas, B. Red Data Book of Dragonflies of Croatia; Franković, M., Ed.; Ministry of Culture, State Institute for Nature Protection: Zagreb, Croatia, 2008. [Google Scholar]
  43. Kalkman, V.J.; Boudot, J.-P.; Bernard, R.; Conze, K.-J.; De Knijf, G.; Dyatlova, E.; Ferreira, S.; Jović, M.; Ott, J.; Riservato, E.; et al. European Red List of Dragonflies; Publications Office of the European Union: Luxembourg, 2010. [Google Scholar]
  44. Vidaković Maoduš, I.; Pozojević, I.; Vilenica, M.; Mihaljević, Z. Longitudinal dynamics of Odonata assemblages in an anthropogenically impacted lotic system. Ann. Limnol. Int. J. Lim. 2022, 58, 1–13. [Google Scholar] [CrossRef]
  45. Horning, C.E.; Pollard, J.E. Macroinvertebrate Sampling Techniques for Streams in Semi-Arid Regions: Comparison of the Surber Method and Unit-Effort Traveling Kick Method; Report Number EPA 600-4-78-040; US EPA: Washington, DC, USA, 1978; p. 37.
  46. Martín, R.; Maynou, X. Dragonflies (Insecta: Odonata) as indicators of habitat quality in Mediterranean streams and rivers in the province of Barcelona (Catalonia, Iberian Peninsula). Int. J. Odonatol. 2016, 19, 107–124. [Google Scholar] [CrossRef]
  47. Samways, M.J.; McGeoch, M.A.; New, T.R. Insect Conservation: Handbook of Approaches and Methods; Oxford University Press: Oxford, UK, 2009; p. 432. [Google Scholar]
  48. Raebel, E.M.; Merckx, T.; Riordan, P.; Macdonald, D.W.; Thompson, D.J. The dragonfly delusion: Why it is essential to sample exuviae to avoid biased surveys. J. Insect Conserv. 2010, 14, 523–534. [Google Scholar] [CrossRef]
  49. Golfieri, B.; Hardersen, S.; Maiolini, B.; Surian, N. Odonates as indicators of the ecological integrity of the river corridor: Development and application of the Odonate River Index (ORI) in northern Italy. Ecol. Indic. 2016, 61, 234–247. [Google Scholar] [CrossRef]
  50. Da Silva Monteiro, C., Jr.; Couceiro, S.R.M.; Hamada, N.; Juen, L. Effect of vegetation removal for road building on richness and composition of Odonata communities in Amazonia. Brazil. Int. J. Odonatol. 2013, 16, 135–144. [Google Scholar] [CrossRef]
  51. Pereira-Moura, L.; de Sena, W.S.; Neiss, U.G.; Marques Couceiro, S.R. Environmental integrity as a modeler of the composition of the Odonata community. Environ. Monit. Assess. 2021, 193, 160. [Google Scholar] [CrossRef] [PubMed]
  52. White, E.R. Minimum time required to detect population trends: The need for long-term monitoring programs. BioScience 2019, 69, 40–46. [Google Scholar] [CrossRef] [Green Version]
  53. Dolný, A.; Pyszko, P.; Šigutová, H. Community changes in Odonate monitoring: Why are long-term studies so relevant? Insect Conserv. Divers. 2021, 14, 597–608. [Google Scholar] [CrossRef]
  54. Dijkstra, K.D.B.; Lewington, R. Field Guide to the Dragonflies of Britain and Europe, 1st ed.; British Wildlife Publishing: Gillingham, UK, 2006; p. 320. [Google Scholar]
  55. Dijkstra, K.-D.B.; Wildermuth, H.; Martens, A. Dataset “Odonata”. The Taxa and Autecology Database for Freshwater Organisms, Version 8.0. Available online: www.freshwaterecology.info (accessed on 4 September 2022).
  56. Balázs, A.; Šipoš, J.; Matúšová, Z.; Hamerlík, L.; Novikmec, M.; Svitok, M. Comparison of conservation values among man-made aquatic habitats using Odonata communities in Slovakia. Biologia 2022. [Google Scholar] [CrossRef]
  57. Holtmann, L.; Juchem, M.; Brüggeshemke, J.; Möhlmeyer, A.; Fartmann, T. Stormwater ponds promote dragonfly (Odonata) species richness and density in urban areas. Ecol. Eng. 2018, 118, 1–11. [Google Scholar] [CrossRef]
  58. Šigutová, H.; Pyszko, P.; Valušák, J.; Dolný, A. Highway stormwater ponds as islands of Odonata diversity in an agricultural landscape. Sci. Total Environ. 2022, 837, 155774. [Google Scholar] [CrossRef]
  59. Carvalho, F.G.; Oliveira Roque, F.; Barbosa, L.; Assis Montag, L.F.; Juen, L. Oil palm plantation is not a suitable environment for most forest specialist species of Odonata in Amazonia. Anim. Conserv. 2018, 21, 526–533. [Google Scholar] [CrossRef]
  60. Calvão, L.B.; Juen, L.; De Oliveira Junior, J.M.B.; Batista, J.D.; Júnior, P.D.M. Land use modifies Odonata diversity in streams of the Brazilian Cerrado. J. Insect Conserv. 2018, 22, 675–685. [Google Scholar] [CrossRef]
  61. Oliveira, D.G.; Vargas, R.R.; Saad, A.R.; Arruda, R.O.M.; Dalmas, F.B.; Azevedo, F.D. Land use and its impacts on the water quality of the Cachoeirinha Invernada Watershed, Guarulhos (SP). Rev. Ambient. Agua 2018, 13, 1–17. [Google Scholar] [CrossRef]
  62. Boeykens, S.P.; Piol, M.N.; Legal, L.S.; Saralegui, A.B.; Vázquez, C. Eutrophication decrease: Phosphate adsorption processes in presence of nitrates. J. Environ. Manag. 2017, 203, 888–895. [Google Scholar] [CrossRef]
  63. Clarke, E.; Baldwin, A.H. Responses of wetland plants to ammonia and water level. Ecol. Eng. 2022, 18, 257–264. [Google Scholar] [CrossRef]
  64. Perron, M.A.C.; Pick, F.R. Stormwater ponds as habitat for Odonata in urban areas: The importance of obligate wetland plant species. Biodivers. Conserv. 2020, 29, 913–931. [Google Scholar] [CrossRef]
  65. Purse, B.V.; Thompson, D.J. Oviposition Site Selection by Coenagrion mercuriale (Odonata: Coenagrionidae). Int. J. Odonatol. 2009, 12, 257–273. [Google Scholar] [CrossRef]
  66. Brito, J.S.; Michelan, T.S.; Juen, L. Aquatic macrophytes are important substrates for Libellulidae (Odonata) larvae and adults. Limnology 2021, 22, 139–149. [Google Scholar] [CrossRef]
  67. Hu, W.; Zheng, R.; Feng, X.; Kuang, F.; Chun, J.; Xu, H.; Chen, T.; Lu, J.; Li, W.; Zhang, N. Emergence inhibition, repellent activity and antifeedant responds of mineral oils against Asian citrus psyllid, Diaphorina citri (Hemiptera: Liviidae). Int. J. Pest Manag. 2020, 1–8. [Google Scholar] [CrossRef]
  68. Mulla, M.S.; Darwazeb, H.A.; Ede, L. Evaluation of new pyrethroids against immature mosquitoes and their effects on nontarget organisms. Mosq. News 1982, 42, 583–590. [Google Scholar]
Figure 1. Map of the study area located along the Drava River in Croatia.
Figure 1. Map of the study area located along the Drava River in Croatia.
Water 14 03119 g001
Figure 2. Odonata taxa richness and abundance at the studied reservoirs and old river channels in the Drava River between 2012 and 2020. Sites of the same habitat groups without Odonata records are not shown.
Figure 2. Odonata taxa richness and abundance at the studied reservoirs and old river channels in the Drava River between 2012 and 2020. Sites of the same habitat groups without Odonata records are not shown.
Water 14 03119 g002
Figure 3. Odonata taxa richness and abundance at the studied tailrace canals and drainage ditches in the Drava River between 2012 and 2020. Sites in the same habitat groups without Odonata records are not shown.
Figure 3. Odonata taxa richness and abundance at the studied tailrace canals and drainage ditches in the Drava River between 2012 and 2020. Sites in the same habitat groups without Odonata records are not shown.
Water 14 03119 g003
Figure 4. Cluster analysis of Odonata assemblages based on the Bray–Curtis similarity coefficient (group average linking) and their log(x + 1)-transformed abundances at anthropogenically impacted habitats in the Drava River between 2012 and 2020. Sites without Odonata records were excluded from the analysis.
Figure 4. Cluster analysis of Odonata assemblages based on the Bray–Curtis similarity coefficient (group average linking) and their log(x + 1)-transformed abundances at anthropogenically impacted habitats in the Drava River between 2012 and 2020. Sites without Odonata records were excluded from the analysis.
Water 14 03119 g004
Figure 5. F1 × F2 plane of redundancy analysis for species–environment relationships for Odonata assemblages in anthropogenically impacted habitats in the Drava River between 2012 and 2020. Taxa abbreviations (blue arrow symbols): P pen = Platycnemis pennipes, Cal spl = Calopteryx splendens, Cal vir = Calopteryx virgo, Coe = Coenagrionidae, I ele = Ischnura elegans, E naj = Erythromma najas, E lin = Erythromma lindenii, E cya = Enallagma cyathigerum. G vul = Gomphus vulgatissimus, O for = Onychogomphus forcipatus, O bru = Orthetrum brunneum, and O can = Orthetrum cancellatum. Environmental factors (red arrow symbols): WT = water temperature (°C), pH = pH value, CON = conductivity (µS/cm), DO = dissolved oxygen (%), BOD = biological oxygen demand (O2 mg/L), COD = chemical oxygen demand (O2 mg/L), NH4 = ammonium concentration (N mg/L), NO2 = nitrite concentration (N mg/L), NO3 = nitrate concentration (N mg/L), TN = total nitrogen concentration (N mg/L), TP = total phosphorous concentration (P mg/L), PO4 = orthophosphate concentration (P mg/L), PES = pesticide concentration (µg/L), MIN-O = mineral oil concentration (mg/L), and COLIF = fecal coliform bacteria (nr/1000 mL).
Figure 5. F1 × F2 plane of redundancy analysis for species–environment relationships for Odonata assemblages in anthropogenically impacted habitats in the Drava River between 2012 and 2020. Taxa abbreviations (blue arrow symbols): P pen = Platycnemis pennipes, Cal spl = Calopteryx splendens, Cal vir = Calopteryx virgo, Coe = Coenagrionidae, I ele = Ischnura elegans, E naj = Erythromma najas, E lin = Erythromma lindenii, E cya = Enallagma cyathigerum. G vul = Gomphus vulgatissimus, O for = Onychogomphus forcipatus, O bru = Orthetrum brunneum, and O can = Orthetrum cancellatum. Environmental factors (red arrow symbols): WT = water temperature (°C), pH = pH value, CON = conductivity (µS/cm), DO = dissolved oxygen (%), BOD = biological oxygen demand (O2 mg/L), COD = chemical oxygen demand (O2 mg/L), NH4 = ammonium concentration (N mg/L), NO2 = nitrite concentration (N mg/L), NO3 = nitrate concentration (N mg/L), TN = total nitrogen concentration (N mg/L), TP = total phosphorous concentration (P mg/L), PO4 = orthophosphate concentration (P mg/L), PES = pesticide concentration (µg/L), MIN-O = mineral oil concentration (mg/L), and COLIF = fecal coliform bacteria (nr/1000 mL).
Water 14 03119 g005
Figure 6. Odonata taxa richness in relation to the (a) ammonium, (b) orthophosphate, (c) nitrate, and (d) mineral oil concentrations in water in anthropogenically impacted habitats in the Drava River between 2012 and 2020.
Figure 6. Odonata taxa richness in relation to the (a) ammonium, (b) orthophosphate, (c) nitrate, and (d) mineral oil concentrations in water in anthropogenically impacted habitats in the Drava River between 2012 and 2020.
Water 14 03119 g006
Figure 7. Odonata abundance in relation to the (a) ammonium, (b) orthophosphate, (c) nitrate, and (d) mineral oil concentrations in water in anthropogenically impacted habitats in the Drava River between 2012 and 2020.
Figure 7. Odonata abundance in relation to the (a) ammonium, (b) orthophosphate, (c) nitrate, and (d) mineral oil concentrations in water in anthropogenically impacted habitats in the Drava River between 2012 and 2020.
Water 14 03119 g007
Table 1. Mean values (and SD, standard deviation) of environmental parameters in anthropogenically impacted habitats in the Drava River measured between 2012 and 2020. Higher values of environmental parameters are in bold.
Table 1. Mean values (and SD, standard deviation) of environmental parameters in anthropogenically impacted habitats in the Drava River measured between 2012 and 2020. Higher values of environmental parameters are in bold.
Habitats/Environmental ParametersReservoirsDrainage DitchesTailrace CanalsOld River Channels
MeanSDMeanSDMeanSDMeanSD
Water temperature (°C)14.982.0913.400.9814.831.2415.311.55
pH8.120.167.980.178.120.158.090.16
Conductivity (µS/cm)28325339192842130729
Dissolved oxygen concentration (O2 mg/L)9.840.626.470.719.640.499.390.90
Oxygen saturation (%)101.565.9364.907.2797.593.8495.919.49
Biological oxygen demand (O2 mg/L)1.590.861.090.911.420.841.480.90
Chemical oxygen demand (O2 mg/L)6.335.905.875.305.605.495.505.27
Ammonium concentration (N mg/L)0.0280.0210.0460.0570.0200.0080.0850.120
Nitrite concentration (N mg/L)0.0380.1420.0460.1470.0450.1700.0440.138
Nitrate concentration (N mg/L)0.7930.1920.5950.1780.7690.1750.7530.213
Total nitrogen concentration (N mg/L)1.2170.4220.8290.4471.1050.2731.1140.352
Total phosphorous concentration (P mg/L)0.0620.0500.0440.0480.0460.0320.0540.040
Orthophosphate concentration (P mg/L)0.0170.0070.0200.0160.0150.0030.0250.024
Phenol concentration (mg/L)0.00610.00450.00530.00300.00730.00620.00820.0104
Pesticide concentration (µg/L)0.00330.00240.00310.00240.00310.00250.00320.0026
Mineral oil concentration (mg/L)0.0160.0190.0170.0230.0200.0290.0160.024
Faecal coliform bacteria (nr/1000 mL)4610592220,728116,6042975445914,44937,149
Table 2. Odonata taxa recorded in anthropogenically impacted habitats in the Drava River between 2012 and 2020. Abundance is shown as the number of individuals per m2.
Table 2. Odonata taxa recorded in anthropogenically impacted habitats in the Drava River between 2012 and 2020. Abundance is shown as the number of individuals per m2.
Habitat TypeReservoirsDrainage DitchesTailrace
Canals
Old River
Channels
SiteC1C1AC2D1D2V1V2C3C4D3D4C5D5V3C6D6V4
Platycnemis pennipes (Pallas, 1771)6.71.7 3.3 3.36.7 30.015.61.7
Calopteryx splendens (Harris, 1782) 1.7
Calopteryx virgo (Linnaeus, 1758)1.7 3.3
Coenagrionidae non det. 6.7 1.7
Ischnura elegans Vander Linden, 1820 2.81.7
Erythromma najas Hansemann, 1823 1.7
Erythromma lindenii (Selys, 1840) 11.7
Enallagma cyathigerum (Charpentier, 1840) 11.7
Gomphus vulgatissimus (Linnaeus, 1758) 1.7
Onychogomphus forcipatus (Linnaeus, 1758)5.0 1.7 5.06.710.015.0
Orthetrum brunneum (Fonscolombe, 1837) 1.7
Orthetrum cancellatum (Linnaeus, 1758) 6.7
Taxa richness31030000011112455
Abundance13.41.70.011.70.00.00.00.00.01.71.73.36.78.360.136.821.8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vilenica, M.; Mihaljević, Z. Odonata Assemblages in Anthropogenically Impacted Habitats in the Drava River—A Long-Term Study. Water 2022, 14, 3119. https://doi.org/10.3390/w14193119

AMA Style

Vilenica M, Mihaljević Z. Odonata Assemblages in Anthropogenically Impacted Habitats in the Drava River—A Long-Term Study. Water. 2022; 14(19):3119. https://doi.org/10.3390/w14193119

Chicago/Turabian Style

Vilenica, Marina, and Zlatko Mihaljević. 2022. "Odonata Assemblages in Anthropogenically Impacted Habitats in the Drava River—A Long-Term Study" Water 14, no. 19: 3119. https://doi.org/10.3390/w14193119

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop