Ecosystem-level effects of re-oligotrophication and N:P imbalances in rivers and estuaries on a global scale

Trends and ecological consequences of phosphorus (P) decline and increasing nitrogen (N) to phosphorus (N:P) ratios in rivers and estuaries are reviewed and discussed. Results suggest that re-oligotrophication is a dominant trend in rivers and estuaries of high-income countries in the last two-three decades, while in low-income countries widespread eutrophication occurs. The decline in P is well documented in hundreds of rivers of United States and the European Union, but the biotic response of rivers and estuaries besides phytoplankton decline such as trends in phytoplankton composition, changes in primary production, ecosystem shifts, cascading effects, changes in ecosystem metabolism, etc., have not been sufficiently monitored and investigated, neither the effects of N:P imbalance. N:P imbalance has significant ecological effects that need to be further investigated. There is a growing number of cases in which phytoplankton biomass have been shown to decrease due to re-oligotrophication, but the potential regime shift from phytoplankton to macrophyte dominance described in shallow lakes has been documented only in a few rivers and estuaries yet. The main reasons why regime shifts are rarely described in rivers and estuaries are, from one hand the scarcity of data on macrophyte cover trends, and from the other hand physical factors such as peak flows or high turbidity that could prevent a general spread of submerged macrophytes as observed in shallow lakes. Moreover, re-oligotrophication effects on rivers may be different compared to lakes (e.g., lower dominance of macrophytes) or estuaries (e.g., limitation of primary production by N instead of P) or may be dependent on river/estuary type. We conclude that river and estuary re-oligotrophication effects are complex, diverse and still little known, and in some cases are equivalent to those described in shallow lakes, but the regime shift is more likely to occur in mid to high-order rivers and shallow estuaries.

Individual body mass and length dataset for over 12,000 fish from Iberian streams

We provide a unique fish individual body size dataset collected from our own sampling and public sources in north-eastern Spain. The dataset includes individual body size measures (fork length and mass) of 12,288 individuals of 24 fish species within 10 families collected at 118 locations in large rivers and small streams. Fish were caught by one-pass electrofishing following European standard protocols. The fish dataset has information on the local instream conditions including climatic variables (i.e., temperature and precipitation), topography (i.e., altitude), nutrient concentration (i.e., total phosphorus and nitrates), and the IMPRESS values (a measure of cumulative human impacts in lotic ecosystems). The potential uses of this new fish dataset are manifold, including developing size-based indices to further estimate the ecological status of freshwater ecosystems, allometric models, and analysis of variation in body size structure along environmental gradients.

Pristine vs. human-altered Ebro Delta habitats display contrasting resilience to RSLR

River deltas are ecologically and economically valuable coastal ecosystems but low elevations make them extremely sensitive to relative sea level rise (RSLR), i.e. the combined effects of sea level rise and subsidence. Most deltas are subjected to extensive human exploitation, which has altered the habitat composition, connectivity and geomorphology of deltaic landscapes. In the Ebro Delta, extensive wetland reclamation for rice cultivation over the last 150 years has resulted in the loss of 65% of the natural habitats. Here, we compare the dynamics of habitat shifts under two departure conditions (a simulated pristine delta vs. the human-altered delta) using the Sea Level Affecting Marshes Model (SLAMM) under the 4.5 and 8.5 RCP (Representative Concentration Pathways) scenarios for evaluating their resilience to RSLR (i.e. resistance to inundation). Results showed lower inundation rates in the human delta (~10 to 22% by the end of the century, depending on RCP conditions), mostly due to ~4.5 times lower initial extension of coastal lagoons compared to the pristine delta. Yet, inundation rates from ~15 to 30% of the total surface represent the worst possible human scenario, assuming no flooding protection measures. Besides, accretion rates within rice fields are disregarded since this option is not available in SLAMM for developed dry land. In the human delta, rice fields were largely shifted to other wetland habitats and experienced the highest reductions, mostly because of their larger surface. In contrast, in the pristine delta most of the habitats showed significant decreases by 2100 (~2 to 32% of the surface). Coastal infrastructures (dykes or flood protection dunes) and reintroduction of riverine sediments through irrigation channels are proposed to minimize impacts of RSLR. In the worst RCP scenarios, promoting preservation of natural habitats by transforming unproductive rice fields into wetlands could be the most sustainable option.

The combined use of chemical and biochemical markers in Rutilus rutilus to assess the effect of dredging in the lower course of the Ebro River

