Perennial frost blisters of the outer Mackenzie Delta,western Arctic coast,Canada

Saturated floodplains in Arctic deltas provide conditions favourable for frost mound growth. Little work has been reported from these settings to determine the origin of frost mounds and the controls on their distribution, to assess the longevity of individual mounds, or to quantify variation of mound distribution over time. A case study is presented on low mounds in low‐centred syngenetic ice‐wedge polygons of Big Lake Delta Plain, outer Mackenzie Delta. In 2008 and 2009, 12 mounds were examined by drilling to describe their morphologic variations and to investigate their growth processes. The mounds, containing a core of ice 15 to 58 cm thick, were less than 1 m high and 3 · 7 to 8 · 5 m in diameter; other mounds were over 10 m long. Organic inclusions in the ice, bubble densities, electrical conductivity profiles, and ice‐crystal structure indicated that the mounds were hydrostatic frost blisters. Up to six frost blisters were found within individual polygons due to the relatively small volume of water needed to create each mound. Frost‐blister densities, of greater than 1700 km^–2, increased toward the wet centres of alluvial islands down gentle topographic gradients. The frost blisters were perennial, with individuals remaining identifiable on aerial photographs and satellite images for up to 10 years. Frost blisters collapsed along dilation cracks opened by hydrostatic uplift and by thawing from their sides caused by snow drifting and water ponding. Cyclical growth and decay of the mounds may degrade the visible polygonal network over time. Copyright © 2013 John Wiley & Sons, Ltd.     

The significance of channel recharge rates for estimating debris‐flow magnitude and frequency

The quantification of debris‐flow hazard requires estimates of debris‐flow frequency and magnitude. Several methods have been proposed to determine the probable volume of future debris flows from a given basin, but most have neglected to account for debris recharge rates over time, which may lead to underestimation of debris‐flow volumes in basins with rare debris flows. This paper deals with the determination of debris recharge rates in debris‐flow channels based on knowledge of debris storage and the elapsed time since the last debris flow. Data are obtained from coastal British Columbia and a relation is obtained across a sample of basins with similar terrain and climatic conditions. For Rennell Sound on the west coast of the Queen Charlotte Islands, the power‐law relation for area‐normalized recharge rate, R_t, versus elapsed time, t_e was R_t = 0·23t_e^?0·58 with an explained variance of 75 per cent. A difference in recharge rates may exist between creeks in logged and unlogged forested terrain. The power function for undisturbed terrain was R_t = 0·20t_e^?0·49, while the function for logged areas was R_t = 0·30t_e^?0·77. This result suggests that for the same elapsed time since the last debris flow, clearcut gullies tend to recharge at a slower rate than creeks in old growth forest. This finding requires verification, particularly for longer elapsed times since debris flow, but would have important implications for forest resource management in steep coastal terrain. This study demonstrates that commonly used encounter probability equations are inappropriate for recharge‐limited debris flow channels. Copyright © 2005 John Wiley & Sons, Ltd.     

Synergistic effect of vegetation and air temperature changes on soil water content in alpine frost meadow soil in the permafrost region of Qinghai‐Tibet

Seasonal changes over 2 years (2004–2006) in soil moisture content (θ_v) of frozen alpine frost meadow soils of the Qinghai‐Tibet plateau permafrost region under three different levels of vegetation cover were investigated. Vegetation cover and air temperature changes had significant effects (synergistic effect) on θ_v and its distribution in the soil profile. During periods of soil freezing or thawing, the less the vegetation cover, the quicker the temperature drop or rise of soil water, and the shorter the duration of the soil water freeze–thaw response in the active soil layer. Under 30% and 65% vegetation cover the amplitude of variation in θ_v during the freezing period was 20–26% greater than that under 93% cover, while during the thawing period, it was 1·5‐ to 40·5‐fold greater. The freezing temperature of the surface soil layer, ^fT_s, was 1·6 °C lower under 30% vegetation cover than under 93% vegetation cover. Changes in vegetation cover of the alpine frost meadow affected θ_v and its distribution, as well as the relationship between θ_v and soil temperature (T_s). As vegetation cover decreased, soil water circulation in the active layer increased, and the response to temperature of the water distribution across the soil profile was heightened. The quantity of transitional soil phase water at different depths significantly increased as vegetation cover decreased. The influence of vegetation cover and soil temperature distribution led to a relatively dry soil layer in the middle of the profile (0·70–0·80 m) under high vegetation cover. Alpine meadow θ_v and its pattern of distribution in the permafrost region were the result of the synergistic effect of air temperature and vegetation cover. Copyright © 2007 John Wiley & Sons, Ltd.     

