Inland navigation and land use interact to impact European freshwater biodiversity

Inland navigation in Europe is proposed to increase in the coming years, being promoted as a low-carbon form of transport. However, we currently lack knowledge on how this would impact biodiversity at large scales and interact with existing stressors. Here we addressed this knowledge gap by analysing fish and macroinvertebrate community time series across large European rivers comprising 19,592 observations from 4,049 sampling sites spanning the past 32 years. We found ship traffic to be associated with biodiversity declines, that is, loss of fish and macroinvertebrate taxonomic richness, diversity and trait richness. Ship traffic was also associated with increases in taxonomic evenness, which, in concert with richness decreases, was attributed to losses in rare taxa. Ship traffic was especially harmful for benthic taxa and those preferring slow flows. These effects often depended on local land use and riparian degradation. In fish, negative impacts of shipping were highest in urban and agricultural landscapes. Regarding navigation infrastructure, the negative impact of channelization on macroinvertebrates was evident only when riparian degradation was also high. Our results demonstrate the risk of increasing inland navigation on freshwater biodiversity. Integrative waterway management accounting for riparian habitats and landscape characteristics could help to mitigate these impacts.

Field-based molecular detection of Batrachochytrium dendrobatidis in critically endangered Atelopus toads and aquatic habitats in Ecuador

Batrachochytrium dendrobatidis (Bd) is a lethal fungal species that parasitizes vertebrates and is associated with the worldwide decline of amphibian populations. The development of sensitive, rapid detection methods, particularly DNA-based techniques, is critical for effective management strategies. This study evaluates the efficacy of DNA extraction and a portable PCR device in a mountable field laboratory setup for detecting Bd near the habitats of three critically endangered Atelopus toad species in Ecuador. We collected skin swabs from Atelopus balios, A. nanay, and A. bomolochos, and environmental DNA (eDNA) samples from streams in Andean and coastal regions of Ecuador. For eDNA, a comparison was made with duplicates of the samples that were processed in the field and in a standard university laboratory. Our findings revealed Bd detection in eDNA and swabs from 6 of 12 water samples and 10 of 12 amphibian swab samples. The eDNA results obtained in the field laboratory were concordant with those obtained under campus laboratory conditions. These findings highlight the potential of field DNA-based monitoring techniques for detecting Bd in amphibian populations and their aquatic habitats, particularly in remote areas. Furthermore, this research aligns with the National Action Plan for the Conservation of Ecuadorian Amphibians and contributes to the global effort to control this invasive and deadly fungus.

Modeling greenhouse gas emissions from riverine systems: A review

Despite the recognized importance of flowing waters in global greenhouse gas (GHG) budgets, riverine GHG models remain oversimplified, consequently restraining the development of effective prediction for riverine GHG emissions feedbacks. Here we elucidate the state of the art of riverine GHG models by investigating 148 models from 122 papers published from 2010 to 2021. Our findings indicate that riverine GHG models have been mostly data-driven models (83%), while mechanistic and hybrid models were uncommonly applied (12% and 5%, respectively). Overall, riverine GHG models were mainly used to explain relationships between GHG emissions and biochemical factors, while the role of hydrological, geomorphic, land use and cover factors remains missing. The development of complex and advanced models has been limited by data scarcity issues; hence, efforts should focus on developing affordable automatic monitoring methods to improve data quality and quantity. For future research, we request for basin-scale studies explaining river and land-surface interactions for which hybrid models are recommended given their flexibility. Such a holistic understanding of GHG dynamics would facilitate scaling-up efforts, thereby reducing uncertainties in global GHG estimates. Lastly, we propose an application framework for model selection based on three main criteria, including model purpose, model scale and the spatiotemporal characteristics of GHG data, by which optimal models can be applied in various study conditions.

Impact of salinity gradient, water pollution and land use types on greenhouse gas emissions from an urbanized estuary

Estuaries have been recognized as one of the major sources of greenhouse gases (GHGs) in aquatic systems; yet we still lack insights into the impact of both anthropogenic and natural factors on the dynamics of GHG emissions. Here, we assessed the spatiotemporal dynamics and underlying drivers of the GHG emissions from the Scheldt Estuary with a focus on the effects of salinity gradient, water pollution, and land use types, together with their interaction. Overall, we found a negative impact of salinity on carbon dioxide (CO2) and nitrous oxide (N2O) emissions which can be due to the decrease of both salinity and water quality when moving upstream. Stronger impact of water pollution on the GHG emissions was found at the freshwater sites upstream compared to saline sites downstream. In particular, when water quality of the sites reduced from good, mainly located in the mouth and surrounded by arable sites, to polluted, mainly located in the upstream and surrounded by urban sites, CO2 emissions from the sites doubled while N2O emissions tripled. Similarly, the effects of water pollution on methane (CH4) emissions became much stronger in the freshwater sites compared to the saline sites. These decreasing effects from upstream to the mouth were associated with the increase in urbanization as sites surrounded by urban areas released on average almost two times more CO2 and N2O than sites surrounded by nature and industry areas. Applied machine learning methods also revealed that, in addition to salinity effects, nutrient and organic enrichment stimulated the GHG emissions from the Scheldt Estuary. These findings highlight the importance of the interaction between salinity, water pollution, and land use in order to understand their influences on GHG emissions from dynamic estuarine systems.

Linking fine-scale behaviour to the hydraulic environment shows behavioural responses in riverine fish

BACKGROUND

Fish migration has severely been impacted by dam construction. Through the disruption of fish migration routes, freshwater fish communities have seen an incredible decline. Fishways, which have been constructed to mitigate the problem, have been shown to underperform. This is in part due to fish navigation still being largely misunderstood. Recent developments in tracking technology and modelling make it possible today to track (aquatic) animals at very fine spatial (down to one meter) and temporal (down to every second) scales. Hidden Markov models are appropriate models to analyse behavioural states at these fine scales. In this study we link fine-scale tracking data of barbel (Barbus barbus) and grayling (Thymallus thymallus) to a fine-scale hydrodynamic model. With a HMM we analyse the fish's behavioural switches to understand their movement and navigation behaviour near a barrier and fishway outflow in the Iller river in Southern Germany.

METHODS

