A rational approach to seepage flow effects on bottom friction beneath random waves

The bottom friction beneath random waves is predicted taking into account the effect of seepage flow. This is achieved by using wave friction factors for rough turbulent, smooth turbulent and laminar flow valid for regular waves together with a modified Shields parameter which includes the effect of seepage flow. Examples using data typical to field conditions are included to illustrate the approach. The analytical results can be used to make assessment of seepage effects on the bottom friction based on available wave statistics. Generally, it is recommended that a stochastic approach should be used rather than using the rms values in an otherwise deterministic approach.     

The permeability of faults within siliciclastic petroleum reservoirs of the North Sea and Norwegian Continental Shelf

Faulting in Middle Jurassic reservoirs occurred at shallow depth during regional extension. Clean sandstones (<15% clay) deformed without significant grain fracturing and permeability reduction. Faults in impure sandstones (15–40% clay) experienced significant syn-deformation compaction and permeability reduction. Enhanced compaction during deeper burial reduced their permeabilities further from an average of 0.05 mD at <2.5 km to 0.001 mD at >4 km. Clay-rich sediments (>40% clay) deformed to produce clay smears with very low permeabilities (<0.001 mD). Faulting in the Rotliegendes occurred at greater depth during both basin extension and inversion. Extensional faulting produced cataclasites with permeability reductions of <10–>10⁶; their permeabilities range from 0.2 to 0.0001 mD and are inversely related to their maximum burial depth. Faults formed or reactivated during inversion experienced permeability increase. These results can be extrapolated to other hydrocarbon reservoirs if differences in stress and temperature history are taken into account.The permeability of most (>80%) faults is not sufficiently low, compared to their wallrock, to retard single-phase fluid flow on a km-scale. Nevertheless, most faults could retard the flow of a non-wetting phase if present at low saturations. It may be necessary to incorporate the two-phase fluid flow properties of fault rocks into reservoir simulators using upscaling or pseudoisation techniques. Fault property data should be calibrated against production data before it can be used confidently.     

Time-averaged sub-layer concentration profiles in oscillatory sheet flows

Although thin on the order of several to tens of millimeters, sheet flows normally comprise a lower pick-up sub-layer and an upper contact-load sub-layer, separated at the bed level (z = 0). The time-averaged concentration profile in the pick-up sub-layer shows a ‘convex upward’ curvature, but ‘concave upward’ shape characterizes the time-averaged concentration profile in the contact-load sub-layer. The time-dependent concentration in the contact-load sub-layer is approximately in-phase with the free stream flow velocity, whereas it is nearly in anti-phase with the free stream flow velocity in the pick-up sub-layer. Two distinct analytical expressions of the time-averaged concentration profiles for the respective sub-layers are proposed. The expressions are validated with detailed observation datasets collected in the Groβer Wellenkanal (GWK) prototype wave flume in Hannover, Germany. The agreement between the predicted and the measured values is excellent. Interparticle collisions in the pick-up sub-layer and convective lifting processes associated with vortex shedding in the contact-load sub-layer are considered responsible for the opposite curvatures and in-phase/anti-phase concentration variations. Both transitional boundary and reference concentration are also elaborated.     

Issues in Eulerian–Lagrangian modeling of sediment transport under saltation regime

