Real - time prediction of rainfall induced instability of residual soil slopes associated with deep cracks

The early warning based on real-time prediction of rain-induced instability of natural residual slopes helps to minimise human casualties due to such slope failures. Slope instability prediction is complicated, as it is influenced by many factors, including soil properties, soil behaviour, slope geometry, and the location and size of deep cracks in the slope. These deep cracks can facilitate rainwater infiltration into the deep soil layers and reduce the unsaturated shear strength of residual soil. Subsequently, it can form a slip surface, triggering a landslide even in partially saturated soil slopes. Although past research has shown the effects of surface-cracks on soil stability, research examining the influence of deep-cracks on soil stability is very limited. This study aimed to develop methodologies for predicting the real-time rain-induced instability of natural residual soil slopes with deep cracks. The results can be used to warn against potential rain-induced slope failures. The literature review conducted on rain induced slope instability of unsaturated residual soil associated with soil crack, reveals that only limited studies have been done in the following areas related to this topic: - Methods for detecting deep cracks in residual soil slopes. - Practical application of unsaturated soil theory in slope stability analysis. - Mechanistic methods for real-time prediction of rain induced residual soil slope instability in critical slopes with deep cracks. Two natural residual soil slopes at Jombok Village, Ngantang City, Indonesia, which are located near a residential area, were investigated to obtain the parameters required for the stability analysis of the slope. A survey first identified all related field geometrical information including slope, roads, rivers, buildings, and boundaries of the slope. Second, the electrical resistivity tomography (ERT) method was used on the slope to identify the location and geometrical characteristics of deep cracks. The two ERT array models employed in this research are: Dipole-dipole and Azimuthal. Next, bore-hole tests were conducted at different locations in the slope to identify soil layers and to collect undisturbed soil samples for laboratory measurement of the soil parameters required for the stability analysis. At the same bore hole locations, Standard Penetration Test (SPT) was undertaken. Undisturbed soil samples taken from the bore-holes were tested in a laboratory to determine the variation of the following soil properties with the depth: - Classification and physical properties such as grain size distribution, atterberg limits, water content, dry density and specific gravity. - Saturated and unsaturated shear strength properties using direct shear apparatus. - Soil water characteristic curves (SWCC) using filter paper method. - Saturated hydraulic conductivity. The following three methods were used to detect and simulate the location and orientation of cracks in the investigated slope: (1) The electrical resistivity distribution of sub-soil obtained from ERT. (2) The profile of classification and physical properties of the soil, based on laboratory testing of soil samples collected from bore-holes and visual observations of the cracks on the slope surface. (3) The results of stress distribution obtained from 2D dynamic analysis of the slope using QUAKE/W software, together with the laboratory measured soil parameters and earthquake records of the area. It was assumed that the deep crack in the slope under investigation was generated by earthquakes. A good agreement was obtained when comparing the location and the orientation of the cracks detected by Method-1 and Method-2. However, the simulated cracks in Method-3 were not in good agreement with the output of Method-1 and Method-2. This may have been due to the material properties used and the assumptions made, for the analysis. From Method-1 and Method-2, it can be concluded that the ERT method can be used to detect the location and orientation of a crack in a soil slope, when the ERT is conducted in very dry or very wet soil conditions. In this study, the cracks detected by the ERT were used for stability analysis of the slope. The stability of the slope was determined using the factor of safety (FOS) of a critical slip surface obtained by SLOPE/W using the limit equilibrium method. Pore-water pressure values for the stability analysis were obtained by coupling the transient seepage analysis of the slope using finite element based software, called SEEP/W. A parametric study conducted on the stability of an investigated slope revealed that the existence of deep cracks and their location in the soil slope are critical for its stability. The following two steps are proposed to predict the rain-induced instability of a residual soil slope with cracks. (a) Step-1: The transient stability analysis of the slope is conducted from the date of the investigation (initial conditions are based on the investigation) to the preferred date (current date), using measured rainfall data. Then, the stability analyses are continued for the next 12 months using the predicted annual rainfall that will be based on the previous five years rainfall data for the area. (b) Step-2: The stability of the slope is calculated in real-time using real-time measured rainfall. In this calculation, rainfall is predicted for the next hour or 24 hours and the stability of the slope is calculated one hour or 24 hours in advance using real time rainfall data. If Step-1 analysis shows critical stability for the forthcoming year, it is recommended that Step-2 be used for more accurate warning against the future failure of the slope. In this research, the results of the application of the Step-1 on an investigated slope (Slope-1) showed that its stability was not approaching a critical value for year 2012 (until 31st December 2012) and therefore, the application of Step-2 was not necessary for the year 2012. A case study (Slope-2) was used to verify the applicability of the complete proposed predictive method. A landslide event at Slope-2 occurred on 31st October 2010. The transient seepage and stability analyses of the slope using data obtained from field tests such as Bore-hole, SPT, ERT and Laboratory tests, were conducted on 12th June 2010 following the Step-1 and found that the slope in critical condition on that current date. It was then showing that the application of the Step-2 could have predicted this failure by giving sufficient warning time.