Fish were tracked with acoustic telemetry as they approached a hydropower facility and were presented with a fishway. Tracking resulted in fish tracks with variable intervals between subsequent fish positions. This variability stems from both a variable interval between tag emissions and missing detections within a track. After track regularisation hidden Markov models were fitted using different parameters. The tested parameters are step length, straightness index calculated over a 3-min moving window, and straightness index calculated over a 10-min window. The best performing model (based on a selection by AIC) was then expanded by allowing flow velocity and spatial velocity gradient to affect the transition matrix between behavioural states.

RESULTS

In this study it was found that using step length to identify behavioural states with hidden Markov models underperformed when compared to models constructed using straightness index. Of the two different straightness indices assessed, the index calculated over a 10-min moving window performed better. Linking behavioural states to the ecohydraulic environment showed an effect of the spatial velocity gradient on behavioural switches. On the contrary, flow velocity did not show an effect on the behavioural transition matrix.

CONCLUSIONS

We found that behavioural switches were affected by the spatial velocity gradient caused by the attraction flow coming from the fishway. Insight into fish navigation and fish reactions to the ecohydraulic environment can aid in the construction of fishways and improve overall fishway efficiencies, thereby helping to mitigate the effects migration barriers have on the aquatic ecosystem.

The recovery of European freshwater biodiversity has come to a halt

Owing to a long history of anthropogenic pressures, freshwater ecosystems are among the most vulnerable to biodiversity loss1. Mitigation measures, including wastewater treatment and hydromorphological restoration, have aimed to improve environmental quality and foster the recovery of freshwater biodiversity2. Here, using 1,816 time series of freshwater invertebrate communities collected across 22 European countries between 1968 and 2020, we quantified temporal trends in taxonomic and functional diversity and their responses to environmental pressures and gradients. We observed overall increases in taxon richness (0.73% per year), functional richness (2.4% per year) and abundance (1.17% per year). However, these increases primarily occurred before the 2010s, and have since plateaued. Freshwater communities downstream of dams, urban areas and cropland were less likely to experience recovery. Communities at sites with faster rates of warming had fewer gains in taxon richness, functional richness and abundance. Although biodiversity gains in the 1990s and 2000s probably reflect the effectiveness of water-quality improvements and restoration projects, the decelerating trajectory in the 2010s suggests that the current measures offer diminishing returns. Given new and persistent pressures on freshwater ecosystems, including emerging pollutants, climate change and the spread of invasive species, we call for additional mitigation to revive the recovery of freshwater biodiversity.

Epigenetic analytical approaches in ecotoxicological aquatic research

Environmental epigenetics has become a key research focus in global climate change studies and environmental pollutant investigations impacting aquatic ecosystems. Specifically, triggered by environmental stress conditions, intergenerational DNA methylation changes contribute to biological adaptive responses and survival of organisms to increase their tolerance towards these conditions. To critically review epigenetic analytical approaches in ecotoxicological aquatic research, we evaluated 78 publications reported over the past five years (2016-2021) that applied these methods to investigate the responses of aquatic organisms to environmental changes and pollution. The results show that DNA methylation appears to be the most robust epigenetic regulatory mark studied in aquatic animals. As such, multiple DNA methylation analysis methods have been developed in aquatic organisms, including enzyme restriction digestion-based and methyl-specific immunoprecipitation methods, and bisulfite (in)dependent sequencing strategies. In contrast, only a handful of aquatic studies, i.e. about 15%, have been focusing on histone variants and post-translational modifications due to the lack of species-specific affinity based immunological reagents, such as specific antibodies for chromatin immunoprecipitation applications. Similarly, ncRNA regulation remains as the least popular method used in the field of environmental epigenetics. Insights into the opportunities and challenges of the DNA methylation and histone variant analysis methods as well as decreasing costs of next generation sequencing approaches suggest that large-scale epigenetic environmental studies in model and non-model organisms will soon become available in the near future. Moreover, antibody-dependent and independent methods, such as mass spectrometry-based methods, can be used as an alternative epigenetic approach to characterize global changes of chromatin histone modifications in future aquatic research. Finally, a systematic guide for DNA methylation and histone variant methods is offered for ecotoxicological aquatic researchers to select the most relevant epigenetic analytical approach in their research.

Freshwater Management Discourses in the Northern Peruvian Andes: The Watershed-Scale Complexity for Integrating Mining, Rural, and Urban Stakeholders

The Peruvian environmental action plan seeks headwaters protection as one of its integrated watershed management objectives. However, heterogeneous social and environmental conditions shape this freshwater management challenge at subnational scales. We have noticed different interpretations of this challenge. To map the debate, understand the diverse interpretations, and frame political choices, we conducted semi-structured interviews with institutional and non-institutional stakeholders for performing discourse analysis in an Andean watershed where mountaintop gold mining, midstream farmers, and the downstream Cajamarca city coexist. One discourse dominates the debate on protecting the freshwater supply and argues the importance of river impoundment, municipal storage capacity, and institutional leadership. The other two discourses revolve around protecting the mountain aquifer. The second discourse does so with a fatalistic view of headwaters protection and rural support. The third discourse partially shifts the debate towards the need for improving rural capacity building and (ground)water inventories. To understand evolutions in society, it is crucial to understand these three discourses, including the types of knowledge that actors present as legitimate, the attributed roles to all stakeholders, and the kinds of worldviews informing each discourse. The interaction among discourses could hinder integrated watershed management at worst or, at best, help inspire multi-stakeholder collaboration.

Widespread triazole pesticide use affects infection dynamics of a global amphibian pathogen

The sixth mass extinction is a consequence of complex interplay between multiple stressors with negative impact on biodiversity. We here examine the interaction between two globally widespread anthropogenic drivers of amphibian declines: the fungal disease chytridiomycosis and antifungal use in agriculture. Field monitoring of 26 amphibian ponds in an agricultural landscape shows widespread occurrence of triazole fungicides in the water column throughout the amphibian breeding season, together with a negative correlation between early season application of epoxiconazole and the prevalence of chytrid infections in aquatic newts. While triazole concentrations in the ponds remained below those that inhibit growth of Batrachochytrium dendrobatidis, they bioaccumulated in the newts' skin up to tenfold, resulting in cutaneous growth-suppressing concentrations. As such, a concentration of epoxiconazole, 10 times below that needed to inhibit fungal growth, prevented chytrid infection in anuran tadpoles. The widespread presence of triazoles may thus alter chytrid dynamics in agricultural landscapes.

