Introduction

An important element of the biological assessment of stream conditions is an assessment of the direct or indirect effects of human activities or disturbances. Most biological assessments are based directly or indirectly on the concept of comparing the current condition to natural conditions (structure, composition, function, diversity) in the absence of human disturbance or alterations (Stoddard et al., 2006). The Reference Condition Approach (RCA; Reynoldson et al., 1997; Bailey et al., 1998; Reynoldson & Wright, 2000; Bailey et al., 2004) has been also defined under this Reference Condition (RC) concept, defined as the condition in the absence of human disturbance which is used to describe the standard, or benchmark, against which the current condition is compared. Hence, the RCA depends on the identification and sampling of reference sites and the result is a database with the biological and environmental attributes that represent this RC (Bailey et al., 2004).

This concept has been adopted by the Water Framework Directive (WFD-European Commission, 2000) since it requires the ecological status assessment, which may be expressed as a deviation from the near-natural reference condition. According to the WFD, the RC should be linked to stream typology and reference sites should present the full range of conditions expected to occur naturally within a given stream type (e.g. Barbour et al., 1996, 2004; Reynoldson & Wright, 2000; Stoddard et al., 2006). The bases for the characterisation of river typology and the identification of the RC for each stream type are provided in Annexe II of the Directive.

Before the WFD, water quality monitoring programmes in most EU Member States were mainly based on chemical and physical variables. Only around half the programmes included biological parameters in their assessment and classification of river quality (Hering et al., 2003). Nowadays, in the WFD context, the ecological status of rivers must be determined by biological quality elements (phytobenthos, macrophytes, macroinvertebrates and fishes) and has to be supported by hydro-morphological and physico-chemical quality elements (European Commission, 2000).

Many authors have defined the term RC (e.g. Davies, 1994; Hughes, 1995; Reynoldson et al., 1997; Chovarec et al., 2000; Bailey et al., 2004; Stoddard et al., 2006), and all of them emphasise that this condition corresponds to a state of very low environmental pressure or degradation. Although several methods have been proposed to establish the RC (extensive spatial survey, predictive modelling, historical data, paleo-reconstruction and expert judgment), the extensive spatial surveys are one of the most widely used approaches in rivers and streams (e.g. Wallin et al., 2003). The use of this method implies that a sufficient number of undisturbed or minimally disturbed sites are available in a specific stream type to establish the RC by statistic measures such as median values or arithmetic means. Regardless of the method followed by the EU Member States, the WFD requires the approaches being applied to provide a high level of confidence.

The recognition and selection of reference sites is a critical step in the design of extensive spatial surveys in establishing the RC (Reynoldson & Wright, 2000; Bailey et al., 2004) because these sites are the baseline for comparing test sites (e.g. Reynoldson et al., 1997; Norris & Thoms, 1999; Bailey et al., 2004). This selection of reference sites is commonly founded on a priori established criteria (Barbour et al., 1996; Reynoldson et al., 1997; Stoddard et al., 2006) based on different pressures derived from human activities that may affect the ecological status and which distinguish a reference site from a pressure-exposed site (Hering et al., 2003). Reference sites must fulfil specific operational criteria that easily indicate the absence of exposure stressors (Bailey et al., 2004). The factors chosen to be included in the a priori criteria attempt to define the least amount of environmental disturbance caused by human activities (Stoddard et al., 2006). Therefore when establishing criteria, the goal is to explicitly define the reference or acceptable healthy ecosystem (Bailey et al., 2004).

After this a priori selection, site validation must be applied to confirm and refine the selection of reference sites (Barbour et al., 1996) since some kind of disturbances are difficult to detect with the common screening methods used (Hering et al., 2003; Nijboer et al., 2004). The validation of a preliminary selection of reference sites can be specially needed in European rivers because they have been affected by multiple pressures such as organic pollution and flow regulation, among others (Hering et al., 2003). Therefore, establishing the RC by the selection method of reference sites in European rivers is expected to be problematic (Bonada et al., 2004; Nijboer et al., 2004), particularly in lowland rivers with a long history of intense flow regulation, agricultural development and pollution (Petts, 1989).

