Science of the Total Environment 409 (2011) 612–624

Contents lists available at ScienceDirect

Science of the Total Environment

j o u r n a l h o m e p a g e : w w w. e l s ev i e r. c o m / l o c a t e / s c i t o t e n v

Assessing land-use effects on water quality, in-stream habitat, riparian ecosystems

and biodiversity in Patagonian northwest streams
María Laura Miserendino a,⁎, Ricardo Casaux a, Miguel Archangelsky a, Cecilia Yanina Di Prinzio a,
Cecilia Brand a, Adriana Mabel Kutschker b
a
CONICET Laboratorio de Investigaciones en Ecología y Sistemática Animal-Universidad Nacional de la Patagonia San Juan Bosco, Argentina
b
Cátedra de Ecología-Facultad de Ciencias Naturales - Universidad Nacional de la Patagonia San Juan Bosco, Argentina

a r t i c l e i n f o a b s t r a c t

Article history: Changes in land-use practices have affected the integrity and quality of water resources worldwide. In
Received 15 April 2010 Patagonia there is a strong concern about the ecological status of surface waters because these changes are
Received in revised form 14 October 2010 rapidly occurring in the region. To test the hypothesis that greater intensity of land-use will have negative
Accepted 16 October 2010
effects on water quality, stream habitat and biodiversity we assessed benthic macroinvertebrates, riparian/
Available online 20 November 2010
littoral invertebrates, fish and birds from the riparian corridor and environmental variables of 15 rivers
(Patagonia) subjected to a gradient of land-use practices (non-managed native forest, managed native forest,
Keywords:
Water quality
pine plantations, pasture, urbanization). A total of 158 macroinvertebrate taxa, 105 riparian/littoral
Macroinvertebrates invertebrate taxa, 5 fish species, 34 bird species, and 15 aquatic plant species, were recorded considering
Fish all sites. Urban land-use produced the most significant changes in streams including physical features,
Bird conductivity, nutrients, habitat condition, riparian quality and invertebrate metrics. Pasture and managed
Pastures native forest sites appeared in an intermediate situation. The highest values of fish and bird abundance and
Forest management diversity were observed at disturbed sites; this might be explained by the opportunistic behavior displayed by
these communities which let them take advantage of increased trophic resources in these environments. As
expected, non-managed native forest sites showed the highest integrity of ecological conditions and also great
biodiversity of benthic communities. Macroinvertebrate metrics that reflected good water quality were
positively related to forest land cover and negatively related to urban and pasture land cover. However, by
offering stream edge areas, pasture sites still supported rich communities of riparian/littoral invertebrates,
increasing overall biodiversity. Macroinvertebrates were good indicators of land-use impact and water quality
conditions and resulted useful tools to early alert of disturbances in streams. Fish and birds having a greater
ability of dispersion and capacity to move quickly from disturbances would reflect changes at a higher scale.
© 2010 Elsevier B.V. All rights reserved.

