A smoothed particle hydrodynamics modelling of soil–water mixing and resulting changes in average strength

Soil–water interaction is a pivotal process in many underwater geohazards such as underwater landslides where soil sediments gradually evolve into turbidity currents after interactions with ambient water. Due to the large deformations, multiphase interactions and phase changes this involves, investigations from numerical modelling of the transition process have been limited so far. This study explores a simple numerical replication of such soil–water mixing with respect to changes in average strength using smoothed particle hydrodynamics (SPH). A uniform viscoplastic model is used for both the solid-like and fluid-like SPH particles. The proposed numerical solution scheme is verified by single-phase dam break tests and multiphase simple shear tests. SPH combinations of solid-like and fluid-like particles can replicate the clay–water mixture as long as the liquidity index of the solid-like particles is larger than unity. The proposed numerical scheme is shown to capture key features of an underwater landslide such as hydroplaning, water entrainment and wave generation and thus shows promise as a tool to simulate the whole process of subaquatic geohazards involving solid–fluid transition during mass transport.     

Elastic models for nonlinear response of rigid passive piles

Recent study indicates that the response of rigid passive piles is dominated by elastic pile–soil interaction and may be estimated using theory for lateral piles. The difference lies in that passive piles normally are associated with a large scatter of the ratio of maximum bending moment over maximum shear force and induce a limiting pressure that is ~1/3 that on laterally loaded piles. This disparity prompts this study. This paper proposes pressure‐based pile–soil models and develops their associated solutions to capture response of rigid piles subjected to soil movement. The impact of soil movement was encapsulated into a power‐law distributed loading over a sliding depth, and load transfer model was adopted to mimic the pile–soil interaction. The solutions are presented in explicit expressions and can be readily obtained. They are capable of capturing responses of model piles in a sliding soil owing to the impact of sliding depth and relative strength between sliding and stable layer on limiting force prior to ultimate state. In comparison with available solutions for ultimate state, this study reveals the 1/3 limiting pressure (of the active piles) on passive piles was induced by elastic interaction. The current models employing distributed pressure for moving soil are more pertinent to passive piles (rather than plastic soil flow). An example calculation against instrumented model piles is provided, which demonstrates the accuracy of the current solutions for design slope stabilising piles. Copyright © 2014 John Wiley & Sons, Ltd.     

Water level observations from unmanned aerial vehicles for improving estimates of surface water–groundwater interaction

Integrated hydrological models are usually calibrated against observations of river discharge and piezometric head in groundwater aquifers. Calibration of such models against spatially distributed observations of river water level can potentially improve their reliability and predictive skill. However, traditional river gauging stations are normally spaced too far apart to capture spatial patterns in the water surface, whereas spaceborne observations have limited spatial and temporal resolution. Unmanned aerial vehicles can retrieve river water level measurements, providing (a) high spatial resolution; (b) spatially continuous profiles along or across the water body, and (c) flexible timing of sampling. A semisynthetic study was conducted to analyse the value of the new unmanned aerial vehicle‐borne datatype for improving hydrological models, in particular estimates of groundwater–surface water (GW–SW) interaction. Mølleåen River (Denmark) and its catchment were simulated using an integrated hydrological model (MIKE 11–MIKE SHE). Calibration against distributed surface water levels using the Differential Evolution Adaptive Metropolis algorithm demonstrated a significant improvement in estimating spatial patterns and time series of GW–SW interaction. After water level calibration, the sharpness of the estimates of GW–SW time series improves by ~50% and root mean square error decreases by ~75% compared with those of a model calibrated against discharge only.     

Dynamic response of a pile considering the interaction of pile variable cross section with the surrounding layered soil