The GUADALMED Project focuses on the development of tools to facilitate the establishment of the RC in Mediterranean streams in Spain (Prat 2004). During the first phase of the project, the results showed that the use of a priori criteria was a better tool for the selection of reference sites than expert judgment because subjectivity was implicit in the experts’ interpretation of the “minimally disturbed site” concept (see Bonada et al., 2004). The main aim of the second phase was to define clear and operational criteria adapted to the particularities of Mediterranean, based on the criteria defined by Sánchez-Montoya et al. (2005), which allow the selection of validated reference sites in different Mediterranean stream types.

The specific objectives of this study were:

  1. 1.

    To review commonly used criteria to select reference sites in different regions.

  2. 2.

    To propose clear and operational criteria to select reference sites in Mediterranean streams in Spain (Mediterranean reference criteria, MRC).

  3. 3.

    To evaluate the selection of reference sites according to the MRC by comparing the values of four macroinvertebrate metrics between reference and disturbed sites for distinct stream types.

  4. 4.

    To corroborate the reference site selection by applying validation criteria.

Materials and methods

Study area and sampling sites

This study was performed in 33 basins located on a latitudinal, thermal and pluviometric gradient along the eastern coast of Spain and the Balearic Islands (Fig. 1). Basins vary in area and discharge regime, among others features. They are all influenced by the Mediterranean climate with significant spring and autumn rainfall (Köppen, 1923). Hence, all the rivers sampled show high variability in annual discharge, with frequent floods and droughts (Gasith & Resh, 1999). More information about the studied basins can be found in Sánchez-Montoya et al. (2007).

Fig. 1
figure 1

Distribution of the 162 sampling sites. Reference and non reference sites according to the final selection process

We examined 162 potential reference sites defined as “least disturbed”, following the criteria defined by Bonada et al. (2004). These sites belong to perennial or temporary streams, which are classified into five different stream types: temporary streams, evaporite calcareous streams at medium altitude, siliceous headwaters streams at high altitude, calcareous headwaters streams at medium to high altitude and large watercourses (see Sánchez-Montoya et al., 2007). Table 1 shows the main characteristics of these five stream types. The 162 studied sites are distributed along the main rivers and their tributaries (Fig. 1).

Table 1 Average and SD (in parenthesis) values of main environmental variables which characterised the five stream types established in the study area by Sánchez-Montoya et al. (2007)

Revision of reference criteria

Many criteria used to select the reference sites in various regions were compiled to not only analyse which were the most frequently used, but to determine the most appropriate criteria for Mediterranean streams. The criteria analysed were those defined by Davies (1994) in Australia, and Hughes (1995) and Barbour et al. (1996) in the United States. We also examined the criteria proposed more recently by Ehlert et al. (2002), Hering et al. (2003), Wallin et al. (2003), Bonada et al. (2004) and Nijboer et al. (2004) for European rivers for the selection of reference sites equivalent to “high status” in the context of WFD.

Selection of reference sites

After carefully reviewing the criteria used by aforementioned authors, 20 criteria were defined to be representative of Mediterranean conditions (Mediterranean Reference Criteria-MRC). More details about the selected criteria are provided in the Results and Discussion sections of this article.

To calculate the percentage of different land uses, firstly a digital terrain model (DTM 30 × 30 m; Centro Geográfico del Ejército, Ministerio de Defensa, Spain, 2005) and the Arc/Info software (version 9.0, ESRI, Redlands, California, USA, 2005) were used to delimit and calculate the water drainage area as a new polygon for all the sampling sites. Secondly, the percentages of intensive and extensive agriculture, urban zones, burnt area and natural land use were calculated by intersecting the Corine Land Cover 2000 Programme (Spain, Scale 1:100.000, Instituto Geográfico Nacional, Centro Nacional de Información Geográfica) with the water drainage area. Information on the remaining criteria was obtained from maps and available data, some of which were from the water agencies of each river basin district. This information was checked during a preliminary site inspection.

The 162 sites were classified into three “exposure to stressors” levels on the basis of the level of disturbance (1–3). These levels were determined by the number of MRC criteria fulfilled by each site. Only those that satisfied all 20 criteria were considered reference sites, which were termed undisturbed sites or very minor disturbed sites (exposure to stressors level 1). The rest which fulfilled fewer criteria were classified as disturbed sites with two degrees of degradation. If sites fulfilled from 16 to 19 criteria, they were considered minor disturbed sites (exposure to stressors level 2). Finally, sites were considered to be disturbed sites when they fulfilled fewer than 16 of the 20 criteria (exposure to stressors level 3).

