Next Article in Journal
Doping Ag in ZnO Nanorods to Improve the Performance of Related Enzymatic Glucose Sensors
Next Article in Special Issue
Adaptation of Dubins Paths for UAV Ground Obstacle Avoidance When Using a Low Cost On-Board GNSS Sensor
Previous Article in Journal
Energy Neutral Wireless Bolt for Safety Critical Fastening
Previous Article in Special Issue
Secure Utilization of Beacons and UAVs in Emergency Response Systems for Building Fire Hazard
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Towards a Transferable UAV-Based Framework for River Hydromorphological Characterization

by
Mónica Rivas Casado
1,*,
Rocío Ballesteros González
2,
José Fernando Ortega
2,
Paul Leinster
1 and
Ros Wright
3
1
School of Water, Energy and Environment, Cranfield University, Cranfield, Bedfordshire MK430AL, UK
2
Regional Centre of Water Research, Universidad de Castilla-La Mancha, Carretera de las Peñas km 3.2, 02071 Albacete, Spain
3
National Fisheries Services, Environment Agency, Threshelfords Business Park, Inworth Road, Feering, Essex CO61UD, UK
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(10), 2210; https://doi.org/10.3390/s17102210
Submission received: 23 June 2017 / Revised: 18 September 2017 / Accepted: 21 September 2017 / Published: 26 September 2017
(This article belongs to the Special Issue UAV or Drones for Remote Sensing Applications)

Abstract

:
The multiple protocols that have been developed to characterize river hydromorphology, partly in response to legislative drivers such as the European Union Water Framework Directive (EU WFD), make the comparison of results obtained in different countries challenging. Recent studies have analyzed the comparability of existing methods, with remote sensing based approaches being proposed as a potential means of harmonizing hydromorphological characterization protocols. However, the resolution achieved by remote sensing products may not be sufficient to assess some of the key hydromorphological features that are required to allow an accurate characterization. Methodologies based on high resolution aerial photography taken from Unmanned Aerial Vehicles (UAVs) have been proposed by several authors as potential approaches to overcome these limitations. Here, we explore the applicability of an existing UAV based framework for hydromorphological characterization to three different fluvial settings representing some of the distinct ecoregions defined by the WFD geographical intercalibration groups (GIGs). The framework is based on the automated recognition of hydromorphological features via tested and validated Artificial Neural Networks (ANNs). Results show that the framework is transferable to the Central-Baltic and Mediterranean GIGs with accuracies in feature identification above 70%. Accuracies of 50% are achieved when the framework is implemented in the Very Large Rivers GIG. The framework successfully identified vegetation, deep water, shallow water, riffles, side bars and shadows for the majority of the reaches. However, further algorithm development is required to ensure a wider range of features (e.g., chutes, structures and erosion) are accurately identified. This study also highlights the need to develop an objective and fit for purpose hydromorphological characterization framework to be adopted within all EU member states to facilitate comparison of results.

1. Introduction

A recognition of the importance of protecting and improving the water environment and the increasing pressures on water resources resulting from population growth, patterns of use and climate change has seen the introduction of major water protection laws in many countries. This includes the development of regulatory frameworks based on an assessment of the ecological quality of freshwater systems such as the European Union Water Framework Directive (EU WFD) [1], which aims to achieve good ecological status or potential of inland and coastal waters.
Each member state has adopted a specific methodology within the overall WFD to classify the ecological status of a water body based on an assessment of the biological, chemical, physico-chemical and supporting elements. The measures derived for individual parameters are based on significantly different methods and are therefore subject to multiple sources of difference. This variation in approach between member states means that it can be difficult to compare the results obtained in the various countries. In the particular case of supporting elements, recent studies have looked at the multiple methods used for hydromorphological characterization. Hydromorphology has a key role in the assessment of hydrology (i.e., the quantity and dynamics of water flow and connection to groundwater bodies), morphology (i.e., reach depth and width variation), structure and substrate of the river, structure of the riparian zone and river continuity [2].
The strengths and limitations in relation to the implementation of the WFD of a total of 139 available (European and non-European) hydromorphological characterization methods were identified under the REFORM project [3]. Overall, the majority of methods relied on the assessment of the physical habitat and did not provide sufficient consideration of physical processes or were pressure-specific. As a result, over the last few years there has been an increasing effort to develop geomorphologically based methods that facilitate the understanding of river functioning and evolution as a basis for assessing current conditions [4]. Examples of such methods are the Morphological Quality Index (MQI), the Morphological Quality Index for Monitoring (MQIm), the Geomorphic Units survey and classification System (GUS) [4,5,6,7,8,9] and the Hydromorphological Evaluation Tool (HYMET) [10].
The variables used for the characterization of hydromorphology depend upon the method being implemented. For example, the GUS system focuses on the identification of areas containing landforms created by erosion and or deposition inside or outside the river channel [6,9]. These can correspond to sedimentary units or include living or dead vegetation. Examples of geomorphic units are alternate bars, abandoned channels, cut-off channel, glides, island and riffles. The HYMET tool focuses on artificiality and sediment budget assessment [10]. The implementation of the method requires (i) the evaluation of the connectivity of the reach to the sediment production in its catchment; (ii) the analysis of the sediment transfer through the river network to the downstream reach and (iii) the quantification of the river sediment budget and artificiality. Variables such as the number of reservoirs present, the total reservoir storage volume and the quantity of produced sediment with free access to reach are estimated. More traditional methods such as the River Habitat Survey (RHS) [11] are based on descriptors of physical attributes (e.g., predominant bank material, predominant substrate within the channel, flow types and channel modification indicators).
Recent studies have looked at the harmonization of existing methods for hydromorphological characterization. Raven et al. [12] compared qualitatively three survey methods namely, the German Länder-Arbeitsgemeinschaft Wasser-Field Survey (LAWA-FS), the United Kingdom River Habitat Survey (RHS) and the French Système d’Evaluation de la Qualité du Milieu Physique (SEQ-MP). Differences in the quality assessment between the methods arose from the survey strategy and scale, data collection and analysis. Raven et al. [13] developed a methodology to test whether an UK hydromorphological assessment method (RHS) was suitable for rivers outside the UK. The aim was to assess whether the sampling scale (500 m) used was adequate and to suggest benchmarking strategies. Although some modifications were needed, the sampling scale was considered effective for the characterization of small rivers with the designation of hydro-ecological regions being identified as a key step towards a benchmarking program. Further work [14] compared the Ukrainian Field-Survey (UA-FS) and German LAWA-FS approaches to identify similarities and differences to the conformity of the outputs with the WFD requirements. Differences between methods included the assessment and interpretation of lateral erosion, sinuosity, type and depth profile, substrate diversity and special structures whereas similarities where observed for parameters describing land use, current diversity and water depth variation. Benjankar et al. [15] compared three hydromorphological characterization methods (i.e., the German LAWA, a special approach for urban rivers and a hydraulic based method), with results showing that assessment approaches developed for a certain geographic region may not be suitable for rivers in different contexts (i.e., ecoregion) and a transferable characterization approach was required. Langhans et al. [16] proposed a seven step procedure to combine assessment methods and harmonize metric outputs to a common scale from 0 to 1. The method was successfully tested for the integration of three hydromorphological methods developed in the USA (the Rapid Bioasessment Protocol-RBP), Switzerland (the Swiss Modular Concept of stream assessment-SMC) and Germany (LAWA) for four river sections. Similarly, in [2] a total of 121 hydromorphological assessment methods currently in use were compared to identify strengths, limitations and the potential to integrate different approaches. The results highlighted the need for an assessment approach based on integrated analysis, where the morphological and hydrological components drive the hydromorphological characterization. Fernandez et al. [17] reviewed more than 50 methods for river habitat characterization used worldwide. Results indicated that the key differences are the objectives for which the methods were designed, the time required for their implementation and whether the methods measure or evaluate characteristics. Overall, the parameters recorded most often in the various approaches were bank stability, channel substrate, artificial structures, riparian vegetation structure, channel dimensions, flow types, adjacent land uses and bars. The study concluded that methods that can be implemented at multiple spatial scales and provide quantitative information (i.e., indices vs. characterization protocols) were the most effective.
The WFD intercalibration exercise between member states [18] seeks to ensure comparability of ecological status boundaries and national assessment methods across Europe [19,20,21]. For this purpose, a set of six geographical intercalibration groups (GIGs) comprising waters of similar biogeophysical types (i.e., common intercalibration types) have been defined and are: Alpine, Eastern Continental, Central-Baltic, Mediterranean, Northern and Very Large Rivers (VLR) [18]. These are equivalent to the standard geographical stratification by ecoregion, where the ecoregions are geographical zones influenced by similar geophysical drivers.
The need to find pragmatic and cost-effective assessment approaches that can be adopted as common assessment methods across all countries or even as a baseline benchmark for the intercalibration exercise has been highlighted by some authors [19]. Recent studies have proposed the use of remote sensing data as a way to help achieve the harmonization of hydromorphological sampling protocols [2,22,23] at multiple spatial scales. Several tools have been developed to estimate key hydromorphological variables from remote sensing data. For example, [24] developed a method to derive valley width, active channel width and channel slope automatically from a 5 m spatial resolution DEM and 0.25–0.5 m spatial resolution orthoimage in a geographic information system (GIS) environment. Demarchi et al. [25] developed a semi-automated methodology for mapping hydromorphology indicators of rivers at regional scale using remote sensing data (i.e., high resolution near-infrared and LIDAR topography). The method enables the delineation of the natural fluvial corridor and primary riverscape units (e.g., water channel, unvegetated sediment bars, riparian densely-vegetated units) and in-stream mesohabitats. However, the resolution of commercially available remote sensing products (e.g., Lidar, aerial photography or synthetic aperture radar (SAR) data) is coarser than 10 cm and may not be sufficiently detailed for accurate hydromorphological characterization. Based on the work by [26], resolutions coarser than 2.5 cm provide biased characterizations of hydromorphological features when using artificial neural networks (ANNs) for their identification. The trade-off between high aerial imagery resolution and survey area coverage needs to be taken into account within this context. Aerial imagery at resolutions finer than 2 cm is difficult to obtain for wide-area characterization due to limitations in the current state-of-the art of the technology. Other limitations exist for commercially available remote sensing products: satellite data do not provide on-demand or temporally continuous data sets and provide oblique imagery which may result in partial spatial coverage when structures obstruct the view angle. In addition, both satellite and aircraft based remote sensing products fail to provide information under low cloud conditions.
Methodologies based on high resolution aerial imagery captured from Unmanned Aerial Vehicles (UAVs) have been proposed by several authors as potential approaches to overcome some of these limitations. For example, in [27] an UAV-based methodology for the automated recognition of hydromorphological features using ANNs was developed. Similarly, other authors [28,29] have identified photogrammetric methodologies based on high resolution UAV aerial imagery as the future tool of choice for reliable hydromorphological (i.e., physical river habitat) assessment. These techniques, albeit within their limitations [30], provide a rapid, inexpensive and increasingly accessible alternative to traditional remote sensing methods for physical habitat assessment whilst enabling quantitative microscale characterization [28]. The findings from these studies [26,27,28,30,31] raise the question as to whether frameworks based on the automated recognition of hydromorphological features from high resolution UAV aerial imagery could provide an objective and unbiased approach that can be adopted as a common assessment method across all EU member states and provide benchmark metrics for the intercalibration exercise.
In [27], we developed an UAV-based framework for the automated recognition of hydromorphological features and demonstrated an 81% average level of accuracy in feature classification for a 1.4 km reach (River Dee, Wales, UK). In [26], the level of resolution required for unbiased identification of hydromorphological features was determined for the same river reach and estimated to be 2.5 cm. Here, the aim is to assess the application of the framework to a range of different fluvial settings. This will be achieved through the following three core objectives:
  • Objective 1. To test the validity of the framework for a range of fluvial settings as identified by the WFD GIGs.
  • Objective 2. To compare the accuracy of the framework in hydromorphological feature identification within and between fluvial settings.
  • Objective 3. To interpret the outputs from (1) and (2) in line with the WFD regulatory framework.

