Water Quality Assessment and Modeling of an Effluent-Dominated Stream, the Notwane River, Botswana

In an effort to assess current and future water quality of the only perennial river in southeastern Botswana, this study presents water quality monitoring and modeling results for the effluent-dependent Notwane River. The water quality along the Notwane River, pre- and post-implementation of secondary wastewater treatment, was compared and results demonstrated that water quality improved after the new wastewater treatment plant (WWTP) went online. However, stream standards for chemical oxygen demand, total dissolved phosphorous, and fecal coliform were exceeded in most locations and the critical dissolved oxygen (DO) reached concentrations of less than 4 mg L⁻¹. High dissolved P concentrations and intense macrophyte growth at the impounding ponds and at sites within 30 km of the effluent waste stream confluence suggest that eutrophication was a function of P release from the ponds. Results of DO modeling demonstrated that an unpolluted inflow at approximately 10 km downstream of the confluence was responsible for raising DO concentrations by 2.3 mg L⁻¹, while SOD was responsible for a decline in DO concentrations of 1.4 mg L⁻¹ at 6 km downstream of the confluence. Simulations also showed higher DO concentrations during winter months, when water temperatures were lower. Simulations, in which the distributed biochemical oxygen demand (BOD) loading from cattle excrement was decreased, produced nominal increases in DO concentrations. An increase in WWTP BOD loadings to projected 2020 values resulted in a 1.3 mg L⁻¹ decrease in the critical DO concentration. Furthermore, a decrease in treatment plant efficiency, from 94% to 70% BOD removal, produced critical DO concentrations and anoxia along much of the modeled reach. This has significant implications for Gaborone, especially if decreased WWTP efficiency occurs as a result of the expected future increase in pollutant loadings.     

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.     

Biogeochemistry of the Penobscot River watershed, Maine, USA: nutrient export patterns for carbon, nitrogen, and phosphorus

Watershed exports of carbon, nitrogen, phosphorus, major solutes, and suspended sediments were examined during five water years in the Penobscot River basin, which forms part of the Gulf of Maine watershed. Mean annual exports of dissolved organic carbon (DOC) in the Penobscot River were 58 kg C ha(-1) year(-1), whereas cumulative yearly watershed flux of DOC during the study period ranged from 8.6 to 16.1 × 10(10) g C year(-1) and averaged 11.7 × 10(10) g C year(-1). Watershed exports of total soluble N (TN) and total soluble P in the Penobscot River averaged 1.9 and 0.02 kg ha(-1) year(-1), respectively. Companion studies in two other major Maine rivers indicated that mean annual exports of DOC and TN in the Androscoggin River were 40 kg C ha(-1) year(-1) and 2.0 kg N ha(-1) year(-1), whereas exports in the Kennebec River were 43 kg C ha(-1) year(-1) and 2.2 kg N ha(-1) year(-1). Extrapolation of results from this investigation and a previous complementary study indicates that estuaries and coastal waters in the Gulf of Maine receive at least 1.0 × 10(10) g N year(-1) and 2.5 × 10(11) g C year(-1) in combined runoff from the four largest Maine river basins. Soluble exports of Ca + Mg + Na minus wet deposition inputs of cations in the Penobscot system were approximately 1,840 mol(c) ha(-1) year(-1), which represents a minimum estimate of cation denudation from the watershed. Based on its low N and P export rates, the Penobscot River watershed represents an example of reference conditions for use as a benchmark in ecological assessments of river water quality restoration or impairment. In addition, the biogeochemical metrics from this study provide an historical baseline for analysis of future trends in nutrient exports from the Penobscot watershed as a function of changing climatic and land use patterns.     

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.     

Spatial and temporal variations of nitrogen pollution in Wen-Rui Tang River watershed, Zhejiang, China