The lower course of the Ebro River is polluted with high concentrations of organochlorine compounds dumped by a chloro-alkali plant during the last century. A remediation plan, including building of a protective wall, removal and disposal of polluted sediments started in 2012. With the aim of assessing the effects of dredging of contaminated sediments and potential alterations of water quality, areas located upstream (RR) and downstream (BE, A) the chemical plant (FL) were monitored prior (October 2012) and during dredging (June 2013) using roach (Rutilus rutilus) as sentinel organisms. Concentrations of organochlorine compounds (OCs) in fish muscle and biliary levels of polycyclic aromatic hydrocarbons (PAHs), galaxolide (HHCB) and alkyphenols (APEs) were determined together with selected enzymatic activities (7-ethoxyresorufin-O-deethylase (EROD), 7-benzyloxy-4-trifluoromethyl-coumarin O-debenzyloxylase (BFCOD) and UDP-glucuronyltransferase (UGT)) in the liver. The obtained results proved the effectiveness of the wall retaining suspended particles and avoiding further contamination of downstream sites as fish sampled at downstream sites showed up to 9-fold higher concentrations of OCs in muscle during wall construction than during dredging. EROD and UGT activities were induced in fish from downstream sites; however, no clear response to the observed pollution gradient was detected.

Effects of Camellia sinensis crude saponin on survival and biochemical markers of oxidative stress and multixenobiotic resistance of the Mediterranean mussel, Mytilus galloprovincialis

The Ebro Delta (NE Spain), formed by two bays, northern and southern hemidelta, is an area mainly devoted to rice farming. As a result shellfish species inhabiting or cultured in nearby bays may become increasingly threatened by exposure to pollutants associated to agriculture. The latest product applied in the rice fields is a plant based molluscicide called saponin, used to control populations of the giant apple snail (GAS). Saponins are known to also affect other organisms, hence research of its toxicity towards non target species is needed. In this study, mussels (Mytilus galloprovincialis) were exposed to three concentrations of suspended solid pellets containing saponin extracts from the plant species Camellia sinensis. Effects across a large set of biochemical markers were investigated in the digestive gland and gill tissue of exposed mussels to saponin for 1 and 7days. In addition, crude saponin was extracted from solid pellets to confirm its correct concentration. Results obtained from this study, showed that the concentration of crude saponin in dry pellets was of 5.5%. Lethal levels of saponin were 3.2 fold higher than the maximal predicted field concentration, and activated respiratory metabolism and expression of transmembrane protein transporters. Mussels, exposed to sub lethal concentrations of saponin showed increase of antioxidant defenses.

Effects of enhanced hydrological connectivity on Mediterranean salt marsh fish assemblages with emphasis on the endangered Spanish toothcarp (Aphanius iberus)

The hydrological connectivity between the salt marsh and the sea was partially restored in a Mediterranean wetland containing isolated ponds resulting from former salt extraction and aquaculture activities. A preliminary assessment provided evidence that ponds farther from the sea hosted very large numbers of the endangered Spanish toothcarp, Aphanius iberus, suggesting that individuals had been trapped and consequently reach unnaturally high densities. In order to achieve both habitat rehabilitation and toothcarp conservation, efforts were made to create a gradient of hydrologically connected areas, including isolated fish reservoirs, semi-isolated, and connected salt marsh-sea areas that could allow migratory movements of fish and provide some protection for A. iberus. The fish community was monitored prior to, and for three years after rehabilitation. Results showed an increase in the number of fish species within semi-isolated areas (Zone A), whereas areas adjacent to the sea (Zone B) increased the number of marine species and decreased that of estuarine species (ES). Yet overall differences in fish assemblages were much higher between zones than among study years. Generalized linear models (GLMs) evidenced that distance to the sea was the most important variable explaining the local diversity of the fish community after restoration, with occasional influence of other factors such as temperature, and depth. The abundance of A. iberus was consistently higher in semi-isolated areas at greater distances from the sea, but a decline occurred in both zones and in isolated reservoir ponds after restoration efforts, which may be attributable to interannual differences in recruitment success and, to a lesser extent, to dispersal into adjacent habitats. A negative effect of restoration works on fish population cannot be excluded, but the final outcome of the intervention likely needs a longer period.

Sea level rise impacts on rice production: The Ebro Delta as an example

