River Cross-section Research Articles

Transport of (Micro)plastic Within a River Cross-Section—Spatio-Temporal Variations and Loads

Despite substantial research, the spatio-temporal dynamics of microplastic fluxes remain underexplored, especially in lower-order rivers. This study aims to quantify microplastic loads using a spatio-temporal sampling approach in a single cross-section of the Lahn River, a typical low-mountain river in Central Germany, over a sampling period from July 2020 to April 2021, covering varying discharge conditions, from low to high flow. A total of 198 plastic particles were detected, averaging 3.67 particles per hour, with a mean microplastic load of 0.03 ± 0.027 particles per cubic metre. Microplastic abundance varied spatially within the river cross-section, with lower concentrations found at deeper sampling positions. The data indicate that higher discharge conditions correlate with increased microplastic loads, predominantly at the water surface, suggesting that hydrological conditions significantly influence plastic transport dynamics. However, it remains unclear whether the microplastics observed at higher discharges originate from additional sources or are reactivated from river sediments. This research highlights the need for further studies to validate model assumptions and better understand the reactivation and transport mechanisms of microplastics in river systems.

River suspended-sand flux computation with uncertainty estimation using water samples and high-resolution ADCP measurements

Abstract. Measuring suspended-sand fluxes in rivers remains a scientific challenge due to their high spatial and temporal variability. To capture the vertical and lateral gradients of concentration in the cross-section, measurements with point samples are performed. However, the uncertainty related to these measurements is rarely evaluated, as few studies of the major sources of error exist. Therefore, the aim of this study is to develop a method to determine the cross-sectional sand flux and estimate its uncertainty. This SDC (for sand discharge computing) method combines suspended-sand concentrations from point samples with ADCP (acoustic Doppler current profiler) high-resolution depth and velocity measurements. The MAP (for multitransect averaged profile) method allows obtaining an average of several ADCP transects on a regular grid, including the unmeasured areas. The suspended-sand concentrations are integrated vertically by fitting a theoretical exponential suspended-sand profile to the data using Bayesian modeling. The lateral integration is based on the water depth as a proxy for the local bed shear stress to evaluate the bed concentration and sediment diffusion along the river cross-section. The estimation of uncertainty combines ISO standards and semi-empirical methods with a Bayesian approach to estimate the uncertainty due to the vertical integration. The new method is applied to data collected in four rivers under various hydro-sedimentary conditions: the Colorado, Rhône, Isère, and Amazon rivers, with computed flux uncertainties ranging between 18 % and 32 %. The relative difference between the suspended-sand flux in 21 cases calculated with the proposed SDC method compared to the International Organization for Standardization (ISO) 4363 standard method ranges between −40 % and +23 %. This method that comes with a flexible, open-source code is the first to propose an applicable uncertainty estimation that could be adapted to other flux computation methods.

Avoiding Flood by Improving Cross-Sectional Capacity through River Normalization

Bekasi is one of the regions in Indonesia that often suffer from flooding. To overcome flooding problems in the Bekasi area, the Ciliwung Cisadane River Basin Center has begun a normalization project along the Bekasi River. This research aimed to evaluate the Bekasi River's cross-sectional capacity both before and after normalization. The analysis focused on the section where the Cileungsi and Cikeas rivers converge to the Bekasi weir. In this research, secondary data were employed, such as national digital elevation model data, detailed engineering designs of Bekasi River flood control activities, and rainfall data from the Cibinong Station that covers the period from 2006 to 2020. The result from planned flood discharge with the Nakayasu synthetic unit hydrograph method for a return period of 25 years on the Bekasi River was obtained at 665.81 m3/s. The results of hydraulic analysis before normalization with the HEC-RAS application show that flooding occurred at 102 out of 116 stations, with capacities of river cross-sections ranging from 453.49 m³/s to 665.71 m³/s at these overflow stations. Following the river normalization process, the cross-sectional capacity at all Bekasi River stations can accommodate flood discharge without any instances of overflow.Keywords: Flood, HEC-RAS, Nakayasu, Normalization, Rainfall.

Hydrological model calibration in data-deficient basins using satellite altimetry and a hydrodynamic model

