Influence of conduit network geometry on solute transport in karst aquifers with a permeable matrix

In karst aquifers with significant matrix permeability, water and solutes are exchanged between the conduits and carbonate matrix. Transport through the matrix increases the spread of solutes and increases travel times. This study numerically evaluates advective solute transport in synthetic karst systems that contain 3D branching conduit networks. Particle tracking is performed to analyze the spatial and temporal transport history of solute that arrives at the conduit outlet. Three measures of transport connectivity are used to quantify the solute migration behavior: the skewness of the particle arrival time distribution, the normalized fifth percentile of arrival times, and the fraction of the total travel time that occurs within conduits. All three of these metrics capture the influence of conduit network geometry on solute transport. A more tortuous network leads to enhanced conduit-matrix mixing, which reduces the transport connectivity and yields a broader distribution of solute arrival times. These results demonstrate that the conduit network geometry is an important control on solute transport in karst systems with a permeable matrix.     

Use of geochemical mass balance modelling to evaluate the role of weathering in determining stream chemistry in five mid-atlantic watersheds on different lithologies

The importance of mineral weathering was assessed and compared for five mid-Atlantic watersheds receiving similar atmospheric inputs but underlain by differing bedrock. Annual solute mass balances and volume-weighted mean solute concentrations were calculated for each watershed for each year of record. In addition, primary and secondary mineralogy were determined for each of the watersheds through analysis of soil samples and thin sections using petrographic, scanning electron microscope, electron microprobe and X-ray diffraction techniques. Mineralogical data were also compiled from the literature. These data were input to NETPATH, a geochemical program that calculates the masses of minerals that react with precipitation to produce stream water chemistry. The feasibilities of the weathering scenarios calculated by NETPATH were evaluated based on relative abundances and reactivities of minerals in the watershed. In watersheds underlain by reactive bedrocks, weathering reactions explained the stream base cation loading. In the acid-sensitive watersheds on unreactive bedrock, calculated weathering scenarios were not consistent with the abundance of reactive minerals in the underlying bedrock, and alternative sources of base cations are discussed. © 1997 John Wiley & Sons, Ltd.     

Comparison of methods for calculating annual solute exports from six forested appalachian watersheds

Six methods were compared for calculating annual stream exports of sulfate, nitrate, calcium, magnesium and aluminum from six small Appalachian watersheds. Approximately 250–400 stream samples and concurrent stream flow measurements were collected during baseflows and storm flows for the 1989 water year at five Pennsylvania watersheds and during the 1989–1992 water years at a West Virginia watershed. Continuous stream flow records were also collected at each watershed. Solute exports were calculated from the complete data set using six different scenarios ranging from instantaneous monthly measurements of stream chemistry and stream flow, to intensive monitoring of storm flow events and multiple regression equations. The results for five of the methods were compared with the regression method because statistically significant models were developed and the regression equations allowed for prediction of solute concentrations during unsampled storm flows. Results indicated that continuous stream flow measurement was critical to producing exports within 10% of regression estimates. For solutes whose concentrations were not correlated strongly with stream flow, weekly grab samples combined with continuous records of stream flow were sufficient to produce export estimates within 10% of the regression method. For solutes whose concentrations were correlated strongly with stream flow, more intensive sampling during storm flows or the use of multiple regression equations were the most appropriate methods, especially for watersheds where stream flows changed most quickly. Concentration–stream flow relationships, stream hydrological response, available resources and required level of accuracy of chemical budgets should be considered when choosing a method for calculating solute exports. © 1997 John Wiley & Sons, Ltd.     

Impacts of spatially and temporally varying snowmelt on subsurface flow in a mountainous watershed: 1. Snowmelt simulation

Abstract

The dominant source of streamflow in many mountainous watersheds is snowmelt recharge through shallow groundwater systems. The hydrological response of these watersheds is controlled by basin structure and spatially distributed snowmelt. The purpose of this series of two papers is to simulate spatially varying snowmelt and groundwater response in a small mountainous watershed. This paper examines the spatially and temporally variable snowmelt to be used as input to the groundwater flow modelling described in the second paper. Snowmelt simulation by the Simultaneous Heat and Water (SHAW) model (a detailed process model of the interrelated heat, water and solute movement through vegetative cover, snow, residue and soil) was validated by applying the model to two years of data at three sites ranging from shallow transient snow cover on a west-facing slope to a deep snow drift on a north-facing slope. The simulated energy balances for several melt periods are presented. Snow depth, density, and the magnitude and timing of snow cover outflow were simulated well for all sites.     

