Influence of spatial resolution on assessing channelization impacts on fish and macroinvertebrate communities in a warmwater stream in the southeastern United States

The choice of spatial and temporal scale used in environmental assessments may influence the observed results. One method of assessing the impact of stream habitat alterations involves the comparison of response variables among treatment categories (i.e., impacted and unimpacted sites). The influence of spatial resolution on patterns of response variables among treatment categories in assessments of stream channelization and other types of habitat alterations has not been evaluated. We examined how patterns of 10 community response variables among channel types and our interpretations of channelization impacts on fish and macroinvertebrate communities differed among three spatial resolutions in a warmwater stream in Mississippi and Alabama. Four fish and three macroinvertebrate community response variables exhibited different patterns among channel types at different spatial resolutions. Our interpretations of the impacts of channelization on fish and macroinvertebrate communities differed among spatial resolutions. Channelization had a negative influence on fish communities either with or without evidence of potential community recovery in one channel type. Channelization impacts on macroinvertebrate communities ranged from a negative influence to no effect. Our results suggest that spatial resolution can influence the observed results and interpretations derived from assessments of stream habitat alterations.     

Correspondence of biological condition models of California streams at statewide and regional scales

We used boosted regression trees (BRT) to model stream biological condition as measured by benthic macroinvertebrate taxonomic completeness, the ratio of observed to expected (O/E) taxa. Models were developed with and without exclusion of rare taxa at a site. BRT models are robust, requiring few assumptions compared with traditional modeling techniques such as multiple linear regression. The BRT models were constructed to provide baseline support to stressor delineation by identifying natural physiographic and human land use gradients affecting stream biological condition statewide and for eight ecological regions within the state, as part of the development of numerical biological objectives for California’s wadeable streams. Regions were defined on the basis of ecological, hydrologic, and jurisdictional factors and roughly corresponded with ecoregions. Physiographic and land use variables were derived from geographic information system coverages. The model for the entire state (n?=?1,386) identified a composite measure of anthropogenic disturbance (the sum of urban, agricultural, and unmanaged roadside vegetation land cover) within the local watershed as the most important variable, explaining 56 % of the variance in O/E values. Models for individual regions explained between 51 and 84 % of the variance in O/E values. Measures of human disturbance were important in the three coastal regions. In the South Coast and Coastal Chaparral, local watershed measures of urbanization were the most important variables related to biological condition, while in the North Coast the composite measure of human disturbance at the watershed scale was most important. In the two mountain regions, natural gradients were most important, including slope, precipitation, and temperature. The remaining three regions had relatively small sample sizes (n?≤?75 sites) and had models that gave mixed results. Understanding the spatial scale at which land use and land cover affect taxonomic completeness is imperative for sound management. Our results suggest that invertebrate taxonomic completeness is affected by human disturbance at the statewide and regional levels, with some differences among regions in the importance of natural gradients and types of human disturbance. The construction and application of models similar to the ones presented here could be useful in the planning and prioritization of actions for protection and conservation of biodiversity in California streams.     

Influence of Riffle and Snag Habitat Specific Sampling on Stream Macroinvertebrate Assemblage Measures in Bioassessment

Stream macroinvertebrate communities vary naturally among types of habitats where they are sampled, which affects the results of environmental assessment. We analyzed macroinvertebrates collected from riffle and snag habitats to evaluate influences of habitat-specific sampling on taxon occurrence, assemblage measures, and biotic indices. We found considerably more macroinvertebrate taxa unique to snags (143 taxa) than to riffles (75 taxa), and the numbers of taxa found in both riffles and snags (149 taxa) were similar to that found in snags. About 64% of the 47 macroinvertebrate measures we tested differed significantly between riffles and snags. Eighty percent intercepts of regressions between biotic indices and urban or agricultural land uses differed significantly between riffles and snags. The Hilsenhoff biotic index calculated from snag samples explained 69% of the variance of riffle samples and classified 66% of the sites into the same stream health group as the riffle samples. However, four multimetric indices for snag samples explained less than 50% of the variance of riffle samples and classified less than 50% of the sites into the same health group as the riffle samples. We concluded that macroinvertebrate indices developed for riffle/run habitat should not be used for snag samples to assess stream impairment. We recommend developing an index of biotic integrity specifically for snags and using snags as an alternate sampling substrate for streams that naturally lack riffles.     