Invasion impacts and dynamics of a European-wide introduced species

Globalization has led to the introduction of thousands of alien species worldwide. With growing impacts by invasive species, understanding the invasion process remains critical for predicting adverse effects and informing efficient management. Theoretically, invasion dynamics have been assumed to follow an "invasion curve" (S-shaped curve of available area invaded over time), but this dynamic has lacked empirical testing using large-scale data and neglects to consider invader abundances. We propose an "impact curve" describing the impacts generated by invasive species over time based on cumulative abundances. To test this curve's large-scale applicability, we used the data-rich New Zealand mud snail Potamopyrgus antipodarum, one of the most damaging freshwater invaders that has invaded almost all of Europe. Using long-term (1979-2020) abundance and environmental data collected across 306 European sites, we observed that P. antipodarum abundance generally increased through time, with slower population growth at higher latitudes and with lower runoff depth. Fifty-nine percent of these populations followed the impact curve, characterized by first occurrence, exponential growth, then long-term saturation. This behaviour is consistent with boom-bust dynamics, as saturation occurs due to a rapid decline in abundance over time. Across sites, we estimated that impact peaked approximately two decades after first detection, but the rate of progression along the invasion process was influenced by local abiotic conditions. The S-shaped impact curve may be common among many invasive species that undergo complex invasion dynamics. This provides a potentially unifying approach to advance understanding of large-scale invasion dynamics and could inform timely management actions to mitigate impacts on ecosystems and economies.

Sampling errors and variability in video transects for assessment of reef fish assemblage structure and diversity

Video monitoring is a rapidly evolving tool in aquatic ecological research because of its non-destructive ability to assess fish assemblages. Nevertheless, methodological considerations of video monitoring techniques are often overlooked, especially in more complex sampling designs, causing inefficient data collection, processing, and interpretation. In this study, we discuss how video transect sampling designs could be assessed and how the inter-observer variability, design errors and sampling variability should be quantified and accounted for. The study took place in the coastal areas of the Galapagos archipelago and consisted of a hierarchical repeated-observations sampling design with multiple observers. Although observer bias was negligible for the assessment of fish assemblage structure, diversity and counts of individual species, sampling variability caused by simple counting/detection errors, observer effects and instantaneous fish displacement was often important. Especially for the counts of individual species, sampling variability most often exceeded the variability of the transects and sites. An extensive part of the variability in the fish assemblage structure was explained by the different transects (13%), suggesting that a sufficiently high number of transects is required to account for the within-location variability. Longer transect lengths allowed a better representation of the fish assemblages as sampling variability decreased by 33% if transect length was increased from 10 to 50 meters. However, to increase precision, including more repeats was typically more efficient than using longer transect lengths. The results confirm the suitability of the technique to study reef fish assemblages, but also highlight the importance of a sound methodological assessment since different biological responses and sampling designs are associated with different levels of sampling variability, precision and ecological relevance. Therefore, besides the direct usefulness of the results, the procedures to establish them may be just as valuable for researchers aiming to optimize their own sampling technique and design.

Impact of heavy metal exposure on biological control of a deadly amphibian pathogen by zooplankton

Despite devastating effects on global biodiversity, efficient mitigation strategies against amphibian chytridiomycosis are lacking. Since the free-living pathogenic zoospores of Batrachochytrium dendrobatidis (Bd), the infective stage of this disease, can serve as a nutritious food source for components of zooplankton communities, these groups may act as biological control agents by eliminating zoospores from the aquatic environment. Such pathogen-predator interaction is, however, embedded in the aquatic food web structure and is therefore affected by abiotic factors interfering with these networks. Heavy metals, released from both natural and anthropogenic sources, are widespread contaminants of aquatic ecosystems and may interfere with planktonic communities and thus pathogen elimination rates. We investigated the interaction between zooplankton communities and chytridiomycosis infections in a Flemish agricultural region. Moreover, we also investigated the impact of heavy metal contamination, that was previously investigated in the region and presented in recent work, on zooplankton assemblages and chytridiomycosis infections. Finally, we tested the effect of sublethal concentrations of copper and zinc on Bd removal rates by Daphnia magna in a laboratory assay. Although zinc, copper, nickel and chromium were widely abundant pollutants, heavy metals were no driving force for zooplankton assemblages at our study locations. Moreover, our field survey did not reveal indirect effects of zooplankton assemblages on chytridiomycosis infections. However, sampling occasions testing negative for Bd showed a higher degree of copper contamination compared to positive sampling occasions, indicating a potential inhibitory effect of copper on Bd prevalence. Finally, whereas D. magna significantly reduced zoospore densities in its environment, sublethal concentrations of copper and zinc showed no interference with pathogen removal in the laboratory assay. Our results provide perspectives for further research on such a biological control strategy against chytridiomycosis by optimizing environmental conditions for pathogen predation.

Greenhouse gas dynamics in an urbanized river system: influence of water quality and land use

Rivers act as a natural source of greenhouse gases (GHGs). However, anthropogenic activities can largely alter the chemical composition and microbial communities of rivers, consequently affecting their GHG production. To investigate these impacts, we assessed the accumulation of CO₂, CH₄, and N₂O in an urban river system (Cuenca, Ecuador). High variation of dissolved GHG concentrations was found among river tributaries that mainly depended on water quality and land use. By using Prati and Oregon water quality indices, we observed a clear pattern between water quality and the dissolved GHG concentration: the more polluted the sites were, the higher were their dissolved GHG concentrations. When river water quality deteriorated from acceptable to very heavily polluted, the mean value of pCO₂ and dissolved CH₄ increased by up to ten times while N₂O concentrations boosted by 15 times. Furthermore, surrounding land-use types, i.e., urban, roads, and agriculture, could considerably affect the GHG production in the rivers. Particularly, the average pCO₂ and dissolved N₂O of the sites close to urban areas were almost four times higher than those of the natural sites while this ratio was 25 times in case of CH₄, reflecting the finding that urban areas had the worst water quality with almost 70% of their sites being polluted while this proportion of nature areas was only 12.5%. Lastly, we identified dissolved oxygen, ammonium, and flow characteristics as the main important factors to the GHG production by applying statistical analysis and random forests. These results highlighted the impacts of land-use types on the production of GHGs in rivers contaminated by sewage discharges and surface runoff.