Water quality has degraded dramatically in Wen-Rui Tang River watershed, Zhejiang, China, especially due to rapid economic development since 1995. This paper aims to assess spatial and temporal variations of the main pollutants (NH??-N, TN, BOD(5), COD(Mn), DO) of water quality in Wen-Rui Tang River watershed, using the geographic information system, cluster analysis (CA) and principal component analysis (PCA). Results showed that concentrations of BOD(5), COD(Mn), NH??-N, and TN were significantly higher in tertiary rivers than in primary and secondary rivers. From April 2006 to March 2007, the concentrations of NH? ?-N (2.25-57.9 mg/L) and TN (3.78-70.4 mg/L) in all samples exceeded Type V national water quality standards (≥2 mg/L), while 5.3% of all COD(Mn) (1.83-27.5 mg/L) and 33.6% of all BOD(5) (0.34-50.4 mg/L) samples exceeded Type V national water quality standards (COD(Mn)?≥ 15 mg/L, BOD(5)?≥ 10 mg/L). Monthly changes of pollutant concentrations did not show a clear pattern, but correlation analysis indicated that NH??-N and TN in tertiary rivers had a significant negative correlation with 5-day cumulative rainfall and monthly rainfall, while there were no significant correlations in primary and secondary rivers. The results of CA and spatial analysis showed that the northern part of Wen-Rui Tang River watershed was the most seriously polluted. This region is characterized by the high population density and industrial and commercial activities. The PCA and spatial analysis indicated that the degraded water quality is caused by anthropogenic activities and poor wastewater management.     

Field assessment of flood event suspended sediment transport from ephemeral streams in the tropical semi-arid catchments

An assessment of suspended sediment transport was carried out in a number of semiarid catchments during flood events in order to quantify the degradation rates. In order to quantify these, a systematic sampling procedure of the episodic flood events was proposed for representative catchments. The procedure allows for an integration over the whole run-off episode using both the rising and falling limbs of the run-off hydrograph to compute the sediment quantities for each individual flood event. Higher sediment concentrations occurred in the rising limb than those at the recession for any stage of flow. The maximum suspended sediment concentration was observed at the peak of the flood hydrograph. An integration of the sediment concentration over its duration gave the total sediment yield from the flood event. For the ephemeral channels, only a small number of flood events were observed over a three-year experimental period each with a duration of the order of 3–6 h. It is notable that high sediment loads were associated with high flow volumes which were effectively the result of the catchment characteristics and incident rainfall causing the flood events in the respective catchments. A large percentage of the annual sediment yield from a catchment is transported by the ephemeral streams during a small number of flood events. The correct determination of the total sediment yield from any of the flood events depends entirely on the accuracy of the measurements. The understanding of run-off and sediment loss for the representative catchments aims at assisting planning, management and control of water and land resources for sustainable development in the semi-arid parts of the tropics. The sediment rates reveal the degradation of catchments which have repercussions on the crop and pasture production and this has a bearing on the soil and water conservation programmes in the delicate ecological balance of the semi-arid areas. Further, these rates will determine the lifespan of the reservoirs planned for the dry river valleys (ephemeral streams) and existing ones which serve livestock and domestic needs. These occasionally will require costly rehabilitation and scooping to increase effective storage unless conservation measures are taken, and these measures are bound to vary from place to place as per the representative catchments output.     

Environmental change monitoring in the arid and semi-arid regions: a case study Al-Basrah Province, Iraq

In recent years, land use/cover dynamic change has become a key subject urgently to be dealt with in the study of global environmental change. This research utilizes the integrated remote sensing and geographic information systems (GIS) in the southern part of Iraq (Basrah Province was taken as a case) to monitor, map, and quantify the environmental change using a 1:250,000 mapping scale. Remote sensing and GIS software were used to classify Landsat TM in 1990 and Landsat ETM+ in 2003 imagery into five land use and land cover (LULC) classes: vegetation land, sand land, urban area, unused land, and water bodies. Supervised classification and normalized difference buildup index, normalized difference vegetation index, normalized difference bare land index, the normalized differential water index, crust index (CI) algorithms, and change detection techniques were adopted in this research and used, respectively, to retrieve its class boundary. An accuracy assessment was performed on the 2003 LULC map to determine the reliability of the map. Finally, GIS software was used to quantify and illustrate the various LULC conversions that took place over the 13-year span of time. The results showed that the urban area, sand lands, and bare lands had increased by the rate of 1.2%, 0.8%, and 0.4% per year, with area expansion from 3,299.1, 4,119.1 km², and 3,201.9 km² in 1990 to 3,794.9, 4,557.7, and 3,351.7 km² in 2003, respectively. While the vegetation cover and water body classes were about 43.5% in 1990, the percentage decreased to about 39.6% in 2003. This study demonstrates the effectiveness of the remote sensing and GIS technologies in detecting, assessing, mapping, and monitoring the environmental changes.     

Concentrations and inventories of polycyclic aromatic hydrocarbons and organochlorine pesticides in watershed soils in the Pearl River Delta,China