Base flow‐driven shifts in tropical stream temperature regimes across a mean annual rainfall gradient

Climate change is expected to affect air temperature and watershed hydrology, but the degree to which these concurrent changes affect stream temperature is not well documented in the tropics. How stream temperature varies over time under changing hydrologic conditions is difficult to isolate from seasonal changes in air temperature. Groundwater and bank storage contributions to stream flow (i.e., base flow [BF]) buffer water temperatures against seasonal and daily fluctuations in solar radiation and air temperature, whereas rainfall‐driven runoff produces flooding events that also influence stream temperature. We used a space‐for‐time substitution to examine how shifts in BF and runoff alter thermal regimes in streams by analyzing hydrological and temperature data collected from similar elevations (400–510 m above sea level) across a 3,500‐mm mean annual rainfall gradient on Hawai'i Island. Sub‐daily water temperature and stream flow gathered for 3 years were analyzed for daily, monthly, and seasonal trends and compared with air temperature measured at multiple elevations. Results indicate that decreases in median BF increased mean, maximum, and minimum water temperatures as well as daily temperature range. Monthly and daily trends in stream temperature among watersheds were more pronounced than air temperature, driven by differences in groundwater inputs and runoff. Stream temperature was strongly negatively correlated to BF during the dry season but not during the wet season due to frequent wet season runoff events contributing to total flow. In addition to projected increases in global air temperature, climate driven shifts in rainfall and runoff are likely to affect stream flow and groundwater recharge, with concurrent influences on BF resulting in shifts in water temperature that are likely to affect aquatic ecosystems.     

Evaluation of infiltration‐based stormwater management to restore hydrological processes in urban headwater streams

Urbanization threatens headwater stream ecosystems globally. Watershed restoration practices, such as infiltration‐based stormwater management, are implemented to mitigate the detrimental effects of urbanization on aquatic ecosystems. However, their effectiveness for restoring hydrologic processes and watershed storage remains poorly understood. Our study used a comparative hydrology approach to quantify the effects of urban watershed restoration on watershed hydrologic function in headwater streams within the Coastal Plain of Maryland, USA. We selected 11 headwater streams that spanned an urbanization–restoration gradient (4 forested, 4 urban‐degraded, and 3 urban‐degraded) to evaluate changes in watershed hydrologic function from both urbanization and watershed restoration. Discrete discharge and continuous, high‐frequency rainfall‐stage monitoring were conducted in each watershed. These datasets were used to develop 6 hydrologic metrics describing changes in watershed storage, flowpath connectivity, or the resultant stream flow regime. The hydrological effects of urbanization were clearly observed in all metrics, but only 1 of the 3 restored watersheds exhibited partially restored hydrologic function. At this site, a larger minimum runoff threshold was observed relative to the urban‐degraded watersheds, suggesting enhanced infiltration of stormwater runoff within the restoration structure. However, baseflow in the stream draining this watershed remained low compared to the forested reference streams, suggesting that enhanced infiltration of stormwater runoff did not recharge subsurface storage zones contributing to stream baseflow. The highly variable responses among the 3 restored watersheds were likely due to the spatial heterogeneity of urban development, including the level of impervious cover and extent of the storm sewer network. This study yielded important knowledge on how restoration strategies, such as infiltration‐based stormwater management, modulated—or failed to modulate—hydrological processes affected by urbanization, which will help improve the design of future urban watershed management strategies. More broadly, we highlighted a multimetric approach that can be used to monitor the restoration of headwater stream ecosystems in disturbed landscapes.     

Effects of forest roads on the hydrological response of a small‐scale mountain watershed in Greece