Agricultural contaminants in amphibian breeding ponds: Occurrence, risk and correlation with agricultural land use

Anthropogenic pressure such as agricultural pollution globally affects amphibian populations. In this study, a total of 178 different compounds from five agrochemical groups (i.e. antimicrobial drugs residues (ADRs), coccidiostats and anthelmintics, heavy metals, mycotoxins and pesticides) were determined monthly, from March until June 2019 in 26 amphibian breeding ponds in Flanders, Belgium. Furthermore, a possible correlation between the number and concentration of selected contaminants that were found and the percentage of arable land within a 200 m radius was studied. Within each group, the highest detected concentrations were obtained for 4-epioxytetracycline (0.422 μg L^-1), levamisole (0.550 μg L^-1), zinc (333.1 μg L^-1), 3-acetyldeoxynivalenol (0.013 μg L^-1), and terbuthylazine (38.7 μg L^-1), respectively, with detection frequencies ranging from 1 (i.e. 3-acetyldeoxynivalenol) to 26 (i.e. zinc) out of 26 ponds. Based on reported acute and chronic ecotoxicological endpoints, detected concentrations of bifenthrin, cadmium, copper, cypermethrin, hexachlorobenzene, mercury, terbuthylazine, and zinc pose a substantial ecological risk to aquatic invertebrates such as Daphnia magna and Ceriodaphnia dubia, which both play a role in the food web and potentially in amphibian disease dynamics. Additionally, the detected concentrations of copper were high enough to exert chronic toxicity in the gray treefrog (Hyla versicolor). The number of detected compounds per pond ranged between 0 and 5 (ADRs), 0 - 2 (coccidiostats and anthelmintics), 1 - 7 (heavy metals), 0 - 4 (mycotoxins), and 0 - 12 (pesticides) across the four months. Furthermore, no significant correlation was demonstrated between the number of detected compounds per pond, as well as the detected concentrations of 4-epioxytetracycline, levamisole, copper, zinc, enniatin B and terbuthylazine, and the percentage of arable land within a 200 m radius. For heavy metals and pesticides, the number of compounds per pond varied significantly between months. Conclusively, amphibian breeding ponds in Flanders were frequently contaminated with agrochemicals, yielding concentrations up to the high μg per liter level, regardless of the percentage surrounding arable land, however showing temporal variation for heavy metals and pesticides. This research also identifies potential hazardous substances which may be added to the European watch list (CD 2018/408/EC) in the future.

A global agenda for advancing freshwater biodiversity research

Global freshwater biodiversity is declining dramatically, and meeting the challenges of this crisis requires bold goals and the mobilisation of substantial resources. While the reasons are varied, investments in both research and conservation of freshwater biodiversity lag far behind those in the terrestrial and marine realms. Inspired by a global consultation, we identify 15 pressing priority needs, grouped into five research areas, in an effort to support informed stewardship of freshwater biodiversity. The proposed agenda aims to advance freshwater biodiversity research globally as a critical step in improving coordinated actions towards its sustainable management and conservation.

Integrated mechanistic and data-driven modeling for risk assessment of greenhouse gas production in an urbanized river system

Surrounded by intense anthropogenic activities, urban polluted rivers have increasingly been reported as a significant source of greenhouse gases (GHGs). However, unlike pollution and climate change, no integrated urban water models have investigated the GHG production in urban rivers due to system complexity. In this study, we proposed a novel integrated framework of mechanistic and data-driven models to qualitatively assess the risks of GHG accumulation in an urban river system in different water management interventions. Particularly, the mechanistic model delivered elaborated insights into river states in four intervention scenarios in which the installation of a new wastewater treatment plant using two different technologies, together with new sewage systems and additional retention tanks, were assessed during dry and rainy seasons. From the insights, we applied fuzzy rule-based models as a decision support tool to predict the GHG accumulation risks and identify their driving factors in the scenarios. The obtained results indicated the important role of new discharge connection and additional storage capacity in decreasing pollutant concentrations, consequently, reducing the risks. Moreover, among the major variables explaining the GHG accumulation in the rivers, DO level was considerably affected by the reaeration capacity of the rivers that was strongly dependent on river slope and flow. Furthermore, river water quality emerged as the most critical variable explaining the pCO₂ and N₂O accumulation that implied that the more polluted and anaerobic the sites were, the higher were their GHG accumulation. Given its simplicity and transparency, the proposed modeling framework can be applied to other river basins as a decision support tool in setting up integrated urban water management plans.

Vertical flow constructed wetlands using expanded clay and biochar for wastewater remediation: A comparative study and prediction of effluents using machine learning

This study evaluated and compared the performance of two vertical flow constructed wetlands (VF) using expanded clay (VF₁) and biochar (VF₂), of which both are low-cost, eco-friendly, and exhibit potentially high adsorption as compared to conventional filter layers. Both VFs achieved relatively high removal for organic matters (i.e. Biological oxygen demand during 5 days, BOD₅) and nitrogen, accounting for 9.5 - 10.5 g^.BOD₅^.m^-2.d^-1 and 3.5 - 3.6 g^.NH₄-N^.m^-2.d^-1, respectively. The different filter materials did not exert any significant discrepancy to effluent quality in terms of suspended solids, organic matters and NO₃-N (P > 0.05), but they did influence NH₄-N effluent as evidenced by the removal rate of that by VF₁ and VF₂ being of 8₂.4 ± 5.7 and 84.6 ± 6.4%, respectively (P < 0.05). The results obtained from the designed systems were further subject to machine learning to clarify the effecting factors and predict the effluents. The optimal algorithms were random forest, generalized linear model, and support vector machine. The values of the coefficient of determination (R²) and the root mean square error (RMSE) of whole fitting data achieved 74.0% and 5.0 mg^.L^-1, 80.0% and 0.3 mg^.L^-1, 90.1% and 2.9 mg^.L^-1, and 48.5% and 0^.5 mg^.L^-1 for BOD₅_VF₁, NH₄^-N_VF₁, BOD₅_VF₂, and NH₄-N_VF₂, respectively.