The concentration levels, source, and inventories of polycyclic aromatic hydrocarbons (PAHs) and organochlorine pesticides (OCPs) in 55 surface vegetable soils in the watershed of the Pearl River Delta (PRD) were analyzed and compared with those of the surface sediments in the Pearl River Estuary (PRE) and northern South China Sea (SCS). The 16 priority PAHs on US EPA list range from 58 to 3,077 microg/kg (average: 315 microg/kg). The concentrations of DDTs and HCHs range from 3.58 to 831 microg/kg (average: 82.1 microg/kg) and from 0.19 to 42.3 microg/kg (average: 4.42 microg/kg). The ratios of DDT/ (DDD+DDE) are higher than 2 in majority of the soil samples, suggesting that DDT contamination still exists. The PAH ratios suggest that the source of PAHs is petroleum, and combustion of fossil fuel, biomass, and coal. The average concentrations of PAHs and the linear regression slope between PAHs and TOC for the soils and the sediments are quite similar. It was estimated that the soil mass inventories at 0-20 cm depth are 1,292 metric tons for PAHs and 356 metric tons for OCPs in the studied region. The average PAHs inventory per unit area for the soil samples investigated in PRD is about 0.86 time that of surface sediments in the Pearl River Estuary, and about 2.43 times that of surface sediments in the northern South China Sea. PAHs in the soils in PRD have similar source to those of the surface sediments in PRE. All of those may suggest that PAHs in PRE and SCS are probably mainly inputted from the soils in PRD via soil erosion and river transport.     

A GIS-based spatial multi-index model for flood risk assessment in the Yangtze River Basin,China

This paper developed a GIS-based spatial multi-index model for large basin-scale flood risk assessment. In terms of the risk definition proposed by the IPCC, the flood risk in the Yangtze River Basin (YRB) was classified into indexes of hazard, vulnerability, and exposure. The model systematically accounts for various flood risk indicators related to the economic, social and ecological environment of the YRB. Using the robust data space analysis and processing capabilities of ArcGIS, these flood risk indicators were superimposed and analyzed to generate an integrated flood risk spatial distribution map for the YRB. The modeling results were verified reasonably well using observation data from the YRB floods in 1998, 2008, and 2016. We found that 24.90% of the study area was at very high and high risk in 1998, and the risk in these areas decreased to 15.95% and 17.61% in 2008 and 2016, respectively. We believe that the GIS-based spatial multi-index model can be applied to other areas where basin-scale flood risk assessment is desired and contribute to further scientific research on flood forecasting and mitigation.     

Soil erosion and sediment fluxes analysis: a watershed study of the Ni Reservoir,Spotsylvania County,VA, USA

Anthropogenic forces that alter the physical landscape are known to cause significant soil erosion, which has negative impact on surface water bodies, such as rivers, lakes/reservoirs, and coastal zones, and thus sediment control has become one of the central aspects of catchment management planning. The revised universal soil loss equation empirical model, erosion pins, and isotopic sediment core analyses were used to evaluate watershed erosion, stream bank erosion, and reservoir sediment accumulation rates for Ni Reservoir, in central Virginia. Land-use and land cover seems to be dominant control in watershed soil erosion, with barren land and human-disturbed areas contributing the most sediment, and forest and herbaceous areas contributing the least. Results show a 7 % increase in human development from 2001 (14 %) to 2009 (21.6 %), corresponding to an increase in soil loss of 0.82 Mg ha^-1 year^-1 in the same time period. ²¹⁰Pb-based sediment accumulation rates at three locations in Ni Reservoir were 1.020, 0.364, and 0.543 g cm^-2 year^-1 respectively, indicating that sediment accumulation and distribution in the reservoir is influenced by reservoir configuration and significant contributions from bedload. All three locations indicate an increase in modern sediment accumulation rates. Erosion pin results show variability in stream bank erosion with values ranging from 4.7 to 11.3 cm year^-1. These results indicate that urban growth and the decline in vegetative cover has increased sediment fluxes from the watershed and poses a significant threat to the long-term sustainability of the Ni Reservoir as urbanization continues to increase.     

Applicability of hydromorphological monitoring data to locate flood risk reduction measures: Blanice River basin, Czech Republic

