Accuracy of the kinematic wave and diffusion wave approximations for time-independent flows with infiltration included

Errors in the kinematic wave and diffusion wave approximation for time-independent (or steady-state) cases of channel flow with infiltration were derived for three types of boundary conditions: zero flow at the upstream end, and critical flow depth and zero depth gradient at the downstream end. The diffusion wave approximation was found to be in excellent agreement with the dynamic wave approximation, with errors of less than 1·4% for KF²₀≥7·5, and up to 14% for KF²₀≤0·75 for the upstream boundary condition of zero discharge and finite depth, where K is the kinematic wave number and F₀ is the Froude number. The kinematic wave approximation was reasonably accurate except at the channel boundaries and for small values of KF²₀ (≤1). The accuracy of these approximations was significantly influenced by the downstream boundary, both in terms of the magnitude of the error and the segment of the channel reach for which these approximations would be applicable. © 1998 John Wiley & Sons, Ltd.     

Errors of kinematic wave and diffusion wave approximations for time‐independent flows with infiltration and momentum exchange included

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.     

ACCURACY OF KINEMATIC WAVE AND DIFFUSION WAVE APPROXIMATIONS FOR TIME-INDEPENDENT FLOW WITH MOMENTUM EXCHANGE INCLUDED

Errors in the kinematic wave and diffusion wave approximation for time-independent (or steady-state) cases of channel flow with momentum exchange included were derived for three types of boundary conditions: zero flow at the upstream end, and critical flow depth and zero depth gradient at the downstream end. The diffusion wave approximation was found to be in excellent agreement with the dynamic wave approximation, with errors of less than 1% for KF²₀≥7·5 and up to 12% for KF²₀≤0·75 for the upstream boundary condition of zero discharge and finite depth, where K is the kinematic wave number and F₀ is the Froude number. The kinematic wave approximation was reasonably accurate except at the channel boundaries and for small values of KF²₀ (≤1). The accuracy of these approximations was significantly influenced by the downstream boundary condition both in terms of the error magnitude and the segment of the channel reach for which these approximations would be applicable. © 1997 by John Wiley & Sons, Ltd.     

Accuracy of kinematic wave and diffusion wave approximations for time-independent flows

Errors in the kinematic wave and diffusion wave approximations for time-independent (or steady-state) cases of channel flow were derived for three types of boundary conditions: zero flow at the upstream end, and critical flow depth and zero depth gradient at the downstream end. The diffusion wave approximation was found to be in excellent agreement with the dynamic wave approximation, with errors in the range 1–2% for values of 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²₀ = 0.75), the errors were typically less than 15%. The accuracy of the diffusion wave approximation was greatly influenced by the downstream boundary condition. The error of the kinematic wave approximation was found to be less than 13% in the region 0.1 ? x ? 0.95 for KF = 7.5 and was greater than 30% for smaller values of KF (? 0.75). This error increased with strong downstream boundary control.     

A simple two‐parameter model for scaling hillslope surface runoff

The objective of this research was to develop and parameterise a physically justified yet low‐parameter model to quantify observed changes in surface runoff ratios with hillslope length. The approach starts with the assumption that a unit of rainfall‐excess runoff generated at a point is a fraction β of precipitation P (m) which travels some variable distance down a slope before reinfiltrating, depending on the local rainfall, climate, soils, etc. If this random distance travelled Y is represented by a distribution, then a survival function will describe the probability of this unit of runoff travelling further than some distance x (m). The total amount of per unit width runoff Q (m²) flowing across the lower boundary of a slope of length λ (m) may be considered the sum of all the proportions of the units of rainfall excess runoff integrated from the lower boundary x = 0 to the upper boundary x = λ of the slope. Using these assumptions we derive a model Q(λ) = βPμλ/(μ + λ), (μ > 0, 0 ≤ β ≤ 1, λ ≥ 0) that describes the change in surface runoff with slope length, where μ (m) is the mean of the random variable Y. Dividing both sides of this equation by Pλ yields a simple two‐parameter equation for the dimensionless hillslope runoff ratio Q_h(λ) = βμ/(μ + λ). The model was parameterised with new rainfall and runoff data collected from three replicates of bounded 2 m wide plots of four different lengths (0.5, 1.0, 2.0 and 4.0 m) for 2 years from a forested SE Australian site, and with 32 slope length–runoff data sets from 12 other published studies undertaken between 1934 and 2010. Using the parameterised model resulted in a Nash and Sutcliffe statistic between observed and predicted runoff ratio (for all data sets combined) of 0.93, compared with –2.1 when the runoff ratio was fixed at the value measured from the shortest plot. Copyright © 2014 John Wiley & Sons, Ltd.     