Assuming that the pile variable cross section interacts with the surrounding soil in the same way as the pile toe does with the bearing stratus, the interaction of pile variable cross section with the surrounding soil is represented by a Voigt model, which consists of a spring and a damper connected in parallel, and the spring constant and damper coefficient are derived. Thus, a more rigid pile–soil interaction model is proposed. The surrounding soil layers are modeled as axisymmetric continuum in which its vertical displacements are taken into account and the pile is assumed to be a Rayleigh–Love rod with material damping. Allowing for soil properties and pile defects, the pile–soil system is divided into several layers. By means of Laplace transform, the governing equations of soil layers are solved in frequency domain, and a new relationship linking the impedance functions at the variable‐section interface between the adjacent pile segments is derived using a Heaviside step function, which is called amended impedance function transfer method. On this basis, the impedance function at pile top is derived by amended impedance function transfer method proposed in this paper. Then, the velocity response at pile top can be obtained by means of inverse Fourier transform and convolution theorem. The effects of pile–soil system parameters are studied, and some conclusions are proposed. Then, an engineering example is given to confirm the rationality of the solution proposed in this paper. Copyright © 2017 John Wiley & Sons, Ltd.     

An analytical solution for deforming twin‐parallel tunnels in an elastic half plane

The interaction between twin‐parallel tunnels affects the tunnelling‐induced ground deformation, which may endanger the nearby structures. In this paper, an analytical solution is presented for problems in determining displacements and stresses around deforming twin‐parallel tunnels in an elastic half plane, on the basis of complex variable theory. As an example, a uniform radial displacement was assumed as the boundary condition for each of the two tunnels. Special attention was paid to the effects of tunnel depth and spacing between the two tunnels on the surface movement to gain deep insight into the effect of the interaction between twin‐parallel tunnels using the proposed analytical approach. It is revealed that the influence of twin tunnel interaction on surface movements diminishes with both the increase of the tunnel depth and the spacing between the two tunnels. The presented analytical solution manifests that, similar to most of the existing numerical results, the principle of superposition can be applied to determine ground deformation of twin‐parallel tunnels with a certain large depth and spacing; otherwise, the interaction effect between the two tunnels should be taken into account for predicting reliable ground movement. Copyright © 2014 John Wiley & Sons, Ltd.     

Scales and magnitude of hyporheic,river–aquifer and bank storage exchange fluxes

Many studies have investigated the exchange processes that occur between rivers and groundwater systems and have successfully quantified the water fluxes involved. Specifically, these exchange processes include hyporheic exchange, river–aquifer exchange (groundwater discharge and river loss) and bank storage exchange. Remarkably, there are relatively few examples of field studies where more than one exchange process is quantified, and as a consequence, the relationships between them are not well understood. To compare the relative magnitudes of these common exchange processes, we have collected data from 54 studies that have quantified one or more of these exchange flux types. Each flux value is plotted against river discharge at the time of measurement to allow the different exchange flux types to be compared. We show that there are positive relationships between the magnitude of each exchange flux type and increasing river discharge across the different studies. For every one order of magnitude increase in river discharge, the hyporheic, river–aquifer and bank storage exchange fluxes increase by factors of 2.7, 2.9 and 2.5, respectively. On average, hyporheic exchange fluxes are almost an order of magnitude greater than river–aquifer exchange fluxes, which are, in turn, approximately four times greater than bank storage exchange fluxes for the same river discharge. Unless measurement approaches that can distinguish between different types of exchange flux are used, there is potential for hyporheic exchange fluxes to be misinterpreted as river–aquifer exchange fluxes, with possible implications for water resource management decisions. Copyright © 2015 John Wiley & Sons, Ltd.     

Improved spatial delineation of streambed properties and water fluxes using distributed temperature sensing

A new method was developed for analysing and delineating streambed water fluxes, flow conditions and hydraulic properties using coiled fibre‐optic distributed temperature sensing or closely spaced discrete temperature sensors. This method allows for a thorough treatment of the spatial information embedded in temperature data by creating a matrix visualization of all possible sensor pairs. Application of the method to a 5‐day field dataset reveals the complexity of shallow streambed thermal regimes. To understand how velocity estimates are affected by violations of assumptions of one‐dimensional, saturated, homogeneous flow and to aid in the interpretation of field observations, the method was also applied to temperature data generated by numerical models of common field conditions: horizontal layering, presence of lateral flow and variable streambed saturation. The results show that each condition creates a distinct signature visible in the triangular matrices. The matrices are used to perform a comparison of the behaviour of one‐dimensional analytical heat‐tracing models. The results show that the amplitude ratio‐based method of velocity calculation leads to the most reliable estimates. The minimum sensor spacing required to obtain reliable velocity estimates with discrete sensors is also investigated using field data. The developed method will aid future heat‐tracing studies by providing a technique for visualizing and comparing results from fibre‐optic distributed temperature sensing installations and testing the robustness of analytical heat‐tracing models. Copyright © 2016 John Wiley & Sons, Ltd.     