Spatial prediction of N2O emissions in pasture : a Bayesian model averaging analysis

Nitrous oxide (N2O) is one of the greenhouse gases that can contribute to global warming. Spatial variability of N2O can lead to large uncertainties in prediction. However, previous studies have often ignored the spatial dependency to quantify the N2O - environmental factors relationships. Few researches have examined the impacts of various spatial correlation structures (e.g. independence, distance-based and neighbourhood based) on spatial prediction of N2O emissions. This study aimed to assess the impact of three spatial correlation structures on spatial predictions and calibrate the spatial prediction using Bayesian model averaging (BMA) based on replicated, irregular point-referenced data. The data were measured in 17 chambers randomly placed across a 271 m(2) field between October 2007 and September 2008 in the southeast of Australia. We used a Bayesian geostatistical model and a Bayesian spatial conditional autoregressive (CAR) model to investigate and accommodate spatial dependency, and to estimate the effects of environmental variables on N2O emissions across the study site. We compared these with a Bayesian regression model with independent errors. The three approaches resulted in different derived maps of spatial prediction of N2O emissions. We found that incorporating spatial dependency in the model not only substantially improved predictions of N2O emission from soil, but also better quantified uncertainties of soil parameters in the study. The hybrid model structure obtained by BMA improved the accuracy of spatial prediction of N2O emissions across this study region.

Blast response and failure analysis of pile foundations subjected to surface explosion

This paper presents the response of pile foundations to ground shocks induced by surface explosion using fully coupled and non-linear dynamic computer simulation techniques together with different material models for the explosive, air, soil and pile. It uses the Arbitrary Lagrange Euler coupling formulation with proper state material parameters and equations. Blast wave propagation in soil, horizontal pile deformation and pile damage are presented to facilitate failure evaluation of piles. Effects of end restraint of pile head and the number and spacing of piles within a group on their blast response and potential failure are investigated. The techniques developed and applied in this paper and its findings provide valuable information on the blast response and failure evaluation of piles and will provide guidance in their future analysis and design.

Soil moisture monitoring at the field scale using neutron probe

Measurement of the moisture variation in soils is required for geotechnical design and research because soil properties and behavior can vary as moisture content changes. The neutron probe, which was developed more than 40 years ago, is commonly used to monitor soil moisture variation in the field. This study reports a full-scale field monitoring of soil moisture using a neutron moisture probe for a period of more than 2 years in the Melbourne (Australia) region. On the basis of soil types available in the Melbourne region, 23 sites were chosen for moisture monitoring down to a depth of 1500 mm. The field calibration method was used to develop correlations relating the volumetric moisture content and neutron counts. Observed results showed that the deepest “wetting front” during the wet season was limited to the top 800 to 1000 mm of soil whilst the top soil layer down to about 550mmresponded almost immediately to the rainfall events. At greater depths (550 to 800mmand below 800 mm), the moisture variations were relatively low and displayed predominantly periodic fluctuations. This periodic nature was captured with Fourier analysis to develop a cyclic moisture model on the basis of an analytical solution of a one-dimensional moisture flow equation for homogeneous soils. It is argued that the model developed can be used to predict the soil moisture variations as applicable to buried structures such as pipes.

Calibration functions for estimating soil moisture from GPR dielectric constant measurements

Ground-penetrating radar (GPR) is widely used for assessment of soil moisture variability in field soils. Because GPR does not measure soil water content directly, it is common practice to use calibration functions that describe its relationship with the soil dielectric properties and textural parameters. However, the large variety of models complicates the selection of the appropriate function. In this article an overview is presented of the different functions available, including volumetric models, empirical functions, effective medium theories, and frequency-specific functions. Using detailed information presented in summary tables, the choice for which calibration function to use can be guided by the soil variables available to the user, the frequency of the GPR equipment, and the desired level of detail of the output. This article can thus serve as a guide for GPR practitioners to obtain soil moisture values and to estimate soil dielectric properties.

Impact of a change in tillage and crop residue management practice on soil chemical and microbiological properties in a cereal-producing red duplex soil in NSW, Australia