Large basins with insufficient hydrological station distribution and lacking sufficient meteorological data, such as water level and discharge, pose significant challenges in developing reliable hydrological models. These models are crucial for addressing water resource-related challenges such as climate change (e.g., floods and droughts) and human-induced activities (e.g., dams) impacts. Recent advancements in satellite altimetry offer an efficient alternative by providing water level data; however, integrating this data into hydrological models remains an ongoing challenge. Here we demonstrate that a properly developed hydrodynamic model can effectively translate the satellite-observed water levels at various river cross-sections (hereafter referred to as virtual stations (VSs)) into the required discharge data, facilitating reliable hydrological modeling for data-deficient basins. Our results, based on the assessment of water level data derived from a hydrodynamic model and satellite altimetry, showed high consistency (NSE > 0.96) in the Lancang-Mekong River (LMR) basin as a case study. We found that the distribution of stations along the mainstream, in terms of their number and distance from each other, can significantly impact the performance of hydrological modeling processes. The results revealed that when a large stretch of the river is not monitored by a station, a large uncertainty occurs in runoff generation, at least for that segment of the river. However, the inclusion of virtual stations in the modeling process improved the overall performance, indicating the advantages of satellite altimetry data in better modeling of data-deficient basins. Our findings revealed that virtual stations can be effectively utilized in physically based and spatially distributed hydrological models to develop a sub-basin modeling approach that addresses the regional contribution of sub-basins to downstream river flow. The integrated modeling framework holds potential for hydrological analysis and water resources management for basins lacking sufficient gauging data and provides insight into how the distribution of stations can influence the calibration of hydrological models.

Development and evaluation of hydraulic geometry relationships for large Brazilian river basins

ABSTRACT River cross-section characteristics are not available for many regions of the world, and they are generally expressed as hydraulic geometry (HG) relationships. These relationships usually represent average characteristics of a region, which can generate biased estimates in smaller basins. Therefore, to improve the representation of river channel geomorphological characteristics, the objective of this study was to use a large database of river cross-section geometries to refine HG relationships for five river basins in Brazil. HG relationships were developed for the entire basins and for hydraulically homogeneous regions, obtained through a geographical convenience method. The relationship between the coefficients of the geometric relationships and morphoclimatic characteristics was assessed, showing a trend between precipitation and slope and the coefficients of the depth relationships. The estimated bankfull width and depth values were also compared to observed data and global datasets, showing an improvement in estimates with the more refined HG relationships.

Potential Legacy of SWOT Mission for the Estimation of Flow–Duration Curves

Flow–duration curves (FDCs) provide a compact view of the historical variability of river flows, reflecting climate conditions and the main hydrologic features of river basins. The Surface Water and Ocean Topography (SWOT) satellite mission will enable the estimation of river flows globally, by sensing rivers wider than 100 m with a sampling recurrence from 3 to 21 days. This study investigated the lifetime mission potential for FDC estimation through the comparison between remotely-sensed and empirical FDCs. We employed the Global Runoff Data Center dataset and derived SWOT-like river flows by selecting gauging stations of rivers wider than 100 m with more than 10-year long daily river flow time series. Overall, 1200 gauged river cross-sections were examined. For each site, we created a set of 24 SWOT-simulated FDCs (i.e., based on different sampling recurrences, mean biases, and random errors) to be compared against their empirical counterparts through the Nash–Sutcliffe efficiency and the mean relative error. Our results show that climate and the sampling recurrence play a key role on the performance of SWOT-based FDCs. Tropical and temperate climates performed the best, whereas arid climates mostly revealed higher uncertainties, especially for high- and low-flows.

Hydroclimatic drivers of at‐a‐station hydraulic geometry of Brazilian rivers

The geometry of river channels defines their discharge conveyance capacity and can alter the flood regime. While most rivers follow the Theory of Hydraulic Geometry at-a-station power laws, the spatial patterns of exponents, coefficients, and hydrological controls beyond the scale of river reaches remain unclear. Here, we analyze 290 Brazilian river cross-sections to test the hypothesis that the hydraulic geometry is related to basin hydroclimatic characteristics, with a significant difference between humid and dry basins. We classify the basins draining to each cross-section according to the aridity index and the baseflow index, and analyze the controls on the exponents and coefficients (i.e., parameters) of the at-a-station hydraulic geometry power-laws. Our results show that the exponents and coefficients follow a spatially coherent pattern. The exponent of the width of AGH power law is related to the aridity of the basin, being larger in drier basins with a low baseflow index. The width and velocity coefficients are strongly dependent on the basin area. Whereas the velocity coefficient decreases in large basins, the width coefficient increases in large basins, and the spatial distribution of width exponent is related to the bankfull cross-section shape. The combined analysis of hydraulic geometry and basin characteristics is a promising tool for discharge estimation in ungauged basins, flood management, and maintenance of natural rivers.