The saltation regime is very important for understanding the sediment transport mechanism. However,there is no consensus on a model for the saltation regime. This study answers several questions raised with respect to the Eulerian-Lagrangian modeling of sediment transport. The first question is why the previous saltation models that use different combinations of hydrodynamic forces yielded acceptable results? The second question is which shear lift model(i.e. a shear lift expression and its coefficient) is more appropriate? Another important question is which hydrodynamic forces have greater contributions to the saltation characteristics of a sediment particle? The last question is what are the contributions of the turbulence fluctuations as well as effects of using two-and three-dimensional(2 D and 3 D) models on the simulation results? In order to fairly answer these questions, a systematic study was done by considering different scenarios. The current study is the first attempt to clearly discuss these issues. A comprehensive 3 D saltation model for non-cohesive sediment was developed that includes all the hydrodynamic forces acting on the particle. The random nature of sediment transport was included using turbulent flow and bed-particle collision models. The eddy interaction model was applied to generate a3 D turbulent flow field. Bed-particle collisions were considered using the concept of a contact zone and a corresponding contact point. The validation of the model was done using the available experimental data for a wide range of sediment size(0.03 to 4.8 cm). For the first question, the results indicated that some of the hydrodynamic effects show opposing trends and some have negligible effects. With these opposing effects it is possible to adjust the coefficients of different models to achieve acceptable agreement with the same experimental data while omitting some aspects of the physics of the process. A suitable model for the shear lift force was developed by linking the lift coefficient to the drag coefficient and the contributions of the hydrodynamic forces and turbulence fluctuations as well as the consequences of using of 2 D and 3 D models were studied. The results indicate that the shear lift force and turbulent flow fluctuations are important factors for the saltation of both sand and gravel, and they cannot be ignored.     

Comparing threshold definition techniques for rainfall‐induced landslides: A national assessment using radar rainfall

Translational landslides and debris flows are often initiated during intense or prolonged rainfall. Empirical thresholds aim to classify the rain conditions that are commonly associated with landslide occurrence and therefore improve understating of these hazards and predictive ability. Objective techniques that are used to determine these thresholds are likely to be affected by the length of the rain record used, yet this is not routinely considered. Moreover, remotely sensed spatially continuous rainfall observations are under‐exploited. This study compares and evaluates the effect of rain record length on two objective threshold selection techniques in a national assessment of Scotland using weather radar data. Thresholds selected by ‘threat score’ are sensitive to rain record length whereas, in a first application to landslides, ‘optimal point’ (OP) thresholds prove relatively consistent. OP thresholds increase landslide detection and may therefore be applicable in early‐warning systems. Thresholds combining 1‐ and 12‐day antecedence variables best distinguish landslide initiation conditions and indicate that Scottish landslides may be initiated by lower rain accumulation and intensities than previously thought. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

A model of secondary upwelling over the shelf break

Motivated by our puzzling high-resolution radar observations of surface vortices in the nearly rectangular Gulf of Eilat/Aqaba, northern Red Sea, we propose and explore the driven cavity approach to this geophysical phenomenon. While the lid-driven cavity has long been considered a benchmark problem in computational fluid dynamics, its oceanographic context has not been considered. Despite the additional effects of rotation and stratification, our modeling demonstrates that when the fluid within a cavity geometrically similar to the Gulf of Eilat is driven by the external current, an interior vortex can develop as in our observations. Furthermore, the Eilat vortices appear only under relatively calm conditions, adding evidence to the intriguing possibility of their simple shear-driven origin.     

Statistical analysis of the effects on overland flow of spatial variability in soil hydraulic conductivity / Analyse statistique des effets de la variabilité spatiale de la conductivité hydraulique du sol sur l'écoulement de surface

Abstract

The effects of spatial variation of the saturated hydraulic conductivity (K _s ) of the soil on the variation of overland flow were tested by analysing 2000 synthetic rainfall—runoff events, all generated from real, observed rainfall events but with runoff modelled by a two-dimensional distributed model using different spatially variable K _s fields in a small (12 ha) agricultural catchment. The purpose is to determine the influence of spatial variation in K _s on runoff generation. The statistical measures used to describe the variation in the generated K _s were its coefficient of variation and correlation length. Both of these had two levels of typical values obtained from field measurements in other studies. The storms were analysed at a general event level, first using simple graphical and statistical methods and then using analysis of variance (ANOVA). The observed scale of the spatial variation of K _s does cause statistically significant variation in overland flow. The graphical analysis showed that the first flow peak in a multi-event storm had the largest variation and that differences were greater in the rising part of the hydrograph than in its recession. The greatest variation in overland flow was produced by the combination of the greater coefficient of variation and the longer correlation lengths. The smallest variation in overland flow was produced by the combination of the smaller coefficient of variation and the shorter correlation lengths. ANOVA showed that the coefficient of variation and correlation length alone did not explain all the variation of the total flow. ANOVA was not very useful due to the many restrictive assumptions that were not satisfied by the nature of the data and therefore analysis methods with less restrictive assumptions need to be tested.     

