Sources of stream sulphate in headwater catchments in Otter Creek Wilderness,West Virginia,USA

Upland forested catchments in the Appalachian Plateau region receive among the greatest rates of atmospheric sulphur (S) deposition in the eastern USA, although coal mines and S‐bearing minerals in bedrock may also contribute to stream acidity in this region. Watershed mass balance and stable S isotopic values (δ³⁴S) of sulphate (SO₄^2?) were used to assess the contributions to stream SO₄^2? from atmospheric and lithogenic sources at Yellow Creek (YC), a headwater catchment on the Appalachian Plateau in West Virginia. Oxygen isotopic values (δ¹⁸O) of water were used to study catchment hydrology. Stream output of SO₄^2? was c. 60% of atmospheric S deposition during a relatively dry year, whereas atmospheric S input was nearly balanced by stream output during a year with above normal amounts of precipitation. The temporal patterns and values of δ³⁴S were similar between bulk precipitation and stream water at two upper elevation sites. At the lowest elevation site, stream δ³⁴S values were similar to bulk precipitation values during the dormant season but were slightly lower than precipitation during the low‐flow summer, probably as the result of a greater proportion of stream water being derived from deep hydrological flowpaths that have contacted S‐bearing minerals with low δ³⁴S values in coal seams. Stream δ³⁴S values at YC were significantly higher than at Coal Run, a catchment containing abandoned coal prospects and having a greater amount of S‐bearing minerals than YC. Results suggested that lithogenic S is a relatively minor source and that atmospheric deposition is the principal source of stream SO₄^2?, and thus stream acidity, at YC. Copyright © 2001 John Wiley & Sons, Ltd.     

A simple index explains annual atrazine transport from surface runoff‐prone watersheds in the north‐central USA

Year‐to‐year dynamics in weather affect both the timing of application and the potential hydrologic transport of pesticides. Further, the most commonly used pesticides dissipate in the environment during the growing season. Interactions among these factors – hydrology, timing of application and dissipation kinetics – hinder the detection of temporal trends in transport. It is increasingly important to be able to discern such trends, to judge effectiveness of management practices or to determine whether observed changes were caused by management or weather. In previous work, a cumulative vulnerability index was developed to account for these three factors. It explained 63% of annual variation in atrazine load in the Goodwater Creek Experimental Watershed (GCEW). The objectives of the current work were (i) to generalize the cumulative vulnerability index to explicitly account for variation in watershed size, area treated with atrazine and average application rate; (ii) to test the overall performance on watersheds showing such variation; and (3) to test whether the generalized index properly accounted for the additional input parameters. The generalized index was tested using data from GCEW (73.7 km²) and seven additional watersheds in the northeast Missouri claypan region that varied in size from 212 to 1180 km² and from 4% to 23% of watershed area planted to corn or sorghum. Across 32 site‐years, the generalized index explained 84% of variation in annual atrazine load. Further, tests of residuals showed no dependence on either watershed area or fraction of area planted to corn and sorghum, indicating that these parameters were properly integrated into the index. The performance of the index supports the conclusion that data obtained from GCEW is representative of the Mark Twain Lake Basin and likely the entire Central Claypan Major Land Resource Area. Published 2012. This article is a U.S. Government work and is in the public domain in the USA.     

Comparison of stormflow responses of surface‐mined and forested watersheds in the Appalachian Mountains,USA

The results of a hydrological analysis that was conducted as part of a larger, multifaceted, collaborative effort to quantify ecosystem functions in watersheds subjected to land‐use and land‐cover change are presented. The primary goal of the study was to determine whether a small watershed in the Appalachian region (USA) that was recently subjected to surface mining and reclamation practices produces stormflow responses to rain events that are different from those produced by a nearby reference watershed covered by young, second‐growth forest. Water balances indicated that runoff yields did not vary significantly between the two watersheds on an annual basis. Statistically significant differences (p?0·05) in runoff responses were observed on an event basis, however, with the mined/reclaimed watershed producing, on average (a) higher storm runoff coefficients (2·5×), (b) greater total storm runoff (3×), and (c) higher peak hourly runoff rates (2×) when compared with the reference watershed. Results of a unit hydrograph analysis also showed, unexpectedly, that the modelled unit responses of the two watersheds to effective rainfall pulses were similar, despite the noted differences in land cover. Differences in stormflow responses were thus largely explained by dramatic reductions in cumulative rates of rainfall abstraction (measured using infiltrometers) attributable to soil compaction during land reclamation. Additional field hydrological measurements on other mined watersheds will be needed to generalize our results, as well as to understand and predict the cumulative hydrological impacts of widespread surface mining in larger watersheds and river basins. Copyright © 2006 John Wiley & Sons, Ltd.     