The effects of land use changes on the ecology and hydrology of natural watersheds have long been debated. However, less attention has been given to the hydrological effects of forest roads. Although less studied, several researchers have claimed that streamflow changes related to forest roads can cause a persistent and pervasive effect on hillslope hydrology and the functioning of the channel system. The main potential direct effects of forest roads on natural watersheds hydrologic response are runoff production on roads surfaces due to reduced infiltration rates, interruption of subsurface flow by road cutslopes and rapid transfer of the produced runoff to the stream network through roadside ditches. The aforementioned effects may significantly modify the total volume and timing of the hillslope flow to the stream network. This study uses detailed field data, spatial data, hydro‐meteorological records, as well as numerical simulation to investigate the effects of forest roads on the hydrological response of a small‐scale mountain experimental watershed, which is situated in the east side of Penteli Mountain, Attica, Greece. The results of this study highlight the possible effects of forest roads on the watersheds hydrological response that may significantly influence direct runoff depths and peak flow rates. It is demonstrated that these effects can be very important in permeable watersheds and that more emphasis should be given on the impact of roads on the watersheds hydrological response. Copyright © 2014 John Wiley & Sons, Ltd.     

Embedding isochronous cells overland flow module into MODFLOW

This study integrates a simple overland flow module (isochronous cells model) with the river module of MODFLOW such that temporal and spatial interactions between stream flow and groundwater can be simulated using net rainfall data of a watershed. The isochronous cells model is an efficient travel time runoff approach based on geographic information system (GIS) that considers both spatial and temporal variations of net rainfall through hill slope of the watershed. This overland module is easily coupled with MODFLOW river routing module. Specifically, the stream flow from the isochronous cells model is directly assigned to both sides of river cells of the MODFLOW model. Such an integration of MODFLOW and isochronous cells model is especially useful in watersheds where river flow data are limited. The feasibility of this integrated model was demonstrated using a case study in the middle and downstream regions of the Yitong River watershed, China. Copyright © 2012 John Wiley & Sons, Ltd.     

Seasonal and event-scale variations in solute chemistry for four Sierra Nevada catchments

Hydrobiogeochemical processes controlling stream water chemistry were examined in four small (<5 km²) catchments having contrasting bedrock lithologies in the western Sierra Nevada foothills of California. The Mediterranean climate with its cool/wet and hot/dry cycle produces strong seasonal patterns in hydrological, biological and geochemical processes. Stream water solutes fall into three general groups according to seasonal fluctuation in concentration: strong, rainy season minimum–dry season maximum (Cl⁻, SO₄²⁻, base cations); weak, rainy season minimum–dry season maximum (Si); and rainy season maximum–dry season minimum (NO₃⁻ and K⁺). Solute dynamics in soil solutions and stream water suggest that mixing of drainage waters from bedrock and soil sources regulate stream water solute concentrations. Patterns are further altered by the leaching of solutes accumulated in the soil over the summer period of desiccation and the temporal discoupling of nutrient cycles that occurs due to differences in the timing between vegetation growth (late spring) and leaching (early winter). Solute concentrations are remarkably similar between watersheds with varying bedrock types, with the exception of nitrate, sulfate and bicarbonate. Three watersheds have nitrogen-bearing metasedimentary bedrock that contributes to elevated nitrate concentrations in stream waters. Watersheds whose bedrock includes mineralized veins of sulfide and carbonate minerals similarly have greater sulfate and bicarbonate concentrations in stream water. Hydrobiogeochemical processes are highly dynamic at the seasonal and storm-event temporal scales and spatially complex at the watershed scale making management of stream water chemical composition, such as nitrate concentrations, very challenging.     

Complex patterns of catchment solute–discharge relationships for coastal plain rivers