This paper presents a new methodological approach, using hydromorphological monitoring data from river networks, to identify stream elements that may negatively influence the course of floods, as well as elements suitable for modifications that will reinforce the natural retention potential of the floodplain. The approach uses data derived from hydromorphological monitoring that conforms to the EN 14614 standard "Water quality--Guidance standard for assessing the hydromorphological features of rivers." Four major aspects of stream modifications that can potentially affect the course and consequences of floods, and that can be derived from hydromorphological monitoring, were identified. These are: (1) modifications that increase the speed of runoff, (2) modifications that restrict the retention potential of the floodplain, (3) potential obstacles to the flood course, and (4) improper alterations of stream route modifications. A GIS is used to identify and locate these critical aspects of stream and floodplain modification. This method was applied in the Blanice River basin, an area in central Europe hit by a severe flood in August 2002. The critical modifications identified were compared with geomorphological mapping to display the different effects of individual types of modifications.     

Spatial heterogeneity of stream environmental conditions and macroinvertebrates community in an agriculture dominated watershed and management implications for a large river (the Liao River,China) basin

Understanding the effects of watershed land uses (e.g., agriculture, urban industry) on stream ecological conditions is important for the management of large river basins. A total of 41 and 56 stream sites (from first to fourth order) that were under a gradient of watershed land uses were monitored in 2009 and 2010, respectively, in the Liao River Basin, Northeast China. The monitoring results showed that a total of 192 taxa belonging to four phyla, seven classes, 21 orders and 91 families were identified. The composition of macroinvertebrate community in the Liao River Basin was dominated by aquatic insect taxa (Ephemeroptera and Diptera), Oligochaeta and Molluscs. The functional feeding group GC (Gatherer/Collector) was dominant in the whole basin. Statistical results showed that sites with less watershed impacts (lower order sites) were characterized by higher current velocity and habitat score, more sensitive taxa (e.g., Ephemeroptera), and the substrate was dominated by high percentage of cobble and pebble. The sites with more impacts from agriculture and urban industry (higher order sites) were characterized by higher biochemical (BOD₅) and chemical oxygen demand (COD), more tolerant taxa (e.g., Chironominae), and the substrate was dominated by silt and sand. Agriculture and urban-industry activities have reduced habitat condition, increased organic pollutants, reduced macroinvertebrate abundance, diversity, and sensitive taxa in streams of the lower Liao River Basin. Restoration of degraded habitat condition and control of watershed organic pollutants could be potential management priorities for the Basin.     

Trace metals solubility in rainwater: evaluation of rainwater quality at a watershed area, Istanbul

In this study, 79 bulk precipitation samples were collected at two sampling sites near Büyükçekmece Lake, one of the important drinking water sources of Istanbul, for the period of October 2001 to July 2002. The study comprised the determination of trace and toxic metals concentrations in rain water. The concentrations of the metals in this study were found to be higher than those reported by other researchers around the world. The solubility of toxic metals was found in the order of Cd > Cu > V > Zn > Ni > Pb > Cr. Solubility of metals under acidic conditions (pH < 5.5) was approximately five times higher than those under neutral conditions with Cd as the most soluble metal (50% soluble). Statistical evaluations including seasonal variations, crustal enrichment factors, and correlation matrix were discussed to identify the possible sources of these pollutants. The study revealed that anthropogenic elements were highly enriched especially for Cd > Cu > Pb which were found to be highly enriched. Significant portion of Cu and Pb could be increased by the effect of local sources like cement industry in the area; however, the rest of the investigated trace metals could be brought to the sampling site by long-range transport to the Büyükçekmece Lake watershed area.     

Land-cover changes in an urban lake watershed in a mega-city, Central China

Urbanization can exert a profound influence on land covers and landscape characteristics. In this study, we characterize the impact of urbanization on land cover and lacustrine landscape and their consequences in a large urban lake watershed, Donghu Lake watershed (the largest urban lake in China), Central China, by using Landsat TM satellite images of three periods of 1987, 1993 and 1999 and ground-based information. We grouped the land covers into six categories: water body, vegetable land, forested land, shrub-grass land, open area and urban land, and calculated patch-related landscape indices to analyze the effects of urbanization on landscape features. We overlaid the land cover maps of the three periods to track the land cover change processes. The results indicated that urban land continuously expanded from 9.1% of the total watershed area in 1987, to 19.4% in 1993, and to 29.6% in 1999. The vegetable land increased from 7.0% in 1987, 11.9% in 1993, to 13.9% in 1999 to sustain the demands of vegetable for increased urban population. Concurrently, continuous reduction of other land cover types occurred between 1987 and 1999: water body decreased from 30.4% to 23.8%, and forested land from 33.6% to 24.3%. We found that the expansion of urban land has at least in part caused a decrease in relatively wild habitats, such as urban forest and lake water area. These alterations had resulted in significant negative environmental consequences, including decline of lakes, deterioration of water and air quality, and loss of biodiversity.     

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.     

