Advantages of Geographically Weighted Regression for Modeling Benthic Substrate in Two Greater Yellowstone Ecosystem Streams

Stream habitat assessments are commonplace in fish management, and often involve nonspatial analysis methods for quantifying or predicting habitat, such as ordinary least squares regression (OLS). Spatial relationships, however, often exist among stream habitat variables. For example, water depth, water velocity, and benthic substrate sizes within streams are often spatially correlated and may exhibit spatial nonstationarity or inconsistency in geographic space. Thus, analysis methods should address spatial relationships within habitat datasets. In this study, OLS and a recently developed method, geographically weighted regression (GWR), were used to model benthic substrate from water depth and water velocity data at two stream sites within the Greater Yellowstone Ecosystem. For data collection, each site was represented by a grid of 0.1 m² cells, where actual values of water depth, water velocity, and benthic substrate class were measured for each cell. Accuracies of regressed substrate class data by OLS and GWR methods were calculated by comparing maps, parameter estimates, and determination coefficient r ². For analysis of data from both sites, Akaike’s Information Criterion corrected for sample size indicated the best approximating model for the data resulted from GWR and not from OLS. Adjusted r ² values also supported GWR as a better approach than OLS for prediction of substrate. This study supports GWR (a spatial analysis approach) over nonspatial OLS methods for prediction of habitat for stream habitat assessments.     

Application of geographically weighted regression model to analysis of spatiotemporal varying relationships between groundwater quantity and land use changes (case study: Khanmirza Plain,Iran)

Understanding the spatiotemporal relationships between land use/cover changes (LUCC) and groundwater resources is necessary for effective and efficient land use management. In this paper, geographically weighted regression (GWR) and ordinary least squares (OLS) models have been expanded to analyze varying spatial relationships between groundwater quantity changes and LUCC for three periods: 1987–2000, 2000–2010, and 1987–2010 in the Khanmirza Plain of southwestern Iran. For this purpose, TM images were used to generate LUCC (rainfed, irrigated, meadow, and bare lands). Groundwater quantity variables, including groundwater level changes (GLC) and groundwater withdrawal differences (GWD), were gathered from piezometric and agricultural wells data. The analysis of spatial autocorrelation (Moran’s I and local indicators of spatial association ) demonstrated that GWR has a better ability to model spatially varying data with very minimal clustering of residuals. The results R ² and corrected Akaike’s Information Criterion parameters revealed that the GWR has the lowest similarity in space and time in neighboring situations and it has the high ability to explain more variance in the LUCC as a function of the groundwater quantity changes. All results of the distribution of local R ² values from GWR confirm our assertion that there is a spatiotemporal relationship between types of land use and each of groundwater quantity variables within the region. According to the t test results from GWR, there are significant differences between the GLC and GWD and the land use types in different places of region in each of the three time series. The GWR results can help decision-makers to make appropriate decisions for future planning.     

Spatial and Seasonal Variations in Chlorophyll-Nutrient Relationships in The Shallow Hypertrophic Lake Manyas, Turkey

Regression and correlation analyses were used to predict responses of phytoplankton biomass (chlorophyll) (μg L⁻¹) to nitrate (NO₃) (mg L⁻¹), phosphate (PO₄) (mg L⁻¹) and ammonium (NH₄) (mg L⁻¹) dynamics in the shallow hypertrophic Lake Manyas, Turkey. Nutrient concentrations showed a descending gradient with distance, while chlorophyll concentrations showed an ascending gradient with the distance from the Sığırcı Inlet to the Karadere Outlet. Higher nutrient concentrations did always not coincide with higher chlorophyll concentrations. The results showed that regression models developed using seasonal data were more accurate in predicting chlorophyll concentrations than those developed using the pooled data from whole year (based on R ² and the difference between the measured and predicted values). The findings also revealed that within a single large shallow lake, chlorophyll-nutrient relationships might show significant variations spatially. The objective of this study was to determine the seasonal and spatial variations in the relationships between chlorophyll, nitrate, phosphate and ammonium in the shallow hypertrophic Lake Manyas, Turkey.     

滇池东南岸农业和富磷区入湖河流地表径流及污染特征