Climate change and sea level rise (SLR) are global impacts threatening the sustainability of coastal territories and valuable ecosystems such as deltas. The Ebro Delta is representative of the vulnerability of coastal areas to SLR. Rice cultivation is the main economic activity in the region. Rice fields occupy most of the delta (ca. 65%) and are vulnerable to accelerated SLR and consequent increase in soil salinity, the most important physical factor affecting rice production. We developed a model to predict the impacts of SLR on soil salinity and rice production under different scenarios predicted by the Fifth Assessment Report of the Intergovernmental Panel on Climate Change by coupling data from Geographic Information Systems with Generalized Linear Models. Soil salinity data were measured in agricultural parcels and rice production from surveys among farmers. The correlation between observed and soil salinity predicted values was high and significant (Pearson's r=0.72, P<0.0001), thus supporting the predictive ability of the model. Soil salinity was directly related to distances to the river, to the delta inner border, and to the river old mouth, while clay presence, winter river flow and surface elevation were inversely related to it. Surface elevation was the most important variable in explaining soil salinity. Rice production was negatively influenced by soil salinity, thus the models predict a decrease from higher elevation zones close to the river to the shoreline. The model predicts a maximum reduction in normalized rice production index from 61.2% in 2010 to 33.8% by 2100 in the worst considered scenario (SLR=1.8m), with a decrease of profit up to 300 € per hectare. The model can be applied to other deltaic areas worldwide, and help rice farmers and stakeholders to identify the most vulnerable areas to SLR impacts.

Benthic macrofaunal dynamics and environmental stress across a salt wedge Mediterranean estuary

The spatial distribution of benthic macroinvertebrate community in relation to environmental factors was studied along the Ebro Estuary (NE Iberian Peninsula), a salt wedge Mediterranean estuary. Both ordination methods and generalized additive models were performed to identify the different benthic assemblages and their relationship to abiotic factors. Our results showed a strong relationship between macrofaunal assemblages and the predominant environmental gradients (e.g. salinity); thus revealing spatial differences in their structure and composition. Two different stretches were identified, namely the upper (UE) and the lower Ebro Estuary (LE). UE showed riverine characteristics and hence was colonized by a freshwater community; whereas LE was influenced by marine intrusion and sustained a complex marine-origin community. However, within each stretch, water and sediment characteristics played an important role in explaining species composition differences among sampling stations. Moreover, outcomes suggested a total species replacement pattern, instead of the nestedness pattern usually associated with well-mixed temperate estuaries. The sharp species turnover together with the estuarine stratification point out that the Ebro Estuary is working, in terms of ecological boundaries, under an ecotone model. Finally, despite obvious differences with well mixed estuaries (i.e. lack of tidal influence, stratification and species turnover), the Ebro Estuary shares important ecological attributes with well-mixed temperate estuaries.

Trophic interactions between native and introduced fish species in a littoral fish community

The trophic interactions between 15 native and two introduced fish species, silverside Odontesthes bonariensis and rainbow trout Oncorhynchus mykiss, collected in a major fishery area at Lake Titicaca were explored by integrating traditional ecological knowledge and stable-isotope analyses (SIA). SIA suggested the existence of six trophic groups in this fish community based on δ(13)C and δ(15)N signatures. This was supported by ecological evidence illustrating marked spatial segregation between groups, but a similar trophic level for most of the native groups. Based on Bayesian ellipse analyses, niche overlap appeared to occur between small O. bonariensis (<90 mm) and benthopelagic native species (31.6%), and between the native pelagic killifish Orestias ispi and large O. bonariensis (39%) or O. mykiss (19.7%). In addition, Bayesian mixing models suggested that O. ispi and epipelagic species are likely to be the main prey items for the two introduced fish species. This study reveals a trophic link between native and introduced fish species, and demonstrates the utility of combining both SIA and traditional ecological knowledge to understand trophic relationships between fish species with similar feeding habits.

Comparability of fish-based ecological quality assessments for geographically distinct Iberian regions

In this work we compare two Iberian and a pan-European fish-based methods to assess ecological quality in rivers: the Fish-based Index of Biotic Integrity for Portuguese Wadeable Streams (F-IBIP), the Mediterranean Index of Biotic Integrity (IBIMED) and the pan-European Fish Index (EFI+). The results presented herein were developed in the context of the 2nd phase of the Intercalibration Exercise (IC), as required by the Water Frame Directive (WFD). The IC is aimed at ensuring comparability of the quality boundaries among the different WFD assessment methods developed by the Member States for each biological quality element. Although the two national assessment methods were developed for very distinct regions of Iberia (Western and Eastern Iberian Peninsula) they share the same methodological background: both are type-specific and guild-based multimetric indices. EFI+ is a multimetric guild-based model, but it is site-specific and uses a predictive modelling approach. The three indices were computed for all sites included in the Iberian Intercalibration database to allow the direct comparison, by means of linear regressions, of the resulting three quality values per site. The quality boundary harmonization between the two Iberian methods was only possible through an indirect comparison between the two indices, using EFI+ as a common metric. The three indices were also shown to be responsive to a common set of human induced pressures. This study highlights the need to develop general assessment methods adapted to wide geographical ranges with high species turnover to help intercalibrating assessment methods tailored for geographically more restricted regions.

Monitoring the effects of floods on submerged macrophytes in a large river

The lower Ebro River (Catalonia, Spain) has recently undergone a regime shift from a phytoplankton to a macrophyte-dominated system. Macrophytes started to spread at the end of the 1990s and since 2002 artificial floods (flushing flows) of short duration (1-2 days) are released from the Riba-roja dam once or twice a year in order to reduce macrophyte density. The aim of this study was to analyse the spatiotemporal trends of the submerged macrophytes in two stretches of the lower Ebro River using high-resolution hydroacoustic methods, in order to elucidate the effects of artificial floods and natural floods on its distribution and abundance. Results showed that the mean cover in the two studied stretches (Móra and Ginestar) was not reduced after a flushing flow (from 36.59% to 55.85% in Móra, and from 21.18% to 21.05% in Ginestar), but it was greatly reduced after the natural flood (down to 9.79% in Móra and 2.04% in Ginestar); surprisingly the cover increased in Móra after the artificial flood. In order to increase the efficiency of floods in controlling macrophyte spreading, the magnitude and frequency of them should largely increase, as well as the suspended sediment load, approaching as much as possible to the original flood pattern before dam construction. Hydroacoustic methods combined with geostatistics and interpolation in GIS can accurately monitor spatiotemporal trends of submerged macrophytes in large rivers. This is the first article to apply this monitoring system to submerged macrophytes in rivers.

Regime shift from phytoplankton to macrophyte dominance in a large river: Top-down versus bottom-up effects

The lower Ebro River (Catalonia, Spain) has recently undergone a regime shift from a phytoplankton-dominated to a macrophyte-dominated system. This shift is well known in shallow lakes but apparently it has never been documented in rivers. Two initial hypotheses to explain the collapse of the phytoplankton were considered: a) the diminution of nutrients (bottom-up); b) the filtering effect due to the colonization of the zebra mussel (top-down). Data on water quality, hydrology and biological communities (phytoplankton, macrophytes and zebra mussel) was obtained both from existing data sets and new surveys. Results clearly indicate that the decrease in phosphorus is the main cause of a dramatic decrease in chlorophyll and large increase in water transparency, triggering the subsequent colonization of macrophytes in the river bed. A Generalized Linear Model analysis showed that the decrease in dissolved phosphorus had a relative importance 14 times higher than the increase in zebra mussel density to explain the variation of total chlorophyll. We suggest that the described changes in the lower Ebro River can be considered a novel ecosystem shift. This shift is triggering remarkable changes in the biological communities beyond the decrease of phytoplankton and the proliferation of macrophytes, such as massive colonization of Simulidae (black fly) and other changes in the benthic invertebrate communities that are currently investigated.

