Next Article in Journal
Climate Change Impact on Hydropower Resources in Gauged and Ungauged Lithuanian River Catchments
Next Article in Special Issue
Fish Rescue during Streamflow Intermittency May Not Be Effective for Conservation of Rio Grande Silvery Minnow
Previous Article in Journal
Uncertainty Relations in Hydrodynamics
Previous Article in Special Issue
Fish Biodiversity Conservation and Restoration, Yangtze River Basin, China, Urgently Needs ‘Scientific’ and ‘Ecological’ Action
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Water
Volume 12
Issue 11
10.3390/w12113264
water-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Factors Influencing Abundances and Population Size Structure of the Threatened and Endemic Cyprinodont Aphanius iberus in Mediterranean Brackish Ponds

by
Serena Sgarzi
1,*,
Sandra Brucet
1,2,*,
Mireia Bartrons
1,
Ignasi Arranz
1,3,
Lluís Benejam
1 and
Anna Badosa
1
1
Aquatic Ecology Group, University of Vic—Central University of Catalonia, 08500 Barcelona, Spain
2
Catalan Institution for Research and Advanced Studies (ICREA), 08010 Barcelona, Spain
3
Laboratoire Evolution et Diversité Biologique (EDB UMR 5174), Université de Tolouse, CNRS, IRD, UPS, 118 Route de Narbonne, F-31062 Tolouse, France
*
Authors to whom correspondence should be addressed.
Water 2020, 12(11), 3264; https://doi.org/10.3390/w12113264
Submission received: 29 October 2020 / Revised: 16 November 2020 / Accepted: 17 November 2020 / Published: 21 November 2020
(This article belongs to the Special Issue Ecology and Conservation of Freshwater Fishes Biodiversity)
Browse Figures
Review Reports Versions Notes

Abstract

:
Aphanius iberus is an endemic cyprinodontoid fish species of Mediterranean ponds in danger of extinction. In this study, we studied some abiotic and biotic factors that can influence A. iberus’s size structure and density in Mediterranean brackish ponds. We sampled fish using fyke nets in 10 ponds of Empordà (Spain) during the spring season. Our results showed that a better ecological status (according to the Water Quality of Lentic and Shallow Ecosystems (QAELS) index), pond’s depth and pond’s isolation (reflected by an increase in total nitrogen) were related to larger individual sizes and more size-diverse populations. Increasing the salinity is known to help the euryhaline A. iberus acting as a refuge from competitors. Nevertheless, our results showed that higher conductivities had a negative effect on A. iberus’s size structure, leading to a decrease in the mean and maximum size of the fish. Fish abundance (expressed as captures per unit of effort (CPUE)) seemed to increase with increasing the pond’s depth and total nitrogen (the latter reflecting pond isolation). In conclusion, our results suggest that achieving a better pond ecological status may be important for the conservation of endangered A. iberus, because better size-structured populations (i.e., larger mean and average lengths) were found at higher water quality conditions.

Graphical Abstract

1. Introduction

Aquatic ecosystems of the Iberian Peninsula are a hotspot for endemic freshwater fish fauna; still, most of the fish species are critically threatened by habitat destruction, intensive agricultural activities or the introduction of exotic species [1,2]. The Spanish toothcarp (Aphanius iberus, Valenciennes, 1846) is a small cyprinodont (up to 6 cm in length) endemic from the Eastern Mediterranean lowland waters of the Iberian Peninsula [3,4,5] and in danger of extinction [6,7]. As other cyprinodonts, it is characterized by fast growth, early maturity, high reproductive effort and multiple spawnings [8,9], which implies a short longevity (age up to two+). A. iberus is an eurytherm and euryhaline species, well-adapted to changes in environmental conditions [10] such as sudden alterations in temperature and salinity due to marine intrusions or freshwater floodings [11]. This cyprinodont originally inhabited a wide range of lowland waterbodies, but now, its geographical distribution is limited to brackish and hypersaline coastal waterbodies [1,5,8,12] due to habitat degradation (e.g., intensive agriculture, water pollution and wetland desiccation) and the introduction of invasive species, which usually act in an additive manner, since habitat degradation facilitates biological invasions [13,14]. The high degree of isolation among the remaining populations also poses a threat to their conservation, as they show higher rates of extinction than populations in well-connected locations [5].
Some studies have shown that the abundance and size structure of A. iberus depends on the ecological status of the ponds, with larger individuals and higher densities found in ponds with a higher water quality [15]. Indeed, in Italy, another species of the Aphanius genus (Aphanius fasciatus) has been proposed as an indicator of the ecological status of salt marshes [16], suggesting that those fish are sensitive to changes in the ecological status of their environments. More confined and less accessible ponds also seem to host populations of A. iberus more abundant and stable over time [17]. Another variable that may influence the density and size structure of A. iberus is the pond morphometry (area and depth), because it has a strong impact on the structural complexity and niche availability, as has been found for other fish species and communities [18,19,20]. However, to our knowledge, there have been not studies in this respect.
The Eastern mosquitofish (Gambusia holbrooki) [21,22] is an invasive species that very often interacts with A. iberus, because both share similar habitats [23] and compete for the same resources [24,25]. Both fish species are zooplanktivorous, but G. holbrooki consumes mainly cladocerans, ostracods and copepods [26,27], and A. iberus prefers harpacticoid copepods, copepod nauplii and detritus [24]. Sometimes, G. holbrooki can act in an aggressive way against A. iberus, and this behavior seems to be inversely proportional to the salinity, as well as its ability to capture prey [28]. Young individuals of A. iberus have been found to capture less prey in the presence of conspecific adults and G. holbrooki, suggesting both strong intraspecific and interspecific competition [12,29]. Currently, A. iberus has disappeared from fresh and oligohaline waters, and its habitat is restricted to salt marshes, coastal lagoons and river mouths [24,30], where the invasion success of the mosquitofish is limited due to the high salinity fluctuations [24,31,32]. Nevertheless, although A. iberus tolerates high salinity conditions, its metabolism may be affected when the salinity levels in the location are high [33,34,35]. Physiological functions such as oxygen consumption, critical swimming speed and routine activity level show a general decrease at the extreme salinity in Aphanius dispar [36], although the spawning efficiency seems not to be significantly affected by the changes in salinity [37].
Identifying the key factors that influence the population structure of A. iberus is relevant to develop efficient conservation and management plans for this endangered species. The body size of A. iberus has been used to assess growth-related parameters, such as age [8], fecundity and sexual maturity [38], as well as ammonia excretion rates [39]. Although, in Mediterranean brackish ponds, trophic interactions are very often body size-dependent [40,41,42], studies about the size structure of A. iberus and the factors that determine it are scarce.
The present study aims to identify the factors influencing the abundances and population size structure of A. iberus in the north-east of the Iberian Peninsula in late-spring (i.e., when this species finishes the first period of annual reproduction). Specifically, we assessed whether abiotic factors (i.e., conductivity, nutrient concentrations and pond morphology); the ponds’ ecological status; food resource availability (zooplankton biomass) and the presence of the main competitor, Gambusia holbrooki, are correlated with the size structure and abundance of this endangered species in 10 coastal brackish and hypersaline ponds. We assessed the size structure using the size diversity index [43], in addition to several size metrics, such as the maximum size, mean size and size range.
We hypothesized that a good ecological status, together with a larger pond dimension (depth and area), would increase the possibilities to find well size-structured populations with a higher size diversity, as well as higher densities of fish, as these two factors are supposed to set good conditions for the fish growth. Concerning conductivity, we expected a decrease in the size-related variables of A. iberus, because the high conductivity negatively affects its metabolism. In contrast, locations at higher conductivity levels could host higher A. iberus densities, because high conductivity may prevent the colonization of invasive species, such as G. holbrooki. We also hypothesized that the presence of G. holbrooki (main competitor of A. iberus) would lead to lower densities of A. iberus, as G. holbrooki have been observed to outcompete A. iberus [30]. Finally, we would expect that the abundance of A. iberus would be negatively correlated with the zooplankton biomass due to fish predation on zooplankton.