Reduced riparian zone width compromises aquatic macroinvertebrate communities in streams of southern Brazil

Recent changes in Brazilian legislation reduced the width of riparian forest buffer needed to be preserved in private properties from 30 to 15 m or less. The consequences of these modifications can be dramatic, mainly because riparian buffer width is an important parameter for riparian forest structure and functioning. Our study assessed whether (1) macroinvertebrate family richness and Ephemeroptera, Plecoptera, and Trichoptera (EPT) family richness decrease with reduced riparian buffer width; (2) taxonomic composition and functional feeding group (FFG) composition of macroinvertebrates vary with a reduced riparian buffer width; and (3) reduced riparian buffer width similarly influence the macroinvertebrate community in different stream substrates. We selected three fragments with different riparian buffer widths (>40, <30, and <15 m) in three streams (fourth and fifth orders) in the Sinos River watershed, southern Brazil. Our results show that on all substrate types, reducing the width of the riparian buffer altered neither the macroinvertebrate richness nor EPT richness. However, EPT richness was greater in the substrates stone and gravel than leaf litter, independent of riparian buffer width. There was a significant difference in macroinvertebrate composition among riparian buffer widths. The macroinvertebrate composition and FFG differed among substrates, independent of riparian buffer width. This study showed that riparian buffer widths <15 m altered the macroinvertebrate community. A width greater than 15 m is necessary to maintain the composition and trophic conditions of macroinvertebrate families similar to those found in reference states of conservation.     

Can Basin Land Use Effects on Physical Characteristics of Streams Be Determined at Broad Geographic Scales?

The environmental setting (e.g., climate, topography, geology) and land use affect stream physical characteristics singly and cumulatively. At broad geographic scales, we determined the importance of environmental setting and land use in explaining variation in stream physical characteristics. We hypothesized that as the spatial scale decreased from national to regional, land use would explain more of the variation in stream physical characteristics because environmental settings become more homogeneous. At a national scale, stepwise linear regression indicated that environmental setting was more important in explaining variability in stream physical characteristics. Although statistically discernible, the amount of variation explained by land use was not remarkable due to low partial correlations. At level II ecoregion spatial scales (southeastern USA plains, central USA plains, and a combination of the western Cordillera and the western interior basins and ranges), environmental setting variables were again more important predictors of stream physical characteristics, however, as the spatial scale decreased from national to regional, the portion of variability in stream physical characteristics explained by basin land use increased. Development of stream habitat indicators of land use will depend upon an understanding of relations between stream physical characteristics and environmental factors at multiple spatial scales. Smaller spatial scales will be necessary to reduce the confounding effects of variable environmental settings before the effects of land use can be reliably assessed.     

Benthic Macroinvertebrate Responses to Increasing Levels of Cattle Grazing in Blue Ridge Mountain Streams, Virginia, USA

The relationship between benthic macroinvertebrate assemblages and cattle density was assessed from fall 2002 through spring 2004 in five small streams that represented a gradient of cattle grazing intensity. All study stream reaches were in pasture with no woody riparian vegetation, but varied in the intensity of cattle grazing (0 cattle ha⁻¹ at site 1 to 2.85 cattle ha⁻¹ at site 5). Regression analysis indicated highly significant and strong macroinvertebrate metric responses to cattle density during most sampling periods. The majority of metrics responded negatively to increased grazing, while a few (total taxa richness, number of sensitive taxa, and % collector filterers) increased along the gradient before declining at the most heavily grazed sites. Total number of sensitive taxa and % Coleoptera had the strongest relationship with cattle density throughout the study period. During some sampling periods, nearly 80% of the variation in these metrics was explained by cattle density. The elmid beetle, Oulimnius, had a particularly strong negative response to the grazing gradient. Study site groupings based on taxa composition, using detrended correspondence analysis (DCA), indicated that benthic samples collected from the reference site and light rotational grazing site were more similar in macroinvertebrate taxa composition than samples collected from the intermediate grazing and heavy grazing sites. Our findings demonstrate that biological integrity, as measured by benthic macroinvertebrate metrics and assemblage composition, is highly related to cattle density in small streams in the Blue Ridge mountains, Virginia, USA. This suggests that the degree of agricultural intensity should be given consideration in stream assessments, as well as land use planning and regulatory decisions.     

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.     