1. Introduction et al., 2007). In this sense, current threats to global biodiversity require
understanding of how progressively or locally human activities might
All over the world the ecological integrity of river systems is being scale up to affect regional biotas (Manel et al., 2000). This phenomenon
endangered by land-use changes. This widespread and ever-increasing has captured the attention of scientists dedicated to the study of
phenomenon involving urbanization, agriculture, pasture conversion, communities living in aquatic ecosystems and those associated with
deforestation, and the replacement of native species by exotic ones with their riparian corridors (Danger and Robson, 2004; Diamond et al.,
commercial value, represents a real and pervasive threat to the 2002; Meador and Goldstein, 2003; Saunders et al., 2002; Thompson
biodiversity and conservation of lotic ecosystems (Allan 2004; Morley and Townsend, 2004).
and Karr, 2001; Paul and Meyer, 2001; Scrimgeour and Kendall, 2003). Widespread removal of basin vegetation has been shown to alter
However, the relationship between land-use change and alteration of flow characteristics, to change the amount of sediment introduced to
stream physical characteristics can be complex and influenced by many stream systems (Walling and Fang, 2003), to decrease infiltration and
factors acting from different spatial scales ranging from local (e.g., increase surface runoff. Frequently, stream concentrations of many
adjacent to the stream) to basin level or to ecoregional scales (Goldstein chemical constituents such as nitrate concentration and phosphorus
compounds are affected both directly and indirectly by human land-
use (Bahar et al., 2008; Quinn et al., 1997). In turn, increases in
⁎ Corresponding author. CONICET Laboratorio de Investigaciones en Ecología y
Sistemática Animal FCN-Universidad Nacional de la Patagonia-Sede Esquel Sarmiento
nutrients can result in overgrowth of algae and aquatic plants, which
849-9200 Esquel-Chubut-Argentina. Tel.: + 54 2945 451165; fax: + 54 2945 452271. potentially alter habitat suitability for endemic fauna (Hall et al.,
E-mail address: [email protected] (M.L. Miserendino). 2001).

0048-9697/$ – see front matter © 2010 Elsevier B.V. All rights reserved.
doi:10.1016/j.scitotenv.2010.10.034
 M.L. Miserendino et al. / Science of the Total Environment 409 (2011) 612–624 613

Channelization and realignment of rivers are also common 2. Materials and methods
practices in urban or agricultural areas, which can produce substrate
modification and disruption of the riffle/pool sequences that provide a 2.1. Study area and site selection
diversity of habitats for aquatic species such as fish and invertebrates
(Jones et al., 1999; Nolan and Guthrie, 1998). The clearing of riparian The study area belongs to the Andean-Humid and Sub-Andean
areas alongside rivers, as modification and/or extirpation of stream- Sub-Humid regions (Del Valle et al., 1998; Paruelo et al., 1998) and is
side vegetation, can also alter the functioning of river ecosystems and located in a transitional mountain and piedmont area in the
disrupt fluxes of organic matter and energy (McCord et al., 2007). Northwest of Chubut province, Argentina between the Chilean border
Moreover this can affect directly or indirectly the river biota. For and the city of Esquel in Northern Patagonia (Fig. 1). Studied basins
example, the loss of oviposition sites and habitat for aerial stages of are similar in size, gradient and geology. From a phytogeographic
insects with aquatic larvae can be critical to preserve “in stream” perspective, the study area is located in the ecotone between the
macroinvertebrate diversity (Bojsen and Jacobsen, 2003). Subantarctic Forest and the Patagonian steppe, and exhibits a marked
There is ample evidence that landscape factors influence aquatic altitudinal gradient. The Subantarctic Forest is constituted by
vertebrate fauna, as for instance landscape characteristics can be perennial (Austrocedrus chilensis, Nothofagus dombeyi and Maitenus
linked to fish species assemblages (Walters et al., 2003). Allan (2004) boaria) and deciduous tree species (N. pumilio, N. antarctica), whereas
has shown significant relationships between some landscape metrics the shrub and herbaceous strata are composed mainly by Chusquea
such as amount of agricultural, forest or urban land cover, and some culeou, Berberis buxifolia, Fuchsia magellanica, Aristotelia chilensis,
biological metrics for fish. At a more local scale, habitat conditions can Oenothera odorata, Fragaria chiloensis and Geranium sp. The valleys are
affect fish distribution as has been observed by Wohl (2000), who dominated by the herbaceous, shrub-like steppe vegetation, Mulinum
reported that some native species in highland stream systems may be spinosum, Stipa speciosa, S. humilis, Corynabutilon bicolor, Verbena
particularly vulnerable to replacement by generalist species in tridens, Nassauvia glomerulosa and Berberis heterophylla (Tell et al.,
response to increased sediment input. 1997).
Fragmentation of continuous habitats (e.g. riparian corridors and The exotic vegetation is composed mainly by forests of Pseudot-
native forest) results in habitat loss and an increase in isolation, which suga menziesii, Pinus ponderosa, P. radiata and P. lambertiana; in many
might affect birds by decreasing colonization, increasing mortality cases these species occupy the riparian corridor replacing most of the
and extinction, and reducing biodiversity, among other potential native species.
effects (Saunders et al., 1991; Wilcox and Murphy, 1985). Changes in Across our study region, we sampled streams impacted by four
agricultural practices also continue to have a wide range of effects on different land-uses: selective harvest of native forests (HNF: managed
terrestrial and aquatic birds species (Duncan et al., 1999; Henderson native forest), plantations of non-native commercially-harvest conifers
et al., 2000), whereas urbanization is considered one of the leading (Pn: pine plantations), grazed pasture (Pas: pasture), and urbanization
drivers of native bird species extirpation at the continental and (Ur: urban). We also sampled two types of reference streams: non-
regional scales (Czech et al., 2000). managed native forests (NF) and reference urban streams (R-ur).
Biodiversity in Patagonian mountains, is globally significant due to Within each of these six categories, we collected samples from three
pronounced endemism, habitat heterogeneity and biogeographic streams. A sampling site was a 100-m reach of the selected stream, for
location (Morrone 2006; Paruelo et al., 1998). Presently there is a bird censuses this area was extended to accomplish the protocol (see
strong concern because land-use practices are changing rapidly the biological sampling).
landscape in the region, being the most widespread activities the Sites on non-managed native Nothofagus forest (NF: Chiquito, Loro
conversion of native forest to pasture, forest exploitation (wood and Comisario) were considered as reference sites to test for the
collection for fuel, carpentry, etc.), and the substitution of native following impacts: conversion of forest to pasture, pine plantation and
forest for pine plantations. These human activities have gained wood extraction. Three sites on streams with native N. antarctica and
attention in recent studies: analyses of water quality and macro- N. pumilio managed forest were selected (HNF: Glyn, Cabeza de Vaca
invertebrate response to impairment (Miserendino and Pizzolón, and Pipo). At these sites the harvest was achieved during the spring
2000; Miserendino et al., 2008); land-use changes (Miserendino, and summer seasons previous to the study. The logging operations
2004; Miserendino and Pizzolón 2004); effect on fish (Di Prinzio et al., implied a selective cut of young to mature trees (dead and live
2009; Lara et al., 2009) and bird communities (Lantschner and Rusch, specimens), which affected no more than 50% of the studied area/
2007; Willson et al., 1994). parcel. Exploitation should exclude specimens from the river corridor
Assessing the impact of land-use activities on river catchments is and operators also have to keep the streams free of debris and
an important subject faced by resource managers. Although stream branches once the tree is logged. However during the survey we
ecologists traditionally have tried to address how land-use changes observed some river corridors having logged trees, and some water-
affect aquatic environments, integrated approaches examining dis- courses obstructed with debris. Most of the properties having wood
turbance through environmental features and biological communities collection permission have a mixed management, thus areas are used
(including aquatic plants, macroinvertebrates, riparian and littoral for extensive livestock which is sustained by the herbaceous stratum.
invertebrates, fish and birds) are not frequent. It is also crucial to To assess pasture impact we selected sites (Pas: Nant y Fall,
regionally calibrate and test existing stream assessment protocols. Manguera and Los Ñires) that had been cleared from the native
To test the hypothesis that greater intensity of land-use will have Nothofagus forest (~60–70 years ago) and had been intensively grazed.
negative effects on water quality, stream habitat and biodiversity, in Pasture sites sustain livestock, mainly cows and sheep. The pine
this work we assess 18 stream sites (north-western Chubut province) plantations sites corresponded to streams running (950–1200 m)
subjected to a gradient of land-use practices (non-managed native through plantations of the exotics P. menziesii, P. ponderosa, P. radiata
forest, managed native forest, pine plantations, pasture, urbaniza- and P. lambertiana (Pn: Golondrinas, Patriada and Ifona) that were
tion). We also evaluate the quality of riparian ecosystems and in mainly established between 1954 and 1970. Three urban rivers, (Ur: Las
stream habitat conditions using two indexes adapted to Patagonian Minas, Carbón and Esquel) plus their respective reference sites (R-ur:
environments and assess the values of these indexes as predictors of Las Minas ref, Carbón ref and Esquel ref) were included in order to assess
the biotic assemblages. Finally, we attempt to identify which specifically the urban impact and make sites comparable in terms of
communities are more responsive to local and large scale changes in environmental features (Fig. 1).
Patagonian mountain streams in order to suggest better management For stream selection we looked at similitude among streams
practices. within a geographical range, but also evaluated the accessibility to the
614 M.L. Miserendino et al. / Science of the Total Environment 409 (2011) 612–624

Fig. 1. Location of sampling sites (site codes in Table 1), Chubut Province, Patagonia, Argentina. Symbols are: ■: urban sites, △: managed native forest, ○: pasture, □: reference urban
sites, ●: non-managed native forest and ▲: pine plantations.

study reach. We also revised existing data sources in different percent coverage of native and exotic vegetation, pasture and urban
governmental administration offices (DPByP Dirección Provincial de settlements.
Bosques y Parques, INTA Instituto Nacional de Tecnología Agrope- At each site, percentages of boulder, cobble, gravel, pebble, and
cuaria). We conducted this study on 18 sites on 15 mountain and sand in the reach were estimated using a 1-m2 grid. Current speed
piedmont streams. Sites were visited in May (autumn), September was measured in mid channel (average of three trials) by timing a
(winter) and December (spring) 2005 and in February (summer) float as it moved over a distance of 10 m (Gordon et al., 1994).
2006, under normal hydrological conditions (avoiding rainstorms or Average depth was estimated from five measurements with a
extremely high discharge events). calibrated stick along one transverse profile across the channel. Wet
and dry widths (from bank to bank) of the channel were also
2.2. Site characterization determined. Discharge was obtained by combining depth, wet width
and current velocity as in Gordon et al. (1994). At each site, air and
In order to identify basin features, vegetation types and coverage, a water temperature were measured with a mercury thermometer.
set of LandSat images ETM (2000 and 2001, mapping resolution 15 m) At each occasion and in the mid channel section of the reach,
were classified and processed using a geographical information specific conductance, pH, turbidity and dissolved oxygen were
system (GIS) (Arc VIEW 3.3 package) by the DPByP. For each surveyed measured with a Horiba U2-probe. For nutrients analyses water
site we obtained sub-basin area above to the sampling site, and the samples were collected below the water surface, kept at 4 °C and
 M.L. Miserendino et al. / Science of the Total Environment 409 (2011) 612–624 615

transported to the laboratory for analysis. Total suspended solids visual inspection over rocks, trees and bushes. Riparian and littoral
(TSS) were measured from separate water samples. Total nitrogen species were identified to the lowest taxonomic level and treated
(TN) and total phosphorus (TP) were determined on unfiltered together in the analyses. Sampling of these communities aimed at
samples digested with persulphate, whereas nitrate plus nitrite increasing the diversity of taxa sampled in the watercourses; they will
nitrogen (NO3 + NO2), ammonia (NH4), and soluble reactive phos- also be used for additional studies not included in this paper (e.g.
phate (SRP) were analyzed using standard methods (APHA, 1994). insects available as prey for fish communities).
Attributes of the riparian vegetation were examined at each site Fish sampling was carried out by electro-fishing with a powered
using an adaptation of the QBR index (Riparian corridor quality index, backpack machine and a hand net. Stream reaches were blocked off
Munné et al., 1998) for Patagonian streams: the QBRp (Kutschker et with a downstream net. Each reach was sampled three times;
al., 2009). This index combines information from total cover, downstream passes of equal effort were done in a slow zigzag
structure, complexity and naturalness of vegetation, and the degree pattern. Specimens were identified to species following López et al.
of channel alteration (e.g. bank modifications, dredging, etc.). The (2003) and counted. Given that an intensive sampling program may
total QBRp score ranges from 0 points (extreme degradation) to 100 affect the fish populations under study (see Casaux and Barrera-Oro,
points (excellent quality, natural riparian forest). 2002), up to 30 individuals per species (selected randomly) were
We evaluated habitat quality (HCI: habitat condition index) using preserved for posterior analyses; the remaining specimens were
the assessment procedure for high gradient streams of Barbour et al. returned alive to the stream. Individuals carried to the laboratory
(1999). This method ranks 10 river channel features (e.g. epifaunal were weighed (W, accuracy 0.01 g). Considering the total of the
substrate availability, frequency of riffles, etc.) from 0 to 20. A score of individuals caught, mean density (ind. m2) and biomass (g. m2) were
200 points indicates that the river is natural and pristine and in its best estimated for each species at each sampling site. For the biomass
possible condition (range: 150–200). This index evaluates the ability estimation, the mass of the individuals returned to the stream was
of the stream's physical habitat to support a given fauna, and then is a assumed to be similar to the mean mass of the individuals carried to
measure of the spatial heterogeneity of the stream (Castela et al., the laboratory.
2008). Taking advantage of the overall sampling effort, we opportunis-
A vegetation sampling was conducted in February of 2006. We tically performed exploratory avian censuses aimed to estimate bird
collected and made vouchers of 60 plants. Species were identified in abundance (number of individuals per point count), richness and
the laboratory using the keys in Correa (1978, 1999). At each reach diversity (Simpson Index) and to look for correlations between these
coverage of aquatic plants and also filamentous algae (eg. Cladophora and the remaining parameters measured throughout the study (see
sp.) were assessed visually as in Hauer and Lamberti (1996). previous discussion). At each stream and during each sampling period
we censused birds at three stations located within the riparian
2.3. Biological sampling corridor and separated by at least 200 m, using 30 m fixed-radius
point counts (0.28 ha) of 10 min duration (the censuses started two
Quantitative macroinvertebrate samples were taken with a minutes after arrival at the station). Due to the overall sampling
Surber sampler (0.09 m2; 250 μm pore size). On each reach, three protocol, censuses were not performed according to the standard
samples from riffles and three from pools (n = 6) were taken. methodology (observations within the three hours after the sunrise or
Samples were fixed in situ with 4% formaldehyde, and sorted in the before the sunset). The individuals observed were identified following
laboratory under at least 5x magnification. Macroinvertebrate the descriptions in Narosky and Yzurieta (1987).
species were identified to the lowest possible taxonomic level
using regional keys (Archangelsky and Manzo, 2007; Fernández 2.4. Statistical analysis
and Domínguez, 2001; Manzo, 2005) and counted.
We calculated a set of macroinvertebrate community descriptors One way ANOVA analyses with land-use as categorical variable in
for each site and sampling date, including richness measures: taxa the treatments followed by post hoc test when significant (Tukey test
richness (SR) and Ephemeroptera, Plecoptera and Trichoptera (EPT) p b 0.05), or non-parametric ANOVA models (Kruskal–Wallis) were
richness; enumerations measures: total density (TD) (Barbour et al., used to assess differences in environmental variables among land-
1999; Resh et al., 1995) and diversity measures: Shannon-Weaver uses. Rivers subjected to certain land-use (pastures, urban, etc.) were
diversity index (H'). A biotic index previously adapted for use in the cases in the models (n = 3). Additionally sampling dates were used as
Patagonian region, BMPS (Biotic Monitoring Patagonian Streams) was cases to improve significance in comparisons, which is acceptable in
calculated (Miserendino and Pizzolón, 1999). The BMPS is an order to contrast multiple independent samples (Kruskal–Wallis)
adaptation of the BMWP (Biological Monitoring Working Party (Sokal and Rohlf, 1995). This included physico-chemical parameters
index; Armitage et al., 1983) and is obtained from a table of 95 (geomorphological measurements, nutrients, etc.), the QBRp and the
families with different degrees of pollution sensitivity (scores 1–10) HCI. The same analysis was performed for community descriptors and
present in Patagonia. The total BMPS score ranges from 0 to N150. metrics. Differences in SR, TD, H´, EPT richness and BMPS were
Riparian and littoral insects were collected using an array of assessed for the benthic community, species richness for riparian and
different techniques. Littoral species are aquatic species associated to littoral invertebrates, fish density and total biomass for the fish
the aquatic vegetation, or found underneath rocks or wood within the community, and species richness and total density for the bird
water; these species usually do not appear in Surber samples. Riparian community.
species are those species that live close to the water, or that can be Prior to analysis data were checked for normality and homogeneity
associated to the marginal vegetation, but are terrestrial. The samples of variances using Kolmogorov–Smirnov and Levene's test, respectively,
were qualitative, and the material was fixed in 80% ethyl alcohol. and log (x+ 1)-transformed when appropriate (Gotelli and Ellison,
Littoral species were collected by passing a D-frame net (Merritt and 2004).
Cummins, 1996) through the aquatic vegetation and also by removing Redundancy Analysis (RDA) was run using CANOCO (ter Braak and
the substrate (rocks, branches, algae, etc.). The content of the net was Smilauer, 1999) to assess relationships between main community
placed in a white pan and inspected visually and the insects caught attributes (riparian and littoral invertebrate richness, macroinverte-
were pulled out with featherweight forceps. Riparian species were brate metrics, density and biomass of fish and bird) and environ-
sampled by washing the substrate in white pans, and also by direct mental variables. RDA was chosen because previous inspection of the
observation. Flying adults of some orders were collected with a data revealed a linear mode rather than a unimodal response in the
sweeping net used through the vegetation on the margins, and also by biotic variables (ter Braak and Smilauer, 1998). All environmental
616 M.L. Miserendino et al. / Science of the Total Environment 409 (2011) 612–624

variables included in Table 1 and 2 were used, initially, to evaluate the 3.2. Aquatic plants
response of community descriptors and sites to environmental
gradients. Variables (except pH) and community attributes (except A total of 15 species of aquatic plants were recorded, with pasture
H´) were transformed (log x + 1), prior to analysis. Variables that and urban sites supporting more species than the rest (Table 4). In
were strongly intercorrelated with others (those with an inflation contrast pine and non-managed native forest sites had fewer species
factor N10) in the initial analysis, were removed and a further analysis mostly represented by Veronica anagallis-aquatica, Juncus balticus and
was carried out with the remaining environmental variables. A Monte Ranunculus sp. Bryophytes were well represented at both managed
Carlo permutation test (9999 permutations) was used to verify the and non-managed native forest sites. Aquatic plant coverage ranged
significance of the models (ter Braak and Smilauer, 1998). between 0 to N50% (Table 1).

3.3. Habitat condition and riparian quality

3. Results
QBRp mean scores showed that urban sites had strong riparian
3.1. Environmental characteristics and land-use ecosystem alteration whereas pasture and pine plantations sites
showed a moderate disturbance (Fig. 2). Sites Manguera (Pas) and
Elevation of sampling sites ranged between 361 and 820 m.a.s.l. Carbón down (Ur) showed the lowest values of riparian vegetation
(Table 1). Stream shade was higher at non-managed native forest sites coverage, having simplified and poor quality communities mostly
(Table 2), which resulted in significantly lower mean values of water represented by; Plantago laneolata, Rumex acestosella, Artemisia
temperature than the rest (ANOVA, p b 0.002, Table 3). Wet channel absinthium and Trifolium repens. Other urban sites showed dominance
width (range 1.1–25 m) varied considerably between sites, with of exotics such as Cytisus scoparia, Rosa rubiginosa and Conium
pasture and urban sites having the widest channels. However maculatum. Non-managed native forest sites showed significantly
differences were only significant when comparing urban and pasture higher values of QBRp than the rest (Kruskal–Wallis, p = 0.017, Table 3).
wet width with all types of forested channels (ANOVA, p b 0.02, HCI index mean values ranged between 111.3 and 182 (Fig. 2).
Table 3). At most sites the substrate particle size was mainly boulder, Urban and pasture sites showed suboptimal habitat quality, possibly
cobble and pebble; however a higher proportion of soft sediments as consequence of channel realignment and dredging which was a
(N30%) was observed on some sites (Table 1), having managed native regular practice at these particular sites. When comparing all values,
significantly higher sand percentage than non-managed native sites scores were significantly higher at non-managed native forest, pine
(Table 3). plantation and reference urban sites than at the rest (Kruskal–Wallis,
Mean conductivity values ranged between 35 and 160.7 μS cm− 1 p = 0.013).
(Table 2), having urban and pasture sites significantly higher
conductivity values than native forest (with and without manage- 3.4. Riparian and littoral invertebrates
ment) (Kruskal–Wallis, p = 0.02).
Except for managed native forest and reference stations for urban We collected and identified 105 riparian and littoral taxa from all
sites, dissolved oxygen contents were significantly higher at non- sites, being Trichoptera (24%), Plecoptera (19%), Ephemeroptera
managed native forest sites (Kruskal–Wallis p b 0.02, Table 3) than at (18%), Diptera (16%) and Coleoptera (16%) the most diverse orders.
the rest. Most chemical variables including SRP, TP, NO3–NO2 and NH4 Riparian and littoral invertebrate richness ranged from 0 to 25 taxa
showed higher mean values at urban use than at the rest of the sites when taking into account absolute values, whereas mean values per
(Table 2). Except for reference stations in urban sites, NO3–NO2 values site ranged from 0 to 15 taxa. Significant differences on taxa richness
were significantly higher at urban sites than at the rest (Kruskal– were found among land-uses (Kruskal–Wallis p b 0.008, Fig. 3a)
Wallis p b0.027, Table 3). However NH4 and SRP values were highly having urban sites significantly lower richness than the rest, and
variable, which might be explained by extreme values at Esquel down pasture sites significantly higher richness than pine, and managed
urban site. native forest. This pattern was mostly explained by the presence of

Table 1
Physical characteristics of 18 sampled sites of streams in Patagonia, Argentina. Land-use codes: R-ur: reference for urban site; Ur: urban site; Pas: pasture; Pn: pine plantation; HNF:
managed native forest; NF: non-managed native forest. Bo: boulder, Co: cobble, Pe: pebble, Gra: gravel, Sa: sand.

Site Site Land- Elevation Catchment size % Native % Pasture % Exotic % Urban % Timberline Substrate % % Aquatic plant
Code use (m.a.s.l) (km2) forest wetlands forest areas and lakes size Sand coverage

Esquel EU R-ur 721 22.9 11.3 0 0.7 0 88.0 Bo/co 10 b5

Las Minas LMU R-ur 573 35 93.1 1.3 0 0 5.6 Bo/co 10 b5
Carbón CU R-ur 384 80.9 97.0 1.4 1.4 0.2 0 Bo/co 10 b5
Esquel ED Ur 491 299 11.0 18.7 5.6 9.0 55.7 Pe/gra 20 N50
Las Minas LMD Ur 555 35.5 92.3 1.3 0 0.7 5.7 Co/pe 5 b5
Carbón CD Ur 403 110 96.6 1.5 1.1 0.4 0.4 Bo/co 10 b5
Nant y Fall NyF Pas 690 161.8 69.7 4.4 0.2 0 25.7 Gra/co 15 25
E. Los Ñires LÑ Pas 667 102.9 81 8 0 0 11 Sa/gra 40 b5
Manguera MG Pas 699 20.6 76 b1 0 0 23 Bo/co 5 b5
Golondrinas GOL Pn 361 5.7 63 0 5.3 0 31.7 Co/sa 30 0
Patriada PAT Pn 495 1.8 79.1 0 7.3 0 13.6 Gra/sa 35 0
IFONA IFO Pn 368 35.5 69.5 0 3.1 0 27.4 Bo/co 5 0
Cab. de vaca CVA HNF 699 88.7 90.4 2.8 0 0 6.8 Co/gra 15 15
Pipo PIP HNF 728 2.2 99.0 0 0 0 1.0 Gra/sa 35 20
Glyn GLY HNF 615 21.7 98.8 0 0 0 1.2 Pe/co 15 20
Comisario COM NF 770 21.8 62.8 0 0 0 37.2 Bo/co 5 15
Loro LO NF 820 8.2 43.7 0 0 0 56.3 Bo/co 5 35
Chiquito CHI NF 670 11.1 97.1 0 0 0 2.9 Co/pe 10 15
 M.L. Miserendino et al. / Science of the Total Environment 409 (2011) 612–624 617

Table 2
Means and ranges (in brackets) of environmental quality variables in the study streams during the one-year study period (n = 12).

Land-use: Ref. urban Urban Pine plantation Pasture Managed native Non-managed
forest native forest
Variable:

Physicochemical variables
Water temp. (°C) 8.6 (5.1–10.2) 10 (5.7–11.9) 8.1 (7.1–11.5) 8.2 (3.4–13.6) 7.1 (4–10.6) 4.5 (1.7–7.4)
Wet width (m) 5.3 (2.1–9.1) 11.5 (3.8–23.5) 4.1 (1.1–9.2) 11.3 (5–25) 4.3 (1.3–11) 4.7 (1.5–8.9)
Depth (cm) 21 (9.2–31.5) 23 (7–38.2) 22 (7.8–41.3) 30 (13.5–33.7) 18 (5.2–28) 25 (10.4–34.3)
W. velocity (m s− 1) 0.84 (0.33–1.36) 0.91 (0.33–1.81) 0.59 (0.20–1.07) 0.94 (0.66–1.11) 0.81 (0.10–1.66) 0.92 (0.10–1.66)
Discharge (m3 s− 1) 1.02 (0.12–2.83) 2.83 (0.15–8.97) 1.09 (0.02–3.8) 3.60 (0.45–8) 1.03 (0.01–3.08) 1.36 (0.02–4.7)
pH 7.33 (6.7–7.7) 7.23 (6.7–7.7) 7.16 (6.8–7.5) 7.38 (6.6–7.9) 7.28 (7–7.6) 7.28 (7–7.6)
Conductivity 20 ºC (μS. cm−¹) 97.6 (44–186) 160.7 (52–390) 79.8 (29–137) 87.3 (61–118) 61.3 (38–97) 35 (9–92)
Dis. oxygen (mg l−¹) 11.5 (9.2–18) 10.5 (7.9–17.1) 11.6 (8.6–15.7) 11.1 (7.7–14.9) 12.7 (8.9–18.8) 14.1 (11–20)
Turbidity (NTU) 25.1 (0–158) 48.2 (1–287) 33.7 (0–240) 9.7 (2–32) 7.6 (0–37) 30.5 (0–192)
TSS (mg l−¹) 7.56 (0.1–52.3) 14.75 (0.7–102.7) 2.15 (0.27–5.59) 5.29 (1.08–16.1) 4.24 (0.67–33.9) 1.97 (0.1–6.1)
SRP (ug l−¹) 0.62 (0.36–1.17) 4.88 (0.29–20.96) 0.50 (0.19–1.24) 0.72 (0.32–1.79) 0.72 (0.21–2.65) 0.56 (0.16–1.88)
Nitrate/nitrite (ug l−¹) 2.34 (0–17.44) 16.47 (0.03–69) 0.42 (0–1.61) 0.16 (0–0.60) 0.22 (0.07–0.62) 0.20 (0–1.55)
Ammonia (ug l−¹) 1.23 (0–3.26) 35.51 (0.49–271.5) 0.96 (0.08–1.58) 1.15 (0.38–2.49) 1.28 (0.14–4) 1.20 (0.17–3.89)

Sigara sp., Notonecta sp., Belostoma sp. (Hemiptera, Heteroptera), 3.6. Fish communities and land-use
Rionaeschna sp. and Cyanallagma interruptum (Odonata), Rhantus
signatus, Liodessus patagonicus and Gymnochthebius sp. (Coleoptera), Five fish species were represented in the samples. Three of them
which were frequently recorded at these sites. (Oncorhynchus mykiss, Salmo trutta and Salvelinus fontinalis) were
exotic species (99.1% of the individuals caught) and the remaining
(Hatcheria macraei and Odontesthes hatcheri) were native ones
3.5. Benthic macroinvertebrates and land-use relationships (Table 5). Whereas exotic species were represented at all of the
land-uses, native fish only occurred in reference urban, pastures and
A total of 158 taxa were recorded in the study; macroinvertebrate managed native sites. Overall, fish density ranged between 0 and 2.69
density and number of taxa ranged from 298 to 38,831 individuals ind. m− 2 (Table 5). Total fish density did not vary in a systematic
m− 2 and from 11 to 54 taxa site− 1 respectively (Fig. 3c and d). manner among land-uses (Kruskal–Wallis test, p = 0.55, Fig. 3g). Fish
Number of taxa differed among land-uses having non-managed native biomass ranged between 0 and 54.0 g m− 2 (Table 5). Total biomass
forest sites higher richness than the rest with the exception of was significantly higher at urban than at pine sites (Tukey test
managed native forest sites (p = 0.002, Fig. 3c); an identical pattern p b 0.05, Fig. 3h).
was displayed by EPT richness (p = 0.003, Fig. 3e). Mean macro-
invertebrate density was significantly higher at urban sites except for 3.7. Bird community and land-use
managed native sites (p = 0.001, Fig. 3d). Non-managed native forest
sites had significantly more macroinvertebrate diversity than the rest A total of 34 species were identified throughout the study, but only
of land-uses, and as expected urban sites showed the lowest diversity, 8 of them (23.5%) were observed actively interacting with the streams
being significantly lower than reference urban and pine sites (Table 6). The highest number of species (26) and mean annual
(p b 0.001). Analysis of BMPS showed that only urban sites fall into diversity (1.46) were observed at urban sites. Maximum density (4.96
lower categories of water quality (Mean score = 84.9). Non-managed individuals per point count) was observed at pasture sites (Table 6).
native sites (Mean score = 163) showed very clean waters, whereas Differences observed in bird diversity and density among land-uses
pine sites (Mean score = 135) displayed clean waters (Fig. 3f). were not statistically significant (Kruskal–Wallis p = 0.16 and
Managed native sites showed clean waters though BMPS values p = 0.13 Fig. 3i and j).
were significantly lower than reference sites (p b 0.0001).
3.8. Ordination analysis
Table 3
Results of the Kruskal–Wallis ANOVA by rank analyses (non-parametric) among First two axes of RDA analysis accounted for 56.1 of the total
environmental variables and land-use classes. n = 3 for QBRp and HCI, n = 12 for the variance in the community data (Fig. 4). The 1st axis was significant
rest of variables. ns: non significant. Land-use codes: R-ur: reference for urban site; Ur: (p b 0.005), and represented an environmental gradient mainly
urban site; Pas: pasture; Pn: pine plantation; HNF: managed native forest; NF: non- defined by land-use coverage and water quality conditions whereas
managed native forest.
the 2nd axis was associated to pasture land cover and oxygen contents
Variable p Relationship observed (Table 7). According to the model landscape variables such as urban
Water temperature 0.002 Ur N HNF,NF Pas, Pn, R-ur, HNF N NF percentage and pasture land cover decreased from the positive to the
Wet width 0.0052 Pas,Ur N rest negative end of the RDA1, whereas percent of native forest cover,
Depth ns increased to the negative end of same axis. Urban sites were
Current velocity ns
positioned in the lower right quadrant whereas non-managed native
Discharge ns
% Sand 0.04 HNF N NF forested sites were in the lower left quadrant (Fig. 4 and Table 7). A
Dissolved oxygen 0.04 NF N Ur,Pas,Pn gradient in water quality conditions including dissolved oxygen,
pH ns nitrate and TSS contents was also highlighted in the model and
Turbidity ns associated with RDA1 (Fig. 4 and Table 7). Sites having better water
TSS ns
Conductivity 0.004 Pas,Ur N HNF N NF
quality conditions were positioned at the negative extreme of RDA1
NH4 ns and polluted sites at the positive extreme of RDA1. Sites having
NO3 0.026 Ur N Pas, HNF, NF, Pn physical disturbances such as dredging, realignment of the channel
SRP ns and vegetation removal (which in turn produced high TSS levels)
QBRp 0.017 NF N rest HNF N Pn N Ur
were positioned in the lower right quadrant. Community descriptors
HCI 0.013 NF,Pn,R-ur N Pas,HNF,Ur
were located along the same gradients, thus H´, invertebrate riparian
618 M.L. Miserendino et al. / Science of the Total Environment 409 (2011) 612–624

Table 4
Aquatic vegetation present in 18 sampling sites in Patagonian streams (Chubut, Argentina) according to selected land-use during the study (2005–2006). Nf: no found.

Reference urban Urban Pine plantations Pasture Managed native forest Non-managed native forest

Veronica anagallis-aquatica V. anagallis-aquatica Nf V. anagallis-aquatica V. serpyllifolia J. arcticus

Juncus microcephalus V. serpyllifolia V. anagallis-aquatica Ranunculus sp.
J. arcticus J. arcticus M. glabratus Bryophytes
Myriophyllum quitense J. burkartii R. hydrophilus
Callitriche lechleri J. diemii C. lechleri
Ranunculus sp. J. microcephalus E. pseudoalbibracteata
Mimulus glabratus Eleocharis albibracteata Bryophytes
Eleocharis sp. E. melanostachys
Cladophora sp. C. lechleri
Ranunculus flagelliformis
Limosella australis
Lilaeopsis macloviana
M. glabratus
Cladophora sp.

richness, BMPS, EPTr and macroinvertebrate richness increased conductivity, nutrients and oxygen. On the other hand, stream-bank
significantly with percentage of native cover and better habitat erosion was an obvious feature in managed native sites as well as in
conditions. Instead, macroinvertebrate density augmented with urban some pasture sites with an extreme situation at Los Ñires; changes in
land cover percentage, whereas fish biomass, bird density and bird the riparian vegetation, together with the direct effects of stock
richness were positively correlated with the percentage of pasture trampling, probably accounted for an increase of fine material in the
cover. streambed. Stocking with sheep or cattle in the absence of streamside
fences or riparian strips can lead to bank erosion, bed disturbance, and
increased supply of fine sediment (Quinn et al., 1997).
4. Discussion
We detected significant increases in water temperature and
conductivity values in harvested forest sites when compared with
4.1. Environmental variables and land-use
non-harvested sites, meaning that wood collection is affecting the
environmental features. Baillie et al. (2005) documented marked
A considerable number of physical, chemical, and environmental
increases in stream light levels and maximum monthly water
changes associated with different land-use practices were documen-
temperatures in streams draining through 25% clear-cut forests in
ted in this study. A reduction of the riparian canopy resulted in an
New Zealand watersheds and Kedzierski and Smock (2001) reported
increase of the water temperature and sedimentation symptoms,
important differences in incident light reaching streambeds when
which diminished the quality of the stream habitat. As reported in
compared with unlogged and logged sections in a lowland creek.
other studies (Collier, 1995; Miserendino et al., 2008) water quality
variables showed that urban sites were the most affected in terms of
4.2. Riparian and littoral invertebrate richness pattern

100 Very good quality Our results suggest that urban sites showed lower richness than
the remaining uses; human impacts observed in urban sites include
75 Good quality afforestation, river canalization and pollution, mostly by domestic
QBRp score

Riparian ecosystem sewage. These impacts help to reduce taxonomic diversity (Maitland
alteration and Morgan, 1997). Our observations are in agreement with what
50
was showed by Subramanian et al. (2005). These authors reported
Strong riparian
ecosystem that generic richness of insects was higher in the streams with
25 alteration natural riparian vegetation than in those modified by humans; they
Extreme degradation also found that the diversity of riparian grasslands is comparable to
of riparian ecosystem
that of streams with riparian paddy fields. Our forested sites had an
0
intermediate richness of riparian and littoral invertebrates; this
R-ur Ur Pas Pn HNF NF
could probably be explained by the absence of marginal areas with
200 slow moving water with macrophytes and algae, which harbor an
important number of species. Instead, our pasture sites supported
Optimal
more species of riparian and littoral invertebrates than other uses.
150
HABITAT score

The taxa occurring at these sites were characteristic of marginal

Suboptimal pools (Merritt and Cummins, 1996), sometimes with poor or slow
100 flow and with macrophytes and algae, biotopes that were present
mostly at our pasture sites. Some of the species in these marginal
Marginal
areas are also opportunistic (e.g. Tropisternus setiger and Rhantus
50
signatus colonize all kind of stagnant or slow moving waters). The
Poor species assemblages of riparian and littoral invertebrates at pasture
0 sites seem to be more diverse because they also incorporate adult
R-ur Ur Pas Pn HNF NF stages of aquatic insects, for example many of the aquatic beetles
and bugs collected in pastures have both larval and adult stages
Fig. 2. Habitat condition index (upper) and riparian ecosystem condition score (lower)
per land-use during the study. R-ur: reference for urban sites; Ur: urban; Pas: pastures;
aquatic (e.g. Tropisternus setiger, Enochrus sp., Rhantus signatus,
Pn: Pine plantation; HNF: managed native forest; NF: non-managed native forest Liodessus patagonicus, Gymnochthebius sp., Notonecta sp., Sigara sp.,
(n = 3). etc.).
 M.L. Miserendino et al. / Science of the Total Environment 409 (2011) 612–624 619

4.3. Macroinvertebrate environmental relationships percentage of urban cover are fully consistent with those reported by
Strayer et al. (2003) and Morley and Karr (2001) who documented a
We found several significant relationships among macroinverte- significant relationship between macroinvertebrate species richness
brate attributes and water quality variables, land-use coverage types, and urban land cover for North American watersheds.
and riparian and habitat condition indexes. Most metrics (SR, EPT Regarding water quality parameters, the strong relationships
richness and BMPS) were positively correlated with native forest observed between macroinvertebrate metrics with dissolved oxygen,
coverage and negatively correlated with urban and pasture land nitrate and TSS confirmed the ability of benthos community to detect
cover. These observations are in agreement with Collier (1995, 2008) stream impairment. Nitrate and ammonia are strong indicators of
findings for New Zealand streams and with Roy et al. (2003) for organic pollution (Hall et al., 2001; Paul and Meyer, 2001),
streams in Georgia (USA). They reported positive relationship consequently significant reductions in SR, EPT richness, H, BMPS
between total richness, EPT richness and the percentage of indigenous occurred with increases of nitrogen concentrations at urban sites. We
forest, but also a negative correlation with the increase of pasture can conclude that these metrics are useful to detect impairment
cover. The relationships between the metrics considered and the especially at urban reaches in Patagonia.

20 3,0
NF>all Ur< R-ur, Pn, NF
RIP./LIT. INVERT. RICHNESS

BENTHOS DIVERSITY (H)

Pas> Pn, HNF, Ur p=0.001
15 Ur<all
p=0.008 2,5

10

2,0
5

0 a 1,5
b
R-ur Ur Pn Pas HNF NF R-ur Ur Pn Pas HNF NF
LAND USE LAND USE

50 20000
NF> R-ur, Ur, Pn, Pas Ur>R-ur, Pn, Pas, NF
HNF> Ur p=0.001
BENTHOS S. RICHNESS

BENTHOS DENSITY

p=0.002 15000
40

10000

30
5000

c d
20 0
R-ur Ur Pn Pas HNF NF R-ur Ur Pn Pas HNF NF
LAND USE LAND USE

25
NF>R-ur, Ur, Pn, Pas NF>R-ur, Ur, Pas, HNF Very
HNF>Ur Pn>R-ur, Ur, Pas p=0.0001 clean
p=0.003 150 water
20
EPT RICHNESS

BMPS

15 Clean
water
100
Incipient pollution or other
10
kinds of perturbation

e f
5 50
R-ur Ur Pn Pas HNF NF R-ur Ur Pn Pas HNF NF
LAND USE LAND USE

Fig. 3. Distribution values of selected metrics for a) riparian and littoral invertebrate richness, b) benthic diversity (Shannon-Weaver), c) macroinvertebrate species
richness, d) macroinvertebrate total density (ind.m− 2 ), e) EPT richness, f) BMPS quality index, g) fish density (ind.m− 2 ), h) fish biomass (g.m− 2 ), i) bird diversity
(S: Simpson index) and j) bird density per land-use (ind. per point count). Range bars show maxima and minima, boxes are interquartile ranges (25–75%), small
squares are medians. Significant relationships (ANOVA Kruskal–Wallis, Tukey test) are shown in bold. Land-uses: R-ur: reference for urban sites; Ur: urban; Pas:
pastures; Pn: Pine plantation; HNF: managed native forest; NF: non-managed native forest. Sample size was n = 12 except for item a) n = 9.
620 M.L. Miserendino et al. / Science of the Total Environment 409 (2011) 612–624