Dynamics of stream nitrate sources and flow pathways during stormflows on urban,forest and agricultural watersheds in central Pennsylvania,USA

Understanding the influence of storm events on nitrate (NO₃^?) dynamics is important for efficiently managing NO₃^? pollution. In this study, five sites representing a downstream progression of forested uplands underlain by resistant sandstone to karst lowlands with agricultural, urban and mixed land‐use were established in Spring Creek, a 201 km² mixed land‐use watershed in central Pennsylvania, USA. At each site, stream water was monitored during six storm events in 2005 to assess changes in stable isotopes of NO₃^? (δ¹⁵N‐NO₃^? and δ¹⁸O‐NO₃^?) and water (δ¹⁸O‐H₂O) from baseflow to peakflow. Peakflow fractions of event NO₃^? and event water were then computed using two‐component mixing models to elucidate NO₃^? flow pathway differences among the five sites. For the forested upland site, storm size appeared to affect NO₃^? sources and flow pathways. During small storms (<35 mm rainfall), greater event NO₃^? fractions than event water fractions indicated the prevalence of atmospheric NO₃^? source contributions at peakflow. During larger storms (>35 mm rainfall), event NO₃^? fractions were less than event water fractions at peakflow suggesting that NO₃^? was flushed from stored sources via shallow subsurface flow pathways. For the urbanized site, wash‐off of atmospheric NO₃^? was an important NO₃^? source at peakflow, especially during short‐duration storms where event water contributions indicated the prevalence of overland flow. In the karst lowlands, very low fractions of event water and even lower fractions of event NO₃^? at peakflow suggested the dominance of ground water flow pathways during storms. These ground water flow pathways likely flushed stored NO₃^? sources into the stream, while deep soils in the karst lowlands also may have promoted NO₃^? assimilation. The results of this study illustrated how NO₃^? isotopes and δ¹⁸O‐H₂O could be combined to show key differences in water and NO₃^? delivery between forested uplands, karst valleys and fully urbanized watersheds. Copyright © 2009 John Wiley & Sons, Ltd.     

Acidification recovery in a changing climate: Observations from thirty-five years of stream chemistry monitoring in forested headwater catchments at the Turkey Lakes watershed,Ontario