Local physical habitat quality cloud the effect of predicted pesticide runoff from agricultural land in Danish streams

The combination of intensive agricultural activities and the close connectivity between land and stream emphasise the potential risk of pesticide exposure in Danish streams. Benthic macroinvertebrates are applied in the assessment of stream ecological status, and some sensitive species have been shown to respond strongly to brief pulses of pesticide contamination. In this study we investigate the impact of agriculturally derived pesticides on stream macroinvertebrate communities in Denmark. As a measure of toxic pressure we apply the Runoff Potential. We investigated a total of 212 streams. These were grouped into distinct classes according to the magnitude of pesticide contamination in the period from 2003-2006. A total of 24 different macroinvertebrate indices were applied to detect effects of pesticide runoff (e.g. the SPEAR-index and the number of EPT taxa). We found high predicted pesticide runoff in 39% of the streams, but we found no significant effect of predicted pesticide exposure on stream macroinvertebrate indices. We, additionally, examined the influence of a series of environmental parameters ranging from site scale to catchment scale on the macroinvertebrate community. Relative proportions of gravel, sand and silt in bed sediments explained most of the variation in macroinvertebrate indices as well as the upstream riparian habitat quality. We suggest that the Runoff Potential model overestimate pesticide runoff contamination in Danish streams due the presence of buffer strips enforced by Danish legislation. When pesticide runoff contamination is low to moderate, poor physical properties (indirectly related to agricultural activity) are the main impediment for the ecological quality of Danish streams.     

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.     

Limitations of habitats as biodiversity surrogates for conservation planning in estuaries

Increasing pressures on global biodiversity and lack of data on the number and abundance of species have motivated conservation planners and researchers to use more readily available information as proxies or surrogates for biodiversity. “Habitat” is one of the most frequently used surrogates but its assumed value in marine conservation planning is not often tested. The present study developed and tested three alternative habitat classification schemes of increasing complexity for a large estuary in south-east Australia and tested their effectiveness in predicting spatial variation in macroinvertebrate biodiversity and selecting estuarine protected areas to represent species. The three habitat classification schemes were: (1) broad-scale habitats (e.g., mangroves and seagrass), (2) subdivision of each broad-scale habitat by a suite of environmental variables that varied significantly throughout the estuary, and (3) subdivision of each broad-scale habitat by the subset of environmental variables that best explained spatial variation in macroinvertebrate biodiversity. Macroinvertebrate assemblages differed significantly among the habitats in each classification scheme. For each classification scheme, habitat richness was significantly correlated with species richness, total density of macroinvertebrates, assemblage dissimilarity, and summed irreplaceability. However, in a reserve selection process designed to represent examples of each habitat, no habitat classification scheme represented species significantly better than a random selection of sites. Habitat classification schemes may represent variation in estuarine biodiversity; however, the results of this study suggest they are inefficient in designing representative networks of estuarine protected areas.     

不同土地利用对溪流大型底栖无脊椎动物群落的影响

2010年4月调查了钱塘江中游区域29个样点的水环境特征和底栖动物。聚类排序将样点分为参照、农业和城镇3组,相似性分析表明不同组间底栖动物群落有显著差异(r=0.863,P=0.001)。环境因子的主成分分析表明,研究区域主要的环境胁迫是农业和城镇用地及其引起的水质变化,并能较好解释组间物种差异的关键环境胁迫因子为农业用地比例、城镇用地比例、溶解氧、总氮和平均底质得分(Rho=0.568,P=0.001)。底栖动物参数(总分类单元数、Shannon-Weaver多样性指数、BI指数和丰富度指数)和k-优势度曲线显示农业和城镇组的生物完整性遭到很大程度的破坏,且农业组较城镇组严重。     

Effects of geomorphology, habitat, and spatial location on fish assemblages in a watershed in Ohio, USA

