Electric imaging and laboratory resistivity testing for geotechnical investigation of Pusan clay deposits

Large-scale reclamation works for new land development in the Nakdong River plain have been extensively carried out on soft clays. Several geotechnical characteristics of the clays could not be well evaluated, partly due to easy disturbance during drilling and sampling. Consequently, geophysical methods, seen as nondestructive testing tools, have been applied in geotechnical investigation of Pusan clays for the first time. In this study, the 2D electric imaging technique was employed to map the thick soft clay deposits in four reclamation sites. The Pusan clay deposit was very well mapped. Electric resistivity of Pusan clays was measured on over 50 core samples in the laboratory, and then correlated with other geotechnical parameters such as salinity, organic content, water content, plasticity, unit weight and sampling depth. Additionally, electric resistivity of about 20 natural clays collected worldwide was measured and compared to that of the Pusan clays as an initial effort in creating a database of clay electric resistivity to help further application of electric imaging in geotechnical investigation of clayey soils.     

Solute transport and water content measurements in clay soils using time domain reflectometry

Abstract

Clayey and saline soils have been shown to be problematic for time domain reflectometry (TDR) measurements. This study presents some of these problems and discusses solutions to them. Thirteen solute transport experiments were carried out in three undisturbed soil columns of swelling clay soil from Tunisia, labelled S1, S2, and S3 respectively. The columns were collected at three different physiographical regions within a catchment. Water fluxes ranged from 1.2 to 7.2 cm day^?1. The large solute transport heterogeneity and large tailing indicated that preferential flow was most pronounced in S1. The preferential flow took place in voids between structural elements and in wormholes. In S3, preferential flow was also evident, but not to the same extent as in S1. In S2, the solute transport was more uniform with little preferential flow. The heterogeneity of the solute transport increased with the water flux in S1 and to a smaller extent in S3, whereas it remained constant in S2. In a previous dye experiment in the field, preferential flow in cracks was observed at those sites where S1 and S3 were collected. In the column experiments, preferential flow in these cracks was less due to the higher initial water content compared to the dye experiments, indicating that the desiccation cracks were closed by the swelling clay.     

The use of laboratory measurements of electric resistivity of rocks for evaluating the magnetotelluric method

Summary A profile model of electric resistivity as a function of depth was compiled, the initial parameters being the values of the electric resistivity determined experimentally for basaltic and eclogitic rocks from the Bohemian Massif. The curves of the apparent resistivity were computed for this model and the measure of information evaluated for the various model layers. The dominant influence of the subsurface layer was proved.     

探地雷达测定土壤含水量的研究进展

结合国内外关于探地雷达测定土壤含水量的最新研究成果,阐述了探地雷达测定土壤含水量的基本原理和计算方法,总结和评价了土壤含水量与土壤介电常数的不同计算模型,针对在实际应用中存在的不足,提出了相关改进措施,并展望发展方向.     

Error distribution modelling of satellite soil moisture measurements for hydrological applications

Satellite‐based soil moisture data accuracies are of important concerns by hydrologists because they could significantly influence hydrological modelling uncertainty. Without proper quantification of their uncertainties, it is difficult to optimize the hydrological modelling system and make robust decisions. Currently, the satellite soil moisture data uncertainty has been limited to summary statistics with the validations mainly from the in situ measurements. This study attempts to build the first error distribution model with additional higher‐order uncertainty modelling for satellite soil moisture observations. The methodology is demonstrated by a case study using the Soil Moisture and Ocean Salinity satellite soil moisture observations. The validation is based on soil moisture estimates from hydrological modelling, which is more relevant to the intended data use than the in situ measurements. Four probability distributions have been explored to find suitable error distribution curves using the statistical tests and bootstrapping resampling technique. General extreme value is identified as the most suitable one among all the curves. The error distribution model is still in its infant stage, which ignores spatial and temporal correlations, and nonstationarity. Further improvements should be carried out by the hydrological community by expanding the methodology to a wide range of satellite soil moisture data using different hydrological models. Copyright © 2016 John Wiley & Sons, Ltd.     