应用聚类分析与因子分析方法,通过8次常规监测,对滇池东南岸10条以农业面源和受磷矿开采区影响的入湖河流的地表径流及其水质污染特征进行了分析,并探讨了其空间差异性。在南岸选取降雨过程相同的3条河流,开展暴雨径流监测,探讨污染物在降雨过程中的流失特征。结果表明,新宝象河的平均流量为2.6 m3/s,占总入湖流量的26.5%;总氮、总磷、化学需氧量、悬浮物是滇池的主要污染指标,许多河流均已严重超标。河流水质在空间上可分为3类,具有明显的空间差异性。总氮、总磷、溶解磷、硝态氮对水质污染的贡献率达到了53.636%,氮、磷含量是河流水质污染的主要贡献因子。降雨条件下化学需氧量、悬浮物浓度增长迅速,流量、悬浮物与大多数水质指标均有相关性,磷矿开采对河流水质的影响在降雨条件下更加明显,其悬浮物浓度在降雨条件下比只受农业面源影响的河流最高高出1.9倍。     

Implications of climate change for evaporation from bare soils in a Mediterranean environment

The purpose of this study was to predict quantitative changes in evaporation from bare soils in the Mediterranean climate region of Turkey in response to the projections of a regional climate model developed in Japan (hereafter RCM). Daily RCM data for the estimation of reference evapotranspiration (ET _r) and soil evaporation were obtained for the periods of 1994–2003 and 2070–2079. Potential evaporation (E _p) from bare soils was calculated using the Penman–Monteith equation with a surface resistance of zero. Simulation of actual soil evaporation (E _a) was carried out using Aydin model (Aydin et al., Ecological Modelling 182:91–105, 2005) combined with Aydin and Uygur (2006, A model for estimating soil water potential of bare fields. In Proceedings of the 18th International Soil Meeting (ISM) on Soils Sustaining Life on Earth, Managing Soil and Technology, Sanliurfa, 477–480pp.) model of predicting soil water potential at the top surface layer of a bare soil, after performances of Aydin model (R ² = 94.0%) and Aydin and Uygur model (R ² = 97.6) were tested. The latter model is based on the relations among potential soil evaporation, hydraulic diffusivity, and soil wetness, with some simplified assumptions. Input parameters of the model are simple and easily obtainable such as climatic parameters used to compute the potential soil evaporation, average diffusivity for the drying soil, and volumetric water content at field capacity. The combination of Aydin and Aydin and Uygur models appeared to be useful in estimating water potential of soils and E _a from bare soils, with only a few parameters. Unlike ET _r and E _p projected to increase by 92 and 69 mm (equivalent to 8.0 and 7.3% increases) due to the elevated evaporative demand of the atmosphere, respectively, E _a from bare soils is projected to reduce by 50 mm (equivalent to a 16.5% decrease) in response to a decrease in rainfall by 46% in the Mediterranean region of Turkey by the 2070s predicted by RCM, and consequently, to decreased soil wetness in the future.     

Uncertainty in watershed response predictions induced by spatial variability of precipitation

Negligence to consider the spatial variability of rainfall could result in serious errors in model outputs. The objective of this study was to examine the uncertainty of both runoff and pollutant transport predictions due to the input errors of rainfall. This study used synthetic data to represent the “true” rainfall pattern, instead of interpolated precipitation. It was conducted on a synthetic case area having a total area of 20 km² with ten subbasins. Each subbasin has one rainfall gauge with synthetic precipitation records. Six rainfall storms with varied spatial distribution were generated. The average rainfall was obtained from all of the ten gauges by the arithmetic average method. The input errors of rainfall were induced by the difference between the actual rainfall pattern and estimated average rainfall. The results show that spatial variability of rainfall can cause uncertainty in modeling outputs of hydrologic, which would be transport to pollutant export predictions, when uniformity of rainfall is assumed. Since rainfall is essential information for predicting watershed responses, it is important to consider the properties of rainfall, particularly spatial rainfall variability, in the application of hydrologic and water quality models.     

Spatial and temporal estimation of soil loss for the sustainable management of a wet semi-arid watershed cluster