Modelling groundwater/surface water interaction in a managed riparian chalk valley wetland

Understanding hydrological processes in wetlands may be complicated by management practices and complex groundwater/surface water interactions. This is especially true for wetlands underlain by permeable geology, such as chalk. In this study, the physically based, distributed model MIKE SHE is used to simulate hydrological processes at the Centre for Ecology and Hydrology River Lambourn Observatory, Boxford, Berkshire, UK. This comprises a 10‐ha lowland, chalk valley bottom, riparian wetland designated for its conservation value and scientific interest. Channel management and a compound geology exert important, but to date not completely understood, influences upon hydrological conditions. Model calibration and validation were based upon comparisons of observed and simulated groundwater heads and channel stages over an equally split 20‐month period. Model results are generally consistent with field observations and include short‐term responses to events as well as longer‐term seasonal trends. An intrinsic difficulty in representing compressible, anisotropic soils limited otherwise excellent performance in some areas. Hydrological processes in the wetland are dominated by the interaction between groundwater and surface water. Channel stage provides head boundaries for broad water levels across the wetland, whilst areas of groundwater upwelling control discrete head elevations. A relic surface drainage network confines flooding extents and routes seepage to the main channels. In‐channel macrophyte growth and its management have an acute effect on water levels and the proportional contribution of groundwater and surface water. The implications of model results for management of conservation species and their associated habitats are discussed. Copyright © 2015 John Wiley & Sons, Ltd.     

Taxonomic analysis of the Quaternary archaeofauna found at The Lapa do Santo site,Lagoa Santa region,Brazil

The Lapa do Santo archaeological site, located in the Lagoa Santa region, Brazil, represents an important hunter-gatherer occupation dated from the Early and Middle Holocene. Prior studies of archaeofauna dating to this period are few and most of them only provide basic faunal identification with limited information on taphonomic processes. The main goal of this study is to identify the archaeofauna, record the taphonomic processes, and make inferences about its natural (interpreted as the death of the animal in the rockshelter area due to natural causes or due to predation) or anthropic origins. No extinct species (including megamastofauna) were identified. Taxonomic analysis indicated that most faunal remains consist of taxa that may be of mixed origin (natural or anthropic), such as microvertebrates and carnivorous mammals. The Cervidae are the main vertebrate family found in the studied material and these are likely the result of hunting. The Cervidae Ozotoceros support the presence of a fauna that is typical of savannah environments, although the assemblage contains mostly generalist fauna that can thrive in a variety of environments. Very few remains showed clear evidence for taphonomic processes related to human interaction, such as burnt bones.     

Soil‐pile interaction in the pile vertical vibration considering true three‐dimensional wave effect of soil

The dynamic response of an end bearing pile embedded in a linear visco‐elastic soil layer with hysteretic type damping is theoretically investigated when the pile is subjected to a time‐harmonic vertical loading at the pile top. The soil is modeled as a three‐dimensional axisymmetric continuum in which both its radial and vertical displacements are taken into account. The pile is assumed to be vertical, elastic and of uniform circular cross section. By using two potential functions to decompose the displacements of the soil layer and utilizing the separation of variables technique, the dynamic equilibrium equation is uncoupled and solved. At the interface of soil‐pile system, the boundary conditions of displacement continuity and force equilibrium are invoked to derive a closed‐form solution of the vertical dynamic response of the pile in frequency domain. The corresponding inverted solutions in time domain for the velocity response of a pile subjected to a semi‐sine excitation force applied at the pile top are obtained by means of inverse Fourier transform and the convolution theorem. A comparison with two other simplified solutions has been performed to verify the more rigorous solutions presented in this paper. Using the developed solutions, a parametric study has also been conducted to investigate the influence of the major parameters of the soil‐pile system on the vertical vibration characteristics of the pile. Copyright © 2013 John Wiley & Sons, Ltd.