In this paper, we evaluate relationships between in-stream habitat, water chemistry, spatial distribution within a predominantly agricultural Midwestern watershed and geomorphic features and fish assemblage attributes and abundances. Our specific objectives were to: (1) identify and quantify key environmental variables at reach and system wide (watershed) scales; and (2) evaluate the relative influence of those environmental factors in structuring and explaining fish assemblage attributes at reach scales to help prioritize stream monitoring efforts and better incorporate all factors that influence aquatic biology in watershed management programs. The original combined data set consisted of 31 variables measured at 32 sites, which was reduced to 9 variables through correlation and linear regression analysis: stream order, percent wooded riparian zone, drainage area, in-stream cover quality, substrate quality, gradient, cross-sectional area, width of the flood prone area, and average substrate size. Canonical correspondence analysis (CCA) and variance partitioning were used to relate environmental variables to fish species abundance and assemblage attributes. Fish assemblages and abundances were explained best by stream size, gradient, substrate size and quality, and percent wooded riparian zone. Further data are needed to investigate why water chemistry variables had insignificant relationships with IBI scores. Results suggest that more quantifiable variables and consideration of spatial location of a stream reach within a watershed system should be standard data incorporated into stream monitoring programs to identify impairments that, while biologically limiting, are not fully captured or elucidated using current bioassessment methods.     

Adjusting the effect of seasonal variability in the bioassessment of streams

Bioassessment tools should distinguish between the effects of anthropogenic degradation in communities and natural temporal changes. The present study tests the influence of natural seasonal variability on macroinvertebrate stream communities assessed by a predictive model (PORTRIV) and a multimetric index (IPtI) calibrated for spring. The scores of PORTRIV decreased significantly between spring and autumn, and between spring and winter (ca. 37 to 53 %, respectively), while those of IPtI did not change significantly between seasons. For non-reference samples, the results of the predictive model also indicate no significant differences. A correction factor (CF) was calculated to adjust the existing differences in the model assessments between seasons, based on the percentage of variation of reference site scores from spring to autumn and winter. After the application of the CF to the OE50 scores of spring reference samples, the differences were no longer significant. Independent reference validation sites confirmed this tendency. This method has the advantage of avoiding large efforts required for the construction of databases from other seasons and the development of new models to allow the assessment of streams in seasons other than spring. Further tests with models developed in regions with more marked seasonal changes should be done to confirm its wider applicability.     

Assessing the ecological condition of streams in a southeastern Brazilian basin using a probabilistic monitoring design

Prompt assessment and management actions are required if we are to reduce the current rapid loss of habitat and biodiversity worldwide. Statistically valid quantification of the biota and habitat condition in water bodies are prerequisites for rigorous assessment of aquatic biodiversity and habitat. We assessed the ecological condition of streams in a southeastern Brazilian basin. We quantified the percentage of stream length in good, fair, and poor ecological condition according to benthic macroinvertebrate assemblage. We assessed the risk of finding degraded ecological condition associated with degraded aquatic riparian physical habitat condition, watershed condition, and water quality. We describe field sampling and implementation issues encountered in our survey and discuss design options to remedy them. Survey sample sites were selected using a spatially balanced, stratified random design, which enabled us to put confidence bounds on the ecological condition estimates derived from the stream survey. The benthic condition index indicated that 62 % of stream length in the basin was in poor ecological condition, and 13 % of stream length was in fair condition. The risk of finding degraded biological condition when the riparian vegetation and forests in upstream catchments were degraded was 2.5 and 4 times higher, compared to streams rated as good for the same stressors. We demonstrated that the GRTS statistical sampling method can be used routinely in Brazilian rain forests and other South American regions with similar conditions. This survey establishes an initial baseline for monitoring the condition and trends of streams in the region.     

European Farmland Bird Distribution Explained by Remotely Sensed Phenological Indices