The ungauged wet semi-arid watershed cluster, Seethagondi, lies in the Adilabad district of Telangana in India and is prone to severe erosion and water scarcity. The runoff and soil loss data at watershed, catchment, and field level are necessary for planning soil and water conservation interventions. In this study, an attempt was made to develop a spatial soil loss estimation model for Seethagondi cluster using RUSLE coupled with ARCGIS and was used to estimate the soil loss spatially and temporally. The daily rainfall data of Aphrodite for the period from 1951 to 2007 was used, and the annual rainfall varied from 508 to 1351 mm with a mean annual rainfall of 950 mm and a mean erosivity of 6789 MJ mm ha^?1 h^?1 year^?1. Considerable variation in land use land cover especially in crop land and fallow land was observed during normal and drought years, and corresponding variation in the erosivity, C factor, and soil loss was also noted. The mean value of C factor derived from NDVI for crop land was 0.42 and 0.22 in normal year and drought years, respectively. The topography is undulating and major portion of the cluster has slope less than 10°, and 85.3 % of the cluster has soil loss below 20 t ha^?1 year^?1. The soil loss from crop land varied from 2.9 to 3.6 t ha^?1 year^?1 in low rainfall years to 31.8 to 34.7 t ha^?1 year^?1 in high rainfall years with a mean annual soil loss of 12.2 t ha^?1 year^?1. The soil loss from crop land was higher in the month of August with an annual soil loss of 13.1 and 2.9 t ha^?1 year^?1 in normal and drought year, respectively. Based on the soil loss in a normal year, the interventions recommended for 85.3 % of area of the watershed includes agronomic measures such as contour cultivation, graded bunds, strip cropping, mixed cropping, crop rotations, mulching, summer plowing, vegetative bunds, agri-horticultural system, and management practices such as broad bed furrow, raised sunken beds, and harvesting available water using farm ponds and percolation tanks. This methodology can be adopted for estimating the soil loss from similar ungauged watersheds with deficient data and for planning suitable soil and water conservation interventions for the sustainable management of the watersheds.     

Spatial Modeling of Mean Annual Temperature in Iran: Comparing Cokriging and Geographically Weighted Regression

Mapping spatial distribution of climatological parameters with a good degree of accuracy is crucial in environmental modeling and planning. Nowadays, there are various models to estimate and predict spatial variables in an area but some such as cokriging and geographically weighted regression (GWR) have got more attention from experts in this field. The objectives of this study are to evaluate and compare GWR with ordinary cokriging (OCK) techniques for estimating the mean annual air temperature (MAT) of Iran using European Centre for Medium-Range Weather Forecasts (ECMWF) data and auxiliary variables (e.g., longitude, latitude and altitude). The MAT-gridded data for Iran was collected in pixels during the time interval of 1987–2015 from the ERA-Interim re-analysis version of ECMWF. Validation results indicate that cokriging model with latitude and altitude for estimating MAT has the lowest MAE (0.0155), MBE (0.00085), RMSE (0.0251), and the highest NS (0.9999) in relation to other cokriging methods. On the other hand, GWR with altitude has better results than those of GWR with other auxiliary variables because of its MAE (0.1271), MBE (0.0124), RMSE (0.1760), and NS (0.9969). By comparing two mentioned methods, cokriging with latitude and altitude has provided the best performance in relation to GWR for prediction of MAT in Iran. To obtain accurate estimation of the spatial distribution of MAT, local residuals were analyzed. Results concluded that residuals of the OCK model have high spatial adaptations between the observed and predicted MAT data compared to the GWR model. Hence, OCK was a relatively optimum method for the estimation of MAT compared with GWR.

     

Land use impacts on river health of Uma Oya,Sri Lanka: implications of spatial scales

Human actions on landscapes are a principal threat to the ecological integrity of river ecosystems worldwide. Tropical landscapes have been poorly investigated in terms of the impact of catchment land cover alteration on water quality and biotic indices in comparison to temperate landscapes. Effects of land cover in the catchment at two spatial scales (catchment and site) on stream physical habitat quality, water quality, macroinvertebrate indices and community composition were evaluated for Uma Oya catchment in the upper Mahaweli watershed, Sri Lanka. The relationship between spatial arrangement of land cover in the catchment and water quality, macroinvertebrate indices and community composition was examined using univariate and multivariate approaches. Results indicate that chemical water quality variables such as conductivity and total dissolved solids are mostly governed by the land cover at broader spatial scales such as catchment scale. Shannon diversity index was also affected by catchment scale forest cover. In stream habitat features, nutrients such as N-NO₃ ^?, macroinvertebrate family richness, %shredders and macroinvertebrate community assemblages were predominantly influenced by the extent of land cover at 200 m site scale suggesting that local riparian forest cover is important in structuring macroinvertebrate communities. Thus, this study emphasizes the importance of services provided by forest cover at catchment and site scale in enhancing resilience of stream ecosystems to natural forces and human actions. Findings suggest that land cover disturbance effects on stream ecosystem health could be predicted when appropriate spatial arrangement of land cover is considered and has widespread application in the management of tropical river catchments.     

Sorption of carbon tetrachloride,ethylene dibromide,and trichloroethylene on soil and clay

