Spatiotemporal variability in hydrochemistry of shallow groundwater in a small pre‐alpine catchment: The importance of landscape elements

Topography and landscape characteristics affect the storage and release of water and, thus, groundwater dynamics and chemistry. Quantification of catchment scale variability in groundwater chemistry and groundwater dynamics may therefore help to delineate different groundwater types and improve our understanding of which parts of the catchment contribute to streamflow. We sampled shallow groundwater from 34 to 47 wells and streamflow at seven locations in a 20‐ha steep mountainous catchment in the Swiss pre‐Alps, during nine baseflow snapshot campaigns. The spatial variability in electrical conductivity, stable water isotopic composition, and major and trace ion concentrations was large and for almost all parameters larger than the temporal variability. Concentrations of copper, zinc, and lead were highest at sites that were relatively dry, whereas concentrations of manganese and iron were highest at sites that had persistent shallow groundwater levels. The major cation and anion concentrations were only weakly correlated to individual topographic or hydrodynamic characteristics. However, we could distinguish four shallow groundwater types based on differences from the catchment average concentrations: riparian zone‐like groundwater, hillslopes and areas with small upslope contributing areas, deeper groundwater, and sites characterized by high magnesium and sulfate concentrations that likely reflect different bedrock material. Baseflow was not an equal mixture of the different groundwater types. For the majority of the campaigns, baseflow chemistry most strongly resembled riparian‐like groundwater for all but one subcatchment. However, the similarity to the hillslope‐type groundwater was larger shortly after snowmelt, reflecting differences in hydrologic connectivity. We expect that similar groundwater types can be found in other catchments with steep hillslopes and wet areas with shallow groundwater levels and recommend sampling of groundwater from all landscape elements to understand groundwater chemistry and groundwater contributions to streamflow.     

On the pore‐scale mechanisms leading to brittle and ductile deformation behavior of crystalline rocks

Deformation mechanisms at the pore scale are responsible for producing large strains in porous rocks. They include cataclastic flow, dislocation creep, dynamic recrystallization, diffusive mass transfer, and grain boundary sliding, among others. In this paper, we focus on two dominant pore‐scale mechanisms resulting from purely mechanical, isothermal loading: crystal plasticity and crofracturing. We examine the contributions of each mechanism to the overall behavior at a scale larger than the grains but smaller than the specimen, which is commonly referred to as the mesoscale. Crystal plasticity is assumed to occur as dislocations along the many crystallographic slip planes, whereas microfracturing entails slip and frictional sliding on microcracks. It is observed that under combined shear and tensile loading, microfracturing generates a softer response compared with crystal plasticity alone, which is attributed to slip weakening where the shear stress drops to a residual level determined by the frictional strength. For compressive loading, however, microfracturing produces a stiffer response than crystal plasticity because of the presence of frictional resistance on the slip surface. Behaviors under tensile, compressive, and shear loading invariably show that porosity plays a critical role in the initiation of the deformation mechanisms. Both crystal plasticity and microfracturing are observed to initiate at the peripheries of the pores, consistent with results of experimental studies. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical and analytical observations on long and infinite slopes

Interest in the mechanics of landslides has led to renewed evaluation of the infinite slope equations, and the need for a more general framework for estimating the factor of safety of long and infinite slopes involving non‐homogeneous soil profiles. The paper describes finite element methods that demonstrate the potential for predicting failure in long slope profiles where the critical mechanism is not necessarily at the base of the soil layer. The influence of slope angle is also examined in long slopes, leading to some counter‐intuitive conclusions about the impact of slope steepness on the factor of safety. Copyright © 2010 John Wiley & Sons, Ltd.     

Impact of climate indicators on continental‐scale potential groundwater recharge in Africa

In the last decades, human activity has been contributing to climate change that is closely associated with an increase in temperatures, increase in evaporation, intensification of extreme dry and wet rainfall events, and widespread melting of snow and ice. Understanding the intricate linkage between climate warming and the hydrological cycle is crucial for sustainable management of groundwater resources, especially in a vulnerable continent like Africa. This study investigates the relationship between climate‐change drivers and potential groundwater recharge (PGR) patterns across Africa for a long‐term record (1960–2010). Water‐balance components were simulated by using the PCR‐GLOBWB model and were reproduced in both gridded maps and latitudinal trends that vary in space with minima on the Tropics and maxima around the Equator. Statistical correlations between temperature, storm occurrences, drought, and PGR were examined in six climatic regions of Africa. Surprisingly, different effects of climate‐change controls on PGR were detected as a function of latitude in the last three decades (1980–2010). Temporal trends observed in the Northern Hemisphere of Africa reveal that the increase in temperature is significantly correlated to the decline of PGR, especially in the Northern Equatorial Africa. The climate indicators considered in this study were unable to explain the alarming negative trend of PGR observed in the Sahelian region, even though the Standardized Precipitation‐Evapotranspiration Index (SPEI) values report a 15% drought stress. On the other hand, increases in temperature have not been detected in the Southern Hemisphere of Africa, where increasing frequency of storm occurrences determine a rise of PGR, particularly in southern Africa. Time analysis highlights a strong seasonality effect, while PGR is in‐phase with rainfall patterns in the summer (Northern Hemisphere) and winter (Southern Hemisphere) and out‐of‐phase during the fall season. This study helps to elucidate the mechanism of the processes influencing groundwater resources in six climatic zones of Africa, even though modelling results need to be validated more extensively with direct measurements in future studies. Copyright © 2016 John Wiley & Sons, Ltd.     

