Nuno Caiola

Patterns in species richness and endemism of European freshwater fish

Aim To analyse the patterns in species richness and endemism of the native European riverine fish fauna, in the light of the Messinian salinity crisis and the Last Glacial Maximum (LGM). Location European continent. Methods After gathering native fish faunistic lists of 406 hydrographical networks, we defined large biogeographical regions with homogenous fish fauna, based on a hierarchical cluster analysis. Then we analysed and compared the patterns in species richness and endemism among these regions, as well as species–area relationships. Results Among the 233 native species present in the data set, the Cyprinidae family was strongly dominant (> 50% of the total number of species). Seven biogeographical regions were defined: Western Peri-Mediterranea, Central Peri-Mediterranea, Eastern Peri-Mediterranea, Ponto-Caspian Europe, Northern Europe, Central Europe and Western Europe. The highest regional species richness was observed for Central Peri-Mediterranea and Ponto-Caspian Europe. The highest endemic richness was found in Central Peri-Mediterranea. Species–area relationships were characterized by high slope values for Peri-Mediterranean Europe and low values for Central and Western Europe. Main conclusions The results were in agreement with the ‘Lago Mare’ hypothesis explaining the specificity of Peri-Mediterranean fish fauna, as well as with the history of recolonization of Central and Western Europe from Ponto-Caspian Europe following the LGM. The results also agreed with the mechanisms of speciation and extinction influencing fish diversity in hydrographical networks. We advise the use of the seven biogeographical regions for further studies, and suggest considering Peri-Mediterranean Europe and Ponto-Caspian Europe as ‘biodiversity hotspots’ for European riverine fish.

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.

Feeding ecology of the endangered Valencia toothcarp, Valencia hispanica (Actinopterygii: Valenciidae)

Patterns of food resource utilisation by Valencia toothcarp were studied in a small coastal lagoon in the Iberian Peninsula, its most northern population. Emphasis was placed on feeding activity, diet composition, seasonal and intraspecific variation in diet and feeding strategy. Important features of the overall feeding pattern of the Valencia toothcarp are decreased feeding activity during winter, concentration on benthic invertebrates and narrow dietary breadth, with a focus on gammarid amphipods. The feeding of V. hispanica was independent of seasonal variation in food availability, because its main prey were much the same throughout the year. There were no differences in seasonal niche width and niche overlap between seasons was high. There are also no differences in diets between sexes, although diet composition, niche width and niche overlap differ between juvenile and adult fish. The Valencia toothcarp population from Santes Creus lagoon was formed by specialist individuals, with narrow niche widths, that fed on four preferential prey types, but also consumed some occasional prey.

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.

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.

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.

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.

Environmental filtering determines metacommunity structure in wetland microcrustaceans.

AbstractMetacommunity approaches are becoming popular when analyzing factors driving species distribution at the regional scale. However, until the popularization of the variation partitioning technique it was difficult to assess the main drivers of the observed patterns (spatial or environmental). Here we propose a new framework linking the emergence of different metacommunity structures (e.g., nested, Gleasonian, Clementsian) to spatial and environmental filters. This is a novel approach that provides a more profound analysis of how both drivers could lead to similar metacommunity structures. We tested this framework on 110 sites covering a strong environmental gradient (i.e., microcrustacean assemblages organized along a salinity gradient, from freshwater to brackish water wetlands). First we identified the metacommunity structure that better fitted these microcrustacean assemblages. Then, we used hierarchical variation partitioning to quantify the relative influences of environmental filters and the distance among wetlands on the identified structure. Our results showed that under strong environmental filtering metacommunity structures were non-random. We also noted that even passive dispersers, that are supposed to be poorly spatially filtered, showed spatial signals at a large geographical scale. However, some difficulties arose when inferring biotic interactions at finer-scale spatial signals. Overall, our study shows the potential of elements of metacommunity structure combined with variation partition techniques to detect environmental drivers and broadscale patterns of metacommunity structure, and that some caution is needed when interpreting finer-scale spatial signals.

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.

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 (Pearsons 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.