Riverine solute versus discharge (C–Q) relationships provide information about the magnitude and dynamics of material fluxes from landscapes. We analysed long‐term patterns of C–Q relationships for 44 rivers in Florida across a suite of geogenic, nutrient, and organic solutes and investigated land cover, watershed size, and surficial geology as controls on these patterns. Solute concentrations generally exhibited far less variability than did discharge, with coherent solute‐specific behaviours repeated across watersheds. Geogenic solutes generally diluted with increasing discharge, whereas organic solutes generally enriched; patterns for nutrients were highly variable across watersheds, but on average exhibited chemostasis. Despite strong evidence of both geologic and land cover controls on solute flow‐weighted concentrations, these variables were poor predictors of C–Q slopes (β) or relative coefficients of variation (CV_C:CV_Q). CV_C:CV_Q generally increased with watershed size, and wetland area appeared to influence C–Q patterns for base cations and organic solutes. Perhaps most importantly, we observed significant slope breaks in C–Q association in approximately half of our observations, challenging the generality of using single power functions to describe catchment solute export patterns. For all solutes except phosphorus (P), C–Q slopes decreased above statistically identified breaks (slopes for P increased), with breaks consistently at or near median flow (i.e., 50% flow exceedance probability). This common pattern significantly impacts solute load estimates; failing to account for slope breaks overestimates nitrate and total organic carbon loads as much as 125% and underestimates P loads as much as 35%. In addition to challenging generic power‐law characterization of C–Q relationships for these coastal plain rivers, and exploring the load estimate consequences thereof, our study supports emerging insights about watershed hydrochemical behaviours across a wide array of solutes.     

Fine scale variations of surface water chemistry in an ephemeral to perennial drainage network

Although temporal variation in headwater stream chemistry has long been used to document baseline conditions and response to environmental drivers, less attention is paid to fine scale spatial variations that could yield clues to processes controlling stream water sources. We documented spatial and temporal variation in water composition in a headwater catchment (41 ha) at the Hubbard Brook Experimental Forest, NH, USA. We sampled every 50 m along an ephemeral to perennial stream network as well as groundwater from seeps and 35 shallow wells across varying flow conditions. Groundwater influences on surface water in this region have not been considered to be important in past studies as relatively coarse soils were assumed to be well drained in steep catchments with flashy runoff response. However, seeps displayed perennial discharge, upslope accumulated areas (UAA) smaller than those for channel initiation sites and higher pH, Ca and Si concentrations than streams, suggesting relatively long groundwater residence time or long subsurface flow paths not bound by topographic divides. Coupled with a large range in groundwater chemistry seen in wells, these results suggest stream chemistry variation reflects the range of connectivity with, and quality of, groundwater controlled by hillslope hydropedological processes. The magnitude of variations of solute concentrations seen in the first order catchment was as broad as that seen at the fifth order Hubbard Brook Valley (3519 ha). Reduction in variation in solute concentrations with increasing UAA suggested a representative elementary area (REA) value of less than 3 ha in the first order catchment, compared with 100 ha for the fifth order basin. Thus, the REA is not necessarily an elementary catchment property. Rather, the partitioning of variation between highly variable upstream sources and relatively homogenous downstream characteristics may have different physical significance depending on the scale and complexity of the catchment under examination. Copyright © 2012 John Wiley & Sons, Ltd.     

Plot and small‐watershed scale runoff from two reclaimed surface‐mined watersheds in Alberta

Infiltration excess overland flow has been identified as the dominant flow pathway in recently reclaimed surface mined watersheds as a result of compaction and sorting during the reclamation procedure. Therefore, there could be a fairly direct relationship between runoff generated from the hillslopes to that measured at the watershed outlet. A 3‐year study was initiated in 1993 to determine how well surface runoff at a watershed scale could be predicted from 1‐m² runoff frames placed on hillslopes in two reclaimed surface‐mined watersheds in central Alberta. Runoff from the hillslope frames suggests outlet discharge should be high from the 3\4‐ha Sandy Subsoil Watershed and much less for the 9\8‐ha West Watershed, but the opposite occurred. Most of the hillslope runoff from the Sandy Subsoil Watershed infiltrated once it reached the channel and depression storage played an insignificant role in determining runoff. In contrast, most of the runoff from the West Watershed originated from rain falling directly on the saturated channel (depression storage) or near‐channel saturated areas, rather than the hillslopes. Neither watershed runoff magnitude nor timing could be predicted from the same parameters for hillslope runoff frames for either reclaimed watershed. Copyright © 2000 John Wiley & Sons, Ltd.     

The physical framework of the dependence between channel flow hydrographs and drainage network morphometry