An R-based integrated method for producing river bathymetry and cross-sections from recreational-grade sonar sensor data

Water bodies’ bathymetry is a crucial information for understanding and sustainably managing water resources. Bathymetric surveys can be expensive due to sonar equipment cost, but low-cost alternatives options exist. We present a methodology that standardize the bathymetric data collection and processing of recreational-grade sonar data. The sonar data postprocessing if fully implemented in R, with ready to use functions able to produce bathymetric maps or extract river cross sections’ metrics with minimal computing efforts. The method robustly produces a variety of outputs; the performance of the equipment adopted and of the interpolation technique allow for high accuracy and low-cost bathymetric reconstruction.•The method implemented allow for a robust and consistent processing of recreational-grade sonar water depth measures.•Through R-based functions the data are postprocessed to obtain bathymetry maps also for complex shape waterbodies.•Further metrics of rivers/channel cross sections can be extracted.

UAV Mapping for Flood Routing in Steep and Densely Vegetated Areas: Insights from the Contok River Basin, Garang Watershed, Indonesia

This research utilizes photogrammetry to assess flood routing dynamics in the Contok river basin, a sub-watershed with a challenging landscape characterized by steep slopes, dense vegetation, and meandering patterns. The objectives are to assess Unmanned Aerial Vehicle (UAV) mapping accuracy, evaluate the river's capacity for design flood volumes, quantify the impact of land cover changes on surface runoff, and provide insights for early warning systems and watershed conservation strategies. The study area, encompassing the Contok River Basin, a sub-watershed of the Garang Watershed, covers 7,413 km² and includes a stream length of 5,274 meters in Semarang City, Central Java, Indonesia. This research employed image processing of aerial photographs and satellite imagery. Aerial photos captured using UAV data were utilized to derive elevation data and cross-sectional profiles of the Contok River, essential for understanding channel morphology and hydraulic characteristics. Concurrently, satellite imagery was used for land cover analysis, identifying vegetation and built-up areas that influence surface runoff dynamics. Hydrological analysis was performed to quantify discharge magnitudes, simulated against river cross-sections to evaluate flood behavior under varying scenarios. Our proposed UAV mapping provides adequate accuracy for small and local areas. Furthermore, it remains reliable for flood routing analysis. We discovered that the capacity of the Contok River channel in the downstream area allows it to convey design flood discharges up to a 50-year return period, contrary to the upstream area, it overflows. Notably, the shift from vegetated to built-up and agricultural areas significantly contributes to the 10.6% increase in surface runoff. This research highlights the role of UAV-based photogrammetry in assessing and mitigating flood hazards amidst evolving land cover patterns. It also enhances the understanding of flood dynamics and thus provides insights that will serve as a reference for flood early warning systems, flood management practices, and watershed conservation.

Kajian Hidrologi dan Analisis Kapasitas Pengaliran Penampang Sungai Way Kuripan terhadap Bencana Banjir Wilayah Bandar Lampung Berbasis Hec Ras

This study aimed to determine flood height and also the area of flood in Way Kuripan River. The first step of this study is used hydrological analysis to obtain the maximum flow rate of Way Kuripan River a return periods of 2 years, 5 years, 10 years, 25 years, 50 years, and 100 years. Furthermore, to know capacities and water level of Way Kuripan River is required the hydraulic simulation process that is made easier by using Hydrologic Engineering Center-River Analysis System (HEC-RAS) program through 8 cross-sectional points. The HEC-RAS program can help be model the river flows in existing conditions with using steady flow options and input datas in the form of maximum discharge data. The output of the modeling using the HEC-RAS program can be seen in figure and tables form that present the characteristics of the Way Kuripan river cross section. The analysis of results of the HEC RAS with simulation till discharge plan when a return period of 100 years have been flood at 3 point cross section that is Sta. 4.6 and 7. In the cross section Sta. 6 has been flood as high as 1.94 meters on the right side. The flood areas include: Pesawahan, Kota Karang, Kuripan, Gedong Pakuon, Talang, Sukarame II, Negeri Olok Gading, Sumur Putri and Perwata. Keywords: Flood, HEC RAS, Way Kuripan River.