To assess the aforementioned classification of sampling sites, according to the number of criteria fulfilled, values of four macroinvertebrate metrics were calculated and compared at each exposure to stressors level for the five stream types. For this purpose, the 162 sites were sampled in the spring, summer and autumn of 2003. A total of 469 macroinvertebrate samples were taken. Sites belonging to stream type 1 (temporary streams) were not sampled on each occasion because they were dry in at least one season, most frequently in summer. A multi-habitat semiquantitative kick-sample, as described by Zamora-Muñoz & Alba-Tercedor (1996), was taken on each sampling occasion using the protocol described by Jáimez-Cuéllar et al. (2004). Macroinvertebrate samples were collected from all habitats, incorporating riffles, runs and pool if these habitats were present in the stream stretch in question, with a kick-net (250–400 μm), and they were preserved in 100% ethanol. Effort was allocated in proportion to the occurrence of each habitat. Samples were examined under a stereoscope in the laboratory and at least 200 individuals in each sample were randomly picked and identified at the family level, except for Ostracoda, Oligochaeta and Hydracarina. Large uncommon individuals were picked individually, as described in Barbour et al. (1999).

The Iberian Biomonitoring Working Party Biotic Index (IBMWP; Alba-Tercedor & Sánchez-Ortega 1988; Alba-Tercedor et al., 2004) and the IASPT (Iberian Average Score per Taxon) were applied as biotic integrity indices and the total number of families (NFAM) and the number of EPT taxa (number of families belonging to Ephemeroptera, Plecoptera and Trichoptera) were applied as single metrics. The mean values of macroinvertebrate metrics and indices of the three sampling occasions were used to perform multiple box-plots to represent the range of variation of each metric within the three exposures to stressors levels for the five stream types considered.

Validation of reference sites and final selection

In order to refine and confirm the a priori selection of reference sites by applying of the MRC, a validation process was applied. The validation method must be preferably an approach which does not depend on the method used to select the reference sites (Nijboer et al., 2004). Macroinvertebrate community data were selected as indicators of biological conditions to corroborate whether a priori reference sites presented undisturbed or near-natural conditions. Biological threshold defined by several authors as a validation criterion of undisturbed sites in Mediterranean streams to confirm that near natural conditions was used. This threshold was the value of the IBMWP > 100 established by Alba-Tercedor & Sánchez-Ortega (1988) as a high quality class indicator. This validation criterion was also applied to Mediterranean streams in Portugal by Chaves et al. (2006). This index was designed mainly to detect organic pollution in permanent streams and therefore this threshold can not be applied for temporary streams where, additionally, large differences in macroinvertebrate assemblages were observed as a result of variability in salinity and hydrology (Sánchez-Montoya et al. 2007). A priori reference sites had to fulfil the biological validation criterion in the three sampling occasions (spring, summer and autumn) to be considered as a reference site in the final selection.

Results

Review of the criteria used to select reference sites

To facilitate further analysis, all the criteria were grouped into three general groups of attributes: (a) biological (Table 2), (b) physical and chemical (Table 3) and (c) hydro-morphological (Table 4). Furthermore, these groups were arranged into 12 categories following some authors′ classifications. These categories were: riparian vegetation and floodplain, introduced species, fisheries and aquaculture, biomanipulation, wildlife, catchment vegetation, macroinvertebrate community, point source pollution, diffuse source pollution and land use, channel morphology, habitat composition, hydrologic conditions and regulation.

Table 2 Criteria related to biological attributes applied to select reference sites in previous studies
Table 3 Criteria related to physical–chemical features applied for the selection of reference sites in previous studies
Table 4 Criteria related to hydromorphological features applied for the selection of reference sites in previous studies

Proposal of criteria for Mediterranean streams in Spain (Mediterranean reference criteria, MRC)

On the basis of the reference criteria used by other authors and considering the particular conditions of Mediterranean streams, we proposed 20 criteria that a site has to fulfil to be considered a reference site (Table 5). The spatial scale of application was sub-basin (upstream of the site) for all criteria. These 20 criteria include a wide range of human uses and disturbances on streams in our study area, as well as some general aspects on naturalness which must be present if a site to be considered a reference site.