2. Materials and Methods

2.1. Study Sites

Three sites corresponding to different GIGs (i.e., Central-Baltic, Mediterranean and VLR were selected for analysis (Figure 1 and Table 1). The sites were identified to maximize the diversity in hydromorphological characteristics along a ≈ 1 km reach and included the key hydromorphological features expected to be found in the named GIGs. From an initial set of reaches considered for analysis, only those that were representative of the spatial river variability within the waterbody where selected for further analysis. To ensure the selection was based on the current hydromorphological configuration of the waterbody, the full length of the river within each selected waterbody was assessed via multiple field-visits. The total reach length (≈1 km) was selected to exceed that used by the majority of the hydromorphological assessment methodologies used for WFD purposes, specifically those used in UK and Spain [2,3].
The Central-Baltic GIG was represented by a reach within the upper catchment of the River Dee, Wales, UK. The reach was located 30 km from the river source and 11.5 km downstream from the Bala reservoir. The river bed along the reach was characterized by gravel substrate with silt deposition in areas with low to non-perceptible flow. Bank erosion was visible along the reach with some sections presenting eroding cliffs. The banks were primarily covered by grazed grassland with few riparian trees scattered along the reach. The river length within the water body was 27.87 km [32], with the selected reach representing 5% of its totality. The overall designation for the reach was “heavily modified” [32].
The reaches representing both the Mediterranean and the VLR GIGs fell within the Jucar river basin in Spain. The Mediterranean GIG reach was located within the midland section of the river, near Motilleja (Albacete) and was characterized by a calcareous river bed with associated high entrenchment ratio. The reach contained small man-made structures that regulated the flow for irrigation purposes.
The majority of the calcareous bed channel had a layer of fine sediments due to intensive agricultural practices. The river length within the water body was 21.89 km [33,34], with the selected reach accounting for 5.5% of the total length. The overall WFD designation of the water body was “natural” [33,34].
The VLR GIG reach was located within the lower catchment of the Jucar, 500 m downstream of the Antella sluice (Antella, Valencia). The sluice diverted more than 25% of the river flow to irrigation [33,34]. Gabions were present along both banks for over 500 m at both the normal channel width and bankfull heights. The substrate within the reach was characterized by gravels with fine sediments intermittently present in some areas. The banks were dominated by invasive species (i.e., Arundo donax) with both submerged and emergent vegetation frequently present along the reach. The river length within the waterbody was 4.54 km [33,34], with the selected reach accounting for 26.4% of the total length. The WFD designation of the overall water body was “natural” [33,34].

2.2. Selection of Hydromorphological Variables

A set of variables that are major components in the WFD hydromorphological assessment in the UK and Spain were selected for automated identification (Table 2). These were specifically derived from: (i) the River Habitat Survey [11], the key method for the WFD hydromorphological assessment of rivers in UK; (ii) the UK methodology used for the designation of the WFD water body as heavily modified or artificial [36] and (iii) the equivalent Spanish methodologies [37,38,39,40,41]. The selection was based on the potential of these variables to contribute to the implementation of other WFD methods [2,3,5,7,8].
The scale at which the automated identification took place was an important consideration for the selection of hydromorphological features. This is because rivers are hierarchically organized systems; at each spatial and temporal scale there are a set of variables that are the most important in determining system behaviors and capacities [7,42]. The scale also determines the remote sensing data required for the characterization of hydromorphology [22]. Within this project, the variables selected focused on reach scale characterization which is the scale at which UAVs show optimal performance. The selected variables were divided into substrate, water, vegetation, shadows and artificial features.
Substrate features (i.e., side bars and erosion) provided information on the contemporary evidence of channel adjustment [7]. From a total of 121 hydromorphological characterization methodologies reviewed in [2], substrate related features were considered by 83 methods. Water features (i.e., riffles, glides, pools, shallow water, chute and major impacts) informed the characteristics and dimensions of the channel type and provide the physical base for habitats [7,9]. The identification of such features (i.e., flow type and change in depth [2]) was carried out within at least 34 of the hydromorphological characterization methods out of the 121 methods reviewed in [2]. Vegetation features (i.e., tree, vegetated side bar, vegetated bank, submerged free floating vegetation, emergent free floating vegetation, grass and nuisance plant species) provided evidence of the vegetation dynamics within the reach and allowed the estimation of derived parameters such as the percentage of riparian corridor under riparian vegetation or the riparian vegetation age structure [7]. They also provided information on the riparian vegetation continuity, species composition, coverage and distribution [2]. Belletti et al. [2] estimated that vegetation features (i.e., in channel vegetation and woody debris) were assessed within at least 54 hydromorphological characterization methods out of the 121 reviewed. Shadows were used to estimate the percentage of the channel with shade present whilst artificial features provided an indication of the constrains forced onto natural channel adjustments [7,9]. The presence of artificial features was recorded in over 80 methods out of the 121 reviewed in [2]. There was no clear indication of how many methods recorded shadows/shade was available.

2.3. Sampling Design and Data Collection

Aerial imagery in the visible spectrum was collected over a 19-month period (21 Apirl 2015–24 November 2016) with a range of platforms and sensors specifically selected for the area to be surveyed (Table 3 and Figure 2 and Figure 3). Ground Control Points (GCPs) were deployed along the river banks to obtain parameters for external orientation. At the Central-Baltic GIG reach, the locations of the GCP centroids were obtained from a Leica GS14 Base and Rover Real Time Kinematic (RTK) GPS (Leica Geosystems AG, Heerbrugg, Switzerland) with a position accuracy of more than 0.02 m in the X, Y and Z dimensions. For all the other reaches, a Leica 1200 GPS linked to a GNSS permanent reference station with a position accuracy of 0.02 m in planimetry and 0.03 m in altimetry was used. Two different GCP designs were used in the study: 1 m × 1 m black squared GCPs with white opposite facing triangles where used for the Central-Baltic reach whereas white circular (0.30 m diameter) GCPs with concentric black rings were used for the Mediterranean and VLR reaches (Figure 2 and Figure 3).
Various vertical-take-off-and-landing UAV platforms (Figure 3) were used for image collection and included a Falcon 8 Trinity (ASCTEC, Krailling, Germany), an IRIS9+ (3DR, Berkeley, CA, USA) and a md4-1000 (Microdrones, Inc., Kreuztal, Germany). The Falcon 8 Trinity platform is a 0.77 m × 0.82 m × 0.12 m octocopter equipped with an uBlox LEA 6S GPS (uBlox, Thalwil, Switzerland). The platform has a vertical-take-off weight of 1.8 kg with the sensor payload (Sony Alpha) only requiring 0.34 kg. Flight endurance under optimal conditions and with fully charged lithium polymer (LiPo) batteries (6250 mAh) is 22 min. The IRIS 9+ is a 0.50 m diameter (from rotor hub to rotor hub) quadcopter equipped with a 3DR u-Blox GPS (uBlox, Thalwil, Switzerland) and a 0.4 kg payload capacity, with the camera requiring only 0.14 kg. Flight endurance under optimal conditions with the payload incorporated and with fully charged LiPo batteries (5100 mAh) is between 16–22 min. The md4-1000 is a 1.00 m diameter (from rotor hub to rotor hub) quadcopter with a 0.224 kg embedded camera, a total payload capacity of 1.2 kg and a maximum take of weight of 6 kg. The md4-1000 is equipped with an mdIMU and GNSS u-Blox6 positioning system. The flight endurance with fully charged 13,000 mAh LiPo batteries is 45 min. The ground sampling distance (GSD) for each reach was calculated to achieve a desired target imagery resolution of 2.5 cm based on the camera focal length and the regulatory airspace constraints within each country (Table 3). All flights were carried out by qualified UAV pilots following the country specific aviation regulation.
The methodology was first used and tested in the Central-Baltic GIG reach [26,27]. The lessons learnt from this deployment were transferred to the other reaches. Data were collected under low flow conditions and with a constant volumetric flow rate (Table 4). The number of flights required was determined by the platform, sensor and reach configuration. Each platform was configured to follow a flight plan that captured images at waypoints that achieved 60% along track and 40% across track image overlap. Longitudinal and cross-sectional multi-passes ensured full spatial coverage.
An ancillary data set was collected to aid visual identification of hydromorphological features. This consisted of a detailed rapid habitat map of the reach, associated photographic evidence of the records, accurate RTK GPS location of key features and measurements of water depth and velocity. The habitat map was obtained by walking along the reach and drawing polygons indicating the extent of the observed features described in Table 2. The RTK measurements and photographs were collected wherever key features were present to indicate their position and extent (subject to access). Both water depth and velocity readings were obtained by use of a radio controlled boat with integrated ADCP or a hand held ADCP and a standard surveying metric rod. All the information for a single reach was collected within a maximum of a five day period to ensure spatio-temporal collocation of multiple-sensor measurements.

2.4. Photogrammetry

The imagery collected at each reach was assessed based on quality and spatial coverage for inclusion in the photogrammetric process (Figure 4, step 3). Blurred and distorted frames were excluded [46] for the generation of the standard geomatic products (i.e., orthoimage, digital terrain model and point cloud) via Photoscan Pro version 1.1.6 (Agisoft LLC, St. Petersburg, Russia). The photogrammetric process requires all the frames to be georeferenced (i.e., scale, translate and rotate) into a target Geodetic System (i.e., the World Geodetic System, WGS84) using the GCPs coordinates to minimize geometric distortion. For this purpose, the centroid of each GCP was manually identified in all frames and assigned the corresponding field RTK GPS coordinates. The processing time required for the photogrammetric process for each reach based on the performance of a computer with an Intel Core i7-5820k 3.30 GHz processor, 32 Gb RAM and two graphic cards (Geoforce GTX 980 and Qadro K2200, NVIDIA, Santaclara, CA, USA) is summarized in Table 3. The coregistration errors in X, Y and Z were automatically derived by Photoscan Agisoft the difference between the positions of the GCP centroids measured through RTK GPS and the coordinates derived from the imagery.

2.5. Automatic Classification

The automated identification of hydromorphological features was carried out via the Leaf Area Index Calculation (LAIC) software [47] (Figure 4, step 4), originally developed to discriminate green canopy cover from ground, stones and shadow background from high resolution aerial imagery [47,48]. LAIC bases the classification on an ANN approach (i.e., supervised classification technique) that segments the spectral domain of the RGB imagery into areas that directly relate to the features of interest. Key to the overall approach is the training process by which representative samples of target hydromorphological features are selected using a k-means clustering algorithm and used as the reference signature for the classification of the entire orthoimage.
The classification process was described in more detail in [27] but in brief, prior to the ANN training process, the RGB image was transformed into CIE L × a × b (CIELAB) space to perform a cluster analysis to identify the different river features. From the three parameters describing the CIELAB space (i.e., lightness (L), green to red scale (a) and blue to yellow scale (b)), only a and b were taken into account by the clustering algorithm. The number of clusters (k) depended on the feature being identified and was determined following an iterative process that increased k by one up to a maximum of ten clusters until visually satisfactory results were obtained. Within this context, visually satisfactory results required the image outputs to show that the feature of interest had been adequately identified. This supervised method was used as a basis to calibrate a Multilayer Perceptron (MP) ANN, which was applied to the remaining image. The ANN was based on three consecutive layers composed of inter-connected nodes. Within the context of this study, the outputs from the software were a classified map of the reach and estimates of the total area allocated to each of the hydromorphological features identified.

2.6. Statistical Analysis

The performance of the automated classification technique was assessed independently for each reach. The hydromorphological features were identified at each point of a 2 m × 2 m regular grid super imposed on each of the sampled reaches. The identification was based on the high resolution aerial imagery, the ancillary data sets encompassing the in-situ maps, the RTK measurements of key hydromorphological features and the photographic evidence of their distribution. The resulting output was considered to be the ground truth data set.
The performance of the ANN was assessed following [26,27] through confusion matrices and the estimation of derived standard metrics (Equations (1)–(5)) known as: overall accuracy (AC), true positive ratio (TPR), true negative ratio (TNR), false positive ratio (FPR) and false negative ratio (FNR). These were calculated as follows:
A C = i = 1 I ( T N i + T P i T N i + T P i + F N i + F P i )
T P R i = T P i F N i + T P i
T N R i = T N i T N i + F P i
F N R i = F N i F N i + T P i
F P R i = F P i T N i + F P i
where TPi (True Positives) is the number of points correctly identified as class i, FNi (False Negatives) is the number of points incorrectly rejected as class i, TNi (True Negatives) is the number of points correctly rejected as class i, FPi (False Positives) is the number of points incorrectly identified as class i and I is the total number of classes identified.
TPRi, TNRi FNRi and FPRi were estimated for each of the features of interest whereas AC was a single value of overall classification performance. AC, as well as all the other ratios, ranged from 0 to 1. TPRi and TNRi quantifed the power of LAIC at classifying features correctly when compared to the ground truth whereas FNRi and FPRi showed misclassification rates.

3. Results

The image coregistration errors in X, Y and Z were below 2.7 cm for all reaches, with the VLR reach presenting the largest errors in X and Y (≈2.7 cm) and the smallest errors in Z (≈1 cm) (Table 5). The smallest coregistration errors in X and Y were found for the Mediterranean (≈1.06 cm) and the Central-Baltic (≈1 cm) reaches, respectively. The dimensions of the sampled reach (i.e., width and length) (Table 1) determined the number of points defined by the 2 m × 2 m regular grid and considered for analysis (Table 5). These varied between 4915 points (VLR) and 13,085 points (Central-Baltic).
Based on the visual classification of the 2 m × 2 m grid, the type and extent of hydromorphological features varied considerably between reaches (Figure 5). The configuration within the Central-Baltic was predominantly dominated by riffles (26%), vegetation (25%), shallow waters (20%), deep water (16%) and bars (11%) with each of the other classes identified accounting for less than 3% of the total number of points. The dominant features within the Mediterranean reach include vegetation (65%), shallow water (19%) and deep water (10%) with all the other classes not accounting for more than 5% of the points each. These results were consistent with the pattern observed for the VLR reach, where vegetation, shallow water and deep water accounted for 39%, 35% and 8% of the points, respectively. Each of the remaining classes identified accounted for less than 8% of the classified points. The presence of shadows was significant along the Central-Baltic (≈2%) and VLR (8%) reaches with no shadows being identified along the Mediterranean reach.
The overall accuracy of the automated classification was above 50% for all the reaches, with the Central-Baltic reach presenting the best performance (81%), followed by the Mediterranean (71%) and the VLR (50%) reaches. The patterns of performance, as shown by the TPR, TNR, FNR and FPR values (Table 6), differed between reaches. Overall, the ANN was able to identify vegetation successfully with TPR values above 74% for all sites. Misclassification of vegetation features primarily occurred with riffles (Central-Baltic) and shallow waters (Central-Baltic, Mediterranean and VLR) (Table 6 and Table 7). Vegetation features were also classified as shadow (Central-Baltic and VLR) and deep waters (Mediterranean). These results were consistent with the performance of the ANN for the detection of riffles and side bars (TPR > 73% and FNR < 27% and FPR ≤ 6%) for the Central-Baltic and Mediterranean reaches. However, the performance in riffle and side bar identification decreased down to 0 for the VLR reach (FNR = 1). Overall, riffles were misclassified as shallow water (Central-Baltic and VLR), shadow (Central-Baltic), vegetation (Central-Baltic, Mediterranean and VLR) and deep water (VLR), whereas side bars were mainly misclassified as vegetation (Central-Baltic and VLR), riffles (Mediterranean) and shallow waters (Mediterranean and VLR) (Table 6 and Table 7).
For the specific case of deep waters, the ANN performance (TPR) was above 55% in all instances, with the best performances being observed in the Central-Baltic (TPR > 92%, FNR < 8%) and the VLR reaches (TPR > 66%, FNR < 34%). Misclassification occurred with riffles (Central-Baltic), shallow water (Central-Baltic and Mediterranean) and vegetation (Central-Baltic, Mediterranean and VLR) (Table 6 and Table 7). The performance in shallow water identification was consistent across the three reaches (TPR > 51%, FNR < 49% and FPR < 14%). Here, misclassification was present in shadows (Central-Baltic and VLR), vegetation (Central-Baltic, Mediterranean and VLR), riffles (Central-Baltic) and deep-water (Mediterranean) (Table 6 and Table 7). Shadows, when present, were successfully identified with TPR, FNR and FPR values >53%, <47% and <8%, respectively with misclassification mainly occurring with vegetation (VLR). Poor results were obtained across the three reaches when the ANN focused on the identification of chutes (TPR < 29%), erosion (TPR < 8%), structures (TPR = 0%) and major impacts (Table 2) such as pollution (TPR = 50%). The number of validation points available for these features was below 20 for the Mediterranean GIG site and 55 for the VLR GIG site (Table 7).

4. Discussion

The work reported in this paper aimed at assessing the transferability of an already tested framework for hydromorphological feature characterization [26,27] to a range of different fluvial settings. The three objectives addressed within this context were: (i) to test the validity of the framework to a range of fluvial settings as identified by the GIGs; (ii) to compare the accuracy of the framework in hydromorphological feature identification within and between fluvial settings and (iii) to interpret the outputs from (i) and (ii) in line with the WFD regulatory framework. The following sections discuss the outcomes within the context of these objectives.

4.1. UAV Framework Performance at the GIGs Sites

The proposed framework enables the identification of hydromorphological features with accuracies above 50% for all the GIGs considered. Overall, the ANN works robustly for the identification of vegetation, deep water, shallow water, side bars, riffles and shadows. For the particular case of the VLR GIGs, the ANN appears to underperform in the detection of both riffles and side bars. Riffles were systematically classified as deep water, shallow water or vegetation (submerged) whereas side bars were classified as vegetation or shallow waters. The underperformance detected for the VLR site may be the result of the configuration of the reach—i.e., transitional (from deep to shallow and vice versa) flowing water along a wide gravel bed channel that generated smooth rippled surfaces undetectable by the ANN. Within this setting, the ANN seem to be unable to detect the difference between riffles and shallow/deep waters. This is because all the riffles were present within the extent of the transitional water and could be classified as both feature categories simultaneously (e.g., riffles and shallow waters); any validation point falling under a riffle and being classified as shallow water (or vice versa) cannot be automatically considered a misclassification error.
The calcareous nature of the VLR reach may also affect the characteristics of the RGB aerial imagery captured; the water presented a low turbidity and a turquoise-like color that made feature identification (i.e., identification of differences in depth, velocity and substrate, amongst other parameters) challenging. This could explain why side bars (and riffles) were confused with vegetation as the color of both was similar under the survey conditions and why the distinction between riffles, deep and shallow waters cannot be accurately identified by the ANN. In [27], we explored in detail where the confusion between shallow and deep waters occurred along the Central-Baltic GIG site; the transition zone from deep to shallow areas were the key sources of feature misclassification. This lack of ability to identify the transition zones applies to all the reaches analyzed in the present study. The transition zone occupied larger proportions of the surveyed area in the Mediterranean and VLR reaches.
The ANN performs poorly on the classification of chutes, artificial structures, erosion and pollution. These features appeared occasionally along some of the reaches and were not representative enough for the ANN to register them as independent clusters (i.e., they accounted for less than 8% of the classified points). In the particular case of pollution, the ANN provides encouraging results (TPR > 53%). However, the limited number of validation points (i.e., 16 points along the Mediterranean reach) available for that particular feature class does not allow for a conclusive statement on the performance of the ANN for pollution identification. Note that pollution appeared in the Mediterranean GIG reach only and along a backwater that collected white foam and plastics derived from upstream activity.

4.2. Potential Technical Improvements

For the proposed framework to be used for large scale (>1 km) hydromorphological characterization, it is necessary to increase its cost-effectiveness via the reduction of the number of GCP to be deployed. Recent work by several authors (e.g., [48,49]) has already contributed to address this gap in knowledge. This, coupled with increased UAV battery performance, data capture and less CPU demanding software makes UAV based frameworks a plausible option for wide-area (>1 km reach) robust and accurate hydromorphological assessment. Further work should compare the performance of the proposed framework with existing classification techniques for geomorphological environments (e.g., [50,51]).
The work presented here focuses on the characterization of hydromorphological variables at reach level. It is not yet certain whether the framework can be applied to larger spatial scales (e.g., catchment [7]) and achieve the same efficiency as current remote sensing methodologies [22]. Woodget et al. [30] identify the spatial coverage as a key challenge for the routine operational use of UAVs and digital photogrammetry for river habitat mapping. The strong trade-off between resolution and spatial coverage [30] suggests that some aspects of the hydromorphological characterization for WFD purposes may be achieved from UAV imagery but that complementary remote sensing methods may be required to address the remaining aspects [22]. For example, the strong reliance on image texture required for the implementation of Structure-from-Motion (SfM) may also compromise the wide-area implementation of the framework. If the texture of the imagery is compromised, SfM will not detect matching features between overlapping images and fail to produce an orthoimage of the surveyed area [30]. Further research is therefore required to address these points before the framework can be adopted for wide-area monitoring. In addition, further consideration needs to be given to the type of sensors to use on the UAV platform. The RGB imagery collected for this study enables the qualitative assessment of hydromorphological features (e.g., deep/shallow waters) but does not facilitate the quantitative estimation of depth and velocity within the channel reach [30]. It has been highlighted by some authors that the RGB imagery presents limited radiometric resolution which obstructs the restitution of topography in darker parts (e.g., shadow and deep water) [30]. Sensors able to provide multispectral and RAW format imagery may be beneficial and should be considered in further developments of the framework. The variability in the quality of the geomatic products generated from high resolution UAV aerial imagery has also been raised as a general concern by some authors [52]. For example, the optimal selection of GCP can improve the variability of the Digital Elevation Model from 37 mm to 16 mm. The magnitude of these changes in the DEM may be relevant when changes in hydromorphological characteristics over time are required [52] but not relevant where regulatory assessments are the primary consideration.

4.3. Framework Output Interpretation in Line with the WFD

The results obtained indicate that the methodology can be successfully transferred within and between river GIG types and countries but that water turbidity, color and underlying river bed substrate may affect classification performance. The hydromorphological features successfully identified (i.e., vegetation, deep water, shallow water, side bars, riffles and shadows) are common to the majority of hydromorphological characterization methods reviewed in [2] and are therefore of great relevance for the transferability of this framework. They will be able to contribute to the implementation of the WFD in UK as they align with the descriptors included in the RHS methodology [11] and the designation of WFD artificial and heavily modified water bodies [36]. Similarly, they will aid the implementation of the Spanish hydromorphological assessments [37,38,39,40,41] and contribute to the implementation of WFD methodologies in other countries [2,3,5,7,8].
The framework outlined in this paper is a step forward towards the accurate identification of in-channel and bank features and could be integrated with complementary frameworks. Recent work by other authors have shown that other riverine features can be identified with similar UAV based methodologies. In [53], supervised machine learning approaches were used to identify different types of macrophytes. In [54], object oriented classification was used to map dead wood presence, whereas in [29], geomatic products derived from UAV aerial imagery were used to characterize channel substrate. Further research should focus on integrating these scientific advances into a single methodological framework that enables the objective and automated characterization of key WFD water elements. This will address the need identified in previous works [2] for the development of an integrated hydromorphological analysis that includes both morphological and hydrological components to evaluate and classify hydromorphological state and quality. In turn, this will contribute to addressing some of the current limitations in the WFD intercalibration exercise.
At present it is not possible to identify fully the consequences of the EU member states using different methodological approaches to assess the ecological status of water bodies. The intercalibration exercises have not provided a clear answer yet [14,21] and as a result, there is still insufficient information about how the different methods should be implemented and integrated to achieve readily comparable assessments [20]. It is possible that the consequences will be both financial and legal in terms of the measures and efforts undertaken to achieve good ecological status or potential under the WFD [12,14]. An objective and fit for purpose hydromorphological characterization framework that is widely adopted within all member states will facilitate comparison of the results obtained across Europe and will provide a level playing field. The UAV based framework presented here would contribute for example, to a more consistent designation and overall hydromorphological assessment by addressing some of the current limitations, namely: inadequate characterization of pressure gradients within the GIGs, difficulty comparing results obtained from methods with different metrics, lack of comparable data sets and implementation of subjective methods [14] through the provision of a standardized approach to understanding the hydromorphological quality in rivers.
An interesting application of the framework described in this paper would be its use for the designation of the waterbody as “natural”, “heavily modified” or “artificial”. Currently such designations can be influenced to a significant extent by management and external influence. For the case study areas, the rivers studied representing the Mediterranean and VLR GIGs had been designated by the Spanish authorities as having “natural” hydromorphological conditions [33]. In the UK, the reach representing the Central-Baltic GIGs was less modified than the Spanish reaches but had been designated by the authorities as “highly modified” [35]. The differences in the approaches taken for the designation of the waterbodies translates into different requirements in terms of the measures to be implemented and the associated level of investment required to improve the ecological status of the water bodies. This mismatch in the implementation of the WFD could be addressed if a common framework for hydromorphological characterization were to be adopted. The techniques developed could also be useful in a number of situations where hydromorphological characterization is required aside from the WFD implementations such as river restoration appraisal [55].

5. Conclusions

This paper reports the transferability of an already validated ANN based framework for hydromorphological characterization of three different fluvial settings as defined by the WFD GIGs. The framework provides overall accuracies greater than 50% for all the reaches but shows some limitations when applied to calcareous rivers falling under the VLR GIG. Deep waters, shallow waters, vegetation, side bars, riffles and shadows are successfully identified within all reaches, with the ANN underperforming on the identification of side bars and riffles for the VLR GIG. Further work is required to develop more reliable algorithms that can be incorporated into the existing ANN for the accurate detection of a wider range of features, including chutes, erosion, structures and major impacts (i.e., pollution). This work also highlights the need to develop objective and reliable hydromorphological assessment frameworks that are widely adopted by all EU member states. The framework presented here is a step forward towards that goal.

Acknowledgments

We would like to thank the Environment Agency and EPSRC for funding part of this project under the EPSRC Industrial Case Studentship voucher number 08002930. The underlying data are confidential and cannot be shared. Special thanks go to Natural Resources Wales for their help and support with data collection. The authors acknowledge financial support from the Castilla-La Mancha Regional Government (Spain) under the Postdoctoral Scholarship Programa Operative 2007–2013 de Castilla-La Mancha. Special thanks go to the journal reviewers and editor; they made this manuscript stronger with their positive criticism and valuable comments.

Author Contributions

M.R.C. is the principal investigator and corresponding author. She led and supervised the overall research and field data collection. M.R.C. structured and wrote the paper in collaboration with R.B.G. R.B.G. helped with the flight planning and overall data collection in UK and Spain. R.B.G. processed the UAV imagery and helped with the interpretation of the results. J.F.O. enabled equipment provision and data collection in Spain. J.F.O. coordinated the financial aspects of data collection in the Spanish reaches. He also provided advice and guidance on the regulatory implication of the WFD in Spain. P.L. reviewed the robustness of the approach. P.L. provided input on the policy and regulatory implications of the research outcomes and coordinated communication with key project partners. R.W. coordinated the financial aspects of the research, enabled equipment provision and facilitated fieldwork arrangements in UK. She highlighted research priorities and provided operational advice as and when required. All authors were key for the interpretation of the results.

Conflicts of Interest

The authors declare no conflict of interest.

References and Notes

  1. European Commission. Directive 2000/60/EC of the European Parliament and of the Council of 23 October 2000 establishing a framework for Community action in the field of water policy. Off. J. Eur. Union 2000, L327, 1–72. [Google Scholar]
  2. Belletti, B.; Rinaldi, M.; Buijse, A.D.; Gurnell, A.M.; Mosselman, E. A review of assessment methods for river hydromorphology. Environ. Earth Sci. 2014, 73, 2079–2100. [Google Scholar] [CrossRef]
  3. Rinaldi, M.; Belleti, B.; Van de Bund, W.; Bertoldi, W.; Gurnell, A.; Buijse, T.M.E.; Mosselman, E. Review on Eco-Hydromorphological Methods. Available online: http://www.reformrivers.eu/d-11-review-eco-hydromorphological-methods (accessed on 26 September 2017).
  4. Rinaldi, M.; Belletti, B.; Bussettini, M.; Comiti, F.; Golfieri, B.; Lastoria, B.; Marchese, E.; Nardi, L.; Surian, N. New tools for the hydromorphological assessment and monitoring of European streams. J. Environ. Manag. 2016, 202, 363–378. [Google Scholar] [CrossRef] [PubMed]
  5. Gurnell, M.; Bussettini, M.; Camenen, B.; González Del Tánago, M.; Grabowski, R.; Hendriks, D.; Henshaw, A.; Latapie, A.; Rinaldi, M. Multi-Scale Framework and Indicators of Hydromorphological Processes and Forms I. Main Report. Available online: http://www.reformrivers.eu/multi-scale-framework-and-indicators-hydromorphological-processes-and-forms-i-main-report (accessed on 26 September 2017).
  6. Belletti, B.; Rinaldi, M.; Bussettini, M.; Comiti, F.; Gurnell, A.M.; Mao, L.; Nardi, L.; Vezza, P. Characterising physical habitats and fluvial hydromorphology: A new system for the survey and classification of river geomorphic units. Geomorphology 2017, 283, 143–147. [Google Scholar] [CrossRef]
  7. Gurnell, A.M.; Rinaldi, M.; Belletti, B.; Bizzi, S.; Blamauer, B.; Braca, G.; Buijse, A.D.; Bussettini, M.; Camenen, B.; Comiti, F.; et al. A multi-scale hierarchical framework for developing understanding of river behaviour to support river management. Aquat. Sci. 2016, 78, 1–16. [Google Scholar] [CrossRef] [Green Version]
  8. Rinaldi, M.; Gurnell, A.M.; Belletti, B.; Berga Cano, M.I.; Bizzi, S.; Bussettini, M.; González del Tánago, M.; Grabowski, R.; Habersak, H.; Klösch, M.; et al. Final Report on Methods, Models, Tools to Asses the Hydromorphology of Rivers. Available online: http://www.reformrivers.eu/deliverables/d62-final-report-methods-models-tools-assess-hydromorphology-rivers (accessed on 26 September 2017).
  9. Rinaldi, M.; Belletti, B.; Comiti, F.; Nardi, L.; Bussettini, M.; Mao, L.; Gurnell, A.M. The Geomorphic Unit Survey and Classification System (GUS). Available online: http://www.reformrivers.eu/geomorphic-units-survey-and-classification-system-gus (accessed on 26 September 2017).
  10. Klösch, M.; Habersack, H. The Hydromorphological Evaluation Tool (HYMET). Geomorphology 2017, 291, 143–158. [Google Scholar] [CrossRef]
  11. Environment Agency. River Habitat Survey in Britain and Ireland; Environment Agency: Bristol, UK, 2003.
  12. Raven, P.J.; Holmes, N.T.H.; Charrier, P.; Dawson, F.H.; Naura, M.; Boon, P.J. Towards a harmonized approach for hydromorphological assessment of rivers in Europe: A qualitative comparison of three survey methods. Aquat. Conserv. Mar. Freshw. Ecosyst. 2002, 12, 405–424. [Google Scholar] [CrossRef]
  13. Raven, P.J.; Holmes, N.T.H.; Vaughan, I.P.; Dawson, F.H.; Scarlett, P. Benchmarking habitat quality: Observations using River Habitat Survey on near-natural streams and rivers in northern and western Europe. Aquat. Conserv. Mar. Freshw. Ecosyst. 2010, 20, S13–S30. [Google Scholar] [CrossRef]
  14. Scheifhacken, N.; Haase, U.; Gram-Radu, L.; Kozovyi, R.; Berendonk, T.U. How to assess hydromorphology? A comparison of Ukrainian and German approaches. Environ. Earth Sci. 2012, 65, 1483–1499. [Google Scholar] [CrossRef]
  15. Benjankar, R.; Koenig, F.; Tonina, D. Comparison of hydromorphological assessment methods: Application to the Boise River, USA. J. Hydrol. 2013, 492, 128–138. [Google Scholar] [CrossRef]
  16. Langhans, S.D.; Lienert, J.; Schuwirth, N.; Reichert, P. How to make river assessments comparable: A demonstration for hydromorphology. Ecol. Indic. 2013, 32, 264–275. [Google Scholar] [CrossRef]
  17. Fernández, D.; Barquín, J.; Raven, P.J. A review of river habitat characterisation methods: Indices vs. characterisation protocols. Limnetica 2011, 30, 217–234. [Google Scholar]
  18. European Commission. Commission Decision of 30 October 2008 establishing, pursuant to Directive 2000/60/EC of the European Parliament and the Council, the values of the Member State monitoring system classifications as a result of the intercalibration exercise 2008/915/EC. Off. J. Eur. Union 2008, L332, 20–44. [Google Scholar]
  19. Poikane, S.; Zampoukas, N.; Borja, A.; Davies, S.P.; van de Bund, W.; Birk, S. Intercalibration of aquatic ecological assessment methods in the European Union: Lessons learned and way forward. Environ. Sci. Policy 2014, 44, 237–246. [Google Scholar] [CrossRef]
  20. Poikane, S.; Birk, S.; Böhmer, J.; Carvalho, L.; de Hoyos, C.; Gassner, H.; Hellsten, S.; Kelly, M.; Lyche Solheim, A.; Olin, M.; et al. A hitchhiker’s guide to European lake ecological assessment and intercalibration. Ecol. Indic. 2015, 52, 533–544. [Google Scholar] [CrossRef]
  21. Poikane, S.; Johnson, R.K.; Sandin, L.; Schartau, A.K.; Solimini, A.G.; Urbanič, G.; Arbačiauskas, K.; Aroviita, J.; Gabriels, W.; Miler, O.; et al. Benthic macroinvertebrates in lake ecological assessment: A review of methods, intercalibration and practical recommendations. Sci. Total Environ. 2016, 543, 123–134. [Google Scholar] [CrossRef] [PubMed]
  22. Bizzi, S.; Demarchi, L.; Grabowski, R.C.; Weissteiner, C.J.; Van de Bund, W. The use of remote sensing to characterise hydromorphological properties of European rivers. Aquat. Sci. 2016, 78, 57–70. [Google Scholar] [CrossRef] [Green Version]
  23. Gurnell, A.M.; Rinaldi, M.; Buijse, A.D.; Brierley, G.; Piégay, H. Hydromorphological frameworks: Emerging trajectories. Aquat. Sci. 2016, 78, 135–138. [Google Scholar] [CrossRef]
  24. Martínez-Fernández, V.; Solana-Gutiérrez, J.; González del Tánago, M.; García de Jalón, D. Automatic procedures for river reach delineation: Univariate and multivariate approaches in a fluvial context. Geomorphology 2016, 253, 38–47. [Google Scholar] [CrossRef]
  25. Demarchi, L.; Bizzi, S.; Piégay, H. Hierarchical Object-Based Mapping of Riverscape Units and in-Stream Mesohabitats Using LiDAR and VHR Imagery. Remote Sens. 2016, 8, 97. [Google Scholar] [CrossRef] [Green Version]
  26. Rivas Casado, M.; Ballesteros González, R.; Wright, R.; Bellamy, P. Quantifying the effect of aerial imagery resolution in automated hydromorphological river characterisation. Remote Sens. 2016, 8, 650. [Google Scholar] [CrossRef]
  27. Rivas Casado, M.; Ballesteros González, R.; Kriechbaumer, T.; Veal, A. Automated identification of river hydromorphological features using UAV high resolution aerial imagery. Sensors 2015, 15, 27969–27989. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Woodget, A.S.; Visser, F.; Maddock, I.P.; Carbonneau, P.E. The accuracy and reliability of traditional surface flow type mapping: Is it time for a new method of characterizing physical river habitat? River Res. Appl. 2016, 39, 1902–1914. [Google Scholar] [CrossRef]
  29. Woodget, A.S.; Carbonneau, P.E.; Visser, F.; Maddock, I.P. Quantifying submerged fluvial topography using hyperspatial resolution UAS imagery and structure from motion photogrammetry. Earth Surf. Process. Landf. 2015, 40, 47–64. [Google Scholar] [CrossRef]
  30. Woodget, A.; Austrums, R.; Maddock, I.P.; Habit, E. Drones and digital photogrammetry: From classifications to continuums for monitoring river habitat and hydromorphology. Wiley Interdiscip. Rev. Water 2017, 4, e1222. [Google Scholar] [CrossRef]
  31. DeBell, L.; Anderson, K.; Brazier, R.E.; King, N.; Jones, L. Water resource management at catchment scales using lightweight UAVs: Current capabilities and future perspectives. J. Unmanned Veh. Syst. 2016, 4, 7–30. [Google Scholar] [CrossRef]
  32. Environment Agency. Catchment Data Explorer. Available online: http://environment.data.gov.uk/catchment-planning/ManagementCatchment/3024 (accessed on 30 December 2016).
  33. Ministerio de Agricultura Alimentación y Medio Ambiente. Confederación Hidrográfica del Júcar. Plan Hidrológico de la Demarcación Hidrográfica del Júcar. Available online: http://www.chj.es/es-es/medioambiente/planificacionhidrologica/Paginas/PHC-2015-2021-Plan-Hidrologico-cuenca.aspx (accessed on 5 March 2017).
  34. Ministerio de Agricultura Alimentación y Medio Ambiente. Real Decreto 1/2016, de 8 de enero, por el que se aprueba la revisión de los Planes Hidrológicos de las demarcaciones hidrográficas del Cantábrico Occidental, Guadalquivir, Ceuta, Melilla, Segura y Júcar, y de la parte española de las demarcaciones hidrográficas del Cantábrico Oriental, Miño-Sil, Duero, Tajo, Guadiana y Ebro. Boletín Of. Estado 2016, 16, 2972–4301. [Google Scholar]
  35. Water Watch Wales Map Gallery. WFD Cycle 2 Rivers and Water-Bodies in Wales. Available online: https://nrw.maps.arcgis.com/apps/webappviewer/index.html?id=2176397a06d64731af8b21fd69a143f6 (accessed on 8 March 2017).
  36. Environment Agency. Technical Assessment Method (Rivers): Morphological Alteration/Identification of pHMWBs and pAWBs (Drainage Channels).
  37. Ministerio de Agricultura Alimentación y Medio Rural y Marino. Orden ARM/2656/2008, de 10 de septiembre, por la que se aprueba la instruccion y planificacion hidrologica. Boletín Of. Estado 2008, 229, 38472–38582. [Google Scholar]
  38. Ministerio de Agricultura Alimentación y Medio Ambiente Real Decreto 817/2015, de 11 de septiembre, por el que se establecen los criterios de seguimiento y evaluacion del estado de las aguas superficiales y las normas de calidad ambiental. Boletín Of. Estado 2015, 219, 80582–80677.
  39. Ministerio de Alimentación Agricultura y Medio Ambiente. Confederación Hidrográfica del Júcar. Plan hidrológico de la demarcación hidrográfica del Júcar. Memoria-Anejo 1. Designación de Masas de Agua Artificiales y muy Modificadas. Available online: http://www.chj.es/descargas/ProyectosOPH/Consulta%20publica/Anejos/PHJ_Anejo01_MAMM.pdf (accessed on 5 March 2017).
  40. Ministerio de Alimentación Agricultura y Medio Ambiente. Confederación Hidrográfica del Júcar. Plan hidrológico de la demarcación hidrográfica del Júcar. Memoria-Anejo 12. Evaluación del estado de las masas de agua superficiales y subterránea. Available online: http://www.chj.es/es-es/medioambiente/planificacionhidrologica/Paginas/PHC-2015-2021-Plan-Hidrologico-cuenca.aspx (accessed on 5 March 2017).
  41. Ministerio de Alimentación Agricultura y Medio Ambiente; Confederación Hidrográfica del Ebro. Protocolos de muestreo y análisis para indicadores hidromorfológicos. In Metodología para establecimiento del estado ecológico según la Directiva Marco del Agua; Ministerio de Alimentación Agricultura y Medio Ambiente: Madrid, Spain, 2013. [Google Scholar]
  42. Frissell, C.A.; Liss, W.J.; Warren, C.E.; Hurley, M.D. A Hierarchical Framework for Stream Habitat Classification: Viewing Streams in a Watershed Context. Environ. Manag. 1986, 10, 199–214. [Google Scholar] [CrossRef]
  43. The National Registry for Adoption (NRFA). Station Mean Flow Data. Available online: http://nrfa.ceh.ac.uk/data/station/meanflow/67001 (accessed on 31 December 2016).
  44. Ministerio de Agricultural y Pesca Alimentación y Medio Ambiente. Sistema de información agroclimática para el regadío. Available online: http://eportal.mapama.gob.es/websiar/SeleccionParametrosMap.aspx?dst=1 (accessed on 1 May 2017).
  45. Ministerio de Alimentación Agricultura y Medio Ambiente. Confederación Hidrográfica del Júcar. Plan hidrológico de la demarcación hidrográfica del Júcar. Memoria-Anejo 5. Régimen de caudales ecológicos. Available online: http://www.chj.es/descargas/ProyectosOPH/Consulta%20publica/Anejos/PHJ_Anejo05_QEco.pdf (accessed on 5 March 2017).
  46. Ribeiro-Gomes, K.; Hernández-Lopez, D.; Ballesteros González, R.; Moreno, M.A. Approximate georeferencing and automatic blurred image detection to reduce the costs of UAV use in environmental and agricultural applications. Biosyst. Eng. 2016, 151, 308–327. [Google Scholar] [CrossRef]
  47. Córcoles, J.I.; Ortega, J.F.; Hernández, D.; Moreno, M.A. Estimation of leaf area index in onion (Allium cepa L.) using an unmanned aerial vehicle. Biosyst. Eng. 2013, 115, 31–42. [Google Scholar] [CrossRef]
  48. Ballesteros, R.; Ortega, J.F.; Hernández, D.; Moreno, M.A. Applications of georeferenced high-resolution images obtained with unmanned aerial vehicles. Part I: Description of image acquisition and processing. Precis. Agric. 2014, 15, 579–592. [Google Scholar] [CrossRef]
  49. Carbonneau, P.E.; Dietrich, J.T. Cost-Effective non-metric photogrammetry from consumer-grade sUAS: Implications for direct georeferencing of structure from motion photogrammetry. Earth Surf. Process. Landf. 2017, 42, 473–486. [Google Scholar] [CrossRef]
  50. Brodu, N.; Lague, D. 3D terrestrial lidar data classification of complex natural scenes using a multi-scale dimensionality criterion: Applications in geomorphology. ISPRS J. Photogramm. Remote Sens. 2012, 68, 121–134. [Google Scholar] [CrossRef] [Green Version]
  51. Lague, D.; Brodu, N.; Leroux, J. Accurate 3D comparison of complex topography with terrestrial laser scanner: Application to the Rangitikei canyon (N-Z). ISPRS J. Photogramm. Remote Sens. 2013, 82, 10–26. [Google Scholar] [CrossRef]
  52. James, M.R.; Robson, S.; d’Oleire-Oltmanns, S.; Niethammer, U. Optimising UAV topographic surveys processed with structure-from-motion: Ground control quality, quantity and bundle adjustment. Geomorphology 2017, 280, 51–66. [Google Scholar] [CrossRef]
  53. Göktoǧan, A.H.; Sukkarieh, S.; Bryson, M.; Randle, J.; Lupton, T.; Hung, C. A rotary-wing unmanned air vehicle for aquatic weed surveillance and management. J. Intell. Robot. Syst. Theory Appl. 2010, 57, 467. [Google Scholar] [CrossRef]
  54. Dunford, R.; Michel, K.; Gagnage, M.; Piégay, H.; Trémelo, M.L. Potential and constraints of Unmanned Aerial Vehicle technology for the characterization of Mediterranean riparian forest. Int. J. Remote Sens. 2009, 30, 4915–4935. [Google Scholar] [CrossRef]
  55. Poppe, M.; Kail, J.; Aroviita, J.; Stelmaszczyk, M.; Giełczewski, M.; Muhar, S. Assessing restoration effects on hydromorphology in European mid-sized rivers by key hydromorphological parameters. Hydrobiologia 2016, 769, 21–40. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram showing the location of the selected study sites within each Geographical Intercalibration Group (GIG) and detailed imagery of the selected reaches. The maps of Spain and UK show the delineation of the main river basins with those basins containing the study sites highlighted in red.
Figure 1. Schematic diagram showing the location of the selected study sites within each Geographical Intercalibration Group (GIG) and detailed imagery of the selected reaches. The maps of Spain and UK show the delineation of the main river basins with those basins containing the study sites highlighted in red.
Sensors 17 02210 g001
Figure 2. Example of classification outputs obtained for each Geographical Intercalibration Group (GIG). From left to right, orthoimage, 2 m × 2 m ground truth grid and classified outputs from the Artificial Neural Network (ANN); (ac) Outputs for the Central-Baltic GIG reach; (df) Outputs for the Mediterranean GIG reach; (gi) Outputs for the Very Large Rivers GIG reach.
Figure 2. Example of classification outputs obtained for each Geographical Intercalibration Group (GIG). From left to right, orthoimage, 2 m × 2 m ground truth grid and classified outputs from the Artificial Neural Network (ANN); (ac) Outputs for the Central-Baltic GIG reach; (df) Outputs for the Mediterranean GIG reach; (gi) Outputs for the Very Large Rivers GIG reach.
Sensors 17 02210 g002
Figure 3. Types of Ground Control Points (GCPs) and Unmanned Aerial Vehicles (UAVs) used to collect the imagery at each reach. (a) 1 m × 1 m Squared GCP used in the Central-Baltic reach; (b) 0.30 m diameter GCP used in the Mediterranean and Very Large Rivers reaches; (c) IRIS9+ UAV (3DR, Berkeley, CA, USA) used at the Mediterranean reach; (d) Falcon 8 Trinity (ASCTEC, Krailling, Germany) used at the Central-Baltic reach; (e) md4-1000 UAV (Microdrones, Inc., Kreuztal, Germany) used at the Very Large Rivers reach.
Figure 3. Types of Ground Control Points (GCPs) and Unmanned Aerial Vehicles (UAVs) used to collect the imagery at each reach. (a) 1 m × 1 m Squared GCP used in the Central-Baltic reach; (b) 0.30 m diameter GCP used in the Mediterranean and Very Large Rivers reaches; (c) IRIS9+ UAV (3DR, Berkeley, CA, USA) used at the Mediterranean reach; (d) Falcon 8 Trinity (ASCTEC, Krailling, Germany) used at the Central-Baltic reach; (e) md4-1000 UAV (Microdrones, Inc., Kreuztal, Germany) used at the Very Large Rivers reach.
Sensors 17 02210 g003
Figure 4. Workflow followed from imagery collection to multiple comparison analysis. ANN, GCP and GSD stand for artificial neural network, ground control point and ground sampling distance, respectively.
Figure 4. Workflow followed from imagery collection to multiple comparison analysis. ANN, GCP and GSD stand for artificial neural network, ground control point and ground sampling distance, respectively.
Sensors 17 02210 g004
Figure 5. Number of points of the 2 m × 2 m ground truth grid allocated to each feature for each of the Geographical Intercalibration Groups (GIGs).
Figure 5. Number of points of the 2 m × 2 m ground truth grid allocated to each feature for each of the Geographical Intercalibration Groups (GIGs).
Sensors 17 02210 g005
Table 1. Characteristics of the selected reaches within each of the Water Framework Directive Geographical Intercalibration Groups (GIGs). WFD D, WB ID, WB L, RL, RW and HMWB stand for Water Framework Directive Designation, Water Body Identification code, Water Body Length, Reach Length, mean wetted Reach Width and Heavily Modified Water Body, respectively. VLR stands for Very Large Rivers GIG. Area (m2) refers to the total area within which hydromorphological features were identified. RW was estimated based on a total of 20 width measurements randomly taken along the reach. 1 [35] and 2 [33].
Table 1. Characteristics of the selected reaches within each of the Water Framework Directive Geographical Intercalibration Groups (GIGs). WFD D, WB ID, WB L, RL, RW and HMWB stand for Water Framework Directive Designation, Water Body Identification code, Water Body Length, Reach Length, mean wetted Reach Width and Heavily Modified Water Body, respectively. VLR stands for Very Large Rivers GIG. Area (m2) refers to the total area within which hydromorphological features were identified. RW was estimated based on a total of 20 width measurements randomly taken along the reach. 1 [35] and 2 [33].
Geographical Intercalibration Group
DescriptorCentral-BalticMediterraneanVLR
RiverDeeJucarJucar
CountryUKSpainSpain
WFD DHMWB 1Natural 2Natural 2
WB IDGB111067052240 1ES080MSPF18.12 2ES080MSPF18.28 2
WB L (km)27.73 121.89 24.54 2
RL (km)1.41.20.96
RW (m)32.6211.0918.78
Area (m2)46,38530,85921,784
Table 2. Description for each of the hydromorphological features identified within the selected study sties. The features are adapted from the River Habitat Survey [11], the key method for the hydromorphological assessment of rivers in UK and for the designation of the water body [36]. This was further supported by the features used by the Spanish methodologies [37,38,39,40,41]. The water features recorded extended for over 5 m or >1% of the channel length following [11]. These features were recorded even if they were the result of an artificial structure.
Table 2. Description for each of the hydromorphological features identified within the selected study sties. The features are adapted from the River Habitat Survey [11], the key method for the hydromorphological assessment of rivers in UK and for the designation of the water body [36]. This was further supported by the features used by the Spanish methodologies [37,38,39,40,41]. The water features recorded extended for over 5 m or >1% of the channel length following [11]. These features were recorded even if they were the result of an artificial structure.
FeatureDescription
SubstrateSide BarConsolidated river bed material along the margins of a reach which is exposed at low flow.
ErosionPredominantly derived from eroding cliffs which are vertical or undercut banks, with a minimum height of 0.5 m and less than 50% vegetation cover.
WaterRiffleArea within the river channel presenting shallow and fast-flowing water. Generally over gravel, pebble or cobble substrate with disturbed (rippled) water surface.
Deep Water (Glides and Pools)Deep glides: deep homogeneous areas with visible flow movement along the surface.
Pools: localized deeper parts of the channel created by scouring.
Both features present fine substrate, non-turbulent and slow flow.
Shallow WaterIncludes any slow flowing and non-turbulent areas.
ChuteLow curving fall in contact with substrate.
Major impacts (pollution)Indicators of water quality pollution (e.g., accumulation of white/sluggish foam, tipping, litter, sewage, abstraction).
VegetationTreeTrees obscuring the aerial view of the river channel. The distinction between perennial and tree in dormant period was made when possible.
Vegetated Side BarSide bar presenting plant cover in more than 50% of its surface.
Vegetated BankBanks not affected by erosion. When possible the difference was made between grass and shrub cover.
Submerged Free Floating VegetationPlants rooted on the river bed with floating leaves.
Emergent Free Floating VegetationPlants rooted on the river bed with floating leaves on the water surface.
GrassPresent along the banks and floodplain as a result of intense grazing regime.
Nuisance plant specieInvasive species covering a large proportion of the banks or river channel.
ShadowsExtent of direct, overhead, tree canopy shade. Includes shading of channel and overhanging vegetation.
ArtificialAny weir, sluices, culverts, bridges, fords, deflectors or equivalent that are not underwater.
Table 3. Key characteristics of the platforms and sensors used to gather the imagery at each river reach. GIG, VLR, GCP, GSD, Mill. effect. pix. and FLA stand for Geographical Intercalibration Group, Very Large Rivers, Ground Control Point, Ground Sampling Distance, Million Effective Pixels and Focal Length Applied, respectively. The time required for the photogrammetric process (PT) is estimated based on the performance of a computer with an Intel Core i7-5820k 3.30 GHz processor, 32 Gb RAM and 2 graphic cards (Geoforce GTX 980 and Qadro K2200, NVIDIA, Santa Clara, CA, USA). 1 Sony Corporation, Tokio, Japan. 2 CanonTM, Tokio, Japan.
Table 3. Key characteristics of the platforms and sensors used to gather the imagery at each river reach. GIG, VLR, GCP, GSD, Mill. effect. pix. and FLA stand for Geographical Intercalibration Group, Very Large Rivers, Ground Control Point, Ground Sampling Distance, Million Effective Pixels and Focal Length Applied, respectively. The time required for the photogrammetric process (PT) is estimated based on the performance of a computer with an Intel Core i7-5820k 3.30 GHz processor, 32 Gb RAM and 2 graphic cards (Geoforce GTX 980 and Qadro K2200, NVIDIA, Santa Clara, CA, USA). 1 Sony Corporation, Tokio, Japan. 2 CanonTM, Tokio, Japan.
GIGCentral-BalticMediterraneanVLR
GCPs60208
GSD2.52.172.21
Flight altitude10077.6120
No. Flights422
PlatformFalcon 8 TrinityIRIS9+md4-1000
CameraSony Alpha 6000 1Canon IXUS 115 HS 2Sony Alpha ILCE-5100 1
Sensor typeCMOS APS-C type ExmorTM HD 1BCI-CMOS 2CMOS APS-C type ExmorTM 1
Mill. effect. pix.24.312.124.3
Pixel size (mm)0.003910.021690.02214
FLA (mm)20520
PT (h)121212
Table 4. Description of weather and reach characteristics during the flight at each of the Geographical Intercalibration Groups (GIGs) study sites. Weather conditions at the Central-Baltic reach were estimated based on the Shawbury (Shropshire, UK) meteorological aerodrome report. Similar information was obtained from the Agroclimatic Information for Irrigation Service weather stations at Motilleja (Albacete, Spain) and Xátiva (Valencia, Spain) for the Mediterranean and Very Large Rivers (VLR) GIGs reaches, respectively. Q stands for flow. 1 [43], 2 [44], 3 [45].
Table 4. Description of weather and reach characteristics during the flight at each of the Geographical Intercalibration Groups (GIGs) study sites. Weather conditions at the Central-Baltic reach were estimated based on the Shawbury (Shropshire, UK) meteorological aerodrome report. Similar information was obtained from the Agroclimatic Information for Irrigation Service weather stations at Motilleja (Albacete, Spain) and Xátiva (Valencia, Spain) for the Mediterranean and Very Large Rivers (VLR) GIGs reaches, respectively. Q stands for flow. 1 [43], 2 [44], 3 [45].
GIGCentral-Baltic 1Mediterranean 2,3VLR 2,3
Date21 Apirl 201528 January 201624 November 2016
Discharge (m3 s−1)4.82.53.4
Percentile Q (m3 s−1)Q80Q80Q80
Surface wind1 m s−1–3 m s−10.46 m s−10.93 m s−1
Wind direction60–350°307°293.5°
Table 5. Parameters describing the coregistration errors and the overall performance in hydromorphological feature identification for each of the Geographical Intercalibration Group (GIG) sites. N and AC stand for the number of points in the 2 m × 2 m grid and the accuracy in feature classification. GCP and VLR stand for Ground Control Point and Very Large Rivers, respectively.
Table 5. Parameters describing the coregistration errors and the overall performance in hydromorphological feature identification for each of the Geographical Intercalibration Group (GIG) sites. N and AC stand for the number of points in the 2 m × 2 m grid and the accuracy in feature classification. GCP and VLR stand for Ground Control Point and Very Large Rivers, respectively.
GIGCentral-BalticMediterraneanVLR
Total GCP error in X (cm)1.11.062.65
Total GCP error in Y (cm)1.01.492.52
Total GCP error in Z (cm)1.61.421.01
N13,08577164915
AC (%)817150
Table 6. Summary of the ANN performance in hydromorphological feature identification per Geographical Intercalibration Group (GIG) and feature. TPR, TNR, FPR and FNR stand for true positive ratio, true negative ratio, false positive ratio and false negative ratio, respectively.
Table 6. Summary of the ANN performance in hydromorphological feature identification per Geographical Intercalibration Group (GIG) and feature. TPR, TNR, FPR and FNR stand for true positive ratio, true negative ratio, false positive ratio and false negative ratio, respectively.
FeatureTPRTNRFPRFNR
Central-Baltic
Side bar0.8220.7650.0000.178
Erosion0.0770.7860.0010.923
Riffle0.8140.7560.0600.074
Deep water0.9260.7410.0080.074
Shallow water0.5880.8150.0510.412
Shadow0.8180.7700.0730.182
Vegetation0.8100.7580.0810.192
Mediterranean
Side bar0.7580.7060.0000.241
Riffle0.7360.7070.0140.263
Deep water0.5500.7240.0440.449
Shallow water0.5150.7530.0930.484
Vegetation0.7850.5650.2990.214
Pollution0.5000.7080.0010.500
Structure0.0000.7090.0001.000
Chute0.0000.7080.0001.000
Very Large Rivers
Side bar0.0000.5240.0021.000
Riffle0.0000.5270.0001.000
Deep water0.6650.4880.0340.334
Shallow water0.5550.4750.1360.444
Shadow0.5310.5000.0730.468
Vegetation0.7430.4810.3640.256
Structure0.0000.5030.0001.000
Chute0.2830.5050.0030.716
Table 7. Confusion matrix obtained per Geographical Intercalibration Group (GIG) site and feature considered. The features identified are described in Table 4 and are as follows: Side Bars (SB), Riffles (RI), Erosion (ER), Deep Water (DW), Shallow Water (SW), Chute (CH), Shadow (SH), Vegetation (VG), Pollution (PL), Structure (ST) and Georeferencing Error (GE). ANN refers to the features identified with the Artificial Neural Network algorithms. Visual refers to the features identified through visual observation and considered to be the ground truth data set.
Table 7. Confusion matrix obtained per Geographical Intercalibration Group (GIG) site and feature considered. The features identified are described in Table 4 and are as follows: Side Bars (SB), Riffles (RI), Erosion (ER), Deep Water (DW), Shallow Water (SW), Chute (CH), Shadow (SH), Vegetation (VG), Pollution (PL), Structure (ST) and Georeferencing Error (GE). ANN refers to the features identified with the Artificial Neural Network algorithms. Visual refers to the features identified through visual observation and considered to be the ground truth data set.
ANNSBRIERDWSWCHSHVGPLSTGETotal
Visual
Central-Baltic
SB10978--2-10214--31334
RI-27171-318-21976--83339
ER-132213-10238---287
DW-60-192754-829--42082
SW-262-801514-493217--72573
CH------------
SH--4-5-18031---220
VG-245117156-1972129--93250
PL------------
ST------------
Total109733053820152052-11173430--3113,085
Meditteranean
SB17619--16-----21232
RI-198--1--65--5269
ER------------
DW---38577--188--47697
SW-25-59752--541--711448
CH----2--1---3
SH------------
VG-61-248495--391914-2534990
PL-------58-316
ST-------11-810
Total17630306921343--472023-4087665
Very Large River
SB----49-5150--22226
RI---104751144--5230
ER------------
DW---2684-12119--25428
SW7--39429137597--1161811
CH---12-15422--50103
SH---136-214152--402805
VG6--3527651711431--5512475
PL------------
ST------35--19
Total13--4231382305472520--11726087

Share and Cite

MDPI and ACS Style

Rivas Casado, M.; González, R.B.; Ortega, J.F.; Leinster, P.; Wright, R. Towards a Transferable UAV-Based Framework for River Hydromorphological Characterization. Sensors 2017, 17, 2210. https://doi.org/10.3390/s17102210

AMA Style

Rivas Casado M, González RB, Ortega JF, Leinster P, Wright R. Towards a Transferable UAV-Based Framework for River Hydromorphological Characterization. Sensors. 2017; 17(10):2210. https://doi.org/10.3390/s17102210

Chicago/Turabian Style

Rivas Casado, Mónica, Rocío Ballesteros González, José Fernando Ortega, Paul Leinster, and Ros Wright. 2017. "Towards a Transferable UAV-Based Framework for River Hydromorphological Characterization" Sensors 17, no. 10: 2210. https://doi.org/10.3390/s17102210

APA Style

Rivas Casado, M., González, R. B., Ortega, J. F., Leinster, P., & Wright, R. (2017). Towards a Transferable UAV-Based Framework for River Hydromorphological Characterization. Sensors, 17(10), 2210. https://doi.org/10.3390/s17102210

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