The physical basis of the linkage between magnitude and timing of channel flow hydrographs and drainage network morphometry is reviewed. Small Hortonian and structurally Hortonian networks are analysed using numerical runoff simulation. For Hortonian networks the variability of the geometry of individual channels and subcatchments within each Strahler order has generally little effect upon the overall character of the hydrograph in channels of higher order. If the network is also structurally Hortonian, the analysis of the simultaneous formation, travel, and concentration of the hydrographs in all channels of the network can be simplified to a sequence of one representative hydrograph per channel order. This approach is used in this study. Three major runoff processes control the flow hydrograph characteristics: the overland flow process which determines the water supply to the drainage network; the channel flow process which translates the hydrograph in space and time; and the drainage network process which concentrates and magnifies the flow at the junctions of the drainage network. Functional relations for the hydrograph peak, timing, and flow velocity are presented. For a given uniform rainfall and infiltration rate, the peak of the channel flow hydrograph is shown to increase geometrically with channel order, and its magnitude is directly related to the bifurcation ratio. The travel time of the peak also increases geometrically with channel order, and it is directly related to the channel length ratio over velocity ratio. The flow velocity of the peak changes in a downstream direction as a function of the bifurcation and slope ratio. It was also found that for negligible channel storage the channel flow and drainage network processes do not contribute significantly to the observed nonlinear response of a watershed to precipitation.     

不同尺度流域地表径流氮、磷浓度比较

选择太湖上游为研究对象,采集了1-400 km2不同尺度小流域产出径流TN、TP浓度实测数据,结合前期开展的地表坡面流人工暴雨实验监测结果,开展不同尺度流域水质监测对水体面源污染产出浓度估算影响的比较研究,探讨流域尺度之间入渗、汇流以及伴随的流域生态系统营养盐调节机制的差异．结果表明,流域监测尺度对土地利用面源污染产出浓度估算有较大影响．地表坡面流由于未经过流域汇流过程伴随的下渗滤过与吸附等过程,产出径流TN、TP浓度一般高于小流域．小流域林地生态系统具有较强的入渗机制、接近自然的生态沟谷汇流网络,对面源污染TN、TP有较强的削减作用．农业生态系统较弱的入渗机制、人工沟渠汇流网络对面源污染TN、TP的削减作用较弱．现代农业造成流域面源污染增加不仅仅是因为人类农业活动对流域局部土体及养分的改变,农业生态系统改变流域自然生态系统整体水文过程及营养盐调节机制也是面源污染增加的重要因素之一,恢复小尺度的生态沟谷网络系统对削减流域面源污染具有重要的意义．     

Rainfall–runoff responses on Arctic hillslopes underlain by continuous permafrost,North Slope,Alaska, USA

The Arctic hydrologic cycle is intensifying, as evidenced by increased rates of precipitation, evapotranspiration, and riverine discharge. However, the controls on water fluxes from terrestrial to aquatic systems in upland Arctic landscapes are poorly understood. Upland landscapes account for one third of the Arctic land surface and are often drained by zero‐order geomorphic flowpath features called water tracks. Previous work in the region attributed rapid runoff response at larger stream orders to water tracks, but models suggest water tracks are hydrologically disconnected from the surrounding hillslope. To better understand the role of water tracks in upland landscapes, we investigated the surface and subsurface hydrologic responses of 6 water tracks and their hillslope watersheds to natural patterns of rainfall, soil thaw, and drainage. Between storms, both water track discharge and the water table in the hillslope watersheds exhibited diel fluctuations that, when lagged by 5 hr, were temporally correlated with peak evapotranspiration rate. Water track soils remained saturated for more of the summer season than soils in their surrounding hillslope watersheds. When rainfall occurred, the subsurface response was nearly instantaneous, but the water tracks took significantly longer than the hillslopes to respond to rainfall, and longer than the responses previously observed in nearby larger order Arctic streams. There was also evidence for antecedent soil water storage conditions controlling the magnitude of runoff response. Based on these observations, we used a broken stick model to test the hypothesis that runoff production in response to individual storms was primarily controlled by rainfall amount and antecedent water storage conditions near the water track outlet. We found that the relative importance of the two factors varied by site, and that water tracks with similar watershed geometries and at similar landscape positions had similar rainfall–runoff model relationships. Thus, the response of terrestrial water fluxes in the upland Arctic to climate change depends on the non‐linear interactions between rainfall patterns and subsurface water storage capacity on hillslopes. Predicting these interactions across the landscape remains an important challenge.     

