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.     

Modelling two‐dimensional steady‐state groundwater flow and flow sensitivity to boundary conditions in blanket peat complexes

This study used a two‐dimensional steady‐state finite‐element groundwater flow model to simulate groundwater flow in two Newfoundland blanket peat complexes and to examine flow system sensitivity to changes in water table recharge and aquifer properties. The modelling results were examined within the context of peat‐forming processes in the two complexes. Modelled flow compared favourably with observed flow. The sensitivity analyses suggested that more highly decomposed bog peat along bog margins probably has/had a positive impact on net peat accumulation within bog interiors. Peat with lower hydraulic conductivity along bog margins effectively impedes lateral drainage, localizes water table drawdown to extreme bog margins, and elevates water tables along bog interiors. Peat formation and elevated water tables in adjacent poor fens/laggs currently rely on placic and ortstein horizons impeding vertical drainage and water flow inputs from adjacent bogs. Modest reductions in atmospheric recharge were found to govern bog‐flow‐system geometries in a way that would adversely affect paludification processes in adjacent fens/laggs. Copyright © 2004 John Wiley & Sons, Ltd.     

Do not only connect: a model of infiltration‐excess overland flow based on simulation

The paper focusses on connectivity in the context of infiltration‐excess overland flow and its integrated response as slope‐base overland flow hydrographs. Overland flow is simulated on a sloping surface with some minor topographic expression and spatially differing infiltration rates. In each cell of a 128 × 128 grid, water from upslope is combined with incident rainfall to generate local overland flow, which is stochastically routed downslope, partitioning the flow between downslope neighbours. Simulations show the evolution of connectivity during simple storms. As a first approximation, total storm runoff is similar everywhere, discharge increasing proportionally with drainage area. Moderate differences in plan topography appear to have only a second‐order impact on hydrograph form and runoff amount. Total storm response is expressed as total runoff, runoff coefficient or total volume infiltrated; each plotted against total storm rainfall, and allowing variations in average gradient, overland flow roughness, infiltration rate and storm duration. A one‐parameter algebraic expression is proposed that fits simulation results for total runoff, has appropriate asymptotic behaviour and responds rationally to the variables tested. Slope length is seen to influence connectivity, expressed as a scale distance that increases with storm magnitude and can be explicitly incorporated into the expression to indicate runoff response to simple events as a function of storm size, storm duration, slope length and gradient. The model has also been applied to a 10‐year rainfall record, using both hourly and daily time steps, and the implications explored for coarser scale models. Initial trails incorporating erosion continuously update topography and suggest that successive storms produce an initial increase in erosion as rilling develops, while runoff totals are only slightly modified. Other factors not yet considered include the dynamics of soil crusting and vegetation growth. Copyright © 2014 John Wiley & Sons, Ltd.     

Poro‐acoustoelasticity of fluid‐saturated rocks

We generalize the classical theory of acoustoelasticity to the porous case (one fluid and a solid frame) and finite deformations. A unified treatment of non‐linear acoustoelasticity of finite strains in fluid‐saturated porous rocks is developed on the basis of Biot’s theory. A strain‐energy function, formed with eleven terms, combined with Biot’s kinetic and dissipation energies, yields Lagrange’s equations and consequently the wave equation of the medium. The velocities and dissipation factors of the P‐ and S‐waves are obtained as a function of the 2nd‐ and 3rd‐order elastic constants for hydrostatic and uniaxial loading. The theory yields the limit to the classical theory if the fluid is replaced with a solid with the same properties of the frame. We consider sandstone and obtain results for open‐pore jacketed and closed‐pore jacketed ‘gedanken’ experiments. Finally, we compare the theoretical results with experimental data.     

Back‐projection stacking of P‐ and S‐waves to determine location and focal mechanism of microseismic events recorded by a surface array

We present an automatic method of processing microseismic data acquired at the surface by a star‐like array. The back‐projection approach allows successive determination of the hypocenter position of each event and of its focal mechanisms. One‐component vertical geophone groups and three‐component accelerometers are employed to monitor both P‐ and S‐waves. Hypocenter coordinates are determined in a grid by back‐projection stacking of the short‐time‐average‐to‐long‐time‐average ratio of absolute amplitudes at vertical components and polarization norm derived from horizontal components of the P‐ and S‐waves, respectively. To make the location process more efficient, calculation is started with a coarse grid and zoomed to the optimum hypocenter using an oct‐tree algorithm. The focal mechanism is then determined by stacking the vertical component seismograms corrected for the theoretical P‐wave polarity of the focal mechanism. The mechanism is resolved in the coordinate space of strike, dip, and rake angles. The method is tested on 34 selected events of a dataset of hydraulic fracture monitoring of a shale gas play in North America. It was found that, by including S‐waves, the vertical accuracy of locations improved by a factor of two and is equal to approximately the horizontal location error. A twofold enhancement of horizontal location accuracy is achieved if a denser array of geophone groups is used instead of the sparse array of three‐component seismometers. The determined focal mechanisms are similar to those obtained by other methods applied to the same dataset.     