Long-term ecosystem studies are valuable for understanding integrated ecosystem response to global changes in atmospheric deposition and climate. We examined trends for a 35-year period (1982/83–2017/18) in concentrations of a range of solutes in precipitation and stream water from nine headwater catchments spanning elevation and surficial geology gradients at the Turkey Lakes watershed (TLW) in northeastern Ontario, Canada. Average annual water year (WY, October to September) concentrations in precipitation significantly declined over the period for sulphate (SO₄²⁻), nitrate (NO₃⁻) and chloride (Cl⁻), while calcium (Ca²⁺) and potassium (K⁺) concentrations increased, resulting in a significant pH increase from 4.2 to 5.7. Trends in stream chemistry through time are generally consistent with expectations associated with acidification recovery. Concentration of many stream water solutes (SO₄²⁻, Cl⁻, calcium [Ca²⁺], magnesium [Mg²⁺] and NH₄⁺ generally decreased, while others (silica [SiO₂] and dissolved organic carbon [DOC]) generally increased. Increases were also observed for alkalinity (six of nine catchments), acid neutralizing capacity ([ANC]; six of nine catchments) and pH (eight of nine catchments), while conductivity declined (six of nine catchments). Variability in trends among catchments are associated with differences in surficial geology and wetland cover. While absolute solute concentrations were generally lower at bedrock dominated high-elevation catchments compared to till dominated lower elevation catchments, the rate of change of concentration was often greater for high elevation catchments. This study confirms continued, but non-linear stream chemistry recovery from acidification, particularly at the less buffered high and moderate elevation sites. The heterogeneity of responses among catchments highlights our incomplete understanding of the relative importance of different mechanisms influencing stream chemistry and the consequences for downstream ecosystems.     

Sensitivity to acidification of subalpine ponds and lakes in north‐western Colorado

Although acidifying deposition in western North America is lower than in many parts of the world, many high‐elevation ecosystems there are extremely sensitive to acidification. Previous studies determined that the Mount Zirkel Wilderness Area (MZWA) has the most acidic snowpack and aquatic ecosystems that are among the most sensitive in the region. In this study, spatial and temporal variability of ponds and lakes in and near the MZWA were examined to determine their sensitivity to acidification and the effects of acidic deposition during and after snowmelt. Within the areas identified as sensitive to acidification based on bedrock types, there was substantial variability in acid‐neutralizing capacity (ANC), which was related to differences in hydrological flowpaths that control delivery of weathering products to surface waters. Geological and topographic maps were of limited use in predicting acid sensitivity because their spatial resolution was not fine enough to capture the variability of these attributes for lakes and ponds with small catchment areas. Many of the lakes are sensitive to acidification (summer and autumn ANC < 100 µeq L^?1), but none of them appeared to be threatened immediately by episodic or chronic acidification. In contrast, 22 ponds had minimum ANC < 30 µeq L^?1, indicating that they are extremely sensitive to acidic deposition and could be damaged by episodic acidification, although net acidity (ANC < 0) was not measured in any of the ponds during the study. The lowest measured pH value was 5·4, and pH generally remained less than 6·0 throughout early summer in the most sensitive ponds, indicating that biological effects of acidification are possible at levels of atmospheric deposition that occurred during the study. The aquatic chemistry of lakes was dominated by atmospheric deposition and biogeochemical processes in soils and shallow ground water, whereas the aquatic chemistry of ponds was also affected by organic acids and biogeochemical processes in the water column and at the sediment–water interface. These results indicate that conceptual and mechanistic acidification models that have been developed for lakes and streams may be inadequate for predicting acidification in less‐understood systems such as ponds. Copyright © 2004 John Wiley & Sons, Ltd.     

The effect of terrace geology on ground‐water movement and on the interaction of ground water and surface water on a mountainside near Mirror Lake,New Hampshire,USA

The west watershed of Mirror Lake in the White Mountains of New Hampshire contains several terraces that are at different altitudes and have different geologic compositions. The lowest terrace (FSE) has 5 m of sand overlying 9 m of till. The two next successively higher terraces (FS2 and FS1) consist entirely of sand and have maximum thicknesses of about 7 m. A fourth, and highest, terrace (FS3) lies in the north‐west watershed directly adjacent to the west watershed. This highest terrace has 2 m of sand overlying 8 m of till. All terraces overlie fractured crystalline bedrock. Numerical models of hypothetical settings simulating ground‐water flow in a mountainside indicated that the presence of a terrace can cause local ground‐water flow cells to develop, and that the flow patterns differ based on the geologic composition of the terrace. For example, more ground water moves from the bedrock to the glacial deposits beneath terraces consisting completely of sand than beneath terraces that have sand underlain by till. Field data from Mirror Lake watersheds corroborate the numerical experiments. The geology of the terraces also affects how the stream draining the west watershed interacts with ground water. The stream turns part way down the mountainside and passes between the two sand terraces, essentially transecting the movement of ground water down the valley side. Transects of water‐table wells were installed across the stream's riparian zone above, between, and below the sand terraces. Head data from these wells indicated that the stream gains ground water on both sides above and below the sand terraces. However, where it flows between the sand terraces the stream gains ground water on its uphill side and loses water on its downhill side. Biogeochemical processes in the riparian zone of the flow‐through reach have resulted in anoxic ground water beneath the lower sand terrace. Results of this study indicate that it is useful to understand patterns of ground‐water flow in order to fully understand the flow and chemical characteristics of both ground water and surface water in mountainous terrain. Copyright © 2007 John Wiley & Sons, Ltd.     

Evaluating relative sensitivity of SWAT‐simulated nitrogen discharge to projected climate and land cover changes for two watersheds in North Carolina,USA

We investigated how projected changes in land cover and climate affected simulated nitrate (NO₃^?) and organic nitrogen discharge for two watersheds within the Neuse River Basin, North Carolina, USA, for years 2010–2070. We applied the Soil and Water Assessment Tool watershed model to predict nitrogen discharge using (1) atmospheric carbon dioxide (CO₂) concentrations predicted by the Intergovernmental Panel on Climate Change, (2) land cover change predicted by the Integrated Climate and Land Use Change project and (3) precipitation and temperature simulated by two statistically downscaled and bias‐corrected Global Circulation Models. We determined the sensitivity of simulated nitrogen discharge to separate changes in each treatment [(1) CO₂, (2) land cover and (3) precipitation and temperature (PT)] by comparing each treatment to a reference condition. Results showed that nitrogen discharges were most sensitive to changes in PT over the 60‐year simulation. Nitrogen discharges had similar sensitivities to the CO₂ and land cover treatments, which were only one‐tenth the influence of the PT treatment. Under the CO₂ treatment, nitrogen discharges increased with increasing ambient CO₂. NO₃^? discharge decreased with increased urbanization; however, organic nitrogen had a varied response. Under the PT treatment, there was high spatial variability in nitrogen discharges. In a single year, certain sub‐basins showed an 80% increase in nitrogen discharge relative to reference, while others showed a 400% decrease. With nitrogen discharge showing high sensitivity to PT change, we suggest that more emphasis should be placed on investigating impacts of PT on nutrient transport in the Neuse River Basin. Published 2015. This article is a U.S. Government work and is in the public domain in the USA     

Use of soil‐streamwater relationships to assess regional patterns of acidic deposition effects in the northeastern USA

Declines of acidic deposition levels by as much as 50% since 1990 have led to partial recovery of surface waters in the northeastern USA but continued depletion of soil calcium through this same period suggests a disconnection between soil and surface water chemistry. To investigate the role of soil‐surface water interactions in recovery from acidification, the first regional survey to directly relate soil chemistry to stream chemistry during high flow was implemented in a 4144‐km² area of the Catskill region of New York, where acidic deposition levels are among the highest in the East. More than 40% of 95 streams sampled in the southern Catskill Mountains were determined to be acidified and had inorganic monomeric aluminum concentrations that exceeded a threshold that is toxic to aquatic biota. More than 80% likely exceeded this threshold during the highest flows, but less than 10% of more than 100 streams sampled were acidified in the northwestern portion of the region. Median Oa horizon soil base saturation ranged from 50% to 80% at 200 sites across the region, but median base saturation in the upper 10 cm of the B horizon was less than 20% across the region and was only 2% in the southern area. Aluminum is likely to be interfering with root uptake of calcium in the mineral horizon in approximately half the sampled watersheds. Stream chemistry was highly variable over the Catskill region and, therefore, did not always reflect the calcium depletion of the B horizon that our sampling suggested was nearly ubiquitous throughout the region. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

Run‐off processes from mountains to foothills: The role of soil stratigraphy and structure in influencing run‐off characteristics across high to low relief landscapes

The critical zone features that control run‐off generation, specifically at the regional watershed scale, are not well understood. Here, we addressed this knowledge gap by quantitatively and conceptually linking regional watershed‐scale run‐off regimes with critical zone structure and climate gradients across two physiographic provinces in the Southeastern United States. We characterized long‐term (~20 years) discharge and precipitation regimes for 73 watersheds with United States Geological Survey in‐stream gaging stations across the Appalachian Mountain and Piedmont physiographic provinces of North Carolina. Watersheds included in this analysis had <10% developed land and ranged in size from 14.1–4,390 km². Thirty‐four watersheds were located in the Piedmont physiographic province, which is typically classified as a low relief landscape with deep, highly weathered soils and regolith. Thirty‐nine watersheds were located in the Appalachian Mountain physiographic province, which is typically classified as a steeper landscape with highly weathered, but shallower soils and regolith. From the United States Geological Survey daily mean run‐off time series, we calculated annual and seasonal baseflow indices (BFI), minimum, mean, and maximum daily run‐off, and Pearson's correlation coefficients between precipitation and baseflow. Our results showed that Appalachian Mountain watersheds systematically had higher minimum daily flows and BFI values. Piedmont watersheds displayed much larger deviations from mean annual BFI in response to year‐to‐year variability in precipitation. A series of linear regression models between 21 landscape metrics and annual BFIs showed non‐linear and complex terrestrial–hydrological relationships across the two provinces. From these results, we discuss how distinct features of critical zone architecture, with specific focus on soil depth and stratigraphy, may be dominating the regulation of hydrological processes and run‐off regimes across these provinces.     

Stream nitrate uptake and transient storage over a gradient of geomorphic complexity,north‐central Colorado,USA

The understanding of nutrient uptake in streams is impeded by a limited understanding of how geomorphic setting and flow regime interact with biogeochemical processing. This study investigated these interactions as they relate to transient storage and nitrate uptake in small agricultural and urban streams. Sites were selected across a gradient of channel conditions and management modifications and included three 180‐m long geomorphically distinct reaches on each of two streams in north‐central Colorado. The agricultural stream has been subject to historically variable cattle‐grazing practices, and the urban stream exhibits various levels of stabilisation and planform alteration. Reach‐scale geomorphic complexity was characterised using highly detailed surveys of channel morphology, substrate, hydraulics and habitat units. Breakthrough‐curve modelling of conservative bromide (Br^?) and nonconservative nitrate (NO₃^?) tracer injections characterised transient storage and nitrate uptake along each reach. Longitudinal roughness and flow depth were positively associated with transient storage, which was related to nitrate uptake, thus underscoring the importance of geomorphic influences on stream biogeochemical processes. In addition, changes in geomorphic characteristics due to temporal discharge variation led to complex responses in nitrate uptake. Copyright © 2011 John Wiley & Sons, Ltd.     

Mass‐balance modeling of mineral weathering rates and CO2 consumption in the forested,metabasaltic Hauver Branch watershed,Catoctin Mountain,Maryland, USA

Mineral weathering rates and a forest macronutrient uptake stoichiometry were determined for the forested, metabasaltic Hauver Branch watershed in north‐central Maryland, USA. Previous studies of Hauver Branch have had an insufficient number of analytes to permit determination of rates of all the minerals involved in chemical weathering, including biomass. More equations in the mass‐balance matrix were added using existing mineralogic information. The stoichiometry of a deciduous biomass term was determined using multi‐year weekly to biweekly stream‐water chemistry for a nearby watershed, which drains relatively unreactive quartzite bedrock. At Hauver Branch, calcite hosts ~38 mol% of the calcium ion (Ca²⁺) contained in weathering minerals, but its weathering provides ~90% of the stream water Ca²⁺. This occurs in a landscape with a regolith residence time of more than several Ka (kiloannum). Previous studies indicate that such old regolith does not typically contain dissolving calcite that affects stream Ca²⁺/Na⁺ ratios. The relatively high calcite dissolution rate likely reflects dissolution of calcite in fractures of the deep critical zone. Of the carbon dioxide (CO₂) consumed by mineral weathering, calcite is responsible for approximately 27%, with the silicate weathering consumption rate far exceeding that of the global average. The chemical weathering of mafic terrains in decaying orogens thus may be capable of influencing global geochemical cycles, and therefore, climate, on geological timescales. Based on carbon‐balance calculations, atmospheric‐derived sulfuric acid is responsible for approximately 22% of the mineral weathering occurring in the watershed. Our results suggest that rising air temperatures, driven by global warming and resulting in higher precipitation, will cause the rate of chemical weathering in the Hauver Branch watershed to increase until a threshold temperature is reached. Beyond the threshold temperature, increased recharge would produce a shallower groundwater table and reduced chemical weathering rates. Copyright © 2013 John Wiley & Sons, Ltd.     

Relative contributions of stream concentration,stream discharge and shoreline load to base cation trends in Red Chalk and Harp lakes,south‐central Ontario,Canada

We evaluated the importance and contributions of changes in stream base cation concentration, stream discharge and lake shoreline load to changes in base cation concentrations in two small inland lakes in south‐central Ontario. The shoreline load from ungauged drainage areas was calculated with a mass balance equation. An evaluation method based on the partial Mann–Kendall test quantified the relative contributions of these three explanatory drivers to the observed trends in lake concentration of Ca, K, Mg and Na. Over a 29‐year period (1978–2006) at Red Chalk and Harp lakes, declines in stream concentration and discharge were correlated with the declines in lake concentrations of all base cations (except for Na in Harp Lake), with contributions of 21–81% from stream concentration and 12–58% from discharge, whereas the shoreline load had little contribution to observed trends. The observed unusual increase in Na concentration at Harp Lake was correlated with the increase in stream Na concentration and additional load from road salts, with the shoreline load contributing a substantial 37%. These results may be applicable to numerous inland lakes with similar site conditions, including lakes found on the Canadian Shield. Copyright © 2015 John Wiley & Sons, Ltd.     

Geospatially based distributed rainfall‐runoff modelling for simulation of internal and outlet responses in a semi‐forested lower Himalayan watershed

This study presents a Geographic Information System (GIS)‐based distributed rainfall‐runoff model for simulating surface flows in small to large watersheds during isolated storm events. The model takes into account the amount of interception storage to be filled using a modified Merriam ( 1960 ) approach before estimating infiltration by the Smith and Parlange ( 1978 ) method. The mechanics of overland and channel flow are modelled by the kinematic wave approximation of the Saint Venant equations which are then numerically solved by the weighted four‐point implicit finite difference method. In this modelling the watershed was discretized into overland planes and channels using the algorithms proposed by Garbrecht and Martz ( 1999 ). The model code was first validated by comparing the model output with an analytical solution for a hypothetical plane. Then the model was tested in a medium‐sized semi‐forested watershed of Pathri Rao located in the Shivalik ranges of the Garhwal Himalayas, India. Initially, a local sensitivity analysis was performed to identify the parameters to which the model outputs like runoff volume, peak flow and time to peak flow are sensitive. Before going for model validation, calibration was performed using the Ordered‐Physics‐based Parameter Adjustment (OPPA) method. The proposed Physically Based Distributed (PBD) model was then evaluated both at the watershed outlet as well as at the internal gauging station, making this study a first of its kind in Indian watersheds. The results of performance evaluation indicate that the model has simulated the runoff hydrographs reasonably well within the watershed as well as at the watershed outlet with the same set of calibrated parameters. The model also simulates, realistically, the temporal variation of the spatial distribution of runoff over the watershed and the same has been illustrated graphically. Copyright © 2011 John Wiley & Sons, Ltd.     

Seasonal dynamics and exports of elements from a first‐order stream to a large inland lake in Michigan

Headwater streams are critical components of drainage systems, directly connecting terrestrial and downstream aquatic ecosystems. The amount of water in a stream can alter hydrologic connectivity between the stream and surrounding landscape and is ultimately an important driver of what constituents headwater streams transport. There is a shortage of studies that explore concentration–discharge (C‐Q) relationships in headwater systems, especially forested watersheds, where the hydrological and ecological processes that control the processing and export of solutes can be directly investigated. We sought to identify the temporal dynamics and spatial patterns of stream chemistry at three points along a forested headwater stream in Northern Michigan and utilize C‐Q relationships to explore transport dynamics and potential sources of solutes in the stream. Along the stream, surface flow was seasonal in the main stem, and perennial flow was spatially discontinuous for all but the lowest reaches. Spring snowmelt was the dominant hydrological event in the year with peak flows an order of magnitude larger at the mouth and upper reaches than annual mean discharge. All three C‐Q shapes (positive, negative, and flat) were observed at all locations along the stream, with a higher proportion of the analytes showing significant relationships at the mouth than at the mid or upper flumes. At the mouth, positive (flushing) C‐Q shapes were observed for dissolved organic carbon and total suspended solids, whereas negative (dilution) C‐Q shapes were observed for most cations (Na⁺, Mg²⁺, Ca²⁺) and biologically cycled anions (NO₃^?, PO₄^3?, SO₄^2?). Most analytes displayed significant C‐Q relationships at the mouth, indicating that discharge is a significant driving factor controlling stream chemistry. However, the importance of discharge appeared to decrease moving upstream to the headwaters where more localized or temporally dynamic factors may become more important controls on stream solute patterns.     

Comparison of AnnAGNPS and SWAT model simulation results in USDA‐CEAP agricultural watersheds in south‐central Kansas

This study was conducted under the USDA‐Conservation Effects Assessment Project (CEAP) in the Cheney Lake watershed in south‐central Kansas. The Cheney Lake watershed has been identified as ‘impaired waters’ under Section 303(d) of the Federal Clean Water Act for sediments and total phosphorus. The USDA‐CEAP seeks to quantify environmental benefits of conservation programmes on water quality by monitoring and modelling. Two of the most widely used USDA watershed‐scale models are Annualized AGricultural Non‐Point Source (AnnAGNPS) and Soil and Water Assessment Tool (SWAT). The objectives of this study were to compare hydrology, sediment, and total phosphorus simulation results from AnnAGNPS and SWAT in separate calibration and validation watersheds. Models were calibrated in Red Rock Creek watershed and validated in Goose Creek watershed, both sub‐watersheds of the Cheney Lake watershed. Forty‐five months (January 1997 to September 2000) of monthly measured flow and water quality data were used to evaluate the two models. Both models generally provided from fair to very good correlation and model efficiency for simulating surface runoff and sediment yield during calibration and validation (correlation coefficient; R², from 0·50 to 0·89, Nash Sutcliffe efficiency index, E, from 0·47 to 0·73, root mean square error, RMSE, from 0·25 to 0·45 m³ s^?1 for flow, from 158 to 312 Mg for sediment yield). Total phosphorus predictions from calibration and validation of SWAT indicated good correlation and model efficiency (R² from 0·60 to 0·70, E from 0·63 to 0·68) while total phosphorus predictions from validation of AnnAGNPS were from unsatisfactory to very good (R² from 0·60 to 0·77, E from ? 2·38 to 0·32). The root mean square error–observations standard deviation ratio (RSR) was estimated as excellent (from 0·08 to 0·25) for the all model simulated parameters during the calibration and validation study. The percentage bias (PBIAS) of the model simulated parameters varied from unsatisfactory to excellent (from 128 to 3). This study determined SWAT to be the most appropriate model for this watershed based on calibration and validation results. Copyright © 2008 John Wiley & Sons, Ltd.     

Effect of climate changes and water‐related human activities on annual stream flows of the Shiyang river basin in arid north‐west China

The change of annual stream flow in the Shiyang river basin, a typical arid‐inland basin in north‐west China, was investigated using hydrological, meteorological and water‐related human activities' data of the past 50 years. The long‐term trends of the hydrological time series were examined by non‐parametric techniques, including the Pettitt and Mann–Kendall tests. Double cumulative curves and multi‐regression methods were used to separate and quantify the effects of climate changes and human activities on the stream flows. The results show that the study area has been experiencing a significant upward warming trend since 1986 and precipitation shows a decreasing trend in the mountainous region but an increasing trend in the plains region. All stream flows in the upper reach and lower reaches of the Shiyang river exhibit decreasing tendencies. Since 1970, human activities, such as irrigation, have had a significant effect on the upstream flow, and account for 60% of total flow decreases in the 1970s. However, climate changes are the main reason for the observed flow decreases in the 1980s and 1990s, with contributions to total flow decrease of 68% and 63%, respectively. Before 1975, flow decreases in the upper reaches were the main factor causing reduced flows in the lower reaches of the Shiyang river. After 1975, the effect of human activities became more pronounced, with contributions of 63%, 68% and 56% to total flow decreases in the lower reaches of the Shiyang river in the periods 1975 to 1980, 1980s and 1990s, respectively. As a result, climate change is responsible for a large proportion of the flow decreases in the upstream section of the catchment during the 1980s and 1990s, while human activities have caused flow decreases downstream during the same period. Copyright © 2007 John Wiley & Sons, Ltd.     

Oxygen‐18 dynamics in precipitation and streamflow in a semi‐arid agricultural watershed,Eastern Washington,USA

Understanding flow pathways and mechanisms that generate streamflow is important to understanding agrochemical contamination in surface waters in agricultural watersheds. Two environmental tracers, δ¹⁸O and electrical conductivity (EC), were monitored in tile drainage (draining 12 ha) and stream water (draining nested catchments of 6‐5700 ha) from 2000 to 2008 in the semi‐arid agricultural Missouri Flat Creek (MFC) watershed, near Pullman Washington, USA. Tile drainage and streamflow generated in the watershed were found to have baseline δ¹⁸O value of ?14·7‰ (VSMOW) year round. Winter precipitation accounted for 67% of total annual precipitation and was found to dominate streamflow, tile drainage, and groundwater recharge. ‘Old’ and ‘new’ water partitioning in streamflow were not identifiable using δ¹⁸O, but seasonal shifts of nitrate‐corrected EC suggest that deep soil pathways primarily generated summer streamflow (mean EC 250 µS/cm) while shallow soil pathways dominated streamflow generation during winter (EC declining as low as 100 µS/cm). Using summer isotopic and EC excursions from tile drainage in larger catchment (4700‐5700 ha) stream waters, summer in‐stream evaporation fractions were estimated to be from 20% to 40%, with the greatest evaporation occurring from August to October. Seasonal watershed and environmental tracer dynamics in the MFC watershed appeared to be similar to those at larger watershed scales in the Palouse River basin. A 0·9‰ enrichment, in shallow groundwater drained to streams (tile drainage and soil seepage), of δ¹⁸O values from 2000 to 2008 may be evidence of altered precipitation conditions due to the Pacific Decadal Oscillation (PDO) in the Inland Northwest. Copyright © 2009 John Wiley & Sons, Ltd.     

Application of the WEPP model to determine sources of run‐off and sediment in a forested watershed

This study investigates critical run‐off and sediment production sources in a forested Kasilian watershed located in northern Iran. The Water Erosion Prediction Project (WEPP) watershed model was set up to simulate the run‐off and sediment yields. WEPP was calibrated and validated against measured rainfall–run‐off–sediment data. Results showed that simulated run‐off and sediment yields of the watershed were in agreement with the measured data for the calibration and validation periods. While low and medium values of run‐off and sediment yields were adequately simulated by the WEPP model, high run‐off and sediment yield values were underestimated. Performance of the model was evaluated as very good and satisfactory during the calibration and validation stages, respectively. Total soil erosion and sediment load of the study watershed during the study period were determined to be 10 108 t yr^?1 and 8735 t yr^?1, respectively. The northern areas of the watershed with dry farming were identified as the critical erosion prone zones. To prioritize the subwatersheds based on their contribution to the run‐off and sediment production at the watershed's main outlet, unit response approach (URA) was applied. In this regard, subwatersheds close to the main outlet were found to have the highest contribution to sediment yield of the whole watershed. Results indicated that depending on the objective of land and water conservation practices, particularly, for controlling sediment yield at the main outlet, critical areas for implementing the best management practices may be identified through conjunctive application of WEPP and URA. Copyright © 2014 John Wiley & Sons, Ltd.     

Lithology‐controlled evolution of stream bed sediment and basin‐scale sediment yields in adjacent mountain watersheds,Idaho, USA

The composition, grain‐size, and flux of stream sediment evolve downstream in response to variations in basin‐scale sediment delivery, channel network structure, and diminution during transport. Here, we document downstream changes in lithology and grain size within two adjacent ~300 km² catchments in the northern Rocky Mountains, USA, which drain differing mixtures of soft and resistant rock types, and where measured sediment yields differ two‐fold. We use a simple erosion–abrasion mass balance model to predict the downstream evolution of sediment flux and composition using a Monte Carlo approach constrained by measured sediment flux. Results show that the downstream evolution of the bed sediment composition is predictably related to changes in underlying geology, influencing the proportion of sediment carried as bedload or suspended load. In the Big Wood basin, particle abrasion reduces the proportion of fine‐grained sedimentary and volcanic rocks, depressing bedload in favor of suspended load. Reduced bedload transport leads to stronger bed armoring, and coarse granitic rocks are concentrated in the stream bed. By contrast, in the North Fork Big Lost basin, bedload yields are three times higher, the stream bed is less armored, and bed sediment becomes dominated by durable quartzitic sandstones. For both basins, the geology‐based mass balance model can reproduce within ~5% root‐mean‐square error the composition of the bed substrate using realistic erosion and abrasion parameters. As bed sediment evolves downstream, bedload fluxes increase and decrease as a function of the abrasion parameter and the frequency and size of tributary junctions, while suspended load increases steadily. Variable erosion and abrasion rates produce conditions of variable bed‐material transport rates that are sensitive to the distribution of lithologies and channel network structure, and, provided sufficient diversity in bedrock geology, measurements of bed sediment composition allow for an assessment of sediment source areas and yield using a simple modeling approach. Copyright © 2016 John Wiley & Sons, Ltd.