Cracking risk of partially saturated porous media—Part I: Microporoelasticity model

Drying of deformable porous media results in their shrinkage, and it may cause cracking provided that shrinkage deformations are hindered by kinematic constraints. This is the motivation to develop a thermodynamics‐based microporoelasticity model for the assessment of cracking risk in partially saturated porous geomaterials. The study refers to 3D representative volume elements of porous media, including a two‐scale double‐porosity material with a pore network comprising (at the mesoscale) 3D mesocracks in the form of oblate spheroids, and (at the microscale) spherical micropores of different sizes. Surface tensions prevailing in all interfaces between solid, liquid, and gaseous matters are taken into account. To establish a thermodynamics‐based crack propagation criterion for a two‐scale double‐porosity material, the potential energy of the solid is derived, accounting—in particular—for mesocrack geometry changes (main original contribution) and for effective micropore pressures, which depend (due to surface tensions) on the pore radius. Differentiating the potential energy with respect to crack density parameter yields the thermodynamical driving force for crack propagation, which is shown to be governed by an effective macrostrain. It is found that drying‐related stresses in partially saturated mesocracks reduce the cracking risk. The drying‐related effective underpressures in spherical micropores, in turn, result in a tensile eigenstress of the matrix in which the mesocracks are embedded. This way, micropores increase the mesocracking risk. Model application to the assessment of cracking risk during drying of argillite is the topic of the companion paper (Part II). Copyright © 2009 John Wiley & Sons, Ltd.     

Distribution of deep near-inertial waves observed in the Kuroshio Extension

The distribution of deep near-inertial waves (NIWs) is investigated using data mainly from an array of 46 near-bottom acoustic current meter sensors spanning a 600 km × 600 km region as part of the Kuroshio Extension System Study during 2004–2006. The deep NIW distribution is interpreted in the context of both upper-layer and near-bottom mapped circulations. The wintertime-mean mixed-layer NIW energy input, modeled from observed wind stress, has the same range of values north and south of the Kuroshio Extension in this region. Yet, the wintertime-mean deep NIW energy distribution reveals a sharp factor-of-5 decrease from north to south of the Kuroshio jet. This direct observational evidence shows that the Kuroshio Extension blocks the equatorward propagation of NIWs. The NIW energy that does reach the sea floor within the subset of wintertime observations in the subtropical gyre arrives with patchy spatial and temporal distribution. Elevated NIW energy in deep water is associated with anticyclones in the deep barotropic flow and unassociated with upper layer eddies.     

The multiple chronological techniques applied to the Lake Suigetsu SG06 sediment core,central Japan

The varved sediment of Lake Suigetsu (central Japan) provides a valuable opportunity to obtain high‐resolution, multi‐proxy palaeoenvironmental data across the last glacial/interglacial cycle. In order to maximize the potential of this archive, a well‐constrained chronology is required. This paper outlines the multiple geochronological techniques being applied – namely varve counting, radiocarbon dating, tephrochronology (including argon–argon dating) and optically stimulated luminescence (OSL) – and the approaches by which these techniques are being integrated to form a single, coherent, robust chronology. Importantly, we also describe here the linkage of the floating Lake Suigetsu (SG06) varve chronology and the absolute (IntCal09 tree‐ring) time scale, as derived using radiocarbon data from the uppermost (non‐varved) portion of the core. This tie‐point, defined as a distinct (flood) marker horizon in SG06 (event layer B‐07–08 at 1397.4 cm composite depth), is thus derived to be 11 255 to 11 222 IntCal09 cal. years BP (68.2% probability range).     

Deriving TDS values in coarse sediments from long normal and electromagnetic logs

This paper presents an empirical relationship of quantitatively linked electromagnetic (EM) borehole recordings of the total dissolved solids (TDS) in pore water in the Quaternary deposits of the Belgian coastal plain. First, the long normal (LN) logs are linked to EM logs, then the already developed relationships between LN resistivity measurements and the TDS values are rewritten for EM recordings. The main parameter in these equations is the formation factor, which is derived from ground water analyses and LN logs through Archie's law. The EM recording has several advantages compared to the LN logs. The EM analysis allows measuring in PVC-cased wells and is not hindered by the invasion zone around the well. Furthermore, it has a high vertical resolution. LN logs can be measured only once, after drilling a well; EM recordings can be repeated several times in monitoring wells, which allows the gathering of time-dependent data over a complete vertical cross section. Such data could be obtained with LN logs only in wells with screens over the full-depth interval, which causes a hydraulic short circuit. This short circuit can result in a large artificial flow through the well between different levels, resulting in a salinity profile, which is no longer representative for the studied site. Remediation against short circuiting is a reduction of the screened interval, which strongly reduces the gathered information. The application of the derived equations is one of setting up a monitoring network along the Belgian coast to monitor the trend in salinity levels and comparing present salinity levels with older LN recordings to investigate the salinity changes in the last 30 years. Deep wells already present in the Belgiancoastal plain can then be used to monitor both the fresh water head changes and the salt water evolution. The technique has also been used for parameter identification for which real concentration measurements were needed.     

Motion Analysis in SAR Images of Unfocused Objects Using Time–Frequency Methods

