A methodology for regionalizing 3‐D effective porosity at watershed scale in crystalline aquifers

An innovative approach for regionalizing the 3‐D effective porosity field is presented and applied to two large, overexploited, and deeply weathered crystalline aquifers located in southern India. The method derives from earlier work on regionalizing a 2‐D effective porosity field in that part of an aquifer where the water table fluctuates, which is now extended over the entire aquifer using a 3‐D approach. A method based on geological and geophysical surveys has also been developed for mapping the weathering profile layers (saprolite and fractured layers). The method for regionalizing 3‐D effective porosity combines water table fluctuation and groundwater budget techniques at various cell sizes with the use of satellite‐based data (for groundwater abstraction), the structure of the weathering profile, and geostatistical techniques. The approach is presented in detail for the Kudaliar watershed (983 km²) and tested on the 730 km² Anantapur watershed. At watershed scale, the effective porosity of the aquifer ranges from 0.5% to 2% in Kudaliar and between 0.3% and 1% in Anantapur, which agrees with earlier works. Results show that (a) depending on the geology and on the structure of the weathering profile, the vertical distribution of effective porosity can be very different and that the fractured layers in crystalline aquifers are not necessarily characterized by a rapid decrease in effective porosity and (b) that the lateral variations in effective porosity can be larger than the vertical ones. These variations suggest that within a same weathering profile, the density of open fractures and/or degree of weathering in the fractured zone may significantly vary from a place to another. The proposed method provides information on the spatial distribution of effective porosity that is of prime interest in terms of flux and contaminant transport in crystalline aquifers. Implications for mapping groundwater storage and scarcity are also discussed, which should help in improving groundwater resource management strategies.     

The effective stress concept and the evaluation of changes in pore air pressure under jacketed isotropic compression tests for dry sand based on the 2‐phase porous theory

The effective stress concept for solid‐fluid 2‐phase media was revisited in this work. In particular, the effects of the compressibility of both the pore fluid and the soil particles were studied under 3 different conditions, i.e., undrained, drained, and unjacketed conditions based on a Biot‐type theory for 2‐phase porous media. It was confirmed that Terzaghi effective stress holds at the moment when soil grains are assumed to be incompressible and when the compressibility of the pore fluid is small enough compared to that of the soil skeleton. Then, isotropic compression tests for dry sand under undrained conditions were conducted within the triaxial apparatus in which the changes in the pore air pressure could be measured. The ratio of the increment in the cell pressure to the increment in the pore air pressure, m, corresponds to the inverse of the B value by Bishop and was obtained during the step loading of the cell pressure. In addition, the m values were evaluated by comparing them with theoretically obtained values based on the solid‐fluid 2‐phase mixture theory. The experimental m values were close to the theoretical values, as they were in the range of approximately 40 to 185, depending on the cell pressure. Finally, it was found that the soil material with a highly compressible pore fluid, such as air, must be analyzed with the multi‐phase porous mixture theory. However, Terzaghi effective stress is practically applicable when the compressibilities of both the soil particles and the pore fluid are small enough compared to that of the soil skeleton.     

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.     

On the dependence of the saturated hydraulic conductivity upon the effective porosity through a power law model at different scales

We study the scale dependence of the saturated hydraulic conductivity K_s through the effective porosity n_e by means of a newly developed power‐law model (PLM) which allows to use simultaneously measurements at different scales. The model is expressed as product between a single PLM (capturing the impact of the dominating scale) and a characteristic function κ^? accounting for the correction because of the other scale(s). The simple (closed form) expression of the κ^?‐function enables one to easily identify the scales which are relevant for K_s. The proposed model is then applied to a set of real data taken at the experimental site of Montalto Uffugo (Italy), and we show that in this case two (i.e. laboratory and field) scales appear to be the main ones. The implications toward an important application (solute transport) in Hydrology are finally discussed. 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.     

A smoothed cap model for variably saturated soils and its robust numerical implementation