Table 5 Proposed criteria for Mediterranean stream assessment (Mediterranean Reference Criteria, MRC) according to the WFD

A priori selection and validation of reference sites selected by means of the MRC

Most of the sites included in the GUADALMED project (92 out of 162; 57%) fulfilled all 20 criteria and were, therefore, considered reference sites in the a priori selection. Of the total sites, 36% were classified as minor disturbed and only 7% as disturbed sites (Table 6). On the basis of the number of reference sites, the least disturbed stream type was siliceous headwaters (84% of all stream type 3 sites were undisturbed), followed by temporary streams (71% of all stream type 1 sites were undisturbed), calcareous median (62% of all stream type 4 sites were undisturbed) and evaporate calcareous (31% of all stream type 2 sites were undisturbed). Therefore, reference sites were found for all stream types, except for stream type 5 (large watercourses).

Table 6 Number of sites according to the number of Mediterranean Reference Criteria fulfilled for each stream types (T1: stream type 1, T2: stream type 2, T3: stream type 3, T4: stream type 4 and T5: stream type 5)

In general, the most limiting criteria were those related to natural land use and to dry land farming, which affected more than 20% of the sites. Additionally, disturbances caused by the presence of reservoirs, retention of sediments, exotic species and irrigated farming affected a considerable percentage of sites (10–14%). Canalisation, adequate substrates and grazing were the criteria which all the sites fulfilled.

In particular, the most limiting criteria for type 1 (temporary streams) were those mostly related to land use and the presence of exotic species (Fig. 2). Land use and hydrologic conditions were the most restraining criteria with the type 2 sites (evaporite calcareous streams at medium altitude). For stream type 3 sites (siliceous headwater streams at high altitude), river morphology and hydrological conditions were the most limiting criteria, while types 4 and 5 were limited by riparian vegetation and land use, and by river morphology and hydrological conditions, respectively.

Fig. 2
figure 2

Percentage of sites which failed criteria for the six categories from the total of non reference sites for each stream type. (RV, Riparian Vegetation zone; IS, Introduced Species; PP, Point sources of Pollution; DP, Diffuse sources of Pollution and land uses; RM, River Morphology and habit conditions and HC, Hydrological Conditions and regulation)

Reference sites presented higher median values of the four metrics studied than those sites that failed some of the criteria (Fig. 3). In general, as the level of exposure to stressors increased, the median values of macroinvertebrate metrics decreased, with some exceptions. The median and 25th percentile values for number of families, EPT and IBMWP for stream type 2 were lower for exposure to stressors level 2 than for exposure to stressors level 3. Type 3 reference sites (siliceous headwaters) presented the highest median values, followed by types 4, 2 and 1. These values for types 2 and 4 were the most similar, and this pattern was also observed in exposure to stressors level 2, but not for the most disturbed sites, probably because of the scarcity of sites at exposure to stressors level 3.

Fig. 3
figure 3

Multiple box-plots for NFAM, EPT, IBMWP and IASPT values for all five stream types for the levels of exposure to stressors considered. Boxes are interquartile ranges (25th percentile to 75th percentile) ranges bars show maximum and minimum values and lines represent medians. (T1: stream type 1, T2: stream type 2, T3: stream type 3, T4: stream type 4 and T5: stream type 5)

The IBMWP > 100 threshold was used as a validation criterion. All the a priori reference sites exceeded this value on the three sampling occasions. Therefore the 92 a priori reference sites were selected as true reference sites in the final selection. Since all the sites were included in the final selection, consistency was observed between the selection and the validation method, which implies that sites were appropriately selected.

Discussion

The application of criteria for the recognition and selection of reference sites has been widely reported (e.g. Davies, 1994; Hughes, 1995; Barbour et al., 1996; Reynoldson & Wright, 2000; Ehlert, 2002; Hering et al., 2003; Wallin et al., 2003; Bonada et al., 2004; Nijboer et al., 2004). A review of the criteria used by the above-mentioned authors showed that, in relation to the number of criteria examined, the most relevant categories are: riparian vegetation, diffuse and point sources of pollution, channel morphology and hydrological conditions and regulations. However other aspects such as introduced species, biomanipulation, presence of wildlife, catchment vegetation, and macroinvertebrate and habitat composition were not considered by many authors.

