Next Article in Journal
Seasonal Variations of Dissolved Organic Matter by Fluorescent Analysis in a Typical River Catchment in Northern China
Next Article in Special Issue
A Planning Tool for Optimizing Investment to Reduce Drinking Water Risk to Multiple Water Treatment Plants in Open Catchments
Previous Article in Journal
Reconstructing Spatiotemporal Dynamics in Hydrological State Along Intermittent Rivers
Previous Article in Special Issue
Impact of Land Use Change on Non-Point Source Pollution in a Semi-Arid Catchment under Rapid Urbanisation in Bolivia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Intensive Agriculture on Benthic Macroinvertebrate Assemblages and Water Quality in the Aconcagua River Basin (Central Chile)

1
Instituto de Ciencias Marinas y Limnológicas, Universidad Austral de Chile, Valdivia 5090000, Chile
2
Departamento de Sistemas Acuáticos, Facultad de Ciencias Ambientales, Universidad de Concepción, y Centro de Ciencias Ambientales (EULA), Universidad de Concepcion, Concepcion 4070386, Chile
3
Departamento de Ecología, Facultad de Ciencias, Universidad Católica de la Santísima Concepción, Concepción 4090541, Chile
4
Departamento de Física, Facultad de Ciencias, Universidad del Bío-Bío, Concepción 4051381, Chile
5
Centro FONDAP de Investigación en Dinámica de Ecosistemas Marinos de Altas Latitudes (IDEAL), Valdivia 5090000, Chile
*
Author to whom correspondence should be addressed.
Water 2021, 13(4), 492; https://doi.org/10.3390/w13040492
Submission received: 30 December 2020 / Revised: 4 February 2021 / Accepted: 5 February 2021 / Published: 14 February 2021
(This article belongs to the Special Issue Impact of Land-Use Changes on Surface Hydrology and Water Quality)

Abstract

:
This study assessed natural variation in the macroinvertebrate assemblages (MIB) and water quality in one of the main basins with the largest agricultural activities in Chile (Aconcagua River Basin). We sampled throughout the annual cycle; nine sampling sites were established along the basin, classifying according to agricultural area coverage as least-disturbed, intermediate, and most-disturbed. We collected 56 macroinvertebrate taxa throughout the entire study area. Multivariate analysis shows significant differences among the three disturbance categories in different seasons, both water quality variables and the MIB structure. Distance-based linear model (DistLM) analysis for all seasons explained more than 95.9% of the macroinvertebrate assemblages, being significantly explained by chemical oxygen demand, pH, total coliforms, nitrites, elevation, and water temperature. ANOVA test revealed significant differences in the proportion of noninsect individuals, macroinvertebrates density, and the number of taxa among the three disturbance categories (p < 0.05). In general, water temperature, conductivity, chemical oxygen demand, ammonium, nitrites, and nitrates increased their values downstream in the basin. Our results indicate that the elevation gradient and increment in agricultural land use in the basin had a strong influence on water quality and MIB. A better understanding of these ecosystems could help conservation and integrated watershed management.

1. Introduction

Limnologists have studied natural changes in the aquatic community composition through altitude gradient for many decades [1]. In this sense, the River Continuum Concept-RCC [2] proposes a taxonomic composition change in macroinvertebrates along the river continuum, from headwaters to mouth, increasing, for example, filter-feeding species downstream due to an increase in dissolved organic matter. However, natural changes in macroinvertebrates can be altered by anthropogenic activities, such as land-use changes [3]. Worldwide intensification of agriculture, one of the main stressors of aquatic ecosystems worldwide, has adverse effects on water quality and therefore changing the structure of streams communities, impacting the biodiversity and functions of freshwater ecosystems negatively [4]. Lawrence et al. [5] found that an increasing agricultural area in the watershed of California Mediterranean streams was negatively related to the diversity of streams communities and water quality. Similarly, Fierro et al. [6] recorded a negative relation among agriculture land-use percentage and habitat water quality and streams communities in Mediterranean Chilean streams. Thus, it is essential to evaluate if models predicted by RCC are fulfilled in streams affected by anthropogenic land use, to use this information in conservation policies that differentiate both natural and anthropogenic changes in stream communities.
Macroinvertebrates have been used for decades as bioindicators because they respond to environmental changes as water quality and land use, increasing or decreasing their abundances, occasionally even disappearing [7,8]. However, because macroinvertebrate assemblages also change temporally by intrinsic factors (e.g., life-history cycles), studies comprising seasonal changes are essential to discriminate natural and anthropogenic effects on aquatic communities [9]. In southern South America, maximum abundances and richness have been reported in summer months (January, February, and March), related to seasonal fluctuations in hydrology, higher water temperatures, lower flows, and most stable conditions. In contrast, lower abundances and richness of aquatic insects occur in winter [10,11].
Agricultural land use often degrades riparian habitat, decreasing Canopy River, and commonly altering the water quality, increasing water temperature, nutrients, and fine sediments in streams [12,13]. Due to the unidirectional nature of the streams, nutrients and pollutants (e.g., fertilizers or agrochemicals) can accumulate downstream with higher concentrations in lower areas of the basin [14]. Moreover, reductions in streamflow associated with water diversions to irrigation canals are common in agricultural Mediterranean regions [15,16]. These environmental features associated with agricultural activities can negatively influence aquatic macroinvertebrates. Some agricultural land-use consequences are the decrease of sensible taxa (such as EPT), increased relative abundances of tolerant taxa such as Chironomids or oligochaetes, and noninsects invertebrates, such as snails [17,18].
Chilean Andean streams draining to the Pacific Ocean are relatively shorts (<400 km), with marked land use in the Mediterranean ecoregions. Uppers zones in these basins are dominated by native vegetation, such as scrublands, while central valleys are generally dominated by exotic forest plantations and agricultural land use. The lowlands of the basins are dominated almost exclusively by agricultural lands [6]. The Aconcagua River Basin is located in central Chile, comprises an elevation up to 6100 m above sea level, and it is one of the main basins that sustain the economy based on agriculture in Chile (supports 12% of Chile national agriculture), livestock and forestry production [19]. Thus, the Aconcagua basin provides an excellent study model to examine longitudinal patterns of macroinvertebrate assemblages associated with altitudinal and land-use agriculture gradients—agricultural areas increase toward lowlands of the basin.
This study analyzes the seasonal structure of macroinvertebrate assemblages and water quality along the Aconcagua River Basin. We predict that landscape variables, such as altitude or catchment area, and environmental variables derived from agricultural activities, such as increased nutrients, are directly related to macroinvertebrate assemblages composition.

2. Material and Methods

2.1. Study Area

The Aconcagua River basin (32°54′ S, 71°30′ W) has a surface area of 7200 km 2 (Figure 1). The climate is typically Mediterranean, characterized by warm, dry summers (between October and March) and wet, cool winters (mainly between May and August) with intense and irregular rainfall [20]. Temperatures and annual precipitation vary along the basin. Lowland average annual temperature is 14.5 °C, and precipitation reaches 395 mm/year; about 15.2 °C and 261 mm/year in the medium sector 14.1 °C and 467 mm/year in the highlands. The Aconcagua is 5th order river according to Strahler and is 142 km long, forming from Juncal and Blanco River’s confluence. The basin has a mixed hydrological regime with rain and snow contributions, an average annual flow of 33.1 m 3 s 1 . The maximum flow peaks in November–January, whereas the minimum flow occurs between March and September [21].
The Aconcagua River Basin is composed mainly of igneous rocks interbedded with marine and continental sediments, whose ages fluctuate between the Upper Triassic and the Upper Miocene [22]. The basin encompasses the three main geomorphological provinces of Chile: Andes mountains, central valley, and coastal mountains and can be divided into three parts according to its geomorphology and distribution of geological units [23]. In the upper section of the basin, high slopes do not allow a large amount of water to infiltrate the underground water tables. In the medium section, the terrain is smooth, and the slopes are covered with a detrital mantle accompanied by basement rocks in some places. Given the little rainfall, erosive effects are limited compared to the upper part of the basin. In this part of the basin, the detrital formations also influence the soils’ permeability capable of containing a large amount of water, constituting a large groundwater reservoir. In the lower section of the basin, the Aconcagua River narrows to only one watercourse.
The Aconcagua valley is characterized by intensive irrigated agriculture. From the 145 inflows in the river, 90% are diverted for irrigation channels [21]. The basin is divided into four sections used to draw water to agricultural practices. The sections are: (i) from first agriculture terrains in the upper part of the basin (975 m a.s.l.) to San Felipe city (648 m a.s.l.), (ii) from San Felipe city to Puntilla Romeral (211 m a.s.l.), (iii) from Puntilla Romeral to Quillota city (127 m a.s.l.), and (iv) from Quillota city to Aconcagua River mouth in Concón bay. Aconcagua River Basin supports about 500,000 inhabitants, mainly from the cities San Felipe (71,000 inhabitants), Los Andes (63,000 inhabitants), and La Calera (50,000 inhabitants), all located in the middle and lower part of the basin. High nitrate concentrations are related to agriculture activities [24], whereas high concentrations of copper and molybdenum have been reported in the upper basin [25].

