Three dimensional elasto‐plastic interface for embedded beam elements with interaction surface for the analysis of lateral loading of piles

This paper presents a numerical formulation for a three dimensional elasto‐plastic interface, which can be coupled with an embedded beam element in order to model its non‐linear interaction with the surrounding solid medium. The formulation is herein implemented for lateral loading of piles but is able to represent soil‐pile interaction phenomena in a general manner for different types of loading conditions or ground movements. The interface is formulated in order to capture localized material plasticity in the soil surrounding the pile within the range of small to moderate lateral displacements. The interface is formulated following two different approaches: (i) in terms of beam degrees of freedoms; and (ii) considering the displacement field of the solid domain. Each of these alternatives has its own advantages and shortcomings, which are discussed in this paper. The paper presents a comparison of the results obtained by means of the present formulation and by other well‐established analysis methods and test results published in the literature. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of channel transmission losses in a dryland river reach in north‐eastern Brazil using streamflow series,groundwater level series and multi‐temporal satellite data

Scarcity of hydrological data, especially streamflow discharge and groundwater level series, restricts the understanding of channel transmission losses (TL) in drylands. Furthermore, the lack of information on spatial river dynamics encompasses high uncertainty on TL analysis in large rivers. The objective of this study was to combine the information from streamflow and groundwater level series with multi‐temporal satellite data to derive a hydrological concept of TL for a reach of the Middle Jaguaribe River (MJR) in semi‐arid north‐eastern Brazil. Based on this analysis, we proposed strategies for its modelling and simulation. TL take place in an alluvium, where river and groundwater can be considered to be hydraulically connected. Most losses certainly infiltrated only through streambed and levees and not through the flood plains, as could be shown by satellite image analysis. TL events whose input river flows were smaller than a threshold did not reach the outlet of the MJR. TL events whose input flows were higher than this threshold reached the outlet losing on average 30% of their input. During the dry seasons (DS) and at the beginning of rainy seasons (DS/BRS), no river flow is expected for pre‐events, and events have vertical infiltration into the alluvium. At the middle and the end of the rainy seasons (MRS/ERS), river flow sustained by base flow occurs before/after events, and lateral infiltration into the alluvium plays a major role. Thus, the MJR shifts from being a losing river at DS/BRS to become a losing/gaining (mostly losing) river at MRS/ERS. A model of this system has to include the coupling of river and groundwater flow processes linked by a leakage approach. Copyright © 2012 John Wiley & Sons, Ltd.     

The influence of local hydrogeologic forcings on near‐stream event water recharge and retention (Upper San Pedro River,Arizona)

The rise in stream stage during high flow events (floods) can induce losing stream conditions, even along stream reaches that are gaining during baseflow conditions. The aquifer response to flood events can affect the geochemical composition of both near‐stream groundwater and post‐event streamflow, but the amount and persistence of recharged floodwater may differ as a function of local hydrogeologic forcings. As a result, this study focuses on how vertical flood recharge varies under different hydrogeologic forcings and the significance that recharge processes can have on groundwater and streamflow composition after floods. River and shallow groundwater samples were collected along three reaches of the Upper San Pedro River (Arizona, USA) before, during and after the 2009 and 2010 summer monsoon seasons. Tracer data from these samples indicate that subsurface floodwater propagation and residence times are strongly controlled by the direction and magnitude of the dominant stream–aquifer gradient. A reach that is typically strongly gaining shows minimal floodwater retention shortly after large events, whereas the moderately gaining and losing reaches can retain recharged floodwater from smaller events for longer periods. The moderately gaining reach likely returned flood recharge to the river as flow declined. These results indicate that reach‐scale differences in hydrogeologic forcing can control (i) the amount of local flood recharge during events and (ii) the duration of its subsurface retention and possible return to the stream during low‐flow periods. Our observations also suggest that the presence of floodwater in year‐round baseflow is not due to long‐term storage beneath the streambed along predominantly gaining reaches, so three alternative mechanisms are suggested: (i) repeated flooding that drives lateral redistribution of previously recharged floodwater, (ii) vertical recharge on the floodplain during overbank flow events and (iii) temporal variability in the stream–aquifer gradient due to seasonally varying water demands of riparian vegetation. Copyright © 2011 John Wiley & Sons, Ltd.     

Sedimentation of overbank floods in the confined complex channel–floodplain system of the Lower Yellow River,China

