Response of biotic communities to salinity changes in a Mediterranean hypersaline stream

Background This study investigates the relationship between salinity and biotic communities (primary producers and macroinvertebrates) in Rambla Salada, a Mediterranean hypersaline stream in SE Spain. Since the 1980's, the mean salinity of the stream has fallen from about 100 g L-1 to 35.5 g L-1, due to intensive irrigated agriculture in the watershed. Furthermore, large dilutions occur occasionally when the water irrigation channel suffers cracks. Results Along the salinity gradient studied (3.5 – 76.4 g L-1) Cladophora glomerata and Ruppia maritima biomass decreased with increasing salinity, while the biomass of epipelic algae increased. Diptera and Coleoptera species dominated the community both in disturbed as in re-established conditions. Most macroinvertebrates species found in Rambla Salada stream are euryhaline species with a broad range of salinity tolerance. Eight of them were recorded in natural hypersaline conditions (~100 g L-1) prior to important change in land use of the watershed: Ephydra flavipes, Stratyomis longicornis, Nebrioporus ceresyi, N. baeticus, Berosus hispanicus, Enochrus falcarius, Ochthebius cuprescens and Sigara selecta. However, other species recorded in the past, such as Ochthebius glaber, O. notabilis and Enochrus politus, were restricted to a hypersaline source or absent from Rambla Salada. The dilution of salinity to 3.5 – 6.8 gL-1 allowed the colonization of species with low salininty tolerance, such as Melanopsis praemorsa, Anax sp., Simulidae, Ceratopogonidae and Tanypodinae. The abundance of Ephydra flavipes and Ochthebius corrugatus showed a positive significant response to salinity, while Anax sp., Simulidae, S. selecta, N. ceresyi, N. baeticus, and B. hispanicus showed significant negative correlations. The number of total macroinvertebrate taxa, Diptera and Coleoptera species, number of families, Margalef's index and Shannon's diversity index decreased with increasing salinity. However, the rest of community parameters, such as the abundance of individuals, evenness and Simpson's index, showed no significant response to changes in salinity. Classification and ordination analysis revealed major differences in macroinvertebrate community structure between hypersaline conditions (76.4 g L-1) and the rest of the communities observed at the lower salinity levels, and revealed that below ~75 g L-1, dissimilarities in the communities were greater between the two habitats studied (runs and pools) than between salinity levels. Conclusion Salinity was the first factor determining community composition and structure in Rambla Salada stream followed by the type of habitat.


Background
Saline inland waters (salinity ≥ 3 g L -1 ) are geographically widespread, especially in arid and semi-arid regions of the world. Such habitats are globally threatened by human activities that induce changes in the natural hydrology and salinity levels, with a corresponding loss of biodiversity [1]. Hydrological alteration and associated changes in salinity in aquatic systems are frequently associated to changes in land use in the corresponding watershed [2,3]. Natural salt lakes have been the focus of most limnological studies of saline waters. Environmental threats to salt lakes have been well documented [1], but those to saline streams, especially hypersaline streams, have been poorly studied. This is a pity because hypersaline running waters are particularly interesting due to their halotolerant/halophilic biota and the rarity and the high number of endemic species that they sustain: for example, some species of aquatic Coleoptera [4][5][6], which often have restricted geographical ranges and occur as highly isolated populations [7].
Although, the effects of increasing salinity (secondary or anthropogenic salinisation) on freshwater ecosystems and saline lakes have been extensively reviewed, especially in Australia, producing a reduction in diversity [8][9][10][11][12][13][14][15], the contrary process, the reduction of salinity in naturally saline systems is less well documented [16]. In the Mediterranean area, the reduction of salinity in naturally saline systems has been an increasingly common phenomenon in recent decades due principally to changes in agricultural practices, such as the expansion of irrigated agriculture in the watersheds. Díaz et al. [16] emphasized the fragility of hypersaline lagoon ecosystems in the face to anthropogenic disturbances, such as increases in freshwater inflow and nutrient inputs.
The purpose of our study was to investigate the relationship between salinity and the biomass of primary producers, and the composition and structure of macroinvertebrate communities in a hypersaline stream, together with its influence on biodiversity at local and regional scales. Our hypothesis is that salinity is the principal factor affecting the structure of the community, and that an extended period of salinity reduction may increase the abundance and richness of species, but reduce the richness or abundance of the most halophilic species, resulting in a possible reduction in biodiversity at regional level.