Multi-residue analysis of 20 mycotoxins including major metabolites and emerging mycotoxins in freshwater using UHPLC-MS/MS and application to freshwater ponds in flanders, Belgium

Mycotoxins are known for their negative impact on human and animal health as they frequently contaminate food and feed products from crop origin that are consumed by humans and animals. Furthermore, mycotoxins can leach out of plant tissue, to be transported through runoff water into nearby ponds where they can exert negative effects on aquatic organisms, such as fish, amphibians and zooplankton. The overall goal of this study was to develop a SPE-UHPLC-MS/MS method for the detection and quantification of multiple mycotoxins in amphibian breeding ponds. The method was validated and yielded acceptable within-run and between-run apparent recoveries and precision, as well as good linearity. Matrix effects (i.e. 75.7-109.6%, ≤ 17.8% RSD) were evaluated using water from 20 different ponds in Flanders, Belgium. By incorporating internal standards, overall results improved and adequate precision values (i.e. ≤ 15%) were obtained according to the EMA guideline. Additionally, extraction recovery (n = 3) was evaluated, yielding good results for all mycotoxins (i.e. 75.3-109.1%, ≤15% RSD), except for AME (i.e. 6.7 ± 0.7%), which implied the need for a matrix-matched calibration curve. Detection sensitivity was in the low nanograms per liter range. Storage stability experiments indicated that sample storage at 4 °C in amber glass bottles and analysis performed within 96 h after sampling was sufficient to avoid loss by degradation for all compounds, excluding β-ZAL and β-ZEL, for which analysis within 24 h is more indicated. The method was successfully applied to water samples originating from 18 amphibian breeding ponds situated across Flanders. Overall, enniatins B, B1 and A1 were most commonly detected at maximum concentrations of 6.9, 3.3 and 2.6 ng L^-1, respectively, followed by detection of beauvericin (1.1 ng L^-1 and < 1 ng L^-1), alternariol monomethyl ether (< 10 ng L^-1), HT2-toxin (< 40 ng L^-1), zearalenone (< 25 ng L^-1) and α-zearalanol (< 10 ng L^-1). We believe that this method will boost further research into the dynamics and ecotoxicological impact of mycotoxins in aquatic environments.

Spatial and temporal variations of greenhouse gas emissions from a waste stabilization pond: Effects of sludge distribution and accumulation

Due to regular influx of organic matter and nutrients, waste stabilization ponds (WSPs) can release considerable quantities of greenhouse gases (GHGs). To investigate the spatiotemporal variations of GHG emissions from WSPs with a focus on the effects of sludge accumulation and distribution, we conducted a bathymetry survey and two sampling campaigns in Ucubamba WSP (Cuenca, Ecuador). The results indicated that spatial variation of GHG emissions was strongly dependent on sludge distribution. Thick sludge layers in aerated ponds and facultative ponds caused substantial CO₂ and CH₄ emissions which accounted for 21.3% and 78.7% of the total emissions from the plant. Conversely, the prevalence of anoxic conditions stimulated the N₂O consumption via complete denitrification leading to a net uptake from the atmosphere, i.e. up to 1.4±0.2 mg-N m^-2 d^-1. Double emission rates of CO₂ were found in the facultative and maturation ponds during the day compared to night-time emissions, indicating the important role of algal respiration, while no diel variation of the CH₄ and N₂O emissions was found. Despite the uptake of N₂O, the total GHG emissions of the WSP was higher than constructed wetlands and conventional centralized wastewater treatment facilities. Hence, it is recommended that sludge management with proper desludging regulation should be included as an important mitigation measure to reduce the carbon footprint of pond treatment facilities.

Polluted paradise: Occurrence of pesticide residues within the urban coastal zones of Santa Cruz and Isabela (Galapagos, Ecuador)

Although pesticides are frequently used for agriculture in the Galapagos Islands (Ecuador), there are, to date, no investigations of pesticide occurrences in its coastal waters. We examined the presence of pesticide residues in the coastal waters of urban areas in two islands of the Galapagos archipelago using a repeated sampling design. Quantification was performed by solid-phase extraction, followed by chemical analysis using liquid chromatography-tandem mass spectrometry (LC-MS/MS) and gas chromatography-electron capture detector (GC-ECD). The diversity and concentration of pesticide residues in Santa Cruz island were higher compared to Isabela island. In total, sixteen pesticides were detected, including three persistent organic pollutants. Carbendazim (23.93 μg·L^-1), cadusafos (4.74 μg·L^-1), DDT (2.99 μg·L^-1), diuron (1.61 μg·L^-1) and aldrin (1.55 μg·L^-1) were detected with the highest concentrations between samples. Repetitions in locations show that concentrations of pesticide residues varied considerably in space and time. Comparison with local products indicated agricultural activities on the islands as a possible source. Furthermore, evaluation through ecological risk quotients showed that the observed concentration levels of seven pesticides pose a relatively high risk for three biotic groups (i.e. algae, invertebrates and fishes). Taken together, this study provides insights into the need to regulate, monitor and assess the presence of pesticides in the islands. At a global scale, this study is moreover valuable for the many islands that are facing the same challenges.

Highly sensitive multi-residue analysis of veterinary drugs including coccidiostats and anthelmintics in pond water using UHPLC-MS/MS: application to freshwater ponds in Flanders, Belgium