The effect of a change of tillage and crop residue management practice on the chemical and micro-biological properties of a cereal-producing red duplex soil was investigated by superimposing each of three management practices (CC: conventional cultivation, stubble burnt, crop conventionally sown; DD: direct-drilling, stubble retained, no cultivation, crop direct-drilled; SI: stubble incorporated with a single cultivation, crop conventionally sown), for a 3-year period on plots previously managed with each of the same three practices for 14 years. A change from DD to CC or SI practice resulted in a significant decline, in the top 0-5 cm of soil, in organic C, total N, electrical conductivity, NH4-N, NO3-N, soil moisture holding capacity, microbial biomass and CO2 respiration as well as a decline in the microbial quotient (the ratio of microbial biomass C to organic C; P <0.05). In contrast, a change from SI to DD or CC practice or a change from CC to DD or SI practice had only negligible impact on soil chemical properties (P >0.05). However, there was a significant increase in microbial biomass and the microbial quotient in the top 0-5 cm of soil following the change from CC to DD or SI practice and with the change from SI to DD practice (P <0.05). Analysis of ester-linked fatty acid methyl esters (EL-FAMEs) extracted from the 0- to 5-cm and 5- to 10-cm layers of the soils of the various treatments detected changes in the FAME profiles following a change in tillage practice. A change from DD practice to SI or CC practice was associated with a significant decline in the ratio of fungal to bacterial fatty acids in the 0- to 5-cm soil (P <0.05). The results show that a change in tillage practice, particularly the cultivation of a previously minimum-tilled (direct-drilled) soil, will result in significant changes in soil chemical and microbiological properties within a 3-year period. They also show that soil microbiological properties are sensitive indicators of a change in tillage practice.

A comparative analysis of the conservative bridge transition design approaches

Differential settlement at the bridge approach between the deck and rail track on ground is often considered as a source of challenging technical and economical problem. This caused by the sudden stiffness changes between the bridge deck and the track on ground, and changes in soil stiffness of backfill and sub-grade with soil moisture content and loading history. To minimise the negative social and economic impacts due to poor performances of railway tracks at bridge transition zones, it is important, a special attention to be given at design, construction and maintenance stages. It is critically challenging to obtain an appropriate design solution for any given site condition and most of the existing conventional design approaches are unable to address the actual on-site behaviour due to their inherent assumptions of continuity and lack of clarifying of the local effects. An evaluation of existing design techniques is considered to estimate their contributions to a potential solution for bridge transition zones. This paper analyses five different approaches: the Chinese Standard, the European Standard with three different approaches, and the Australian approach. Each design approach is used to calculate the layer thicknesses, accounting critical design features such as the train speed, the axle load, the backfill subgrade condition, and the dynamic loading response. Considering correlation between track degradation and design parameters, this paper concludes that there is still a need of an optimised design approach for bridge transition zones.

Using the Gini coefficient with BIOLOG substrate utilisation data to provide an alternative quantitative measure for comparing bacterial soil communities

A measure quantifying unequal use of carbon sources, the Gini coefficient (G), has been developed to allow comparisons of the observed functional diversity of bacterial soil communities. This approach was applied to the analysis of substrate utilisation data obtained from using BIOLOG microtiter plates in a study which compared decomposition processes in two contrasting plant substrates in two different soils. The relevance of applying the Gini coefficient as a measure of observed functional diversity, for soil bacterial communities is evaluated against the Shannon index (H) and average well colour development (AWCD), a measure of the total microbial activity. Correlation analysis and analysis of variance of the experimental data show that the Gini coefficient, the Shannon index and AWCD provided similar information when used in isolation. However, analyses based on the Gini coefficient and the Shannon index, when total activity on the microtiter plates was maintained constant (i.e. AWCD as a covariate), indicate that additional information about the distribution of carbon sources being utilised can be obtained. We demonstrate that the Lorenz curve and its measure of inequality, the Gini coefficient, provides not only comparable information to AWCD and the Shannon index but when used together with AWCD encompasses measures of total microbial activity and absorbance inequality across all the carbon sources. This information is especially relevant for comparing the observed functional diversity of soil microbial communities.

Capacity of fatty acid profiles and substrate utilization patterns to describe differences in soil microbial communities associated with increased salinity or alkalinity at three locations in South Australia