1,2 1,2
Ur>Pn
p=0.55 p=0.039

FISH BIOMASS (log x+1)

FISH DENSITY

0,6 0,6

g h
0,0 0,0
R-ur Ur Pn Pas HNF NF R-ur Ur Pn Pas HNF NF
LAND USE LAND USE
3 8

p=61.0 p=31.0
BIRD DIVERSITY (S)

BIRD DENSITY
2

1
2

i j
0
R-ur Ur Pn Pas HNF NF R-ur Ur Pn Pas HNF NF
LAND USE LAND USE

Fig. 3 (continued).

Even though impairment was only evident at urban sites when heterogeneity of the stream (Castela et al., 2008) and macroinverte-
using the BMPS biotic index, this revealed that several sensitive brates diversity (H´) was also evidenced in the ordination. Subsequently,
families (e.g. Ameletopsidae, Diamphipnoidea, Nesameletidae, Ecno- native forest sites having the best ecological conditions (highest QBRp
midae, etc.), well represented in native forest and pine sites, were and HCI scores) supported the most diverse communities whereas sites
practically absent at pasture, harvested native and urban sites. As with incipient problems of sedimentation had less diverse communities
anticipated by Miserendino and Pizzolón (1999) BMPS appears composed by tolerant taxa. Physical disturbances such as dredging,
promising to assess disturbances involving sedimentation processes realignment of the channel and vegetation removal probably accounted
and habitat impoverishment. for the high TSS values and low HCI scores at sites with extreme positions
As indicated by Roy et al. (2003), the present paper highlights the role in RDA (Las Minas down, Los Ñires, Pipo).
of the riparian buffers in preserving water quality. This was evidenced in In our study, the reduction in H' and EPT richness was marked and
the RDA analysis performed, where a clear gradient in riparian significant (pair wise comparisons) when contrasting harvested
ecosystem conditions, habitat quality, NO3−NO2 and TSS contents was native forest with native forest. Bojsen and Jacobsen (2003) found
shown. A high relationship between HCI which is a measure of the spatial that diversity index values were greater in forested sites compared to