Spatial and temporal variability of stable isotopes (δ18O and δ2H) in surface waters of arid,mountainous Central Asia

Characterization of spatial and temporal variability of stable isotopes (δ¹⁸O and δ²H) of surface waters is essential to interpret hydrological processes and establish modern isotope–elevation gradients across mountainous terrains. Here, we present stable isotope data for river waters across Kyrgyzstan. River water isotopes exhibit substantial spatial heterogeneity among different watersheds in Kyrgyzstan. Higher river water isotope values were found mainly in the Issyk‐Kul Lake watershed, whereas waters in the Son‐Kul Lake watershed display lower values. Results show a close δ¹⁸O–δ²H relation between river water and the local meteoric water line, implying that river water experiences little evaporative enrichment. River water from the high‐elevation regions (e.g., Naryn and Son‐Kul Lake watershed) had the most negative isotope values, implying that river water is dominated by snowmelt. Higher deuterium excess (average d = 13.9‰) in river water probably represents the isotopic signature of combined contributions from direct precipitation and glacier melt in stream discharge across Kyrgyzstan. A significant relationship between river water δ¹⁸O and elevation was observed with a vertical lapse rate of 0.13‰/100 m. These findings provide crucial information about hydrological processes across Kyrgyzstan and contribute to a better understanding of the paleoclimate/elevation reconstruction of this region.     

Impacts of urban headwater disturbances on downstream sediment loading where streamside management zones for forest harvesting activities are present downstream

Watershed disturbances caused by land use/land cover(LU/LC) change, intensive forestry practices, and other forms of watershed disturbance are known to negatively affect water quality, if best management practices(BMPs) are not present. Disturbances near or around headwaters where BMPs are absent or are inadequate may impact downstream water quality and stream ecology regardless of the presence of downstream BMPs. In this study, the recovery of two small forest-dominated watersheds(Ww and Ew) in Auburn, Alabama, USA, was monitored 6-years after a partial clearcut was done. The streamside management zone(SMZ) in Ww was also partially harvested shortly after the partial clearcut to improve its sediment trapping effectiveness over the long-term. The hydrology and sediment transport in these watersheds were expected to stabilize over this period, however, new disturbances in the headwaters of both watersheds changed the watershed dynamics. Streamflow and sediment load were examined over a15-month period from February 2014 to May 2015. The clearcut in Ww and Ew, and the partially harvested SMZ in Ww seemed to have recovered; however, substantial amount of sediment load were measured across all monitoring sites in both watersheds as a result of upstream disturbances. This study emphasizes the necessity of watershed management as a whole, and the importance of BMP applications within the entire watershed. It also highlights the fact that watershed management and BMP implementation is a dynamic process, thus, strategies and techniques need to be continuously adapted in response to environmental changes.     

Hydrologic comparison between a forested and a wetland/lake dominated watershed using SWAT