Study area
Rambla Salada is a Mediterranean hypersaline stream in the sedimentary Fortuna basin which belongs to the Segura River watershed (southeast of the Iberian Peninsula, Figure 1). The climate in the basin is semiarid with a mean annual temperature of 18°C and annual precipitation of around 250 mm. A long warm and dry summer season is interrupted by spring and autumn rains, the latter followed by a temperate winter. The watershed, of 4470 ha, is characterised by one main channel 11.6 km long and temporary tributaries that flow only during rainy periods. The main channel has permanent flow in its middle and lower sections, but flow is intermittent at the head. The natural cover on the watershed is open Mediterranean scrub, although much is dedicated to citrus and horticultural crops (intensive agriculture).
The high salinity in Rambla Salada is principally due to Miocene gypsiferous marls in the basin [17]. The ionic composition of the streamwater is dominated by Na + and Clions, followed by SO 4 2-and Ca + [18].
Since 1979 the Tajo-Segura river diversion has brought water for irrigation, and agriculture in the watershed has changed from the cultivation of extensive dry crops to intensively irrigated crops. There has been an increase in the input of freshwater, nutrients, pesticides and other pollutants derived from the crops into the stream. As consequence, since the 1980s, the salinity of the stream has fallen from close to 100 g L -1 [19] to a mean of 35.5 g L -1 . The diversion channel that crosses the stream suffers cracks at least once at year, which produce substantial discharges of freshwater into the stream. One such accident occurred in 2003 with the result that on 2 October the Rambla Salada stream received a discharge of 400 L.s -1 , reaching its lowest recorded salinity level (3.5 g L -1 ) as the channel was emptied for repair ( Figure 2). Figure 1 Geographic location of the Rambla Salada stream and land uses in the basin.

Geographic location of the Rambla Salada stream and land uses in the basin
The reference macroinvertebrate community dates from 1983, in the most hypersaline conditions before the important land use changes took place, but is only based on dispersed data of coleoptera and heteroptera species presence [20][21][22][23][24]. In addition, there is no control site upstream, because of the different conditions (lower salinity, temporary flow and human impact). Thus, taking into account the rapid community recovery in saline streams after flash flood disturbance [25], we have considered the re-established conditions of Rambla Salada stream one year after the crack disturbance as the state most similar to natural hypersaline conditions to investigate the effect of salinity changes on the composition and structure of biological communities.