Mediterranean Reference Criteria

In this study, six categories were used to establish criteria to select the reference sites in Mediterranean streams (see Table 5). Stoddard et al. (2006) explained that these criteria, which are developed for the purpose of establishing the least amount of human ambient disturbance, vary across ecological regions as the characteristics of the landscape and the human use of the landscape vary. In this sense, we believe that the conventional criteria developed in more temperate areas may not always be applicable or sufficient to identify a reference condition in Mediterranean regions. Hence, the particular features of Mediterranean streams and the most frequent pressures have been taken account to define and apply these criteria and, for that reason, some consideration about these criteria must be emphasised.

Riparian vegetation cover in Mediterranean streams is highly variable, depends on the system type (permanent, ephemeral or intermittent) and reflects the prevailing hydrologic regimen (Salinas et al., 2000). In fact, natural riparian vegetation of some typical ramblas, intermittent or ephemeral streams (Williams, 2006) in Mediterranean regions may be scarce or even absent as a result of water stress and floods (Suárez et al., 2004; Gómez et al., 2005). Therefore, these differences in stream types must be taken account when applying the criterion related to the presence of natural cover and the composition of riparian vegetation.

Regarding to the criterion related to introduced species, the fact that water basins in Spain have been affected for decades by the introduction of exotic species must be considered. These additions produce a biotic homogenisation of freshwater fauna, like the increased similarity of biota over time caused by the replacement of native species with nonindigenous species (Rahel, 2002). This uniformity is of special concern for fish fauna as 38% of the freshwater species found in Spain are exotic. Native fish fauna is currently threatened by predation and competition with these exotic species (Elvira, 1998; Elvira & Almodóvar, 2001; Vila-Gispert et al., 2005; Clavero & García-Berthou, 2006). Moreover in some basins examined in this study, both the red swamp crayfish (Procambarus clarkii) and the red-eared slider (Trachemys scripta) were present. In this context the presence of some exotic fauna produces highly disturbed stream communities regardless of upstream land use activity.

The absence of urban and industrial effluent discharges upstream of potential reference sites must be considered. These effluent discharges upstream of potential reference sites have a greater effect on ephemeral streams than on other types of streams. At very a low flow or no flow, there can be little effluent dilution, which results in high concentrations of contaminants, particularly nutrients.

Mediterranean streams, like other watercourses in semi-arid regions, are particularly susceptible to water diversion (directly or via groundwater withdrawal) and flow regulation (Gasith & Resh, 1999). The flow fluctuations below reservoirs are highly variable (Petts, 1984; Casado et al., 1989; Bernez et al., 2002) and depend on the inherent properties of the stream, as well as oscillations in the demand of water. As a result, it is difficult to define a minimal downstream distance of recovery (McAllister et al., 2001) of hydrological and geomorphological conditions for the macroinvertebrate community and other biological components. As specific studies are lacking the affected stream reaches in the study area that help to define that minimal distance, we considered that a reference site cannot be located downstream of a dam. In addition, a criterion related to inter-basin water transfers, which are common in Mediterranean areas and lead to serious alterations in both the donor and the receptor basins (Castella et al., 1995; Prat & Ibañez, 1995; Davies et al., 2000; Gibbins et al., 2000), was included, seeing as that sites affected by the flow regulation from this activity cannot be considered as presented in reference conditions.

Specific but very frequent pressures that occur in Mediterranean areas have been considered in this work by the inclusion of novel criteria which, to our knowledge, have not been previously considered. One such case is the criterion related to the effects of wildfire. Effects of fire on stream ecosystems are attributed mainly to physical disturbances associated with increased runoff and sediment scouring, both of which produce long-term changes in abiotic (i.e. abrupt alterations in water chemistry) and biotic properties of lotic ecosystems, such as alterations in the availability of food resources and stream retention capacity (e.g. Minshall et al., 1989; Mishall et al., 1997; Gresswell, 1999). The number of fires occurring in eastern Spain has dramatically increased since the mid 70s, and it is attributed to land use and climate changes (Pausas, 2004). These observations support the inclusion of this criterion when selecting reference sites, because a reference sites cannot be affected by a significant burnt area. The threshold established for the percentage of burnt area has been determined according to studies conducted in Mediterranean areas (Inbar et al., 1998; Mayor et al., 2007).