Modelling the effects of permafrost loss on discharge from a wetland‐dominated,discontinuous permafrost basin

Permafrost degradation in the peat‐rich southern fringe of the discontinuous permafrost zone is catalysing substantial changes to land cover with expansion of permafrost‐free wetlands (bogs and fens) and shrinkage of forest‐dominated permafrost peat plateaux. Predicting discharge from headwater basins in this region depends upon understanding and numerically representing the interactions between storage and discharge within and between the major land cover types and how these interactions are changing. To better understand the implications of advanced permafrost thaw‐induced land cover change on wetland discharge, with all landscape features capable of contributing to drainage networks, the hydrological behaviour of a channel fen sub‐basin in the headwaters of Scotty Creek, Northwest Territories, Canada, dominated by peat plateau–bog complexes, was modelled using the Cold Regions Hydrological Modelling platform for the period of 2009 to 2015. The model construction was based on field water balance observations, and performance was deemed adequate when evaluated against measured water balance components. A sensitivity analysis was conducted to assess the impact of progressive permafrost loss on discharge from the sub‐basin, in which all units of the sub‐basin have the potential to contribute to the drainage network, by incrementally reducing the ratio of wetland to plateau in the modelled sub‐basin. Simulated reductions in permafrost extent decreased total annual discharge from the channel fen by 2.5% for every 10% decrease in permafrost area due to increased surface storage capacity, reduced run‐off efficiency, and increased landscape evapotranspiration. Runoff ratios for the fen hydrological response unit dropped from 0.54 to 0.48 after the simulated 50% permafrost area loss with a substantial reduction of 0.47 to 0.31 during the snowmelt season. The reduction in peat plateau area resulted in decreased seasonal variability in discharge due to changes in the flow path routing, with amplified low flows associated with small increases in subsurface discharge, and decreased peak discharge with large reductions in surface run‐off.     

Modelling a cohesive‐frictional debris flow: an experimental,theoretical, and field‐based study

The rheology of debris flows is difficult to characterize owing to the varied composition and to the uneven distribution of the components that may range from clay to large boulders, in addition to water. Few studies have addressed debris flow rheology from observational, experimental, and theoretical viewpoints in conjunction. We present a coupled rheological‐numerical model to characterize the debris flows in which cohesive and frictional materials are both present. As a first step, we consider small‐scale artificial debris flows in a flume with variable percentages of clay versus sand, and measure separately the rheological properties of sand–clay mixtures. A comparison with the predictions of a modified version of the numerical model BING shows a reasonable agreement between measurements and simulations. As application to a field case, we analyse a recent debris flow that occurred in Fjærland (Western Norway) for which much information is now available. The event was caused by a glacial lake outburst flood (GLOF) originating from the failure of a moraine ridge. In a previous contribution (Breien et al., Landslides, 2008 , 5: 271–280) we focused on the hydrological and geomorphological aspects. In particular we documented the marked erosion and reported the change in sediment transport during the event. In contrast to the laboratory debris flows, the presence of large boulders and the higher normal pressure inside the natural debris flow requires the introduction of a novel rheological model that distinguishes between mud‐to–clast supported material. We present simulations with a modified BING model with the new cohesive‐frictional rheology. To account for the severe erosion operated by the debris flow on the colluvial deposits of Fjærland, we also suggest a simple model for erosion and bulking along the slope path. Numerical simulations suggest that a self‐sustaining mechanism could partly explain the extreme growth of debris flows running on a soft terrain.     

Modelling surface‐water depression storage in a Prairie Pothole Region