Heat tracing of heterogeneous hyporheic exchange adjacent to in‐stream geomorphic features

There are many field techniques used to quantify rates of hyporheic exchange, which can vary in magnitude and direction spatially over distances of only a few metres, both within and between morphological features. We used in‐stream mini‐piezometers and heat transport modelling of stream and streambed temperatures to quantify the rates and directions of water flux across the streambed interface upstream and downstream of three types of in‐stream geomorphic features: a permanent dam, a beaver dam remnant and a stream meander. We derived hyporheic flux estimates at three different depths at six different sites for a month and then paired those flux rates with measurements of gradient to derive hydraulic conductivity (K) of the streambed sediments. Heat transport modelling provided consistent daily flux estimates that were in agreement directionally with hydraulic gradient measurements and also identified vertical heterogeneities in hydraulic conductivity that led to variable hyporheic exchange. Streambed K varied over an order of magnitude (1·9 × 10^?6 to 5·7 × 10^?5 m/s). Average rates of hyporheic flux ranged from static (q < ±0·02 m/day) to 0·42 m/day. Heat transport modelling results suggest three kinds of flow around the dams and the meander. Copyright © 2010 John Wiley & Sons, Ltd.     

The Carbonate Radical (HCO3·/CO3–·) as a Reactive Intermediate in Water Chemistry: Kinetics and Modelling

Results from kinetic laboratory studies of reactions of the carbonate radical anion (CO₃^–·) with aromatic compounds in aqueous solution at T = 298 K are presented. Data were obtained in using a laser photolysis laser long-path absorption (LP-LPLA) apparatus which was designed for direct time-resolved studies of radical reactions. For the reactions of CO₃^–· with hydroquinone dimethyl ether (2), methyl anisole (3), benzene (4), p-xylene (5), toluene (6), chlorobenzene (7), nitrobenzene (8), and benzonitrile (9), rate coefficients of k₂ = (3.0 ± 0.6)·10⁷ M^–1 s^–1, k₃ = (9.7 ± 1.7)·10⁵ M^–1 s^–1, k₄ = (3.2 ± 0.7)·10⁵ M^–1 s^–1, k₅ = (3.8 ± 0.9)·10⁴ M^–1 s^–1, k₆ = (6.8 ± 2.3)·10⁴ M^–1 s^–1, k₇ = (2.7 ± 0.6)·10⁵ M^–1 s^–1, k₈ = (1.4 ± 0.5)·10⁴ M^–1 s^–1, and k₉ < 1.3·10² M^–1 s^–1 were obtained. In further studies the effect of temperature on the reactions (2), (4), and (5) has been studied. The kinetic data obtained for the reaction of the carbonate radical anion with aromatic compounds were compared to the corresponding reactions of the hydroxyl radical. Finally, these kinetic data were used within a simple model system to investigate the implications of carbonate radical anion kinetics within water treatment processes. It is shown that the degradation of organic pollutants in ^·OH-radical based water treatment may proceed via the CO₃^–·/HCO_3· radical under certain conditions.     

Interactions between groundwater and surface water in a virginia coastal plain watershed. 1. Hydrological flowpaths

A field study of surface water and groundwater interactions during baseflow and stormflow conditions was performed at the Reedy Creek watershed in the Virginia Coastal Plain. Three estimates of the average saturated hydraulic conductivity (K_s) of the unconfined aquifer were in reasonable agreement (ranging from 0.0033 to 0.010 cm/s), indicating that baseflow in the creek is entirely from the drainage of shallow groundwater from the relatively thin (1–6 m thick) unconfined aquifer. This relatively permeable surficial aquifer was found to be underlain by dark, olive grey, clay-silt and diatomaceous Miocene deposits of low permeability known as the Calvert Formation, which is believed to function as a confining bed in the area. A chemical hydrograph separation technique was used to resolve the contributions of [old] (pre-event) and [new] (event) water to stormflow. Results from a major rainstorm indicated that old water dominated the stormflow response of the watershed, although the new water contribution approached 40% at the hydrograph peak. Stormflow at Reedy Creek appears to result from saturation overland flow from variable source areas which include the stream channels and a significant part of the riparian wetland area. This response appears to be attributable to the transient dynamics of the shallow groundwater flow system and to the formation of localized groundwater mounds which raise the water-table to the wetland surface.     