2.2. Land Use and Sampling Sites

The Aconcagua Valley is one of the most important agricultural regions of Chile [26]. We sampled streams highly impacted by agriculture on 9 activities, i.e., streams drain land that has been transformed mainly for agriculture use, close to the streams along the basin. A total of 12 types of agriculture can be recognized in the basin (Figure 2, Appendix A), being the majority of the land used for vineyards, fruit trees, horticulture, and avocado plantations. The highest elevations are dominated by the typical evergreen sclerophyll forests and hillsides by dry xerophytic thorn scrub and forest shrublands [27].
We selected a total of nine sampling sites in the Aconcagua River Basin (Table 1), according to a natural elevation gradient occurring in the basin (Figure 1) and a gradient of agricultural land use (Figure 2). Sample stations located at a higher elevation (>1096 m a.s.l.) and with agricultural land use < 0.2% (respect to the entire basin) corresponded to stations JU and JU10 (Juncal river), BL20 (Blanco river), and AC10 (Aconcagua river) and were categorized as least-disturbed (LD). Stations AC20 (Aconcagua river) and PU10 (Putaendo river) located at the middle part of the basin (650–604 m a.s.l.), categorized as intermediate-disturbed (ID) to have an agricultural land use between 6.8 and 7.5%. Finally, stations PO10 (Pocuro river), AC30, and AC40 (Aconcagua river), located at a lower elevation (<600 m a.s.l.), had the higher agricultural land use >10.8% being categorized as the most-disturbed (MD). BL20 sampled site was located downstream of the Blanco River watershed, where a mining plant (copper exploitation) is located in the upper part. AC20 sampled site was located approximately 1 km downstream of the sewage outflow of the San Felipe city.
Land use cover (%) for each site was estimated by screening digitalized satellite images (1:30,000 scale, 2013 Landsat 8 sensor) freely available from the Satellite Information System, IDE Agriculture (http//www.ide.cl (accessed on 1 February 2021)). Land-use types were determined using ArcGis10 [28]. We obtained the subcatchment area above the sampling site and the coverage (%) of land uses types.

2.3. Field Sampling and Data Collection

Macroinvertebrates and water samples were sampled during winter 2014 (24–29 August), spring 2014 (24–29 November), summer 2015 (9–14 March), and fall 2015 (15–20 June). Macroinvertebrates at each site were sampled from the riffles, the most common habitat. For each season, macroinvertebrates samples consisting of 9 randomly located Surber samples (250 µm mesh netting; sampling area = 0.09 m 2 ) were collected. The samples were fixed in situ with 70% ethanol and then taken to the Biodiversity Laboratory of the Universidad de Concepción in Concepción, Chile.
The nine macroinvertebrate samples collected at each site were analyzed separately. A total of 315 replicates were analyzed. Station PU10 was not sampled in June because the stream was dried up. Samples were washed, sorted, and all organisms were separated from the detritus and preserved in 90% ethanol. Then organisms (larvae/pupae and adults in the case of riffle beetles) were identified and counted under stereo microscope (Zeiss, model Stemi Dv4, 32X magnification) to the lowest practical taxonomic resolution using dichotomous keys and taxonomic descriptions [29]. In Chile, there is insufficient information to identify all organisms to species. Most insects of the order Ephemeroptera were identified to species, whereas Trichoptera, Coleoptera, Odonata, Hemiptera, Lepidoptera, and Diptera were identified mostly to genus or family. Except for some rare taxa, most noninsects were identified to genus. The first author performed all identifications to maintain consistency among sample sets. For each site, macroinvertebrate count data were converted to densities by square meter.
At each site, 28 environmental and habitat variables were assessed (Table 1 and Table 2). Water temperature (°C), electrical conductivity (µS/cm), pH, and turbidity were measured in situ with a Hanna multiparameter instrument model HI 9828. Water quality samples were collected in duplicate from the center of the active channel below the water surface, deposited in prewashed glass bottles (1 L), and transported in coolers to 4 °C to the laboratory for chemical analyses. In the laboratory, the following water parameters were measured: DO-Dissolved Oxygen (mg/L), BDO 5 -Biological Oxygen Demand(mg/L), COD-Chemical Oxygen Demand (mg/L), Organic suspended solids (mg/L), Inorganic suspended solids (mg/L), total phosphorous (µg/L), ammonium (µg/L), nitrite (µg/L), nitrate (mg/L), Al, As, Cu, Fe, Mn, Mb, and Zn dissolved (µg/L), sulfate (mg/L), total hardness (mg/L), alkalinity (mg/L), chlorides (mg/L), fecal coliforms, and total coliforms (Nmp/100 mL). All analyses were carried out using standard methods for water and wastewater [30].

2.4. Data Analyses

The following macroinvertebrate assemblages metrics descriptors were calculated using the level taxa (such as [6,16,31]), averaging all seasons for each sampling site: insects in the disturbance sensitive orders Ephemeroptera, Plecoptera, and Trichoptera (% EPT individuals and number of EPT taxa) and % noninsects individuals, abundance descriptors (i.e., total density macroinvertebrates), species richness (i.e., total number taxa), and Shannon–Weaver diversity index (H’). These metric values for low, medium, and highly disturbed agricultural sites were compared using an analysis of variance (ANOVA), followed by a Tukey test to identify different groups. The differences were considered to be significantly different at p < 0.05.
The environmental variables were elevation and catchment area (described in Table 1) and the 26 variables described in Table 2. We used canonical analysis of principal coordinates— CAP [32]—to test for differences in macroinvertebrates taxonomic composition and environmental variables among the land-use types. CAP is a variant of principal coordinates analysis (PCOA), which can find axes through the multivariate cloud of points to discriminate among a priori groups [33]. The analysis can be based on any resemblance matrix. We used Bray–Curtis dissimilarity for biological data and Euclidean distances for environmental data. Biological data were square-root transformed before performing the CAP. The diagnostic test in CAP to allocate the sites to correct classification groups and test for agriculture intensification differences was run in macroinvertebrate composition using random permutations. We also tested the null hypothesis of no differences among group centroids (agriculture intensification type) using 999 permutations [34].
Relationships between macroinvertebrate community composition and environmental variables were examined using distance-based redundancy analysis (dbRDA) and distance-based linear models (DistLM) [33] ( α = 0.05; 999 permutations). These analyses were based on Bray–Curtis dissimilarity matrices of macroinvertebrate data and the whole set of environmental variables. Before dbRDA and DistLM analysis, the environmental features (except pH) were square-root transformed, and highly correlated variables (Pearson’s r > 0.9) were removed. In winter, conductivity, turbidity, Zn, and alkalinity were removed; spring: conductivity, turbidity, SSO, SSI, P, Mn, Mb, sulfate, fecal coliforms, and elevation; summer: turbidity, BOD 5 , nitrite, sulfates, total hardness; fall: turbidity, BOD 5 , SSO, Zn, Mn, and Clorures. The dbRDA and DistLM were run for the macroinvertebrates and environmental data among the 4 seasons. We reported the adjusted R² and cumulative proportion for environmental variables from the constrained ordinations. We ran ANOVAs using the R software [35], whereas CAPs, DistLMs, and dbRDAs were conducted in PERMANOVA+ for PRIMER [33].

3. Results

3.1. Environmental Variables and Water Quality

The elevation of sampling sites ranged from 12 m a.s.l. to 1657 m a.s.l., and Strahler order ranged from 3 to 5, reaching this last value in the central valley. The main substrate types were sand, gravel, cobble, and boulders, with the mostly size substrates in the higher altitude sampling sites (BL20, JU, JU10) (Table 1). Agricultural land use characterized the intermedium and most-disturbed sites, ranging from 6.8% to 19.1% of basin land, recorded agricultural sites at altitudes below 604 m a.s.l. (Table 1).
Generally, water temperature, conductivities, COD, ammonium, nitrites, and nitrates were lower in LD stations (stations higher altitude), increasing downstream. Conversely, metal concentrations were higher in LD stations, specifically at BL20 (mouth of Blanco River), decreasing downstream (Al, Cu, Mn, Mb). Sulfates, total hardness, alkalinity, chlorides, fecal coliforms, and total coliforms followed the same pattern: increasing their concentrations downstream, reaching maximum values near the river mouth, and minimum values around headwaters. The pH, turbidity, DO, BOD 5 , organic suspended solids, and inorganic suspended solids varied through the basin without a clear altitudinal trend (Table 2, Appendix B).

3.2. Macroinvertebrate Assemblages

We collected 56 macroinvertebrate taxa in the study area across all the year from riffle habitats, of which 36 were insects (mostly Diptera taxa) and 20 noninsects. Among insects, Diptera was the most diverse order (11 taxa), followed by Ephemeroptera (7), Trichoptera (6), and Coleoptera (5). Numerically Chironomids were the taxa most abundant and found at 100% of sites throughout the year. Noninsects were represented by crustaceans, spiders, snails, worms, free-living flatworms, leeches, and nematodes, of which snails Physa chilensis, Litorridina sp., leeches Glossiphonidae and worms Tubifex sp. were found in greater abundance in intermedium and most-disturbed sites.
Seasonally, 39 taxa were reported in winter with 23,782 individuals, 36 in spring with 29,885 individuals, 34 in summer with 17,512 individuals, and 38 in fall with 15,644 individuals. Most-disturbed agricultural sites had the highest richness of macroinvertebrates, both aquatic insects and noninsects yearlong. Surprisingly, juvenile life stages of Plecoptera were not recorded in the study area, whereas other disturbance-sensitive orders, Ephemeroptera and Trichoptera were uncommon in least-disturbed sites. In contrast, the Trichopteran Metrichia sp. and Smicridea sp., the Ephemeropteran Meridialaris laminata, Meridialaris chilooense, Massartellopsis irarrazavali, and the Dipteran Blephabericeridae, were taxa founded only in least-disturbed sites.
A clear pattern was observed in the % noninsect individuals, densities, and the number of taxa, which increased significantly from least-disturbed agricultural sites to most-disturbed sites (Figure 3). Statistical analysis (ANOVA, Tukey’s test, p < 0.05) showed statistically significant differences in % of noninsect individuals, densities, and the number of taxa among least and most agricultural disturbed sites. Despite %EPT individuals and Shannon–Weaver diversity shown their highest median values in most perturbed sites, these were not significant (p < 0.05). A similar situation occurred with EPT richness, which did not vary between levels of agricultural sites (Figure 3).

3.3. Multivariate Analysis: Environmental Variables and Macroinvertebrate Assemblages

The CAPs based on environmental features showed that sampling sites were not statistically different among the three agriculture disturbed types across the year (Table 3). Despite this, the percentage of correct classification of sites to their parental groups reached 100% in least and intermediate agricultural streams yearlong. Sites classified as most-disturbed reached 66.7% of correct classification because site PO10 should have been classified as least-disturbed.
The CAPs based on macroinvertebrate assemblages showed that sampling sites were statistically different among the three agriculture intervention types, both in spring and summer (Table 3). In winter, the percentage of correct site classification to their parental groups reached 50%, and none of them was correctly classified as an intermediate-disturbed site. In the fall, the percentage of correct classification of sites to their parental groups reached <75% (Table 3).
Environment-assemblage relationships varied across the seasons, as shown by differences of both the significant predictor variables and the predicted power (adjusted R²), recorded both natural (e.g., temperature, elevation, catchment area) and anthropogenic variables (e.g., total coliforms, nitrates, nitrites) (Table 4). DistLM explained 96.9% of the variation in the macroinvertebrate assemblages in winter, and contributions from seven environmental variables achieved the best solution, in which COD and temperature were significant (Table 4). DistLM for spring explained 97.7% with contribution from seven variables in which temperature and pH were significant. DistLM for summer explained 97% with contribution from seven variables to the best solution in which elevation and total coliforms were significant. Finally, DistLM for fall explained 95.9%, and six environmental variables achieved the best solution. The elevation, total coliforms, and nitrites were significant (Table 4).
The dbRDA model (Forward) indicated that the least and most-disturbed sites were strongly separated along axis 1 and 2, and intermediate-disturbed sites (depending on the seasons) were closer to the least or most-disturbed sites. For winter (Figure 4), vectors for Mn concentration, COD, and fecal coliforms were closely aligned parallel to axis 1 and intermediate and most-disturbed sites were scattered in this direction. The vector for nitrates was oriented along axis 2, also close to intermediate and most-disturbed sites. On the dbRDA plot for spring (Figure 4), vectors for catchment area, nitrates, and total coliforms were closely aligned parallel with axis 1, and vectors for pH, temperature, Fe, and Nitrites were in orientation to axis 2. Intermediate and most-disturbed sites were close to these vectors. On the dbRDA plot for summer (Figure 4), vectors for catchment area, fecal coliforms, elevation, and SSO were closely aligned parallel with axis 1, and total coliforms, COD, and Zn were in orientation to axis 2. Intermediate disturbance sites were partially intermingled with least and most-disturbed sites and scattered almost equally along axis 1 and 2. The dbRDA plot for fall (Figure 4) vector for pH was closely aligned parallel with axis 1, whereas elevation, fecal coliforms, total coliforms, ammonium, and nitrites were principally aligned with axis 2. Vectors for fecal and total coliforms were linked to most-disturbed sites.

4. Discussion

Land-use changes, such as those associated with agricultural activities, degrade the habitat and the water quality in various rivers around the world [36,37]. However, such changes can be masked by the natural changes that occur in the basins, such as the altitudinal gradients that occur along a river, influencing the water’s quality and the aquatic communities [38,39]. We evaluated these effects in the Aconcagua River Basin since the increase in the proportion of agricultural activity is related to an increase in the basin’s size and a decrease in altitude. We found that water quality changes were mainly reflected in increased temperature and nutrients (e.g., nitrates, nitrites) downstream. The increase (decrease) in agricultural activity (altitude) was also associated with the increase in noninsect individuals and the increase in macroinvertebrates’ density and richness.
The increase in agricultural activity impacted the macroinvertebrate assemblages, reflected in the CAP analysis for the spring and summer seasons, and significant differences for three of the six evaluated metrics. These results are consistent with numerous studies related to the impact of agriculture [17], which relate an increase in the proportion of noninsect individuals and an increase in the density of macroinvertebrates. On the other hand, macroinvertebrates richness was higher in most-disturbed sites by agriculture. These results are contrary to those reported in the literature, where lower species richness is expected in places with vast agricultural land [6,40]. However, this increase in richness is due to the decrease in altitude. Studies in rivers covering a wide altitudinal gradient, such as the Colorado River in the USA, have shown that as altitude decreases, there is an increase in the richness of macroinvertebrates [41].
The richness of macroinvertebrate taxa was similar throughout the year. Simultaneously, the density increased notably in spring, resulting from the seasonal changes of life cycles of insects such as Trichoptera Metrichia, Smicridea and the diptera Chironomidae in intermedium and most-disturbed stations. The EPT index also remained stable along the continuous river, so that this index would be a poor indicator of the environmental status of the Aconcagua River Basin. The nonvariation in the index values results from the high abundance in all sampling stations of some ephemeropterans (Andesiops torrens and A. peruvianus) and trichopterans (Hydroptilidae and Smicridea sp.). Both species of Andesiops are widely distributed in southern South America, from high-altitude rivers up to river mouths at sea level [42]. These results suggest that A. peruvianus and A. torrens seem to be tolerant to the increase in nitrites (ranged from 22.5–852.5 µg/L), nitrates (ranged from 1–17.7 mg/L), total phosphorous (ranged from 10–200 µg/L), water temperature (ranged from 6.3–17 °C), and other environmental variables in the Aconcagua River Basin.
Although the Plecoptera order is characteristic of mountain rivers, especially in Andean rivers, no species were reported in the study area. Similar results found Scheibler et al. [43] in a river at the same latitude as the Aconcagua river but on the other side of the Andes, in Argentina, which only recorded one Plecoptera (Gripopterygidae) taxa. Probably the absence of stoneflies nymphs in the upper part of the Aconcagua River is due to the small entry of woody residues or leaf-litter from the riverine vegetation, the main food of these insects, which is naturally scarce in the upper part of the basin, where the use of bare soil and rocky outcrops (see Figure 2). In the upper part of the basin (>1096 m a.s.l.), where there is a low percentage of agricultural land use, the macroinvertebrate assemblages were represented by chironomids. Furthermore, the presence of insects sensitive to pollution, only present in this sector of the basin, were represented by ephemeropterans Meridialaris laminata, Meridialaris chilooense, Massartellopsis irarrazavali, and the Dipteran Blephabericeridae. These taxa are known to prefer mountain rivers with cold, well-oxygenated waters [10,44]. On the other hand, in the middle sector of the basin, sites moderately disturbed (660–604 m a.s.l.), and in the lower sector with sites most-disturbed by agricultural land use (<600 m a.s.l.), the assemblage was again represented by Chironomids and a high abundance of the insects Smicridea and Hydroptilidae. Both taxa are characteristic for being found in sites with a high organic material content, low DO, high conductivities, and warm temperatures [45,46].
Our results showed statistical differences among the three agricultural disturbed types based on macroinvertebrates assemblages in spring and summer. On the other hand, non statistical differences were founded based on environmental features. These results indicate that ecological indicators are more sensitive that environmental features to detect different degrees of disturbance based on agricultural disturbances, as described in other Chilean streams [47]. Notice that site PO10 was incorrectly classified in its parental group (most-disturbed) based on environmental features, which should have been classified as leas-disturbed according to CAP analysis. Pocuro stream (where PO10 site was located) has the highest proportion of agricultural land use (19.1%), principally grapes and parron-vid (Figure 2, Appendix A). However, it also has the highest proportion of shrublands land use (55.9%). Thus, it seems that this stream’s elevated agricultural activities are not enough to classify this sampling station as a most-disturbed site based on environmental features. From the 28 environmental variables included in the DistLM for each station, six to seven variables contributed to the model’s best solution (Table 4). The significant variables with a high percentage of contribution (>26.8%) throughout the year were the DOC, temperature, pH, altitude, total coliforms, and nitrites. Several studies have shown that high levels of DOC (an indicator of organic pollution) and high pH values can affect the structure of the macroinvertebrate assemblages in rivers associated with the use of agricultural land [48,49]. Thus, our significant variables follow the general models of agricultural basins, having strong effects on the structure of aquatic biodiversity. Altitude is an environmental factor that affects macroinvertebrate assemblages and water quality, for example, by altering temperature and dissolved oxygen [50]. Vannote et al. [2] suggests that the highest species richness occurs in the basins middle reaches. However, we found the highest richness in the lower reaches of the basin. Thus, the expected assemblage of macroinvertebrates by the RCC would be modified by the effects of agriculture in the Aconcagua River Basin.
The Aconcagua River Basin decreased its flow in recent years, being declared to have water scarcity by the Chilean government [51]. The dramatic decrease in river flow is associated with climate change and activities related to agriculture, such as excessive water diversions for irrigation canals in the middle and lower part of the basin, which implies a decrease in the main river flow and consequent increase in its temperature. The water temperature is a critical variable in the macroinvertebrate assemblage structuring, which naturally varies along the continuous river related to the basin’s altitude and size. In our study, the highest altitude station had an average of 6.3 °C throughout the year, whereas the lowest stations had an average temperature range between 14.4 and 17 °C. Rivers with the highest proportion of agricultural soils generally have little riparian vegetation. This lack of vegetation causes rivers to be more exposed to solar radiation, and therefore, an increase in water temperature, nutrients, fine sediments, and pesticides is expected [52,53,54]. In our study, the total and fecal coliform variables, nitrates, nitrites, increased their concentrations in intermedium and most-disturbed areas. For example, the average values of phosphates and ammonium reached 100 µg/L and 1067 µg/L, respectively, closest to the river mouth. Higher values of the total N-content in the water have been recorded in rainy seasons in the study area, associated with agricultural sources [24]. Nutrients play an essential role in structuring aquatic communities [55]. One possible mechanism is that the increase in nutrients recorded in this study in the middle and lower areas of the basin (i.e., phosphorus, ammonium nitrites, nitrates) promotes an increase in the biomass of the periphyton, the main food for browsing macroinvertebrates [6]. Indeed, a large number of gastropod snails were recorded in these areas of the basin. This would also explain the lack of statistical differences in most-disturbed sites based on the Shannon–Wiever diversity index. Heino et al. [56] found that taxonomic richness indices could not differentiate disturbed sites from reference sites based on macroinvertebrate assemblages. The replacement of species from least- to most-disturbed sites was reflected by the increase in the richness and density of noninsect taxa in the middle and lower basin. These results are consistent with studies carried out in the lower area of the Aconcagua River Basin, where sites affected by agriculture dominated flatworms, hirudineans, worms, snails, and chironomids [57].
The high concentration of heavy metals in the upper basin is consistent with the natural sources of metals in the environment (include metallic minerals and parent rocks) and the presence of a mining company of copper and molybdenum in the Blanco River (BL20 station) [58]. The decrease in the concentration of metals downstream may result from their dissolution by contributing other tributaries to the main river and the natural dispersion with an increase in distance from the source [59]. Heavy metals As, Mn, Fe, and Zn were nonsignificant variables but were included in the model throughout the year. These parameters were generally high in the upper basin, and their concentrations decreased downstream. Our results are similar to those reported by Alvial et al. [60] in a nearby basin with heavy metals naturally present and variables that change naturally through a continuum river such as altitude and dissolved oxygen explaining the structure of the macroinvertebrate assemblages.
Given the limited knowledge of altitude and agriculture’s effects on these Mediterranean ecosystems, our findings are relevant for the basin’s management and regulation. Agricultural activities affect water quality parameters and benthic communities. However, their impacts can be underestimated because of the river’s natural changes due to altitude gradients. The Aconcagua river basin is one of the main basins of central Chile affected by climate change, which is reflected in drought conditions. Thus, we suggest that decision-makers regulate agricultural activity in the region, emphasizing the river’s flow and its tributaries to preserve aquatic biodiversity.

5. Conclusions

The altitude gradient and the increase in agricultural activity downstream in the Aconcagua River Basin affected several physicochemical parameters, which varied throughout the year and led to increased nutrient levels in the areas mostly disturbed by agricultural activities. These changes were also reflected in the macroinvertebrate assemblages structure, associated with an increase in noninsect taxa and greater abundance and richness of macroinvertebrates downstream. Considering that a large part of the basins of the Mediterranean area of Chile follows the same land use pattern (the upper basin has little human intervention, while agricultural activities dominate the middle and lower areas), we expect similar results in other rivers. Integrated watershed management must consider the effect of various anthropogenic stressors on water quality and aquatic biota, including the changes that occur naturally by altitude gradients and the seasonal variability of water quality and biota.

Author Contributions

Conceptualization, P.F. and C.V.; methodology and data curation, P.F., C.V., G.S.S., and C.L.; writing—original draft preparation, P.F and C.V.; writing—review and editing, P.F., C.V., G.S.S., and C.L.; funding acquisition, P.F and C.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by FONDECYT grant number 11190631 and DAND Codelco-Andina. CL has been funded by UPWELL Millenium Nucleus (grant ICM NCN 19153) and red GeoLIBERO–CYTED.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sets are contained within the article. Further data and materials requests should be addressed to [email protected].

Acknowledgments

We thank W. San Martin, A. Jara, C. Cerna, J. González, and A. Manosalva for his assistance in sampling. Three anonymous reviewers and an editor provided comments that improved our manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Map of land uses and sampling sites on the Aconcagua River Basin. Direct land use corresponds to land use adjacent to watercourses, principally anthropogenic land use. Indirect land use corresponds to land use in a wider catchment area. Name of sampling stations in Figure 1, black circles, least-disturbed sites; lead squares, intermedium-disturbed sites; white triangles, most-disturbed sites by agricultural activities.
Figure A1. Map of land uses and sampling sites on the Aconcagua River Basin. Direct land use corresponds to land use adjacent to watercourses, principally anthropogenic land use. Indirect land use corresponds to land use in a wider catchment area. Name of sampling stations in Figure 1, black circles, least-disturbed sites; lead squares, intermedium-disturbed sites; white triangles, most-disturbed sites by agricultural activities.
Water 13 00492 g0a1

Appendix B

Table A1. Environmental features at 9 sampling sites on the Aconcagua River basin, Chile. Values represent average of four seasons (August and November 2014, March and June 2015). Minimal and maximum (in brackets) are consigned.
Table A1. Environmental features at 9 sampling sites on the Aconcagua River basin, Chile. Values represent average of four seasons (August and November 2014, March and June 2015). Minimal and maximum (in brackets) are consigned.
JUBL 20JU 10AC 10AC 20PU 10PO 10AC 30AC 40
Temperature (°C)6.3 (3–10.3)8.5 (4.5–14.9)10.2 (6.9–12.6)13.4 (10–16.2)17 (10.9–21.7)15.6 (12.1–17.7)15.2 (7.3–22.8)16.7 (10.2–20.4)14.4 (5.4–18.9)
Electrical conductivity (S/cm)367.8 (204.2–560)720.8 (362–1122)303.2 (235.2–375)326.5 (257.2–432)643 (452–794)653.7 (558–721)493 (379–585)621.3 (594–663)938.5 (813–1089)
pH7.1 (6.3–8)6.9 (6–7.5)7.4 (6.6–8)8 (7.5–8.4)7.9 (7–8.6)7.6 (7.4–7.8)7.5 (6.9–8)8.2 (7.9–8.6)7.6 (7.3–8.3)
Turbidity (NTU)32 (11–55)36.8 (5–100)14.5 (2–47)47.8 (5–170)50.5 (6–111)3 (1–5)85.3 (45–132)5.3 (4–9)16.3 (3–44)
Dissolved oxygen (mg/L)10.4 (9.6–11.6)9.3 (8.4–11.8)9.9 (9.5–10.3)9.9 (8.8–10.8)10.6 (9.6–11.2)9.7 (7.9–10.8)9.8 (8.5–11.8)11.3 (10.4–12)10.2 (6–13.2)
BDO5 (mg/L)1.1 (1–1.3)5.8 (2.6–7.8)1.1 (1–1.4)1 (1–1.1)2.9 (2.2–3.4)1 (1–1.1)2 (1.2–2.6)1.4 (1–1.6)1.5 (1.1–1.9)
COD (mg/L)5.6 (2.5–8.5)15.5 (9.7–24.1)5.5 (2.5–7.6)6.3 (2.5–11.5)14.7 (3.3–23.9)5.2 (2.7–8.5)9.8 (4–12.3)8.6 (5.4–12.3)10.5 (6.8–14.2)
Organic suspended solids (mg/L)1.6 (1.2–2.5)22 (1–3.8)2.5 (1–7)1.8 (1–3.5)7.1 (3.1–11.8)111.8 (6.3–24)1.1 (1–1.3)2.5 (1–4.2)
Inorganic suspended solids (mg/L)18.9 (10.1–35.8)8 (1–15.7)4.7 (1–11.8)16.5 (2.6–57.4)37.3 (3.1–79.2)1 (1–1.1)48.3 (31.3–62.8)2.2 (1.1–3.9)15.2 (2–42.2)
Total phosphorus (g/L)67.5 (20–120)62.5 (10–150)35 (10–70)132.5 (10–420)200 (70–270)10187.5 (120–240)50 (10–90)65 (20–120)
Ammonium (g/L)2032.5 (20–70)25 (20–40)201067.5 (60–1920)20142.5 (20–320)22.5 (20–30)57.5 (30–90)
Nitrite (g/L)22.5 (15–45)133.5 (15–478)24.5 (15–53)26.3 (15–60)852.5 (114–1889)33.3 (15–0)87.3 (36–174)36.5 (15–59)153.5 (24–366)
Nitrate (mg/L)1 (0.9–1.2)4.4 (1.4–11.9)1.4 (0.9–1.8)1.4 (1–2.2)3.1 (1.8–5.9)17.7 (11.9–23.8)2.6 (1.2–4.4)11.7 (7.7–17)5.7 (0.4–13.7)
Al dissolved (g/L)55 (13–128)18 (5–57)76.5 (5–195)63 (30–90)52.8 (8–140)47.7 (27–71)47.8 (22–93)50.8 (5–81)26.8 (5–83)
As dissolved (g/L)3.5 (1–5)3.7 (1.4–7.1)4.2 (3–5.7)4.3 (2.2–7.5)1.8 (1–2.7)0.7 (0.2–1)3.2 (2.7–3.7)1.9 (0.5–3)2.3 (1.4–3.5)
Cu dissolved (g/L)182.4 (0.5–617.4)38.2 (13.1–63.8)3.5 (0.5–8.9)8.6 (0.5–13.8)10.7 (5.7–16.2)3.7 (0.5–6.9)10.2 (4.3–17.5)10.2 (4.2–15.3)6.5 (0.5–13.8)
Fe dissolved (g/L)21.5 (3–35)16 (1–38)30.3 (18–45)18.3 (1–27)17 (5–36)12.7 (1–22)28.5 (17–38)11.8 (1–32)27.5 (5–70)
Mn dissolved (g/L)27.5 (6–47)235.3 (68–570)4.5 (1–7)36 (4–94)30.3 (20–46)8 (4–11)42.5 (3–129)20.3 (5–43)31.8 (3–63)
Mb dissolved (g/L)3.3 (1–5)20.3 (1–65)2.3 (1–6)2.3 (1–5)2.5 (1–6)1.3 (1–2)3.8 (1–7)1.8 (1–3)1.8 (1–3)
Zn dissolved (g/L)18.5 (6.7–34.9)28–9 (7.2–67.6)11.5 (2.5–24.3)7.7 (5–13.8)13.5 (0.9–34.5)9.8 (0.2–22.5)10 (5.1–17.1)9.9 (3.1–21.9)13 (0.2–39.5)
Sulfate (mg/L)70.7 (51–102)186 (93.5–337.6)60.4 (50–66.7)89.2 (55.4–116.8)176.2 (127–234.3)167.4 (152.7–176.4)161.3 (116.9–223.9)155.5 (148.2–161.4)248.6 (215.2–266.3)
Total hardness (mg/L)580.3 (103.6–1875.5)231.9 (177–324.1)160.4 (125.1–187.1)157.9 (106–217.3)288.6 (200.4–367)372.5 (330.9–404.7)233.1 (175.4–315)325.8 (279.6–362)800.2 (388.3–1961.5)
Alkalinity (mg/L)52.3 (45–65)46 (36–55)69.5 (65–73)51 (43–56)124.5 (95–145)137.8 (97–172.5)59.4 (50–72.5)138.8 (115–155)123.4 (52.5–159)
Chlorides (mg/L)18.3 (10–24.8)64.3 (24.1–125.1)19.9 (15.9–25.1)15.7 (9.1–23.9)49.1 (25.8–79.5)26.4 (26.2–26.6)42.1 (20.6–54.8)20.4 (20–20.8)53.6 (39.4–66.9)
Fecal coliforms (Nmp/100 mL)25.1 (4.5–46)190.8 (33–330)90.2 (7.8–220)44 (2–130)16.1 (4.5–26)1776.6 (240–4600)3025 (1100–7900)512.5 (170–920)892.3 (79–1700)
Total coliforms (Nmp/100 mL)65.4 (4.5–130)1780.8 (33–3300)1760 (280–4600)116.5 (17–330)493.3 (13–1300)5410 (330–11,000)39475 (7900–79,000)493.3 (13–1300)5280 (330–17,000)

References

  1. Ormerod, S.; Rundle, S.; Wilkinson, S.; Daly, G.; Dale, K.; Juttner, I. Altitudinal trends in the diatoms, bryophytes, macroinvertebrates and fish of a Nepalese river system. Freshw. Biol. 1994, 32, 309–322. [Google Scholar] [CrossRef]
  2. Vannote, R.L.; Minshall, G.W.; Cummins, K.W.; Sedell, J.R.; Cushing, C.E. The river continuum concept. Can. J. Fish. Aquat. Sci. 1980, 37, 130–137. [Google Scholar] [CrossRef]
  3. Fu, L.; Jiang, Y.; Ding, J.; Liu, Q.; Peng, Q.Z.; Kang, M.Y. Impacts of land use and environmental factors on macroinvertebrate functional feeding groups in the Dongjiang River basin, southeast China. J. Freshw. Ecol. 2016, 31, 21–35. [Google Scholar] [CrossRef] [Green Version]
  4. Allan, J.D. Landscapes and riverscapes: The influence of land use on stream ecosystems. Annu. Rev. Ecol. Evol. Syst. 2004, 35, 257–284. [Google Scholar] [CrossRef] [Green Version]
  5. Lawrence, J.E.; Lunde, K.B.; Mazor, R.D.; Bêche, L.A.; McElravy, E.P.; Resh, V.H. Long-term macroinvertebrate responses to climate change: Implications for biological assessment in mediterranean-climate streams. J. N. Am. Benthol. Soc. 2010, 29, 1424–1440. [Google Scholar] [CrossRef] [Green Version]
  6. Fierro, P.; Valdovinos, C.; Arismendi, I.; Díaz, G.; Jara-Flores, A.; Habit, E.; Vargas-Chacoff, L. Examining the influence of human stressors on benthic algae, macroinvertebrate, and fish assemblages in Mediterranean streams of Chile. Sci. Total Environ. 2019, 686, 26–37. [Google Scholar] [CrossRef] [PubMed]
  7. Martins, F.S.; Formigo, N.; Antunes, S.C. Can be the environmental and biotic factors responsible for macroinvertebrate communities’ alterations in Portuguese alpine ponds? Limnologica 2020, 83, 125782. [Google Scholar] [CrossRef]
  8. Firmiano, K.R.; Castro, D.M.; Linares, M.S.; Callisto, M. Functional responses of aquatic invertebrates to anthropogenic stressors in riparian zones of Neotropical savanna streams. Sci. Total Environ. 2021, 753, 141865. [Google Scholar] [CrossRef]
  9. Chen, K.; Hughes, R.M.; Xu, S.; Zhang, J.; Cai, D.; Wang, B. Evaluating performance of macroinvertebrate-based adjusted and unadjusted multi-metric indices (MMI) using multi-season and multi-year samples. Ecol. Indic. 2014, 36, 142–151. [Google Scholar] [CrossRef]
  10. Miserendino, M.; Brand, C.; Di Prinzio, C. Assessing urban impacts on water quality, benthic communities and fish in streams of the Andes Mountains, Patagonia (Argentina). Water Air Soil Pollut. 2008, 194, 91–110. [Google Scholar] [CrossRef]
  11. Fierro, P.; Bertrán, C.; Tapia, J.; Hauenstein, E.; Peña-Cortés, F.; Vergara, C.; Cerna, C.; Vargas-Chacoff, L. Effects of local land-use on riparian vegetation, water quality, and the functional organization of macroinvertebrate assemblages. Sci. Total Environ. 2017, 609, 724–734. [Google Scholar] [CrossRef]
  12. Hunke, P.; Mueller, E.N.; Schröder, B.; Zeilhofer, P. The Brazilian Cerrado: Assessment of water and soil degradation in catchments under intensive agricultural use. Ecohydrology 2015, 8, 1154–1180. [Google Scholar] [CrossRef]
  13. Brysiewicz, A.; Czerniejewski, P.; Bonisławska, M. Effect of Diverse Abiotic Conditions on the Structure and Biodiversity of Ichthyofauna in Small, Natural Water Bodies Located on Agricultural Lands. Water 2020, 12, 2674. [Google Scholar] [CrossRef]
  14. Atapattu, S.S.; Kodituwakku, D.C. Agriculture in South Asia and its implications on downstream health and sustainability: A review. Agric. Water Manag. 2009, 96, 361–373. [Google Scholar] [CrossRef]
  15. Bonada, N.; Prat, N.; Resh, V.H.; Statzner, B. Developments in aquatic insect biomonitoring: A comparative analysis of recent approaches. Annu. Rev. Entomol. 2006, 51, 495–523. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Novoa, V.; Ahumada-Rudolph, R.; Rojas, O.; Munizaga, J.; Sáez, K.; Arumí, J.L. Sustainability assessment of the agricultural water footprint in the Cachapoal River basin, Chile. Ecol. Indic. 2019, 98, 19–28. [Google Scholar] [CrossRef]
  17. Gerth, W.J.; Li, J.; Giannico, G.R. Agricultural land use and macroinvertebrate assemblages in lowland temporary streams of the Willamette Valley, Oregon, USA. Agric. Ecosyst. Environ. 2017, 236, 154–165. [Google Scholar] [CrossRef] [Green Version]
  18. Horak, C.N.; Assef, Y.A.; Grech, M.G.; Miserendino, M.L. Agricultural practices alter function and structure of macroinvertebrate communities in Patagonian piedmont streams. Hydrobiologia 2020, 847, 3659–3676. [Google Scholar] [CrossRef]
  19. Webb, M.J.; Winter, J.M.; Spera, S.A.; Chipman, J.W.; Osterberg, E.C. Water, agriculture, and climate dynamics in central Chile’s Aconcagua River Basin. Phys. Geogr. 2020, 1, 1–21. [Google Scholar] [CrossRef]
  20. Amigo, J.; Ramírez, C. A bioclimatic classification of Chile: Woodland communities in the temperate zone. Plant Ecol. 1998, 136, 9–26. [Google Scholar] [CrossRef]
  21. DGA. Diagnóstico y Clasificación de los Cursos y Cuerpos de agua Según Objetivos de Calidad: Cuenca del río Aconcagua; DGA: Valparaiso, Indiana, 2004. [Google Scholar]
  22. Carter, D.; Aguirre, L. Structural Geology of Aconcagua province and its relationship to the central Valley Graben, Chile. GSA Bull. 1965, 76, 651–664. [Google Scholar] [CrossRef]
  23. Wilgat, T.; Kosik, A.; Wojcirchowski, K. Hidrología de la cuenca Aconcagua selección de ejemplos. Inv. Geogr. 1968, 1–36. [Google Scholar] [CrossRef] [Green Version]
  24. Ribbe, L.; Delgado, P.; Salgado, E.; Flügel, W.A. Nitrate pollution of surface water induced by agricultural non-point pollution in the Pocochay watershed, Chile. Desalination 2008, 226, 13–20. [Google Scholar] [CrossRef]
  25. Gaete, H.; Aránguiz, F.; Cienfuegos, G.; Tejos, M. Metales pesados y toxicidad de aguas del río Aconcagua en Chile. Quim. Nova 2007, 30, 885–891. [Google Scholar] [CrossRef] [Green Version]
  26. Olivares, Y.; Gaete, H.; Neaman, A. Evaluación de la fitotoxicidad y la genotoxicidad de suelos agrícolas de zonas con actividades mineras de cobre de la cuenca del río Aconcagua (Chile central). Rev. Int. Contam. Ambient. 2015, 31, 237–243. [Google Scholar]
  27. Hernández, A.; Miranda, M.D.; Arellano, E.C.; Dobbs, C. Landscape trajectories and their effect on fragmentation for a Mediterranean semi-arid ecosystem in Central Chile. J. Arid Environ. 2016, 127, 74–81. [Google Scholar] [CrossRef]
  28. ESRI. GIS and Mapping Software. 2007. Available online: http://www.esri.com/ (accessed on 9 February 2021).
  29. Domínguez, E.; Fernández, H.R. Macroinvertebrados bentónicos sudamericanos. Sistemática y biología. Fundación Miguel Lillo, Tucumán, Argentina 2009, 103, 695. [Google Scholar]
  30. APHA. Standard Methods for the Examination of Water and Wastewater; American Public Health Association (APHA): Washington, DC, USA, 2005. [Google Scholar]
  31. Elgueta, A.; Górski, K.; Thoms, M.; Fierro, P.; Toledo, B.; Manosalva, A.; Habit, E. Interplay of geomorphology and hydrology drives macroinvertebrate assemblage responses to hydropeaking. Sci. Total Environ. 2021, 768, 144262. [Google Scholar] [CrossRef]
  32. Anderson, M.J.; Robinson, J. Generalized discriminant analysis based on distances. Aust. NZ J. Stat. 2003, 45, 301–318. [Google Scholar] [CrossRef]
  33. Anderson, M.; Gorley, R.; Clarke, K.P. For PRIMER: Guide to Software and Statistical Methods; PRIMER-E: Plymouth, UK, 2008. [Google Scholar]
  34. Heino, J. Taxonomic surrogacy, numerical resolution and responses of stream macroinvertebrate communities to ecological gradients: Are the inferences transferable among regions? Ecol. Indic. 2014, 36, 186–194. [Google Scholar] [CrossRef]
  35. Team, R.C. R: A Language and Environment for Statistical Computing [Computer Software Manual]; Team, R.C.: Vienna, Austria, 2016. [Google Scholar]
  36. Alavaisha, E.; Lyon, S.W.; Lindborg, R. Assessment of water quality across irrigation schemes: A case study of wetland agriculture impacts in Kilombero Valley, Tanzania. Water 2019, 11, 671. [Google Scholar] [CrossRef] [Green Version]
  37. Burdon, F.J.; Munz, N.A.; Reyes, M.; Focks, A.; Joss, A.; Räsänen, K.; Altermatt, F.; Eggen, R.I.L.; Stamm, C. Agriculture versus wastewater pollution as drivers of macroinvertebrate community structure in streams. Sci. Total Environ. 2019, 659, 1256–1265. [Google Scholar] [CrossRef]
  38. Usseglio-Polatera, P.; Beisel, J.N. Longitudinal changes in macroinvertebrate assemblages in the Meuse River: Anthropogenic effects versus natural change. River Res. Appl. 2002, 18, 197–211. [Google Scholar] [CrossRef]
  39. Jacobsen, D.; Wiberg-Larsen, P.; Brodersen, K.P.; Hansen, S.B.; Lindegaard, C.; Friberg, N.; Dall, P.C.; Kirkegaard, J.; Skriver, J.; Toman, M. Macroinvertebrate communities along the main stem and tributaries of a pre-Alpine river: Composition responds to altitude, richness does not. Limnologica 2020, 84, 125816. [Google Scholar] [CrossRef]
  40. Ruaro, R.; Gubiani, É.A.; Cunico, A.M.; Moretto, Y.; Piana, P.A. Comparison of fish and macroinvertebrates as bioindicators of Neotropical streams. Environ. Monit. Assess. 2016, 188, 45. [Google Scholar] [CrossRef]
  41. Harrington, R.A.; Poff, N.L.; Kondratieff, B.C. Aquatic insect β-diversity is not dependent on elevation in Southern Rocky Mountain streams. Freshw. Biol. 2016, 61, 195–205. [Google Scholar] [CrossRef]
  42. Molineri, C.; Nieto, C.; Dos Santos, D.A.; Emmerich, D.; Zúñiga, M.d.C.; Fierro, P.; Pessacq, P.; Gomez, D.; Márquez, J.A.; Príncipe, R.E.; et al. Do mayflies (Ephemeroptera) support a biogeographic transition zone in South America? J. Biogeogr. 2020, 47, 1980–1993. [Google Scholar] [CrossRef]
  43. Scheibler, E.E.; Fernández Campón, F.; Lagos Silnik, S.; Wellnitz, T. Seasonal influence and local factors affecting macroinvertebrate structure in a high-altitude Andean stream. J. Mt. Sci. 2020, 17, 1374–1386. [Google Scholar] [CrossRef]
  44. Dos Santos, D.A.; Molineri, C.; Nieto, C.; Zuñiga, M.C.; Emmerich, D.; Fierro, P.; Pessacq, P.; Rios-Touma, B.; Márquez, J.; Gomez, D.; et al. Cold/Warm stenothermic freshwater macroinvertebrates along altitudinal and latitudinal gradients in Western South America: A modern approach to an old hypothesis with updated data. J. Biogeogr. 2018, 45, 1571–1581. [Google Scholar] [CrossRef]
  45. Gil, M.; Tripole, S.; Vallania, E. Feeding habits of Smicridea (Rhyacophylax) dythyra Flint, 1974 (Trichoptera: Hydropsychidae) larvae in the Los Molles stream (San Luis–Argentina). Acta Limnol. Bras. 2008, 20, 1–4. [Google Scholar]
  46. Príncipe, R.E.; Gualdoni, C.M.; Raffani, G.B.; Corigliano, M.C. Spatial-temporal patterns of functional feeding groups in mountain streams of Córdoba, Argentina. Ecol. Austral 2010, 20, 257–268. [Google Scholar]
  47. Fierro, P.; Arismendi, I.; Hughes, R.M.; Valdovinos, C.; Jara-Flores, A. A benthic macroinvertebrate multimetric index for Chilean Mediterranean streams. Ecol. Indic. 2018, 91, 13–23. [Google Scholar] [CrossRef]
  48. Rodrigues, V.; Estrany, J.; Ranzini, M.; de Cicco, V.; Martín-Benito, J.M.T.; Hedo, J.; Lucas-Borja, M.E. Effects of land use and seasonality on stream water quality in a small tropical catchment: The headwater of Córrego Água Limpa, São Paulo (Brazil). Sci. Total Environ. 2018, 622, 1553–1561. [Google Scholar] [CrossRef] [Green Version]
  49. Wang, X.; Zhang, Y.; Tan, X.; Zheng, Y.; Zhang, Q. Do water quality, land use, or benthic diatoms drive macroinvertebrate functional feeding groups in a subtropical mountain stream? Inland Waters 2020, 1, 1–11. [Google Scholar] [CrossRef]
  50. Monaghan, K.A.; Peck, M.R.; Brewin, P.; Masiero, M.; Zarate, E.; Turcotte, P.; Ormerod, S. Macroinvertebrate distribution in Ecuadorian hill streams: The effects of altitude and land use. Arch. Hydrobiol. 2000, 149, 421–440. [Google Scholar] [CrossRef]
  51. Crespo, S.A.; Lavergne, C.; Fernandoy, F.; Muñoz, A.A.; Cara, L.; Olfos-Vargas, S. Where does the chilean aconcagua river come from? Use of natural tracers for water genesis characterization in glacial and periglacial environments. Water 2020, 12, 2630. [Google Scholar] [CrossRef]
  52. Climent, M.J.; Sánchez-Martín, M.J.; Rodríguez-Cruz, M.S.; Pedreros, P.; Urrutia, R.; Herrero-Hernández, E. Determination of pesticides in river surface waters of central Chile using SPE-GC-MS multi-residue method. J. Chil. Chem. Soc. 2018, 63, 4023–4031. [Google Scholar] [CrossRef] [Green Version]
  53. Cornejo, A.; Pérez, J.; López-Rojo, N.; Tonin, A.M.; Rovira, D.; Checa, B.; Jaramillo, N.; Correa, K.; Villarreal, A.; Villarreal, V.; et al. Agriculture impairs stream ecosystem functioning in a tropical catchment. Sci. Total Environ. 2020, 745, 140950. [Google Scholar] [CrossRef] [PubMed]
  54. Shuman, T.C.; Smiley, P.C.; Gillespie, R.B.; Gonzalez, J.M. Influence of Physical and Chemical Characteristics of Sediment on Macroinvertebrate Communities in Agricultural Headwater Streams. Water 2020, 12, 2976. [Google Scholar] [CrossRef]
  55. Wang, L.; Robertson, D.M.; Garrison, P.J. Linkages between nutrients and assemblages of macroinvertebrates and fish in wadeable streams: Implication to nutrient criteria development. Environ. Manag. 2007, 39, 194–212. [Google Scholar] [CrossRef]
  56. Heino, J.; Mykrä, H.; Hämäläinen, H.; Aroviita, J.; Muotka, T. Responses of taxonomic distinctness and species diversity indices to anthropogenic impacts and natural environmental gradients in stream macroinvertebrates. Freshw. Biol. 2007, 52, 1846–1861. [Google Scholar] [CrossRef]
  57. Córdova, S.; Gaete, H.; Aránguiz, F.; Figueroa, R. Evaluación de la calidad de las aguas del estero Limache (Chile central), mediante bioindicadores y bioensayos. Lat. Am. J. Aquat. Res. 2009, 37, 199–209. [Google Scholar] [CrossRef]
  58. Pizarro, J.; Vergara, P.M.; Rodríguez, J.A.; Valenzuela, A.M. Heavy metals in northern Chilean rivers: Spatial variation and temporal trends. J. Hazard. Mater. 2010, 181, 747–754. [Google Scholar] [CrossRef] [PubMed]
  59. Suresh, G.; Ramasamy, V.; Meenakshisundaram, V.; Venkatachalapathy, R.; Ponnusamy, V. Influence of mineralogical and heavy metal composition on natural radionuclide concentrations in the river sediments. Appl. Radiat. Isot. 2011, 69, 1466–1474. [Google Scholar] [CrossRef]
  60. Alvial, I.E.; Orth, K.; Durán, B.C.; Alvarez, E.; Squeo, F.A. Importance of geochemical factors in determining distribution patterns of aquatic invertebrates in mountain streams south of the Atacama Desert, Chile. Hydrobiologia 2013, 709, 11–25. [Google Scholar] [CrossRef]
Figure 1. Map of the study area and sampling sites on the Aconcagua River Basin. (a,b) Show the location in a geographical context in South America. (c) Show the basin with the elevations. Least-disturbed sites by agricultural activities are: JU, JU10, BL20, AC10; intermedium-disturbed sites by agricultural activities AC20, PU10; most-disturbed sites by agricultural activities PO10, AC30, AC40.
Figure 1. Map of the study area and sampling sites on the Aconcagua River Basin. (a,b) Show the location in a geographical context in South America. (c) Show the basin with the elevations. Least-disturbed sites by agricultural activities are: JU, JU10, BL20, AC10; intermedium-disturbed sites by agricultural activities AC20, PU10; most-disturbed sites by agricultural activities PO10, AC30, AC40.
Water 13 00492 g001
Figure 2. Map of land uses and sampling sites on the Aconcagua River Basin. Name of sampling stations are shown in Figure 1. The black circles represent the least-disturbed sites; lead squares the intermedium-disturbed sites; and white triangles the most-disturbed sites by agricultural activities.
Figure 2. Map of land uses and sampling sites on the Aconcagua River Basin. Name of sampling stations are shown in Figure 1. The black circles represent the least-disturbed sites; lead squares the intermedium-disturbed sites; and white triangles the most-disturbed sites by agricultural activities.
Water 13 00492 g002
Figure 3. Biological indicators based on macroinvertebrate assemblages. Boxplots show the maximum and minimum as well as the interquartile ranges (25–75%) with solids lines representing median values. p indicates level of significance among the levels of agricultural disturbance in the one-way ANOVA.
Figure 3. Biological indicators based on macroinvertebrate assemblages. Boxplots show the maximum and minimum as well as the interquartile ranges (25–75%) with solids lines representing median values. p indicates level of significance among the levels of agricultural disturbance in the one-way ANOVA.
Water 13 00492 g003
Figure 4. Ordination plots of Distance-based Redundancy Analysis (dbRDA) of macroinvertebrate assemblages to environmental variables between seasons of the year. Vectors indicate the direction of the effects of quantitative variables in the ordination plot. Abbreviations for environmental variables are described in Table 2.
Figure 4. Ordination plots of Distance-based Redundancy Analysis (dbRDA) of macroinvertebrate assemblages to environmental variables between seasons of the year. Vectors indicate the direction of the effects of quantitative variables in the ordination plot. Abbreviations for environmental variables are described in Table 2.
Water 13 00492 g004
Table 1. Catchment and habitat characteristics and land use at 9 sampling sites on the Aconcagua River basin, Chile. Elev = Elevation, CA = Catchment area, Bo = Boulder, Co = Cobble, Gra = Gravel.
Table 1. Catchment and habitat characteristics and land use at 9 sampling sites on the Aconcagua River basin, Chile. Elev = Elevation, CA = Catchment area, Bo = Boulder, Co = Cobble, Gra = Gravel.
Site CodeElev (m asl)CA (ha)Stream OrderSubstrate Type% Agriculture% Native Forest% Shrublands% Forest Plantation% Urban-and-Industrial% Water Body% No Vegetation
JU1657157793Bo/Co0.00.419.10.00.017.063.5
BL201521387293Bo/Co0.01.210.30.04.625.158.9
JU101448173123Bo/Co0.01.210.30.04.625.158.9
AC101096449454Gra/Co0.24.016.10.01.718.160.0
AC206501641625Gra/Co6.87.329.60.01.814.640.0
PU106041038164Gra7.510.529.00.01.110.741.2
PO10600438144Gra/Co19.117.655.90.02.61.23.5
AC301883985405Gra10.815.930.00.11.610.131.6
AC40125389085Gra13.021.328.90.42.48.225.8
Table 2. Environmental features at 9 sampling sites on the Aconcagua River basin, Chile. Values represent average of four seasons (August and November 2014, March and June 2015).
Table 2. Environmental features at 9 sampling sites on the Aconcagua River basin, Chile. Values represent average of four seasons (August and November 2014, March and June 2015).
VariableAbbreviationJUBL20JU 10AC 10AC 20PU 10PO 10AC 300AC 40
Temperature (°C)6.38.51713.417.015.615.216.714.4
Conductivity (µS/cm)Cond367.8720.8303.2326.5643.0653.4493.0621.3938.5
pHpH7.16.97.48.07.97.67.58.27.6
Turbidity (NTU)Tur32.036.814.547.850.53.085.35.316.3
Dissolved oxygen (mg/L)DO10.49.39.99.910.69.79.811.310.2
BDO 5 (mg/L)BOD1.15.81.11.02.91.02.01.41.5
COD (mg/L)COD5.615.55.56.314.75.29.811.310.2
OSS (mg/L)SSO1.622.02.51.87.11.011.81.12.5
ISS (mg/L)SSI18.98.04.716.537.31.048.32.215.2
Total phosphorus (µg/L)TP67.562.535.0132.5200.010.0187.550.065.0
Ammonium (µg/L)Am20.032.525.020.01067.520.0142.522.557.5
Nitrite (µg/L)Nitri22.5133.524.526.3852.533.387.336.5153.5
Nitrate (mg/L)Nitra1.04.41.41.43.117.72.611.75.7
Al dissolved (µg/L)Al55.018.076.563.052.847.747.850.826.8
As dissolved (µg/L)As3.53.74.24.31.80.73.21.92.3
Cu dissolved (µg/L)Cu182.438.23.58.610.73.710.210.26.5
Fe dissolved (µg/L)Fe21.516.030.318.317.012.728.511.827.5
Mn dissolved (µg/L)Mn27.5235.34.536.030.38.042.520.331.8
Mb dissolved (µg/L)Mb3.320.32.32.32.51.33.81.81.8
Zn dissolved (µg/L)Zn18.528.911.57.713.59.810.09.913.0
Sulfate (mg/L)Sulf70.7186.060.489.2176.2167.4161.3155.5248.6
Total hardness (mg/L)TH580.3231.9160.4157.9288.6372.5233.1225.8800.2
Alkalinity (mg/L)Alk52.346.069.551.0124.5137.859.4138.8123.4
Chlorides (mg/L)Cl18.364.319.915.749.126.442.120.453.6
Fecal coliforms (Nmp/100 mL)ColFe25.1190.890.244.016.11776.63025512.5892.3
Total coliforms (Nmp/100 mL)ColTo65.41780.81760116.5493.35410.039,475.0493.35280.0
Table 3. Summary of Canonical Analysis of Principal Coordinates (CAP) for the three agriculture intervention types, environmental datasets, and macroinvertebrate assemblages in Aconcagua River Basin.
Table 3. Summary of Canonical Analysis of Principal Coordinates (CAP) for the three agriculture intervention types, environmental datasets, and macroinvertebrate assemblages in Aconcagua River Basin.
Percentage Correct ClassificationsAmong Land Use Types Differences
CAPLeastIntermediateMostTraceP
Environmental variables
Winter10010066.670.5870.064
Spring10010066.670.5660.055
Summer10010066.670.5660.057
Fall10010066.670.5830.079
Macroinvertebrate assemblages
Winter5001001.4440.080
Spring1001001001.9430.002
Summer100066.671.4900.008
Fall75066.671.2160.094
Table 4. DistLM sequential test results giving significant environmental variables explain macroinvertebrates assemblage variation to winter and spring 2014 and summer and fall 2015. Variables selection was based on the option “Forward”. Adj.R 2 = cumulative adjusted R 2 .
Table 4. DistLM sequential test results giving significant environmental variables explain macroinvertebrates assemblage variation to winter and spring 2014 and summer and fall 2015. Variables selection was based on the option “Forward”. Adj.R 2 = cumulative adjusted R 2 .
VariableAdjust.R 2 Pseudo-FPCumulative Proportion
Winter
CDOM0.2892.8560.0130.289
Temperature0.5283.0290.0120.528
Fecal coliforms0.6571.8840.0520.657
As0.751.6460.130.757
Cl0.8431.6510.2110.843
Nitrates0.9111.5460.2670.912
Mn0.9691.9020.3220.970
Spring
Temperature0.4255.1740.0050.425
pH0.5892.4010.0180.589
Nitrates0.6971.7760.1390.697
Catchment area0.7931.8680.1520.793
Total coliforms0.8571.9710.1440.875
Fe0.9311.6340.2740.931
Nitrites0.9772.0910.330.977
Summer
Elevation0.2682.5690.010.268
Total coliforms0.4431.8880.0240.443
ColFe0.5861.7260.080.586
SSO0.7061.6320.1710.706
Catchment area0.8121.7040.1760.812
Zn0.8911.4390.3010.891
COD0.972.6350.2770.97
Fall
Elevation0.28102.34530.02100.2810
Total coliforms0.52112.50500.03100.5211
Nitrites0.70832.56700.04100.7083
pH0.81931.84300.12200.8193
Ammonium0.89821.54800.26200.8982
Fecal coliforms0.95991.54000.40500.9599
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fierro, P.; Valdovinos, C.; Lara, C.; Saldías, G.S. Influence of Intensive Agriculture on Benthic Macroinvertebrate Assemblages and Water Quality in the Aconcagua River Basin (Central Chile). Water 2021, 13, 492. https://doi.org/10.3390/w13040492

AMA Style

Fierro P, Valdovinos C, Lara C, Saldías GS. Influence of Intensive Agriculture on Benthic Macroinvertebrate Assemblages and Water Quality in the Aconcagua River Basin (Central Chile). Water. 2021; 13(4):492. https://doi.org/10.3390/w13040492

Chicago/Turabian Style

Fierro, Pablo, Claudio Valdovinos, Carlos Lara, and Gonzalo S. Saldías. 2021. "Influence of Intensive Agriculture on Benthic Macroinvertebrate Assemblages and Water Quality in the Aconcagua River Basin (Central Chile)" Water 13, no. 4: 492. https://doi.org/10.3390/w13040492

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

Article Metrics

Back to TopTop