Birds are important components of biodiversity conservation since they are capable of indicating changes in the general status of wildlife and of the countryside. The Pan-European Common Bird Monitoring Scheme (PECBM) has been launched by the BirdLife Partnership in Europe, where the European Bird Census Council has been collecting data from 20 independent breeding bird survey programs across Europe over the last 25 years. These data show dramatic declines in European farmland birds. We suggest that seasonal characteristics of vegetation cover derived from high temporal resolution remote sensing images could facilitate the monitoring the suitability of farmland bird habitats, and that these indicators may be a better choice for monitoring than climate data. We used redundancy analysis to link the PECBM data of the estimated number of farmland birds in Europe to a set of phenological and climatic indicators and to the biogeographic regions of Europe. Variance partitioning was used to account for the variation explained by the phenological and climate variables and by the area of the environmental strata individually, to define the pure effect of the variables, and to extract the total explained variance. The analysis revealed high statistical significance (p < 0.001) of the correlations between species and environment. Phenological indices explained 38% of the variance in community composition of the 23 farmland bird species, whereas climate explained 30% of the variance. After partitioning the other variables as covariables, the pure effect of phenology, climate, and environmental strata were 16%, 8%, and 16%, respectively. Based on the probability results, we suggest that phenological indicators derived from remote sensing may supply better indicators for continental scale biodiversity studies than climate only. In addition, these indicators are cost and time effective, are on continuous scale, and are readily repeatable on a large spatial coverage while supplying standardized results.     

Assessment of wadeable stream resources in the driftless area ecoregion in Western Wisconsin using a probabilistic sampling design

The Wisconsin Department of Natural Resources (WDNR), with support from the U.S. EPA, conducted an assessment of wadeable streams in the Driftless Area ecoregion in western Wisconsin using a probabilistic sampling design. This ecoregion encompasses 20% of Wisconsin’s land area and contains 8,800 miles of perennial streams. Randomly-selected stream sites (n = 60) equally distributed among stream orders 1–4 were sampled. Watershed land use, riparian and in-stream habitat, water chemistry, macroinvertebrate, and fish assemblage data were collected at each true random site and an associated “modified-random” site on each stream that was accessed via a road crossing nearest to the true random site. Targeted least-disturbed reference sites (n = 22) were also sampled to develop reference conditions for various physical, chemical, and biological measures. Cumulative distribution function plots of various measures collected at the true random sites evaluated with reference condition thresholds, indicate that high proportions of the random sites (and by inference the entire Driftless Area wadeable stream population) show some level of degradation. Study results show no statistically significant differences between the true random and modified-random sample sites for any of the nine physical habitat, 11 water chemistry, seven macroinvertebrate, or eight fish metrics analyzed. In Wisconsin’s Driftless Area, 79% of wadeable stream lengths were accessible via road crossings. While further evaluation of the statistical rigor of using a modified-random sampling design is warranted, sampling randomly-selected stream sites accessed via the nearest road crossing may provide a more economical way to apply probabilistic sampling in stream monitoring programs.     

Establishment of a structural equation model for ground-level ozone: a case study at an urban roadside site

This study established a cause–effect relationship between ground-level ozone and latent variables employing partial least-squares analysis at an urban roadside site in four distinct seasons. Two multivariate analytic methods, factor analysis, and cluster analysis were adopted to cite and identify suitable latent variables from 14 observed variables (i.e., meteorological factors, wind and primary air pollutants) in 2008–2010. Analytical results showed that the first six components explained 80.3 % of the variance, and eigenvalues of the first four components were greater than 1. The effectiveness of this model was empirically confirmed with three indicators. Except for surface pressure, factor loadings of observed variables were 0.303–0.910 and reached statistical significance at the 5 % level. Composite reliabilities for latent variables were 0.672–0.812 and average variances were 0.404–0.547, except for latent variable “primary” in spring; thus, discriminant validity and convergent validity were marginally accepted. The developed model is suitable for the assessment of urban roadside surface ozone, considering interactions among meteorological factors, wind factors, and primary air pollutants in each season.     

Community analysis in stream biomonitoring: what we measure and what we don't