Application of kinematic wave theory for predicting flash flood hazards on coupled alluvial fan–piedmont plain landforms

A rainstorm that caused a severe flash flood on the piedmont plain at the toe positions of two alluvial fans located to the west of the Organ Mountains in Dona Ana County, New Mexico, USA, is analysed. The space–time distributions of rainfall are evaluated from the Next Generation Weather Radar (NEXRAD) and overland flow is modelled as kinematic wave. The spatial distribution of rainfall shows a topographic control. The greatest rainfall depth, duration, and intensity occurred at the higher elevation mountain slopes and decreased with decreasing elevation from the alluvial fans to the piedmont plain. The alluvial fan–piedmont plain system is modelled by coupling divergent and rectangular overland flow planes. Explicit finite difference approximations, hybridized with the analytical method of characteristics, are made to the kinematic wave equations to account for the spatial and temporal distribution of the rainfall and variable boundary conditions. Simulation results indicate that sheet‐flow floodwater elevations rise (1) in a nonlinear fashion from the apex to toe positions of the alluvial fans, and (2) near linearly from the toe positions of the alluvial fans onto the piedmont plains with the formation of kinematic shocks near the middle to the upstream end of the plane at times between the initiation of the rainstorm and the time of concentration of the plane. Thus, the maximum flooding occurs at the middle or upstream sections of the piedmont plains regardless of the pattern of space–time variability of rainfall. These results are in agreement with observed geomorphologic features suggesting that piedmont plains are naturally flood‐prone areas. This case study demonstrates that flood hazards on piedmont plains can exceed those on alluvial fans. The models presented in this study suggest that the flood hazard zones on coupled alluvial fan–piedmont plain landforms should be delineated transverse to the flow directions, as opposed to the flood hazard zones with boundaries in the longitudinal direction of the axis of an alluvial fan. Copyright © 2003 John Wiley & Sons, Ltd.     

Vortical and wave modes in 3D rotating stratified flows: random large-scale forcing

Utilizing an eigenfunction decomposition, we study the growth and spectra of energy in the vortical (geostrophic) and wave (ageostrophic) modes of a three-dimensional (3D) rotating stratified fluid as a function of ε = f/N, where f is the Coriolis parameter and N is the Brunt–Vaisala frequency. Throughout, we employ a random large-scale forcing in a unit aspect ratio domain and set these parameters such that the Froude and Rossby numbers are roughly comparable and much less than unity. Working in regimes characterized by moderate Burger numbers, i.e. Bu = 1/ε² < 1 or Bu ≥ 1, our results indicate profound change in the character of vortical and wave mode interactions with respect to Bu = 1. Indeed, previous analytical work concerning the qualitatively different nature of these interactions has been in limiting conditions of rotation or stratification domination (i.e. when Bu ? 1 or Bu ? 1, respectively). As with the reference state of ε = 1, for ε < 1 the wave mode energy saturates quite quickly and the ensuing forward cascade continues to act as an efficient means of dissipating ageostrophic energy. Further, these saturated spectra steepen as ε decreases: we see a shift from k ^?1 to k ^?5/3 scaling for k _f < k < k _d (where k _f and k _d are the forcing and dissipation scales, respectively). On the other hand, when ε > 1 the wave mode energy never saturates and comes to dominate the total energy in the system. In fact, in a sense the wave modes behave in an asymmetric manner about ε = 1. With regard to the vortical modes, for ε ≤ 1, the signatures of 3D quasigeostrophy are clearly evident. Specifically, we see a k ^?3 scaling for k _f < k < k _d and, in accord with an inverse transfer of energy, the vortical mode energy never saturates but rather increases for all k < k _f . In contrast, for ε > 1 and increasing, the vortical modes contain a progressively smaller fraction of the total energy indicating that the 3D quasigeostrophic subsystem, though always present, plays an energetically smaller role in the overall dynamics. Combining the vortical and wave modes, the total energy for k > k _f and ε ≤ 1 shows a transition as k increases wherein the vortical modes contain a large portion of the energy at large scales, while the wave modes dominate at smaller scales. There is no such transition when ε > 1 and the wave modes dominate the total energy for all k > k _f .     

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.     