Frequency analysis of nonstationary annual maximum flood series using the time‐varying two‐component mixture distributions

The most popular practice for analysing nonstationarity of flood series is to use a fixed single‐type probability distribution incorporated with the time‐varying moments. However, the type of probability distribution could be both complex because of distinct flood populations and time‐varying under changing environments. To allow the investigation of this complex nature, the time‐varying two‐component mixture distributions (TTMD) method is proposed in this study by considering the time variations of not only the moments of its component distributions but also the weighting coefficients. Having identified the existence of mixed flood populations based on circular statistics, the proposed TTMD was applied to model the annual maximum flood series of two stations in the Weihe River basin, with the model parameters calibrated by the meta‐heuristic maximum likelihood method. The performance of TTMD was evaluated by different diagnostic plots and indexes and compared with stationary single‐type distributions, stationary mixture distributions and time‐varying single‐type distributions. The results highlighted the advantages of TTMD with physically‐based covariates for both stations. Besides, the optimal TTMD models were considered to be capable of settling the issue of nonstationarity and capturing the mixed flood populations satisfactorily. Copyright © 2016 John Wiley & Sons, Ltd.     

A discontinuous finite element suspended sediment transport model for water quality assessments in river networks

Suspended sediment plays an important role in the distribution and transport of many pollutants (such as radionuclides) in rivers. Pollutants may adsorb on fine suspended particles (e.g. clay) and spread according to the suspended sediment movement. Hence, the simulation of the suspended sediment mechanism is indispensable for realistic transport modelling. This paper presents and tests a simple mathematical model for predicting the suspended sediment transport in river networks. The model is based on the van Rijn suspended load formula and the advection–diffusion equation with a source or sink term that represents the erosion or deposition fluxes. The transport equation is solved numerically with the discontinuous finite element method. The model evaluation was performed in two steps, first by comparing model simulations with the measured suspended sediment concentrations in the Grote Nete–Molse Nete River in Belgium, and second by a model intercomparison with the sediment transport model NST MIKE 11. The simulations reflect the measurements with a Nash‐Sutcliffe model efficiency of 0.6, while the efficiency between the proposed model and the NST MIKE 11 simulations is 0.96. Both evaluations indicate that the proposed sediment transport model, that is sufficiently simple to be practical, is providing realistic results.     

Analysis of shaft resistance of jacked piles in sands

In recent years, pile jacking has become a viable alternative installation method for displacement piles. Pile jacking produces minimal noise, vibration and air pollution during installation. In addition, it is possible, at the end of jacking, to have a good estimate of the ultimate static capacity of the pile. In this paper, the shaft resistance of piles jacked into sand is studied using one‐dimensional finite element analysis. The finite element simulations, using a two‐surface plasticity model, demonstrate the effects of relative density and confinement on the unit shaft resistance of piles jacked in sand. The impact of the number of jacking strokes on the unit shaft capacity is also assessed. Based on the numerical results, we developed equations for shaft resistance quantifying the effects of relative density, initial confinement and number of jacking strokes. Predictions using these equations are compared with data obtained from centrifuge tests and field tests. Copyright © 2010 John Wiley & Sons, Ltd.     

Steady‐state analytical solutions of flow and deformation coupling due to a point sink within a finite fluid‐saturated poroelastic layer

An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.     

Extended finite element framework for fault rupture dynamics including bulk plasticity

We present an explicit extended finite element framework for fault rupture dynamics accommodating bulk plasticity near the fault. The technique is more robust than the standard split‐node method because it can accommodate a fault propagating freely through the interior of finite elements. To fully exploit the explicit algorithmic framework, we perform mass lumping on the enriched finite elements that preserve the kinetic energy of the rigid body and enrichment modes. We show that with this technique, the extended FE solution reproduces the standard split‐node solution, but with the added advantage that it can also accommodate randomly propagating faults. We use different elastoplastic constitutive models appropriate for geomaterials, including the Mohr–Coulomb, Drucker–Prager, modified Cam‐Clay, and a conical plasticity model with a compression cap, to capture off‐fault bulk plasticity. More specifically, the cap model adds robustness to the framework because it can accommodate various modes of deformation, including compaction, dilatation, and shearing. Copyright © 2013 John Wiley & Sons, Ltd.     

Review and enhancement of 3D concrete models for large‐scale numerical simulations of concrete structures

The present paper focuses on selected plasticity and damage‐plasticity models for describing the 3D material behavior of concrete. In particular, a plasticity model and a damage‐plasticity model are reviewed and evaluated. Based on the results of the evaluation, enhancements are proposed, aiming at improving the correspondence between predicted and observed material behavior and aiming at implementing a robust and efficient stress update algorithm in a finite element program for performing large‐scale 3D numerical simulations of concrete structures. The capabilities of the concrete models are demonstrated by 3D numerical simulations of benchmark tests with combined bending and torsional loading and combined compression and shear loading and by a large‐scale 3D finite element analysis of a model test of a concrete arch dam. Copyright © 2011 John Wiley & Sons, Ltd.