The sorption of carbon tetrachloride, ethylene dibromide, and trichloroethylene in two silty clay loam soils and aluminium (Al³⁺) or calcium (Ca²⁺) saturated montmorillonite clay was studied. When the adsorbents were exposed to environmental levels of these chemicals (10 to 1000 ppb in water) the amounts of each of the chemicals sorbed were 6% or less of that available except for a 17% sorption of trichloroethylene by Al-saturated clay. In the case of the Ca-saturated clay, there was no apparent sorption of carbon tetrachloride or trichloroethylene. When soil sorption was normalized based on the soil organic carbon content (K _oc) a correlation was found between the K _oc, water solubility, and octanol/water partitioning coefficients of the chemicals. However, carbon tetrachloride did not behave according to with the predicted relationships.     

Field scale variability of cadmium and zinc in soil and barley

CaCl₂-extractable soil Cd and Zn contents have been suggested as a measure of bioavailability. To investigate the ability of this measure to reflect spatial patterns of Cd and Zn concentrations in barley (Hordeum vulgare L.) in an arable field, plant and soil samples were taken from a 0.5 ha area sandy soil contaminated with Cd and Zn. Cd and Zn contents in barley and yield were spatially variable. Yield was low, which may have been caused by Zn toxicity or atrazine turnover. For Cd, CaCl₂-extractable soil contents explained only 17% of the variation in Cd contents in grain, and for Zn no significant correlation was observed. Nevertheless, surface plots of CaCl₂-extractable soil contents and contents of barley grain illustrated their corresponding spatial patterns. Despite the poor linear correlation between CaCl₂-extractable soil-Cd and grain-Cd, a stochastic model for long term behaviour of Cd in field soils predicted observed variability in Cd contents of barley grain well from spatial variability of soil pH and organic matter content. The probabilistic model predicted behaviour of Cd in terms of probability, and was more appropriate than the deterministic approach.     

Integral assessment of climate impact on the transboundary Mesta River flow formation in Bulgaria

The present work considers the section of the Mesta River on Bulgarian territory using the integral method for evaluation of climate and anthropogenic impact on the river flow. The level of this impact is determined by the index K _i (flow module), the coefficient C _i for the deviation of the average annual water volume Q _i from the flow norm Q _o and the index h _i for the deviation of the average annual rainfall volume H _i from the average multi-annual rainfall volume H _o . The dynamics of the average annual flow Q _i at two typical hydrometric stations – Yakoruda and Khadzhidimovo, as well as the dynamics of the average annual rainfall for the Yakoruda station was examined for the period 1955–2003. The data for the considered period 1955–2003 exhibit a decreasing trend of the average annual water volumes dynamics for both stations due to the impact of climate changes in the Mesta River catchment.     

遂宁城乡集中式饮用水水源地钡分布特征及健康风险评价

为考察遂宁市辖区内集中式饮用水水源地污染物钡的分布特征和健康风险水平,通过电感耦合等离子体原子发射光谱法对研究区域内市级、县级和乡镇级所有在用的56个集中式饮用水水源地钡的浓度进行分析检测,借助空间分析与统计分析的结果,探讨了其空间分布和浓度差异,并利用环境健康风险评价模型,对不同类型水源地钡的健康风险进行了评价。结果表明,38个地表水水源地钡的浓度范围为0.065~0.180 mg/L,均值为0.110 mg/L;18个地下水水源地钡的浓度范围为0.027~0.370 mg/L,均值为0.130 mg/L。地表水与地下水水源地间钡的浓度差异具有统计学意义(P0.05),钡的空间分布也存在不同程度的差异性。各水源地中的钡经饮用和皮肤暴露两种途径对成人和儿童所引起的非致癌风险值为1.34×10~(-8)~1.62×10~(-8),远低于推荐的最大可接受风险水平(1.0×10~(-6)),各水源地因污染物钡导致的非致癌风险极低。     

Monitoring Water Quality and Quantity of National Watersheds in Turkey