Fatty acid methyl ester (FAME) profiles, together with Biolog substrate utilization patterns, were used in conjunction with measurements of other soil chemical and microbiological properties to describe differences in soil microbial communities induced by increased salinity and alkalinity in grass/legume pastures at three sites in SE South Australia. Total ester-linked FAMEs (EL-FAMEs) and phospholipid-linked FAMEs (PL-FAMEs), were also compared for their ability to detect differences between the soil microbial communities. The level of salinity and alkalinity in affected areas of the pastures showed seasonal variation, being greater in summer than in winter. At the time of sampling for the chemical and microbiological measurements (winter) only the affected soil at site 1 was significantly saline. The affected soils at all three sites had lower organic C and total N concentrations than the corresponding non-affected soils. At site 1 microbial biomass, CO 2-C respiration and the rate of cellulose decomposition was also lower in the affected soil compared to the non-affected soil. Biomarker fatty acids present in both the EL- and PL-FAME profiles indicated a lower ratio of fungal to bacterial fatty acids in the saline affected soil at site 1. Analysis of Biolog substrate utilization patterns indicated that the bacterial community in the affected soil at site 1 utilized fewer carbon substrates and had lower functional diversity than the corresponding community in the non-affected soil. In contrast, increased alkalinity, of major importance at sites 2 and 3, had no effect on microbial biomass, the rate of cellulose decomposition or functional diversity but was associated with significant differences in the relative amounts of several fatty acids in the PL-FAME profiles indicative of a shift towards a bacterial dominated community. Despite differences in the number and relative amounts of fatty acids detected, principal component analysis of the EL- and PL-FAME profiles were equally capable of separating the affected and non-affected soils at all three sites. Redundancy analysis of the FAME data showed that organic C, microbial biomass, electrical conductivity and bicarbonate-extractable P were significantly correlated with variation in the EL-FAME profiles, whereas pH, electrical conductivity, NH 4-N, CO 2-C respiration and the microbial quotient were significantly correlated with variation in the PL-FAME profiles. Redundancy analysis of the Biolog data indicated that cation exchange capacity and bicarbonate-extractable K were significantly correlated with the variation in Biolog substrate utilization patterns.

Effect of soil properties on the response of pile to underground explosion

This paper develops and presents a fully coupled non-linear finite element procedure to treat the response of piles to ground shocks induced by underground explosions. The Arbitrary Lagrange Euler coupling formulation with proper state material parameters and equations are used in the study. Pile responses in four different soil types, viz, saturated soil, partially saturated soil and loose and dense dry soils are investigated and the results compared. Numerical results are validated by comparing with those from a standard design manual. Blast wave propagation in soils, horizontal pile deformations and damages in the pile are presented. The pile damage presented through plastic strain diagrams will enable the vulnerability assessment of the piles under the blast scenarios considered. The numerical results indicate that the blast performance of the piles embedded in saturated soil and loose dry soil are more severe than those in piles embedded in partially saturated soil and dense dry soil. Present findings should serve as a benchmark reference for future analysis and design.

Blast response and sensitivity analysis of segmental tunnel

This research treated the response of underground transportation tunnels to surface blast loads using advanced computer simulation techniques. The influences of important parameters, such as tunnel material, geometrical configuration of segments and surrounding soil were investigated. The findings of this research offer significant new information on the blast performance of underground tunnels and will contribute towards future civil engineering applications.

Performance-Based Quality Assurance/Quality Control (QA/QC) Acceptance Procedures for In-Place Soil Testing Phase 3