The effects of rock uplift and rock resistance on river morphology in a subduction zone forearc,Oregon, USA

Relationships between riverbed morphology, concavity, rock type and rock uplift rate are examined to independently unravel the contribution of along-strike variations in lithology and rates of vertical deformation to the topographic relief of the Oregon coastal mountains. Lithologic control on river profile form is reflected by convexities and knickpoints in a number of longitudinal profiles and by general trends of concavity as a function of lithology. Volcanic and sedimentary rocks are the principal rock types underlying the northern Oregon Coast Ranges (between 46°30′ and 45°N) where mixed bedrock–alluvial channels dominate. Average concavity, θ, is 0·57 in this region. In the alluviated central Oregon Coast Ranges (between 45° and 44°N) values of concavity are, on average, the highest (θ = 0·82). South of 44°N, however, bedrock channels are common and θ = 0·73. Mixed bedrock–alluvial channels characterize rivers in the Klamath Mountains (from 43°N south; θ = 0·64). Rock uplift rates of ≥0·5 mm a⁻¹, mixed bedrock–alluvial channels, and concavities of 0·53–0·70 occur within the northernmost Coast Ranges and Klamath Mountains. For rivers flowing over volcanic rocks θ = 0·53, and θ = 0·72 for reaches crossing sedimentary rocks. Whereas channel type and concavity generally co-vary with lithology along much of the range, rivers between 44·5° and 43°N do not follow these trends. Concavities are generally greater than 0·70, alluvial channels are common, and river profiles lack knickpoints between 44·5° and 44°N, despite the fact that lithology is arguably invariant. Moreover, rock uplift rates in this region vary from low, ≤0·5 mm a⁻¹, to subsidence (<0 mm a⁻¹). These observations are consistent with models of transient river response to a decrease in uplift rate. Conversely, the rivers between 44° and 43°N have similar concavities and flow on the same mapped bedrock unit as the central region, but have bedrock channels and irregular longitudinal profiles, suggesting the river profiles reflect a transient response to an increase in uplift rate. If changes in rock uplift rate explain the differences in river profile form and morphology, it is unlikely that rock uplift and erosion are in steady state in the Oregon coastal mountains. Copyright © 2006 John Wiley & Sons, Ltd.     

Seasonal rainfall–runoff relationships in a lowland forested watershed in the southeastern USA

Hydrological processes of lowland watersheds of the southern USA are not well understood compared to a hilly landscape due to their unique topography, soil compositions, and climate. This study describes the seasonal relationships between rainfall patterns and runoff (sum of storm flow and base flow) using 13 years (1964–1976) of rainfall and stream flow data for a low‐gradient, third‐order forested watershed. It was hypothesized that runoff–rainfall ratios (R/P) are smaller during the dry periods (summer and fall) and greater during the wet periods (winter and spring). We found a large seasonal variability in event R/P potentially due to differences in forest evapotranspiration that affected seasonal soil moisture conditions. Linear regression analysis results revealed a significant relationship between rainfall and runoff for wet (r² = 0·68; p < 0·01) and dry (r² = 0·19; p = 0·02) periods. Rainfall‐runoff relationships based on a 5‐day antecedent precipitation index (API) showed significant (r² = 0·39; p < 0·01) correspondence for wet but not (r² = 0·02; p = 0·56) for dry conditions. The same was true for rainfall‐runoff relationships based on 30‐day API (r² = 0·39; p < 0·01 for wet and r² = 0·00; p = 0·79 for dry). Stepwise regression analyses suggested that runoff was controlled mainly by rainfall amount and initial soil moisture conditions as represented by the initial flow rate of a storm event. Mean event R/P were higher for the wet period (R/P = 0·33), and the wet antecedent soil moisture condition based on 5‐day (R/P = 0·25) and 30‐day (R/P = 0·26) prior API than those for the dry period conditions. This study suggests that soil water status, i.e. antecedent soil moisture and groundwater table level, is important besides the rainfall to seasonal runoff generation in the coastal plain region with shallow soil argillic horizons. Copyright © 2011 John Wiley & Sons, Ltd.     