The channel boundary conditions along the Lower Yellow River (LYR) have been altered significantly since the 1950s with the continual reinforcement and construction of both main and secondary dykes and river training works. To evaluate how the confined complex channel–floodplain system of the LYR responds to floods, this study presents a detailed investigation of the relationship between the tempo‐spatial distribution of sedimentation/erosion and overbank floods occurred in the LYR. For large overbank floods, we found that when the sediment transport coefficient (ratio of sediment concentration of flow to flow discharge) is less than 0.034, the bankfull channel is subject to significant erosion, whereas the main and secondary floodplains both accumulate sediment. The amount of sediment deposited on the main and secondary floodplains is closely related to the ratio of peak discharge to bankfull discharge, volume of water flowing over the floodplains, and sediment concentration of overbank flow, whereas the degree of erosion in the bankfull channel is related to the amount of sediment deposited on the main and secondary floodplains, water volume, and sediment load in flood season. The significant increase in erosion in the bankfull channel is due to the construction of the main and secondary dykes and river training works, which are largely in a wide and narrow alternated pattern along the LYR such that the water flowing over wider floodplains returns to the channel downstream after it drops sediment. For small overbank floods, the bankfull channel is subject to erosion when the sediment transport coefficient is less than 0.028, whereas the amount of sediment deposited on the secondary floodplain is associated closely with the sediment concentration of flow. Over the entire length of the LYR, the situation of erosion in the bankfull channel and sediment deposition on the main and secondary floodplains occurred mainly in the upper reach of the LYR, in which a channel wandering in planform has been well developed.     

The life cycle of the intense IOP‐14 storm during CASPII. Part I: Analysis and simulations

Abstract

In this study, a 5‐day life‐cycle of the IOP‐14 storm during CASP II is examined using conventional observations and numerical simulations with a mesoscale version of the Canadian Regional Finite‐Element (RFE) model. Observational analysis reveals that the IOP‐14 storm forms from a lee trough, occurring along a strong baroclinic zone with an intense frontogenetic deformation, that interacts with an upper‐level travelling short‐wave trough across the Canadian Rockies. Then the storm experiences a slow, but nearly steady, growth while traversing the North American continent. It deepens explosively as it moves into the Atlantic Ocean. It appears that i) the enhanced large‐scale baroclinicity due to land‐sea temperature contrasts, ii) the tremendous latent heat release due to the transport of high‐θ_e air from the marine boundary layer, Hi) the decrease of surface drag and iv) the favourable westward tilt of the low with an amplifying trough all contribute to the explosive deepening of the storm.

Two consecutive simulations covering a total of 102 h during the storm development are carried out with a grid size of 50 km. The RFE model reproduces very well the formation of the surface low on the lee side of the Rockies, the track and deepening rates, the explosive development and decay of the storm, and various mesoscale phenomena (e.g., a “bent‐back” warm front, a “T‐bone” thermal pattern, a cold frontal “fracture”, an upper‐level “eye” and warm‐core structures), as verified by conventional observations, satellite imagery, flight‐level and dropsonde data from a research aircraft. It is found from potential vorticity (PV) analysis that the storm reaches its peak intensity as the upper‐level dry PV anomaly, the low‐level moist PV anomaly and surface thermal warmth are vertically superposed. PV inversions reveal that these anomalies contribute about 60%, 30% and 10%, respectively, to the 900‐hPa negative height perturbation. It is shown that the warm‐core structure near the cyclone centre is produced by advection of warmer air ahead of the cold front, rather than by adiabatic warming associated with subsidence.     

On the testing of fully integrated surface–subsurface hydrological models

Studies employing integrated surface–subsurface hydrological models (ISSHMs) have utilized a variety of test cases to demonstrate model accuracy and consistency between codes. Here, we review the current state of ISSHM testing and evaluate the most popular ISSHM test cases by comparing the hydrodynamic processes simulated in each case to the processes found in well‐characterized, real‐world catchments and by comparing their general attributes to those of successful benchmark problems from other fields of hydrogeology. The review reveals that (1) ISSHM testing and intercode comparison have not adopted specific test cases consistently; (2) despite the wide range of ISSHM metrics available for model testing, only two model performance diagnostics are typically adopted: the catchment outflow hydrograph and the catchment water balance; (3) in intercode comparisons, model performance is usually judged by evaluating only one performance diagnostic: the catchment outflow hydrograph; and (4) ISSHM test cases evaluate a small number of hydrodynamic processes that are largely uniform across the model domain, representing a limited selection of the processes of interest in well‐characterized, real‐world catchments. ISSHM testing would benefit from more intercode comparisons using a consistent set of test cases, aimed at evaluating more catchment processes (e.g. flooding) and using a wider range of simulation diagnostics (e.g. pressure head distributions). To achieve this, a suite of test case variations is required to capture the relevant catchment processes. Finally, there is a need for additional ISSHM test problems that compare model predictions with hydrological observations from intensively monitored field sites and controlled laboratory experiments. Copyright © 2012 John Wiley & Sons, Ltd.     

