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.

Fractal dimension of dried foods : a correlation between microstructure and porosity

The quality of dried food is affected by a number of factors including quality of raw material, initial microstructure, and drying conditions. The structure of the food materials goes through deformations due to the simultaneous effect of heat and mass transfer during the drying process. Shrinkage and changes in porosity, microstructure and appearance are some of the most remarkable features that directly influence overall product quality. Porosity and microstructure are the important material properties in relation to the quality attributes of dried foods. Fractal dimension (FD) is a quantitative approach of measuring surface, pore characteristics, and microstructural changes [1]. However, in the field of fractal analysis, there is a lack of research in developing relationship between porosity, shrinkage and microstructure of different solid food materials in different drying process and conditions [2-4]. Establishing a correlation between microstructure and porosity through fractal dimension during convective drying is the main objective of this work.

Effect of cell wall properties of plant tissue on porosity and shrinkage of dried apple

In this study cell wall properties; moisture distribution, stiffness, thickness and cell dimension have been taken into consideration. Cell wall stiffness dependent on complex combination of plant cell microstructures, composition and water holding capacity of the cell. In this work, some preliminary steps taken by investing cell wall properties of apple in order to predict change of porosity and shrinkage during drying. Two different types of apple cell wall characteristic were investigated to correlate with porosity and shrinkage after convective drying. A scanning electron microscope (SEM), 2N Intron, a pyncometer and image J software were used in order to measure and analyze cell characteristics, water dynamics, porosity and shrinkage. Cell stiffness of red delicious apple was found higher than granny smith apples. A significant relationship has found between cell wall characteristics and both heat and mass transfer. Consequently, evolution of porosity and shrinkage noticeably influenced during convective drying by the nature of cell wall. This study has brought better understanding of porosity and shrinkage of dried food stuff in microscopic (cell) level and would provide better insight to attain energy effective drying process and quality food stuff.

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.

Two-scale computational modelling of unsaturated flow in soils exhibiting small scale heterogeneities

Unsaturated water flow in soil is commonly modelled using Richards’ equation, which requires the hydraulic properties of the soil (e.g., porosity, hydraulic conductivity, etc.) to be characterised. Naturally occurring soils, however, are heterogeneous in nature, that is, they are composed of a number of interwoven homogeneous soils each with their own set of hydraulic properties. When the length scale of these soil heterogeneities is small, numerical solution of Richards’ equation is computationally impractical due to the immense effort and refinement required to mesh the actual heterogeneous geometry. A classic way forward is to use a macroscopic model, where the heterogeneous medium is replaced with a fictitious homogeneous medium, which attempts to give the average flow behaviour at the macroscopic scale (i.e., at a scale much larger than the scale of the heterogeneities). Using the homogenisation theory, a macroscopic equation can be derived that takes the form of Richards’ equation with effective parameters. A disadvantage of the macroscopic approach, however, is that it fails in cases when the assumption of local equilibrium does not hold. This limitation has seen the introduction of two-scale models that include at each point in the macroscopic domain an additional flow equation at the scale of the heterogeneities (microscopic scale). This report outlines a well-known two-scale model and contributes to the literature a number of important advances in its numerical implementation. These include the use of an unstructured control volume finite element method and image-based meshing techniques, that allow for irregular micro-scale geometries to be treated, and the use of an exponential time integration scheme that permits both scales to be resolved simultaneously in a completely coupled manner. Numerical comparisons against a classical macroscopic model confirm that only the two-scale model correctly captures the important features of the flow for a range of parameter values.

Monitoring phase behavior of sub- and supercritical CO2 confined in porous fractal silica with 85% porosity

Phase behavior of CO2 confined in porous fractal silica with volume fraction of SiO2 φs = 0.15 was investigated using small-angle neutron scattering (SANS) and ultrasmall-angle neutron scattering (USANS) techniques. The range of fluid densities (0<(FCO2)bulk<0.977 g/cm3) and temperatures (T=22 °C, 35 and 60 °C) corresponded to gaseous, liquid, near critical and supercritical conditions of the bulk fluid. The results revealed formation of a dense adsorbed phase in small pores with sizes D<40 A° at all temperatures. At low pressure (P <55 bar, (FCO2)bulk <0.2 g/cm3) the average fluid density in pores may exceed the density of bulk fluid by a factor up to 6.5 at T=22 °C. This “enrichment factor” gradually decreases with temperature, however significant fluid densification in small pores still exists at temperature T=60°C, i.e., far above the liquid-gas critical temperature of bulk CO2 (TC=31.1 °C). Larger pores are only partially filled with liquid-like adsorbed layer which coexists with unadsorbed fluid in the pore core. With increasing pressure, all pores become uniformly filled with the fluid, showing no measurable enrichment or depletion of the porous matrix with CO2.

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.

Hydrolysis of triasulfuron, metsulfuron-methyl and chlorsulfuron in alkaline soil and aqueous solutions