National data from the hydrological network for 38 rivers out of 25 watersheds were used to detect spatial and temporal trends in water quality and quantity characteristics between 1995 and 2002. Assessment of water quality and quantity included flow rate, water temperature, pH, electrical conductivity, sodium adsorption rate, Na, K, Ca+Mg, CO₃, HCO₃, Cl, SO₄, and boron. Among the major ions assessed on a watershed basis, Turkish river waters are relatively high in Ca+Mg, Na and HCO₃, and low in K and CO₃. The watersheds in Turkey experienced a general trend of 16% decrease in flow rates between 1995 and 2002 at a mean annual rate of about 4 m³ s^?1, with a considerable spatial variation. Similarly, there appeared to be an increasing trend in river water temperature, at a mean annual rate of about 0.2°C. A substantial proportion of watersheds experienced an increase in pH, in particular, after 1997, with a maximum increase from 8.1 to 8.4 observed in Euphrates (P?R ² values in accounting for variations of pH and water temperature only. The findings of the study can provide a useful assessment of controls over water quality and quantity and assist in devising integrated and sustainable management practices for watersheds at the regional scale in Turkey.     

The relative influence of nutrients and habitat on stream metabolism in agricultural streams

Stream metabolism was measured in 33 streams across a gradient of nutrient concentrations in four agricultural areas of the USA to determine the relative influence of nutrient concentrations and habitat on primary production (GPP) and respiration (CR-24). In conjunction with the stream metabolism estimates, water quality and algal biomass samples were collected, as was an assessment of habitat in the sampling reach. When data for all study areas were combined, there were no statistically significant relations between gross primary production or community respiration and any of the independent variables. However, significant regression models were developed for three study areas for GPP (r ² = 0.79–0.91) and CR-24 (r ² = 0.76–0.77). Various forms of nutrients (total phosphorus and area-weighted total nitrogen loading) were significant for predicting GPP in two study areas, with habitat variables important in seven significant models. Important physical variables included light availability, precipitation, basin area, and in-stream habitat cover. Both benthic and seston chlorophyll were not found to be important explanatory variables in any of the models; however, benthic ash-free dry weight was important in two models for GPP.     

Effects of land use types on surface water quality across an anthropogenic disturbance gradient in the upper reach of the Hun River, Northeast China

Surface water quality is vulnerable to pollution due to human activities. The upper reach of the Hun River is an important water source that supplies 52 % of the storage capacity of the Dahuofang Reservoir, the largest reservoir for drinking water in Northeast China, which is suffering from various human-induced changes in land use, including deforestation, reclamation/farming, urbanization and mine exploitation. To investigate the impacts of land use types on surface water quality across an anthropogenic disturbance gradient at a local scale, 11 physicochemical parameters (pH, dissolved oxygen [DO], turbidity, oxygen redox potential, conductivity, biochemical oxygen demand [BOD₅], chemical oxygen demand [COD], total nitrogen [TN], total phosphorus [TP], NO ₃ ^? -N, and NH ₄ ⁺ -N) of water from 12 sampling sites along the upper reach of the Hun River were monitored monthly during 2009–2010. The sampling sites were classified into four groups (natural, near-natural, more disturbed, and seriously disturbed). The water quality exhibited distinct spatial and temporal characteristics; conductivity, TN, and NO ₃ ^? -N were identified as key parameters indicating the water quality variance. The forest and farmland cover types played significant roles in determining the surface water quality during the low-flow, high-flow, and mean-flow periods based on the results of a stepwise linear regression. These results may provide incentive for the local government to consider sustainable land use practices for water conservation.     

The temporal changes in road stormwater runoff quality and the implications to first flush control in Chongqing, China

This study investigates the quality of stormwater runoff from a driveway in the southwest mountainous urban area of Chongqing, China, from 2010 to 2011. The results showed that the mean concentrations of chemical oxygen demand (COD), total nitrogen (TN), and total phosphorus (TP) were 4.1, 2.4, and 2.2 times the grade V levels of the national surface water standard of China. The pollutant concentration peak preceded or synchronized with the rainfall intensity peak and occurred 10 min after the runoff started. The significant high pollutant concentration in the initial stage of the rainfall suggested that first flush control is necessary, especially for the most polluted constitutes, such as total suspended solids, COD, and TN. Three potential pollution sources were identified: the atmospheric dry and wet deposition (TN, NO₃ ^?-N, NH₄ ⁺-N, and DCu), the road sediment and materials (total suspended solids, COD, and TP), and the vehicle emissions (DPb and DZn). Therefore, this study indicates that reductions in road sediments and material pollution and dry and wet deposition should be the priority factors for pollution control of road stormwater runoff in mountainous urban areas.     

Fluctuations of Climatic Conditions, Elemental Cycling and Forest Growth at the Watershed Scale

The relationships between fluctuations in climatic conditions,forest productivity and elemental cycling were studied from 1994 to 1997 in a headwater catchment of the southern Laurentians dominated by sugar maple (Acer saccharumMarsh.) growing on podzolic soils. Annual budgets show that H⁺, K, and NO₃ were retained in the watershed whileCa, Mg and Na were lost. The magnitude of the net annual budget for Ca, Mg and Na was correlated to annual variationsin precipitation with the absolute budget value decreasing during dry years. Stemwood (r² = 0.85) and total tree biomass production (r² = 0.99) were correlated with mean annual temperature but fine roots and leaf litter werenot. During the growing season, the pH of the organic horizons(FH) decreased as the volumetric water content of soildecreased. A positive association was also found between airtemperature and H₂O-soluble (r² = 0.88) and PO₄-extractable (r² = 0.99) SO₄ in the upper B horizon. On a multi-year scale, we suspect that the decreasein the storage of inorganic SO₄ in the soil results from the cumulative effects of annual variations in climatic conditions superimposed on the long-term decrease in SO₄deposition from the atmosphere. These soil changes were reflected by a decline in SO₄, Ca and Mg concentrationsin the stream. The generalisation of the observed short-term patterns to longer time scales must be approached with caution. Yet, our results indicate that the associations between climatic variations and the biogeochemistry of the ecosystem occur at different spatial and temporal scales and integrate abroad range of chemical components and ecosystem compartments. This reflects the inherent complexity of natural systems and offers a vast palette of indicators of the response of terrestrial ecosystems to variations in the intensity of environmental factors such as climatic conditions.     

Does seston size structure reflect fish-mediated effects on water quality in a degraded semiarid wetland?

This study evaluates whether the size structure of seston (the sum of living and nonliving particles in the water column) reflects the effects of fish on wetland water quality. Using enclosures, we measured water quality and zooplankton community structure in the presence and absence of two fish species with distinct foraging strategies [benthivorous carp (Cyprinus carpio) and planktivorous mosquitofish (Gambusia holbrooki)]. Seston collected from the enclosures was counted and sized automatically with a Coulter counter, and the size structure in the range of 1–60 μm was modelled using the underlying Pareto distribution of particles. Only Cyprinus contributed to a loss of water quality (increased nutrient levels, algal and non-algal turbidity, hypoxia), while both fish species changed zooplankton community composition compared to fishless controls. However, these changes were not reflected in parameters [goodness of fit (r ²) and parameter c (slopes)] of Pareto models. Multivariate statistics suggest that seston size structure responded more to environmental gradients related to water depth but the relationships were also weak. Our overall result contrasts with the regulation of size spectra constructed from living organisms. Although seston integrates many structural and functional features of the water column, the lack of strong relationships between Pareto model parameters and water quality suggests that a size-based approach to characterise seston structure has a limited potential for assessing biota-mediated effects in wetlands in a straightforward manner.     

Impact of water quality on bacterioplankton assemblage along Cértima River Basin (central western Portugal) assessed by PCR–DGGE and multivariate analysis

The information on bacterial community composition (BCC) in Portuguese water bodies is very scarce. Cértima River (central western Portugal) is known to have high levels of pollution, namely organic. In the present work, the BCC from a set of 16 water samples collected from Cértima River Basin and its main tributaries was characterized using 16S rDNA–denaturing gradient gel electrophoresis, a culture-independent molecular approach. Molecular data were related to environmental parameters through multivariate analysis to investigate potential impact of water pollution along the river. Principal component analysis using environmental data showed a water quality gradient from more pristine waters (at the mountain tributaries) to waters with increasingly eutrophic potential (such as Fermentelos Lake). This gradient was mainly defined by factors such as organic and inorganic nutrient sources, electrical conductivity, hydrogen carbonate concentration, and pH. Molecular results showed variations in BCC along Cértima River Basin but in the main river section, a Bacteroidetes phylotype (Flavobacterium sp.) proved to be dominant throughout the river course. Multivariate analysis suggests that spatial variation of BCC along the Cértima River Basin depended mainly on parameters such as Chl a, total suspended solid (TSS), total organic carbon, electrical conductivity, and HCO _boxclose^-_{3}^{-} levels. Bacteroidetes phylotypes were all related to higher electrical conductivity and HCO₃^-_{3}^{-} levels although some of these were also correlated with high SO₄^2-_{4}^{2-} and others with high soluble reactive phosphorus, nitrate, TN, and Kjeld-N levels. The Gammaproteobacteria occurrence was correlated with high SO₄^2-_{4}^{2-} levels. One of the Betaproteobacteria phylotypes showed to correlate with low redox potential (E_h) and high temperature, pH, TSS, and Chl a levels while another one showed a negative correlation with Chl a values.