Veterinary drugs, such as coccidiostats and anthelmintics are routinely administered in extensive animal husbandry, finding their way into the aquatic environment through urine and/or feces of treated animals kept outdoors or by the application of contaminated liquid manure on agricultural fields and subsequent mechanisms of surface run-off, leaching and drift. Several of these compounds are known to exert acute and chronic toxicity effects on aquatic organisms, and can lead to changes in biodiversity and ecosystem functioning. The overall objective of this research was to develop, validate and apply a highly sensitive, multi-residue SPE-UHPLC-MS/MS method for the determination of 12 coccidiostats, registered as a feed supplement or veterinary medicine in Europe and three regularly used anthelmintics, in pond water, often functioning as amphibian habitat. Sample extraction was optimized using a fractional factorial resolution design. Pond water filtration efficiency (i.e. 80-118%, ≤25% RSD) and matrix effects (i.e. 72-119%, ≤39% RSD) were evaluated using water from respectively 3 and 20 different ponds in Flanders. By incorporating internal standards, overall results improved and adequate precision values (i.e.≤15%) were obtained according to the EMA guidelines. Acceptable within-run and between-run apparent recoveries, satisfactory precision as well as good linearity were demonstrated according to the CD 2002/657/EC, SANTE/12682/2019 and VICH 49 guidelines, except for robenidine for which the between-day precision was between 21.0 and 34.5%. Sample storage stability studies indicated that storage at 4 °C and analysis performed within 96 hours after sampling was sufficient to avoid loss by degradation for all compounds, excluding robenidine. Values for the limit of detection (LOD) and quantification (LOQ) were in nanograms per liter, which was essential for the environmental application of this novel method. The method was successfully applied on grab water samples from the water surface of 18 different ponds across Flanders, Belgium, detecting amprolium and levamisole at concentrations below the LOQ of 2.5 ng L^-1 and at 250.0 ng L^-1 or below the LOQ of 250.0 ng L^-1, respectively. In conclusion, our newly developed method may provide insights about the contamination status of amphibian breeding ponds.

Multi-class analysis of 46 antimicrobial drug residues in pond water using UHPLC-Orbitrap-HRMS and application to freshwater ponds in Flanders, Belgium

Increasing anthropogenic pressure and agricultural pollution raises concerns regarding antimicrobial resistance and biodiversity loss in aquatic environments. In order to protect and restore water resources and biodiversity, antimicrobial drug residues should be monitored in all aquatic environments including pond water. Consequently, the objective of this research was to develop and validate a novel multi-residue method for the simultaneous quantification of 46 targeted human and veterinary antimicrobial drugs in pond water. A suitable extraction method based on solid-phase extraction (SPE) was developed, assisted by a fractional factorial design. A broad polarity range of compounds was covered (log P from -4.05 to 4.38), including major representatives of the following classes: sulfonamides, tetracyclines, quinolones, macrolides, lincosamides, nitrofurans, penicillins, cephalosporins, diaminopyrimidines, pleuromutilins and phenicols. All analytes were separated using ultra-high performance liquid chromatography (UHPLC) and detected in full-scan by Orbitrap high resolution mass spectrometry (Orbitrap-HRMS). Good linearity was obtained for all compounds with R² ≥ 0.993 and goodness-of-fit coefficient (g) ≤ 11.56%. Method detection limits ranged from 10 to 50 ng L^-1 and method quantification limits were 50 ng L^-1 for all compounds. Acceptable values were obtained for within-day and between-day apparent recoveries (i.e. between 50 and 120%), precision (< 30% and < 45%) and measurement uncertainty (< 50%). Targeted analysis of 18 freshwater ponds throughout Flanders was performed to demonstrate the applicability of the newly developed UHPLC-HRMS method. Overall, 20 antimicrobial drugs were detected with highest concentrations observed for tetracyclines and their transformation products ranging between 51 and 248 ng L^-1. Finally, suspect screening was performed suggesting the presence of 14 additional pharmaceuticals including 3 antimicrobial degradation products (e.g. apo-oxytetracycline, amoxicillin penicilloic acid and penilloic acid) and 11 pesticides.

Integrated ecological modelling for evidence-based determination of water management interventions in urbanized river basins: Case study in the Cuenca River basin (Ecuador)

The growth of urbanization worldwide has contributed to the deterioration of the ecological status of water bodies. Efforts at improving the ecological status have been made either in isolated form or by means of integrated measures by stakeholders, but in many cases, these measures have not been evaluated to determine their benefit. In this study, we implemented a scenario analysis to restore the ecological water quality in the Cuenca River and its tributaries, which are located in the southern Andes of Ecuador. For this analysis, an integrated ecological model (IEM) was developed. The IEM linked an urban wastewater system (IUWS) model, which gave satisfactory results in its calibration and validation processes, with ecological models. The IUWS is a mechanistic model that incorporated the river water quality model, a wastewater treatment plant (WWTP) with activated sludge technology, and discharges from the sewage system. The ecological status of the waterways was evaluated with the Andean Biotic Index (ABI), which was predicted using generalized linear models (GLMs). The GLMs were calculated with physicochemical results from the IUWS model. Four scenarios that would enhance the current ecological water quality were analyzed. In these scenarios, the inclusion of a new WWTP with carbon, and with carbon and nitrogen removal as well as the addition of retention tanks before the discharges of combined sewer overflows (CSOs) were assessed. The new WWTP with carbon and nitrogen removal would bring about a better restoration of the ecological water quality due to better nitrogen removal. The retention tanks would help to enhance the ecological status of the rivers during rainy seasons. The integrated model implemented in this study was shown to be an essential tool to support decisions in the Cuenca River basin management.

Cost-effectiveness and screening performance of ECG handheld machine in a population screening programme: The Belgian Heart Rhythm Week screening programme

AIMS

Overall, 40% of patients with atrial fibrillation are asymptomatic. The usefulness and cost-effectiveness of atrial fibrillation screening programmes are debated. We evaluated whether an atrial fibrillation screening programme with a handheld electrocardiogram (ECG) machine in a population-wide cohort has a high screening yield and is cost-effective.

METHODS

We used a Markov-model based modelling analysis on 1000 hypothetical individuals who matched the Belgian Heart Rhythm Week screening programme. Subgroup analyses of subjects ≥65 and ≥75 years old were performed. Screening was performed with one-lead ECG handheld machine Omron® HeartScan HCG-801.

RESULTS