Table 5
Mean annual values (± SD) of fish species density (ind. m− 2) and biomass (g. m− 2) per land-use in the study streams during the one-year study period (n = 12). Empty cells indicate
that the species was not recorded. N: cumulative number of individuals in the entire study. D: density B: biomass.

Land-use: Reference urban Urban Pine plantation Pasture Managed native forest Non-managed native forest N

Salmonidae:
Oncorhynchus mykiss D 2.34 ± 2.45 1.64 ± 2.05 0.13 ± 0.22 0.37 ± 0.19 0.63 ± 0.58 0.97 ± 1.22 2113
B 17.01 ± 17.83 32.67 ± 35.82 0.44 ± 0.76 15.86 ±16.31 9.37 ± 9.19 18.99 ± 25.58
Salmo trutta D 0.13 ± 0.22 0.20 ± 0.35 0.02 ± 0.04 0.25 ± 0.43 247
B 1.29 ± 1.74 2.62 ± 4.54 0.33 ± 0.58 2.89 ± 5.00
Salvelinus fontinalis D 0.25 ± 0.43 31
B 6.38 ± 11.05
Trichomycteridae
Hatcheria macraei D 0.01 ± 0.02 0.02 ± 0.02 0.02 ± 0.02 31
B 0.01 ± 0.02 0.08 ± 0.09 0.12 ± 0.11
Atherinidae
Odontesthes hatcheri D 0.002 ± 0.002 4
B 0.003 ± 0.01
Total density (min–max) TD 0.01–2.69 0.02–2.63 0–0.28 0.01–0.54 0–0.69 0–1.17
Total biomass (min–max) TB 0.13–21.83 1.42–54.00 0–8.37 0.06–31.14 0–13.86 0–20.85
 M.L. Miserendino et al. / Science of the Total Environment 409 (2011) 612–624 621