Estimation of soil water content in watershed using artificial neural networks

Soil water content (SWC) is an important factor in transfer processes between soil and air, contributing to water and energy balances, and quantifying it remains a challenge. This study uses artificial neural networks (ANNs) to analyse spatial and temporal variation of SWC in a Brazilian watershed, based on climate information, soil physical properties and topographic variables. Thirty eight input variables were tested in 200 models. The outputs were compared with 650 gravimetric moisture measurements collected at 26 points (25 field studies). The results show that it is possible to estimate SWC efficiently (Nash-Sutcliffe statistic, NS = 0.77) using topographic data, soil physical properties and rainfall. If only climate information is considered, modelling is less efficient (NS = 0.28). Using many variables does not necessarily improve performance. Alternatively, SWC can be estimated by simplified models using rainfall and topographic maps information, although the performance is less good (NS = 0.65).     

The effect of clay mineral content on the dynamic response of reconstituted fine grained soil

Fine grained soils with considerable amount of silt may exhibit sand-like or clay-like behavior depending on several factors such as the amount of fines and clay content, as well as the consistency limits, other variables being kept unchanged. This unpredictable behavior makes silts highly problematic, especially under seismic conditions. This paper describes the laboratory behavior of low plasticity Adapazari silt, known to be highly sensitive to cyclic loading. In the first phase of the basic study reported herein, Adapazari silt was mixed with different percentages of bentonite and kaolin and the behavior of these reconstituted mixtures was investigated in cyclic triaxial and dynamic simple shear tests. The purpose was to identify basic index properties and their threshold values to delineate sand- and clay-like behavior. Such a distinction may make it possible to complement field penetration resistance with appropriate adjustment factors to evaluate the pore pressure development potential, thus the risk of ground failure during an earthquake. The results show that there is a range of liquid limit and plasticity index values above which cyclic failure is significantly mitigated. It can now be stated that silts of intermediate and high plasticity may be deemed of relatively low potential for ground failure during seismic loading.     

Two-dimensional resistivity imaging and modelling in areas of complex geology

A system is described for the automatic measurement of electrical resistivity pseudo-sections. This comprises a linear array of up to 32 electrodes connected through a multicore cable to a computer controlled switching module and a resistivity meter. The processing of the measured sections to produce two-dimensional true resistivity images of the subsurface is briefly described. Some account is given of the capabilities and limitations of the technique. This is illustrated by a series of computed constant separation traverses for models of simple subsurface structures. Examples of processed images derived from sections measured in areas of relatively complex geology follow, a comparison being made of the interpretations obtained using an automatic imaging method and a manual iterative approach. It is concluded that with the equipment and software so far developed, in areas of modest subsurface geological complexity where some control is available and where the structures are essentially two-dimensional, then good approximations to the true geoelectric sections can be obtained down to depths of between 100 and 200 m.     

Effects of soil resistivity on currents induced on pipelines

The goal of cathodic protection is to prevent corrosion by maintaining buried pipelines at a constant potential with respect to the surrounding soil. In practice, however, the implementation is very complicated since many factors can contribute to the current flowing off the pipe. Design requires characterization of the parameters impacting the corrosion process, such as soil resistivity, size of the pipe and quality of the coating.In the present paper, we have studied the effect of geomagnetic fields on the pipe-induced currents considering it as an additional cause of corrosion. A theoretical method implemented to model the induced currents was tested in a previous work and the effect during disturbed days was quantified. This theoretical model indicated that the intensity of the current induced in a pipeline by the varying geomagnetic field depends on the intensity and rate of change of the field and the electrical resistivity of the soil. This induced current is in equilibrium with the host current and there is no current drainage between the pipeline and the host until, along the length of the pipeline, the host resistivity becomes different. At that point, current must flow between the pipe and host in order to establish a new equilibrium. It is this drainage current, flowing between the pipeline and the host, which causes corrosion problems.Following these results, experimental tests were performed in Tierra del Fuego. In this zone, a geophysical study was made to determine the discontinuities in soil resistivities and simultaneous measurements of the geomagnetic field and the drainage of current were recorded at different sites. The results obtained from the correlation of the data are consistent with the theoretical predictions.     