The importance of wave-tide-surge interaction to coastal material transport and diffusion

This study is to combine a coastal high-resolution (2′×2′) two-way coupled wave-tide-surge numerical model (including 3 main physical mechanisms) with a material transport/diffusion model for understanding the law of material transport/diffusion. Results show that the law of material trans- port/diffusion driven by background current field simulated by the coupled wave-tide-surge model is dif- ferent from that simulated by pure tide-surge, and more different from traditional ones driven by tidal current. The coupled background current should be taken into account for the simulation.     

Response evaluation of axially loaded fixed‐head pile groups in clayey soils

The aim of this paper is to investigate the interaction between the piles in a group with a rigid head and correlate the response of a group of piles to that of a single pile. For this purpose, a computationally intensive study using 3‐D nonlinear numerical analysis was carried out for different pile group arrangements in clayey soils. The responses of the groups of piles were compared with that of a single pile and the variation of the settlement amplification factor R_a was then quantified. The influence of the number of piles, the spacing, and the settlement level on the group response is discussed. A previously proposed relationship for predicting the response of a pile group, based on its configuration and the response of a single pile, has been modified to extend its applicability for any pile spacing. The modified relationship provides a reasonable prediction for various group configurations in clayey soils. Copyright © 2009 John Wiley & Sons, Ltd.     

An extended Greenwood‐Williamson modelbased normal interaction law for discrete element modelling of spherical particles with surface roughness

The current work aims to develop an improved random normal interaction law based on an extended Greenwood‐Williamson (GW) model for spherical particles with surface roughness in the discrete element modelling of particle systems. The extended GW model overcomes some theoretical defects of the classic GW model when incorporated into the discrete element framework. Based on 2 nondimensional forms in which only 2 surface roughness parameters are involved, an empirical formula of the improved interaction law is derived by the curve‐fitting technique. The resulting interaction law is incorporated into discrete element modelling to investigate the mechanical response of particle systems with different surface roughness. Numerical simulations are performed to model 1‐dimensional and 3‐dimensional compression tests to explore the macro and micro characteristics of granular particles with surface roughness. The results show that surface roughness makes the initial packing of a particle assembly looser and has a greater influence on looser packed samples as expected, but an assembly with moderate roughness may exhibit a higher strength. The limitations of the current development are also highlighted.     

The importance of incorporating diurnally fluctuating stream discharge in stream temperature energy balance models

Although stream temperature energy balance models are useful to predict temperature through time and space, a major unresolved question is whether fluctuations in stream discharge reduce model accuracy when not exactly represented. However, high‐frequency (e.g., subdaily) discharge observations are often unavailable for such simulations, and therefore, diurnal streamflow fluctuations are not typically represented in energy balance models. These fluctuations are common due to evapotranspiration, snow pack or glacial melt, tidal influences within estuaries, and regulated river flows. In this work, we show when to account for diurnally fluctuating streamflow. To investigate how diurnal streamflow fluctuations affect predicted stream temperatures, we used a deterministic stream temperature model to simulate stream temperature along a reach in the Quilcayhuanca Valley, Peru, where discharge varies diurnally due to glacial melt. Diurnally fluctuating streamflow was varied alongside groundwater contributions via a series of computational experiments to assess how uncertainty in reach hydrology may impact simulated stream temperature. Results indicated that stream temperatures were more sensitive to the rate of groundwater inflow to the reach compared with the timing and amplitude of diurnal fluctuations in streamflow. Although incorporating observed diurnal fluctuations in discharge resulted in a small improvement in model RMSE, we also assessed other diurnal discharge signals and found that high amplitude signals were more influential on modelled stream temperatures when the discharge peaked at specific times. Results also showed that regardless of the diurnal discharge signal, the estimated groundwater flux to the reach only varied from 1.7% to 11.7% of the upstream discharge. However, diurnal discharge fluctuations likely have a stronger influence over longer reaches and in streams where the daily range in discharge is larger, indicating that diurnal fluctuations in stream discharge should be considered in certain settings.