Trace-metal contamination in the glacierized Rio Santa watershed,Peru

The objective of this research is to characterize the variability of trace metals in the Rio Santa watershed based on synoptic sampling applied at a large scale. To that end, we propose a combination of methods based on the collection of water, suspended sediments, and riverbed sediments at different points of the watershed within a very limited period. Forty points within the Rio Santa watershed were sampled between June 21 and July 8, 2013. Forty water samples, 36 suspended sediments, and 34 riverbed sediments were analyzed for seven trace metals. The results, which were normalized using the USEPA guideline for water and sediments, show that the Rio Santa water exhibits Mn concentrations higher than the guideline at more than 50% of the sampling points. As is the second highest contaminating element in the water, with approximately 10% of the samples containing concentrations above the guideline. Sediments collected in the Rio Santa riverbed were heavily contaminated by at least four of the tested elements at nearly 85% of the sample points, with As presenting the highest normalized concentration, at more than ten times the guideline. As, Cd, Fe, Pb, and Zn present similar concentration trends in the sediment all along the Rio Santa.The findings indicate that care should be taken in using the Rio Santa water and sediments for purposes that could affect the health of humans or the ecosystem. The situation is worse in some tributaries in the southern part of the watershed that host both active and abandoned mines and ore-processing plants.     

Urban stream syndrome in a small,lightly developed watershed: a statistical analysis of water chemistry parameters,land use patterns,and natural sources

The relationships among land use patterns, geology, soil, and major solute concentrations in stream water for eight tributaries of the Kayaderosseras Creek watershed in Saratoga County, NY, were investigated using Pearson correlation coefficients and multivariate regression analysis. Sub-watersheds corresponding to each sampling site were delineated, and land use patterns were determined for each of the eight sub-watersheds using GIS. Four land use categories (urban development, agriculture, forests, and wetlands) constituted more than 99 % of the land in the sub-watersheds. Eleven water chemistry parameters were highly and positively correlated with each other and urban development. Multivariate regression models indicated urban development was the most powerful predictor for the same eleven parameters (conductivity, TN, TP, NO $_{3}^-$ , Cl^?, HCO $_{3}^-$ , SO $_{4}^{2-}$ , Na⁺, K⁺, Ca²⁺, and Mg²⁺). Adjusted R ² values, ranging from 19 to 91 %, indicated that these models explained an average of 64 % of the variance in these 11 parameters across the samples and 70 % when Mg²⁺ was omitted. The more common R ², ranging from 29 to 92 %, averaged 68 % for these 11 parameters and 72 % when Mg²⁺ was omitted. Water quality improved most with forest coverage in stream watersheds. The strong associations between water quality variables and urban development indicated an urban source for these 11 water quality parameters at all eight sampling sites was likely, suggesting that urban stream syndrome can be detected even on a relatively small scale in a lightly developed area. Possible urban sources of Ca²⁺ and HCO $_{3}^-$ are suggested.     

Adapting the RUSLE and GIS to model soil erosion risk in a mountains karst watershed, Guizhou Province, China

Soil erosion is a serious environmental problem in Guizhou Province, which is located in the centre of the karst areas of southwestern China. Unfortunately, Guizhou Province suffers from a lack of financial resources to research, monitor and model soil erosion at large watershed. In order to assess the soil erosion risk, soil erosion modeling at the watershed scale are urgently needed to be undertaken. This study integrated the Revised Universal Soil Loss Equation (RUSLE) with a Geographic Information System (GIS) to estimate soil loss and identify the risk erosion areas in the Maotiao River watershed, which is a typical rural watershed in Guizhou Province. All factors used in the RUSLE were calculated for the watershed using local data. It was classified into five categories ranging from minimal risk to extreme erosion risk depending on the calculated soil erosion amount. The soil erosion map was linked to land use, elevation and slope maps to explore the relationship between soil erosion and environmental factors and identify the areas of soil erosion risk. The results can be used to advice the local government in prioritizing the areas of immediate erosion mitigation. The integrated approach allows for relatively easy, fast, and cost-effective estimation of spatially distributed soil erosion. It thus indicates that RUSLE-GIS model is a useful and efficient tool for evaluating and mapping soil erosion risk at a large watershed scale in Guizhou Province.