Research Article |
Corresponding author: Dušanka Cvijanović ( dusanka.cvijanovic@dbe.uns.ac.rs ) Academic editor: Mathias Scholz
© 2025 Dušanka Cvijanović, Maja Novković, Djuradj Milošević, Milica Stojković Piperac, Laszlo Galambos, Dubravka Čerba, Olivera Stamenković, Bojan Damnjanović, Minučer Mesaroš, Dragoslav Pavić, Vladica Simić, Ivana Trbojević, Ana Anđelković, Nusret Drešković, Barbara Stammel, Bernd Cyffka, Snežana Radulović.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Cvijanović D, Novković M, Milošević D, Stojković Piperac M, Galambos L, Čerba D, Stamenković O, Damnjanović B, Mesaroš M, Pavić D, Simić V, Trbojević I, Anđelković A, Drešković N, Stammel B, Cyffka B, Radulović S (2025) Conservation and ecological screening of small water bodies in temperate riverine wetlands using UAV Photogrammetry (Middle Danube). Nature Conservation 58: 61-82. https://doi.org/10.3897/natureconservation.58.116663
|
Aquatic ecosystems in riverine wetlands are important refuges and nurseries for freshwater biota. Given the significant global loss and degradation of wetlands, regular conservation assessments of these habitats, even in not easily accessible regions, are crucial for implementing effective management. Thus, developing cost-effective approaches for rapid ecological and conservation screening of water bodies in floodplains, such as the Danube, is a priority. One potential solution is the use of UAV-based (Unmanned Aerial Vehicle) ecological indicators to complement existing monitoring frameworks. This paper aims to explore whether UAV-based macrophyte data can provide a more precise indication of the trophic state and conservation indices (assessed through fish and macroinvertebrate communities) of temperate wetland lentic ecosystems, compared to traditional field surveys. The fieldwork was conducted during the summer months of 2019 at 23 sampling sites within eight lentic water bodies located in three wetland areas along the Middle Danube in Serbia. Data on aquatic vegetation, fish, and macroinvertebrate communities, and samples for water quality analysis were collected simultaneously. UAV images were acquired using an RGB camera. Orthomosaics were processed using supervised object-based image (OBIA) classification to obtain a single vector layer with macrophyte functional groups and taxa. Macrophyte cover metrics obtained during the fieldwork and UAV data processing were correlated against water quality parameters and conservation indices calculated for fish and macroinvertebrate assemblages. The study demonstrated that UAV photogrammetry can provide relatively precise measurements of macrophyte cover characteristics compared to traditional plot-based monitoring methods, making it effective for assessing aquatic ecosystems. The analysis revealed that sites with high values of fish and macroinvertebrate conservation indices, optimal oxygen conditions, and mesotrophic states were associated with UAV orthomosaic polygons showing relatively high macrophyte functional diversity and a presence of floating-rooted species. Conversely, sites experiencing eutrophication and a poor oxygen regime with species-poor fish assemblages correlated positively with a higher cover of amphibian and free-floating vegetation, as well as filamentous algae. In conclusion, UAV photogrammetry offers a cost-effective method to monitor aquatic habitats along large river floodplains, including those that are not easily accessible.
Aquatic vegetation, fish, macroinvertebrates, ponds, riverine wetlands, UAV indicators
Aquatic ecosystems in riverine wetlands represent important refuges and nurseries for freshwater biota (
Aquatic vegetation is widely used in the conservation assessment of freshwater ecosystems, serving as surrogates for diversity indices of macroinvertebrate and fish communities (
The presence and characteristics of aquatic vegetation, which are considered robust ecological or conservation indicators in dynamic wetland landscapes (
The aim of this paper is to explore if macrophyte cover data derived from UAV images can provide a more precise indication of the trophic state and conservation index values (assessed through fish and macroinvertebrate communities) of the Danube wetland lentic ecosystems in Serbia, compared to the traditional field survey. In order to achieve the primary objective of the study, the following tasks were established: 1) to determine water quality (dissolved oxygen, orthophosphate, and total organic carbon) and conservation indices of lentic water bodies based on fish and macroinvertebrate assemblages; 2) to compare the sensitivity of macrophyte metrics derived by UAV monitoring (UAVM) and by field monitoring (FM) for the conservation assessment; 3) to identify macrophyte metrics obtained by UAV monitoring, which significantly indicate sites having high conservation index values.
The large floodplains along the Middle Danube are recognized as sites of high conservation value and importance at the national and international levels (
Field work was carried out during the summer months of 2019 on 23 sampling sites within the 8 lentic water bodies (Fig.
Vegetation data were collected within the circle polygons of 2.5 m radius using the species relative abundance DAFOR scale following standard method EN 15460:2007 (
Macroinvertebrate communities were collected in transects - one per each sampling site. At each transect, three benthic samples were collected with a 15 × 15 cm Ekman grab. Transects were distributed to cover all mesohabitats, starting from the shoreline, towards the increasingly deeper water. The benthic samples were preserved in 70% ethyl alcohol and individuals were sorted in the laboratory. All macroinvertebrates were identified to the lowest possible taxonomic level (mostly species or genus) using the relevant taxonomic keys (
The fish were sampled along transects from a boat using a DC Aquatech IG 1300 electro-fisher (2.6 kW, 80–470 V). For each selected transect, the constant catch-per-unit-effort (CPUE) of time (10 min) was provided. Each fish was identified at the species level.
UAV images were acquired by Phantom 4 FC330 (12.5MP) RGB camera on summer sunny days in August 2019, between 7:10 a.m. and 12:17 p.m. to correspond with an in-situ field survey (
Conservation indices, relevant for the conservation management of fluvial lentic ecosystems were obtained for fish and macroinvertebrate assemblages for each sampling site: Shannon diversity index (SD) (
Conservation scores C and Csp (
For each sampling site, the cover of a single species was summarized to the following macrophyte metrics, considered as explanatory variables in the further analyses: the total macrophyte cover, the total cover of emergent macrophytes, the total cover of rooted floating-leaved macrophytes; total cover of free-floating macrophytes; and the total cover of submerged macrophytes. The floating-leaved (rooted) macrophyte group was also differentiated into the cover value of waterlily species, the cover value of Nymphoides peltata, and the cover value of Trapa natans.
For each water body, UAV-based geotagged images were block adjusted and stitched into individual georeferenced orthomosaics using default settings of the Adjust and Orthomosaic wizard tools within the ArcGIS Pro 2.6.0 software (Perform Camera Calibration checked, Blunder Point Threshold 5, Image Resolution Factor 8× Source Resolution) (
Image segmentation was performed using LargeScaleMeanShift algorithm (Spatial Radius: 30; Range Radius: 10; Minimum Segment Size: 50 px and Tile Size: 1024 × 1024 px). During the segmentation process, orthomosaic features were partitioned into discrete entities – segments based on the similarity of their spectral characteristics and spatial distribution. In order to increase classification accuracy, a set of classification attributes (RGB spectral and texture indices, Suppl. material
For each orthomosaic, a training data set was created by selecting 50 representative reference polygons representing a specific macrophyte image feature class (macrophyte functional groups or stands of particular macrophyte taxa; for further details please see Suppl. material
A validation data set for each orthomosaic included independent and unbiased polygons selected using the Random points tool in QGIS. A different number of random polygons were selected depending on the waterbody size (<1.5 ha – 100; 1.5–2.5 ha – 200; 2.5–3.5 ha – 300 and > 3.5 h – 400 points) (for further details please see Suppl. material
Object-based classification of orthomosaic segments was further performed using a Random Forest (RF) classifier (TrainVectorClassifier and VectorClassifier tools of the Orfeo Toolbox (OTB) 7.2.0), and the following parameters: Maximum depth of trees: 10; Minimum number of samples in each node: 7; Maximum number of trees in forest: 225; OBBerror: 0.01. TrainVectorClassifier tool performs training of the RF algorithm using training and validation dataset, while VectorClassifier tool performs classification of the orthomosaic segments using model file obtained in the previous step.
After the initial classification, orthomaps were visually evaluated. In the orthomosaic areas which were poorly classified additional training segments were assigned to the misclassified image feature categories and added to the training data set. The training process was repeated. This allowed the lowest size of the training data set to be considered in the analysis and to target and address challenging areas of the orthomosaics. This allowed the lowest size of the training data set to be considered in the analysis, while targeting the classification challenging areas of the orthomosaics and image feature categories. As a result of classification, each orthomosaic segment was assigned to a specific macrophyte functional group or macrophyte taxa.
Two approaches of accuracy analysis were applied to macrophyte metrics, Per-Pixel and Per-Polygon. Per-Polygon analysis was performed with the TrainVectorClassifier tool based on the Kappa index. Per-Pixel accuracy was estimated using the Accuracy tool from Semi-Automatic classification plugin, which includes Kappa-hat index, Standard error, Standard error area, Users accuracy, Producers accuracy and Kappa-hat index.
Therefore, the result of the OBIA classification was a single vector layer with different macrophyte image feature classes (macrophyte functional groups or taxa) (for further details please see Suppl. material
An orthomosaic of the Doktor Pumpa fluvial lake - the wetland area located near Bačko Novo Selo (a) with classified main macrophyte metrics; (b) enlarged fraction of an RGB orthomosaic (c) and image classification showing location of field sampling point and also GIS approach for extraction of macrophyte metric areas within the polygons of different radii (d).
Macrophyte metrics obtained during fieldwork and UAV data processing were further correlated against the conservation metrics and water quality (dissolved oxygen content, total organic carbon and orthophosphates) using the non-parametric Spearman’s rank in STATISTICA 14 software (
In total, 43 macrophyte taxa were recorded in the study area, forming vegetation stands of free-floating duckweeds, occasionally submerged anchored ceratophyllids, and rooted aquatic vegetation (for further details please see Suppl. material
Dissolved oxygen concentrations varied among water bodies and within individual sampling points in a single water body—from low oxygen levels (1.37 mg/L) to good water quality (>9 mg/L), indicating conditions ranging from eutrophic to oligo-mesotrophic (
The relatively lower OBIA classification accuracy was observed for the cover of amphibian vegetation (Producer’s acc. 49.35%, Users acc. 30.43%, Kappa-hat index 0.29) and Trapa natans species (Producer’s acc. 93.45%, Users acc. 44.22%, Kappa-hat index 0.42), compared to other macrophyte metrics (maximal Producer’s acc. 87–100%, maximal Users acc. 90–100%, maximal Kappa-hat index 0.89–1) (Suppl. material
Non-parametric Spearman’s rank values (P < 0.05) for significant correlations of macrophyte cover classes obtained during the fieldwork against the conservation indices for fish and macroinvertebrate assemblages and water quality parameters.
Fish Species richness | Fish Csp value | Fish C value | Fish Shannon-Wiener index | Macroinvertebrate Species richness | Macroinvertebrate Csp | Macroinvertebrate Shannon-Wiener index | Orthophosphates | Dissolved oxygen | |
---|---|---|---|---|---|---|---|---|---|
Total cover of free-floating macrophytes | -0.50 | ||||||||
Total cover of floating rooted macrophytes | 0.50 | 0.52 | -0.44 | ||||||
Total cover of emerged macrophytes | -0.56 | -0.54 | -0.52 | -0.50 | 0.46 | ||||
Total cover of submerged macrophytes | 0.45 | 0.45 | 0.58 | -0.58 | -0.50 | 0.51 | |||
Cover of Trapa natans | -0.52 | ||||||||
Cover of Nymphoides peltata | 0.65 | 0.43 | 0.50 | -0.55 |
Non-parametric Spearman’s rank values for significant correlation (P < 0.05) of macrophyte metrics obtained using the UAV data processing against the conservation indices for fish and macroinvertebrate assemblages.
Radius (m) | Fish Species richness | Fish Csp | Fish C | Fish Shannon-Wiener index | Macroinvertebrate Species richness | Macroinvertebrate Csp | Macroinvertebrates C | Macroinvertebrate Shannon-Wiener index | Total Organic Carbon | Orthophosphates | Dissolved oxygen | |
---|---|---|---|---|---|---|---|---|---|---|---|---|
Total cover of aquatic vegetation | 2.5 | 0.43 | 0.42 | |||||||||
5 | 0.43 | |||||||||||
30 | 0.48 | -0.52 | ||||||||||
2.5 | 0.50 | |||||||||||
Total cover of free-floating macrophytes | 5 | 0.53 | 0.53 | -0.44 | ||||||||
10 | 0.48 | 0.48 | -0.47 | |||||||||
15 | 0.48 | 0.48 | ||||||||||
20 | 0.49 | -0.44 | ||||||||||
30 | -0.44 | |||||||||||
Total cover of floating rooted vegetation | 2.5 | 0.56 | 0.42 | 0.65 | 0.50 | 0.42 | ||||||
5 | 0.47 | 0.53 | ||||||||||
10 | 0.48 | |||||||||||
15 | 0.48 | |||||||||||
20 | 0.49 | |||||||||||
30 | -0.45 | |||||||||||
Total cover of submerged vegetation | 5 | -0.44 | ||||||||||
15 | 0.48 | |||||||||||
30 | 0.64 | |||||||||||
Total cover of amphibian vegetation | 15 | -0.49 | -0.49 | -0.49 | -0.55 | -0.48 | ||||||
20 | -0.49 | -0.49 | -0.49 | -0.55 | -0.48 | |||||||
30 | -0.49 | -0.49 | -0.49 | -0.55 | -0.48 | |||||||
Total cover of emergent vegetation | -0.52 | 0.46 | ||||||||||
Filamentous algae | 2.5 | 0.69 | 0.45 | |||||||||
5 | 0.56 | 0.51 | ||||||||||
10 | 0.53 | 0.51 | ||||||||||
15 | 0.53 | 0.51 | ||||||||||
20 | 0.53 | 0.51 | ||||||||||
30 | 0.70 | 0.49 | 0.46 | |||||||||
Trapa natans | 15 | 0.47 | ||||||||||
20 | 0.50 | |||||||||||
Nymphoides peltata | 2.5 | -0.68 | 0.56 | |||||||||
5 | -0.62 | 0.65 | ||||||||||
10 | 0.49 | -0.65 | 0.69 | |||||||||
15 | 0.46 | -0.49 | -0.65 | 0.68 | ||||||||
20 | 0.46 | -0.49 | -0.67 | 0.69 | ||||||||
30 | 0.47 | -0.51 | -0.67 | 0.70 | ||||||||
Nymphaea alba and Nuphar luteum | 2.5 | 0.54 | 0.42 | |||||||||
5 | 0.48 | |||||||||||
10 | 0.60 | |||||||||||
Nymphaea alba and Nuphar luteum | 15 | 0.49 | 0.57 | |||||||||
20 | 0.51 | 0.53 | ||||||||||
30 | 0.49 | -0.44 | ||||||||||
Number of macrophyte communities / dominant macrophyte species | 5 | 0.64 | ||||||||||
10 | 0.70 | |||||||||||
15 | 0.72 | |||||||||||
20 | 0.72 | |||||||||||
30 | 0.46 | 0.46 | 0.56 | 0.50 | 0.43 |
When compared with FM equivalents, UAVM metrics (i.e. free-floating macrophytes, floating-rooted macrophytes, amphibian vegetation and Nymphoides peltata) showed more significant relationships with conservation metrics and water quality (Tables
In summary, fish conservation metrics were positively correlated with floating vegetation types and species in both datasets, in polygons of 2.5 m radius (Tables
The highest correlation coefficients were obtained between water quality parameters and both types of macrophyte metrics, compared to conservation metrics (Table
Our study has shown that total organic carbon and C value (calculated for macroinvertebrate assemblages) are captured exclusively by UAVM metrics. These results were obtained for 5–30 m radius polygons, which are larger than the traditional vegetation plots (16–20 m2) (
Another important finding of this study is that UAVM macrophyte communities serve as effective indicators for fish, macroinvertebrates, and water quality variables. The number of UAVM macrophyte communities was expressed in our study as a number of dominant macrophyte species, aligning with the formal definitions of Level 4 of the European Nature Information System (EUNIS) aquatic habitat classification (
Overall, in our study UAVM metrics demonstrated stronger relationships with conservation metrics and water quality attributes compared to FM metrics. This was expected as UAV photogrammetry allows for precise measurements of aquatic vegetation stands, especially floating ones (
For floating rooted and submerged macrophytes, UAV-based metrics were found to be reliable indicators of ecosystem conservation indices and moderate trophic states. Meanwhile, for the free-floating macrophytes, UAV metrics were effective predictors of eutrophic conditions. As expected, submerged vegetation was shown to be a better indicator in the FCM data set compared to the UAVM data set. This is likely because traditional field monitoring provides a more comprehensive view and inspection of the entire submerged layer than UAV imaging. Submerged vegetation can be detected with high accuracy up to 1.2 m depth in conditions of relatively high turbidity (Secchi depth > 2 m) using UAV RGB imagery (
Although the scoring of the water body conservation indices is usually performed at the ecosystem level (
While the methodologies deployed in this study have shown significant potential, there are noteworthy limitations associated with the UAV photogrammetry and field sampling protocols employed. To enhance the georeferencing accuracy of orthomosaics, it would be beneficial to use Real-Time Kinematic (RTK) correction or integrate ground control points data, collected during field operations, as a corrective measure in the orthomosaic generation process. Nevertheless, non-RTK UAVs equipped with consumer-grade GPS systems produce orthoimages with horizontal accuracy suitable for our study design (
Furthermore, this study was designed to align with the practical constraints often faced by wetland stakeholders, such as limited resources for monitoring and restricted temporal windows for fieldwork. Consequently, a single field sampling and UAV imaging session was conducted immediately following the Danube flooding event, during the peak of the aquatic vegetation season. This timing aligns with typical periods used for the conservation assessment of water bodies in the Middle Danube wetlands, as noted in previous research (
This study demonstrates the importance of relatively precise measuring of macrophyte cover metrics using UAV photogrammetry compared to traditional plot-based monitoring methods in aquatic monitoring. The cost-effective conservation and ecological screening of aquatic habitats along the great river floodplains can be performed using UAV photogrammetry of macrophyte vegetation. A specific combination of macrophyte functional groups and taxa within the complex wetland mosaics showed good surrogates for water trophic state and conservation indices derived for fish and macroinvertebrate assemblages. In summary, sites of a potentially high conservation indices and mesotrophic conditions can be tracked by the presence of floating rooted species and high macrophyte functional diversity. On the other hand, sites subjected to eutrophication and low dissolved oxygen concentration, with species-poor fish assemblages can be detected based on the cover of amphibian and free-floating vegetation and filamentous algae. While this study was conducted on a temperate large river floodplain, the developed methodological approach should be easily upscaled to other catchments and regions. This is especially important for hard-to-reach wetland biodiversity hot spots, including war-affected and mined areas (
The authors have declared that no competing interests exist.
No ethical statement was reported.
This study was funded by the European Union Horizon EU project no. 101112736 –Restoration of wetland complexes as life supporting systems in the Danube Basin (Restore4Life); by the project of the Provincial Secretariat for Higher Education and Scientific Research, Autonomous Province of Vojvodina, Republic of Serbia, no 142-451-2095/2022-01 -"Trophic status assessment of the Danube floodplain using UAV photogrammetry"; by the Rufford grant, No 28388-1, "Toward Cost-Effective UAV-Assisted Multimetric System for Detection of Freshwater Patches of High Conservation Value within the Danube Floodplain in Serbia"; by the Ministry of Science, Technological Development and Innovation of the Republic of Serbia (Grants No. 451-03-66/2024-03/ 200125 & 451-03-65/2024-03/200125, 451-03-66/2024-03/200010 and 451-03-66/2024-03/200124), and funded by the European Union project no.101131220 -“Engaging approaches and services for meaningful climate actions (ClearClimate), under the Marie Skłodowska-Curie Actions (MSCA).
D. Cvijanović: Conceptualization, Formal analysis, Investigation, Software, Methodology, Writing - original draft, Funding acquisition; M. Novković: Formal analysis, Investigation, Software, Writing - review & editing, Funding acquisition; Dj. Milošević: Investigation, Methodology, Writing - review & editing, Funding acquisition; M. Stojković Piperac: Investigation, Methodology, Writing - review & editing; L. Galambos: Investigation, Methodology, Writing - review & editing; D. Čerba: Investigation, Methodology, Writing - review & editing, Funding acquisition; O. Stamenković, B. Damnjanović, M. Mesaroš, D. Pavić, V. Simić, I. Trbojević: Investigation, Writing - review & editing; A. Anđelković, N. Drešković: Writing - review & editing; B. Stammel, B. Cyffka, S. Radulović: Writing - review & editing.
Dušanka Cvijanović https://orcid.org/0000-0001-8732-1269
Maja Novković https://orcid.org/0000-0002-7646-3998
Dubravka Čerba https://orcid.org/0000-0003-2563-8695
Olivera Stamenković https://orcid.org/0000-0001-5438-8870
Minučer Mesaroš https://orcid.org/0000-0003-2505-5633
Dragoslav Pavić https://orcid.org/0000-0002-7113-0887
Ivana Trbojević https://orcid.org/0000-0002-6715-8422
Ana Anđelković https://orcid.org/0000-0001-6616-1710
Barbara Stammel https://orcid.org/0000-0003-3208-4571
All of the data that support the findings of this study are available in the main text or Supplementary Information.
Average parameters per pond
Data type: xlsx
Explanation note: Object based classification parameters for macrophyte cover variables.
RGB-based and texture indices
Data type: docx
Macrophyte OBIA parameters
Data type: xlsx
Explanation note: Object based classification parameters for macrophyte cover variables.
Physical and chemical parameters
Data type: xlsx
Explanation note: Physical and chemical parameters and macrophyte metrics.