Estimation of soil hydraulic properties and their uncertainty: comparison between laboratory and field experiment

Often the soil hydraulic parameters are obtained by the inversion of measured data (e.g. soil moisture, pressure head, and cumulative infiltration, etc.). However, the inverse problem in unsaturated zone is ill‐posed due to various reasons, and hence the parameters become non‐unique. The presence of multiple soil layers brings the additional complexities in the inverse modelling. The generalized likelihood uncertainty estimate (GLUE) is a useful approach to estimate the parameters and their uncertainty when dealing with soil moisture dynamics which is a highly non‐linear problem. Because the estimated parameters depend on the modelling scale, inverse modelling carried out on laboratory data and field data may provide independent estimates. The objective of this paper is to compare the parameters and their uncertainty estimated through experiments in the laboratory and in the field and to assess which of the soil hydraulic parameters are independent of the experiment. The first two layers in the field site are characterized by Loamy sand and Loamy. The mean soil moisture and pressure head at three depths are measured with an interval of half hour for a period of 1 week using the evaporation method for the laboratory experiment, whereas soil moisture at three different depths (60, 110, and 200 cm) is measured with an interval of 1 h for 2 years for the field experiment. A one‐dimensional soil moisture model on the basis of the finite difference method was used. The calibration and validation are approximately for 1 year each. The model performance was found to be good with root mean square error (RMSE) varying from 2 to 4 cm³ cm^?3. It is found from the two experiments that mean and uncertainty in the saturated soil moisture (θ_s) and shape parameter (n) of van Genuchten equations are similar for both the soil types. Copyright © 2010 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.     

3D modelling and sensitivity in DC resistivity using charge density

A three‐dimensional (3D) electrical resistivity modelling code is developed to interpret surface and subsurface data. Based on the integral equation, it calculates the charge density caused by conductivity gradients at each interface of the mesh, allowing the estimation of the potential everywhere without the need to interpolate between nodes. Modelling generates a huge matrix, made up of Green's functions, which is stored by using the method of pyramidal compression. The potential is compared with the analytical and the numerical solutions obtained by finite‐difference codes for two models: the two‐layer case and the vertical contact case. The integral method is more accurate around the source point and at the limits of the domain for the potential calculation using a pole‐pole array. A technique is proposed to calculate the sensitivity (Jacobian) and Hessian matrices in 3D. The sensitivity is based on the derivative with respect to the block conductivity of the potential computed using the integral equation; it is only necessary to compute the electrical field at the source location. A direct extension of this technique allows the determination of the second derivatives. The technique is compared with the analytical solutions and with the calculation of the sensitivity according to the method using the inner product of the current densities calculated at the source and receiver points. Results are very accurate when the Green's function that includes the source image is used. The calculation of the three components of the electric field on the interfaces of the mesh is carried out simultaneously and quickly, using matrix compression.     

Ultrasonic wave propagation in dry and brine-saturated sandstones as a function of effective stress: laboratory measurements and modelling1

Compressional and shear-wave velocities have been measured and a novel approach using digital processing employed to study wave attenuation for brine- and gas- saturated sandstones, over a range of effective stresses from 5 to 60 MPa. Also measured were the complex conductivity in the brine-saturated state and permeability in the gas-saturated state over the same range of stresses as for the velocity measurements. Broadband ultrasonic pulses of P- and orthogonally polarized S-waves in the frequency range 0.3–0.8 MHz are transmitted through the specimen to be characterized for comparison with a reference (aluminium) having low attenuation. The attenuation is calculated in terms of the quality factor Q from the Fourier spectral ratios, using the frequency spectral ratios technique. The corrections necessary for the effects of diffraction due to the finite size of the ultrasonic transducers have been carried out for the case of measurements under lower confining stress. To interpret the laboratory measured velocity and attenuation data under the physical conditions of this study and to estimate the effects of pore structure, numerical modelling of velocities and attenuation as functions of the confining stress have been performed, based on the MIT model. Theoretical models based on several hypothesized attenuation mechanisms are considered in relation to laboratory data of the effects of confining pressure, fluid saturation and pore structure on attenuation. Numerical calculations using these models with the experimental data indicate that friction on thin cracks and grain boundaries is the dominant attenuation mechanism for dry and brine-saturated sandstones at low effective stresses for the frequencies tested. However, for brine-saturated sandstones at moderately high effective stresses, fluid flow could play a more important role in ultrasonic S-wave attenuation, depending on the pore structure of the sample.     