In this study, the Precipitation‐Runoff Modelling System (PRMS) was used to simulate changes in surface‐water depression storage in the 1,126‐km² Upper Pipestem Creek basin located within the Prairie Pothole Region of North Dakota, USA. The Prairie Pothole Region is characterized by millions of small water bodies (or surface‐water depressions) that provide numerous ecosystem services and are considered an important contribution to the hydrologic cycle. The Upper Pipestem PRMS model was extracted from the U.S. Geological Survey's (USGS) National Hydrologic Model (NHM), developed to support consistent hydrologic modelling across the conterminous United States. The Geospatial Fabric database, created for the USGS NHM, contains hydrologic model parameter values derived from datasets that characterize the physical features of the entire conterminous United States for 109,951 hydrologic response units. Each hydrologic response unit in the Geospatial Fabric was parameterized using aggregated surface‐water depression area derived from the National Hydrography Dataset Plus, an integrated suite of application‐ready geospatial datasets. This paper presents a calibration strategy for the Upper Pipestem PRMS model that uses normalized lake elevation measurements to calibrate the parameters influencing simulated fractional surface‐water depression storage. Results indicate that inclusion of measurements that give an indication of the change in surface‐water depression storage in the calibration procedure resulted in accurate changes in surface‐water depression storage in the water balance. Regionalized parameterization of the USGS NHM will require a proxy for change in surface‐storage to accurately parameterize surface‐water depression storage within the USGS NHM.     

Controlled‐source electromagnetic monitoring of reservoir oil saturation using a novel borehole‐to‐surface configuration

To advance and optimize secondary and tertiary oil recovery techniques, it is essential to know the areal propagation and distribution of the injected fluids in the subsurface. We investigate the applicability of controlled‐source electromagnetic methods to monitor fluid movements in a German oilfield (Bockstedt, onshore Northwest Germany) as injected brines (highly saline formation water) have much lower electrical resistivity than the oil within the reservoir. The main focus of this study is on controlled‐source electromagnetic simulations to test the sensitivity of various source–receiver configurations. The background model for the simulations is based on two‐dimensional inversion of magnetotelluric data gathered across the oil field and calibrated with resistivity logs. Three‐dimensional modelling results suggest that controlled‐source electromagnetic methods are sensitive to resistivity changes at reservoir depths, but the effect is difficult to resolve with surface measurements only. Resolution increases significantly if sensors or transmitters can be placed in observation wells closer to the reservoir. In particular, observation of the vertical electric field component in shallow boreholes and/or use of source configurations consisting of combinations of vertical and horizontal dipoles are promising. Preliminary results from a borehole‐to‐surface controlled‐source electromagnetic field survey carried out in spring 2014 are in good agreement with the modelling studies.     

A method for coupled arch dam‐foundation‐reservoir seismic behaviour analysis

This paper presents a method for coupled arch dam–foundation–reservoir seismic behaviour analysis. The dam is discretized by finite elements (FE) and the foundation and reservoir are discretized by boundary elements (BE). The opening of contraction joints and the spatial variability of the seismic action is taken into account. The study of Pacoima dam by this method is also presented. The computed results show that no cracks were to be expected due to the vibrations induced during the Feb. 9, 1971 San Fernando earthquake. Copyright © 2000 John Wiley & Sons, Ltd.     

Morphology and sedimentology of a complex debris flow in a clay‐shale basin

Coupling morphological, sedimentological, and rheological studies to numerical simulations is of primary interest in defining debris‐flow hazard on alluvial fans. In particular, numerical runout models must be carefully calibrated by morphological observations. This is particularly true in clay‐shale basins where hillslopes can provide a large quantity of poorly sorted solid materials to the torrent, and thus change both the mechanics of the debris flow and its runout distance. In this context, a study has been completed on the Faucon stream (southeastern French Alps), with the objectives of (1) defining morphological and sedimentological characteristics of torrential watersheds located in clay‐shales, and (2) evaluating through a case study the scouring potential of debris flows affecting a clay‐shale basin. Morphological surveys, grain‐size distributions and petrographic analyses of the debris‐flow deposits demonstrate the granular character of the flow during the first hectometre, and its muddy character from there to its terminus on the debris fan. These observations and laboratory tests suggest that the contributing areas along the channel have supplied the bulk of the flow material. Copyright © 2005 John Wiley & Sons, Ltd.     

2.5D direct‐current resistivity forward modelling and inversion by finite‐element–infinite‐element coupled method