Table 6
Mean annual diversity (± SD) and density (overall and by species; individuals per point count ± SD) per land-use of the bird species observed in the studied streams during the one-
year study period (n = 12). The species with (*) corresponds to those observed actively interacting with the streams. S: Simpson index. Code of uses in Table 1.

R-ur Ur Pn Pas HNF NF

Diversity (S) 1.20 ± 0.7 1.46 ± 1.1 0.91 ± 0.8 1.39 ± 1.1 1.31 ± 0.7 0.58 ± 0.5

Density 2.01 ± 1.2 4.94 ± 4.3 1.67 ± 1.7 4.96 ± 5.8 2.37 ± 1.4 0.91 ± 0.9

Threskiornithidae
Theristicus caudatus (Black-faced Ibis) – 0.03 ± 0.1 – 0.03 ± 0.1 – –
Anatidae
Chloephaga picta * (Upland Goose) – – – 0.11 ± 0.4 – –
Anas flavirostris * (Speckled Teal) 0.06 ± 0.2 – – 0.17 ± 0.6 0.11 ± 0.3 –
Accipitridae
Circus cinereus (Cinereous Harrier) – 0.03 ± 0.1 – – – –
Falconidae
Milvago chimango (Chimango Caracara) – 0.03 ± 0.1 0.03 ± 0.1 0.08 ± 0.2 0.03 ± 0.1 –
Falco sparverius (American Kestrel) – 0.03 ± 0.1 – 0.03 ± 0.10 – –
Charadriidae
Vanellus chilensis (Southern Lapwing) 0.06 ± 0.2 0.06 ± 0.1 – 0.33 ± 0.6 – –
Scolopacidae
Gallinago gallinago (Common Snipe) – 0.03 ± 0.1 – 0.03 ± 0.1 – –
Columbidae
Zenaida auriculata (Eared Dove) 0.03 ± 0.1 – – 0.06 ± 0.2 – –
Trochilidae
Sephanoides galeritus (Green-backed Firecrown) – 0.03 ± 0.1 – – – –
Alcedinidae
Ceryle torquata * (Ringed Kingfisher) 0.03 ± 0.1 0.06 ± 0.1 – – 0.14 ± 0.2 –
Picidae
Colaptes pitius (Chilean Flicker) 0.06 ± 0.2 – – 0.06 ± 0.2 – 0.03 ± 0.1
Picoides lignarius (Striped Woodpecker) – – – – 0.03 ± 0.1 –
Furnariidae
Cinclodes patagonicus * (Dark-bellied Cinclodes) 0.11 ± 0.3 0.17 ± 0.6 – 0.08 ± 0.29 0.03 ± 0.1 –
Cinclodes fuscus * (Bar-winged Cinclodes) – 0.03 ± 0.1 – – 0.03 ± 0.1 –
Cinclodes oustaleti * (Grey-flanked Cinclodes) – 0.14 ± 0.4 – – – –
Aphrastura spinicauda (Thorn-tailed Rayadito) – 0.11 ± 0.4 0.50 ± 0.9 0.03 ± 0.1 0.42 ± 0.5 0.11 ± 0.3
Rhinocryptidae
Scelorchilus rubecula (Chucao Tapaculo) 0.06 ± 0.2 – 0.14 ± 0.2 – 0.03 ± 0.1 –
Tyrannidae
Xolmis pyrope (Fired-eyed Diucon) 0.03 ± 0.1 0.17 ± 0.3 – 0.03 ± 0.1 0.06 ± 0.1 –
Lessonia rufa (Rufous-backed Negrito) – 0.06 ± 0.2 – 0.28 ± 0.5 – –
Anairetes parulus (Tufted Tit-tyrant) 0.08 ± 0.3 0.06 ± 0.2 0.03 ± 0.1 – – 0.03 ± 0.1
Elaenia albiceps (White-crested Elaenia) 0.25 ± 0.5 0.19 ± 0.5 0.17 ± 0.4 0.17 ± 0.5 0.33 ± 0.6 0.11 ± 0.3
Hirundinidae
Tachycineta leucopyga * (Chilean Swallow) 0.17 ± 0.6 1.61 ± 3.8 – 1.94 ± 4.3 0.17 ± 0.6 0.03 ± 0.1
Notiochelidon cyanoleuca * (Blue-and-White Swallow) – 0.42 ± 1.4 – – – –
Troglodytidae
Troglodytes aedon (House Wren) 0.17 ± 0.2 0.08 ± 0.3 0.17 ± 0.3 0.11 ± 0.2 0.19 ± 0.5 0.42 ± 0.6
Turdidae
Turdus falcklandii (Austral Thrush) 0.14 ± 0.3 0.50 ± 0.9 0.03 ± 0.1 – 0.09 ± 0.1 0.08 ± 0.3
Motacillidae
Anthus correndera (Correndera Pipit) – – – 0.08 ± 0.3 – –
Vireonidae
Passer domesticus (House sparrow) 0.03 ± 0.1 0.06 ± 0.2 – – – –
Emberizidae
Diuca diuca (Common Diuca-Finch) 0.06 ± 0.2 – – – 0.03 ± 0.1 –
Phrygilus patagonicus (Patagonian Sierra-Finch) 0.08 ± 0.2 0.11 ± 0.4 – – – –
Zonotrichia capensis (Rufous-collared Sparrow) – – – 0.06 ± 0.2 – –
Fringilidae
Carduelis barbata (Black–chinned Siskin) 0.17 ± 0.6 0.36 ± 0.8 0.12 ± 0.4 – – –
Icteridae
Molothrus bonariensis (Shiny Cowbird) – – – – 0.11 ± 0.3 –
Sturnella loyca (Long-tailed Meadowlark) – 0.08 ± 0.3 – 0.03 ± 0.1 – –
Passer. Unidentified * 0.03 ± 0.1 0.11 ± 0.4 – – – 0.03 ± 0.1

