Green''''s functions for uniformly distributed loads acting on an inclined line in a poroelastic layered site

Based on one type of practical Biot's equation and the dynamic-stiffness matrices of a poroelastic soil layer and half-space, Green's functions were derived for uniformly distributed loads acting on an inclined line in a poroelastic layered site. This analysis overcomes significant problems in wave scattering due to local soil conditions and dynamic soil-structure interaction. The Green's functions can be reduced to the case of an elastic layered site developed by Wolf in 1985. Parametric studies are then carried out through two example problems.     

A practical method for estimating dynamic soil stiffness on surface of multi‐layered soil

It is important to estimate the influence of layered soil in soil–structure interaction analyses. Although a great number of investigations have been carried out on this subject, there are very few practical methods that do not require complex calculations. In this paper, a simple and practical method for estimating the horizontal dynamic stiffness of a rigid foundation on the surface of multi‐layered soil is proposed. In this method, waves propagating in the soil are traced using the conception of the cone model, and the impulse response function can be calculated directly and easily in the time domain with a good degree of accuracy. The characteristics of the impedance, that is the transformed value to the frequency domain of the obtained impulse response, are studied using two‐ to four‐layered soil models. The cause of the fluctuation of impedance is expressed clearly from its relation to reflected waves from the lower layer boundary in the model. Copyright © 2005 John Wiley & Sons, Ltd.     

Pressure sensitivity of olivine slip systems: first-principle calculations of generalised stacking faults

We have used a first-principle approach based on the calculation of generalised stacking faults (GSF) to study the influence of pressure on the mechanical properties of forsterite. Six cases corresponding to [100] glide over (010), (021) and (001), and [001] glide over (100), (010) and (110) have been considered. The relaxed energy barriers associated with plastic shear have been calculated by constraining the Si atoms to move perpendicular to the fault plane and allowing Mg and O atoms to move in every direction. These conditions, which preserve dilations as a relaxation process, introduce Si–O tetrahedral tilting as an additional relaxation mechanism. Relaxed GSF show little plastic anisotropy of [100] glide over different planes and confirms that [001] glide is intrinsically easier than [100] glide. The GSF are affected by the application of a 10 GPa confining pressure with a different response for each slip system that cannot be explained by sole elastic effect. In particular, [100](010) is found to harden significantly under pressure compared to [001](010). Our results give the first theoretical framework to understand the pressure-induced change of dominant slip systems observed by Couvy et al. (in Eur J Mineral 16(6):877–889, 2004) and P. Raterron et al. (in GRL, submitted). It appears necessary to account for the influence of pressure on the mechanical properties of silicates in the context of the deep Earth.     

9: Processes of tectonism, magmatism and mineralization: Lessons from Europe

Metallogenic provinces in Europe range in age from the Archaean to the Neogene. Deposit types include porphyry copper and epithermal Cu–Au, volcanic-hosted massive sulphide (VMS), orogenic gold, Fe-oxide–Cu–Au, anorthosite Fe–Ti-oxide and sediment-hosted base-metal deposits. Most of them formed during short-lived magmatic events in a wide range of tectonic settings; many can be related to specific tectonic processes such as subduction, hinge retreat, accretion of island arcs, continental collision, lithosphere delamination or slab tear. In contrast, most sediment-hosted deposits in Europe evolved in extensional, continental settings over significant periods of time. In Europe, as elsewhere, ore formation is an integral part of the geodynamic evolution of the Earth's crust and mantle. Many tectonic settings create conditions conducive to the generation of water-rich magma, but the generation of ore deposits appears to be restricted to locations and short periods of change in temperature and stress, imposed by transitory plate motions. Crustal influence is evident in the strong structural controls on the location and morphology of many ore deposits in Europe. Crustal-scale fault–fracture systems, many involving strike-slip elements, have provided the fabric for major plumbing systems. Rapid uplift, as in metamorphic core complexes, and hydraulic fracturing can generate or focus magmatic–hydrothermal fluid flow that may be active for time spans significantly less than a million years. Once a hydrologically stable flow is established, ore formation is strongly dependent on the steep temperature and pressure gradients experienced by the fluid, particularly within the upper crust. In Europe, significant fracture porosity deep in the crystalline basement (1%) is not only important for magmatic–hydrothermal systems, but allows brines to circulate down through sedimentary basins and then episodically upward, expelled seismically to produce sediment-hosted base-metal deposits and Kupferschiefer copper deposits. Emerging research, stimulated by GEODE, can improve the predicting power of numerical simulations of ore-forming processes and help discover the presence of orebodies beneath barren overburden.     

Numerical simulation of percussive drilling

This paper presents a novel dynamical model to analyze the long‐term response of a percussive drilling system. This departs from existing approaches that usually consider a single activation and bit/rock interaction cycle for the analysis of the process performance. The proposed model integrates the axial dynamics of an elastic piston and an elastic drill bit, a motion‐dependent pressure law to drive the piston, and a generalized bit/rock interaction law representative of the dynamic indentation taking place at the bit/rock interface. It applies to down‐the‐hole percussive drilling as well as top‐hole, with minor modifications. The model does not account for the angular motion or the hole cleaning, however. The model is first formulated mathematically; then, a finite‐dimensional approximation is proposed for computations. Numerical analyses of the model response, for a low‐size down‐the‐hole percussive system, follow. The period‐1 stationary response for the reference configuration is studied in detail, and parametric analyses assessing the influence on the rate of penetration of the bit/rock interaction parameters, the feed force, and the percussive activation parameters are conducted. These analyses reveal that the multiscale nature of the process is well captured by the model and recover expected trends for the influence of the parameters. They also suggest that a significant increase of the penetration rate can be achieved by increasing the percussive frequency. Copyright © 2015 John Wiley & Sons, Ltd.     

Integration of aquifer geology,groundwater flow and arsenic distribution in deltaic aquifers – A unifying concept

Groundwater arsenic (As) presents a public health risk of great magnitude in densely populated Asian delta regions, most acutely in the Bengal Basin (West Bengal, India and Bangladesh). Research has focused on the sources, mobilisation, and heterogeneity of groundwater As, but a consistent explanation of As distribution from local to basin scale remains elusive. We show for the Bengal Aquifer System that the numerous, discontinuous silt‐clay layers together with surface topography impose a hierarchical pattern of groundwater flow, which constrains As penetration into the aquifer and controls its redistribution towards discharge zones, where it is re‐sequestered to solid phases. This is particularly so for the discrete periods of As release to groundwater in the shallow subsurface associated with sea level high‐stand conditions of Quaternary inter‐glacial periods. We propose a hypothesis concerning groundwater flow ( S ilt‐clay layers I mpose H ierarchical groundwater flow patterns constraining A rsenic progression [SIHA]), which links consensus views on the As source and history of sedimentation in the basin to the variety of spatial and depth distributions of groundwater As reported in the literature. SIHA reconciles apparent inconsistencies between independent, in some cases contrasting, field observations. We infer that lithological and topographic controls on groundwater flow, inherent to SIHA, apply more generally to deltaic aquifers elsewhere. The analysis suggests that groundwater As may persist in the aquifers of Asian deltas over thousands of years, but in certain regions, particularly at deeper levels, As will not exceed low background concentrations unless groundwater flow systems are short‐circuited by excessive pumping.     

Elastic lateral dynamic impedance functions for a rigid cylindrical shell type foundation

Elastic lateral dynamic impedance functions are defined as the ratio of the lateral dynamic force/moment to the corresponding lateral displacement/rotation at the top ending of a foundation at very small strains. Elastic lateral dynamic impedance functions have a defining influence on the natural frequencies of offshore wind turbines supported on cylindrical shell type foundations, such as suction caissons, bucket foundations, and monopiles. This paper considers the coupled horizontal and rocking vibration of a cylindrical shell type foundation embedded in a fully saturated poroelastic seabed in contact with a seawater half‐space. The formulation of the coupled seawater–shell–seabed vibration problem is simplified by treating the shell as a rigid one. The rigid shell vibration problem is approached by the integral equation method using ring‐load Green's functions for a layered seawater‐seabed half‐space. By considering the boundary conditions at the shell–soil interface, the shell vibration problem is reduced to Fredholm integral equations. Through an analysis of the corresponding Cauchy singular equations, the intrinsic singular characteristics of the problem are rendered explicit. With the singularities incorporated into the solution representation, an effective numerical method involving Gauss–Chebyshev method is developed for the governing Fredholm equations. Selected numerical results for the dynamic contact load distributions, displacements of the shell, and lateral dynamic impedance functions are examined for different shell length–radius ratio, poroelastic materials, and frequencies of excitation. Copyright © 2016 John Wiley & Sons, Ltd.     

Tensionless–frictionless interaction of flexible annular foundation with a transversely isotropic multi‐layered half‐space

A transversely isotropic multi‐layered half‐space, with axis of material symmetry perpendicular to the free surface, supports a flexible either annular or solid circle foundation. The contact area of the foundation and the half‐space is considered to be both frictionless and tensionless. The foundation is assumed to be affected by a vertical static axisymmetric load. Detailed analysis of the interaction of these two systems with different thickness of layers is the target of this paper. With the use of ring load Green's functions for both the foundation and the continuum half‐space, an integral equation accompanied with some inequalities is introduced to model the complex BVP. With the incorporation of ring‐shape FEM, we are capable of capturing both regular and singular solution smoothly. The validity of the combination of the analytical and numerical method is proved with comparing the results of this paper with a number of benchmark cases of both linear and nonlinear interaction of circular and annular foundation with half‐space. Some new illustrations are presented to portray the aspect of the anisotropy and layering of the half‐space. Copyright © 2014 John Wiley & Sons, Ltd.     

Analytical solutions for advective–dispersive solute transport in double‐layered finite porous media

Analytical solutions for advection and dispersion of a conservative solute in a one‐dimensional double‐layered finite porous media are presented. The solutions are applicable to five scenarios that have various combinations of fixed concentration, fixed flux and zero concentration gradient conditions at the inlet and outlet boundaries that provide a wide number of options. Arbitrary initial solute concentration distributions throughout the media can be considered via explicit formulations or numerical integration. The analytical solutions presented have been verified against numerical solutions from a finite‐element‐based approach and an existing closed‐form solution for double‐layered media with an excellent correlation being found in both cases. A practical application pertaining to advective transport induced by consolidation of underlying sediment layers on contaminant movement within a capped contaminated sediment system is presented. Comparison of the calculated concentrations and fluxes with alternative approaches clearly illustrates the need to consider advection processes. Consideration of the different features of contaminant transport due to varying pore‐water velocity fields in primary consolidation and secondary consolidation stages is achieved via the use of non‐uniform initial concentration distributions within the proposed analytical solutions. Copyright © 2010 John Wiley & Sons, Ltd.