Development of the instantaneous unit hydrograph using stochastic differential equations

Recognizing that simple watershed conceptual models such as the Nash cascade ofn equal linear reservoirs continue to be reasonable means to approximate the Instantaneous Unit Hydrograph (IUH), it is natural to accept that random errors generated by climatological variability of data used in fitting an imprecise conceptual model will produce an IUH which is random itself. It is desirable to define the random properties of the IUH in a watershed in order to have a more realistic hydrologic application of this important function. Since in this case the IUH results from a series of differential equations where one or more of the uncertain parameters is treated in stochastic terms, then the statistical properties of the IUH are best described by the solution of the corresponding Stochastic Differential Equations (SDE's). This article attempts to present a methodology to derive the IUH in a small watershed by combining a classical conceptual model with the theory of SDE's. The procedure is illustrated with the application to the Middle Thames River, Ontario, Canada, and the model is verified by the comparison of the simulated statistical measures of the IUH with the corresponding observed ones with good agreement.     

VORTEX FLOW FIELD IN A SCOUR HOLE AROUND ABUTMENTS

The three-dimensional flow field in a scour hole around different abutments under a clear water regime was experimentally measured in a laboratory flume, using the Acoustic Doppler Velocimeter (ADV). Three types of abutments used in the experiments were vertical-wall (rectangular section), 45° wing-wall (45° polygonal section) and semicircular. The three-dimensional time-averaged velocity components were detected at different vertical planes for vertical-wall abutment and azimuthal planes for wing-wall and semicircular abutments. The velocity components were also measured at different horizontal planes. In the upstream, presentations of flow field through vector plots at vertical / azimuthal and horizontal planes show the existence of a primary vortex associated with the downflow inside the scour hole. On the other hand, in the downstream, the flow field shows a reversed flow near the abutments having a subsequent recovery with a passage of flow as a part of the main flow. The data presented in this p     

OUTBURST FLOODS FROM GLACIER-DAMMED LAKES: THE EFFECT OF MODE OF LAKE DRAINAGE ON FLOOD MAGNITUDE

Published accounts of outburst floods from glacier-dammed lakes show that a significant number of such floods are associated not with drainage through a tunnel incised into the basal ice—the process generally assumed—but rather with ice-marginal drainage, mechanical failure of part of the ice dam, or both. Non-tunnel floods are strongly correlated with formation of an ice dam by a glacier advancing from a tributary drainage into either a main river valley or a pre-existing body of water (lake or fiord). For a given lake volume, non-tunnel floods tend to have significantly higher peak discharges than tunnel-drainage floods. Statistical analysis of data for floods associated with subglacial tunnels yields the following empirical relation between lake volume V and peak discharge Q_p : Q_p = 46V^0.66 (r² = 0.70), when Q_p is expressed in metres per second and V in millions of cubic metres. This updates the so-called Clague–Mathews relation. For non-tunnel floods, the analogous relation is Q_p = 1100V^0.44 (r² = 0.58). The latter relation is close to one found by Costa (1988) for failure of constructed earthen dams. This closeness is probably not coincidental but rather reflects similarities in modes of dam failure and lake drainage. We develop a simple physical model of the breach-widening process for non-tunnel floods, assuming that (1) the rate of breach widening is controlled by melting of the ice, (2) outflow from the lake is regulated by the hydraulic condition of critical flow where water enters the breach, and (3) the effect of lake temperature may be dealt with as done by Clarke (1982). Calculations based on the model simulate quite well outbursts from Lake George, Alaska. Dimensional analysis leads to two approximations of the form Q_p ∝ V^qf(h_i, θ₀), where q = 0.5 to 0.6, h_i is initial lake depth, θ₀ is lake temperature, and the form of f (h_i, θ₀) depends on the relative importance of viscous dissipation and the lake's thermal energy in determining the rate of breach opening. These expressions, along with the regression relations, should prove useful for assessing the probable magnitude of breach-type outburst floods.     

Spatial rainfall model using a pattern classifier for estimating missing daily rainfall data

Missing data in daily rainfall records are very common in water engineering practice. However, they must be replaced by proper estimates to be reliably used in hydrologic models. Presented herein is an effort to develop a new spatial daily rainfall model that is specifically intended to fill in gaps in a daily rainfall dataset. The proposed model is different from a convectional daily rainfall generation scheme in that it takes advantage of concurrent measurements at the nearby sites to increase the accuracy of estimation. The model is based on a two-step approach to handle the occurrence and the amount of daily rainfalls separately. This study tested four neural network classifiers for a rainfall occurrence processor, and two regression techniques for a rainfall amount processor. The test results revealed that a probabilistic neural network approach is preferred for determining the occurrence of daily rainfalls, and a stepwise regression with a log-transformation is recommended for estimating daily rainfall amounts.     

Analytical investigation of the exact groundwater divide between rivers beyond the Dupuit–Forchheimer approximation

The groundwater divide is a key feature of river basins and significantly influenced by subsurface hydrological processes. For an unconfined aquifer between two parallel rivers or ditches, it has long been defined as the top of the water table based on the Dupuit–Forchheimer approximation. However, the exact groundwater divide is subject to the interface between two local flow systems transporting groundwater to rivers from the infiltration recharge. This study contributes a new analytical model for two-dimensional groundwater flow between rivers of different water levels. The flownet is delineated in the model to identify groundwater flow systems and the exact groundwater divide. Formulas with two dimensionless parameters are derived to determine the distributed hydraulic head, the top of the water table and the groundwater divide. The locations of the groundwater divide and the top of the water table are not the same. The distance between them in horizontal can reach up to 8.9% of the distance between rivers. Numerical verifications indicate that simplifications in the analytical model do not significantly cause misestimates in the location of the groundwater divide. In contrast, the Dupuit–Forchheimer approximation yields an incorrect water table shape. The new analytical model is applied to investigate groundwater divides in the Loess Plateau, China, with a Monte Carlo simulation process taking into account the uncertainties in the parameters.     

A new multiscale scheme for computing statistical moments in single phase flow in heterogeneous porous media

In this work we develop a new multiscale procedure to compute numerically the statistical moments of the stochastic variables which govern single phase flow in heterogeneous porous media. The technique explores the properties of the log-normally distributed hydraulic conductivity, characterized by power-law or exponential covariances, which shows invariance in its statistical structure upon a simultaneous change of the scale of observation and strength of heterogeneity. We construct a family of equivalent stochastic hydrodynamic variables satisfying the same flow equations at different scales and strengths of heterogeneity or correlation lengths. Within the new procedure the governing equations are solved in a scaled geology and the numerical results are mapped onto the original medium at coarser scales by a straightforward rescaling. The new procedure is implemented numerically within the Monte Carlo algorithm and also in conjunction with the discretization of the low-order effective equations derived from perturbation analysis. Numerical results obtained by the finite element method show the accuracy of the new procedure to approximated the two first moments of the pressure and velocity along with its potential in reducing drastically the computational cost involved in the numerical modeling of both power-law and exponential covariance functions.     

Stochastic dynamic stiffness of a surface footing for offshore wind turbines: Implementing a subset simulation method to estimate rare events

The purpose of this study, which concerns the stochastic dynamic stiffness of foundations for large offshore wind turbines, is to quantify uncertainties related to the first natural frequency of a turbine supported by a surface footing and to estimate the low event probabilities. Herein, a simple model of a wind turbine structure with equivalent coupled springs at the base is calibrated with the mean soil property values. A semianalytical solution, based on the Green׳s function for a layered half-space is utilized for estimation of foundation responses. Soil elastic modulus and layer depth are considered as random variables with lognormal distributions. The uncertainties are quantified, and the estimation of rare events of the first natural frequency is discussed through an advanced reliability approach based on subset simulation. This analysis represents a first step in the estimation of the safety with respect to the failure of a turbine in the fatigue limit state.     

A new approach to rapid,desktop-level,environmental flow assessments for rivers in South Africa

Abstract

An approach is presented for desktop-level environmental flow requirement (EFR) determination that is aligned with the Habitat Flow–Stressor Response (HFSR) method which evolved in South Africa over recent years. The HFSR method integrates hydrological, hydraulic and ecological habitat data, involves ecological and hydraulic specialists and is data-intensive and time-consuming. The revised desktop method integrates hydrological information with estimates of channel hydraulic cross-sectional characteristics to generate habitat-type frequencies under changing flow conditions. This information is used with the expected natural habitat requirements to determine acceptable habitat availability under different levels of ecological protection, which is then used with the hydraulic data to define flow regime characteristics that meet the ecological objectives. The paper describes the model components, discusses the assumptions, data requirements and limitations and presents some example results. The revised desktop approach uses approaches that are aligned with the more complex methods and generates results that are similar.

Editor D. Koutsoyiannis; Guest editor M. Acreman

Citation Hughes, D.A., Desai, A.Y., Birkhead, A.L., and Louw, D., 2014. A new approach to rapid, desktop-level, environmental flow assessments for rivers in South Africa. Hydrological Sciences Journal, 59 (3–4), 673–687.     

Stochastic modelling of unresolved eddy fluxes

A subgrid-scale parameterization scheme motivated by statistical closure theory, but employing statistics obtained from high-resolution direct numerical simulations, is applied to large eddy simulations of two-level quasigeostrophic turbulence on the sphere. It is shown that these parameterizations are consistent with the phenomenology of quasigeostrophic turbulence. The parameterizations consist of 2 × 2 dissipation and stochastic forcing covariance matrices at each wavenumber, with the off-diagonal elements of the matrices representing vertical mixing. Two flow regimes, characterized by their deformation scales, are considered, namely atmospheric and oceanic. In the former, the deformation scale is fully resolved, and the truncation scale is within the enstrophy cascading interial range. In the latter, the deformation scale is not fully resolved, and the truncation scale is within the energy cascading inertial range. It is demonstrated through numerical experiments that both stochastic and deterministic variants of the scheme give comparable results for the energy spectra in the atmospheric regime. In the oceanic regime, the stochastic variant again gives excellent results, but the deterministic variant is found to be numerically unstable.