The presence of sand and gravel extraction activities is also a novel criterion. These activities, together with the presence of dams, are the main causes of sediment deficiency and the altered granulometric balance noted in many rivers of Spain (Batalla, 2003). Therefore, we believe that the sites affected by these activities cannot present reference conditions.

Drainage channels from traditional agricultural floodwater systems are a typical and frequent pressure in Mediterranean streams (Hooke & Mant, 2002). These channels may increase nutrient and pesticide concentrations in streams because they return the water use to irrigate agricultural land to streams. Hence, a criterion to avoid the presence of these drainage channels upstream of a site was included.

Finally, we believe that the quantitative state of groundwater must be included as an indicator of ecosystem health. Groundwater use is especially intense and important in Mediterranean areas, where human uses of water often compromise the natural function of stream ecosystems (Gasith & Resh, 1999). Therefore we consider that near natural groundwater level is an essential criterion to take account of in Mediterranean streams.

A priori selection of reference sites

Of the 162 sites studied as potential reference sites, 92 were identified as a priori reference sites by applying of the Mediterranean reference criteria, which showed that approximately half of these potential reference sites had no pressures and, consequently, presented natural conditions. Although many sites were influenced by multiple stressors, reference sites were found for all stream types except for large watercourses at medium and low altitude (stream type 5). This stream type included sites allocated at a lower altitude that the rest of stream types in the main axes of rivers and consequently the presence of multiple stressors, that are more common in lowlands, made the selection of undisturbed sites impossible. This problem has been described by other authors in European streams (e.g. Ehlert et al., 2002; Hering et al., 2003; Nijboer et al., 2004; Lorenz et al., 2004; Dodkins et al., 2005), particularly in Mediterranean streams (e.g. Prat & Munné, 2000; Bonada et al., 2004). The analysis of the pressures identified in each stream type in this study shows that one of the most limiting criteria to find reference sites in all the stream type was the presence of a significant percentage of dry land use, except for the siliceous headwater streams which include sites allocated in the highest regions of the study area (see Table 1). This finding was expected because traditionally the agriculture has been a very important productive sector in the Mediterranean Europe (Allen, 2001), where dry and irrigated farm lands are significantly higher than in other European countries (MIMAN, 2000). Another important pressure presented in all the stream types was the alteration of the natural hydrological conditions caused by the presence of reservoirs, except in temporary streams due the low mean annual discharges that characterise this type of streams. This is in the agreement with the knowledge that Mediterranean basins in European countries have a high level of flow regulation given the presence of reservoirs whose storage capacity represents approximately two-thirds of the total capacity of European countries (MIMAM, 2000).

In situations where undisturbed sites are not available to generate a database of reference sites, as is the case of large watercourse at medium and low altitude stream in this study, the use of historical data has been suggested as an additional method to describe reference communities (Ehlert et al., 2002; Nijboer et al., 2004). However, this approach is impossible to apply in this study area because in Spain, as well as in other Mediterranean countries, data of biotic communities covering large catchment areas are recent (e.g. Pinto et al., 2004; Chaves et al., 2006). In a similar context, it has been suggested that those sites which have received less intense impacts fulfilling most of the reference criteria, or meet them closely, might be justifiably used to establish the RC (e.g. Davies, 1994; Hering et al., 2003; Bailey et al., 2004). In this sense, we believe that minor disturbed sites (sites at level 2 which fulfilled more than 15 of the 20 criteria) belonging to large watercourses (stream type 5), where it is not possible to find undisturbed sites, could be considered suitable to establish RC in the sense of the Beast Attainable Condition (BAC) as defined by Stoddard et al. (2006). This term indicates that sites are located in areas where the impact of inevitable land use on biota is minimised. In conclusion, this permissiveness regarding the presence of some pressures in reference sites allows environmental and management goals in order to apply the WFD since other aforementioned approaches are not possible in this specific stream type.

We hypothesised that the values of the macroinvertebrate metrics and indices studied for each type-specific reference site would be higher than those from disturbed sites. This hypothesis was confirmed because undisturbed sites for all the stream types presented higher medians than disturbed sites at different exposure of stressors levels, proving that the analysis of macroinvertebrate communities independently confirm the a priori site selection. This result demonstrates that the criteria proposed for Mediterranean streams in this study reflect a true gradient of disturbance where the values of metrics and indices decreased as site fulfilled lower number of criteria.