One of the objectives of this study was to evaluate soil testing equipment based on its capability of measuring in-place stiffness or modulus values. As design criteria transition from empirical to mechanistic-empirical, soil test methods and equipment that measure properties such as stiffness and modulus and how they relate to Florida materials are needed. Requirements for the selected equipment are that they be portable, cost effective, reliable, a ccurate, and repeatable. A second objective is that the selected equipment measures soil properties without the use of nuclear materials.The current device used to measure soil compaction is the nuclear density gauge (NDG). Equipment evaluated in this research included lightweight deflectometers (LWD) from different manufacturers, a dynamic cone penetrometer (DCP), a GeoGauge, a Clegg impact soil tester (CIST), a Briaud compaction device (BCD), and a seismic pavement analyzer (SPA). Evaluations were conducted over ranges of measured densities and moistures.Testing (Phases I and II) was conducted in a test box and test pits. Phase III testing was conducted on materials found on five construction projects located in the Jacksonville, Florida, area. Phase I analyses determined that the GeoGauge had the lowest overall coefficient of variance (COV). In ascending order of COV were the accelerometer-type LWD, the geophone-type LWD, the DCP, the BCD, and the SPA which had the highest overall COV. As a result, the BCD and the SPA were excluded from Phase II testing.In Phase II, measurements obtained from the selected equipment were compared to the modulus values obtained by the static plate load test (PLT), the resilient modulus (MR) from laboratory testing, and the NDG measurements. To minimize soil and moisture content variability, the single spot testing sequence was developed. At each location, test results obtained from the portable equipment under evaluation were compared to the values from adjacent NDG, PLT, and laboratory MR measurements. Correlations were developed through statistical analysis. Target values were developed for various soils for verification on similar soils that were field tested in Phase III. The single spot testing sequence also was employed in Phase III, field testing performed on A-3 and A-2-4 embankments, limerock-stabilized subgrade, limerock base, and graded aggregate base found on Florida Department of Transportation construction projects. The Phase II and Phase III results provided potential trend information for future research—specifically, data collection for in-depth statistical analysis for correlations with the laboratory MR for specific soil types under specific moisture conditions. With the collection of enough data, stronger relationships could be expected between measurements from the portable equipment and the MR values. Based on the statistical analyses and the experience gained from extensive use of the equipment, the combination of the DCP and the LWD was selected for in-place soil testing for compaction control acceptance. Test methods and developmental specifications were written for the DCP and the LWD. The developmental specifications include target values for the compaction control of embankment, subgrade, and base materials.

Determination of imidacloprid in paddy water and soil by liquid chromatography electrospray ionization-tandem mass spectrometry

A method for the determination of imidacloprid in paddy water and soil was developed using liquid chromatography electrospray ionization-tandem mass spectrometry (LC/ESI-MS/MS). Separation of imidacloprid was carried out on a Shimadzu C18 column (150 mm × 4.6 mm, 4.6 μm) with an acetonitrile?water (50 : 50, v/v) mobile phase containing 0.1% of acetic acid. The flow rate was 0.3 mL/min in isocratic mode. The product ion at 209 m/z was selected for quantification in multiple-reaction monitoring scan mode. Imidacloprid residues in soil were extracted by a solid-liquid extraction method with acetonitrile. Water samples were filtered and directly injected for analysis without extraction. Detection limits of 0.5 μg/kg and 0.3 μg/L were achieved for soil and water samples, respectively. The method had recoveries of 90 ± 2% (n = 4) for soil samples and 100 ± 2% (n = 4) for water samples. A linear relationship was observed throughout the investigated range of concentrations (1-200 μg/L), with the correlation coefficients ranging from 0.999 to 1.000. © Pleiades Publishing, Ltd., 2010.

Analysis of parameter uncertainty and sensitivity in PCPF-1 modeling for predicting concentrations of rice herbicides

This paper demonstrates the procedures for probabilistic assessment of a pesticide fate and transport model, PCPF-1, to elucidate the modeling uncertainty using the Monte Carlo technique. Sensitivity analyses are performed to investigate the influence of herbicide characteristics and related soil properties on model outputs using four popular rice herbicides: mefenacet, pretilachlor, bensulfuron-methyl and imazosulfuron. Uncertainty quantification showed that the simulated concentrations in paddy water varied more than those of paddy soil. This tendency decreased as the simulation proceeded to a later period but remained important for herbicides having either high solubility or a high 1st-order dissolution rate. The sensitivity analysis indicated that PCPF-1 parameters requiring careful determination are primarily those involve with herbicide adsorption (the organic carbon content, the bulk density and the volumetric saturated water content), secondary parameters related with herbicide mass distribution between paddy water and soil (1st-order desorption and dissolution rates) and lastly, those involving herbicide degradations. © Pesticide Science Society of Japan.

Blast response and failure analysis of a segmented buried tunnel

Underground tunnels are vulnerable to terrorist attacks which can cause collapse of the tunnel structures or at least extensive damage, requiring lengthy repairs. This paper treats the blast impact on a reinforced concrete segmental tunnel buried in soil under a number of parametric conditions; soil properties, soil cover, distance of explosive from the tunnel centreline and explosive weight and analyses the possible failure patterns. A fully coupled Fluid Structure Interaction (FSI) technique incorporating the Arbitrary Lagrangian-Eulerian (ALE) method is used in this study. Results indicate that the tunnel in saturated soil is more vulnerable to severe damage than that buried in either partially saturated soil or dry soil. The tunnel is also more vulnerable to surface explosions which occur directly above the centre of the tunnel than those that occur at any equivalent distances in the ground away from the tunnel centre. The research findings provide useful information on modeling, analysis, overall tunnel response and failure patterns of segmented tunnels subjected to blast loads. This information will guide future development and application of research in this field.