The hydrolysis of triasulfuron, metsulfuron-methyl and chlorsulfuron in aqueous buffer solutions and in soil suspensions at pH values ranging from 5.2 to 11.2 was investigated. Hydrolysis of all three compounds in both aqueous buffer and soil suspensions was highly pH-sensitive. The rate of hydrolysis was much faster in the acidic pH range (5.2-6.2) than under neutral and moderately alkaline conditions (8.2-9.4), but it increased rapidly as the pH exceeded 10.2. All three compounds degraded faster at pH 5.2 than at pH 11.2. Hydrolysis rates of all three compounds could be described well with pseudo-first-order kinetics. There were no significant differences (P =0.05) in the rate constants (k, day-1) of the three compounds in soil suspensions from those in buffer solutions within the pH ranges studied. A functional relationship based on the propensity of nonionic and anionic species of the herbicides to hydrolyse was used to describe the dependence of the 'rate constant' on pH. The hydrolysis involving attack by neutral water was at least 100-fold faster when the sulfonylurea herbicides were undissociated (acidic conditions) than when they were present as the anion at near neutral pH. In aqueous buffer solution at pH > 11, a prominent degradation pathway involved O-demethylation of metsulfuron-methyl to yield a highly polar degradate, and hydrolytic opening of the triazine ring. It is concluded that these herbicides are not likely to degrade substantially through hydrolysis in most agricultural (C) 2000 Society of Chemical Industry.

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.

Statistical analysis of reduction in tensile strength of cotton strips as a measure of soil microbial activity

The cotton strip assay (CSA) is an established technique for measuring soil microbial activity. The technique involves burying cotton strips and measuring their tensile strength after a certain time. This gives a measure of the rotting rate, R, of the cotton strips. R is then a measure of soil microbial activity. This paper examines properties of the technique and indicates how the assay can be optimised. Humidity conditioning of the cotton strips before measuring their tensile strength reduced the within and between day variance and enabled the distribution of the tensile strength measurements to approximate normality. The test data came from a three-way factorial experiment (two soils, two temperatures, three moisture levels). The cotton strips were buried in the soil for intervals of time ranging up to 6 weeks. This enabled the rate of loss of cotton tensile strength with time to be studied under a range of conditions. An inverse cubic model accounted for greater than 90% of the total variation within each treatment combination. This offers support for summarising the decomposition process by a single parameter R. The approximate variance of the decomposition rate was estimated from a function incorporating the variance of tensile strength and the differential of the function for the rate of decomposition, R, with respect to tensile strength. This variance function has a minimum when the measured strength is approximately 2/3 that of the original strength. The estimates of R are almost unbiased and relatively robust against the cotton strips being left in the soil for more or less than the optimal time. We conclude that the rotting rate X should be measured using the inverse cubic equation, and that the cotton strips should be left in the soil until their strength has been reduced to about 2/3.

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.

Model for internal porosity development of different shaped foods

Three particular geometrical shapes of foods were prepared from food materials. Cuboidal (aspect ratio = 1:1, 2:1, 3:1) , cylindrical (length: dameter = 1:1, 2:1, 3:1) and spheres were selected from potato, beans and peas respectively. Internal porosity was determined from solid density (theoretical)and particle density (experimental) during fluidised bed drying at different moisture contents. Solid density was calculated using formulae (conservation of mass and volume) already published in the literature by previous researchers. Determined porosity values were correlated with moisture ratio for different geometrical shapes.

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.

Impact of nitrification inhibitor (DMPP) on soil nitrous oxide emissions from an intensive broccoli production system in sub-tropical Australia

Vegetable cropping systems are often characterised by high inputs of nitrogen fertiliser. Elevated emissions of nitrous oxide (N2O) can be expected as a consequence. In order to mitigate N2O emissions from fertilised agricultural fields, the use of nitrification inhibitors, in combination with ammonium based fertilisers, has been promoted. However, no data is currently available on the use of nitrification inhibitors in sub-tropical vegetable systems. A field experiment was conducted to investigate the effect of the nitrification inhibitor 3,4-dimethylpyrazole phosphate (DMPP) on N2O emissions and yield from broccoli production in sub-tropical Australia. Soil N2O fluxes were monitored continuously (3 h sampling frequency) with fully automated, pneumatically operated measuring chambers linked to a sampling control system and a gas chromatograph. Cumulative N2O emissions over the 5 month observation period amounted to 298 g-N/ha, 324 g-N/ha, 411 g-N/ha and 463 g-N/ha in the conventional fertiliser (CONV), the DMPP treatment (DMPP), the DMMP treatment with a 10% reduced fertiliser rate (DMPP-red) and the zero fertiliser (0N), respectively. The temporal variation of N2O fluxes showed only low emissions over the broccoli cropping phase, but significantly elevated emissions were observed in all treatments following broccoli residues being incorporated into the soil. Overall 70–90% of the total emissions occurred in this 5 weeks fallow phase. There was a significant inhibition effect of DMPP on N2O emissions and soil mineral N content over the broccoli cropping phase where the application of DMPP reduced N2O emissions by 75% compared to the standard practice. However, there was no statistical difference between the treatments during the fallow phase or when the whole season was considered. This study shows that DMPP has the potential to reduce N2O emissions from intensive vegetable systems, but also highlights the importance of post-harvest emissions from incorporated vegetable residues. N2O mitigation strategies in vegetable systems need to target these post-harvest emissions and a better evaluation of the effect of nitrification inhibitors over the fallow phase is needed.