To reduce the numerical errors arising from the improper enforcement of the artificial boundary conditions on the distant surface that encloses the underground part of the subsurface, we present a finite‐element–infinite‐element coupled method to significantly reduce the computation time and memory cost in the 2.5D direct‐current resistivity inversion. We first present the boundary value problem of the secondary potential. Then, a new type of infinite element is analysed and applied to replace the conventionally used mixed boundary condition on the distant boundary. In the internal domain, a standard finite‐element method is used to derive the final system of linear equations. With a novel shape function for infinite elements at the subsurface boundary, the final system matrix is sparse, symmetric, and independent of source electrodes. Through lower upper decomposition, the multi‐pole potentials can be swiftly obtained by simple back‐substitutions. We embed the newly developed forward solution to the inversion procedure. To compute the sensitivity matrix, we adopt the efficient adjoint equation approach to further reduce the computation cost. Finally, several synthetic examples are tested to show the efficiency of inversion.     

Flow behaviour and runout modelling of a complex debris flow in a clay‐shale basin

Identification of debris‐flow hazard areas necessitates the knowledge of the flow thickness and the runout distance. Both have been investigated using a numerical runout model. On the Faucon stream (South French Alps), representative of clay‐shale basins, results of various rheological tests and numerical experiments are presented and discussed. The calibration of the model was undertaken using the results of both geomorphological surveys and sedimentological analyses. Rheological tests using either a parallel‐plate rheometer, a coaxial rheometer, slump tests, and an inclined plane were carried out on several samples. Results have shown that the flow behaviour could be described by an Herschel‐Bulkley constitutive equation. The rheological responses of several natural suspensions collected from surficial deposits (sandstones, moraines, weathered black marls) were also investigated. In order to model the runout of the flow, the model BING was used. The model describes well the influence of each type of sediment on the behaviour (runout distance, deposit thickness) of the flow, although the velocities were significantly overestimated. Different risk scenarios are tested and discussed. Copyright © 2005 John Wiley & Sons, Ltd.     

Effect of non‐uniformity of flow on velocity and turbulence intensities over a cobble‐bed

Non‐uniform flows encompassing both accelerating and decelerating flows over a cobble‐bed flume have been experimentally investigated in a flume at a scale of intermediate relative submergence. Measurements of mean longitudinal flow velocity u, and determinations of turbulence intensities u′, v′, w′, and Reynolds shear stress ?u_fw_f have been made. The longitudinal velocity distribution was divided into the inner zone close to the bed and the outer zone far from the bed. In the inner zone of the boundary layer (near the bed) the velocity profile closely followed the ‘Log Law’; however, in the outer zone the velocity distribution deviated from the Log Law consistently for both accelerating and decelerating flows and the changes in bed slopes ranging from ?2% to + 2% had no considerable effect on the outer zone. For a constant bed slope (S = ±0·015), the larger the flow rate, the smaller the turbulence intensities. However, no detectable pattern has been observed for u′, v′ and w′ distributions near the bed. Likewise, for a constant flow rate (Q = 0·040 m³/s), with variation in bed slope the longitudinal turbulent intensity profile in the longitudinal direction remained concave for both accelerating and decelerating flows; whereas vertical turbulent intensity (w′) profile presented no specific form. The results reveal that the positions of maximum values of turbulence intensities and the Reynolds shear stress depend not only on the flow structure (accelerating or decelerating) but also on the intermediate relative submergence scale. Copyright © 2009 John Wiley & Sons, Ltd.     

Estimation of pore‐pressure change in a compacting reservoir from time‐lapse seismic data

An approach is developed to estimate pore‐pressure changes in a compacting chalk reservoir directly from time‐lapse seismic attributes. It is applied to data from the south‐east flank of the Valhall field. The time‐lapse seismic signal of the reservoir in this area is complex, despite the fact that saturation changes do not have an influence. This complexity reflects a combination of pressure depletion, compaction and stress re‐distribution throughout the reservoir and into the surrounding rocks. A simple relation is found to link the time‐lapse amplitude and time‐shift attributes to variations in the key controlling parameter of initial porosity. This relation is sufficient for an accurate estimation of pore‐pressure change in the inter‐well space. Although the time‐lapse seismic estimates mostly agree with reservoir simulation, unexplained mismatches are apparent at a small number of locations with lower porosities (less than 38%). The areas of difference between the observations and predictions suggest possibilities for simulation model updating or a better understanding of the physics of the reservoir.     

Frequency‐domain waveform modelling and inversion for coupled media using a symmetric impedance matrix