deforested sites, and Hall et al. (2001) reported lower macroinverte- some large epibenthic taxa (Molineri, 2008), reducing individuals
brate diversity in pasture streams. Baillie et al. (2005) mentioned that body size and modifying the functional structure (Buria et al. 2007).
harvesting effects can vary depending on stream type and harvesting Thus, it is possible that some of the changes detected in benthic
practices, particularly post-harvesting wood treatment, but usually communities in our study can be related to some extent with trout
there is a decrease in relative abundance of sensitive taxa (EPT) as is impact on receiving ecosystems.
reported here. A combination of environmental modifications such as
higher temperatures, more sediment in the bottom channel, and 4.4. Usefulness of habitat assessment and riparian conditions indexes
changes in the organic material supply, could account for the changes
observed in these sensitive insect communities (Richardson, 2008). In relation to the QBRp index, site conditions varied from poor to
Introduced salmonids, which were present at most of the sampling very good. The lower QBRp values were registered at urban and
sites, usually affect benthic communities by diminishing the density of pasture sites, confirming that these land-uses negatively impact rivers
622 M.L. Miserendino et al. / Science of the Total Environment 409 (2011) 612–624

successful expansion of salmonids in Patagonia which, at least for

Argentina, is supported by a continuous stocking. Oncorhynchus
mykiss was the most abundant fish both in density and biomass. Its
highest values of density and biomass where observed at reference
urban and urban sites respectively. This agrees with studies reporting
that salmonids can occupy and compete against other species in
altered ecosystems (Boët et al., 1999; McKinney, 2002) and the ability
of O. mykiss to grow faster and to larger sizes at urban sites,
presumably due to higher water temperatures and food availability
(Paul and Meyer, 2001).
Differently to the observed in salmonids, native species were
restricted to high altitude sites located in the ecotone region between
the Andes and the steppe. This scenario seems to be evidencing a
process of geographical segregation between exotic and native
species (Di Prinzio et al., 2009; Soto et al., 2006). However, given
Fig. 4. Triplot of the sampling sites and community descriptors in relation to the first
ordination axes of RDA based on land cover percentages, water quality parameters and
that the sampling sites located in the ecotone region were affected by
reach scale variables as environmental variables. Symbols are described in Fig. 1. two of the less aggressive land-uses considered in this study (pasture
and managed native forest) whereas most of the sites where native
fish were absent were associated to more disturbing land-uses (urban
in Northwest Chubut. The QBRp allowed us to assess not only the and pine plantations), both processes (i.e. changes in land-uses and
ecological integrity of riparian buffer zones, but also to identify strong competence between native and exotic fish) seem to be interacting
relationships between the status of those zones and some biotic and more studies are required to better understand their real effects
elements (macroinvertebrate metrics). By using the QBR index and on fish communities (Di Prinzio et al., 2009).
the fluvial habitat index Castela et al. (2008) detected that the benthic
community responded negatively to increases in nutrient concentra- 4.6. Bird patterns
tion, changes in the riparian vegetation and to the reduction of habitat
quality. Riparian plant cover can control local processes such as Only 23.5% of the species represented in the censuses interacted
shading, loading of leaf litter and wood into the channel, provision of actively (feeding or swimming while looking for food) with the water
egg-laying or adult emergence sites (Danger and Robson, 2004; Hall courses, which might be explained by the low order of the streams
et al., 2001; Richardson, 2008). studied. Consequently, overall bird density and diversity were not
In our work, the HCI did not show the same picture as the QBRp; directly correlated with the stream physical and chemical variables
when using HCI exotic pine plantations did not differ from non- considered. Bird density and diversity only responded in a significant
managed native forest sites. The presence of riparian forest, whether way to the percentage of pasture cover in the riparian corridor. Birds
native or exotic, would help to structure the in-stream habitat. display a high dispersal capacity thus their communities at the study
However, the quality of the litter supply (Pinus vs Nothofagus) can lead streams seem to be structured by environmental characteristics at a
to changes in the functional attributes of macroinvertebrate commu- scale larger than the one here considered.
nities as for example the shredder richness (Miserendino and Masi, The density and diversity were higher in urban and pasture sites
2010). where the riparian forest was strongly altered and almost absent, and
lower at pine and non-managed native forest sites where the riparian
4.5. Fish patterns corridor was densely forested. This might be explained by the space
available to transit within the sites given that this parameter
Exotic fish dominated in density and biomass at all the sampling positively interacts with bird abundance. This pattern is also explained
sites. This is in line with the reports for Patagonia, both in Chile (Soto by an increase in the presence of steppe and ecotone birds as the
et al., 2006) and Argentina (Pascual et al., 2002) and reflects the disturbance in the riparian forest increases, a phenomenon also
observed by Lantschner and Rusch (2007) in North Patagonian forests.
In this sense, studies carried out elsewhere reported that bird species
Table 7
Eigenvalues and correlation of standardized environmental variables with the first two
richness and diversity are strongly and positively correlated with the
RDA axes. The species-environmental correlations scale the strength of the relationship structural complexity of vegetation (MacArthur, 1964; Pearson and
between community descriptors and environment for the axes. F-ratio statistics are Ralph, 1978). However, our results are fully in line with those obtained
listed for the first axis and for all the axes combined. at Northern Patagonia where the species diversity and richness was
Axis RDA1 RDA 2 higher at structurally simple environments (Christie et al., 2004; Ralph,
1985; Vuilleumier, 1972).
Eigenvalue 0.463 0.098
Species-environmental correlations 0.95 0.77 The overall bird abundance, richness and diversity do not help to
Cumulative percentage variance assess land-use effects on stream systems at the study area. For such
Of species data 46.3 56.1 purposes, future studies should identify specific indicator species,
Of species-environmental relation 67.6 81.9 something not possible with our sampling protocol.
Correlations
% Urban 0.87 − 0.15
% Native − 0.37 0.29 5. Conclusions
% Pasture 0.78 0.43
Habitat condition − 0.57 0.10 According to our integrated approach urban land-use produced
QBRp − 0.46 0.21
the most significant changes in aquatic ecosystems including physical
TSS 0.69 − 0.16
NO3 0.80 − 0.05 features, nutrients, habitat conditions, riparian quality, riparian and
Dissolved oxygen − 0.74 − 0.36 littoral invertebrate and macroinvertebrate communities. This is
F-value and significance, Montecarlo test.
significant considering that percentage of urban areas found in this
First axis: 7.768 p b 0.005. study was very low. In an intermediate situation were pasture, pine
All axis: 2.451, p b 0.001. plantations and managed native forest land-uses. As expected, native
 M.L. Miserendino et al. / Science of the Total Environment 409 (2011) 612–624 623