This paper deals with the numerical implementation of a cap model for unsaturated soils. It provides a brief review of existing cap model approaches, based on which an improved model formulated in terms of generalised effective stress and matric suction is derived and described in detail. Although the proposed model is a multisurface plasticity model, it can efficiently be implemented using only single‐surface projections because of the smoothness of the model, which is obtained by construction. Numerical algorithms are provided for these single‐surface stress projections, using a single‐equation approach whenever possible. The robustness of the utilised single‐equation approaches is enhanced by proposing problem‐fitted start‐up procedures based on investigations of the nonlinear projection equations. A comparison of the model response with extensive material test data is used to validate the model and to demonstrate the robust application of the approach to silty sands and low to medium plasticity clays. Copyright © 2015 John Wiley & Sons, Ltd.     

Predicting the dependency of a degree of saturation on void ratio and suction using effective stress principle for unsaturated soils

The paper presents an approach to predicting variation of a degree of saturation in unsaturated soils with void ratio and suction. The approach is based on the effective stress principle for unsaturated soils and several underlying assumptions. It focuses on the main drying and wetting processes and does not incorporate the effects of hydraulic hysteresis. It leads to the dependency of water retention curve (WRC) on void ratio, which does not require any material parameters apart from the parameters specifying WRC for the reference void ratio. Its validity is demonstrated by comparing predictions with the experimental data on four different soils taken over from the literature. Good correlation between the measured and predicted behaviour indirectly supports applicability of the effective stress principle for unsaturated soils. Copyright © 2009 John Wiley & Sons, Ltd.     

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.     

History of faulting in the Northern Hastings Block,southern New England Orogen

Brittle failure is common in the Devonian to Permian rocks in the Northern Hastings Block (NHB) and is manifested by faults of different orientation and kinematic histories, but the timing of fault movement is not well defined. In this study, faults in the NHB were analysed with the map pattern of cross-cutting faults used to estimate the relative time of movement and relationship to other faults. We defined five episodes of faulting or fault reactivation that affected the NHB. The Yarras Fault System on the southwestern side of the NHB and the Parrabel Fault and related faults on the eastern side of the NHB are the two major fault systems responsible for transporting and rotating the NHB in the late Carboniferous. Faults on the eastern, northeastern and northern part of Parrabel Dome started and stopped moving after emplacement of the Hastings Block and before the intrusion of the Werrikimbe Triassic granitoids. We suggested that the movement on the major bounding faults is related to the accommodation of the NHB to the folding and cleavage development in the adjoining Nambucca Block, and is associated with the earliest part of the Hunter–Bowen Orogeny. Limited dextral movement on the extensions of the Taylors Arm Fault System caused minor displacements in the northeastern part of the NHB during the Late Triassic. Some small faults cut the Triassic granitoids or Triassic Lorne Basin sediments indicating tectonic activity continued post-Triassic.     

A modeling study of effective radiative forcing and climate response due to tropospheric ozone

This study simulates the effective radiative forcing(ERF) of tropospheric ozone from 1850 to 2013 and its effects on global climate using an aerosol–climate coupled model, BCC AGCM2.0.1 CUACE/Aero, in combination with OMI(Ozone Monitoring Instrument) satellite ozone data. According to the OMI observations, the global annual mean tropospheric column ozone(TCO) was 33.9 DU in 2013, and the largest TCO was distributed in the belts between 30°N and 45°N and at approximately 30°S; the annual mean TCO was higher in the Northern Hemisphere than that in the Southern Hemisphere;and in boreal summer and autumn, the global mean TCO was higher than in winter and spring. The simulated ERF due to the change in tropospheric ozone concentration from 1850 to 2013 was 0.46 W m~(-2), thereby causing an increase in the global annual mean surface temperature by 0.36℃, and precipitation by 0.02 mm d~(-1)(the increase of surface temperature had a significance level above 95%). The surface temperature was increased more obviously over the high latitudes in both hemispheres, with the maximum exceeding 1.4?C in Siberia. There were opposite changes in precipitation near the equator,with an increase of 0.5 mm d~(-1)near the Hawaiian Islands and a decrease of about-0.6 mm d~(-1)near the middle of the Indian Ocean.