Seasonal rainfall partitioning under runon and runoff conditions on sandy soil in Niger. On-farm measurements and water balance modelling

In sparsely cropped farming systems in semi-arid tropics, rainfall partitioning can be complex due to various interactions between vertical and horizontal water flows, both in the atmosphere and in the soil. Despite this, quantifying the seasonal rainfall partitioning is essential, in order to identify options for increased yields. Results are presented on water flow components, based on field measurements and water balance modelling, for three years (1994–96) in a farmer's field cultivated with pearl millet [Pennisetum glaucum (L.) Br.] in the Sahel (Niger). Water balance modelling was carried out for three common infiltration categories: runoff producing surfaces, surfaces receiving inflow of runon water from upstream zones, and a reference surface with zero runoff and runon. Runoff was calculated to 25%–30% of annual rainfall (which ranged from 488 to 596 mm), from crust observations, rainfall, soil wetness data, and infiltration estimates. Inflow of runon was estimated from field observations to 8%–18% of annual rainfall. The parameters in the functions for soil surface and canopy resistances were calibrated with field measurements of soil evaporation, stomatal conductance and leaf area. The model estimates of soil water contents, which were validated against neutron probe measurements, showed a reasonable agreement with observed data, with a root mean square error (RMSE) of approximately 0.02 m³ m⁻³ for 0–160 cm soil depth. Estimated productive water flow as plant transpiration was low, amounting to 4%–9% of the available water for the non-fertilised crop and 7%–24% for the fertilised crop. Soil evaporation accounted for 31%–50% of the available water, and showed a low variation for the observed range of leaf area (LAI <1 m² m⁻²). Deep percolation was high, amounting to 200–330 mm for the non-crusted surfaces, which exceeded soil evaporation losses, for 1994–95 with relatively high annual rainfall (517–596 mm). Even a year with lower rainfall (488 mm) and a distinct dry spell during flowering (1996), resulted in an estimated deep percolation of 160 mm for the non-fertilised crop. The crop did not benefit from the additional inflow of runon water, which was partitioned between soil water storage and deep percolation. The only exception to this was the fertilised crop in 1996, where runon somewhat compensated for the limited rainfall and the higher water demand as a result of a larger leaf area than the non-fertilised crop. The effects of rainfall erraticness, resulting in episodic droughts, explain why a crop that uses such a small proportion of the available water, in an environment with substantial deep percolation, still suffers from water scarcity. Application of small levels of phosphorus and nitrogen roughly doubled yields, from 380 to 620 kg ha⁻¹, and plant transpiration, from 33 to 78 mm. Evapotranspirational water use efficiency (WUE_ET) was low, 6500–8300 m³ ton⁻¹ grain for non-fertilised crop, which is an effect of the low on-farm yields and high non-productive water losses. The estimated seasonal rainfall partitioning indicates the possibility of quantifying vertical water flows in on-farm environments in the Sahel, despite the presence of surface overland flow.     

Model tank experiments for resistivity measurements on non-conducting and conducting sheets

Summary Model tank experiments for conducting and non-conducting sheets have been described and the results compared with those of other workers.     

不同测向地电阻率测值及映震能力初探