The Soil and Water Assessment Tool (SWAT) is a physically‐based hydrologic model developed for agricultural watersheds, which has been infrequently validated for forested watersheds, particularly those with deep overwinter snow accumulation and abundant lakes and wetlands. The goal of this study was to determine the applicability of SWAT for modelling streamflow in two watersheds of the Ontonagon River basin of northern Michigan which differ in proportion of wetland and lake area. The forest‐dominated East Branch watershed contains 17% wetland and lake area, whereas the wetland/lake‐dominated Middle Branch watershed contains 26% wetland and lake area. The specific objectives were to: (1) calibrate and validate SWAT models for the East Branch and Middle Branch watersheds to simulate monthly stream flow, and (2) compare the effects of wetland and lake abundance on the magnitude and timing of streamflow. Model calibration and validation was satisfactory, as determined by deviation of discharge D and Nash and Sutcliffe coefficient values E that compared simulated monthly mean discharge versus measured monthly mean discharge. Streamflow simulation discrepancies occurred during summer and fall months and dry years. Several snow melting parameters were found to be critical for the SWAT simulation: TIMP (snow temperature lag factor) and SMFMX and SMFMN (melting factors). Snow melting parameters were not transferable between adjacent watersheds. Differences in seasonal pattern of long‐term monthly streamflow were found, with the forest‐dominated watershed having a higher peak flow during April but a lower flow during the remainder of the year in comparison to the wetland and lake‐dominated watershed. The results suggested that a greater proportion of wetland and lake area increases the capacity of a watershed to impound surface runoff and to delay storm and snow melting events. Representation of wetlands and lakes in a watershed model is required to simulate monthly stream flow in a wetland/lake‐dominated watershed. Copyright © 2007 John Wiley & Sons, Ltd.     

Changes in stream water quality due to logging of the boreal forest in the Montmorency Forest,Québec

Summer stream water quality was monitored before and following the logging of 50% of the boreal forest within three small watersheds (<50 ha) nested in the ‘Ruisseau des Eaux‐Volées’ Experimental Watershed, Montmorency Forest (Québec, Canada). Logging was conducted in winter, on snow cover according to recommended best management practices (BMPs) to minimize soil disturbance and protect advance growth. A 20‐m forest buffer was maintained along perennial streams. In watershed 7·2, cut‐blocks were located near the stream network and logging was partially allowed within the riparian buffer zone. In watersheds 7·5 and 7·7, logging occurred farther away from the stream network. Observations were also made for watershed 7·3 that collected the runoff from watersheds 7·2 and 7·5, and watershed 7·6, the uproad portion of watershed 7·7. The control watershed 0·2 was contiguous to the impacted watersheds and remained undisturbed. Following clearcutting, changes in summer daily maximum and minimum stream temperatures remained within ± 1 °C while changes in diurnal variation did not decrease by more than 0·5 °C. Concentrations of NO₃^? greatly increased by up to 6000% and concentrations of K⁺ increased by up to 300% during the second summer after logging. Smaller increases were observed for Fe_total (up to 71%), specific conductance (up to 26%), and Mg²⁺ (up to 19%). Post‐logging pH decreased slightly by no more than 7% while PO₄^3? concentration remained relatively constant. Suspended sediment concentrations appeared to increase during post‐logging, but there was not enough pre‐logging data to statistically confirm this result. Logging of moderate intensity and respecting established BMPs may account for the limited changes of water quality parameters and the low exceedances of the criteria for the protection of aquatic life. The proximity of the cutover to the stream network and logging within the riparian zone did not appear to affect water quality. Copyright © 2008 John Wiley & Sons, Ltd.     

DRAINAGE NETWORKS AFTER WILDFIRE

Predicting runoff and erosion from watersheds burned by wildfires requires an understanding of the three-dimensional structure of both hillslope and channel drainage networks. We investigate the small- and large-scale structures of drainage networks using field studies and computer analysis of 30- m digital elevation model. Topologic variables were derived from a composite 30-m DEM, which included 14 order 6 watersheds within the Pikes Peak batholith. Both topologic and hydraulic variables were measured in the field in two smaller burned watersheds （3.7 and 7.0 hectares） located within one of the order 6 watersheds burned by the 1996 Buffalo Creek Fire in Central Colorado. Horton ratios of topologic variables （stream number, drainage area, stream length, and stream slope） for small-scale and large-scale watersheds are shown to scale geometrically with stream order （i.e., to be scale invariant）. However, the ratios derived for the large-scale drainage networks could not be used to predict the rill and gully drainage network structure. Hydraulic variables （width, depth, cross- sectional area, and bed roughness） for small-scale drainage networks were found to be scale invariant across 3 to 4 stream orders. The relation between hydraulic radius and cross-sectional area is similar for fills and gullies, suggesting that their geometry can be treated similarly in hydraulic modeling. Additionally, the rills and gullies have relatively small width-to-depth ratios, implying sidewall friction may be important to the erosion and evolutionary process relative to main stem channels.