Understanding changes in annual runoff following land use changes: a systematic data‐based approach

A simple modelling framework for assessing the response of ungauged catchments to land use change in South‐Western Australia is presented. The framework uses knowledge of transpiration losses from native vegetation and pasture and then partitions the ‘excess’ water (resulting from reduced transpiration after land use change) between runoff and deep storage. The simple partitioning is achieved by using soft information (satellite imagery, previous mapping and field assessment) to delimit the spread of the permanently saturated area close to the stream. Runoff is then assumed to increase in proportion to the saturated area, with the residual difference going to deep storage. The model parameters to describe the annual water yield are obtained a priori and no calibration to streamflow is required. We tested the model using gauged records over 25 years from paired catchment experiments in South‐Western Australia. Very good estimates of runoff were obtained from high rainfall (>1100 mm yr⁻¹) catchments (R² > 0·9) and for low rainfall (<900 mm yr⁻¹) catchments after clearing (R² = 0·96) but results were poorer (R² = 0·55) for an uncleared low rainfall catchment, although overall balances were reasonable. In the drier uncleared catchments, the within‐year distributions of rainfall may exert a substantial influence on runoff response that is not completely captured by the presented model. Copyright © 2005 John Wiley & Sons, Ltd.     

Flood routing based on diffusion wave equation using mixing cell method

A one-dimensional non-linear diffusion wave equation is derived from the Saint Venant equations with neglect of the inertia terms. This non-linear equation has no general analytical solution. Numerical schemes are therefore employed to discretize the space and time axes and convert the differential equation to difference form. In this study, the mixing cell method is used to convert the diffusion wave equation to difference form, in which the difference term can be eliminated by selecting an optimal space step size Δx when time step size Δt is given. When the time step size Δt→0, the space step size Δx=Q/(2S₀BC]_k) where Q is discharge, S₀ is bed slope, B is channel width and C_k is kinematic wave celerity, which is the same as the characteristic length proposed by Kalinin and Milyukov. The results of application to two cases show that the mixing cell and linear channel flow routing methods produce hydrographs that are in agreement with the observed flood hydrographs. © 1997 John Wiley & Sons, Ltd.     

Analytical solution of transient flow in a sloping soil layer with recharge

Abstract

An analytical solution of planar flow in a sloping soil layer described by the linearized extended Boussinesq equation is presented. The solution consists of the sum of steady-state and transient-series solutions, the latter in a separation-of-variables form, and can satisfy an arbitrary initial condition via collocation; this feature reduces the number of series terms, making the solution efficient. Key parameter is the dimensionless linearization depth η _o (R), R being the dimensionless recharge. The variable η _o (R), not the slope, characterizes the flow as kinematic or diffusive, and R ≈ 0.2 demarcates the two regimes. The transient series converges rapidly for large η _o (large R, near-diffusive flow) and slowly as η _o → 0 (kinematic flow). The quasi-steady (QS) state method of Verhoest & Troch is also analysed and it is shown that the QS depth profiles approximate the transient ones well, only if Δt exceeds a system-dependent transition time between flow states (possibly >>1 day). In an application example for a 30-day recharge series, the QS solution differs from the transient one by as much as 20% (RMSE = 15%), does not track recharge changes as well and fails to conserve mass.     

Coupling the Xinanjiang model to a kinematic flow model based on digital drainage networks for flood forecasting

The Xinanjiang model, which is a conceptual rainfall‐runoff model and has been successfully and widely applied in humid and semi‐humid regions in China, is coupled by the physically based kinematic wave method based on a digital drainage network. The kinematic wave Xinanjiang model (KWXAJ) uses topography and land use data to simulate runoff and overland flow routing. For the modelling, the catchment is subdivided into numerous hillslopes and consists of a raster grid of flow vectors that define the water flow directions. The Xinanjiang model simulates the runoff yield in each grid cell, and the kinematic wave approach is then applied to a ranked raster network. The grid‐based rainfall‐runoff model was applied to simulate basin‐scale water discharge from an 805‐km² catchment of the Huaihe River, China. Rainfall and discharge records were available for the years 1984, 1985, 1987, 1998 and 1999. Eight flood events were used to calibrate the model's parameters and three other flood events were used to validate the grid‐based rainfall‐runoff model. A Manning's roughness via a linear flood depth relationship was suggested in this paper for improving flood forecasting. The calibration and validation results show that this model works well. A sensitivity analysis was further performed to evaluate the variation of topography (hillslopes) and land use parameters on catchment discharge. Copyright © 2009 John Wiley & Sons, Ltd.     

