Owing to a long history of anthropogenic pressures, freshwater ecosystems are among the most vulnerable to biodiversity loss¹. Mitigation measures, including wastewater treatment and hydromorphological restoration, have aimed to improve environmental quality and foster the recovery of freshwater biodiversity². 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.

Local biodiversity trends over time are likely to be decoupled from global trends, as local processes may compensate or counteract global change. We analyze 161 long-term biological time series (15–91 years) collected across Europe, using a comprehensive dataset comprising ~6,200 marine, freshwater and terrestrial taxa. We test whether (i) local long-term biodiversity trends are consistent among biogeoregions, realms and taxonomic groups, and (ii) changes in biodiversity correlate with regional climate and local conditions. Our results reveal that local trends of abundance, richness and diversity differ among biogeoregions, realms and taxonomic groups, demonstrating that biodiversity changes at local scale are often complex and cannot be easily generalized. However, we find increases in richness and abundance with increasing temperature and naturalness as well as a clear spatial pattern in changes in community composition (i.e. temporal taxonomic turnover) in most biogeoregions of Northern and Eastern Europe.

The riverine dissolved organic carbon (DOC) flux is of similar magnitude to the terrestrial sink for atmospheric CO ₂ , but the factors controlling it remain poorly determined and are largely absent from Earth system models (ESMs). Here, we show, for a range of European headwater catchments, that electrolyte solubility theory explains how declining precipitation ionic strength (IS) has increased the dissolution of thermally moderated pools of soluble soil organic matter (OM), while hydrological conditions govern the proportion of this OM entering the aquatic system. Solubility will continue to rise exponentially with declining IS until pollutant ion deposition fully flattens out under clean air policies. Future DOC export will increasingly depend on rates of warming and any directional changes to the intensity and seasonality of precipitation and marine ion deposition. Our findings provide a firm foundation for incorporating the processes dominating change in this component of the global carbon cycle in ESMs.

Freshwater acidification continues to be a major problem affecting large areas of Europe, and while there is evidence for chemical recovery, similar evidence for biological recovery of freshwaters is sparse. The need for a methodology to identify waterbodies impacted acidification and to assess the extent of biological recovery is relevant to the EU Water Framework Directive, which requires methods to quantify differences in biology between impacted and unimpacted or reference sites. This study presents a new WFD-compliant metric based on diatoms (Diatom Acidification Metric: DAM) for assessing the acidification status of rivers. A database of 558 benthic diatom samples and associated water chemistry data was assembled. Diatom taxa were assigned to one of 5 indicator classes on the basis of their pH optimum, assessed using Gaussian logistic regression, and these indicator values used to calculate a DAM score for each site using weighted averaging. Reference sites were selected on the basis of their acid neutralising capacity (ANC) and calcium concentration, and a regression model developed to predict expected DAM for each site using pH and total organic carbon (TOC) concentration. Site-specific DAM scores were used to calculate ecological quality ratios ranging from ≥1, where the diatom assemblage showed no impact, to (theoretically) 0, when the diatom assemblage was indicative of major anthropogenic activities. The boundary between 'high' and 'good' status was defined as the 25th percentile of Ecological Quality Ratios (EQRs) of all reference sites. The boundary between 'good' and 'moderate' status was set at the point at which nutrient-sensitive and nutrient-tolerant taxa were present in equal relative abundance. The methodology was evaluated using long-term data from 11 sites from the UK Uplands Waters Monitoring Network and is shown to perform well in discriminating naturally acid from acidified sites.