Validation of reference sites

Identifying the RC in areas that are strongly affected by human activities is difficult because of the occurrence of natural and human generated stress and the interaction between stressor-type and the capacity of individual ecosystems (or types of systems) to resist or to recover from disturbances. In this situation, site validation is the next key step to confirm and refine the previous selection of reference sites (Barbour et al., 1996); regardless of the method used to establish the RC, this process must be applied preferably using independent approaches (Nijboer et al., 2004). In this study the macroinvertebrate community was used for validation purposes by applying the reliable biological threshold defined by Alba-Tercedor & Sánchez-Ortega (1988) and using high quality class indicators in Mediterranean streams (e.g. Chaves et al., 2006).

Although the structure of the macroinvertebrate community should be not used for the a priori selection of reference sites (e.g. Barbour et al., 1996; Bailey et al., 2004), to avoid preconceived ideas about the composition or diversity of the community in typical sites, this community has been recommended to corroborate the RC and, therefore, the high ecological status in the WFD context (e.g. Wallin et al., 2003; Bailey et al., 2004; Nijboer et al., 2004; Chaves et al., 2006). The use of macroinvertebrate communities as a validation criterion is justified since the knowledge that they reflect the influence of all the stressors on its environment (Karr, 1999; Karr & Chu, 2000), and are capable of detecting the presence of certain types of disturbances that are difficult to recognise with the common screening methods used. In this study the validation process let to a 100% confirmation of the a priori reference sites, showing that the applying criteria proposed acts as an appropriate tool to select references sites in Mediterranean streams.

The validation of temporary stream by applying the biological threshold was not possible given the knowledge that macroinvertebrate community in reference sites included in this stream type differs greatly from the rest of permanent stream types (Sánchez-Montoya et al., 2007). It must be also emphasised that temporary streams in the study area includes diverse macroinvertebrate communities, and that a further division into subtypes has been suggested because of the large variability in salinity and hydrology conditions (Sánchez-Montoya et al., 2007). Thus, the application of only one value of the IBMWP threshold to identify the RC in temporary stream is not likely to be suitable to ensure a proper validation of this stream type.

The aforementioned large differences in the reference macroinvertebrate communities between temporary streams and the rest of permanent stream types have been confirmed in this study since the values of the studied metrics in reference sites for this stream type are lower than in the rest of seasonal permanent stream types (see Fig. 3). The temporary stream type in the study area includes natural intermittent or ephemeral sites which underwent intense drought periods, and which lead to a predominance of pool habitats frequently disconnected from riffles. The degree of habitat connectivity may influence assemblage richness and diversity. Low taxa richness should be expected in disconnected habitat (Bonada et al., 2006) because low connectivity produces low richness given the exchanges of matter, energy and organism that are constrained between parches (Ward et al., 1999). Based on the knowledge that OCH families taxa are predominant in pools, because disconnected pools could be a refuge from desiccation for some tolerant lentic macroinvertebrates (i.e. OCH) but probably not for rheophilic macroinvertebrates (i.e. EPT) (e.g. Brown & Brussock, 1991), lower values of both taxa richness (family level) and EPT richness (relative to OCH richness) in temporary streams is expected (Bonada et al., 2006; Sánchez-Montoya et al., 2007). Besides, the lower IBMWP values in temporary streams in the RC is also accounted by the score of EPT families in the IBMWP (mean value = 8.4) which is significantly higher than in OCH families (mean = 4.9). Additionally, spatial heterogeneity increases from ephemeral to intermittent and to permanent sites because higher flows result in an increase of available refuges and substrate diversity (Bonada et al., 2007) and higher richness taxa is expected with increases of this heterogeneity. In conclusion, lower metric values in the RC in temporary streams, if compared to permanent streams, must be attributed to natural differences in macroinvertebrate communities.

The 25th percentile values of the IBMWP for the distribution of reference sites selected in this study for each stream type are higher than those obtained previously by Alba-Tercedor et al. (2004) in similar stream types. This finding corroborates that the methodology to recognise and select reference sites proposed in this study is an appropriate method to identify undisturbed sites in Mediterranean streams.