Analisis Kapasitas Efektif Tampungan Kali Pancir

Kali Pancir is one of the rivers in Jombang Regency which often experiences flooding. The cause was heavy rain and the critical condition of the embankments. In this research, flood discharge analysis was carried out using synthetic hydrographs and hydraulic modeling for river flow. The data needed is rain data from rain stations around the Pancir River and river profiles. The method used to calculate flood discharge is the Nakayasu Synthetic Unit Hydrograph for 6 different return periods. Simulation using HEC-RAS software is used to simulate flow profiles to determine river cross-sections that are critical to flooding. Based on the results of the calculation analysis, it was found that the peak flood discharge occurred at 3.5 hours. The greater the return period, the greater the river cross section that experiences overflow. At a return period of 1.01 years, the critical cross section starts at STA 45, while at a return period of 50 years, the critical cross section occurs starting at STA 6.

KNN Local Linear Regression for Demarcating River Cross-Sections with Point Cloud Data from UAV Photogrammetry URiver-X

Aerial surveying with unmanned aerial vehicles (UAVs) has been popularly employed in river management and flood monitoring. One of the major processes in UAV aerial surveying for river applications is to demarcate the cross-section of a river. From the photo images of aerial surveying, a point cloud dataset can be abstracted with the structure from the motion technique. To accurately demarcate the cross-section from the cloud points, an appropriate delineation technique is required to reproduce the characteristics of natural and manmade channels, including abrupt changes, bumps and lined shapes. Therefore, a nonparametric estimation technique, called the K-nearest neighbor local linear regression (KLR) model, was tested in the current study to demarcate the cross-section of a river with a point cloud dataset from aerial surveying. The proposed technique was tested with synthetically simulated trapezoidal, U-shape and V-shape channels. In addition, the proposed KLR model was compared with the traditional polynomial regression model and another nonparametric technique, locally weighted scatterplot smoothing (LOWESS). The experimental study was performed with the river experiment center in Andong, South Korea. Furthermore, the KLR model was applied to two real case studies in the Migok-cheon stream on Hapcheon-gun and Pori-cheon stream on Yecheon-gun and compared to the other models. With the extensive applications to the feasible river channels, the results indicated that the proposed KLR model can be a suitable alternative for demarcating the cross-section of a river with point cloud data from UAV aerial surveying by reproducing the critical characteristics of natural and manmade channels, including abrupt changes and small bumps as well as different shapes. Finally, the limitation of the UAV-driven demarcation approach was also discussed due to the penetrability of RGB sensors to water.

The Impact of Beaver Dams on the Dynamic of Groundwater Levels at Łąki Soleckie

Areas excluded from agricultural production are susceptible to the presence of beaver families. The most significant changes occur during the initial period, when agricultural utilization is abandoned and beavers establish their presence on the land. During this period, some parcels remain uncultivated, while agricultural activities persist in neighboring areas. This situation is accompanied by the destruction of beaver dams, especially during periods of abundant water resources, and notably during intensive fieldwork. The article presents field studies aimed at determining the extent to which constructed and operational beaver dams contribute to changes in groundwater levels in drained peatland areas. In order to protect and sustainably use peat soils, it is necessary to maintain their high moisture content by ensuring a high groundwater level elevation. This can be achieved through the use of existing damming structures in the area (levees, weirs). Beaver dams can also serve a similar function, blocking the outflow of water from peat lands by raising the water level and consequently retaining it naturally. The specific objective was to develop principles for verifying factors influencing the effects of beaver dam construction on groundwater levels in fields within their range of influence. The water table levels within the study area during rainless periods were influenced by water levels in ditches, dependent on beaver activity in the nearby river. Beaver activities, manifested through dam construction, were influenced by periodic water resources in the river, defined by the cumulative monthly precipitation. Factors affecting groundwater levels in rainless periods on the plots also included the distance from the river cross-section and the permeability of soils expressed by the filtration coefficient of the active layer. Beaver dams had the greatest impact on stabilizing the water table in the soil profile closest to the river.

Efficient river hydrodynamics modelling in realistic river systems using a Fourier neural operator-based network