A laboratory study of formative conditions for characteristic ripple patterns associated with a change in wave conditions

A wave‐?ume experiment was conducted to examine the formative condition for three types of distinctive bedforms that emerged through deformation of existing ripples due to waning wave power. They were ripple marks with: (1) a single secondary crest, (2) double secondary crests, and (3) a rounded crest. Data were analysed using two parameters, kh and d₀/λ_*, where k is the wave number, h is the water depth, d₀ is the near‐bottom orbital diameter, and λ_* is the spacing of existing ripples. The former quantity, kh, was employed as a surrogate for the degree of ?ow asymmetry. The result showed that ripples with secondary crests developed under a rather symmetrical ?ow ?eld with kh ? 0·7, if d₀/λ_* ? 1·2, whereas rounded‐crest ripples emerged under asymmetrical ?ow ?eld with kh ? 0·7, if d₀/λ_* ? ?2·9 kh + 3·2. The number of secondary crests, which initially occurred in each trough, was single if d₀/λ_* ? 0·8, or double if d₀/λ_* ? 0·8. Copyright © 2004 John Wiley & Sons, Ltd.     

Influence of upstream inflow on wave celerity and time to equilibrium on an overland plane

Abstract

Based on the kinematic wave equations, formulae for the wave celerity along an overland plane subject to uniform rainfall excess and with a constant upstream inflow together with the corresponding average wave celerity and time to equilibrium for the entire plane are derived. The formulae are further developed in terms of both the Darcy-Weisbach resistance coefficient and the Manning resistance coefficient. By comparing the wave celerities, the average wave celerities and the time to equilibrium for planes with and without upstream inflow show that the upstream inflow causes the wave celerity and the average wave celerity to be faster and the times to equilibrium to be shorter. The effect of upstream inflow is greater with increasing inflow, but the marginal effect decreases with increasing inflow. The effect is greatest for laminar flow and least for turbulent flow. For the wave celerity, the effect is also greatest at the upstream end of the plane and least at the downstream end of the plane.     

Experimental verification of the fracture density and shear‐wave splitting relationship using synthetic silica cemented sandstones with a controlled fracture geometry

We present laboratory ultrasonic measurements of shear‐wave splitting from two synthetic silica cemented sandstones. The manufacturing process, which enabled silica cementation of quartz sand grains, was found to produce realistic sandstones of average porosity 29.7 ± 0.5% and average permeability 29.4 ± 11.3 mD. One sample was made with a regular distribution of aligned, penny‐shaped voids to simulate meso‐scale fractures in reservoir rocks, while the other was left blank. Ultrasonic shear waves were measured with a propagation direction of 90° to the coincident bedding plane and fracture normal. In the water saturated blank sample, shear‐wave splitting, the percentage velocity difference between the fast and slow shear waves, of <0.5% was measured due to the bedding planes (or layering) introduced during sample preparation. In the fractured sample, shear‐wave splitting (corrected for layering anisotropy) of 2.72 ± 0.58% for water, 2.80 ± 0.58% for air and 3.21 ± 0.58% for glycerin saturation at a net pressure of 40 MPa was measured. Analysis of X‐ray CT scan images was used to determine a fracture density of 0.0298 ± 0.077 in the fractured sample. This supports theoretical predictions that shear‐wave splitting (SWS) can be used as a good estimate for fracture density in porous rocks (i.e., SWS = 100ε_f, where ε_f is fracture density) regardless of pore fluid type, for wave propagation at 90° to the fracture normal.     

An enhanced rainfall–runoff model with coupled canopy interception