对南京地震台、新沂地震台和郫县地震台地电阻率观测资料变化特征及邻区4次MS 5．0以上地震反应进行研究,发现：同一台站、不同测向地电阻率的测值大小可能存在较大差异,测值大的测道地电阻率年变化较规则,往往具有更强的映震能力。     

Estimation of regional material yield from coastal landslides based on historical digital terrain modelling

High‐resolution historical (1942) and recent (1994) digital terrain models were derived from aerial photographs along the Big Sur coastline in central California to measure the long‐term volume of material that enters the nearshore environment. During the 52‐year measurement time period, an average of 21 000 ± 3100 m³ km^?1 a^?1 of material was eroded from nine study sections distributed along the coast, with a low yield of 1000 ± 240 m³ km^?1 a^?1 and a high of 46 700 ± 7300 m³ km^?1 a^?1. The results compare well with known volumes from several deep‐seated landslides in the area and suggest that the processes by which material is delivered to the coast are episodic in nature. In addition, a number of parameters are investigated to determine what influences the substantial variation in yield along the coast. It is found that the magnitude of regional coastal landslide sediment yield is primarily related to the physical strength of the slope‐forming material. Coastal Highway 1 runs along the lower portion of the slope along this stretch of coastline, and winter storms frequently damage the highway. The California Department of Transportation is responsible for maintaining this scenic highway while minimizing the impacts to the coastal ecosystems that are part of the Monterey Bay National Marine Sanctuary. This study provides environmental managers with critical background data on the volumes of material that historically enter the nearshore from landslides, as well as demonstrating the application of deriving historical digital terrain data to model landscape evolution. Published in 2005 by John Wiley & Sons, Ltd.     

Estimation of particle velocity based on blast event measurements at different rock units

This paper analyses the results of the peak particle velocity measurements, at the trial, construction and quarry blasting carried out in sediment rock deposits comprising mainly limestone and dolomite. The tests are divided into 3 groups depending on different geological strength index (GSI) values in these rock groups. Based on the results of seismic measurements the empirical relationships between peak particle velocity and scaled distance were established for each group. In order to establish a useful relationship between peak particle velocity and scaled distance, a simple regression analysis was conducted with the Blastware software program from Instantel. The established relationships and the results of those analyses are presented in this paper.     

基于介电特性获取污染土壤中污染物含量的研究

从污染土壤复介电常数与污染成分的关系模型、污染土壤复介电常数测量技术两个方面综述了目前国内外学者利用介电特性和电磁波信号开展的土壤污染探测研究.最后指出,基于污染土壤的介电特性获取污染土壤中污染物含量是一项非常有前景的技术,其目的是通过获取土壤复介电常数的变化信息,定性和定量的对污染场地做出评价.     

基于多种度量的电离层TEC混沌预测分析

利用东经120°,北纬45°上空的2008年年积日101～150 d时间段内共600个电离层格网TEC数据为时间序列样本,分析了该点上空电离层TEC参数的混沌特性,发现其关联维数为2.2632,嵌入维数为5,最大Lyapunov指数为λ=0.0833,表明该点上空电离层TEC时间序列具有混沌的特征,存在混沌现象.采用加权一阶局域法对TEC时间序列进行预测时,在有效预报尺度内,利用各种度量进行预报得到的中误差都随着预报天数的增加而逐渐增大,但预报的相对误差却是在一定幅度范围内上下波动.采用四种距离度量以及皮尔森相关系数和Jaccard相似系数进行预报,其预测结果或者平均中误差较大,或者平均相对误差较大,结果不理想.其中,采用标准化欧氏距离得到的平均中误差最大,达到5.107TECU,而采用皮尔森相关系数得到的平均中误差最小,为5.078TECU;对于平均相对误差而言,由切比雪夫距离得到的平均相对误差最大,为0.185,而Jaccard相似系数得到的平均相对误差最小,为0.166.而采用余弦相似度得到的预测结果,其平均中误差和平均相对误差都很小,分别为5.079TECU和0.167.因此,在一阶局域混沌预测时,采用余弦相似度作为衡量相空间轨迹的差异性是较为合适的.