Effective river management relies heavily on the accurate simulation of river flow dynamics to develop scenario-based strategies and inform decision making. Deep learning (DL) models, as alternatives to conventional simulation techniques, offer a compelling blend of precision and efficiency. This study introduces FNOCL (Fourier Neural Operator with ConvLSTM units), a novel DL neural network for modelling unsteady natural river flows. The network is adept at handling the complexities of spatial and temporal variations. FNOCL builds upon the Fourier neural operator (FNO) by incorporating convolutional long short-term memory (ConvLSTM) layers, combining robust performance with high proficiency in terms of capturing spatiotemporal patterns. Notably, FNOCL accommodates complex river cross-section profiles, thus enabling the accurate modelling of flow dynamics within intricate river morphologies. Our study applies FNOCL to two distinct river networks with varying morphology complexity levels. With reference to the HEC-RAS model, FNOCL showcases its ability to attain high accuracy in model prediction tasks and superior generalization with regard to addressing the uncertainty associated with inflow variations and cross-sectional changes. In comparative analyses with benchmark models (i.e., FNO and ConvLSTM), FNOCL performs most efficiently during the training process and obtains the most accurate water level and river flow predictions in the spatiotemporal domain. Furthermore, FNOCL exhibits a substantial speedup over HEC-RAS, notably reaching up to 50 times in a complex river network. These results highlight FNOCL as an efficient surrogate for conventional simulators, facilitating precise analyses across various scenarios. It provides efficient computations for system responses within complex river flow networks, thereby enhancing applications in river management, such as flood risk assessment, flood control, and the planning and design of hydraulic control structures.

Morphological Evolution of Euphrates River Cross-Sections in Central Iraq: Response to Streamflow and Sediment Load Variability

Morphological Evolution of Euphrates River Cross-Sections in Central Iraq: Response to Streamflow and Sediment Load Variability

Overbank Flow, Sediment Transport, and Channel Morphology in the Lower Yellow River: A Review

As a prerequisite and foundation for studying the evolution mechanism of river channels, an in-depth understanding of the cross-sectional morphology adjustment is required. As a starting point, it is crucial to systematically summarize and generalize the research findings on channel morphological adjustment obtained to date, particularly in the context of the significant changes in the water and sediment conditions of large rivers that have occurred worldwide. This paper provides a comprehensive review of the research findings on the three following aspects of the Lower Yellow River: the transverse distribution of overbank flow velocity, the transverse distribution of suspended sediment concentration, and the morphological adjustment of the river cross-section. There are various equations available to predict the lateral depth–average flow velocity distribution. These equations are classified into the two following categories: empirical and theoretical formulas. Theoretical formulas are obtained through consideration of the cross-sectional morphology, accounting for inertial force terms caused by secondary flow, and momentum transfer between the main channel and its floodplain. Similarly, empirical equations and theoretical formulas for sediment concentration transverse distribution are also summarized, given the different influencing factors and assumptions. We also discuss the morphological adjustment of river cross-sections based on the analysis of measured data, mathematical model calculation, and the physical model test. In particular, we propose the idea of revealing channel cross-section morphology evolution mechanisms from the theoretical level of water and sediment movement and distribution. This review aims to enhance understanding of overbank flow, sediment transport, and channel morphology in the Lower Yellow River and may also serve to some extent as a reference for the evolution and management of channels in other rivers.

A Numerical Expedition through the Mathematical Representation of Complex Braided Morphometry—A Case Study of Brahmaputra River in India

The present work explores the process of mathematical representation for the complex geometry of a wide alluvial river with high braiding intensities. It primarily focuses on an approach to developing a numerical solution algorithm for representing the complex channel geometry of the braided Brahmaputra River. Traditional elliptic PDEs with boundary-fitted coordinate transformation were deployed, converting the non-uniform physical plane into a transformed uniform orthogonal computational plane. This study was conducted for the river channel reach with upstream and downstream nodes at Pandu and Jogighopa (reach length ~100 km), respectively, within the Assam flood plain in India, with fourteen measured river cross-sections for the year of 1997. The geo-referenced image covering the river stretch in 1997 was delineated using a ArcGIS software 9.0 tool by digitizing the bank lines. Stream bed interpolation was conducted by interpolating bed elevation from a bathymetrical database onto code-generated mesh nodes. Discretization of the domain was performed through the developed computer code, and the bed-level matrix was generated by the IDW method as well as the MATLAB tool using the nearest neighborhood technique. A mathematical representation of a digital terrain model was thus developed. This generated model was employed as a geometrical data input to simulate secondary flow utilizing 2D depth-averaged equations with the flow dispersion stress tensor as an extra source component, coming from curvilinear flow patterns caused by severe river braiding. The developed model may further be useful in mathematically representing the geometrical complexities of braided rivers with a relatively realistic assessment of the various parameters involved if deployed with improved river modeling with morphometric evolution.