Coarse bed load transport in an alluvial gravel bed stream,Dupuyer Creek,Montana

Coarse bed load was sampled in a gravel/cobble bed stream during two major floods in the snowmelt runoff season. The channel is characterized by high rates of bank erosion and, therefore, high rates of sediment supply and bed load flux. Peak discharge reached four times bank‐full, and bed load was sampled at flows 0·7–1·7 times bank‐full. A large aperture bed load sampler (1 m by 0·45 m) captured the largest particles in motion, and specifically targeted the coarse bed load size distribution by using a relatively large mesh (32 mm or D₂₅ of streambed surface size distribution). Bed load flux was highly variable, with a peak value of 0·85 kg/s/m for the coarse fraction above 38 mm. Bed load size distribution and maximum particle size was related to flow strength. Entrainment was size selective for particles D₇₀ and larger (88–155 mm), while particles in the range D₃₀–D₇₀ (35–88 mm) ceased to move at essentially the same flow. Bed load flux was size selective in that coarse fractions of the streambed surface were under‐represented in or absent from the bed load. Painted tracer particles revealed that the streambed surface in the riffles could remain stable even during high rates of bed load transport. These observations suggest that a large proportion of bed load sediments was sourced from outside the riffles. Repeat surveys confirmed major scour and fill in pools (up to 0·75 m), and bank erosion (>2 m), which together contributed large volumes of sediment to the bed load. Copyright © 2007 John Wiley & Sons, Ltd.     

Persistence of road runoff generation in a logged catchment in Peninsular Malaysia

Measurements of saturated hydraulic conductivity (K_s) and diagnostic model simulations show that all types of logging road/trail in the 14·4 ha Bukit Tarek Experimental Catchment 3 (BTEC3) generate substantial Horton overland flow (HOF) during most storms, regardless of design and level of trafficking. Near‐surface K_s(0–0·05 m) on the main logging road, skid trails and newly constructed logging terraces was less than 1, 2 and 34 mm h^?1, respectively. Near‐surface K_s on an abandoned skid trail in an adjacent basin was higher (62 mm h^?1), owing to the development of a thin organic‐rich layer on the running surface over the past 40 years. Saturated hydraulic conductivity measured at 0·25 m below the surface of all roads was not different (all <6 mm h^?1) and corresponded to the K_s of the adjacent hillslope subsoil, as most roads were excavated into the regolith more than 0·5–1 m. After 40 years, only limited recovery in near‐surface K_s occurred on the abandoned skid trail. This road generated HOF after the storage capacity of the upper near‐surface layer was exceeded during events larger than about 20 mm. Thus, excavation into low‐K_s substrate had a greater influence on the persistence of surface runoff production than did surface compaction by machinery during construction and subsequent use during logging operations. Overland flow on BTEC3 roads was also augmented by the interception of shallow subsurface flow traveling along the soil–saprolite/bedrock interface and return flow emerging from the cutbank through shallow biogenic pipes. The most feasible strategy for reducing long‐term road‐related impacts in BTEC3 is limiting the depth of excavation and designing a more efficient road network, including minimizing the length and connectivity of roads and skid trails. Copyright © 2007 John Wiley & Sons, Ltd.     

Salinograph trends as indicators of the recession characteristics of stream components