forest sites showed the highest integrity of ecological conditions, Christie M, Ramilo E, Bettinelli M. Aves del Noroeste Patagónico: Atlas y Guía. Buenos
Aires: L.O.L.A; 2004.
which in turn resulted in higher values of QBRp and HCI, and also in a Collier KJ. Environmental factors affecting the taxonomic composition of aquatic
greater biodiversity of the local benthic communities (richness, EPT macroinvertebrate communities in lowland waterways of Northland, New Zealand.
richness). However, by offering stream edge areas, pasture sites still N Z J Mar Fresh 1995;29:453–65.
Collier K. Temporal patterns in the stability, persistence and condition of stream
supported richer communities of riparian and littoral invertebrates, macroinvertebrate communities: relationships with catchment land-use and
which increased overall biodiversity. regional climate. Freshwater Biol 2008;53:603–16.
Of all the biological communities used to evaluate the impact of Correa MN. Flora Patagónica, Tomo VIII, parte III: Gramineae. Buenos Aires: INTA; 1978.
Correa MN. Flora Patagónica, Tomo VIII, parte VI: Dicotiledóneas Gamopétalas
different land-uses, aquatic macroinvertebrates showed the most (Ericaceae a Calyceraceae). Buenos Aires: INTA; 1999.
obvious and significant responses. This was also the only community Czech B, Krausman P, Devers P. Economic associations among causes of species
that allowed detection of differences between managed native and enlargement in the United States. Bioscience 2000;50:593–601.
Danger AR, Robson BJ. The effects of land-use on leaf-litter processing by macro-
non-managed native forest sites. On the other hand, fish and birds
invertebrates in an Australian temperate coastal stream. Aquat Sci 2004;66:1–9.
having a greater ability of dispersion and capacity to move quickly del Valle HF, Elissalde NO, Gagliardini DA, Milovich J. Status of desertification in the
from disturbances would reflect changes at a higher scale. Patagonian region: assessment and mapping from satellite imagery. Arid Soil Res
Pine plantations affected in different ways the communities Rehabil 1998;12:1-27.
Diamond JM, Dressier DW, Serveiss VB. Assessing relationships between human land-
considered. Exotic plantations produced an intermediate impact on uses and the decline of native mussels, fish, and macroinvertebrates in the Clinch
benthic communities whereas this use resulted in impoverished fish and Powell river watershed, USA. Environ Toxicol Chem 2002;21:1147–55.
and bird communities both in diversity and abundance. The analysis Di Prinzio CY, Casaux RJ, Miserendino ML. Effects of land-use on fish assemblages in
Patagonian low order streams. Ann Limnol 2009;45:267–77.
of different variables (environmental, biotic) with the percentage of Duncan P, Hewison AJM, Houte S, Rosoux R, Tournebize T, Dubs F, et al. Long term
exotic cover did not show any consistent pattern. changes in agricultural practices and wildfowling in an internationally important
Changes in land-use practices are emerging issues concerning the wetland, and their effects on the guild of wintering ducks. J Appl Ecol 1999;36:
11–23.
integrity and conservation of water resources in Patagonia, and this Fernández HR, Domínguez E. Guía para la determinación de los Artrópodos bentónicos
paper proves that stream ecosystems are particularly vulnerable; sudamericanos. Tucumán, Argentina: EUDET; 2001.
moreover it expands understanding to a poorly researched part of the Goldstein RM, Carlisle DM, Meador MR, Short TM. Can basin land-use effects on physical
characteristics of streams be determined at broad geographic scales? Environ
world with distinctive assemblages of species. Monit Assess 2007;130:495–510.
Gordon ND, McMahon TA, Finlayson BL. Stream hydrology, an introduction for
ecologists. New York: Wiley and Sons; 1994.
Acknowledgements Gotelli NJ, Ellison AM. A primer of ecological statistics. Sunderland, Massachusetts, U.S.
A: Sinauer Associates Inc. Publishers; 2004.
Hall MJ, Closs P, Riley RH. Relationships between land-use and stream invertebrate
This paper was supported by the CONICET (PIP 5733), the PADI community structure in a South Island, New Zealand, coastal stream catchment. N Z
Foundation (Grant 19/2006), and Project A.W.A.R.E. Foundation J Mar Fresh 2001;35:591–603.
(P-001036). We also thank Estancia Tecka for the loan of the electro Hauer FR, Lamberti GA. Methods in Stream Ecology. San Diego, California: Academic
press; 1996.
fishing equipment, and the staffs of estancias Los Ñires and Río Frío, Henderson IG, Cooper J, Fuller RJ, Vickery J. The relative abundance of birds on set-aside
and Mr. Cirilo Ñancual who facilitated the access to some of the and neighbouring fields in summer. J Appl Ecol 2000;37:335–47.
sampling sites. We greatly appreciate the reviews of Dr. John Jones EBD, Helfman GS, Harper JO, Bolstad PV. Effects of riparian forest removal on fish
assemblages in southern Appalachian streams. Conserv Biol 1999;13:1454–65.
Richardson for his thoughtful and constructive suggestions. Thanks Kedzierski WM, Smock LA. Effects of logging on macroinvertebrate production in a
to anonymous reviewers for valuable comments that greatly sand-bottomed, low gradient stream. Freshwater Biol 2001;46:821–33.
improved the manuscript. This is Scientific Contribution No. 52 from Kutschker AM, Brand C, Miserendino ML. Un índice de calidad de ribera para ambientes
lóticos de cordillera patagónica sometidos a diferentes usos de la tierra. Ecologia
LIESA.
Aust 2009;19:19–34.
Lara A, Little C, Urrutia R, McPhee J, Alvarez-Garretón C, Oyarzún C, et al. Assessment of
ecosystem services as an opportunity for the conservation and management of
References native forests in Chile. For Ecol Manage 2009;258:415–24.
Lantschner V, Rusch V. Impacto de diferentes disturbios antrópicos sobre las
Allan JD. Landscapes and riverscapes: the influence of land-use on stream ecosystems. comunidades de aves de bosques y matorrales de Nothofagus antarctica en el NO
Annu Rev Ecol Evol Syst 2004;35:257–84. Patagónico. Ecología Aust 2007;17:99-112.
APHA. Standard Methods for the Examination of Water and Wastewater. Hanover, López HL, Miquelarena AM, Menni RC. Lista comentada de los peces continentales de la
Maryland, USA: American Public Health Association; 1994. Argentina. La Plata. División Zoología de Vertebrados. Museo de La Plata. 2003.
Archangelsky M, Manzo V. Descripción de las larvas maduras de los géneros Stethelmis MacArthur R. Environmental factors affecting species diversity. Am Nat 1964;98:387–97.
Hinton y Luchoelmis Spangler and Staines (Insecta: Coleoptera, Elmidae). Rev Mus McCord S, Grippo RS, Eagle DM. Effects of silviculture using best management practices
Argent Cienc Nat 2007;9:79–87. on stream macroinvertebrate communities in three ecoregions of Arkansas, USA.
Armitage PD, Moss D, Wright JF, Furse MT. The performance of a new biological water Water Air Soil Pollut 2007;184:299–311.
quality score system based on macroinvertebrates over a range of unpolluted McKinney ML. Urbanization, biodiversity, and conservation. Bioscience 2002;52:883–90.
running water sites. Water Res 1983;17:333–47. Maitland PS, Morgan NC. Conservation Management of Freshwater Habitats. London:
Bahar MM, Ohmori H, Yamamuro M. Relationship between river water quality and Chapman and Hall; 1997.
land-use in a small river basin running through the urbanizing area of Central Manel S, Buckton ST, Ormerod SJ. Testing large-scale hypotheses using surveys: the
Japan. Limnology 2008;9:19–26. effects of land-use on the habitat, invertebrates and birds of Himalayan rivers. J
Baillie BB, Collier KJ, Nagels J. Effects of forest harvesting and woody debris removal on Appl Ecol 2000;37:756–70.
two Northland streams, New Zealand. NZ J Mar Fresh 2005;39:1-15. Manzo MV. Key to the South American genera of Elmidae (Insecta: Coleoptera) with
Barbour MT, Gerritsen J, Snyder BD, Stribling JB. Rapid Bioassessment Protocols for Use distributional data. Stud Neotrop Fauna Environ 2005;40:201–8.
in Streams and Wadeable Rivers: Periphyton, Benthic Macroinvertebrates and Fish. Meador MR, Goldstein RM. Assessing water quality at large geographic scales: relations
second edition. Washinton. D.C.: U.S. Environmental Protection Agency; 1999. 841- among land-use, water physicochemistry, riparian condition, and fish community
b-99-002. structure. J Environ Manage 2003;31:504–17.
Boët P, Duvoux B, Allardi J, Belliard J. Incidence des orages estivaux sur les peuplements Merritt RW, Cummins KW. An Introduction to the Aquatic Insects of North America.
piscicoles de la Seine á l´aval de l´agglomération parisiense (bief Andrésy- Dubuque, Iowa: Kendall-Hunt; 1996.
Méricourt). Houille Blanche 1999;1–2:141–7. Miserendino ML. Effects of landscape and desertification on the macroinvertebrate
Bojsen BH, Jacobsen D. Effects of deforestation on macroinvertebrate diversity and assemblages of rivers in Andean Patagonia. Arch Hydrobiol 2004;159:185–209.
assemblage structure in Ecuadorian Amazon streams. Arch Hydrobiol 2003;158: Miserendino ML, Brand C, Di Prinzio C. Assessing urban impacts on water quality,
317–42. benthic communities and fish in streams of the Andes Mountains, Patagonia
Buria L, Albariño R, Díaz Villanueva V, Modenutti B, Balseiro E. Impact of exotic rainbow (Argentina). Water Air Soil Pollut 2008;194:91-110.
trout on the benthic macroinvertebrate community from Andean–Patagonian Miserendino ML, Masi CI. The effects of land-use on environmental features and
headwater streams. Arch Hydrobiol 2007;168:145–54. functional organization of macroinvertebrate communities in Patagonian low order
Casaux R, Barrera-Oro E. Effect of a shore-based sampling program on Notothenia streams. Ecolog Ind 2010;10:311–9.
coriiceps populations. Antarct Sci 2002;14:221–4. Miserendino ML, Pizzolón LA. Rapid assessment of river water quality using
Castela J, Ferreira V, Graça M. Evaluation of stream ecological integrity using litter macroinvertebrates: a family level biotic index for the Patagonic Andean zone.
decomposition and benthic invertebrates. Environ Pollut 2008;153:440–9. Acta Limnol Bras 1999;11:137–48.
624 M.L. Miserendino et al. / Science of the Total Environment 409 (2011) 612–624