Sampling and processing of the samples
To remove the seasonal variation component, sampling was performed on four dates: June and October 2003 during the dilution disturbance, and the same months in 2004 when the system had been reestablished. The four dates pointed to an in-site salinity gradient (3.5 -76.4 g L -1 ). The study was made in the middle section of the stream, near an abandoned salt-pan exploitation (lat. 38°7'28"N, long. 1°6'46"W) (Figures 3 and 4). Pools and runs constitute the principal aquatic habitats in the reach ( Figure 5). Shallow pools (≤ 60 cm) occupy more than 80% of the channel surface. They are characterised by large deposit of silts covered by a biofilm with diatoms and cyanobacteria as primary producers. Runs are characterised by gravel and sand as dominant substrates, are <10 cm deep and have a maximum water velocity of 0.37 m.s -1 . This habitat was covered principally by the filamentous algae Cladophora glomerata.
On each date, water temperature and dissolved oxygen were measured in situ at 15 minute intervals during 24 hours using a multi-parametric recorder (WTW, MultiLine P4). Conductivity was measured in the morning with an ECmeter (TetraConR 325) that automatically calculates salinity. Discharge was estimated from measurements of depth and current velocity along a cross-section in the run. Cover percentage of the different primary producers was visually estimated in each habitat.
In the run, three water samples were taken to measure suspended solids, chlorophyll a and dissolved nutrients (ammonium, nitrate, nitrite, and phosphate concentrations). The samples were kept cool until their arrival in the laboratory, where they were filtered onto preashed and preweighed GF/F glass-fiber filters. The filters were ovendried at 60°C to constant dry weight and then ashed at 450°C for 4 h to estimate the input of suspended sediments (SS) and particulate organic matter (POM). Chlorophyll a concentration was determined by spectrophotometry, following extraction in 90% acetone. Ammonium was converted to ammonia by adding 10 M NaOH solution and measured with an ammonia electrode connected to a pH/mV meter. The rest of the dissolved inorganic nutrients were determined according to standard methods [26]: nitrate by the cadmium reduction method, nitrite by sulfanylic acid colorimetry, and phosphate by ascorbic acid colorimetry.
In both habitats, benthic macroinvertebrates were collected from the bottom using a 0.014 m 2 core. Substrate was removed by hand and a sieve (250 μm mesh size) was passed across the bottom and through the water column until no individuals were found. In lentic habitat, nektonic macroinvertebrates were also collected from the water column using a big core (0.1 m 2 ). Five replicate cores per habitat were taken randomly on each sampling date, and material was preserved and stored in 75% ethanol until separation and identification in the laboratory.
Five benthic organic matter samples of sediments were collected by cores of 5.3 cm diameter and 24.5 cm length. Sediment samples were dried at 60°C to constant dry weight and then ashed at 450°C for 4 h to estimate the content of ash free dry weight (AFDW).
The biomass of the different primary producers in each habitat was estimated from five replicate samples using for macrophytes species the same core type as for sediments. For epipelic algae, we used a minicore of 2 cm diameter and 2 cm length. Samples were stored on ice and frozen until processing. Sixty milligrams of fresh macrophyte were taken from each sample to determine Chl a concentration by spectrophotometry, following extraction in acetone. The rest of the macrophyte was dried at 60°C to constant dry weight (DW). Epipelic samples were filtered onto Whatman A filters to determine Chl a concentration and AFDW.

Data analyses
Macroinvertebrate data from each habitat were analysed both separately and combined to create a composite sample, representing one date.
The following structural macroinvertebrate community parameters were calculated: density, richness (number of total taxa, Diptera and Coleoptera species, and families). In addition, several diversity indices (Margalef's index, evenness, Shannon's and Simpson's diversity indices) were calculated using PRIMER (V.9). Comparisons were made using nonparametric Mann-Whitney U paired tests to compare differences between disturbed and re-established conditions in a given month. Spearman rank correlation was used to analyse the relationships of each community parameter with the salinity. Statistical analyses were performed using STATISTICA (V. 6). An analysis of similarity (UPGMA, Bray-Curtis distance metric) was used to compare differences in macroinvertebrate communities between habitats and under the different salinity conditions. Finally, we used canonical correspondence General aspect of the study area in Rambla Salada stream Figure 4 General aspect of the study area in Rambla Salada stream. October 2004.
Aerial view of the Rambla Salada stream, showing the location of the study area, the abandoned salt-pans, and the Tajo-Segura river diversion channel Figure 3 Aerial view of the Rambla Salada stream, showing the location of the study area, the abandoned salt-pans, and the Tajo-Segura river diversion channel.
analysis (CCA) with a down weight of the rare species to extract patterns of variation in macroinvertebrate communities in relation to environmental variables. Multivariate analyses were performed using MVSP (V.3.1) after log (x+1) data transformation.

Physical and chemical environment
Water temperatures were slighter higher during 2004 than 2003 for the same months. The highest discharge value was reached in October 2003 due to large freshwater inputs from the Tajo-Segura channel, while the minimum discharge was reached in October 2004 after a long drought period. Salinity ranged from 3.5 to 76.4 g L -1 , coinciding with the maximum and minimum discharge values (Table 1). Water was well oxygenated in the morning on all dates due to high benthic primary production, although night values of dissolved oxygen were low due to high ecosystem respiration (unpublished data). Mean daily percentage of oxygen saturation exceed 100%.
The availability of total dissolved inorganic nitrogen (DIN) was higher than phosphate availability throughout the study period. The dominant form of nitrogen varied with salinity. High amounts of reduced forms (NH 4 + and NO 2 -) were registered in October 2004, when the highest salinity value was reached; while the oxidized form (NO 3 -) predominated at all other times. Phosphate concentrations were higher on the disturbed dates than on the reestablished dates. The nutrient levels measured in July 1983 [19], considered as the natural hypersaline reference conditions, were lower than those observed during the studied period, although very close to that registered in June 2004 (Table 1) pool and run habitats ( Figure 6a). There were significant differences between both June samplings but not between October samplings (Table 2). No significant correlation was found between FPOM and salinity (Table 3).

Primary producers
During the study period, primary producers in the stream were represented by an extended biofilm containing diatoms and cyanobacteria that covered the sediments, the filamentous algae Cladophora glomerata and Enteromorpha intestinalis, and the widgeon-grass Ruppia maritima.Enteromorpha intestinalis covered less than 5% of the channel and was not sampled for the estimation of biomass. The three variables studied (coverage, biomass (gAFDW.m -2 and gChl a m -2 ) showed high correlation coefficients for each primary producer.
In general, the biomass of macrophytes was higher in the pool than in the run habitat (Figure 6b-c). The biomass of R. maritima was significantly different between dates (Table 2), Chl a concentration being higher in 2003 (disturbed conditions) than in 2004 (reestablished conditions). Cladophora glomerata showed significant differences in biomass between the October months, being much higher in the more dilute conditions. However, the biomass of epipelic algae was significantly higher in October 2004 than in October 2003, but no significant differences were found between the June months ( Figure  6d and Table 2). Cladophora glomerata and R. maritima showed a significant decrease in Clh a content with increasing salinity, while the biomass of epipelic algae increased (Table 3 and Figure 6b-d).

Macroinvertebrates
The macroinvertebrate community was characterized by low species diversity. A total of 25 taxa representing 16 families were recorded in the Rambla Salada stream during the study period (Table 4). Diptera and Coleoptera were the richest insect orders, with 10 and 9 species, respectively. Total taxon richness ranged from 11 to 15 in the run, and from 9 to 17 in the pool.

Individual taxon response
Eight of the species collected were known to be essential components of the natural hypersaline reference community of Rambla Salada from 1983 and were abundant during the studied period, both in disturbed and re-  Nebrioporus baeticus and S. selecta showed significant differences in abundance between the October months, being higher in the more dilute conditions (Table 4) Table 5). The rest of species showed no significant abundance response to the salinity changes.

Community response
The number of taxa, Margalef's index and Shannon's index differed significantly between the extreme conditions of the salinity gradient ( Figure 7 and Table 2). Values of these indices were significantly higher in disturbed than in reestablished conditions. Significant changes in the number of Diptera and total species were also observed between both summer months. Macroinvertebrate communities showed a significant increase in richness, Margalef's index and Shannon's diversity index with decreasing salinity (Table 6). However, the rest of the community parameters, such as the abundance of individuals, evenness and Simpson's index, showed no significant response to the salinity changes.
Variation in run and pool habitats of the content of fine particulate organic matter of sediments (a) and biomass of Ruppia mar-itima (b), Cladophora glomerata (c) and epipelic algae (d) in relation to the salinity gradient The analysis of similarity, which compared communities by habitat, showed the greatest differences in macroinvertebrate assemblages between October 2004 and the rest of dates (41% of dissimilarity distance, Figure 8). In October 2004, the distance between both habitats was 39.5 %. On the disturbed dates, the differences between habitats were greater than between salinity levels (3.5 -6.8 g L -1 ). Samples from June 2004, showed an intermediate position, reflecting their position in the salinity gradient, with 27.4% of distance between habitats.
Canonical axes from the CCA analysis pointed to a high percentage of variability in the samples, with 28 % on the first axis and 27 % on the second. The first axis was related to the salinity gradient and differentiated samples from October 2004, when the highest salt water concentration was reached, from the rest of samples (Figure 9a). Besides salinity, this axis was positively related to ammonium, nitrite and suspended sediments, and negatively associated to dissolved oxygen, water Chl a concentration and R. maritima biomass. The second axis differentiated between samples from the run (positive scores) and pool (negative scores) habitats. The second axis was associated positively with discharge and biomass of R. maritima.
Species with a wide range of salinity tolerance and present in both habitats (pools and runs) appeared near the centre of the variable scores plot (Figure 9b

Discussion
Salinity changes ecosystem functioning by altering abiotic and biotic processes [13]. Changes in salinity affect aquatic biota directly or indirectly. Salinity can influence, for example, pH, dissolved oxygen and the nutrient balance of plants. It may be that biota respond as much to these indirect effects as to the direct ionic effects of salt concentration [27]. In our study, salinity had a strong effect on the biomass of primary producers in Rambla Salada stream, the biomass of C. glomerata and R. maritima increasing as salinity levels fell, while the biomass of epipelic algae decreased. Dilution changed the ecosystem's state from epipelic algae dominated to filamentous algae dominated. This result was expected because macrophytes are often absent or present in insignificant amounts at high salinity values [28]. In saline lakes, Wolheim and   Lovvorn [29] found that aerial cover and both the species and structural diversity of macrophytes often decline as salinity increases. Although both the macrophyte species studied are halotolerant and are frequently found in mesosaline waters, they were adversely affected when salinity exceeded 75 g L -1 . Ruppia maritima has specialized features, such as epidermal leaf cells modified to absorb both cations and anions for osmoregulation [30,31], which enable it to survive in a variety of saline conditions and at high temperature beyond those tolerated by other submerged angiosperms [32]. Widgeon-grass on the other hand, turns white and dies when exposed to high temperatures and widely varying salinity levels, because of the additional energy required for increased osmoregulation [31]. The same plant adapts poorly to turbid waters or anaerobic sediments [32], conditions prevalent in Rambla Salada stream in October 2004, when its biomass decreased significantly. Changes in macrophytic production with decreasing salinity might result from a combination of lower osmoregulation costs and greater nutrient availability. The high discharge and large input of freshwater would have mixed the sediments and associated nutrients within the water column and increased phosphorous availability. In the disturbed diluted conditions, high nitrate and phosphate concentrations favoured C. glomerata and R. maritima growth. Mason [33] found that the high productivity of C. glomerata was directly associ-ated with high nitrate and high biochemical oxygen demand in farm ponds. However, ammonia and nitrite dominated in the hypersaline conditions of Rambla Salada stream and could have been the principal N source for epipelic algae. Although a diversity study of the phototrophic microbial community was not undertaken, in hypersaline conditions, high species diversity is expected, as Major et al. [34] found in Salt Plains of Oklahoma (USA).
On the other hand, because the form of primary production can change with salinity, involving shifts in the relative abundance of open-water and macrophyte/epiphyte habitats, indirect changes in macroinvertebrate communities are also to be expected [35]. These new habitats can support additional invertebrate community components [36]. Cladophora glomerata provide attachment sites for epiphytes, which can serve as food source for invertebrate grazers [37], besides supplying a refuge from predation [38].
The direct effects of salinity changes on macroinvertebrate community depend to a great extend on the salinity tolerance range of the species concerned [39]. Due to the absence of studies of osmotic regulation and salinity tolerance of the macroinvertebrate species occurring in Rambla Salada stream, and because there is no salt sensitivity database of Paleartic species as there is for Australian biota [27], salinity tolerance must be inferred from the maximum salinity at which species were collected. For Coleoptera and Heteroptera species, the upper salinity levels were obtained from a database of field collections carried out since 1980's in the Segura River basin, which comprises aquatic systems across a wide salinity gradient ( Figure 10). For the rest of the taxa, we obtained information of the upper salinity levels of the family, genera or species (if possible) from published literature.
Dipterans and coleopterans showed the greatest diversity in the saline stream studied like as in saline lakes [40]. Among the dipterans showing the highest salinity tolerance, Ephydridae and Dolichopodidae families are known to be present in saline lakes up to 118 g L -1 of salt concentration [41]. Chironomids have been collected at salinities greater than seawater [40]: for example, C. salinarius and H. varians that are found frequently in brackish waters in coastal lagoons and salt marshes [42,43]. Among Coleoptera, Hydraenids, Dytiscids and Hydrophilids contain species that lives in broad ranges of salinity.
The genus Ochthebius has a great number of halophilic species, O. glaber exhibiting the highest salinity tolerance, being present throughout the entire salinity range and able to inhabit saturated brine ( Figure 10). This species is followed in decreasing order of salinity tolerance by O. notabilis, O. corrugatus, O. cuprescens, O. delgadoi and O. tudmirensis. Among ditiscids, N. ceresyi and N. baeticus larvae have been recorded at salinities up to 100 and the adults to up to 150 g L -1 . Among the hydrophilids, E. falcarius and B. hispanicus are the two more abundant species in mesosaline waters (~20 -~50 g L -1 ), although have been found living at 81 g L -1 . Enochrus bicolor, a halobiont species that can be found sporadically in Rambla Salada stream, is confined to saltmarsh/brackish water in the UK in a salinity range of 4.7-62.6 g L -1 [44]. The corixid S. selecta, typical species of both brackish and inland saline waters [45], was found in salinities ranging from 3.5 to 134 g L -1 , completing its life cycle up to 55 g L -1 .
Mayflies are generally restricted to freshwaters, but some Baetidae species are common in moderately saline lakes [40]. Cloeon schoenemundi has been found both in fresh and hyposaline streams of semiarid areas of the Segura River watershed [46] and in Rambla Salada stream it occurred up to 75 g L -1 , the maximum registered for any mayfly.
The New Zealand mud snail, P. antipodarum, is highly tolerant to salinity and resistant to desiccation [47]. These characteristics, together with its asexual reproduction and Variation of different composition and structural macroinvertebrate community parameters in relation to the salinity gradient Figure 7 Variation of different composition and structural macroinvertebrate community parameters in relation to the salinity gradient. The conceptual model developed by Herbst [48] to explain the abundance of Ephydra hians along a salinity gradient is useful for understanding species response to salinity dilution. This model proposes that abundance is maximal at salinity levels midway between the physiological limitations of high-salinity stress and the ecological limitations imposed by biotic interactions (e.g., predation and competition) at low salinity. On the other hand, the ability of many halophiles and halobionts to tolerate lower salinities in the laboratory than those to which they are exposed in the field suggests that biotic factors exclude them from habitats of lower salinity [40].
Several of the diversity indices changed along the salinity gradient. Total richness, the number of Diptera and Coleoptera species, the number of families, Margalef's index and Shannon's diversity index increased with decreasing salinity. However, the rest of community parameters, such as the abundance of individuals, evenness and Simpson's index, showed no significant response to changes in salinity. Higher densities and lower richness, diversity and evenness were found by Bunn and Davies [11] for benthic biota during river salinization when salinity exceeded 3 g L -1 .
Our results were to be expected, since the macroinvertebrate diversity follows a trend parallel to that shown by other animals and plants: i) richness and diversity indices fall as conditions become more extreme [49], and ii) at lower salinities, most species are halotolerant while at higher salinities most of them are halophilic [50].
Although assemblages changed gradually along the salinity gradient, at a given level of salinity there was an abrupt and major change in the composition and structure of the macroinvertebrate community. The analysis of similarity confirms that significant changes toock place in the community structure once salinity exceeded 75 g L -1 . Bellow this salinity level, difference in flow velocity, the organic content of sediments and the biomass of macrophytes between pool and run habitats had a greater effect on macroinvertebrate communities than salinity changes. Similarly, some authors have found that small changes in salinity are not well related to macroinvertebrate [51,52] and microbial [34] communities, because of the broad ranges of salinity and temperature tolerance of inhabitants of saline systems in arid and semi-arid environments. Streams in Mediterranean climate regions are influenced by a sequence of regular and often extreme flooding and drying periods and their inhabitants are adapted to extreme, often rapid, environmental change [53]. However, prolonged dilution or long-term effects of decreasing salinity by anthropogenic disturbances might result in substantial biological changes, and the most halophilic biota will be adversely affected as salinity falls below 75 g L -1 . For example, O. glaber, an endemic halophilic species of high genetic diversity and great conservation interest [54,55] is known to have disappeared in recent decades from the study area, where it remains only in a small hypersaline source (~150 g L -1 ) in the Rambla Salada stream. Ochthebius notabilis, another hydraenid species proposed as vulnerable species in the red list of Spanish invertebrates [55], frequent in salt-pans and hypersaline pools habitats [56], has not been found in Rambla Salada stream since 1991 (unpublished data). Our results indicate that although diversity at local level increases with dilution, with the introduction of estenohaline species with low to moderate salinity tolerance and wide geo-Bray-Curtis dendrogram based on macroinvertebrate assem-blages Figure 8 Bray-Curtis dendrogram based on macroinvertebrate assemblages. Samples are labelled with three digits: month, year, and habitat (R: run, P: pool). Canonical correspondence analysis plots representing samples and environmental variables (a), and macroinvertebrate species and environmental variables (b) on first and second axis Figure 9 Canonical correspondence analysis plots representing samples and environmental variables (a), and macroinvertebrate species and environmental variables (b) on first and second axis. Samples are labelled with three digits: month, year, and habitat (R run, P pool). Species are labelled with codes (see Table 4). Codes of environmental variables: BCL, biomass of Cladophora; BRU, biomass of Ruppia; BEP, biomass of epipelon; SSE, suspended sediments; OMS, organic matter of sediments; DIS, discharge; OXI, dissolved oxygen; and CLA, concentration of chlorophyll a in water. a b graphical ranges, there is a risk that the aquatic diversity at regional and global levels will be reduced as the most halophilic and vulnerable species in the Mediterranean area, which often have restricted geographical ranges and occur as highly isolated populations [7], are eliminated.

Conclusion
This study contributes to our understanding of the response of hypersaline stream biota to dilution and adds information concerning the saline thresholds of primary producers and macroinvertebrates. Our results support the initial hypothesis that dilution causes an increase in richness and biotic diversity, but a reduction in the abundance or elimination of the most halophilic species. There is a critical level of ~75 g L -1 when an abrupt and major change in the structure of the macroinvertebrate community occurs, but below which occasional changes in salinity do not greatly affect community structure. Nevertheless, more studies of saline streams are needed to refine and confirm this assertion. Changes in macroinvertebrate communities due to other factors, such as habitat type and biotic interactions will reduce the relevance of changes due to salinity. Chronic dilution may produce a loss of biodiversity at regional level.
Salinity tolerance for Coleoptera and Heteroptera species recorded at past and recent dates in the Rambla Salada stream, obtained from distributional data for the aquatic sys-tems of the Segura River basin Figure 10 Salinity tolerance for Coleoptera and Heteroptera species recorded at past and recent dates in the Rambla Salada stream, obtained from distributional data for the aquatic systems of the Segura River basin.