Compared to hydrograph recession analysis, which is widely applied in engineering hydrology, the quantitative assessment of stream salinity with time (i.e. the salinograph) has received significantly less attention. In particular, while in many previous hydrological studies an inverse relationship between hydrograph and salinograph responses is apparent, the concept of salinity accession (the inversely related salinity counterpart to hydrograph recession) has not been introduced nor quantitatively evaluated in previous literature. In this study, we conduct a mathematical analysis of salinograph accession, and determine new quantitative relationships between salinity accession and hydrograph recession parameters. An equation is formulated that reproduces the general trend in salinity accession. A salinity accession parameter k_c is then introduced and is shown to be the ratio of direct runoff to total stream flow recession parameters: k_r/k. The groundwater recession parameter k_g was estimated using a simple and rapid method that uses both salinograph and hydrograph data. Salinity accession type‐curves illustrate that under certain conditions, the relative steepness of individual salinographs is dependent upon the ratio of groundwater salinity to direct runoff salinity: C_g/C_r. The salinity accession algorithms are applied to two contrasting field settings: Scott Creek, South Australia and Sandy Creek, northern Queensland, Australia. It was found that k_g > k during periods of obvious stream flow recession, for the events analysed. Salinograph accession behaviour was fairly similar for both sites, despite contrasting environments. Using assumed end‐member salinities for groundwater and direct runoff based upon field observations, the behaviour of k_c from the Scott Creek site was approximately reproduced by varying the initial groundwater to runoff flow ratio: Q_g0/Q_r0, within reasonable parameter ranges. The use of salinograph information when used in addition to standard hydrograph analyses provided useful information on recession characteristics of stream components. Copyright © 2008 John Wiley & Sons, Ltd.     

Environmental isotopic and hydrochemical characteristics of groundwater systems in Daying and Qicun geothermal fields,Xinzhou Basin,Shanxi, China

The conceptual hydrogeological model of the low to medium temperature Daying and Qicun geothermal fields has been proposed, based on hydrochemical characteristics and isotopic compositions. The two geothermal fields are located in the Xinzhou basin of Shanxi, China and exhibit similarities in their broad‐scale flow patterns. Geothermal water is derived from the regional groundwater flow system of the basin and is characterized by Cl·SO₄‐Na type. Thermal water is hydrochemically distinct from cold groundwater having higher total dissolved solids (TDS) (>0·8 g/l) and Sr contents, but relatively low Ca, Mg and HCO₃ contents. Most shallow groundwater belongs to local flow systems which are subject to evaporation and mixing with irrigation returns. The groundwater residence times estimated by tritium and ¹⁴C activities indicate that deep non‐thermal groundwater (130–160 m) in the Daying region range from modern (post‐1950s) in the piedmont area to more than 9·4 ka BP (Before Present) in the downriver area and imply that this water belong to an intermediate flow system. Thermal water in the two geothermal fields contains no detectable active ¹⁴C, indicating long residence times (>50 ka), consistent with this water being part of a large regional flow system. The mean recharge elevation estimated by using the obtained relationship Altitude (m) = ? 23·8 × δ²H (‰ ) ? 121·3, is 1980 and 1880 m for the Daying and Qicun geothermal fields, respectively. The annual infiltration rates in the Daying and Qicun geothermal fields can be estimated to be 9029 × 10³ and 4107 × 10³ m³/a, respectively. The variable ⁸⁶Sr/⁸⁷Sr values in the thermal and non‐thermal groundwater in the two fields reflect different lithologies encountered along the flow path(s) and possibly different extents of water‐rock interaction. Based on the analysis of groundwater flow systems in the two geothermal fields, hydrogeochemical inverse modelling was performed to indicate the possible water‐rock interaction processes that occur under different scenarios. Copyright © 2010 John Wiley & Sons, Ltd.     

Characteristics of preferential flow and groundwater discharge to Shingobee Lake,Minnesota, USA

Small‐scale heterogeneities and large changes in hydraulic gradient over short distances can create preferential groundwater flow paths that discharge to lakes. A 170 m² grid within an area of springs and seeps along the shore of Shingobee Lake, Minnesota, was intensively instrumented to characterize groundwater‐lake interaction within underlying organic‐rich soil and sandy glacial sediments. Seepage meters in the lake and piezometer nests, installed at depths of 0·5 and 1·0 m below the ground surface and lakebed, were used to estimate groundwater flow. Statistical analysis of hydraulic conductivity estimated from slug tests indicated a range from 21 to 4·8 × 10^?3 m day^?1 and small spatial correlation. Although hydraulic gradients are overall upward and toward the lake, surface water that flows onto an area about 2 m onshore results in downward flow and localized recharge. Most flow occurred within 3 m of the shore through more permeable pathways. Seepage meter and Darcy law estimates of groundwater discharge agreed well within error limits. In the small area examined, discharge decreases irregularly with distance into the lake, indicating that sediment heterogeneity plays an important role in the distribution of groundwater discharge. Temperature gradients showed some relationship to discharge, but neither temperature profiles nor specific electrical conductance could provide a more convenient method to map groundwater–lake interaction. These results suggest that site‐specific data may be needed to evaluate local water budget and to protect the water quality and quantity of discharge‐dominated lakes. Copyright © 2002 John Wiley & Sons, Ltd.     

Wetland nitrogen dynamics in an Adirondack forested watershed

Wetlands often form the transition zone between upland soils and watershed streams, however, stream–wetland interactions and hydrobiogeochemical processes are poorly understood. We measured changes in stream nitrogen (N) through one riparian wetland and one beaver meadow in the Archer Creek watershed in the Adirondack Mountains of New York State, USA from 1 March to 31 July 1996. In the riparian wetland we also measured changes in groundwater N. Groundwater N changed significantly from tension lysimeters at the edge of the peatland to piezometer nests within the peatland. Mean N concentrations at the peatland perimeter were 1·5, 0·5 and 18·6 µmol L^?1 for NH₄⁺, NO₃^? and DON (dissolved organic nitrogen), respectively, whereas peatland groundwater N concentration was 56·9, 1·5 and 31·6 µmol L^?1 for NH₄⁺, NO₃^? and DON, respectively. The mean concentrations of stream water N species at the inlet to the wetlands were 1·5, 10·1 and 16·9 µmol L^?1 for NH₄⁺, NO₃^? and DON, respectively and 1·6, 28·1 and 8·4 µmol L^?1 at the wetland outlet. Although groundwater total dissolved N (TDN) concentrations changed more than stream water TDN through the wetlands, hydrological cross‐sections for the peatland showed that wetland groundwater contributed minimally to stream flow during the study period. Therefore, surface water N chemistry was affected more by in‐stream N transformations than by groundwater N transformations because the in‐stream changes, although small, affected a much greater volume of water. Stream water N input–output budgets indicated that the riparian peatland retained 0·16 mol N ha^?1 day^?1 of total dissolved N and the beaver meadow retained 0·26 mol N ha^?1 day^?1 during the study period. Nitrate dominated surface water TDN flux from the wetlands during the spring whereas DON dominated during the summer. This study demonstrates that although groundwater N changed significantly in the riparian peatland, those changes were not reflected in the stream. Consequently, although in‐stream changes of N concentrations were less marked than those in groundwater, they had a greater effect on stream water chemistry—because wetland groundwater contributed minimally to stream flow. Copyright © 2004 John Wiley & Sons, Ltd.     

In situ quantification of spatial and temporal variability of hyporheic exchange in static and mobile gravel‐bed rivers

Seepage meters modified for use in flowing water were used to directly measure rates of exchange between surface and subsurface water in a gravel‐ and cobble bed river in western Pennsylvania, USA (Allegheny River, Q_mean = 190 m³/s) and a sand‐ and gravel‐bed river in Colorado, USA (South Platte River, Q_mean = 9·7 m³/s). Study reaches at the Allegheny River were located downstream from a dam. The bed was stable with moss, algae, and river grass present in many locations. Median seepage was + 0·28 m/d and seepage was highly variable among measurement locations. Upward and downward seepage greatly exceeded the median seepage rate, ranging from + 2·26 (upward) to ? 3·76 (downward) m/d. At the South Platte River site, substantial local‐scale bed topography as well as mobile bedforms resulted in spatial and temporal variability in seepage greatly in exceedence of the median groundwater discharge rate of 0·24 m/d. Both upward and downward seepage were recorded along every transect across the river with rates ranging from + 2·37 to ? 3·40 m/d. Despite a stable bed, which commonly facilitates clogging by fine‐grained or organic sediments, seepage rates at the Allegheny River were not reduced relative to those at the South Platte River. Seepage rate and direction depended primarily on measurement position relative to local‐ and meso‐scale bed topography at both rivers. Hydraulic gradients were small at nearly all seepage‐measurement locations and commonly were not a good indicator of seepage rate or direction. Therefore, measuring hydraulic gradient and hydraulic conductivity at in‐stream piezometers may be misleading if used to determine seepage flux across the sediment‐water interface. Such a method assumes that flow between the well screen and sediment‐water interface is vertical, which appears to be a poor assumption in coarse‐grained hyporheic settings. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical simulation of three‐dimensional flow hydraulics in a braided channel

Results are presented from a numerical simulation of three‐dimensional flow hydraulics around a mid‐channel bar carried out using the FLUENT/UNS computational fluid dynamics (CFD) software package. FLUENT/UNS solves the three‐dimensional Reynolds‐averaged form of the Navier–Stokes equations. Turbulence closure is achieved using a RNG k–ϵ model. Simulated flow velocities are compared with measured two‐dimensional velocities (downstream and cross‐stream) obtained using an electromagnetic current meter (ECM). The results of the simulation are qualitatively consistent with the flow structures observed in the field. Quantitative comparison of the simulated and measured velocity magnitudes indicates a strong positive correlation between the two (r=0·88) and a mean difference of 0·09 m s⁻¹. Deviations between simulated and measured velocities may be identified that are both random and systematic. The former may reflect a number of factors including subgrid‐scale natural spatial variability in flow velocities associated with local bed structures and measurement uncertainty resulting from problems of ECM orientation. Model mesh configuration, roughness parameterization and inlet boundary condition uncertainty may each contribute to systematic differences between simulated and measured flow velocities. These results illustrate the potential for using CFD software to simulate flow hydraulics in natural channels with complex configurations. They also highlight the need for detailed spatially distributed datasets of three‐dimensional flow variables to establish the accuracy and applicability of CFD software. Copyright © 1999 John Wiley & Sons, Ltd.     

A modelling study of hyporheic exchange pattern and the sequence,size, and spacing of stream bedforms in mountain stream networks,Oregon, USA

Studies of hyporheic exchange flows have identified physical features of channels that control exchange flow at the channel unit scale, namely slope breaks in the longitudinal profile of streams that generate subsurface head distributions. We recently completed a field study that suggested channel unit spacing in stream longitudinal profiles can be used to predict the spacing between zones of upwelling (flux of hyporheic water into the stream) and downwelling (flux of stream water into the hyporheic zone) in the beds of mountain streams. Here, we use two‐dimensional groundwater flow and particle tracking models to simulate vertical and longitudinal hyporheic exchange along the longitudinal axis of stream flow in second‐, third‐, and fourth‐order mountain stream reaches. Modelling allowed us to (1) represent visually the effect that the shape of the longitudinal profile has on the flow net beneath streambeds; (2) isolate channel unit sequence and spacing as individual factors controlling the depth that stream water penetrates the hyporheic zone and the length of upwelling and downwelling zones; (3) evaluate the degree to which the effects of regular patterns in bedform size and sequence are masked by irregularities in real streams. We simulated hyporheic exchange in two sets of idealized stream reaches and one set of observed stream reaches. Idealized profiles were constructed using regression equations relating channel form to basin area. The size and length of channel units (step size, pool length, etc.) increased with increasing stream order. Simulations of hyporheic exchange flows in these reaches suggested that upwelling lengths increased (from 2·7 m to 7·6 m), and downwelling lengths increased (from 2·9 m to 6·0 m) with increase in stream order from second to fourth order. Step spacing in the idealized reaches increased from 5·3 m to 13·7 m as stream size increased from second to fourth order. Simulated downwelling lengths increased from 4·3 m in second‐order streams to 9·7 m in fourth‐order streams with a POOL–RIFFLE–STEP channel unit sequence, and increased from 2·5 m to 6·1 m from second‐ to fourth‐order streams with a POOL–STEP–RIFFLE channel unit sequence. Upwelling lengths also increased with stream order in these idealized channels. Our results suggest that channel unit spacing, size, and sequence are all important in determining hyporheic exchange patterns of upwelling and downwelling. Though irregularities in the size and spacing of bedforms caused flow nets to be much more complex in surveyed stream reaches than in idealized stream reaches, similar trends emerged relating the average geomorphic wavelength to the average hyporheic wavelength in both surveyed and idealized reaches. This article replaces a previously published version (Hydrological Processes, 19 (17), 2915–2929 (2005) [ DOI:10.1002/hyp.5790 ]. See also retraction notice DOI:10.1002/hyp.6350 Copyright © 2006 John Wiley & Sons, Ltd.     

Hydrological processes controlling ground and surface water flow from a hypermaritime forest–peatland complex,Diana Lake Provincial Park,British Columbia,Canada

The proposed harvesting of previously undeveloped forests in north coastal British Columbia requires an understanding of hydrological responses. Hydrometric and isotopic techniques were used to examine the hydrological linkages between meteoric inputs to the surface‐groundwater system and runoff response patterns of a forest‐peatland complex. Quickflow accounted for 72–91% of peak storm discharge. The runoff ratio was lowest for open peatland areas with thick organic horizons (0·02–0·05) due to low topographic gradients and many surface depressions capable of retaining surface water. Runoff ratio increased comparatively for ephemeral surface seep flows (0·06–0·40) and was greatest in steeply sloping forest communities with more permeable soils (0·33–0·69). The dominant mechanism for runoff generation was saturated shallow subsurface flow. Groundwater fluxes from the organic horizon of seeps (1·70–1·72 m³ day^?1 m^?1) were an important component of quickflow. The homogeneous δ²H? δ¹⁸O composition of groundwater indicated attenuation of the seasonal rainfall signal by mixing during recharge. The positive correlation (r² = 0·64 and 0·38, α = 0·05) between slope index and δ¹⁸O values in groundwater suggests that the spatial pattern in the δ¹⁸O composition along the forest‐peatland complex is influenced by topography and provides evidence that topographic indices may be used to predict groundwater residence time. Copyright © 2006 John Wiley & Sons, Ltd.     

Flow resistance in steep mountain streams

Resistance to flow at low to moderate stream discharge was examined in five small (12–77 km² drainage area) tributaries of Chilliwack River, British Columbia, more than half of which exhibit planar bed morphology. The resulting data set is composed of eight to 12 individual estimates of the total resistance to flow at 61 cross sections located in 13 separate reaches of five tributaries to the main river. This new data set includes 625 individual estimates of resistance to flow at low to moderate river stage. Resistance to flow in these conditions is high, highly variable and strongly dependent on stage. The Darcy–Weisbach resistance factor (ff) varies over six orders of magnitude (0·29–12 700) and Manning's n varies over three orders of magnitude (0·047–7·95). Despite this extreme range, both power equations at the individual cross sections and Keulegan equations for reach‐averaged values describe the hydraulic relations well. Roughness is divided into grain and form (considered as all non‐grain sources) components. Form roughness is the dominant component, accounting for about 90% of the total roughness of the system (i.e., form roughness is on average 8.6 times as great as grain roughness). Of the various quantitative and qualitative form‐roughness indicators observed, only the sorting coefficient (σ = D₈₄/D₅₀) correlates well with form roughness. Copyright © 2008 John Wiley & Sons, Ltd.     

Controls on effective settling velocities of mud flocs in the Changjiang Estuary,China

Theoretical and experimental studies were undertaken to gain insight into physical parameters controlling the flocculation and settling properties of mud flocs in the Changjiang Estuary, China. The Rouse equation is applied to vertical profiles of suspended sediment concentration to determine the bulk mean settling velocity (w_s) of sediment suspended in the Changjiang Estuary. Both in situ point‐sampled and acoustically measured profiles of suspended mud concentrations were fit selectively. The calculated settling velocities w_s mainly ranged from 0·4 to 4·1 mm s^?1 for the point‐sampled data set, and from 1·0 to 3·0 mm s^?1 for the acoustically measured data set. Furthermore, the settling velocities of mud flocs increased with mean concentration (C?) of mud flocs in suspension and were proportional to increasing bottom shear stress (τ_b) of tidal flow. The best equation for the field settling velocity of mud flocs in the Changjiang Estuary can be expressed by the power law: w_s = mC?ⁿ (m, 1·14–2·37; n, 0·84–1·03). It is suggested that C? and τ_b were the dominant physical parameters controlling the flocculation and w_s of mud flocs in suspension. Copyright © 2004 John Wiley & Sons, Ltd.