Synthetic aperture radar (SAR) image formation processing assumes that the scene is stationary, and to focus an object, one coherently sums a large number of independent returns. Any target motion introduces phases that distort and/or translate the target's image. Target motion produces a smear primarily in the azimuth direction of the SAR image. Time-frequency (TF) modeling is used to analyze and correct the residual phase distortions. An interactive focusing algorithm based on TF modeling demonstrates how to correct the phase and to rapidly focus the mover. This is demonstrated on two watercraft observed in a SAR image. Then, two time-frequency representations (TFRs) are applied to estimate the motion parameters of the movers or refocus them or both. The first is the short-time Fourier transform, from which a velocity profile is constructed based on the length of the smear. The second TFR is the time-frequency distribution series, which is a robust derivative of the Wigner-Ville distribution that works well in this SAR environment. The smear is a modulated chirp, from which a velocity profile is plotted and the phase corrections are integrated to focus the movers. The relationship between these two methods is discussed. Both methods show good agreement on the example.     

南海西北部沉积物捕获器中颗石藻通量和属种组成变化记录

本文利用沉积物捕获器回收的连续时间序列沉降颗粒物样品,研究了南海西北部海南岸外陆坡外缘海域颗石藻的通量和属种组成变化及其环境意义。研究发现,东亚夏季风引起的海南岸外上升流以及冬季风导致的混合层深度增加都会刺激该区域颗石藻的生长,其中夏季上升流的影响尤为显著。从颗石粒属种组成来看,Florisphaera profunda是南海西北部中深层(1 000m)水体中主要的颗石粒优势种,此外Emiliania huxleyi,Gephyrocapsa oceanica和Gephyrocaps ericsonii等属种含量也颇高。然而在将颗石粒折算成颗石球数量后,Emiliania huxleyi则成为主要优势种。沉积物捕获器中F.profunda的相对百分含量与颗石粒总通量呈显著的负相关,表明F.profunda的相对含量与总的颗石藻生产力有着密切联系。本结果对于利用沉积物中的颗石藻记录,尤其是F.profunda相对含量变化,恢复和重建古海洋环境尤其是古海洋生产力的研究提供了直接的理论证据。     

Comparative study of tracer diffusion of HTO, Na and Cl in compacted kaolinite, illite and montmorillonite

The through-diffusion of HTO, ²²Na⁺ and ³⁶Cl⁻ in kaolinite, homo-ionic Na-illite and homo-ionic Na-montmorillonite was measured at a high degree of compaction as a function of the salt concentration in the ‘external solution’, i.e. in the solution in contact with the clay sample. The clays were chosen for this study because of their differences in the number and nature of ion exchange sites leading to different proportions of interlayer-, inter-particle and free pore water. It was found that the diffusive mass transfer of Na⁺ in Na-montmorillonite and Na-illite increased with decreasing external salt concentration, while the opposite trend was observed for the diffusion of Cl⁻. These trends are more pronounced in the case of Na-montmorillonite than in Na-illite, while almost no salt effect was observed for kaolinite. Similarly no salt effect was observed for the diffusion of HTO through all of the clays tested. These observations are in agreement with a conceptual model where it is assumed that cations diffuse preferentially in the interlayer or diffuse double-layer porosity, while anions are almost completely excluded from these regions. In the case of Na⁺ diffusion, the salt effects can be explained by an influence on the concentration gradient of diffusing cations, while in the case of Cl⁻ the external salt concentration has an effect on the accessible porosity. Effective diffusion coefficients of Cl⁻ fulfil the same relationship to porosity as those of the uncharged HTO, when using accessible porosities for such a comparison. Furthermore it is shown that pore diffusion coefficients for the three tracers are fairly well correlated with the respective diffusion coefficients in bulk water, if the effective diffusion coefficients for Na⁺ are derived from calculated tracer concentration gradients in the interlayer or diffuse double-layer porosities.     

Spatial modelling and uncertainty assessment of pyroclastic cover thickness in the Sorrento Peninsula

Assessing spatial variability of soil thickness is a critical issue for understanding and predicting slope processes. The present work was aimed at estimating the spatial scales at which the variation of pyroclastic cover thickness occurs in a sample area in the Sorrento Peninsula (Italy). Stochastic simulation was used to understand the spatial variability of pyroclastic cover thickness on Mount Pendolo and to assess its spatial uncertainty. In the study area, covering about 0.7 km², thickness measurements were collected using electrical resistivity tomography profiles, continuous core drillings and steel rod penetrometric tests. Variographic analysis revealed the occurrence of an anisotropic behaviour along the N50 and N140 directions. In the latter anisotropic direction, a nested variogram was fitted including (1) a long-range component which could be related to large-scale factors, like the curvature of the slope and contributing area and (2) a shorter scale variation which is probably associated with the occurrence of denudation processes or to the articulate cover/bedrock interface. To assess the spatial variability and uncertainty of pyroclastic cover thickness, a stochastic simulation algorithm was used and 500 equally probable images of cover thickness were yielded. The results showed that a better thickness distribution map can be drawn by simulating the data collected on the slope and at the footslope separately. The approach also allowed delineating the areas characterized by greater uncertainty, suggesting supplementary measurements to further improve the cover thickness distribution model, thus reducing the uncertainty.