Miserendino ML, Pizzolón LA. Macroinvertebrates of a fluvial system in Patagonia: Saunders D, Hobbs R, Margules C. Biological consequences of ecosystem fragmentation:
altitudinal zonation and functional structure. Arch Hydrobiol 2000;150:55–83. a review. Conserv Biol 1991;5:18–32.
Miserendino ML, Pizzolón LA. Interactive effects of basin features and land-use change Saunders DL, Meeuwig JJ, Vincent CJ. Freshwater protected areas: strategies for
on macroinvertebrate communities of headwater streams in the Patagonian Andes. conservation. Conserv Biol 2002;16:30–41.
River Res Appl 2004;20:967–83. Scrimgeour GJ, Kendall S. Effects of livestock grazing on benthic invertebrates from a
Molineri C. Impact of rainbow trout on aquatic invertebrate communities in subtropical native grassland ecosystem. Freshwater Biol 2003;48:347–62.
mountain streams of northwest Argentina. Ecol Aust 2008;18:101–17. Sokal RR, Rohlf FJ. Biometry. 3rd edition. New York: W.H. Freeman and Company; 1995.
Morley SA, Karr JR. Assessing and restoring the health of urban streams in the Puget Soto D, Arismendi I, Gonzalez J, Sanzana J, Jara F, Jara C, et al. Southern Chile, trout and
Sound Basin. Conserv Biol 2001;16:1498–509. salmon country: invasion patterns and threats for native species. Rev Chil Hist Nat
Morrone JJ. Biogeographic areas and transition zones of Latin America and the 2006;79:97-117.
Caribbean Islands based on Panbiogeographic and Cladistic analyses of the Strayer DL, Beighley RE, Thompson LC, Brooks S, Nilsson Ch, Pinay G, et al. Effects of land
entomofauna. Annu Rev Entomol 2006;51:467–94. cover on stream ecosystems: roles of empirical models and scaling issues.
Munné A, Solá C, Prat N. QBR Un índice rápido para la evaluación de la calidad de los Ecosystems 2003;6:407–23.
ecosistemas de ribera. Tecnol Agua 1998;175:20–37. Subramanian KA, Sivaramaskrishnan KG, Gadgil M. Impact of riparian land-use on
Narosky T, Yzurieta D. Guía para la identificación de las aves de Argentina y Uruguay. stream insects of Kudremukh National Park, Karnataka State, India. J Insect Sci
Buenos Aires: Asociación Ornitológica del Plata; 1987. 2005;49:1-10.
Nolan PA, Guthrie N. River rehabilitation in an urban environment: examples from the Tell G, Izaguirre I, Quintana R. Flora y Fauna Patagónicas. Bariloche, Argentina:
Mersey basin, North West England. Aquat Conserv 1998;8:685–700. Ediciones Caleuche; 1997.
Paul MJ, Meyer JL. Streams in the urban landscape. Annu Rev Ecol Evol 2001;32:333–65. ter Braak CJF, Smilauer P. CANOCO Reference Manual and User's guide to Canoco for
Pascual M, Macchi P, Urbanski J, Marcos F, Riva Rossi C, Novara M, et al. Evaluating Windows: Software for Canonical Community Ordination (version 4). Microcom-
potential effects of exotic freshwater fish from incomplete species presence puter power; 1998. (Ithaca, NY, USA).
absence data. Biol Invasions 2002;4:101–13. ter Braak CJF, Smilauer P. CANOCO for Windows (version 4.02) a FORTRAN program for
Paruelo JM, Jobbagy EG, Sala OE. Biozones of Patagonia (Argentina). Ecol Aust 1998;8: canonical community ordination. Centre for biometry, Wageningen, The Netherlands.
170–8. 1999.
Pearson O, Ralph C. The diversity and abundance of vertebrates along an altitudinal Thompson RM, Townsend CR. Land-use influences on New Zealand stream commu-
gradient in Peru. Mem Mus Hist Nat “Javier Prado”, Nº 18. ; 1978. Lima, Perú. nities effects on species composition, functional organization, and food-web
Quinn JM, Cooper AB, Davies-Colliey RJ, Rutherford JC, Williamson RB. Land-use effects structure. N Z J Mar Fresh 2004;38:595–608.
on habitat, water quality, periphyton, and benthic invertebrates in Waikato, New Vuilleumier F. Bird species diversity in Patagonia (temperate South America). Am Nat
Zealand, hill-country streams. N Z J Mar Fresh 1997;31:579–97. 1972;106:266–71.
Ralph C. Habitat association patterns of forest and steppe birds of Northern Patagonia, Walling DE, Fang D. Recent trends in the suspended sediment loads of the world's
Argentina. Condor 1985;87:471–83. rivers. Global Planet Change 2003;39:111–26.
Resh VH, Norris R, Barbour M. Design and implementation of rapid assessment Walters DM, Leigh DS, Bearden AB. Urbanization, sedimentation, and the homogenization
approaches for water resource monitoring using benthic macroinvertebrates. Aust J of fish assemblages in the Etowah River Basin, USA. Hydrobiologia 2003;494:5-10.
Ecol 1995;20:108–21. Wilcox B, Murphy D. Conservation strategy: the effects of fragmentation on extinction.
Richardson JS. Aquatic arthropods and forestry: effects of large-scale land-use on Am Nat 1985;125:879–87.
aquatic systems in Nearctic temperate regions. Can Entomol 2008;140:495–509. Willson MF, De Santo T, Sabag C, Armesto JJ. Avian communities of fragmented South
Roy AH, Rosemond AD, Paul MJ, Length DS, Wallace B. Stream macroinvertebrate temperate rainforest in Chile. Conserv Biol 1994;8:508–20.
response to catchment urbanization (Georgia, USA). Freshwater Biol 2003;48: Wohl E. Mountain Rivers, Water Resources Monograph 14. Washington, D.C.: American
329–46. Geophysical Union; 2000.