The Effect of Changes in Flood Discharge on Sedimentation Rates in Middle Segment of Batang Kuranji Check Dam using HEC-RAS 6.0.0

The magnitude of the longitudinal slope of river gratly determines the magnitude of sediment transport. This is due to the large flow velocity. This shape causes gradation of the river bed, especially upstream of the river. This phenomenon is seen in the middle segment of Batang Kuranji. where a sediment control buildg or check dam has been built. However, sediment transport and changes in river cross-section upstream still occur. Based on this phenomenon, a simulation of changes in flood discharge on sediment rate was carried out. This simulation was carried out using the Meyer Peter Muller method found in HEC-RAS 6.0.0. The simulation was carried out using 25 year planned flood discharge with smaller variations in discharge changes, namely Q1, Q2 and Q3. In general, the simulation results show that changes in flood discharge determine the magnitude of changes in sediment rates, the greater the flow discharge, the higher the sediment rate, both with and without check dams. The simulation results also show that the sediment rate in conditions with a check dam has decreased compared to conditions without a check dam, this shows that the presence of a check dam can reduce the sediment rate that occurred before the check dam was built.

Microplastics in the Water Column of the Rhine River Near Basel: 22 Months of Sampling.

Measured microplastic concentrations in river surface waters fluctuate greatly. This variability is affected by season and is codriven by factors, such as sampling methodologies, sampling site, or sampling position within site. Unfortunately, most studies comprise single-instance measurements, whereas extended sampling periods are better suited to assessing the relevance of such factors. Moreover, microplastic concentrations in riverine water column remain underexplored. Similar to the oceans, however, this compartment likely holds significant amounts of microplastics. By representatively sampling the entire Rhine River cross-section near Basel through five sampling points over 22 months, we found a median microplastic (50-3000 μm) concentration of 4.48 n m-3, and estimated a widely ranging load between 4.04 × 102 n s-1 and 3.57 × 105 n s-1. We also show that the microplastic concentration in the water column was not well explained by river discharge. This suggests that although high discharge events as observed here can over short time periods lead to peak microplastic concentrations (e.g., 1.23 × 102 n m-3), microplastic load variance was not dominated by discharge in the study area.

Environmental DNA time series analysis of a temperate stream reveals distinct seasonal community and functional shifts

AbstractEnvironmental DNA (eDNA) extracted from water is routinely used in river biodiversity research, and via metabarcoding eDNA can provide comprehensive taxa lists with little effort and cost. However, eDNA‐based species detection in streams and rivers may be influenced by sampling season and other key factors such as water temperature and discharge. Research linking these factors and also informing on the potential of eDNA metabarcoding to detect shifts in ecological signatures, such as species phenology and functional feeding groups across seasons, is missing. To address this gap, we collected water samples every 2 weeks for 15 months at a long‐term ecological research (LTER) site and at three different positions in the river's cross section, specifically the water surface, riverbed, and riverbank. For these 102 samples, we analyzed macroinvertebrate species and molecular operational taxonomic unit (OTU) richness and temporal community turnover across seasons based on cytochrome c oxidase subunit I (COI) metabarcoding data. Using Generalized Additive Models, we found a significant influence of sampling season on species richness. Community turnover followed a cyclic pattern, reflecting the continuous change of the macroinvertebrate community throughout the year (“seasonal clock”). Although water temperature had no influence on the inferred species richness, higher discharge reduced the number of Annelida and Ephemeroptera species detectable with eDNA. Most macroinvertebrate taxa showed the highest species richness in spring, in particular merolimnic species with univoltine life cycles. Further, we detected an increase in the proportion of shredders in winter and parasites in summer. Our results show the usefulness of highly resolved eDNA metabarcoding time series data for ecological research and biodiversity monitoring in streams and rivers.