Diatom assemblages from 83 epilithic samples taken from the Mesta River, Bulgaria, were regressed against three sets of predictor variables, i.e. environmental, spatial, and temporal. Redundancy analysis (RDA) of species and environmental data explained 36% of the diatom variance and extracted several important gradients of species distribution, associated with a downstream increase in nutrient levels, pH, temperature, and organic pollution. The inclusion of spatial and temporal variables in the RDA model captured additional 24% of the diatom variance and revealed three more gradients, a spatial gradient represented by higher order polynomial terms of latitude and longitude, and two temporal gradients of annual and seasonal variation. Partial RDAs demonstrated that the unique contribution of each predictor set to the explained diatom variance was the highest in the spatial dataset (16%), followed by the environmental (9%), and the temporal (7%) datasets. The remaining 28% of the variance was explained by the covariance of the predictor sets. This suggests that in biomonitoring of single stream basins, the cheap and simple account of space and time would explain most of the variance in assemblage composition obviating the necessity of expensive and time-consuming environmental assessments. The nature of the underlying environmental mechanisms can be easily inferred from the diatom composition itself.     

A risk assessment study of water quality,biota, and legacy sediment prior to small dam removal in a tributary to the Delaware River

The proposed removal of three run-of-river dams (all ≤5-m height) in eastern Pennsylvania along lower Bushkill Creek, a tributary of the Delaware River, has provided a valuable opportunity for multidisciplinary research involving the collection of more than 5 years of pre-removal monitoring data, analysis of heavy metals in legacy sediment cores, and associated toxicity assays to determine the singular and interactive effects of lead, copper, and cadmium on survival and behavior of a common macroinvertebrate found in Bushkill Creek. Monitoring data were collected from sites approximately 35 m upstream and downstream of dams and reference sites located approximately 5 km upstream of all dams. Results indicate that oxygen levels, macroinvertebrate diversity, and proportion of sensitive taxa were significantly lower upstream and downstream of dams in comparison with upstream reference reaches. The strong correlation between water quality and macroinvertebrates in this system implies that removal of the lower three dams would lead to improvements in water quality, biotic integrity, and resilience in lower Bushkill Creek. Sediment analyses and toxicity assays suggest that dam removal and sediment mobilization may route contaminated sediments downstream at concentrations that may harm more sensitive biota. However, macroinvertebrate mortality and behavior were not significantly different from clean water controls for the large majority of toxicity assays. All together, these results suggest that dams 1–3 are good candidates for successful stream restoration but that the removals would best be planned in a way that mitigates potential impacts of contaminated legacy sediment.     

Land use impact on soil quality in eastern Himalayan region of India

Quantitative assessment of soil quality is required to determine the sustainability of land uses in terms of environmental quality and plant productivity. Our objective was to identify the most appropriate soil quality indicators and to evaluate the impact of six most prevalent land use types (natural forestland, cultivated lowland, cultivated upland terrace, shifting cultivation, plantation land, and grassland) on soil quality in eastern Himalayan region of India. We collected 120 soil samples (20 cm depth) and analyzed them for 29 physical, chemical, and biological soil attributes. For selection of soil quality indicators, principal component analysis (PCA) was performed on the measured attributes, which provided four principal components (PC) with eigenvalues >1 and explaining at least 5 % of the variance in dataset. The four PCs together explained 92.6 % of the total variance. Based on rotated factor loadings of soil attributes, selected indicators were: soil organic carbon (SOC) from PC-1, exchangeable Al from PC-2, silt content from PC-3, and available P and Mn from PC-4. Indicators were transformed into scores (linear scoring method) and soil quality index (SQI) was determined, on a scale of 0–1, using the weighting factors obtained from PCA. SQI rating was the highest for the least-disturbed sites, i.e., natural forestland (0.93) and grassland (0.87), and the lowest for the most intensively cultivated site, i.e., cultivated upland terrace (0.44). Ratings for the other land uses were shifting cultivation (0.60)?>?cultivated low land (0.57)?>?plantation land (0.54). Overall contribution (in percent) of the indicators in determination of SQI was in the order: SOC (58 %)?>?exch. Al (17.1 %)?>?available P (8.9 %)?>?available Mn (8.2 %)?>?silt content (7.8 %). Results of this study suggest SOC and exch. Al as the two most powerful indicators of soil quality in study area. Thus, organic C and soil acidity management holds the key to improve soil quality under many exploitatively cultivated land use systems in eastern Himalayan region of India.