In both overall population and subgroups, the use of the screening procedure diagnosed a consistently higher number of diagnosed atrial fibrillation than not screening. In the base-case scenario, the screening procedure resulted in 106.6 more atrial fibrillation patient-years, resulting in three fewer strokes, 10 more life years and five more quality-adjusted life years (QALYs). The number needed-to-screen (NNS) to avoid one stroke was 361. In subjects ≥65 years old, we found 80.8 more atrial fibrillation patient-years, resulting in three fewer strokes, four more life-years and five more QALYs. The NNS to avoid one stroke was 354. Similar results were obtained in subjects ≥75 years old, with a NNS to avoid one stroke of 371. In the overall population, the incremental cost-effectiveness ratio for any gained QALY showed that the screening procedure was cost-effective in all groups.

CONCLUSIONS

In a population-wide screening cohort, the use of a handheld ECG machine to identify subjects with newly diagnosed atrial fibrillation was cost-effective in the general population, as well as in subjects ≥65 and subjects ≥75 years old.

An integrated mechanistic modeling of a facultative pond: Parameter estimation and uncertainty analysis

Imitating natural lakes, pond treatment systems inherit a high complexity with interconnected web of biochemical reactions and complex hydraulic processes. As such, its simulation requires a large and integrated model, which has been a challenge for pond engineers. In this study, we develop an all-encompassing model to gain a quantitative and comprehensive understanding of the hydraulic, physicochemical and microbiological conversion processes in the most common pond, a facultative pond. Moreover, to deal with an evitable issue of large mechanistic models as overparameterization leading to poor identifiability, a systematic parameter estimation was implemented. The application of sensitivity analysis reveals the most influential parameters on pond performance. Particularly, physical parameters, such as vertical eddy diffusivity, water temperature, and maximum growth rate of heterotrophs induce the most changes of organic matters while microbial assimilation and ammonia volatilization appear to be main processes for nutrient removal. In contrast, the efficiency of phosphate precipitation and nutrient biological removal via polyphosphate accumulating organisms and denitrifying bacteria is limited. Identifiability problems are addressed mainly by the characterization of light dependence of algal growth, interaction between water temperature and its coefficient, and the growth of autotrophic bacteria while based on the determinant measures, the most important parameter subsets affecting model outputs are related to physical processes and algal activity. After the establishment of the influential and identifiable parameter subset, an automatic calibration with the data collected from Ucubamba pond system (Ecuador) demonstrates the effect of high-altitude climatic conditions on pond behaviors. An aerobic prevailing condition is observed as a result of high light intensity causing accelerated algal activities, hence, leading to the limitation of hydrolysis, anaerobic processes, and the growth of anoxic heterotrophs for denitrification. Furthermore, the output of uncertainty analysis indicates that a large avoidable uncertainty as a result of vast complexity of the applied model can be reduced greatly via a systematic approach for parameter estimation.

Distribution of agricultural pesticides in the freshwater environment of the Guayas river basin (Ecuador)

The rapid increase and transition to more intensive agricultural activities in developing nations are often leading to misuse and overuse of pesticides, making their environment vulnerable for pesticide accumulation. In the present study, the Guayas river basin was taken as a representative case study to evaluate pesticide contamination of the Ecuadorean freshwater environment. Pesticide contamination was determined at 181 sampling sites by a multi-residue method using solid phase extraction (SPE) and pesticide residues were linked with agricultural land use activities to identify the main pollution sources. Moreover, the biological water quality status based on macroinvertebrate communities was determined at every location and the relation with the occurrence of pesticide residues was further investigated. Results showed that pesticide contamination of the freshwater environment was widely present in the Guayas river basin with detections at 108 sampling sites (60%). A total of 26 pesticide products were identified. Most frequently detected pesticides included cadusafos (62 locations), butachlor (21 locations) and pendimethalin (21 locations), with concentrations up to 0.081, 2.006 and 0.557 μg·L^-1 respectively. Pesticide residues detected in this study did not significantly influence the biological water quality (p = 0.69), but were observed to be positively correlated with ammonium concentrations, supporting the assumed combined application of chemical fertilizers and pesticides in agriculture. These pesticide residues were also associated with one or more agricultural crops, with in particular the banana and rice industries identified as major pollution sources. Both high consumption rates and non-specific application methods, such as aerial spraying of banana plantations and application directly into the water layer of irrigated rice fields, may attribute to pesticide contamination of the freshwater environment of the Guayas river basin. It is therefore suggested that measures, e.g. legal regulations and awareness campaigns, taken to prevent environmental pollution and accumulation of pesticides primarily focus on these industries.

Unimodal head-width distribution of the European eel (Anguilla anguilla L.) from the Zeeschelde does not support disruptive selection

Since the early 20th century, European eels (Anguilla anguilla L.) have been dichotomously classified into 'narrow' and 'broad' heads. These morphs are mainly considered the result of a differential food choice, with narrow heads feeding primarily on small/soft prey and broad heads on large/hard prey. Yet, such a classification implies that head-width variation follows a bimodal distribution, leading to the assumption of disruptive selection. We investigated the head morphology of 272 eels, caught over three consecutive years (2015-2017) at a single location in the Zeeschelde (Belgium). Based on our results, BIC favored a unimodal distribution, while AIC provided equal support for a unimodal and a bimodal distribution. Notably, visualization of the distributions revealed a strong overlap between the two normal distributions under the bimodal model, likely explaining the ambiguity under AIC. Consequently, it is more likely that head-width variation followed a unimodal distribution, indicating there are no disruptive selection pressures for bimodality in the Zeeschelde. As such, eels could not be divided in two distinct head-width groups. Instead, their head widths showed a continuum of narrow to broad with a normal distribution. This pattern was consistent across all maturation stages studied here.

Implications of movement for species distribution models - Rethinking environmental data tools

Movement is considered an essential process in shaping the distributions of species. Nevertheless, most species distribution models (SDMs) still focus solely on environment-species relationships to predict the occurrence of species. Furthermore, the currently used indirect estimates of movement allow to assess habitat accessibility, but do not provide an accurate description of movement. Better proxies of movement are needed to assess the dispersal potential of individual species and to gain a more practical insight in the interconnectivity of communities. Telemetry techniques are rapidly evolving and highly capable to provide explicit descriptions of movement, but their usefulness for SDMs will mainly depend on the ability of these models to deal with hitherto unconsidered ecological processes. More specifically, the integration of movement is likely to affect the environmental data requirements as the connection between environmental and biological data is crucial to provide reliable results. Mobility implies the occupancy of a continuum of space, hence an adequate representation of both geographical and environmental space is paramount to study mobile species distributions. In this context, environmental models, remote sensing techniques and animal-borne environmental sensors are discussed as potential techniques to obtain suitable environmental data. In order to provide an in-depth review of the aforementioned methods, we have chosen to use the modelling of fish distributions as a case study. The high mobility of fish and the often highly variable nature of the aquatic environment generally complicate model development, making it an adequate subject for research. Furthermore, insight into the distribution of fish is of great interest for fish stock assessments and water management worldwide, underlining its practical relevance.

