Observations and simulations of an unsteady island wake in the Firth of Forth, Scotland

Beamer Rock, a 50-m-wide island in the Firth of Forth, produces a distinctive von Kármán vortex street wake, the characteristics of which depend on the speed and direction of the tidal flow. ADCP (acoustic Doppler current profiler), CTD (conductivity–temperature–depth) and aerial photograph data were collected from the region during flood and ebb tidal flow under neap and spring conditions. Good agreement was found between the observations and the results from a fixed-grid depth-averaged numerical tidal model. The island wake parameter correctly predicted the unsteady nature of the wake, and the Strouhal number (defined in terms of flow past a circular cylinder) was found to give excellent predictions of the wake wavelength when scaled on the island width. Contrary to published results, the study shows that it is possible to accurately simulate an unsteady island wake using a relatively coarse fixed-grid numerical model.Responsible Editor: Jens Kappenberg     

第二松花江流域水电厂综合中长期水文预报的分析研究

基于对第二松花江流域上游小山、松山、两江水电站中长期水文预报研究成果，论述了综合中长期水文预报的研究思路和定性预报、定量预报的分析研究方法。该研究通过2002年实践检验，具有较高的预报精度。     

A Finite Volume Scheme for the Transport of Radionucleides in Porous Media

The COUPLEX1 Test case (Bourgeat et al., 2003) is devoted to the comparison of numerical schemes on a convection–diffusion–reaction problem. We first show that the results of the simulation can be mainly predicted by a simple analysis of the data. A finite volume scheme, with three different treatments of the convective term, is then shown to deliver accurate and stable results under a low computational cost.     

Turbulent flow over a dune: Green River,Colorado

Detailed echo‐sounder and acoustic Doppler velocimeter measurements are used to assess the temporal and spatial structure of turbulent flow over a mobile dune in a wide, low‐gradient, alluvial reach of the Green River. Based on the geometric position of the sensor over the bedforms, measurements were taken in the wake, in transitional flow at the bedform crest, and in the internal boundary layer. Spatial distributions of Reynolds shear stress, turbulent kinetic energy, turbulence intensity, and correlation coefficient are qualitatively consistent with those over fixed, two‐dimensional bedforms in laboratory flows. Spectral and cospectral analysis demonstrates that energy levels in the lee of the crest (i.e. wake) are two to four times greater than over the crest itself, with minima over the stoss slope (within the developing internal boundary layer). The frequency structure in the wake is sharply defined with single, dominant peaks. Peak and total spectral and cross‐spectral energies vary over the bedform in a manner consistent with wave‐like perturbations that ‘break’ or ‘roll up’ into vortices that amalgamate, grow in size, and eventually diffuse as they are advected downstream. Fluid oscillations in the lee of the dune demonstrate Strouhal similarity between laboratory and field environments, and correspondence between the peak frequencies of these oscillations and the periodicity of surface boils was observed in the field. Copyright © 2005 John Wiley & Sons, Ltd.     

The effective number of cycles of earthquake ground motion

The seismic response of any system that accumulates damage under cyclic loading is dependent not only on the maximum amplitude of the motion but also its duration. This is explicitly recognized in methods for estimating the liquefaction potential of soil deposits. Many researchers have proposed that the effective number of cycles of the ground motion is a more robust indicator of the destructive capacity of the shaking than the duration. However, as is the case with strong‐motion duration, there is no universally accepted approach to determining the effective number of cycles of motion, and the different methods that have been proposed can give widely varying results for a particular accelerogram. Definitions of the effective number of cycles of motion are reviewed, classified and compared. Measurements are found to differ particularly for accelerograms with broad‐banded frequency content, which contain a significant number of non‐zero crossing peaks. The key seismological parameters influencing the number of cycles of motion and associated equations for predicting this quantity for future earthquakes are identified. Correlations between cycle counts and different duration measures are explored and found to be rather poor in the absence of additional parameters. Copyright © 2004 John Wiley & Sons, Ltd.     

An alternative Boussinesq equation considering the effect of hysteresis on coastal groundwater waves

This paper presents an alternative Boussinesq equation considering hysteresis effect via a third‐order derivative term. By introducing an improved moisture–pressure retention function, this equation describes, with reasonable precision, groundwater propagation in coastal aquifers subject to Dirichlet boundary condition of different oscillation frequencies. Test results confirmed that it is necessary to consider horizontal and vertical flows in unsaturated zone, because of their variable influences on hysteresis. Hysteresis in unsaturated zone can affect the water table wave number of groundwater wave motion, such as wave damping rate and phase lag. Oscillations with different periods exert different hysteresis effect on wave propagation. Truncation/shrinkage of unsaturated zones also affects the strength of hysteresis. These impacts can be reflected in the alternative Boussinesq equation by adjusting the parameter representing the variation rate of moisture associated with pressure change, as opposed to traditional computationally expensive hysteresis algorithms. The present Boussinesq equation is simple to use and can provide feasible basis for future coupling of groundwater and surface water models. Copyright © 2016 John Wiley & Sons, Ltd.     

