The effective number of cycles of earthquake ground motion

The seismic response of any system that accumulates damage under cyclic loading is dependent not only on the maximum amplitude of the motion but also its duration. This is explicitly recognized in methods for estimating the liquefaction potential of soil deposits. Many researchers have proposed that the effective number of cycles of the ground motion is a more robust indicator of the destructive capacity of the shaking than the duration. However, as is the case with strong‐motion duration, there is no universally accepted approach to determining the effective number of cycles of motion, and the different methods that have been proposed can give widely varying results for a particular accelerogram. Definitions of the effective number of cycles of motion are reviewed, classified and compared. Measurements are found to differ particularly for accelerograms with broad‐banded frequency content, which contain a significant number of non‐zero crossing peaks. The key seismological parameters influencing the number of cycles of motion and associated equations for predicting this quantity for future earthquakes are identified. Correlations between cycle counts and different duration measures are explored and found to be rather poor in the absence of additional parameters. Copyright © 2004 John Wiley & Sons, Ltd.     

Penny-shaped fractures revisited

All methods of seismic characterization of fractured reservoirs are based on effective media theories that relate geometrical and material properties of fractures and surrounding rock to the effective stiffnesses. In exploration seismology, the first-order theory of Hudson is the most popular. It describes the effective model caused by the presence of a single set of thin, aligned vertical fractures in otherwise isotropic rock. This model is known to be transversely isotropic with a horizontal symmetry axis (HTI). Following the theory, one can invert the effective anisotropy for the crack density and type of fluid infill of fractures, the quantities of great importance for reservoir appraisal and management.Here I compute effective media numerically using the finite element method. I deliberately construct models that contain a single set of vertical, ellipsoidal, non-intersecting and non-interconnected fractures to check validity of the first-order Hudson’s theory and establish the limits of its applicability. Contrary to conventional wisdom that Hudson’s results are accurate up to crack density e ≈ 0.1, I show that they consistently overestimate the magnitudes of all effective anisotropic coefficients ε^(V), δ^(V), and γ^(V). Accuracy of theoretically derived anisotropy depends on the type of fluid infill and typically deteriorates as e grows. While the theory gives | ε^(V)|, |δ^(V)|, |γ^(V)| and close to the upper bound of the corresponding numerically obtained values for randomly distributed liquid-filled fractures, theoretical predictions of ε^(V), δ^(V) are not supported by numerical computations when the cracks are dry. This happens primarily because the first-order Hudson’s theory makes no attempt to account for fracture interaction which contributes to the final result much stronger for gas- than for liquid-filled cracks. I find that Mori-Tanaka’s theory is superior to Hudson’s for all examined crack densities and both types of fluid infill.The paper was presented at the 11th International Workshop on Seismic Anisotropy (11IWSA) held in St. John’s, Canada in 2004.     

Did Paleogene North Atlantic rift-related eruptions drive early Eocene climate cooling?

The delivery of volcanogenic sulphur into the upper atmosphere by explosive eruptions is known to cause significant temporary climate cooling. Therefore, phreatomagmatic and phreatoplinian eruptions occurring during the final rifting stages of active flood basalt provinces provide a potent mechanism for triggering climate change.

During the early Eocene, the northeast Atlantic margin was subjected to repeated ashfall for 0.5 m.y. This was the result of extensive phreatomagmatic activity along 3000 km of the opening northeast Atlantic rift. These widespread, predominantly basaltic ashes are now preserved in marine sediments of the Balder Formation and its equivalents, and occur over an area extending from the Faroe Islands to Denmark and southern England. These ash-bearing sediments also contain pollen and spore floras derived from low diversity forests that grew in cooler, drier climates than were experienced either before or after these highly explosive eruptions. In addition, coeval plant macrofossil evidence from the Bighorn Basin, Wyoming, USA, also shows a comparable pattern of vegetation change. The coincidence of the ashes and cooler climate pollen and spore floras in northwest Europe identifies volcanism as the primary cause of climate cooling. Estimates show that whilst relatively few phreatomagmatic eruptive centres along the 3000 km opening rift system could readily generate 0.5–1 °C cooling, on an annual basis, only persistent or repeated volcanic phases would have been able to achieve the long-term cooling effect observed in the floral record. We propose that the cumulative effect of repeated Balder Formation eruptions initiated a biodiversity crisis in the northeast Atlantic margin forests. Only the decline of this persistent volcanic activity, and the subsequent climatic warming at the start of the Eocene Thermal Maximum allowed the growth of subtropical forests to develop across the region.     

关于岩石圈有效弹性厚度的地质理解

简要回顾了岩石圈均衡理论的发展及岩石圈区域均衡和挠曲理论在岩石圈动力学研究中起的作用,阐述岩石圈有效弹性厚度(Te)的概念和特征。强调Te的研究是地质学和地球物理学的紧密结合,即通过岩石圈挠曲理论和区域均衡原理,对地形和重力资料进行谱分析计算,来获取岩石圈的物理性质信息。计算的Te(和相应的挠曲刚度)是岩石圈等效的强度,与爆破地震、地震层析成像和大地电磁测深等方法观测到的岩石圈和地壳厚度不同,它们之间只有通过岩石圈的屈服刚度包络面(YSE)才能比较。大洋和大陆岩石圈YSE的理论计算,表明Te值显著小于地震学的地壳和岩石圈厚度。尤其对于大陆岩石圈,地壳厚度、热年龄和应变率均可显著影响岩石圈的强度。本文还以滇西为实例介绍了对相干值曲线计算的新认识和当前岩石圈Te研究的最新趋势。     

Effect of suction on the mechanical behaviour of iron ore rock

The effect of suction on the behaviour of iron ore has been studied from both physical and mechanical points of view. The porosity and the suction phenomena have been analysed using different experimental techniques. Uniaxial compressive tests on partially saturated samples have shown that the suction is responsible for strength and cohesion improvement. Considering the theory of partially saturated porous soils of Coussy and Dangla (Mécanique des sols non saturés (2002 edn). Hermès Science: 2002; 390), we have proposed a constitutive law for partially saturated iron ore. The real increase in the apparent cohesion due to the capillary attraction forces is overestimated if the yield function is written in terms of effective stresses. The effect of the capillary cohesion has been modelled with a function in the expression of the apparent cohesion of the yield function. The effect of suction on the mechanical behaviour has been represented in the effective stresses space and in the total stresses space like the Alonso model (Géotechnique 1990; 40 :405–430). Copyright © 2005 John Wiley & Sons, Ltd.     

Precessional forcing of tropical vegetation carbon storage

Since the Mid Pleistocene Revolution, which occurred about one million years ago, global temperatures have fluctuated with a quasi‐periodicity of ca. 100 ka. The pattern of past change in the extent of woodlands, and therefore by inference vegetation carbon storage, has been demonstrated to have a strong positive link with this global temperature change at high and mid latitudes. However, understanding of climate systems and ecosystem function indicates that the pattern of woodland change at low latitudes may follow a fundamentally different pattern. We present output from the intermediate complexity model GENIE‐1, comprising a single transient simulation over the last 800 ka and a 174‐member ensemble of 130 ka transient simulations over the last glacial cycle. These simulations suggest that while vegetation carbon storage in mid–high northern latitudes robustly follows the characteristic ca. 100 ka cycle, this signal is not a robust feature of tropical vegetation, which is subject to stronger direct forcing by the precessional (21 ka) orbital cycle (albeit with a highly uncertain response). We conclude that the correlation of palaeoenvironmental records from low latitudes with global temperature change must be done with caution. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

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.     

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.