Effect of the dilution rate on microbial competition: r-strategist can win over k-strategist at low substrate concentration

The conditions present in both in vitro and in vivo ecosystems determine the microbial population harbouring it. One commonly accepted theory is that a species with a high substrate affinity and low growth rate (k-strategist) will win the competition against a second species with a lower substrate affinity and higher growth rate (r-strategist) if both species are subjected to low substrate concentrations. In this study two nitrite oxidizing bacteria (NOB), Nitrospira defluvii (k-strategist) and Nitrobacter vulgaris (r-strategist), were cultivated in a continuous reactor systems. The minimal hydraulic retention time (HRT) required for maintaining the slower growing Nitrospira was first determined. A reactor containing Nitrobacter was set to the same HRT and Nitrospira was injected to evaluate the effect of the dilution rate on the competition between both species. By following the microbial population dynamics with qPCR analysis, it was shown that not only the substrate affinity drives the competition between k- and r-strategists but also the dilution rate. Experimental data and numerical simulations both revealed that the washout of Nitrobacter was significantly delayed at dilution rates close to the μmax of Nitrospira. The competition could be even reverted towards Nitrobacter (r-strategist) despite of low nitrite concentrations and dilution rates lower than the μmax of Nitrospira.

A population screening programme for atrial fibrillation: a report from the Belgian Heart Rhythm Week screening programme

AIMS

Despite the increased prevalence of atrial fibrillation (AF), data for the implementation of nationwide screening programmes are limited. The aim of this national screening study was to increase nationwide awareness about AF and stroke risk, to determine the prevalence of AF in Belgian general population using an ECG handheld machine and its feasibility to identify new AF cases.

METHODS AND RESULTS

We analysed data obtained from 5 years of the 'Belgian Heart Rhythm Week' screening programme. All subjects were screened using a one-lead ECG handheld machine. Among 65 747 subjects screened, AF was recorded in 911, with an overall prevalence of 1.4% [95% confidence interval (CI) 1.2-1.6%]. High thrombo-embolic risk, as assessed by CHA₂DS₂-VASc score ≥2, was recorded in 69% of AF subjects. In subjects with high thrombo-embolic risk, only 5.4% were treated with oral anticoagulant (OAC) and 5.8% were treated with OAC and antiplatelet drugs. Among recorded AF cases, the use of the ECG handheld machine allowed identification of 603 new AF patients (1.1%, 95% CI 0.9-1.3%). Factors associated with incident AF were chronic heart failure (P < 0.001), age (P < 0.001), diabetes mellitus (P < 0.001), previous stroke (P < 0.001), vascular disease (P < 0.001), and male sex (P < 0.001).

CONCLUSION

In this Belgian national screening programme, prevalence of AF was 1.4%. The use of an ECG handheld machine is feasible to identify a significant number of new AF cases, most with a high thrombo-embolic risk. Given the low OAC use recorded, greater efforts in AF detection and treatment are urgently needed to reduce the burden of stroke associated with this common arrhythmia.

The application of predictive modelling for determining bio-environmental factors affecting the distribution of blackflies (Diptera: Simuliidae) in the Gilgel Gibe watershed in Southwest Ethiopia

Blackflies are important macroinvertebrate groups from a public health as well as ecological point of view. Determining the biological and environmental factors favouring or inhibiting the existence of blackflies could facilitate biomonitoring of rivers as well as control of disease vectors. The combined use of different predictive modelling techniques is known to improve identification of presence/absence and abundance of taxa in a given habitat. This approach enables better identification of the suitable habitat conditions or environmental constraints of a given taxon. Simuliidae larvae are important biological indicators as they are abundant in tropical aquatic ecosystems. Some of the blackfly groups are also important disease vectors in poor tropical countries. Our investigations aim to establish a combination of models able to identify the environmental factors and macroinvertebrate organisms that are favourable or inhibiting blackfly larvae existence in aquatic ecosystems. The models developed using macroinvertebrate predictors showed better performance than those based on environmental predictors. The identified environmental and macroinvertebrate parameters can be used to determine the distribution of blackflies, which in turn can help control river blindness in endemic tropical places. Through a combination of modelling techniques, a reliable method has been developed that explains environmental and biological relationships with the target organism, and, thus, can serve as a decision support tool for ecological management strategies.

CFD study to determine the optimal configuration of aerators in a full-scale waste stabilization pond

Aerated lagoons (ALs) are important variants of the pond wastewater treatment technology that have not received much attention in the literature. The hydraulic behaviour of ALs and especially the Facultative aerated lagoons (FALs) is very complex since the aeration in these systems is designed for oxygen transfer but not necessarily to create complete mixing. In this work, the energy expenditure of the aerators was studied by means of a scenario analysis. 3D CFD models (one phase and multiphase) of a 3 ha FAL in a waste stabilization pond system in Cuenca (Ecuador) were built for different configurations of aerators. The thrust produced by the aerators was modelled by an external momentum source applied as velocity vectors into the pond fluid. The predictions of a single phase model were in satisfactory agreement with experimental results. Subsequently, a scenario analysis assessing several aeration schemes with different numbers of aerators in operation were tested with respect to velocity profiles and residence time distribution (RTD) curves. This analysis showed that the aeration scheme with all 10 aerators switched on produces a similar hydraulic behaviour compared to using only 6 or 8 aerators. The current operational schemes comprise of switching off some aerators during the peak hours of the day and operating all 10 aerators during night. This current practice could be economically replaced by continuously operating 4 or 6 aerators without significantly affecting the overall mixing. Furthermore, a continuous mixing regime minimises the sediment oxygen demand enhancing the oxygen levels in the pond.