2. Materials and Methods

2.1. Study Area

The studied ponds are located in two protected areas of the Empordà coastal wetlands (Figure 1), between 42°01′42″ N–3°11′18″ E and 42°15′58″ N–3°08′17″ E of the Ter River Basin (NE Iberian Peninsula). Eight of the ten ponds were located in the “Aiguamolls de l’Empordà” Natural Park, and two were located south in the “El Montgrí, Illes Medes i el Baix Ter” Natural Park (some examples of individual ponds are shown in Figure S1. The climate is Mediterranean, with hot, dry summers and mild, wet winters. The hydrology of these Mediterranean coastal wetlands is characterized by a prolonged confinement period during warm seasons without water inputs, followed by irregular flooding events (i.e., rainfall or marine intrusions during sea storms, the latter no more often than twice a year) [44]. Hydrological connections among ponds and/or to rivers and the sea take place only during such irregular flooding events [45]. The studied ponds are characterized by their shallowness and cover a range of morphometry (e.g., pond area and depth) and conductivity [46]. Concentrations of inorganic nutrients (nitrates and phosphates) in late-spring are low due to the scarceness of water inputs, but concentrations of total nutrients, especially total nitrogen, are high due to a concentration effect produced by a high evaporation rate [44,45,47]. In these confined coastal environments, nitrogen, rather than phosphorous, usually limits the primary production [41,48,49]. For more details about the environmental characteristics and planktonic composition of the studied ponds, see [46].
Fish community in the studied ponds is mainly composed by Aphanius iberus and the invasive fish species Gambusia holbrooki [46]. Whereas the former is, overall, more abundant in ponds of higher salinity, the latter is more abundant in oligohaline ponds [46]. The rest of the community is composed by the Atherina boyeri (Risso, 1810), Pomatoschistus sp. (Gill, 1863), Mugil cephalus (Linnaeus, 1758), the marine fish Solea solea (Queusel, 1806), Anguila anguila (Linnaeus, 1758) and the invasive Lepomis gibbosus (Linnaeus, 1758), although the last three species are very scarce (see [46] for more details on the relative abundance of the fish species).

2.2. Field Sampling and Analysis

The 10 ponds studied were sampled once during the end of the spring season (from May to early June 2016). Conductivity (mS·cm−1) was measured using a multiparameter probe (Hanna Instruments, Woonsocket, RI, USA). Total area (m2) of each pond was estimated by using the “Google Maps Area Calculator Tool” [50] (reference 58, old paper), while the mean water column depth (cm) was calculated from in situ repeated measures obtained with a two-meter rule. Total nitrogen (mg·L−1) was measured according to Koroleff, 1973 [51], adapted by Seal Analytical to an integrated system of a CFA (Continuous Flow Analysis) digester.
Two ecological indices related to the pond’s ecological status were used in each pond: (1) the ECELS (Conservation Status of Lentic and Shallow Ecosystems) estimates the conservation status of lentic ecosystems based on morphological aspects, type of aquatic vegetation and human impacts [52]. The ECELS categories range from bad (0–30 out of 100), deficient (30–50 out of 100), mediocre (50–70 out of 100), good (70–90 out of 100) and very good (90–100 out of 100), and (2) the QAELS (Water Quality of Lentic and Shallow Ecosystems) index evaluates the water quality based on the composition of microcrustacean assemblages and taxonomic richness of aquatic insects and crustaceans in the Mediterranean wetlands [53]. The QAELS categories range from bad (<0.46), deficient (0.46–0.55), mediocre (0.55–0.62), good (0.62–0.72) and very good (≥0.72). The QAELS index was calculated after the observation of macroinvertebrate samples under optic microscope and a stereoscope. Samples were obtained through a dip net (mesh size 250 μm) following standard protocols [53].
Zooplankton samples were taken from each pond by mixing subsamples from five different sites in order to overcome the expected patchy distribution of plankton. Five liters of mixed water samples were filtered through a 50-μm mesh size net and preserved in 4% Lugol’s acid solution. Zooplankton individuals collected (including rotifers, copepods and cladocerans) were counted, identified and measured using a stereoscope and an inverted microscope (Utermöhl method), as was described in [46]. To estimate the zooplankton biomass, the total length (μm) of the first 100 individuals (when possible) was measured assuming that all individuals were equally distributed in the observed sample. Individual biomasses were then calculated using approximation to shape formulas.
Fish were caught by fyke nets set for 24 h, a common and widely used method in coastal lagoons [17,54,55]. Fyke nets consisted of a semicircular entrance ring followed by three smaller circular rings surrounded by a net (3.5-mm mesh) with two consecutive funnels (120 mm of funnel diameter, 1050 cm2 of interception area, 98 cm of length, 30 cm of height and 95 cm of wing length). The total number of fyke nets set in each pond varied according to its area and depth (Table S1). A total of 49 fyke nets were set in all the ponds. All captured fish were sexed (except juveniles < 13 mm), measured for total length (mm) and released. We measured all the individuals of Aphanius iberus in each sample, in order to minimize the error estimation.

2.3. Aphanius iberus Abundance and Size Structure

In each pond, the Aphanius iberus abundance was calculated by dividing the total captures by the number of fyke nets set in each pond (captures per unit of effort; CPUE). The size structure of A. iberus in each fyke net was assessed using four size-based metrics: (1) the maximum size, (2) mean size (computed as the geometric mean), (3) size range (the difference between the maximum and the minimum size) and (4) size diversity index (µ). For each fyke net, the size diversity was calculated using individual size (i.e., length) measurements, as proposed by [43]. Size diversity is based on the Shannon–Wiener diversity index [56] adapted for a continuous variable, such as body size. This index is the continuous analog of the taxonomic Shannon diversity index, and it produces values in a similar range to those of the Shannon index. In our case, it integrates the amplitude of the length range and relative abundance of the different lengths. Thus, the high values of the size diversity would indicate a high diversity of sizes with an equitable numerical frequency of sizes along the distributions [57,58]. In contrast, the low values of size diversity (rarely taking negative values) would indicate a low diversity of fish sizes with an inequitable numerical frequency of sizes along the distribution [43].

2.4. Data Analysis

We used mixed linear models (MLMs) to test the effects of abiotic and biotic factors on the abundance and size metrics of A. iberus. We considered captures of each fyke net as an observation unit (N = 49), and “pond” was introduced as a random effect to deal with pseudoreplication. As predictor variables, we considered conductivity, pond area and mean water depth, total nitrogen, ECELS and QAELS indexes (as estimates of ecological status) and zooplankton biomass (as food resource availability). Pearson’s r index revealed correlations among some of those variables. We applied Bonferroni correction to counteract the multiple comparisons issue and, finally, removed the variables that were highly correlated (>0.6). As G. holbrooki was not present in all the ponds (it was absent in 7 ponds out of 10), we could not include its abundance in our MLMs. Instead, we performed an ANOVA, comparing the situations “with” and “without” the presence of G. holbrooki for the abundance and for each of the size-related variables of A. iberus considered in this study.
We looked for the most parsimonious model from the full models by using a stepwise (backward) selection. The most parsimonious model was chosen using the Akaike information criteria (lowest AIC), which represents the best at explaining the data with the lowest combination of variables. We also calculated the standardized (beta) coefficients for the significant predictors included in the best models by using the R package “QuantPsyc” version 1.5 [59]. Predictors were previously checked for normality and homogeneity of variance, and, if variables did not meet the assumptions, base 10 logarithmic transformations were applied. Additionally, a visual inspection of the residual plots was done to detect any violation of the regression assumptions. In order to improve homoscedasticity, we used the function “varPower” of the package “nlme” [60]. For the creation of the boxplots, we used the “ggplot2” package [61]. All analyses were done with the software R version 3.4.2 (R core Team, 2017, Boston, MA, USA).

3. Results

3.1. Description of the Local Characteristics in the Mediterranean Ponds

There were wide ranges of environmental and ecological conditions across the study area (Table 1). Abiotic and biotic factors measured, as well as ECELS and QAELS indexes of ecological status, are shown in Table 1. The 10 studied ponds during spring showed conductivity values ranging from 10.7 mS·cm−1 to 69.10 mS·cm−1 and mean water column depths ranging from 16 cm up to 150 cm. The ponds differed quite a lot in their areas, with values ranging from 147.90 m2 to 68,150 m2. Total nitrogen, which includes organic and inorganic nitrogen compounds, showed the lowest value at 55.58 μmol·L−1 and the highest value at 234.40 μmol·L−1. The zooplankton biomass ranged from 1.13 μg·L−1 up to 4840.48 μg·L−1. Regarding the ponds’ ecological status, ECELS index values ranged from 43 to 98, indicating a “bad” and “mediocre” status, respectively. Whereas the QAELS index ranged from 0.25 to 0.56, also indicating a “bad” and a “mediocre” status, respectively (Table S1). As an average, ponds showed a “deficient” status according to the mean index value.

3.2. Variation of the Population Structure of A. iberus across Mediterranean Ponds

The mean, standard deviation and minimum and maximum values calculated for the different size metrics of A. iberus, as well as the abundance (expressed as CPUE), are shown in Table 2. Maximum length values ranged from 19 mm to 54 mm, while the mean length had a minimum of 16 mm and a maximum of 41 mm. This last value coincided with the maximum value of the length range, whereas the minimum length range value was 6 mm. Size diversity showed a wide range of values, from 0.27 to 2.34. A. iberus abundance (in CPUE) also varied largely among fyke nets, from two individuals to 525 individuals.

3.3. Main Drivers Affecting the Fish Population Size Structure and Density

The MLMs identified the most important drivers influencing the fish population size structure and density across all ponds. The most parsimonious significant models for each fish metric mentioned above (Table 2) as dependent variables are shown in Table 3. The results showed that the maximum length was negatively related to the conductivity and zooplankton biomass but positively related to the total nitrogen and QAELS index. The zooplankton biomass was the predictor with the strongest effect on the maximum length of A. iberus (Table 3).
Concerning the mean length of the fish, the MLM model showed similar results as when considering the maximum length of A. iberus as a response variable. The mean length of A. iberus significantly decreased with the increasing conductivity and zooplankton biomass but showed a positive relation with the ecological quality index QAELS.
The length range of A. iberus was found to be positively related to the pond mean depth and total nitrogen. Beta coefficients of this model (0.42 and 0.48) showed similar effects of both predictor variables on the length range. The same results were found for fish density, with the CPUE positively related to the mean depth and total nitrogen (Table 3), suggesting that a higher number of individuals inhabited more nutrient-rich and larger ponds. However, the pond mean depth exhibited a stronger correlation on the fish density than the total nitrogen (beta coefficients of 0.55 and 0.36, respectively; Table 3).
Concerning the size diversity, it was positively related to the total nitrogen and only slightly related with the conductivity (Table 3). In this case, the beta coefficients for the two main drivers were similar (0.29 for conductivity and 0.30 for total nitrogen; Table 3). Finally, the ECELS index and pond area were the only variables not retained in any of the models selected (see the full models in Supplementary Table S2).

3.4. Influence of the Presence of G. holbrooki on the Size Structure and Density of A. iberus

Boxplots showed the CPUE and size metrics of A. iberus in the presence and absence of G. holbrooki in the pond (Figure 2). Overall, the presence of G. holbrooki in the pond did not significantly modify the size structure and density of A. iberus (p-values > 0.39; Figure 2). However, the mean size of A. iberus was significantly higher when G. holbrooki was present in the pond (Figure 2), indicating an unexpected increase of body size with the presence of the main competitor.

4. Discussion

Our results suggested that both the maximum and mean sizes of A. iberus increased with the increasing pond’s water quality (QAELS index) and decreased with the increasing conductivity and zooplankton biomass. The size range, maximum size, size diversity and CPUE of A. iberus were positively related to the nutrient concentration (i.e., total nitrogen), while the size range and CPUE were also larger in deeper ponds. In contrast to our hypothesis, the presence of G. holbrooki seemed not to affect negatively the population structure of A. iberus.
We found larger maximum and mean lengths of A. iberus in locations with better water quality (i.e., the QAELS index). These results support a previous study on multiple water bodies in the southernmost distribution area of A. iberus, showing its preference for ponds of better ecological status [15]. Although the results showed that the average of the studied locations had a deficient ecological status, higher values of the QAELS index are usually associated with the predominance of large zooplankton (such as big copepods) over small rotifers (that are more linked to eutrophic and hypoxic conditions; [53,54,62]. Adults of A. iberus (with larger body lengths) are usually associated with glasswort habitats (highly productive and occasionally inundated environments) where big zooplankton is more abundant, and this may support the positive relationship between A. iberus size and the ecological status of the ponds found in our study. In contrast, younger and smaller individuals positively select more eutrophic algal mats, associated with a bad ecological status, where small rotifers dominate [24].
The total nitrogen was found to be related with the CPUE, as well as with all size metrics, expect the mean length. In Mediterranean salt marshes, the concentration of total nitrogen in the water is an indicator of a pond’s confinement or isolation [44], and during the late-spring and summer season, the total nitrogen is more concentrated because of evaporation processes [47]. Our results are in accordance with previous studies that showed that, in more confined and less accessible ponds, the populations of A. iberus are more abundant and stable over time, probably due to a lack of external perturbations, such as isolation from invasive species [17] or fewer entries of freshwater inputs. Since nitrogen is typically the limiting nutrient in these ponds, a higher total nitrogen concentration could also be associated with higher production rates of A. iberus, which, in turn, may favor larger populations.
Our results also showed, as expected, a negative relationship among conductivity and the maximum and mean sizes of A. iberus. In high saline habitats, conductivity can act as a “refuge” for the A. iberus to avoid the colonization of less salt-tolerant fish species, such as G. holbrooki [63]. However, high salinity levels may also have negative effects on the metabolism of cyprinodontoids [5,36], because the energy used for osmoregulation is not available for their growth performance and survival [33,34,35]. This could explain the presence of smaller fish in ponds at higher levels of conductivity.
In our study, the A. iberus CPUE, along with the A. iberus size range, were positively related to the pond depth. Similarly, studies from European wetlands and lakes found that wider and deeper waterbodies hosted greater biomass and sizes of fish [17,64], with the consequent higher probability to find a wider range of fish sizes [20]. Additionally, individuals of A. iberus trapped in brackish ponds due to competition exclusion and habitat degradation can, in some cases, reach unnaturally high densities [17].
Our results also showed that the zooplankton biomass is negatively correlated to the A. iberus maximum and mean sizes. Other studies on A. iberus observed that both juveniles and adults of this species have similar food preferences, as they mainly feed on harpacticoids, copepods and nauplii, detritus and diptera larvae [24]. Still, smaller individuals prefer feeding on small-sized prey, while larger fish show a greater preference for large-sized prey [24]. Larger individuals have higher feeding rates [34,65]. Thus, the presence of larger fish (expressed by higher mean and maximum sizes) may imply a lower zooplankton biomass, as it increases the consumption rates with the fish body sizes. In addition, previous studies observed that, when the potential resource availability is low, the fish size distribution tends to be more diverse, suggesting that competitive interactions for resources promote diversification by size [19,66,67].
As for the influence of G. holbrooki on A. iberus abundance and size structure, the results suggested that both the CPUE and size metrics of A. iberus were not affected by the presence of this allochthone fish. Only the mean size seemed to be slightly positively affected by the presence of the competitor. This result differed from our expectations, in which A. iberus would be smaller and less abundant in the presence of G. holbrooki. This apparent inconsistency could be explained by the fact that G. holbrooki was found just in few of the studied ponds, the ones with lower conductivity, and what we observed could be an indirect effect of environmental conditions that favor A. iberus development more than the effect of direct competition. Thus, the low number of ponds with G. holbrooki in our study did not unable us to derive strong conclusions about the influence of the G. holbrooki presence on A. iberus abundance and size structure.
In conclusion, our results suggest that the ponds’ ecological status (as shown by the QAELS index), depth, conductivity and nutrient concentrations are key variables that determine the variations of the size structure and abundance of A. iberus in Mediterranean brackish ponds. Achieving a better pond ecological status seems to be important for the conservation of endangered A. iberus, because better size-structured populations (i.e., larger mean and average lengths) are found at higher water quality conditions. In addition, a pond’s isolation may also be an advantage to preserve A. iberus populations.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4441/12/11/3264/s1: Figure S1: title, Table S1: Main geographic and morphometric characteristics of the studied ponds, along with the number of traps used and the ecological status values (QAELS and ECELS) for each pond, Table S2: Results of the MLMs (N = 49) showing the predictor variables that affect size-related variables and abundance (CPUE) of Aphanius iberus. Both Full models and Best models are presented. For each one, intercept (estimate and standard error, S.E.), Beta coefficients (standardized), t-value, significance (p-value), and df are shown.

Author Contributions

Conceptualization: S.B., A.B. and S.S.; sampling: S.S., A.B., L.B., M.B., I.A. and S.B.; statistical analysis: S.S.; writing—original draft preparation: S.S.; methodology: S.B., S.S., M.B., I.A., L.B., A.B. and writing—review and editing: S.B., S.S., M.B., I.A., L.B., A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was performed within the framework of the “SOS Fartet Project” supported by Fundación Biodiversidad and is also funded by a grant from the Spanish Ministry of Science, Innovation and Universities (grant no. RTI2018-095363-B-I00). S.B. was partially supported by a grant from the Deutsche Forschungsgemeinschaft (DFG) (grant no. Me 1686/7-1).

Acknowledgments

We thank Sergi Romero, director of “Aiguamolls de l’Empordà” Natural Park, and Marc Marí from “El Montgrí, les Illes Medes i el Baix Ter” Natural Park for the facilities and authorizations to develop the project. We are also thankful to the researchers and undergraduate students who helped during sampling: S. Carrasco, P. Ortega, M. Carol, P. Antoni, M. Omer, O. Vigil, and E. Corella. Thanks also to Z. Ersoy for her support with the chlorophyll-a and statistical analyses.

Conflicts of Interest

The authors declare no conflict of interest. The funding sponsors had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript or in the decision to publish the results.

References

  1. Doadrio, I. Atlas y Libro Rojo de los Peces Continentals de España; Ministerio del medio ambiente: Madrid, Spain, 2001. [Google Scholar]
  2. Vilà, M.; García-Berthou, E. Monitoring biological invasions in freshwater habitats. In Conservation Monitoring in Freshwater Habitats a Practical Guide and Case Studies; Springer: Berlin/Heidelberg, Germany, 2010; ISBN 9781402092770. [Google Scholar]
  3. Moreno-Amich, R. Els peixos d’aigua dolça i salabrosa. In Els Aiguamolls de l’Empordà Aspectes Ecologics, Històrics i Socials del Parc Natural; Caries, V., Sargatal, J., Felix, J., Eds.; Editorial ART-3: Figueres, Spain, 1989; pp. 241–250. [Google Scholar]
  4. De Sostoa, A. Las Comunidades de Peces del Delta del Ebro. Ph.D. Thesis, University of Barcelona, Barcelona, Spain, 1983. [Google Scholar]
  5. Oliva-Paterna, F.J.; Torralva, M.; Fernández-Delgado, C. Threatened fishes of the world: Aphanius iberus (Cuvier & Valenciennes, 1846) (Cyprinodontidae). Environ. Biol. Fishes 2006. [Google Scholar] [CrossRef]
  6. Doadrio, I.; Carmona, J.A.; Fernández-Delgado, C. Morphometric study of the Iberian Aphanius (Actinopterygii, Cyprinodontiformes), with description of a new species. Folia Zool. 2002, 51, 67–79. [Google Scholar]
  7. Crivelli, A.J. Aphanius Iberus; The IUCN Red List of Threatened Species. 2006. Available online: http://www.iucnredlist.org (accessed on 25 May 2018).
  8. García-Berthou, E.; Moreno-Amich, R. Age and growth of an Iberian cyprinodont, Aphanius iberus (Cuv. & Val.), in its most northerly population. J. Fish. Biol. 1992. [Google Scholar] [CrossRef]
  9. García-Berthou, E.; Moreno-Amich, R. Multivariate analysis of covariance in morphometric studies of the reproductive cycle. Can. J. Fish. Aquat. Sci. 1993. [Google Scholar] [CrossRef]
  10. García-Berthou, E.; Moreno-Amich, R. Ecología y conservación del fartet (Lebias ibera) en las marismas del Ampurdán (Cataluña). In Peces Ciprinodóntidos Ibéricos. Fartet y Samaruc; Generalitat Valenciana, Conselleria del medio Ambiente: Valencia, Spain, 1999; pp. 151–161. [Google Scholar]
  11. Demestre, M.; Roig, A.; De Sostoa, A.; De Sostoa, F. Contribució a l’estudi de la ictiofauna continental del delta de I’Ebre. Treb. Inst. Cat. Hist. Nat. 1977, 8, 145–226. [Google Scholar]
  12. Rincón, P.A.; Correas, A.M.; Morcillo, F.; Risueño, P.; Lobón-Cerviá, J. Interaction between the introduced eastern mosquitofish and two autochthonous Spanish toothcarps. J. Fish. Biol. 2002. [Google Scholar] [CrossRef]
  13. Mack, R.N.; Simberloff, D.; Lonsdale, W.M.; Evans, H.; Clout, M.; Bazzaz, F.A. Biotic invasions: Causes, epidemiology, global consequences, and control. Ecol. Appl. 2000. [Google Scholar] [CrossRef]
  14. Dudgeon, D.; Arthington, A.H.; Gessner, M.O.; Kawabata, Z.I.; Knowler, D.J.; Lévêque, C.; Naiman, R.J.; Prieur-Richard, A.H.; Soto, D.; Stiassny, M.L.J.; et al. Freshwater biodiversity: Importance, threats, status and conservation challenges. Biol. Rev. Camb. Philos. Soc. 2006, 81, 163–182. [Google Scholar] [CrossRef]
  15. Casas, J.J.; Sánchez-Oliver, J.S.; Sanz, A.; Furné, M.; Trenzado, C.; Juan, M.; Paracuellos, M.; Suárez, M.D.; Fuentes, F.; Gallego, I.; et al. The paradox of the conservation of an endangered fish species in a Mediterranean region under agricultural intensification. Biol. Conserv. 2011. [Google Scholar] [CrossRef]
  16. Franco, A.; Torricelli, P.; Franzoi, P. A habitat-specific fish-based approach to assess the ecological status of Mediterranean coastal lagoons. Mar. Pollut. Bull. 2009. [Google Scholar] [CrossRef]
  17. Prado, P.; Alcaraz, C.; Jornet, L.; Caiola, N.; Ibańẽz, C. Effects of enhanced hydrological connectivity on mediterranean salt marsh fish assemblages with emphasis on the endangered Spanish toothcarp(Aphanius iberus). Peer J. 2017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Holmgren, K.; Appelberg, M. Size structure of benthic freshwater fish communities in relation to environmental gradients. J. Fish. Biol. 2000. [Google Scholar] [CrossRef]
  19. Emmrich, M.; Brucet, S.; Ritterbusch, D.; Mehner, T. Size spectra of lake fish assemblages: Responses along gradients of general environmental factors and intensity of lake-use. Freshw. Biol. 2011, 56, 2316–2333. [Google Scholar] [CrossRef]
  20. Arranz, I.; Mehner, T.; Benejam, L.; Argillier, C.; Holmgren, K.; Jeppesen, E.; Lauridsen, T.L.; Volta, P.; Winfield, I.J.; Brucet, S.; et al. Density-dependent effects as key drivers of intraspecific size structure of six abundant fish species in lakes across Europe 1. Can. J. Fish. Aquat. Sci. 2016, 73, 519–534. [Google Scholar] [CrossRef] [Green Version]
  21. Lowe, S.; Browne, M.; Boudjelas, S.; De Poorter, M. 100 of the World Worst Invasive Alien Species. A Selection from the Global Invasive Species Database; The Invasive Species Specialist Group: Auckland, New Zeland, 2020. [Google Scholar]
  22. Pyke, G.H. A review of the biology of Gambusia affinis and G. holbrooki. Rev. Fish. Biol. Fish. 2005, 15, 339–365. [Google Scholar] [CrossRef]
  23. Vargas, M.J. Interaccion entre Aphanius Iberus y Gambusia Holbrooki en el Helta del Ebro: Sus Ciclos Biologicos y Ecologıas Troficas. Ph.D. Thesis, Universitat de Barcelona, Barcelona, Spain, 1993. [Google Scholar]
  24. Alcaraz, C.; García-Berthou, E. Life history variation of invasive mosquitofish (Gambusia holbrooki) along a salinity gradient. Biol. Conserv. 2007. [Google Scholar] [CrossRef]
  25. Haiahem, D.; Touati, L.; Baaziz, N.; Samraoui, F.; Alfarhan, A.H.; Samraoui, B. Impact of eastern mosquitofish, Gambusia holbrooki, on temporary ponds: Insights on how predation may structure zooplankton communities. Zool. Ecol. 2017. [Google Scholar] [CrossRef]
  26. Crivelli, A.; Boy, V. The diet of the mosquitofish Gambusia affinis (Baird & Girard) (Poecilliidae) in Mediterranean France. Rev. d’Ecologie 1987, 42, 421–435. [Google Scholar]
  27. Garcia-Berthou, E. Food of introduced mosquitofish: Ontogenetic diet shift and prey selection. J. Fish. Biol. 1999. [Google Scholar] [CrossRef]
  28. Alcaraz, C.; Bisazza, A.; García-Berthou, E. Salinity mediates the competitive interactions between invasive mosquitofish and an endangered fish. Oecologia 2008, 155, 205–213. [Google Scholar] [CrossRef]
  29. Caiola, N.; De Sostoa, A. Possible reasons for the decline of two native toothcarps in the Iberian Peninsula: Evidence of competition with the introduced Eastern mosquitofish. J. Appl. Ichthyol. 2005, 21, 358–363. [Google Scholar] [CrossRef]
  30. Alcaraz, C.; Pou-Rovira, Q.; García-Berthou, E. Use of a flooded salt marsh habitat by an endangered cyprinodontid fish (Aphanius iberus). Hydrobiologia 2008. [Google Scholar] [CrossRef]
  31. Magellan, K.; García-Berthou, E. Experimental evidence for the use of artificial refugia to mitigate the impacts of invasive Gambusia holbrooki on an endangered fish. Biol. Invasions 2016. [Google Scholar] [CrossRef]
  32. Prado, P.; Vergara, C.; Caiola, N.; Ibáñez, C. Influence of salinity regime on the food-web structure and feeding ecology of fish species from Mediterranean coastal lagoons. Estuar. Coast. Shelf Sci. 2014. [Google Scholar] [CrossRef]
  33. Brett, J.R.; Groves, T.D.D. Physiological energetics. Fish Physiol. 1979, 8, 279–352. [Google Scholar]
  34. Wootton, R.J. Ecology of Teleost Fishes; Springer: Berlin/Heidelberg, Germany, 1989. [Google Scholar]
  35. Yildirim, Ö.; Karacuha, A. Effect of salinity on growth performance and survival rate of Aphanius chantrei (Galliard, 1895). J. Appl. Ichthyol. 2008. [Google Scholar] [CrossRef]
  36. Plaut, I. Resting metabolic rate, critical swimming speed, and routine activity of the euryhaline cyprinodontid, Aphanius dispar, acclimated to a wide range of salinities. Physiol. Biochem. Zool. 2000. [Google Scholar] [CrossRef]
  37. Oltra, R.; Todolí, R. Reproduction of the endangered killifish Aphanius iberus at different salinities. Environ. Biol. Fishes 2000. [Google Scholar] [CrossRef]
  38. Vargas, M.J.; De Sostoa, A. Life-history pattern of the Iberian toothcarp Aphanius iberus (Pisces, Cyprinodontidae) from a Mediterranean estuary, the Ebro Delta (Spain). Netherlands J. Zool. 1997. [Google Scholar] [CrossRef] [Green Version]
  39. Oliva-Paterna, F.J.; García-Alonso, J.; Cardozo, V.; Torralva, M. Field studies of ammonia excretion in Aphanius iberus (Pisces; Cyprinodontidae): Body size and habitat effects. J. Appl. Ichthyol. 2007. [Google Scholar] [CrossRef]
  40. Quintana, X.D.; Boix, D.; Badosa, A.; Brucet, S.; Compte, J.; Gascon, S.; López-Flores, R.; Sala, J.; Moreno-Amich, R. Community structure in mediterranean shallow lentic ecosystems: Size-based vs. taxon-based approaches. Limnetica 2006, 25, 303–320. [Google Scholar]
  41. Badosa, A.; Boix, D.; Brucet, S.; López-Flores, R.; Gascón, S.; Quintana, X.D. Zooplankton taxonomic and size diversity in Mediterranean coastal lagoons (NE Iberian Peninsula): Influence of hydrology, nutrient composition, food resource availability and predation. Estuar. Coast. Shelf Sci. 2007, 71, 335–346. [Google Scholar] [CrossRef]
  42. Compte, J.; Gascón, S.; Quintana, X.D.; Boix, D. The effects of small fish presence on a species-poor community dominated by omnivores: Example of a size-based trophic cascade. J. Exp. Mar. Bio. Ecol. 2012, 418–419, 1–11. [Google Scholar] [CrossRef]
  43. Quintana, X.D.; Brucet, S.; Boix, D.; López-Flores, R.; Gascón, S.; Badosa, A.; Sala, J.; Moreno-Amich, R.; Egozcue, J.J. A nonparametric method for the measurement of size diversity with emphasis on data standardization. Limnol. Oceanogr. Methods 2008, 6, 75–86. [Google Scholar] [CrossRef] [Green Version]
  44. Quintana, X.D.; Moreno-Amich, R.; Comín, F.A. Nutrient and plankton dynamics in a Mediterranean salt marsh dominated by incidents of flooding. Part 1: Differential confinement of nutrients. J. Plankton Res. 1998. [Google Scholar] [CrossRef]
  45. Quintana, X.D.; Comín, F.A.; Moreno-Amich, R. Biomass-size spectra in aquatic communities in shallow fluctuating Mediterranean salt marshes (Empordà wetlands, NE Spain). J. Plankton Res. 2002. [Google Scholar] [CrossRef] [Green Version]
  46. Sgarzi, S.; Badosa, A.; Leiva-Presa, À.; Benejam, L.; López-Flores, R.; Brucet, S. Plankton taxonomic and size diversity of Mediterranean brackish ponds in spring: Influence of abiotic and biotic factors. Water 2019, 11, 106. [Google Scholar] [CrossRef] [Green Version]
  47. Badosa, A.; Boix, D.; Brucet, S.; López-Flores, R.; Quintana, X.D. Nutrients and zooplankton composition and dynamics in relation to the hydrological pattern in a confined Mediterranean salt marsh (NE Iberian Peninsula). Estuar. Coast. Shelf Sci. 2006, 66, 513–522. [Google Scholar] [CrossRef]
  48. Quintana, X.D.; Baix, D.; Badosa, A.; Brucet, S.; Gascón, S.; Gesti, J.; Trobajo, R. Limnología de los humedales costeros mediterráneos. El caso de Aiguamolls de l ’Empordà. Ecosistemas 2004, 13, 30–36. [Google Scholar]
  49. López-Flores, R.; Boix, D.; Badosa, A.; Brucet, S.; Quintana, X.D. Environmental factors affecting bacterioplankton and phytoplankton dynamics in confined Mediterranean salt marshes (NE Spain). J. Exp. Mar. Bio. Ecol. 2009, 369, 118–126. [Google Scholar] [CrossRef]
  50. Google. (n.d.). Google Maps Area Calculator Tool. 2014. Available online: https://www.daftlogic.com/projects-google-maps-area-calculator-tool.htm (accessed on 8 January 2019).
  51. Sen Gupta, R.; Koroleff, F. A quantitative study of nutrient fractions and a stoichiometric model of the Baltic. Estuar. Coast. Mar. Sci. 1973, 1, 335–360. [Google Scholar] [CrossRef]
  52. Sala, J.; Gascón, S.; Boix, D.; Gesti, J.; Quintana, X.D. Proposal of a rapid methodology to assess the conservation status of Mediterranean wetlands and its application in Catalunya (NE Iberian Peninsula). Arch. Sci. 2004, 57, 141–152. [Google Scholar]
  53. Boix, D.; Gascón, S.; Sala, J.; Martinoy, M.; Gifre, J.; Quintana, X.D. A new index of water quality assessment in Mediterranean wetlands based on crustacean and insect assemblages: The case of Catalunya (NE Iberian Peninsula). Aquat. Conserv. Mar. Freshw. Ecosyst. 2005, 15, 635–651. [Google Scholar] [CrossRef]
  54. Blanco, S.; Romo, S.; Villena, M.J.; Martínez, S. Fish communities and food web interactions in some shallow Mediterranean lakes. Hydrobiologia 2003, 506, 473–480. [Google Scholar] [CrossRef]
  55. Franco, A.; Pérez-Ruzafa, A.; Drouineau, H.; Franzoi, P.; Koutrakis, E.T.; Lepage, M.; Verdiell-Cubedo, D.; Bouchoucha, M.; López-Capel, A.; Riccato, F.; et al. Assessment of fish assemblages in coastal lagoon habitats: Effect of sampling method. Estuar. Coast. Shelf Sci. 2012. [Google Scholar] [CrossRef]
  56. Pielou, E.C. An Introduction to Mathematical Ecology; Wiley-Interscience: New York, NY, USA, 1969. [Google Scholar]
  57. Brucet, S.; Tavşanoğlu, Ü.N.; Özen, A.; Levi, E.E.; Bezirci, G.; Çakıroğlu, A.İ.; Jeppesen, E.; Svenning, J.C.; Ersoy, Z.; Beklioğlu, M. Size-based interactions across trophic levels in food webs of shallow Mediterranean lakes. Freshw. Biol. 2017, 62, 1819–1830. [Google Scholar] [CrossRef]
  58. Brucet, S.; Arranz, I.; Mehner, T.; Argillier, C.; Beklioğlu, M.; Benejam, L.; Boll, T.; Holmgren, K.; Lauridsen, T.L.; Svenning, J.C.; et al. Size diversity and species diversity relationships in fish assemblages of Western Palearctic lakes. Ecography 2018. [Google Scholar] [CrossRef]
  59. Fletcher, T.D. QuantPsyc: Quantitative Psychology Tools; R Package Version 1.5. 2012. Available online: http://cran.r-project.org/package=QuantPsyc (accessed on 15 July 2019).
  60. Pinheiro, J.C.; Bates, D.M. Extending the Basic Linear. In Mixed Effect Models in S & S-Plus; Pinheiro, J.C., Bates, D.M., Eds.; Springer Verlag: Rensselaer, NY, USA, 2000; pp. 201–270. [Google Scholar]
  61. Wickham, H. Ggplot2 Elegant Graphics for Data Analysis (Use R!); R Foundation for Statistical Computing: Vienna, Austria, 2016; ISBN 9780387981406. [Google Scholar]
  62. Brucet, S.; Boix, D.; López-Flores, R.; Badosa, A.; Moreno-Amich, R.; Quintana, X.D. Zooplankton structure and dynamics in permanent and temporary Mediterranean salt marshes: Taxon-based and size-based approaches. Arch. Hydrobiol. 2005, 162, 535–555. [Google Scholar] [CrossRef]
  63. Leonardos, I.; Sinis, A. Reproductive strategy of Aphanius fasciatus Nardo, 1827 (Pisces: Cyprinodontidae) in the Mesolongi and Etolikon lagoons (W. Greece). Fish. Res. 1998. [Google Scholar] [CrossRef]
  64. Søndergaard, M.; Jeppesen, E.; Jensen, J.P. Pond or lake: Does it make any difference? Arch. Hydrobiol. 2005. [Google Scholar] [CrossRef]
  65. Peters, R.H.; Downing, J.A. Empirical analysis of zooplankton filtering and feeding rates. Limnol. Oceanogr. 1984. [Google Scholar] [CrossRef] [Green Version]
  66. Brucet, S.; Boix, D.; López-Flores, R.; Badosa, A.; Quintana, X.D. Size and species diversity of zooplankton communities in fluctuating Mediterranean salt marshes. Estuar. Coast. Shelf Sci. 2006, 67, 424–432. [Google Scholar] [CrossRef]
  67. Quintana, X.D.; Cañedo-Argüelles, M.; Nebra, A.; Gascón, S.; Rieradevall, M.; Caiola, N.; Sala, J.; Ibàñez, C.; Sánchez-Millaruelo, N.; Boix, D. New Tools to Analyse the Ecological Status of Mediterranean Wetlands and Shallow Lakes. In Experiences from Surface Water Quality Monitoring; Springer: Cham, Switzerland, 2015; pp. 171–199. [Google Scholar]
Water 12 03264 g001 550
Figure 1. Location of the study region: Empordà coastal wetlands (Ter River Basin in the NE of the Iberian Peninsula) (C), with the detailed geographical position on the studied ponds in red (A,B). Eight of the studied ponds were located in the “Aiguamolls de l’Empordà” Natural Park (A), and two were located south in the “El Montgrí, Illes Medes i el Baix Ter” Natural Park (B). This map was produced with the online software ArcGIS (version 10.5.1, 2017, ESRI Environmental System Research Institute, Redlands, CA, USA) (https://www.arcgis.com).
Figure 1. Location of the study region: Empordà coastal wetlands (Ter River Basin in the NE of the Iberian Peninsula) (C), with the detailed geographical position on the studied ponds in red (A,B). Eight of the studied ponds were located in the “Aiguamolls de l’Empordà” Natural Park (A), and two were located south in the “El Montgrí, Illes Medes i el Baix Ter” Natural Park (B). This map was produced with the online software ArcGIS (version 10.5.1, 2017, ESRI Environmental System Research Institute, Redlands, CA, USA) (https://www.arcgis.com).
Water 12 03264 g001
Water 12 03264 g002 550
Figure 2. Boxplots showing the distribution of Aphanius iberus capture per unit of effort (CPUE) and size-related variables according to the presence (“with Gambusia”) or absence (“without Gambusia”) of Gambusia holbrooki. Significant differences are marked with the asterisk symbol (*).
Figure 2. Boxplots showing the distribution of Aphanius iberus capture per unit of effort (CPUE) and size-related variables according to the presence (“with Gambusia”) or absence (“without Gambusia”) of Gambusia holbrooki. Significant differences are marked with the asterisk symbol (*).
Water 12 03264 g002
Table
Table 1. Mean, standard deviation (SD), minimum and maximum values of the abiotic and biotic factors and ecological status indexes measured in the study ponds (N = 10). ECELS: Conservation Status of Lentic and Shallow Ecosystems and QAELS: Water Quality of Lentic and Shallow Ecosystems.
Table 1. Mean, standard deviation (SD), minimum and maximum values of the abiotic and biotic factors and ecological status indexes measured in the study ponds (N = 10). ECELS: Conservation Status of Lentic and Shallow Ecosystems and QAELS: Water Quality of Lentic and Shallow Ecosystems.
MeanSDMinimumMaximum
Mean water column depth (cm)59.3041.9316.00150.00
Pond area (m2)10,873.6020,463.50147.9068,150.00
Conductivity (mS·cm−1)46.1417.6310.0769.10
Total nitrogen (μmol·L−1)92.0151.6358.55234.40
Zooplakton biomass (µg·L−1)498.171313.921.134840.48
ECELS index77.0018.7043.0098.00
QAELSe 2010 index0.470.120.250.56
Table
Table 2. A. iberus size metrics and abundance (CPUE) obtained/computed per each sample (N = 49). The descriptive statistics are the mean, standard deviation (SD), minimum and maximum.
Table 2. A. iberus size metrics and abundance (CPUE) obtained/computed per each sample (N = 49). The descriptive statistics are the mean, standard deviation (SD), minimum and maximum.
MeanSDMinimumMaximum
Aphanius iberus maximum length (mm)41.088.0619.0054.00
Aphanius iberus mean length (mm)29.215.3616.0041.00
Aphanius iberus length range (mm)21.677.386.0041.00
Aphanius iberus size diversity1.150.360.272.34
Aphanius iberus capture per effort unit (CPUE)59.4991.682.00525.00
Table
Table 3. Results of the linear mixed models (N = 49 fyke nets) showing the predictor variables that significantly relate with Aphanius iberus size metrics and abundance (expressed as CPUE). Only the most parsimonious significant models were shown for each response variable. For each model, the intercept (estimated and standard error, S.E.), beta coefficients (standardized), t-value, significance (p-value) and degrees of freedom (df) are also reported.
Table 3. Results of the linear mixed models (N = 49 fyke nets) showing the predictor variables that significantly relate with Aphanius iberus size metrics and abundance (expressed as CPUE). Only the most parsimonious significant models were shown for each response variable. For each model, the intercept (estimated and standard error, S.E.), beta coefficients (standardized), t-value, significance (p-value) and degrees of freedom (df) are also reported.
Response VariablePredictorEstimateS.E.Beta Coefficientst-Valuep-Valuedf
Aphanius iberus MAXIMUM lengthConductivity−19.494.69−0.47−4.15<0.015
Log Total Nitrogen 2325.09564.200.454.12<0.015
QAELS index3894.26690.270.595.64<0.015
Zooplankton biomass−0.220.05−0.69−4.51<0.015
Aphanius iberus MEAN lengthConductivity−0.150.05−0.42−3.350.026
QAELS index31.946.760.554.73<0.016
Zooplankton biomass−0.01−0.01−0.59−4.010.016
Aphanius iberus length RANGELog Pond Mean Depth14.455.360.422.700.037
Log Total Nitrogen36.888.930.484.13<0.017
Aphanius iberus SIZE DIVERSITYConductivity<0.01<0.010.292.210.067
Log Total Nitrogen0.930.370.302.510.047
Aphanius iberus capture per effort unit (CPUE)Log Pond Mean Depth1.460.570.552.560.047
Log Total Nitrogen2.150.970.362.230.067
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sgarzi, S.; Brucet, S.; Bartrons, M.; Arranz, I.; Benejam, L.; Badosa, A. Factors Influencing Abundances and Population Size Structure of the Threatened and Endemic Cyprinodont Aphanius iberus in Mediterranean Brackish Ponds. Water 2020, 12, 3264. https://doi.org/10.3390/w12113264

AMA Style

Sgarzi S, Brucet S, Bartrons M, Arranz I, Benejam L, Badosa A. Factors Influencing Abundances and Population Size Structure of the Threatened and Endemic Cyprinodont Aphanius iberus in Mediterranean Brackish Ponds. Water. 2020; 12(11):3264. https://doi.org/10.3390/w12113264

Chicago/Turabian Style

Sgarzi, Serena, Sandra Brucet, Mireia Bartrons, Ignasi Arranz, Lluís Benejam, and Anna Badosa. 2020. "Factors Influencing Abundances and Population Size Structure of the Threatened and Endemic Cyprinodont Aphanius iberus in Mediterranean Brackish Ponds" Water 12, no. 11: 3264. https://doi.org/10.3390/w12113264

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop