Home About us Contact | |||
Flow Depth (flow + depth)
Selected AbstractsEstimation of temporal variation in splash detachment in two Japanese cypress plantations of contrasting ageEARTH SURFACE PROCESSES AND LANDFORMS, Issue 9 2010Y. Wakiyama Abstract To elucidate splash erosion processes under natural rainfall conditions, temporal variations in splash detachment were observed using a piezoelectric saltation sensor (H11B; Sensit Co., Portland, ND, USA). Preliminary laboratory tests of Sensit suggested that they were suitable for field observations. Field observations were conducted between July and September 2006 in 21- and 36-year-old Japanese cypress (Chamaecyparis obtusa) plantations with mean stand heights of 9·2,m and 17·4,m, respectively. Splash detachment (in g m,2) was measured seven times using splash cups, and raindrop kinetic energy (in J,m,2,mm,1) in both stands was measured using laser drop-sizing (LD) gauges. Sensit was installed to record saltation counts, which were converted to temporal data of splash detachment (splash rate; in g m,2 10,min,1) using the relationship between splash detachment and saltation counts. Surface runoff was monitored using runoff plots of 0·5,m width and 2·0,m length to obtain temporal data of flow depth (in millimeters). Both total splash detachment and raindrop kinetic energy were larger in the older stand. Increased splash rates per unit throughfall were found in both stands after rainless durations longer than approximately one day in both stands. However, a lower splash rate was found in the 21-year stand after rainfall events. During extreme rainstorms, the 21-year stand showed a low runoff rate and a decline in the splash rate, while the 36-year stand showed a higher splash rate and increased flow depth. The piezoelectric sensor proved to be a useful means to elucidate splash erosion processes in field conditions. Copyright © 2010 John Wiley & Sons, Ltd. [source] Reduced-complexity flow routing models for sinuous single-thread channels: intercomparison with a physically-based shallow-water equation modelEARTH SURFACE PROCESSES AND LANDFORMS, Issue 5 2009A. P. Nicholas Abstract Reduced-complexity models of fluvial processes use simple rules that neglect much of the underlying governing physics. This approach is justified by the potential to use these models to investigate long-term and/or fundamental river behaviour. However, little attention has been given to the validity or realism of reduced-complexity process parameterizations, despite the fact that the assumptions inherent in these approaches may limit the potential for elucidating the behaviour of natural rivers. This study presents two new reduced-complexity flow routing schemes developed specifically for application in single-thread rivers. Output from both schemes is compared with that from a more sophisticated model that solves the depth-averaged shallow water equations. This comparison provides the first demonstration of the potential for deriving realistic predictions of in-channel flow depth, unit discharge, energy slope and unit stream power using simple flow routing schemes. It also highlights the inadequacy of modelling unit stream power, shear stress or sediment transport capacity as a function of local bed slope, as has been common practice in a number of previous reduced-complexity models. Copyright © 2009 John Wiley & Sons, Ltd. [source] Meltwater discharge through the subglacial bed and its land-forming consequences from numerical experiments in the Polish lowland during the last glaciationEARTH SURFACE PROCESSES AND LANDFORMS, Issue 4 2009Jan A. Piotrowski Abstract Numerical experiments suggest that the last glaciation severely affected the upper lithosphere groundwater system in NW Poland: primarily its flow pattern, velocities and fluxes. We have simulated subglacial groundwater flow in two and three spatial dimensions using finite difference codes for steady-state and transient conditions. The results show how profoundly the ice sheet modifies groundwater pressure heads beneath and some distance beyond the ice margin. All model runs show water discharge at the ice forefield driven by ice-sheet-thickness-modulated, down-ice-decreasing hydraulic heads. In relation to non-glacial times, the transient 3D model shows significant changes in the groundwater flow directions in a regionally extensive aquifer ca. 90 m below the ice,bed interface and up to 40 km in front of the glacier. Comparison with empirical data suggests that, depending on the model run, only between 5 and 24% of the meltwater formed at the ice sole drained through the bed as groundwater. This is consistent with field observations documenting abundant occurrence of tunnel valleys, indicating that the remaining portion of basal meltwater was evacuated through a channelized subglacial drainage system. Groundwater flow simulation suggests that in areas of very low hydraulic conductivity and adverse subglacial slopes water ponding at the ice sole was likely. In these areas the relief shows distinct palaeo-ice lobes, indicating fast ice flow, possibly triggered by the undrained water at the ice,bed interface. Owing to the abundance of low-permeability strata in the bed, the simulated groundwater flow depth is less than ca. 200 m. Copyright © 2009 John Wiley & Sons, Ltd. [source] Flow energy and channel adjustments in rills developed in loamy sand and sandy loam soilsEARTH SURFACE PROCESSES AND LANDFORMS, Issue 1 2009Jovan R. Stefanovic Abstract The storms usually associated with rill development in nature are seldom prolonged, so development is often interrupted by interstorm disturbances, e.g. weathering or tillage. In laboratory simulated rainfall experiments, active rill development can be prolonged, and under these conditions typically passes through a period of intense incision, channel extension and bifurcation before reaching quasi-stable conditions in which little form change occurs. This paper presents laboratory experiments with coarse textured soils under simulated rainfall which show how channel adjustment processes contribute to the evolution of quasi-stability. Newly incised rills were stabilized for detailed study of links between rill configuration and flow energy. On a loamy sand, adjustment towards equilibrium occurred due to channel widening and meandering, whereas on a sandy loam, mobile knickpoints and chutes, pulsations in flow width and flow depth and changes in stream power and sediment discharge occurred as the channel adjusted towards equilibrium. The tendency of rill systems towards quasi-stability is shown by changes in stream power values which show short-lived minima. Differences in energy dissipation in stabilized rills indicate that minimization of energy dissipation was reached locally between knickpoints and at the downstream ends of rills. In the absence of energy gradients in knickpoints and chutes, stabilized rill sections tended toward equilibrium by establishing uniform energy expenditure. The study confirmed that energy dissipation increased with flow aspect ratio. In stabilized rills, flow acceleration reduced energy dissipation on the loamy sand but not on the sandy loam. On both soils flow deceleration tended to increase energy dissipation. Understanding how rill systems evolve towards stability is essential in order to predict how interruptions between storms may affect long-term rill dynamics. This is essential if event-based physical models are to become effective in predicting sediment transport on rilled hillslopes under changing weather and climatic conditions. Copyright © 2008 John Wiley and Sons, Ltd. [source] Organic litter: dominance over stones as a source of interrill flow roughness on low-gradient desert slopes at Fowlers Gap, arid western NSW, AustraliaEARTH SURFACE PROCESSES AND LANDFORMS, Issue 1 2003David Dunkerley Abstract Thirty-six runoff plot experiments provide data on flow depths, speeds, and Darcy,Weisbach friction coefficients (f) on bare soil surfaces, and surfaces to which were added sufficient extra plant litter or surface stones to provide projected cover of 5, 10 and 20 per cent. Precision flow depth data were derived with a computer-controlled gantry and needle gauge for two different discharges for each plot treatment. Taking a fixed flow intensity (Reynolds number, Re = 150) for purposes of comparison shows means of f = 17·7 for bare soil surfaces, f = 11·4 for added stone treatments, and f = 23·8 for added litter treatments. Many individual values of f for stone treatments are lower than for the bare soil surface, but all litter treatments show increases in fcompared to bare soil. The lowering of f in stone treatments relates to the submerged volume that the stones occupied, and the associated concentration of flow onto a smaller part of the plot surface. This leads to locally higher flow intensities and lower frictional drag along threads of flow that the obstacles create. Litter causes higher frictional drag because the particles are smaller, and, for the same cover fraction, are 100 times more numerous and provide 20 times the edge or perimeter length. Along these edges, which in total exceed 2·5 m g,1 (equivalent to 500 m m,2 for a loading of 2 t ha,1), surface tension draws up water from between the litter particles. This reduces flow depth there, and as a consequence of the lower flow intensity, frictional drag rises. Furthermore, no clear passage remains for the establishment of flow threads. These findings apply to shallow interrill flows in which litter is largely immobile. The key new result from these experiments is that under these conditions, a 20 per cent cover of organic litter can generate interrill frictional retardation that exceeds by nearly 41 per cent that of a bare soil surface, and twice that contributed by the same cover fraction of surface stones. Even greater dominance by litter can be anticipated at the many dryland sites where litter covers exceed those tested here. Copyright © 2002 John Wiley & Sons, Ltd. [source] Estimating the mean speed of laminar overland flow using dye injection-uncertainty on rough surfacesEARTH SURFACE PROCESSES AND LANDFORMS, Issue 4 2001David Dunkerley Abstract A common method for estimating mean flow speeds in studies of surface runoff is to time the travel of a dye cloud across a measured flow path. Motion of the dye front reflects the surface flow speed, and a correction must be employed to derive a value for the profile mean speed, which is always lower. Whilst laminar flow conditions are widespread in the interrill zone, few data are available with which to establish the relationship linking surface and profile mean speeds, and there are virtually none for the flow range 100,<,Re,<,500 (Re,=,Reynolds number) which is studied here. In laboratory experiments on a glued sand board, mean flow speeds were estimated from both dye speeds and the volumetric flow relation v,=,Q/wd with d measured using a computer-controlled needle gauge at 64 points. In order to simulate conditions applicable to many dryland soils, the board was also roughened with plant litter and with ceramic tiles (to simulate surface stone cover). Results demonstrate that in the range 100,<,Re,<,500, there is no consistent relation between surface flow speeds and the profile mean. The mean relationship is v,=,0·56 vsurf, which departs significantly from the theoretical smooth-surface relation v,=,0·67 vsurf, and exhibits a considerable scatter of values that show a dependence on flow depth. Given the inapplicability of any fixed conversion factor, and the dependence on flow depth, it is suggested that the use of dye timing as a method for estimating v be abandoned in favour of precision depth measurement and the use of the relation v,=,Q/wd, at least within the laminar flow range tested. Copyright © 2001 John Wiley & Sons, Ltd. [source] Hydrodynamics and geomorphic work of jökulhlaups (glacial outburst floods) from Kverkfjöll volcano, IcelandHYDROLOGICAL PROCESSES, Issue 6 2007Jonathan L. Carrivick Abstract Jökulhlaups (glacial outburst floods) occur frequently within most glaciated regions of the world and cause rapid landscape change, infrastructure damage, and human disturbance. The largest jökulhlaups known to have occurred during the Holocene within Iceland drained from the northern margin of Vatnajökull and along the Jökulsá á Fjöllum. Some of these jökulhlaups originated from Kverkfjöll volcano and were routed through anastomosing, high gradient and hydraulically rough channels. Landforms and sediments preserved within these channels permit palaeoflow reconstructions. Kverkfjöll jökulhlaups were reconstructed using palaeocompetence (point measurements), slope,area (one-dimensional), and depth-averaged two-dimensional (2D) hydrodynamic modelling techniques. The increasing complexity of 2D modelling required a range of assumptions, but produced information on both spatial and temporal variations in jökulhlaup characteristics. The jökulhlaups were volcanically triggered, had a linear-rise hydrograph and a peak discharge of 50 000,100 000 m3 s,1, which attenuated by 50,75% within 25 km. Frontal flow velocities were ,2 m s,1; but, as stage increased, velocities reached 5,15m s,1. Peak instantaneous shear stress and stream power reached 1 × 104 N m,2 and 1 × 105 W m,2 respectively. Hydraulic parameters can be related to landform groups. A hierarchy of landforms is proposed, ranging from the highest energy zones (erosional gorges, scoured bedrock, cataracts, and spillways) to the lowest energy zones (of valley fills, bars, and slackwater deposits). Fluvial erosion of bedrock occurred in Kverkfjallarani above ,3 m flow depth, ,7m s,1 flow velocity, ,1 × 102 N m,2 shear stress, and 3 × 102 W m,2 stream power. Fluvial deposition occurred in Kverkfjallarani below ,8 m flow depth, 11 m s,1 flow velocity, 5 × 102 N m,2 shear stress, and 3 × 103 W m,2 stream power. Hence, erosional and depositional ,envelopes' have considerable overlap, probably due to transitional flow phenomena and the influence of upstream effects, such as hydraulic ponding and topographic constrictions, for example. Holocene Kverkfjöll jökulhlaups achieved geomorphic work comparable to that of other late Pleistocene ,megafloods'. This work was a result of steep channel gradients, topographic channel constrictions, and high hydraulic roughness, rather than to extreme peak discharges. The Kverkfjöll jökulhlaups have implications for landscape evolution in north-central Iceland, for water-sediment inputs into the North Atlantic, and for recognizing jökulhlaups in the rock record. 2D hydrodynamic modelling is likely to be important for hazard mitigation in similar landscapes and upon other glaciated volcanoes, because it only requires an input hydrograph and a digital elevation model to run a model, rather than suites of geomorphological evidence and field-surveyed valley cross-sections, for example. Copyright © 2006 John Wiley & Sons, Ltd. [source] Errors of kinematic wave and diffusion wave approximations for time-independent flows with infiltration and momentum exchange includedHYDROLOGICAL PROCESSES, Issue 9 2005V. P. Singh Abstract Error equations for kinematic wave and diffusion wave approximations were derived for time-independent flows on infiltrating planes and channels under one upstream boundary and two downstream boundary conditions: zero flow at the upstream boundary, and critical flow depth and zero depth gradient at the downstream boundary. These equations specify error in the flow hydrograph as a function of space. The diffusion wave approximation was found to be in excellent agreement with the dynamic wave approximation, with errors below 2% for values of KF (e.g. KF , 7·5), where K is the kinematic wave number and F is the Froude number. Even for small values of KF (e.g. KF = 2·5), the errors were typically less than 3%. The accuracy of the diffusive approximation was greatly influenced by the downstream boundary condition. For critical flow depth downstream boundary condition, the error of the kinematic wave approximation was found to be less than 10% for KF , 7·5 and greater than 20% for smaller values of KF. This error increased with strong downstream boundary control. The analytical solution of the diffusion wave approximation is adequate only for small values of K. Copyright © 2005 John Wiley & Sons, Ltd. [source] Suspended sediment concentration and the ripple,dune transitionHYDROLOGICAL PROCESSES, Issue 17 2004Robert J. Schindler Abstract Flume experiments were conducted in order to monitor changes in flow turbulence intensity and suspended sediment concentration at seven stages across the ripple,dune transition and at three different positions above the bed surface. Three-dimensional velocity measurements were obtained using an acoustic Doppler velocimeter (ADV). Suspended sediment concentration (SSC) was monitored indirectly using ADV signal amplitude. Although limited to time-averaged parameters, the analysis reveals that SSC varies significantly with stage across the transition and with sampling height. The statistical analysis also reveals an apparent uniformity of suspended sediment concentration with height above the bed in the lower half of the flow depth at the critical stage in the transition from ripples to dunes. This is also the stage at which turbulence intensity is maximized. Statistically significant correlations were also observed between suspended sediment concentrations and root-mean-square values of vertical velocity fluctuations. These correlations reflect the various levels of shear-layer activity and the distinct turbulent flow regions across the transition. Conversely, time-averaged values of Reynolds shear stress exhibit a very weak relationship with suspended sediment concentrations. Copyright © 2004 John Wiley & Sons, Ltd. [source] Numerical simulation of overbank processes in topographically complex floodplain environmentsHYDROLOGICAL PROCESSES, Issue 4 2003A. P. Nicholas Abstract This article presents results from an investigation of the hydraulic characteristics of overbank flows on topographically-complex natural river floodplains. A two-dimensional hydraulic model that solves the depth-averaged shallow water form of the Navier,Stokes equations is used to simulate an overbank flow event within a multiple channel reach of the River Culm, Devon, UK. Parameterization of channel and floodplain roughness by the model is evaluated using monitored records of main channel water level and point measurements of floodplain flow depth and unit discharge. Modelled inundation extents and sequences are assessed using maps of actual inundation patterns obtained using a Global Positioning System, observational evidence and ground photographs. Simulation results suggest a two-phase model of flooding at the site, which seems likely to be representative of natural floodplains in general. Comparison of these results with previous research demonstrates the complexity of overbank flows on natural river floodplains and highlights the limitations of laboratory flumes as an analogue for these environments. Despite this complexity, frequency distributions of simulated depth, velocity and unit discharge data closely follow a simple gamma distribution model, and are described by a shape parameter (,) that exhibits clear systematic trends with changing discharge and floodplain roughness. Such statistical approaches have the potential to provide the basis for computationally efficient flood routing and overbank sedimentation models. Copyright © 2002 John Wiley & Sons, Ltd. [source] Depth-integrated, non-hydrostatic model for wave breaking and run-upINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 5 2009Yoshiki Yamazaki Abstract This paper describes the formulation, verification, and validation of a depth-integrated, non-hydrostatic model with a semi-implicit, finite difference scheme. The formulation builds on the nonlinear shallow-water equations and utilizes a non-hydrostatic pressure term to describe weakly dispersive waves. A momentum-conserved advection scheme enables modeling of breaking waves without the aid of analytical solutions for bore approximation or empirical equations for energy dissipation. An upwind scheme extrapolates the free-surface elevation instead of the flow depth to provide the flux in the momentum and continuity equations. This greatly improves the model stability, which is essential for computation of energetic breaking waves and run-up. The computed results show very good agreement with laboratory data for wave propagation, transformation, breaking, and run-up. Since the numerical scheme to the momentum and continuity equations remains explicit, the implicit non-hydrostatic solution is directly applicable to existing nonlinear shallow-water models. Copyright © 2008 John Wiley & Sons, Ltd. [source] A physical, mechanistic and fully coupled hillslope hydrology modelINTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 11 2005David A. Raff Abstract We present the mathematical development and numerical solution of a new model of flow processes on an infiltrating hillslope. We also present validation and sample applications. The model is a distributed, mechanistic, physically based hillslope hydrologic model. The model describes the small-scale processes associated with overland flow, erosion, and sediment transport on an infiltrating surface and is capable of capturing small-scale variations in flow depth, flow velocities, interactive infiltration, erosion rates, and sediment transport. The model couples the fully two-dimensional hydrodynamic equations for overland flow, the one-dimensional Richards equation for infiltration, and a sediment detachment and transport model. Two simulations are presented highlighting the model's ability to capture and describe the interaction between precipitation, overland flow, erosion and infiltration at very small scales. Results of the two-dimensional simulations indicate the system of equations produces hillslopes possessing characteristics of self-organization as observed in real world systems. Copyright © 2005 John Wiley & Sons, Ltd. [source] Spatially varied flow in non-prismatic channels.IRRIGATION AND DRAINAGE, Issue 1 2002I: dynamic equation canal d'irrigation; modélisation; écoulement spatialement varié; équation de transfert; analyse numérique; expérimentation Abstract The bottom width of channels carrying spatially varied flow with increasing discharge is usually flared in the flow direction. This produces a non-prismatic section. This paper, based on the law of linear momentum conservation, presents a new form of the governing dynamic equation for flow in such circumstances. In addition to the computed water surface profile, the proposed equation affects the position of the control section and its flow depth as well. Copyright © 2002 John Wiley & Sons, Ltd. RÉSUMÉ La largeur du lit de cours d'eau, variée spatialement par l'augmentation du débit, influence la direction de l'écoulement. Cette procédure crée une coupe irréguliére de cours d'eau. Cet article présente une nouvelle forme de l'équation dynamique dans cette situation de l'écoulement. La méthodologie est basée sur la loi de conservation linéaire de momentum. L'équation proposée est capable de calculer le profil de niveau de l'eau. Elle affecte aussi sur la situation de coupe de contrôle et la profondeur de l'écoulement. Copyright © 2002 John Wiley & Sons, Ltd. [source] Indian Ocean tsunamis: environmental and socio-economic impacts in Langkawi, MalaysiaTHE GEOGRAPHICAL JOURNAL, Issue 2 2007MICHAEL BIRD We report the results of a study of the physical characteristics and socio-economic impacts of the Indian Ocean Tsunami of 26 December 2004 on the tourist island of Langkawi, Malaysia. In comparison with many other locations struck by the tsunami, the immediate physical and socio-economic impacts in Langkawi were relatively minor. A detailed survey of the watermark and ground elevations was undertaken in the worst affected area between Sungei Kuala Teriang and Sungei Kuala Melaka. Here, the tsunami reached a maximum elevation of 4.29 m as it crossed the coast, with a maximum flow depth of 2.0 m and a very consistent run-up elevation relative to mean sea level of 300 ± 10 cm. The tsunami inundated inshore areas for 300 m and penetrated inland along creeks for 500,1000 m. Structural damage to buildings was confined to within 50,150 m of the shoreline where about 10% of the houses were completely destroyed and 60,70% suffered significant structural damage. Damage was particularly severe in areas where there was no engineered coastal protection, but while coastal revetments did provide enhanced protection for houses at the waterfront, the coastline in the study area appeared to be more heavily impacted than elsewhere in Langkawi because wave energy was focused on the area by offshore breakwaters built to protect the Langkawi port and airport. Emergency response after the tsunami was rapid and efficient but would have been improved if the local police station had not been rendered inoperative by the first wave, and if a mechanism had been in place to ensure that informal advance warnings transmitted between Phuket (Thailand), Langkawi and Penang (Malaysia) by tourist operators could have been more widely disseminated. [source] Organic litter: dominance over stones as a source of interrill flow roughness on low-gradient desert slopes at Fowlers Gap, arid western NSW, AustraliaEARTH SURFACE PROCESSES AND LANDFORMS, Issue 1 2003David Dunkerley Abstract Thirty-six runoff plot experiments provide data on flow depths, speeds, and Darcy,Weisbach friction coefficients (f) on bare soil surfaces, and surfaces to which were added sufficient extra plant litter or surface stones to provide projected cover of 5, 10 and 20 per cent. Precision flow depth data were derived with a computer-controlled gantry and needle gauge for two different discharges for each plot treatment. Taking a fixed flow intensity (Reynolds number, Re = 150) for purposes of comparison shows means of f = 17·7 for bare soil surfaces, f = 11·4 for added stone treatments, and f = 23·8 for added litter treatments. Many individual values of f for stone treatments are lower than for the bare soil surface, but all litter treatments show increases in fcompared to bare soil. The lowering of f in stone treatments relates to the submerged volume that the stones occupied, and the associated concentration of flow onto a smaller part of the plot surface. This leads to locally higher flow intensities and lower frictional drag along threads of flow that the obstacles create. Litter causes higher frictional drag because the particles are smaller, and, for the same cover fraction, are 100 times more numerous and provide 20 times the edge or perimeter length. Along these edges, which in total exceed 2·5 m g,1 (equivalent to 500 m m,2 for a loading of 2 t ha,1), surface tension draws up water from between the litter particles. This reduces flow depth there, and as a consequence of the lower flow intensity, frictional drag rises. Furthermore, no clear passage remains for the establishment of flow threads. These findings apply to shallow interrill flows in which litter is largely immobile. The key new result from these experiments is that under these conditions, a 20 per cent cover of organic litter can generate interrill frictional retardation that exceeds by nearly 41 per cent that of a bare soil surface, and twice that contributed by the same cover fraction of surface stones. Even greater dominance by litter can be anticipated at the many dryland sites where litter covers exceed those tested here. Copyright © 2002 John Wiley & Sons, Ltd. [source] Spring temperatures in the Sagehen Basin, Sierra Nevada, CA: implications for heat flow and groundwater circulationGEOFLUIDS (ELECTRONIC), Issue 3 2009MARIA BRUMM Abstract Heat flow in the Sierra Nevada, CA, is low despite its young geologic age. We investigate the possibility that advective heat transport by groundwater flow leads to an underestimate of heat flow in the Sierras based purely on borehole measurements. Using temperature and discharge measurements at springs in Sagehen Basin, we find that groundwater removes the equivalent of approximately 20,40 mW m,2 of geothermal heat from the basin. This is comparable with other heat flow measurements in the region and indicates that, in this basin, at least, groundwater does transport a significant amount of geothermal heat within the basin. Additionally, we use estimates of the mean residence time of water discharged at the springs along with hourly temperature records in springs to provide constraints on groundwater flow depths within the basin. An analytical model based on these constraints indicates that the heat removed by groundwater may represent 20% to >90% of the total heat flow in the basin. Without better constraints on the regional hydrogeology and the depth of circulation, we cannot determine whether the heat discharged at the springs represents a change in the mode of heat transfer, i.e. from conduction to advection at shallow depths (<100 m) or whether this is a component of heat transfer that should be added to measured conductive values. If the latter is true, and Sagehen Basin is representative of the Sierras, basal heat flow in the Sierra Nevada may be higher than previously thought. [source] |