Bioaccumulation of pollutants in the zebra mussel from hazardous industrial waste and evaluation of spatial distribution using GAMs

In the Flix Reservoir (Ebro River, Spain), ca. 300,000tons of industrial waste were dumped because of the activity of a factory plant in Flix. Within the recovery program implemented, this exceptional situation provides a unique opportunity to test the value of zebra mussel as sentinel organism. Ten metal concentrations were measured in mussels from different sites to assess spatial redistribution of metals and bioavailability to the food web. Our results showed an important metal uptake by mussels; metal concentrations (except As) measured in impacted sites were up to 10 times higher than in control sites, and Mn and Hg exceeded several times the levels previously reported for polluted waters. Concentrations increased downstream showing the metal mobilization from polluted sediments in Flix Reservoir. The higher metal concentrations measured in zebra mussel individuals clearly indicated their bioavailability to the food web, allowing the toxics transfer to predators and occasionally to humans. Thus, zebra mussel is a valuable sentinel organism to identify highly polluted waters, transport routes and trophic transfer.

Life-history traits of the stone loach Barbatula barbatula

The life-history tactics of the stone loach Barbatula barbatula were studied in a Mediterranean-type climate stream (Matarranya River) located in the Ebro River basin (north-east Spain). Maximum observed ages were 2+ years in both sexes (1% of individuals), although only 0+ and 1+ year age groups were well represented. It is the lowest longevity reported for this species in its entire distribution. The seasonal growth period started in June and continued until November, but the pattern observed was different to northern populations. Barbatula barbatula in the Matarranya River was a multiple spawner, releasing small batches of oocytes between April and June. The fecundity of females was higher and the size of oocytes smaller in 1984 than in 1985. The relative fecundity (number of ripening and ripe oocytes g(-1) of fish) was lower than in northern European populations. The role of the particular environmental conditions of a Mediterranean stream was discussed in relation to the life-history tactics of B. barbatula.