The translation of rainfall to runoff is significantly affected by canopy interception. Therefore, a realistic representation of the role played by vegetation cover when modelling the rainfall–runoff system is essential for predicting water, sediment, and nutrient transport on hillslopes. Here, we developed a new mathematical model to describe the dynamics of interception, infiltration, and overland flow on canopy-covered sloping land. Based on the relationship between rainfall intensity and the maximum interception rate, the interception process was modelled under two simplified scenarios (i.e., r_e ≤ Int_m and r_e > Int_m). Parameterization of the model was based on consideration of both vegetation condition and soil properties. By analysing the given examples, we found that Int_m reflects the capacity of the canopy to store the precipitation, k reveals the ability of the canopy to retain the intercepted water, and the processes of infiltration and runoff generation are impacted dramatically by Int_m and k. To evaluate the model, simulated rainfall experiments were conducted in 2 years (2016 and 2017) across six cultivation plots at Changwu State Key Agro-Ecological Experimental Station of the Chinese Loess Plateau. The parameters were obtained by fitting the unit discharge (simulated rainfall experiments in 2016) using the least squares method, and estimation formulas for parameters pertaining to vegetation/soil factors (measured in 2016) were constructed via multiple nonlinear regressions. By matching the simulated results and unit discharge (simulated rainfall experiments in 2017), the validity of the model was verified, and a reasonable precision (average R² = .86 and average root mean square error = 6.45) was obtained. The model developed in this research creatively incorporates the canopy interception process to complement the modelling of rainfall infiltration and runoff generation during vegetation growth and offers an improved hydrological basis to analyse matter transport during rainfall events.     

Lateral flow thresholds for aspen forested hillslopes on the Western Boreal Plain,Alberta, Canada

To predict the long‐term sustainability of water resources on the Boreal Plain region of northern Alberta, it is critical to understand when hillslopes generate runoff and connect with surface waters. The sub‐humid climate (P ≤ ET) and deep glacial sediments of this region result in large available soil storage capacity relative to moisture surpluses or deficits, leading to threshold‐dependent rainfall‐runoff relationships. Rainfall simulation experiments were conducted using large magnitude and high intensity applications to examine the thresholds in precipitation and soil moisture that are necessary to generate lateral flow from hillslope runoff plots representative of Luvisolic soils and an aspen canopy. Two adjacent plots (areas of 2·95 and 3·4 m²) of contrasting antecedent moisture conditions were examined; one had tree root uptake excluded for two months to increase soil moisture content, while the second plot allowed tree uptake over the growing season resulting in drier soils. Vertical flow as drainage and soil moisture storage dominated the water balances of both plots. Greater lateral flow occurred from the plot with higher antecedent moisture content. Results indicate that a minimum of 15–20 mm of rainfall is required to generate lateral flow, and only after the soils have been wetted to a depth of 0·75 m (C‐horizon). The depth and intensity of rainfall events that generated runoff > 1 mm have return periods of 25 years or greater and, when combined with the need for wet antecendent conditions, indicate that lateral flow generation on these hillslopes will occur infrequently. Copyright © 2008 John Wiley & Sons, Ltd.     

An experimental study on the evolution law of roll wave parameters in overland flow

Roll waves commonly occur in overland flow and have an important influence on the progress of soil erosion on slopes. This study aimed to explore the evolution and mechanism of roll waves on steep slopes. The potential effects of flow rate, rainfall intensity and bed roughness on the laws controlling roll wave parameters were investigated. The flow rates, rainfall intensities and bed roughness varied from 5 to 30 L/min, 0 to 150 mm/h, and 0.061 to 1.700 mm, respectively. The results indicate that roll waves polymerize significantly along the propagation path, and bed roughness and rainfall affect the generation and evolution of roll waves. The wave velocity, length and height decreased with bed roughness, whereas the wave frequency increased with increasing bed roughness under fixed flow rate and rainfall intensity conditions. Rainfall increased the wave velocity and wavelength and decreased the wave frequency. The wave velocity, height and wavelength tended to increase with an increasing flow rate. Rainfall promoted the generation of roll waves, whereas bed roughness had the opposite effect. The generation of roll waves is closely related to the Froude number (Fr) and flow resistance. In this experiment, the range of the Reynolds number for the roll waves generated in the laminar region was 142–416, and the range of the flow resistance coefficient was 0.64–4.85. The critical value of the Fr for flow instability in the laminar region was approximately 0.57. Exploring the generation and evolution law of roll waves is necessary for understanding the processes and dynamic mechanisms of slope soil erosion.