To simulate the seismic signals that are obtained in a marine environment, a coupled system of both acoustic and elastic wave equations is solved. The acoustic wave equation for the fluid region simulates the pressure field while minimizing the number of degrees of freedom of the impedance matrix, and the elastic wave equation for the solid region simulates several elastic events, such as shear waves and surface waves. Moreover, by combining this coupled approach with the waveform inversion technique, the elastic properties of the earth can be inverted using the pressure data obtained from the acoustic region. However, in contrast to the pure acoustic and elastic cases, the complex impedance matrix for the coupled media does not have a symmetric form because of the boundary (continuity) condition at the interface between the acoustic and elastic elements. In this study, we propose a manipulation scheme that makes the complex impedance matrix for acoustic–elastic coupled media to take a symmetric form. Using the proposed symmetric matrix, forward and backward wavefields are identical to those generated by the conventional approach; thus, we do not lose any accuracy in the waveform inversion results. However, to solve the modified symmetric matrix, LDLT factorization is used instead of LU factorization for a matrix of the same size; this method can mitigate issues related to severe memory insufficiency and long computation times, particularly for large‐scale problems.     

Post‐processing of ensemble forecasts in low‐flow period

For water supply, navigational, ecological protection or water quality control purposes, there is a great need in knowing the likelihood of the river level falling below a certain threshold. Ensemble streamflow prediction (ESP) based on simulations of deterministic hydrologic models is widely used to assess this likelihood. Raw ESP results can be biased in both the ensemble means and the spreads. In this study, we applied a modified general linear model post‐processor (GLMPP) to correct these biases. The modified GLMPP is built on the basis of regression of simulated and observed streamflow calculated on the basis of canonical events, instead of the daily values as is carried out in the original GLMPP. We conducted the probabilistic analysis of post‐processed ESP results falling below pre‐specified low‐flow levels at seasonal time scale. Raw ESP forecasts from the 1980 to 2006 periods by four different land surface models (LSMs) in eight large river basins in the continental USA are included in the analysis. The four LSMs are Noah, Mosaic, variable infiltration capacity and Sacramento models. The major results from this study are as follows: (1) a modified GLMPP was proposed on the basis of canonical events; (2) post‐processing can improve the accuracy and reduce the uncertainty of hydrologic forecasts; (3) post‐processing can help deal with the effect of human activity; and (4) raw simulation results from different models vary greatly in different basins. However, post‐processing can always remove model biases under different conditions. Copyright © 2014 John Wiley & Sons, Ltd.     

Some early Russian studies of subsurface storm‐flow processes

Subsurface flow is an important component of the catchment hydrological cycle. Various mechanisms of this process and their role in storm‐flow generation attracted the attention of many researchers throughout the twentieth century. The results of most of these studies are well documented. However, similar studies conducted in the past by many Russian hydrologists have never been made available to the English speaking hydrological community. This paper attempts to fill this gap and briefly review some of these investigations, focusing on their main results. It starts with the review of the early experimental studies (after 1930s–1950s), which allow the main characteristics of subsurface storm flow to be established. This is followed by a review of the research conducted in 1960s, which resulted in some conceptualization of the subsurface flow mechanisms. The paper also draws some parallels between this, mostly unknown, subsurface flow research in Russia and the better‐known contemporary studies of this process. Copyright © 2002 John Wiley & Sons, Ltd.     

Empirical prediction of debris‐flow mobility and deposition on fans

A new method to predict the runout of debris flows is presented. A data base of documented sediment‐transporting events in torrent catchments of Austria, Switzerland and northern Italy has been compiled, using common classification techniques. With this data we test an empirical approach between planimetric deposition area and event volume, and compare it with results from other studies. We introduce a new empirical relation to determine the mobility coefficient as a function of geomorphologic catchment parameters. The mobility coefficient is thought to reflect some of the flow properties during the depositional part of the debris‐flow event. The empirical equations are implemented in a geographical information system (GIS) based simulation program and combined with a simple flow routing algorithm, to determine the potential runout area covered by debris‐flow deposits. For a given volume and starting point of the deposits, a Monte‐Carlo technique is used to produce flow paths that simulate the spreading effect of a debris flow. The runout zone is delineated by confining the simulated potential spreading area in the down slope direction with the empirically determined planimetric deposition area. The debris‐flow volume is then distributed over the predicted area according to the calculated outflow probability of each cell. The simulation uses the ARC‐Objects environment of ESRI© and is adapted to run with high resolution (2·5 m × 2·5 m) digital elevation models, generated for example from LiDAR data. The simulation program called TopRunDF is tested with debris‐flow events of 1987 and 2005 in Switzerland. Copyright © 2009 John Wiley & Sons, Ltd.