Glacial flood pulse effects on benthic fauna in equatorial high‐Andean streams

Equatorial glacier‐fed streams present unique hydraulic patterns when compared to glacier‐fed observed in temperate regions as the main variability in discharge occurs on a daily basis. To assess how benthic fauna respond to these specific hydraulic conditions, we investigated the relationships between flow regime, hydraulic conditions (boundary Reynolds number, Re*), and macroinvertebrate communities (taxon richness and abundance) in a tropical glacier‐fed stream located in the high Ecuadorian Andes (> 4000 m). Both physical and biotic variables were measured under four discharge conditions (base‐flow and glacial flood pulses of various intensities), at 30 random points, in two sites whose hydraulic conditions were representative to those found in other streams of the study catchment. While daily glacial flood pulses significantly increased hydraulic stress in the benthic habitats (appearance of Re* > 2000), low stress areas still persisted even during extreme flood events (Re* < 500). In contrast to previous research in temperate glacier‐fed streams, taxon richness and abundance were not significantly affected by changes in hydraulic conditions induced by daily glacial flood pulses. However, we found that a few rare taxa, in particular rare ones, preferentially occurred in highly stressed hydraulic habitats. Monte‐Carlo simulations of benthic communities under glacial flood reduction scenarios predicted that taxon richness would be significantly reduced by the loss of high hydraulic stress habitats following glacier shrinking. This pioneer study on the relationship between hydraulic conditions and benthic diversity in an equatorial glacial stream evidenced unknown effects of climate change on singular yet endangered aquatic systems. Copyright © 2013 John Wiley & Sons, Ltd.     

The effect of expansion ratio on the critical Reynolds number in single fracture flow with sudden expansion

Flow in a single fracture (SF) is an important research subject in groundwater hydrology, hydraulic engineering, radioactive nuclear waste repository and geotechnical engineering. An abruptly changing aperture is a unique type of SF. This study discusses the relation between the values of the critical Reynolds number (Re_c) for the onset of symmetry breaking of flow and the expansion ratio (E) of SF, which is defined as the ratio between the outlet (D) and inlet (d) apertures. This study also investigates the effect of inlet aperture d on Re_c for flow in an SF with abruptly changing apertures (SF‐ACA) using the finite volume method. Earlier numerical and experimental results showed that flow is symmetric in respect to the central plane of the SF‐ACA at small Reynolds number (Re) but becomes asymmetric when Re is sufficiently large. Our simulations show that the value of Re_c decreases with the increasing E, and the relationship between the logarithm of Re_c and E can be described accurately using either a quadratic polynomial function or a logarithmic function. However, the relationship of Re_c and d for a given E value is vague, and Re_c becomes even less sensitive to d when E increases. This study also reveals that the hydraulic gradient (J) and flow velocity (v) follow a super‐linear relationship that can be fitted almost perfectly by the Forchheimer equation. The inertial component (J_i) of J increases monotonically with Re, whereas the viscous component (J_v) of J decreases monotonically with Re. The Re value corresponding to equal inertial and viscous components of J (named as the transitional point Re) decreases when E increases, and such a transitional point Re should be closely related to the critical Reynolds number Re_c, although a rigorous theoretical proof is not yet available. Copyright © 2015 John Wiley & Sons, Ltd.     

Revised runoff curve number for runoff prediction in the Loess Plateau of China

The soil conservation service (now Natural Resources Conservation Service) Curve Number (SCS-CN), one of the most commonly used methods for surface runoff prediction. The runoff calculated by this method was very sensitive to CN values. In this study, CN values were calculated by both arithmetic mean (CN_C) and least square fit method (CN_F) using observed rainfall-runoff data from 43 sites in the Loess Plateau region, which are considerably different from the CN₂ values obtained from the USDA-SCS handbook table (CN_T). The results showed that using CN_C instead of CN_T for each watershed produce little improvement, while replacing CN_T with CN_F improves the performance of the original SCS-CN method, but still performs poorly in most study sites. This is mainly due to the SCS-CN method using a constant CN value and discounting of the temporal variation in rainfall-runoff process. Therefore, three factors—soil moisture, rainfall depth and intensity—affecting the surface runoff variability are considered to reflect the variation of CN in each watershed, and a new CN value was developed. The reliability of the proposed method was tested with data from 38 watersheds, and then applied to the remaining five typical watersheds using the optimized parameters. The results indicated that the proposed method, which boosted the model efficiencies to 81.83% and 74.23% during calibration and validation cases, respectively, performed better than the original SCS-CN and the Shi and Wang (2020b) method, a modified SCS-CN method based on tabulated CN value. Thus, the proposed method incorporating the influence of the temporal variability of soil moisture, rainfall depth, and intensity factors suggests an accurate runoff prediction for general applications under different hydrological and climatic conditions on the Loess Plateau region.