Investigating hydrological regimes and processes in a set of catchments with temporary waters in Mediterranean Europe

Abstract

Seven catchments of diverse size in Mediterranean Europe were investigated in order to understand the main aspects of their hydrological functioning. The methods included the analysis of daily and monthly precipitation, monthly potential evapotranspiration rates, flow duration curves, rainfall—runoff relationships and catchment internal data for the smaller and more instrumented catchments. The results showed that the catchments were less “dry” than initially considered. Only one of them was really semi-arid throughout the year. All the remaining catchments showed wet seasons when precipitation exceeded potential evapotrans-piration, allowing aquifer recharge, “wet” runoff generation mechanisms and relevant baseflow contribution. Nevertheless, local infiltration excess (Hortonian) overland flow was inferred during summer storms in some catchments and urban overland flow in some others. The roles of karstic groundwater, human disturbance and low winter temperatures were identified as having an important impact on the hydrological regime in some of the catchments.     

Response of the hydrological regime of the Yellow River to the changing monsoon intensity and human activity

Abstract

In the past 50 years, influenced by global climate change, the East Asian summer monsoon intensity (SMI) changed significantly, leading to a response by the water cycle of the Yellow River basin. The variation in SMI has three stages: (1) 1951–1963, SMI increased; (2) 1963–1965, SMI declined sharply, a feature that may be regarded as an abrupt change; and (3) 1965–2000, SMI remained at low levels and showed a tendency to decline slowly. The decreased SMI led to a reduction in water vapour transfer from the ocean to the Yellow River basin, and thus precipitation decreased and the natural river runoff of the Yellow River also decreased. Due to the increase in population and therefore in irrigated land area, the ratio of net water diversion to natural river runoff increased continuously. Comparison of the ratio of net water diversion to natural river runoff before and after the abrupt change in SMI indicates some discontinuity in the response of the man-induced lateral branch of the water cycle to the abrupt change in SMI. The frequently occurring flow desiccation in the lower Yellow River can be regarded as a response of the water cycle system to the decreasing summer monsoon intensity and increasing population. When the ratio of net water diversion exceeded the ratio of natural runoff of the low-flow season to the annual total natural runoff, flow desiccation in the lower Yellow River would occur. When the ratio of net water diversion is 0.3 larger than the ratio of the natural runoff of the low-flow season to the annual total natural runoff, an abrupt increase in the number of flow desiccation events is likely to occur.     

Integrating lysimeter drainage and eddy covariance flux measurements in a groundwater recharge model

Abstract

Field-scale water balance is difficult to characterize because controls exerted by soils and vegetation are mostly inferred from local-scale measurements with relatively small support volumes. Eddy covariance flux and lysimeters have been used to infer and evaluate field-scale water balances because they have larger footprint areas than local soil moisture measurements. This study quantifies heterogeneity of soil deep drainage (D) in four 12.5-m² repacked lysimeters, compares evapotranspiration from eddy covariance (ET_EC) and mass balance residuals of lysimeters (ETwb_Lys), and models D to estimate groundwater recharge. Variation in measured D was attributed to redirection of snowmelt infiltration and differences in lysimeter hydraulic properties caused by surface soil treatment. During the growing seasons of 2010, 2011 and 2012, ETwb_Lys (278, 289 and 269 mm, respectively) was in good agreement with ET_EC (298, 301 and 335 mm). Annual recharge estimated from modelled D was 486, 624 and 613 mm for three calendar years 2010, 2011 and 2012, respectively. In summary, lysimeter D and ET_EC can be integrated to estimate and model groundwater recharge.

Editor D. Koutsoyiannis